Expression of antisense small RNAs in response to stress in Pseudomonas aeruginosa

A Novel Regulation of K-antigen Capsule Synthesis in Porphyromonas gingivalis Is Driven by the Response Regulator PG0720-Directed Antisense RNA

The periodontal pathogen Porphyromonas gingivalis strain W83 displays at least three different surface glycans, specifically two types of lipopolysaccharides (O-LPS and A-LPS) and K-antigen capsule. Despite the importance of K-antigen capsule to the virulence of P. gingivalis, little is known as to how expression of genes involved in the synthesis of this surface glycan is regulated. The genes required for K-antigen capsule synthesis are located in a locus that encodes a number of transcripts, including an operon (PG0104 to PG0121, generating ~19.4-kb transcript) which contains a non-coding 77-bp inverted repeat (77 bpIR) region near the 5'-end. Previously, we identified a 550-nucleotide antisense RNA molecule (designated asSuGR for antisense Surface Glycan Regulator) encoded within the 77-bpIR element that influences the synthesis of surface glycans. In this study, we demonstrate that the DNA-binding response regulator PG0720 can bind the promoter region of asSuGR and activate expression of asSuGR, indicating that PG0720 may indirectly influence transcript levels of the K-antigen capsule operon expressed from the sense strand. The data show that deletion of the PG0720 gene confers a defect in the presentation of surface polysaccharides compared with the parent strain and quantitative RT-PCR (qPCR) analysis determined that the overall expression of genes involved in K-antigen capsule synthesis were down-regulated in the PG0720 mutant. Furthermore, the defects of the PG0720 deletion mutant were restored by complementation. Importantly, the PG0720 deletion mutant showed reduced virulence. Altogether, our data show that the response regulator PG0720 regulates expression of asSuGR, a trans-acting antisense RNA molecule involved in modulating the production of surface polysaccharides in P. gingivalis strain W83. The data provide further evidence that surface glycans are key virulence determinants and significantly advances our understanding of the molecular mechanisms controlling the synthesis of P. gingivalis K-antigen capsule, a key virulence determinant.

Direct RNA nanopore sequencing of Pseudomonas aeruginosa clone C transcriptomes

The transcriptomes of Pseudomonas aeruginosa clone C isolates NN2 and SG17M during the mid-exponential and early stationary phase of planktonic growth were evaluated by direct RNA sequencing on the nanopore platform and compared with established short-read cDNA sequencing on the Illumina platform. Fifty to ninety percent of the sense RNAs turned out to be rRNA molecules followed by similar proportions of mRNA transcripts and non-coding RNAs. Both platforms detected similar proportions of uncharged tRNAs and 29 yet undescribed antisense tRNAs. For example, the rarest arginine codon was paired with the most abundant tRNA^Arg, and the tRNA^Arg gene is missing for the most frequent arginine codon. More than 90% of the antisense RNA molecules were complementary to a coding sequence. The antisense RNAs were evenly distributed in the genomes. Direct RNA sequencing identified more than 4,000 distinct non-overlapping antisense RNAs during exponential and stationary growth. Besides highly expressed small antisense RNAs less than 200 bases in size, a population of longer antisense RNAs was sequenced that covered a broad range of a few hundred to thousands of bases and could be complementary to a contig of several genes. In summary, direct RNA sequencing identified yet undescribed RNA molecules and an unexpected composition of the pools of tRNAs, sense and antisense RNAs. IMPORTANCE Genome-wide gene expression of bacteria is commonly studied by high-throughput sequencing of size-selected cDNA fragment libraries of reverse-transcribed RNA preparations. However, the depletion of ribosomal RNAs, enzymatic reverse transcription and the fragmentation, size selection and amplification during library preparation lead to inevitable losses of information about the initial composition of the RNA pool. We demonstrate that direct RNA sequencing on the nanopore platform can overcome these limitations. Nanopore sequencing of total RNA yielded novel insights into the Pseudomonas aeruginosa transcriptome that - if replicated in other species - will change our view of the bacterial RNA world. The discovery of sense - antisense pairs of tmRNA, tRNAs and mRNAs indicates a further and unknown level of gene regulation in bacteria.

Oxidative Stress Response in Pseudomonas aeruginosa

Pseudomonas aeruginosa is a Gram-negative environmental and human opportunistic pathogen highly adapted to many different environmental conditions. It can cause a wide range of serious infections, including wounds, lungs, the urinary tract, and systemic infections. The high versatility and pathogenicity of this bacterium is attributed to its genomic complexity, the expression of several virulence factors, and its intrinsic resistance to various antimicrobials. However, to thrive and establish infection, P. aeruginosa must overcome several barriers. One of these barriers is the presence of oxidizing agents (e.g., hydrogen peroxide, superoxide, and hypochlorous acid) produced by the host immune system or that are commonly used as disinfectants in a variety of different environments including hospitals. These agents damage several cellular molecules and can cause cell death. Therefore, bacteria adapt to these harsh conditions by altering gene expression and eliciting several stress responses to survive under oxidative stress. Here, we used PubMed to evaluate the current knowledge on the oxidative stress responses adopted by P. aeruginosa. We will describe the genes that are often differently expressed under oxidative stress conditions, the pathways and proteins employed to sense and respond to oxidative stress, and how these changes in gene expression influence pathogenicity and the virulence of P. aeruginosa. Understanding these responses and changes in gene expression is critical to controlling bacterial pathogenicity and developing new therapeutic agents.

Global identification of RsmA/N binding sites in Pseudomonas aeruginosa by in vivo UV CLIP-seq

Pseudomonas aeruginosa harbours two redundant RNA-binding proteins RsmA/RsmN (RsmA/N), which play a critical role in balancing acute and chronic infections. However, in vivo binding sites on target transcripts and the overall impact on the physiology remains unclear. In this study, we applied in vivo UV crosslinking immunoprecipitation followed by RNA-sequencing (UV CLIP-seq) to detect RsmA/N-binding sites at single-nucleotide resolution and mapped more than 500 binding sites to approximately 400 genes directly bound by RsmA/N in P. aeruginosa. This also verified the ANGGA sequence in apical loops skewed towards 5'UTRs as a consensus motif for RsmA/N binding. Genetic analysis combined with CLIP-seq results suggested previously unrecognized RsmA/N targets involved in LPS modification. Moreover, the RsmA/N-titrating RNAs RsmY/RsmZ may be positively regulated by the RsmA/N-mediated translational repression of their upstream regulators, thus providing a possible mechanistic explanation for homoeostasis of the Rsm system. Thus, our study provides a detailed view of RsmA/N-RNA interactions and a resource for further investigation of the pleiotropic effects of RsmA/N on gene expression in P. aeruginosa.

A High-Throughput Method for Identifying Novel Genes That Influence Metabolic Pathways Reveals New Iron and Heme Regulation in Pseudomonas aeruginosa

The ability to simultaneously and more directly correlate genes with metabolite levels on a global level would provide novel information for many biological platforms yet has thus far been challenging. Here, we describe a method to help address this problem, which we dub “Met-Seq” (metabolite-coupled Tn sequencing).

Heme is an essential metabolite for most life on earth. Bacterial pathogens almost universally require iron to infect a host, often acquiring this nutrient in the form of heme. The Gram-negative pathogen Pseudomonas aeruginosa is no exception, where heme acquisition and metabolism are known to be crucial for both chronic and acute infections. To unveil unknown genes and pathways that could play a role with heme metabolic flux in this pathogen, we devised an omic-based approach we dubbed “Met-Seq,” for metabolite-coupled transposon sequencing. Met-Seq couples a biosensor with fluorescence-activated cell sorting (FACS) and massively parallel sequencing, allowing for direct identification of genes associated with metabolic changes. In this work, we first construct and validate a heme biosensor for use with P. aeruginosa and exploit Met-Seq to identify 188 genes that potentially influence intracellular heme levels. Identified genes largely consisted of metabolic pathways not previously associated with heme, including many secreted virulence effectors, as well as 11 predicted small RNAs (sRNAs) and riboswitches whose functions are not currently understood. We verify that five Met-Seq hits affect intracellular heme levels; a predicted extracytoplasmic function (ECF) factor, a phospholipid acquisition system, heme biosynthesis regulator Dnr, and two predicted antibiotic monooxygenase (ABM) domains of unknown function (PA0709 and PA3390). Finally, we demonstrate that PA0709 and PA3390 are novel heme-binding proteins. Our data suggest that Met-Seq could be extrapolated to other biological systems and metabolites for which there is an available biosensor, and provides a new template for further exploration of iron/heme regulation and metabolism in P. aeruginosa and other pathogens.

IMPORTANCE The ability to simultaneously and more directly correlate genes with metabolite levels on a global level would provide novel information for many biological platforms yet has thus far been challenging. Here, we describe a method to help address this problem, which we dub “Met-Seq” (metabolite-coupled Tn sequencing). Met-Seq uses the powerful combination of fluorescent biosensors, fluorescence-activated cell sorting (FACS), and next-generation sequencing (NGS) to rapidly identify genes that influence the levels of specific intracellular metabolites. For proof of concept, we create and test a heme biosensor and then exploit Met-Seq to identify novel genes involved in the regulation of heme in the pathogen Pseudomonas aeruginosa. Met-Seq-generated data were largely comprised of genes which have not previously been reported to influence heme levels in this pathogen, two of which we verify as novel heme-binding proteins. As heme is a required metabolite for host infection in P. aeruginosa and most other pathogens, our studies provide a new list of targets for potential antimicrobial therapies and shed additional light on the balance between infection, heme uptake, and heme biosynthesis.

Genome-wide identification of Azospirillum brasilense Sp245 small RNAs responsive to nitrogen starvation and likely involvement in plant-microbe interactions

BACKGROUND

Small RNAs (sRNAs) are non-coding RNAs known to regulate various biological functions such as stress adaptation, metabolism, virulence as well as pathogenicity across a wide range of bacteria, mainly by controlling mRNA stabilization or regulating translation. Identification and functional characterization of sRNAs has been carried out in various plant growth-promoting bacteria and they have been shown to help the cells cope up with environmental stress. No study has been carried out to uncover these regulatory molecules in the diazotrophic alpha-proteobacterium Azospirillum brasilense Sp245 to date.

RESULTS

Expression-based sRNA identification (RNA-seq) revealed the first list of ~ 468 sRNA candidate genes in A. brasilense Sp245 that were differentially expressed in nitrogen starvation versus non-starved conditions. In parallel, in silico tools also identified 2 of the above as candidate sRNAs. Altogether, putative candidates were stringently curated from RNA-seq data based on known sRNA parameters (size, location, secondary structure, and abundance). In total, ~ 59 significantly expressed sRNAs were identified in this study of which 53 are potentially novel sRNAs as they have no Rfam and BSRD homologs. Sixteen sRNAs were randomly selected and validated for differential expression, which largely was found to be in congruence with the RNA-seq data.

CONCLUSIONS

Differential expression of 468 A. brasilense sRNAs was indicated by RNA-seq data, a subset of which was confirmed by expression analysis. Four of the significantly expressed sRNAs were not observed in nitrogen starvation while 16 sRNAs were found to be exclusively expressed in nitrogen depletion. Putative candidate sRNAs identified have potential mRNA targets primarily involved in stress (abiotic and biotic) adaptability; regulation of bacterial cellular, biological and molecular pathways such as nitrogen fixation, polyhydroxybutyrate synthesis, chemotaxis, biofilm formation and transcriptional regulation. In addition to directly influencing bacteria, some of these sRNAs also have targets influencing plant-microbe interactions through adhesion of bacteria to plant roots directly, suppressing host response, inducing plant defence and signalling.

Full Transcriptomic Response of Pseudomonas aeruginosa to an Inulin-Derived Fructooligosaccharide

Pseudomonas aeruginosa is an ubiquitous gram-negative opportunistic human pathogen which is not considered part of the human commensal gut microbiota. However, depletion of the intestinal microbiota (Dysbiosis) following antibiotic treatment facilitates the colonization of the intestinal tract by Multidrug-Resistant P. aeruginosa. One possible strategy is based on the use of functional foods with prebiotic activity. The bifidogenic effect of the prebiotic inulin and its hydrolyzed form (fructooligosaccharide: FOS) is well established since they promote the growth of specific beneficial (probiotic) gut bacteria such as bifidobacteria. Previous studies of the opportunistic nosocomial pathogen Pseudomonas aeruginosa PAO1 have shown that inulin and to a greater extent FOS reduce growth and biofilm formation, which was found to be due to a decrease in motility and exotoxin secretion. However, the transcriptional basis for these phenotypic alterations remains unclear. To address this question we conducted RNA-sequence analysis. Changes in the transcript level induced by inulin and FOS were similar, but a set of transcript levels were increased in response to inulin and reduced in the presence of FOS. In the presence of inulin or FOS, 260 and 217 transcript levels, respectively, were altered compared to the control to which no polysaccharide was added. Importantly, changes in transcript levels of 57 and 83 genes were found to be specific for either inulin or FOS, respectively, indicating that both compounds trigger different changes. Gene pathway analyses of differentially expressed genes (DEG) revealed a specific FOS-mediated reduction in transcript levels of genes that participate in several canonical pathways involved in metabolism and growth, motility, biofilm formation, β-lactamase resistance, and in the modulation of type III and VI secretion systems; results that have been partially verified by real time quantitative PCR measurements. Moreover, we have identified a genomic island formed by a cluster of 15 genes, encoding uncharacterized proteins, which were repressed in the presence of FOS. The analysis of isogenic mutants has shown that genes of this genomic island encode proteins involved in growth, biofilm formation and motility. These results indicate that FOS selectively modulates bacterial pathogenicity by interfering with different signaling pathways.

A rhlI 5' UTR-Derived sRNA Regulates RhlR-Dependent Quorum Sensing in Pseudomonas aeruginosa

The opportunistic human pathogen Pseudomonas aeruginosa possesses multiple quorum sensing systems that regulate and coordinate production of virulence factors and adaptation to different environments. Despite extensive research, the regulatory elements that play a role in this complex network are still not fully understood. By using several RNA sequencing techniques, we were able to identify a small regulatory RNA we named RhlS. RhlS increases translation of RhlI, a key enzyme in the quorum sensing pathway, and represses the fpvA mRNA encoding one of the siderophore pyoverdine receptors. Our results highlight a new regulatory layer of P. aeruginosa quorum sensing and contribute to the growing understanding of the role regulatory RNAs play in bacterial physiology.

N-Acyl homoserine lactone (AHL) quorum sensing (QS) controls expression of over 200 genes in Pseudomonas aeruginosa. There are two AHL regulatory systems: the LasR-LasI circuit and the RhlR-RhlI system. We mapped transcription termination sites affected by AHL QS in P. aeruginosa, and in doing so we identified AHL-regulated small RNAs (sRNAs). Of interest, we noted that one particular sRNA was located within the rhlI locus. We found that rhlI, which encodes the enzyme that produces the AHL N-butanoyl-homoserine lactone (C4-HSL), is controlled by a 5′ untranslated region (UTR)-derived sRNA we name RhlS. We also identified an antisense RNA encoded opposite the beginning of the rhlI open reading frame, which we name asRhlS. RhlS accumulates as wild-type cells enter stationary phase and is required for the production of normal levels of C4-HSL through activation of rhlI translation. RhlS also directly posttranscriptionally regulates at least one other unlinked gene, fpvA. The asRhlS appears to be expressed at maximal levels during logarithmic growth, and we suggest RhlS may act antagonistically to the asRhlS to regulate rhlI translation. The rhlI-encoded sRNAs represent a novel aspect of RNA-mediated tuning of P. aeruginosa QS.

Within-Host Adaptation Mediated by Intergenic Evolution in Pseudomonas aeruginosa

Bacterial pathogens evolve during the course of infection as they adapt to the selective pressures that confront them inside the host. Identification of adaptive mutations and their contributions to pathogen fitness remains a central challenge. Although mutations can either target intergenic or coding regions in the pathogen genome, studies of host adaptation have focused predominantly on molecular evolution within coding regions, whereas the role of intergenic mutations remains unclear. Here, we address this issue and investigate the extent to which intergenic mutations contribute to the evolutionary response of a clinically important bacterial pathogen, Pseudomonas aeruginosa, to the host environment, and whether intergenic mutations have distinct roles in host adaptation. We characterize intergenic evolution in 44 clonal lineages of P. aeruginosa and identify 77 intergenic regions in which parallel evolution occurs. At the genetic level, we find that mutations in regions under selection are located primarily within regulatory elements upstream of transcriptional start sites. At the functional level, we show that some of these mutations both increase or decrease transcription of genes and are directly responsible for evolution of important pathogenic phenotypes including antibiotic sensitivity. Importantly, we find that intergenic mutations facilitate essential genes to become targets of evolution. In summary, our results highlight the evolutionary significance of intergenic mutations in creating host-adapted strains, and that intergenic and coding regions have different qualitative contributions to this process.

Pseudomonas aeruginosa Alters Its Transcriptome Related to Carbon Metabolism and Virulence as a Possible Survival Strategy in Blood from Trauma Patients

While a considerable body of knowledge regarding sepsis in trauma patients is available, the potential influence of trauma-induced changes in the blood of these patients on the pathogenesis of Pseudomonas aeruginosa is basically an unexplored area. Rather than using standard laboratory media, we grew P. aeruginosa in whole blood from either healthy volunteers or trauma patients. The specific changes in the P. aeruginosa transcriptome in response to growth in blood from trauma patients reflect the adaptation of this organism to the bloodstream environment. This knowledge is vital for understanding the strategies this pathogen uses to adapt and survive within the host during systemic infection. Such information will help researchers and clinicians to develop new approaches for treatment of sepsis caused by P. aeruginosa in trauma patients, especially in terms of recognizing the effects of specific therapies (e.g., iron, zinc, or mannitol) on the organism. Further, this information can most likely be extrapolated to all patients with P. aeruginosa septicemia.

Trauma patients (TPs) are highly susceptible to infections, which often lead to sepsis. Among the numerous causative agents, Pseudomonas aeruginosa is especially important, as P. aeruginosa sepsis is often fatal. Understanding the mechanism of its pathogenesis in bloodstream infections is imperative; however, this mechanism has not been previously described. To examine the effect of trauma-induced changes in blood on the expression of P. aeruginosa genes, we grew strain UCBPP-PA14 (PA14) in blood samples from eight TPs and seven healthy volunteers (HVs). Compared with its growth in blood from HVs, the growth of PA14 in blood from TPs significantly altered the expression of 285 genes. Genes whose expression was significantly increased were related to carbon metabolism, especially malonate utilization and mannitol uptake, and efflux of heavy metals. Genes whose expression was significantly reduced included genes of the type VI secretion system, genes related to uptake and metabolism of amino acids, and genes related to biosynthesis and transport of the siderophores pyoverdine and pyochelin. These results suggest that during systemic infection in trauma patients, and to adapt to the trauma-induced changes in blood, P. aeruginosa adjusts positively and negatively the expression of numerous genes related to carbon metabolism and virulence, respectively.

IMPORTANCE While a considerable body of knowledge regarding sepsis in trauma patients is available, the potential influence of trauma-induced changes in the blood of these patients on the pathogenesis of Pseudomonas aeruginosa is basically an unexplored area. Rather than using standard laboratory media, we grew P. aeruginosa in whole blood from either healthy volunteers or trauma patients. The specific changes in the P. aeruginosa transcriptome in response to growth in blood from trauma patients reflect the adaptation of this organism to the bloodstream environment. This knowledge is vital for understanding the strategies this pathogen uses to adapt and survive within the host during systemic infection. Such information will help researchers and clinicians to develop new approaches for treatment of sepsis caused by P. aeruginosa in trauma patients, especially in terms of recognizing the effects of specific therapies (e.g., iron, zinc, or mannitol) on the organism. Further, this information can most likely be extrapolated to all patients with P. aeruginosa septicemia.

Author Video: An author video summary of this article is available.

Coherent Feedforward Regulation of Gene Expression by Caulobacter σT and GsrN during Hyperosmotic Stress

GsrN is a conserved small RNA that is under transcriptional control of the general stress sigma factor, σ^T, and that functions as a posttranscriptional regulator of Caulobacter crescentus survival under multiple stress conditions. We have defined features of GsrN structure that determine survival under hyperosmotic stress, and we have applied transcriptomic and proteomic methods to identify regulatory targets of GsrN under hyperosmotic conditions. The 5' end of GsrN, which includes a conserved cytosine-rich stem-loop structure, is necessary for cell survival after osmotic upshock. GsrN both activates and represses gene expression in this stress condition. Expression of an uncharacterized open reading frame predicted to encode a glycine zipper protein, osrP, is strongly activated by GsrN. Our data support a model in which GsrN physically interacts with osrP mRNA through its 5' C-rich stem-loop to enhance OsrP protein expression. We conclude that sigT, gsrN, and osrP form a coherent feedforward loop in which σ^T activates gsrN and osrP transcription during stress and GsrN activates OsrP protein expression at the posttranscriptional level. This study delineates transcriptional and posttranscriptional layers of Caulobacter gene expression control during hyperosmotic stress, uncovers a new regulatory target of GsrN, and defines a coherent feedforward motif in the Caulobacter general stress response (GSR) regulatory network.IMPORTANCE Bacteria inhabit diverse niches and must adapt their physiology to constant environmental fluctuations. A major response to environmental perturbation is to change gene expression. Caulobacter and other alphaproteobacteria initiate a complex gene expression program known as the general stress response (GSR) under conditions including oxidative stress, osmotic stress, and nutrient limitation. The GSR enables cell survival in these environments. Understanding how bacteria survive stress requires that we dissect gene expression responses, such as the GSR, at the molecular level. This study is significant, as it defines transcriptional and posttranscriptional layers of gene expression regulation in response to hyperosmotic stress. We further provide evidence that a coherent feedforward motif influences the system properties of the Caulobacter GSR pathway.

cAMP and Vfr Control Exolysin Expression and Cytotoxicity of Pseudomonas aeruginosa Taxonomic Outliers

The two-partner secretion system ExlBA, expressed by strains of Pseudomonas aeruginosa belonging to the PA7 group, induces hemorrhage in lungs due to disruption of host cellular membranes. Here we demonstrate that the exlBA genes are controlled by a pathway consisting of cAMP and the virulence factor regulator (Vfr). Upon interaction with cAMP, Vfr binds directly to the exlBA promoter with high affinity (equilibrium binding constant [K_eq] of ≈2.5 nM). The exlB and exlA expression was diminished in the Vfr-negative mutant and upregulated with increased intracellular cAMP levels. The Vfr binding sequence in the exlBA promoter was mutated in situ, resulting in reduced cytotoxicity of the mutant, showing that Vfr is required for the exlBA expression during intoxication of epithelial cells. Vfr also regulates function of type 4 pili previously shown to facilitate ExlA activity on epithelial cells, which indicates that the cAMP/Vfr pathway coordinates these two factors needed for full cytotoxicity. As in most P. aeruginosa strains, the adenylate cyclase CyaB is the main provider of cAMP for Vfr regulation during both in vitro growth and eukaryotic cell infection. We discovered that the absence of functional Vfr in the reference strain PA7 is caused by a frameshift in the gene and accounts for its reduced cytotoxicity, revealing the conservation of ExlBA control by the CyaB-cAMP/Vfr pathway in P. aeruginosa taxonomic outliers.IMPORTANCE The human opportunistic pathogen Pseudomonas aeruginosa provokes severe acute and chronic human infections associated with defined sets of virulence factors. The main virulence determinant of P. aeruginosa taxonomic outliers is exolysin, a membrane-disrupting pore-forming toxin belonging to the two-partner secretion system ExlBA. In this work, we demonstrate that the conserved CyaB-cAMP/Vfr pathway controls cytotoxicity of outlier clinical strains through direct transcriptional activation of the exlBA operon. Therefore, despite the fact that the type III secretion system and exolysin are mutually exclusive in classical and outlier strains, respectively, these two major virulence determinants share similarities in their mechanisms of regulation.

Antisense transcription in Pseudomonas aeruginosa

A large number of antisense transcripts have been detected in diverse microbial genomes and considerable effort has been devoted to elucidating the functional role of antisense transcription. In this study, we reanalysed extensive RNA sequencing data from the opportunistic pathogen Pseudomonas aeruginosa and found that the majority of genes have a propensity for antisense transcription. Although antisense transcripts were found in more than 80 % of the genes of the P. aeruginosa genome, the majority of sequencing reads were mapping sense and only a minority (<2 %) were mapping antisense to genes. Similarly to the sense expression levels, the antisense expression levels varied under different environmental conditions, with the sense and antisense expression levels often being inversely regulated and modulated by the activity of alternative sigma factors. Environment-modulated antisense transcription showed a bias towards being antisense to genes within regions of genomic plasticity and to those encoding small regulatory RNAs. In the future, the validation and functional characterization of antisense transcripts, and novel transcripts that are antisense to small regulatory RNAs in particular, have the potential to contribute to our understanding of the various levels of transcriptional regulation and its dynamics in the bacterial pathogen P. aeruginosa.

Identification of sRNA mediated responses to nutrient depletion in Burkholderia pseudomallei

The Burkholderia genus includes many species that are known to survive in diverse environmental conditions including low nutrient environments. One species, Burkholderia pseudomallei is a versatile pathogen that can survive in a wide range of hosts and environmental conditions. In this study, we investigated how a nutrient depleted growth environment evokes sRNA mediated responses by B. pseudomallei. Computationally predicted B. pseudomallei D286 sRNAs were mapped to RNA-sequencing data for cultures grown under two conditions: (1) BHIB as a nutrient rich media reference environment and (2) M9 media as a nutrient depleted stress environment. The sRNAs were further selected to identify potentially cis-encoded systems by investigating their possible interactions with their flanking genes. The mappings of predicted sRNA genes and interactions analysis to their flanking genes identified 12 sRNA candidates that may possibly have cis-acting regulatory roles that are associated to a nutrient depleted growth environment. Our approach can be used for identifying novel sRNA genes and their possible role as cis-mediated regulatory systems.

Global Transcriptional Responses to Osmotic, Oxidative, and Imipenem Stress Conditions in Pseudomonas putida

Bacteria cope with and adapt to stress by modulating gene expression in response to specific environmental cues. In this study, the transcriptional response of Pseudomonas putida KT2440 to osmotic, oxidative, and imipenem stress conditions at two time points was investigated via identification of differentially expressed mRNAs and small RNAs (sRNAs). A total of 440 sRNA transcripts were detected, of which 10% correspond to previously annotated sRNAs, 40% to novel intergenic transcripts, and 50% to novel transcripts antisense to annotated genes. Each stress elicits a unique response as far as the extent and dynamics of the transcriptional changes. Nearly 200 protein-encoding genes exhibited significant changes in all stress types, implicating their participation in a general stress response. Almost half of the sRNA transcripts were differentially expressed under at least one condition, suggesting possible functional roles in the cellular response to stress conditions. The data show a larger fraction of differentially expressed sRNAs than of mRNAs with >5-fold expression changes. The work provides detailed insights into the mechanisms through which P. putida responds to different stress conditions and increases understanding of bacterial adaptation in natural and industrial settings.IMPORTANCE This study maps the complete transcriptional response of P. putida KT2440 to osmotic, oxidative, and imipenem stress conditions at short and long exposure times. Over 400 sRNA transcripts, consisting of both intergenic and antisense transcripts, were detected, increasing the number of identified sRNA transcripts in the strain by a factor of 10. Unique responses to each type of stress are documented, including both the extent and dynamics of the gene expression changes. The work adds rich detail to previous knowledge of stress response mechanisms due to the depth of the RNA sequencing data. Almost half of the sRNAs exhibit significant expression changes under at least one condition, suggesting their involvement in adaptation to stress conditions and identifying interesting candidates for further functional characterization.

Systems-level understanding of ethanol-induced stresses and adaptation in E. coli

Understanding ethanol-induced stresses and responses in biofuel-producing bacteria at systems level has significant implications in engineering more efficient biofuel producers. We present a computational study of transcriptomic and genomic data of both ethanol-stressed and ethanol-adapted E. coli cells with computationally predicated ethanol-binding proteins and experimentally identified ethanol tolerance genes. Our analysis suggests: (1) ethanol damages cell wall and membrane integrity, causing increased stresses, particularly reactive oxygen species, which damages DNA and reduces the O2 level; (2) decreased cross-membrane proton gradient from membrane damage, coupled with hypoxia, leads to reduced ATP production by aerobic respiration, driving cells to rely more on fatty acid oxidation, anaerobic respiration and fermentation for ATP production; (3) the reduced ATP generation results in substantially decreased synthesis of macromolecules; (4) ethanol can directly bind 213 proteins including transcription factors, altering their functions; (5) all these changes together induce multiple stress responses, reduced biosynthesis, cell viability and growth; and (6) ethanol-adapted E. coli cells restore the majority of these reduced activities through selection of specific genomic mutations and alteration of stress responses, ultimately restoring normal ATP production, macromolecule biosynthesis, and growth. These new insights into the energy and mass balance will inform design of more ethanol-tolerant strains.

Comparative evaluation of rRNA depletion procedures for the improved analysis of bacterial biofilm and mixed pathogen culture transcriptomes

Global transcriptomic analysis via RNA-seq is often hampered by the high abundance of ribosomal (r)RNA in bacterial cells. To remove rRNA and enrich coding sequences, subtractive hybridization procedures have become the approach of choice prior to RNA-seq, with their efficiency varying in a manner dependent on sample type and composition. Yet, despite an increasing number of RNA-seq studies, comparative evaluation of bacterial rRNA depletion methods has remained limited. Moreover, no such study has utilized RNA derived from bacterial biofilms, which have potentially higher rRNA:mRNA ratios and higher rRNA carryover during RNA-seq analysis. Presently, we evaluated the efficiency of three subtractive hybridization-based kits in depleting rRNA from samples derived from biofilm, as well as planktonic cells of the opportunistic human pathogen Pseudomonas aeruginosa. Our results indicated different rRNA removal efficiency for the three procedures, with the Ribo-Zero kit yielding the highest degree of rRNA depletion, which translated into enhanced enrichment of non-rRNA transcripts and increased depth of RNA-seq coverage. The results indicated that, in addition to improving RNA-seq sensitivity, efficient rRNA removal enhanced detection of low abundance transcripts via qPCR. Finally, we demonstrate that the Ribo-Zero kit also exhibited the highest efficiency when P. aeruginosa/Staphylococcus aureus co-culture RNA samples were tested.

Understanding butanol tolerance and assimilation in Pseudomonas putida BIRD-1: an integrated omics approach

Pseudomonas putida  BIRD‐1 has the potential to be used for the industrial production of butanol due to its solvent tolerance and ability to metabolize low‐cost compounds. However, the strain has two major limitations: it assimilates butanol as sole carbon source and butanol concentrations above 1% (v/v) are toxic. With the aim of facilitating BIRD‐1 strain design for industrial use, a genome‐wide mini‐Tn5 transposon mutant library was screened for clones exhibiting increased butanol sensitivity or deficiency in butanol assimilation. Twenty‐one mutants were selected that were affected in one or both of the processes. These mutants exhibited insertions in various genes, including those involved in the TCA cycle, fatty acid metabolism, transcription, cofactor synthesis and membrane integrity. An omics‐based analysis revealed key genes involved in the butanol response. Transcriptomic and proteomic studies were carried out to compare short and long‐term tolerance and assimilation traits. Pseudomonas putida initiates various butanol assimilation pathways via alcohol and aldehyde dehydrogenases that channel the compound to central metabolism through the glyoxylate shunt pathway. Accordingly, isocitrate lyase – a key enzyme of the pathway – was the most abundant protein when butanol was used as the sole carbon source. Upregulation of two genes encoding proteins PPUBIRD1_2240 and PPUBIRD1_2241 (acyl‐CoA dehydrogenase and acyl‐CoA synthetase respectively) linked butanol assimilation with acyl‐CoA metabolism. Butanol tolerance was found to be primarily linked to classic solvent defense mechanisms, such as efflux pumps, membrane modifications and control of redox state. Our results also highlight the intensive energy requirements for butanol production and tolerance; thus, enhancing TCA cycle operation may represent a promising strategy for enhanced butanol production.

A computational strategy for the search of regulatory small RNAs in Actinobacillus pleuropneumoniae

Bacterial regulatory small RNAs (sRNAs) play important roles in gene regulation and are frequently connected to the expression of virulence factors in diverse bacteria. Only a few sRNAs have been described for Pasteurellaceae pathogens and no in-depth analysis of sRNAs has been described for Actinobacillus pleuropneumoniae, the causative agent of porcine pleuropneumonia, responsible for considerable losses in the swine industry. To search for sRNAs in A. pleuropneumoniae, we developed a strategy for the computational analysis of the bacterial genome by using four algorithms with different approaches, followed by experimental validation. The coding strand and expression of 17 out of 23 RNA candidates were confirmed by Northern blotting, RT-PCR, and RNA sequencing. Among them, two are likely riboswitches, three are housekeeping regulatory RNAs, two are the widely studied GcvB and 6S sRNAs, and 10 are putative novel trans-acting sRNAs, never before described for any bacteria. The latter group has several potential mRNA targets, many of which are involved with virulence, stress resistance, or metabolism, and connect the sRNAs in a complex gene regulatory network. The sRNAs identified are well conserved among the Pasteurellaceae that are evolutionarily closer to A. pleuropneumoniae and/or share the same host. Our results show that the combination of newly developed computational programs can be successfully utilized for the discovery of novel sRNAs and indicate an intricate system of gene regulation through sRNAs in A. pleuropneumoniae and in other Pasteurellaceae, thus providing clues for novel aspects of virulence that will be explored in further studies.

A Transcriptome Map of Actinobacillus pleuropneumoniae at Single-Nucleotide Resolution Using Deep RNA-Seq

Actinobacillus pleuropneumoniae is the pathogen of porcine contagious pleuropneumoniae, a highly contagious respiratory disease of swine. Although the genome of A. pleuropneumoniae was sequenced several years ago, limited information is available on the genome-wide transcriptional analysis to accurately annotate the gene structures and regulatory elements. High-throughput RNA sequencing (RNA-seq) has been applied to study the transcriptional landscape of bacteria, which can efficiently and accurately identify gene expression regions and unknown transcriptional units, especially small non-coding RNAs (sRNAs), UTRs and regulatory regions. The aim of this study is to comprehensively analyze the transcriptome of A. pleuropneumoniae by RNA-seq in order to improve the existing genome annotation and promote our understanding of A. pleuropneumoniae gene structures and RNA-based regulation. In this study, we utilized RNA-seq to construct a single nucleotide resolution transcriptome map of A. pleuropneumoniae. More than 3.8 million high-quality reads (average length ~90 bp) from a cDNA library were generated and aligned to the reference genome. We identified 32 open reading frames encoding novel proteins that were mis-annotated in the previous genome annotations. The start sites for 35 genes based on the current genome annotation were corrected. Furthermore, 51 sRNAs in the A. pleuropneumoniae genome were discovered, of which 40 sRNAs were never reported in previous studies. The transcriptome map also enabled visualization of 5'- and 3'-UTR regions, in which contained 11 sRNAs. In addition, 351 operons covering 1230 genes throughout the whole genome were identified. The RNA-Seq based transcriptome map validated annotated genes and corrected annotations of open reading frames in the genome, and led to the identification of many functional elements (e.g. regions encoding novel proteins, non-coding sRNAs and operon structures). The transcriptional units described in this study provide a foundation for future studies concerning the gene functions and the transcriptional regulatory architectures of this pathogen.

Bacterial antisense RNAs are mainly the product of transcriptional noise

cis-Encoded antisense RNAs (asRNAs) are widespread along bacterial transcriptomes. However, the role of most of these RNAs remains unknown, and there is an ongoing discussion as to what extent these transcripts are the result of transcriptional noise. We show, by comparative transcriptomics of 20 bacterial species and one chloroplast, that the number of asRNAs is exponentially dependent on the genomic AT content and that expression of asRNA at low levels exerts little impact in terms of energy consumption. A transcription model simulating mRNA and asRNA production indicates that the asRNA regulatory effect is only observed above certain expression thresholds, substantially higher than physiological transcript levels. These predictions were verified experimentally by overexpressing nine different asRNAs in Mycoplasma pneumoniae. Our results suggest that most of the antisense transcripts found in bacteria are the consequence of transcriptional noise, arising at spurious promoters throughout the genome.

Identification and functional analysis of a bacteriocin, pyocin S6, with ribonuclease activity from a Pseudomonas aeruginosa cystic fibrosis clinical isolate

S‐type pyocins are bacteriocins produced by Pseudomonas aeruginosa isolates to antagonize or kill other strains of the same species. They have a modular organization comprising a receptor‐binding domain recognizing a surface constituent of the target bacterium, a domain for translocation through the periplasm, and a killing or toxic domain with DNase, tRNase, or pore‐forming activity. Pyocins S2, S3, S4, and S5 recognize TonB‐dependent ferri‐siderophore receptors in the outer membrane. We here describe a new nuclease bacteriocin, pyocin S6, encoded in the genome of a P. aeruginosa cystic fibrosis (CF) clinical isolate, CF_PA39. Similarly to pyocins S1 and S2, the S6 toxin–immunity gene tandem was recruited to the genomic region encoding exotoxin A. The pyocin S6 receptor‐binding and translocation domains are identical to those of pyocin S1, whereas the killing domain is similar to the 16S ribonuclease domain of Escherichia coli colicin E3. The cytotoxic activity was abolished in pyocin S6 forms with a mutation in the colicin E3‐equivalent catalytic motif. The CF_PA39 S6 immunity gene displays a higher expression level than the gene encoding the killing protein, the latter being only detected when bacteria are grown under iron‐limiting conditions. In the S1‐pyocinogenic strain P. aeruginosa ATCC 25324 and pyocin S2 producer P. aeruginosa PAO1, a remnant of the pyocin S6 killing domain and an intact S6‐type immunity gene are located downstream of their respective pyocin operons. Strain PAO1 is insensitive for pyocin S6, and its S6‐type immunity gene provides protection against pyocin S6 activity. Purified pyocin S6 inhibits one‐fifth of 110 P. aeruginosa CF clinical isolates tested, showing clearer inhibition zones when the target cells are grown under iron limitation. In this panel, about half of the CF clinical isolates were found to host the S6 genes. The pyocin S6 locus is also present in the genome of some non‐CF clinical isolates.

Antibacterial activity of silver nanoparticles target sara through srna-teg49, a key mediator of hfq, in staphylococcus aureus

Attributed to its antimicrobial effect, Silver nanoparticles (AgNPs) is widely used in various fields, such as biomedicine, textiles, health care products and food, etc. However, the antibacterial mechanism of AgNPs in staphylococcus aureus (S. aureus) by regulating sRNA expression remains largely unknown.

OBJECTIVES

This study was performed to investigate the involvement of the antibacterial mechanism of AgNPs through sRNA-TEG49, a key mediator of Hfq, in S. aureus.

METHODS

Through the antimicrobial tests of AgNPs, its antibacterial laps and minimum inhibitory concentration was measured. A hierarchical cluster analysis of the differentially expressed sRNA in S. aureus was performed to investigate the relationship between AgNPs and sRNA. Expression of genes was analyzed by real-time PCR.

RESULTS

In the present study we found that at the concentrations higher than 1 mg/L, AgNPs could completely restrain bacteria growth, and the antibacterial activity of AgNPs apparently declined at the concentrations lower than 1 mg/L. S. aureus exposure to AgNPs, the expression of sRNA-TEG49, Hfq and sarA was significantly up-regulated in wild-type S. aureus. Moreover, Hfq loss-of-function inhibited the expression of sRNA-TEG49 in mutant-type S. aureus. Furthermore, sRNA-TEG49 loss-of-function associated with down-regulation the expression of sarA in mutant-type S. aureus.

CONCLUSIONS

It was reasonable that Hfq regulated a distinct underlying molecular and antibacterial mechanism of AgNPs by forming a positive feedback loop with sRNA-TEG49. These observations suggested that Hfq plays an important role in the antibacterial mechanism of AgNPs by regulating sRNA-TEG49 expression, via its target sarA.

Diversity of small RNAs expressed in Pseudomonas species

RNA sequencing (RNA-seq) has revealed several hundreds of previously undetected small RNAs (sRNAs) in all bacterial species investigated, including strains of Pseudomonas aeruginosa, Pseudomonas putida and Pseudomonas syringae. Nonetheless, only little is known about the extent of conservation of expressed sRNAs across strains and species. In this study, we have used RNA-seq to identify sRNAs in P. putida DOT-T1E and Pseudomonas extremaustralis 14-3b. This is the first strain of P. extremaustralis and the second strain of P. putida to have their transcriptomes analysed for sRNAs, and we identify the presence of around 150 novel sRNAs in each strain. Furthermore, we provide a comparison based on sequence conservation of all the sRNAs detected by RNA-seq in the Pseudomonas species investigated so far. Our results show that the extent of sRNA conservation across different species is very limited. In addition, when comparing the sRNAs expressed in different strains of the same species, we observe that numerous sRNAs exhibit a strain-specific expression pattern. These results support the idea that the evolution of most bacterial sRNAs is rapid, which limits the extent of both interspecies and intraspecies conservation.