Identification and characterization of small RNAs in Yersinia pestis

Yersiniomics, a Multi-Omics Interactive Database for Yersinia Species

The genus Yersinia includes a large variety of nonpathogenic and life-threatening pathogenic bacteria, which cause a broad spectrum of diseases in humans and animals, such as plague, enteritis, Far East scarlet-like fever (FESLF), and enteric redmouth disease. Like most clinically relevant microorganisms, Yersinia spp. are currently subjected to intense multi-omics investigations whose numbers have increased extensively in recent years, generating massive amounts of data useful for diagnostic and therapeutic developments. The lack of a simple and centralized way to exploit these data led us to design Yersiniomics, a web-based platform allowing straightforward analysis of Yersinia omics data. Yersiniomics contains a curated multi-omics database at its core, gathering 200 genomic, 317 transcriptomic, and 62 proteomic data sets for Yersinia species. It integrates genomic, transcriptomic, and proteomic browsers, a genome viewer, and a heatmap viewer to navigate within genomes and experimental conditions. For streamlined access to structural and functional properties, it directly links each gene to GenBank, the Kyoto Encyclopedia of Genes and Genomes (KEGG), UniProt, InterPro, IntAct, and the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and each experiment to Gene Expression Omnibus (GEO), the European Nucleotide Archive (ENA), or the Proteomics Identifications Database (PRIDE). Yersiniomics provides a powerful tool for microbiologists to assist with investigations ranging from specific gene studies to systems biology studies. IMPORTANCE The expanding genus Yersinia is composed of multiple nonpathogenic species and a few pathogenic species, including the deadly etiologic agent of plague, Yersinia pestis. In 2 decades, the number of genomic, transcriptomic, and proteomic studies on Yersinia grew massively, delivering a wealth of data. We developed Yersiniomics, an interactive web-based platform, to centralize and analyze omics data sets on Yersinia species. The platform allows user-friendly navigation between genomic data, expression data, and experimental conditions. Yersiniomics will be a valuable tool to microbiologists.

CsrA Enhances Cyclic-di-GMP Biosynthesis and Yersinia pestis Biofilm Blockage of the Flea Foregut by Alleviating Hfq-Dependent Repression of the hmsT mRNA

Plague-causing Yersinia pestis is transmitted through regurgitation when it forms a biofilm-mediated blockage in the foregut of its flea vector. This biofilm is composed of an extracellular polysaccharide substance (EPS) produced when cyclic-di-GMP (c-di-GMP) levels are elevated. The Y. pestis diguanylate cyclase enzymes HmsD and HmsT synthesize c-di-GMP. HmsD is required for biofilm blockage formation but contributes minimally to in vitro biofilms. HmsT, however, is necessary for in vitro biofilms and contributes to intermediate rates of biofilm blockage. C-di-GMP synthesis is regulated at the transcriptional and posttranscriptional levels. In this, the global RNA chaperone, Hfq, posttranscriptionally represses hmsT mRNA translation. How c-di-GMP levels and biofilm blockage formation is modulated by nutritional stimuli encountered in the flea gut is unknown. Here, the RNA-binding regulator protein CsrA, which controls c-di-GMP-mediated biofilm formation and central carbon metabolism responses in many Gammaproteobacteria, was assessed for its role in Y. pestis biofilm formation. We determined that CsrA was required for markedly greater c-di-GMP and EPS levels when Y. pestis was cultivated on alternative sugars implicated in flea biofilm blockage metabolism. Our assays, composed of mobility shifts, quantification of mRNA translation, stability, and abundance, and epistasis analyses of a csrA hfq double mutant strain substantiated that CsrA represses hfq mRNA translation, thereby alleviating Hfq-dependent repression of hmsT mRNA translation. Additionally, a csrA mutant exhibited intermediately reduced biofilm blockage rates, resembling an hmsT mutant. Hence, we reveal CsrA-mediated control of c-di-GMP synthesis in Y. pestis as a tiered, posttranscriptional regulatory process that enhances biofilm blockage-mediated transmission from fleas.

Participation of two sRNA RyhB homologs from the fish pathogen Yersinia ruckeri in bacterial physiology

Small noncoding RNAs (sRNAs) are important regulators of gene expression and physiology in bacteria. RyhB is an iron-responsive sRNA well characterized in Escherichia coli and conserved in other Enterobacteriaceae. In this study, we identified and characterized two RyhB homologs (named RyhB-1 and RyhB-2) in the fish pathogen Yersinia ruckeri. We found that, as in other Enterobacteriaceae, both RyhB-1 and RyhB-2 are induced under iron starvation, repressed by the Fur regulator, and depend on Hfq for stability. Despite these similarities in expression, the mutant strains of Y. ruckeri lacking RyhB-1 (ΔryhB-1) or RyhB-2 (ΔryhB-2) exhibited differential phenotypes. In comparison with the wild type, the ΔryhB-1 strain showed a hypermotile phenotype, reduced biofilm formation, increased replication rate, faster growth, and increased ATP levels in bacterial cultures. By contrast, in salmon cell cultures, the ΔryhB-1 strain exhibited an increased survival. On the other hand, the ΔryhB-2 strain was non-motile and showed augmented biofilm formation as compared to the wild type. The expression of a subset of RyhB conserved targets, selected from different bacterial species, was analyzed by quantitative RT-PCR in wild type, ΔryhB-1, ΔryhB-2, and ΔryhB-1 ΔryhB-2 strains cultured in iron-depleted media. RyhB-1 negatively affected the expression of most analyzed genes (sodB, acnA, sdhC, bfr, fliF, among others), whose functions are related to metabolism and motility, involving iron-containing proteins. Among the genes analyzed, only sdhC and bfr appeared as targets for RyhB-2. Taken together, these results indicate that Y. ruckeri RyhB homologs participate in the modulation of the bacterial physiology with non-redundant roles.

Small RNA Mcr11 requires the transcription factor AbmR for stable expression and regulates genes involved in the central metabolism of Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb), the etiologic agent of tuberculosis, must adapt to host-associated environments during infection by modulating gene expression. Small regulatory RNAs (sRNAs) are key regulators of bacterial gene expression, but their roles in Mtb are not well understood. Here, we address the expression and function of the Mtb sRNA Mcr11, which is associated with slow bacterial growth and chronic infections in mice. We found that stable expression of Mcr11 requires multiple factors specific to TB-complex bacteria, including the AbmR transcription factor. Bioinformatic analyses used to predict regulatory targets of Mcr11 identified 7-11 nucleotide regions with potential for direct base-pairing with Mcr11 immediately upstream of Rv3282, fadA3, and lipB. mcr11-dependent regulation of these genes was demonstrated using qRT-PCR and found to be responsive to the presence of fatty acids. Mutation of the putative Mcr11 base-pairing site upstream of lipB in a promoter reporter strain resulted in significant de-repression of lipB expression, similar to that observed in mcr11-deleted Mtb. These studies establish Mcr11's roles in regulating growth and central metabolism in Mtb. Our finding that multiple TB-complex-specific factors are required for production of stable Mcr11 also emphasizes the need to better understand mechanisms of sRNA expression and stability in TB.

Hfq Globally Binds and Destabilizes sRNAs and mRNAs in Yersinia pestis

Hfq is a ubiquitous Sm-like RNA-binding protein in bacteria involved in physiological fitness and pathogenesis, while its in vivo binding nature remains elusive. Here we reported genome-wide Hfq-bound RNAs in Yersinia pestis, a causative agent of plague, by using cross-linking immunoprecipitation coupled with deep sequencing (CLIP-seq) approach. We show that the Hfq binding density is enriched in more than 80% mRNAs of Y. pestis and that Hfq also globally binds noncoding small RNAs (sRNAs) encoded by the intergenic, antisense, and 3' regions of mRNAs. An Hfq U-rich stretch is highly enriched in sRNAs, while motifs partially complementary to AGAAUAA and GGGGAUUA are enriched in both mRNAs and sRNAs. Hfq-binding motifs are enriched at both terminal sites and in the gene body of mRNAs. Surprisingly, a large fraction of the sRNA and mRNA regions bound by Hfq and those downstream are destabilized, likely via a 5'P-activated RNase E degradation pathway, which is consistent with a model in which Hfq facilitates sRNA-mRNA base pairing and the coupled degradation in Y. pestis These results together have presented a high-quality Hfq-RNA interaction map in Y. pestis, which should be important for further deciphering the regulatory role of Hfq-sRNAs in Y. pestis IMPORTANCE Discovered in 1968 as an Escherichia coli host factor that was essential for replication of the bacteriophage Qβ, the Hfq protein is a ubiquitous and highly abundant RNA-binding protein in many bacteria. With the assistance of Hfq, small RNAs in bacteria play important roles in regulating the stability and translation of mRNAs by base pairing. In this study, we want to elucidate the Hfq-assisted sRNA-mRNA regulation in Yersinia pestis A global map of Hfq interaction sites in Y. pestis was obtained by sequencing cDNAs converted from the Hfq-bound RNA fragments using UV cross-linking coupled immunoprecipitation technology. We demonstrate that Hfq could bind to hundreds of sRNAs and the majority of mRNAs in Y. pestis The enriched binding motifs in sRNAs and mRNAs are complementary to each other, suggesting a general base-pairing mechanism for sRNA-mRNA interaction. The Hfq-bound sRNA and mRNA regions were both destabilized. The results suggest that Hfq binding facilitates sRNA-mRNA base pairing and coordinates their degradation, which might enable Hfq to surveil the homeostasis of most mRNAs in bacteria.

Comparative Transcriptomic Profiling of Yersinia enterocolitica O:3 and O:8 Reveals Major Expression Differences of Fitness- and Virulence-Relevant Genes Indicating Ecological Separation

Yersinia enterocolitica is a major diarrheal pathogen and is associated with a large range of gut-associated diseases. Members of this species have evolved into different phylogroups with genotypic variations. We performed the first characterization of the Y. enterocolitica transcriptional landscape and tracked the consequences of the genomic variations between two different pathogenic phylogroups by comparing their RNA repertoire, promoter usage, and expression profiles under four different virulence-relevant conditions. Our analysis revealed major differences in the transcriptional outputs of the closely related strains, pointing to an ecological separation in which one is more adapted to an environmental lifestyle and the other to a mostly mammal-associated lifestyle. Moreover, a variety of pathoadaptive alterations, including alterations in acid resistance genes, colonization factors, and toxins, were identified which affect virulence and host specificity. This illustrates that comparative transcriptomics is an excellent approach to discover differences in the functional output from closely related genomes affecting niche adaptation and virulence, which cannot be directly inferred from DNA sequences.

Yersinia enterocolitica is a zoonotic pathogen and an important cause of bacterial gastrointestinal infections in humans. Large-scale population genomic analyses revealed genetic and phenotypic diversity of this bacterial species, but little is known about the differences in the transcriptome organization, small RNA (sRNA) repertoire, and transcriptional output. Here, we present the first comparative high-resolution transcriptome analysis of Y. enterocolitica strains representing highly pathogenic phylogroup 2 (serotype O:8) and moderately pathogenic phylogroup 3 (serotype O:3) grown under four infection-relevant conditions. Our transcriptome sequencing (RNA-seq) approach revealed 1,299 and 1,076 transcriptional start sites and identified strain-specific sRNAs that could contribute to differential regulation among the phylogroups. Comparative transcriptomics further uncovered major gene expression differences, in particular, in the temperature-responsive regulon. Multiple virulence-relevant genes are differentially regulated between the two strains, supporting an ecological separation of phylogroups with certain niche-adapted properties. Strong upregulation of the ystA enterotoxin gene in combination with constitutive high expression of cell invasion factor InvA further showed that the toxicity of recent outbreak O:3 strains has increased. Overall, our report provides new insights into the specific transcriptome organization of phylogroups 2 and 3 and reveals gene expression differences contributing to the substantial phenotypic differences that exist between the lineages.

IMPORTANCEYersinia enterocolitica is a major diarrheal pathogen and is associated with a large range of gut-associated diseases. Members of this species have evolved into different phylogroups with genotypic variations. We performed the first characterization of the Y. enterocolitica transcriptional landscape and tracked the consequences of the genomic variations between two different pathogenic phylogroups by comparing their RNA repertoire, promoter usage, and expression profiles under four different virulence-relevant conditions. Our analysis revealed major differences in the transcriptional outputs of the closely related strains, pointing to an ecological separation in which one is more adapted to an environmental lifestyle and the other to a mostly mammal-associated lifestyle. Moreover, a variety of pathoadaptive alterations, including alterations in acid resistance genes, colonization factors, and toxins, were identified which affect virulence and host specificity. This illustrates that comparative transcriptomics is an excellent approach to discover differences in the functional output from closely related genomes affecting niche adaptation and virulence, which cannot be directly inferred from DNA sequences.

Enteropathogens: Tuning Their Gene Expression for Hassle-Free Survival

Enteropathogenic bacteria have been the cause of the majority of foodborne illnesses. Much of the research has been focused on elucidating the mechanisms by which these pathogens evade the host immune system. One of the ways in which they achieve the successful establishment of a niche in the gut microenvironment and survive is by a chain of elegantly regulated gene expression patterns. Studies have shown that this process is very elaborate and is also regulated by several factors. Pathogens like, enteropathogenic Escherichia coli (EPEC), Salmonella Typhimurium, Shigellaflexneri, Yersinia sp. have been seen to employ various regulated gene expression strategies. These include toxin-antitoxin systems, quorum sensing systems, expression controlled by nucleoid-associated proteins (NAPs), several regulons and operons specific to these pathogens. In the following review, we have tried to discuss the common gene regulatory systems of enteropathogenic bacteria as well as pathogen-specific regulatory mechanisms.

Discovering RNA-Based Regulatory Systems for Yersinia Virulence

The genus Yersinia includes three human pathogenic species, Yersinia pestis, the causative agent of the bubonic and pneumonic plague, and enteric pathogens Y. enterocolitica and Y. pseudotuberculosis that cause a number of gut-associated diseases. Over the past years a large repertoire of RNA-based regulatory systems has been discovered in these pathogens using different RNA-seq based approaches. Among them are several conserved or species-specific RNA-binding proteins, regulatory and sensory RNAs as well as various RNA-degrading enzymes. Many of them were shown to control the expression of important virulence-relevant factors and have a very strong impact on Yersinia virulence. The precise targets, the molecular mechanism and their role for Yersinia pathogenicity is only known for a small subset of identified genus- or species-specific RNA-based control elements. However, the ongoing development of new RNA-seq based methods and data analysis methods to investigate the synthesis, composition, translation, decay, and modification of RNAs in the bacterial cell will help us to generate a more comprehensive view of Yersinia RNA biology in the near future.

ANNOgesic: a Swiss army knife for the RNA-seq based annotation of bacterial/archaeal genomes

To understand the gene regulation of an organism of interest, a comprehensive genome annotation is essential. While some features, such as coding sequences, can be computationally predicted with high accuracy based purely on the genomic sequence, others, such as promoter elements or noncoding RNAs, are harder to detect. RNA sequencing (RNA-seq) has proven to be an efficient method to identify these genomic features and to improve genome annotations. However, processing and integrating RNA-seq data in order to generate high-resolution annotations is challenging, time consuming, and requires numerous steps. We have constructed a powerful and modular tool called ANNOgesic that provides the required analyses and simplifies RNA-seq-based bacterial and archaeal genome annotation. It can integrate data from conventional RNA-seq and differential RNA-seq and predicts and annotates numerous features, including small noncoding RNAs, with high precision. The software is available under an open source license (ISCL) at https://pypi.org/project/ANNOgesic/.

Primary transcriptome analysis reveals importance of IS elements for the shaping of the transcriptional landscape of Bordetella pertussis

Bordetella pertussis is the causative agent of whooping cough, a respiratory disease still considered as a major public health threat and for which recent re-emergence has been observed. Constant reshuffling of Bordetella pertussis genome organization was observed during evolution. These rearrangements are essentially mediated by Insertion Sequences (IS), a mobile genetic elements present in more than 230 copies in the genome, which are supposed to be one of the driving forces enabling the pathogen to escape from vaccine-induced immunity. Here we use high-throughput sequencing approaches (RNA-seq and differential RNA-seq), to decipher Bordetella pertussis transcriptome characteristics and to evaluate the impact of IS elements on transcriptome architecture. Transcriptional organization was determined by identification of transcription start sites and revealed also a large variety of non-coding RNAs including sRNAs, leaderless mRNAs or long 3' and 5'UTR including seven riboswitches. Unusual topological organizations, such as overlapping 5'- or 3'-extremities between oppositely orientated mRNA were also unveiled. The pivotal role of IS elements in the transcriptome architecture and their effect on the transcription of neighboring genes was examined. This effect is mediated by the introduction of IS harbored promoters or by emergence of hybrid promoters. This study revealed that in addition to their impact on genome rearrangements, most of the IS also impact on the expression of their flanking genes. Furthermore, the transcripts produced by IS are strain-specific due to the strain to strain variation in IS copy number and genomic context.

The evolutionary impact of intragenic FliA promoters in proteobacteria

In Escherichia coli, one sigma factor recognizes the majority of promoters, and six 'alternative' sigma factors recognize specific subsets of promoters. The alternative sigma factor FliA (σ28 ) recognizes promoters upstream of many flagellar genes. We previously showed that most E. coli FliA binding sites are located inside genes. However, it was unclear whether these intragenic binding sites represent active promoters. Here, we construct and assay transcriptional promoter-lacZ fusions for all 52 putative FliA promoters previously identified by ChIP-seq. These experiments, coupled with integrative analysis of published genome-scale transcriptional datasets, strongly suggest that most intragenic FliA binding sites are active promoters that transcribe highly unstable RNAs. Additionally, we show that widespread intragenic FliA-dependent transcription may be a conserved phenomenon, but that specific promoters are not themselves conserved. We conclude that intragenic FliA-dependent promoters and the resulting RNAs are unlikely to have important regulatory functions. Nonetheless, one intragenic FliA promoter is broadly conserved and constrains evolution of the overlapping protein-coding gene. Thus, our data indicate that intragenic regulatory elements can influence bacterial protein evolution and suggest that the impact of intragenic regulatory sequences on genome evolution should be considered more broadly.

Genome-wide identification and characterization of small RNAs in Rhodobacter capsulatus and identification of small RNAs affected by loss of the response regulator CtrA

Small non-coding RNAs (sRNAs) are involved in the control of numerous cellular processes through various regulatory mechanisms, and in the past decade many studies have identified sRNAs in a multitude of bacterial species using RNA sequencing (RNA-seq). Here, we present the first genome-wide analysis of sRNA sequencing data in Rhodobacter capsulatus, a purple nonsulfur photosynthetic alphaproteobacterium. Using a recently developed bioinformatics approach, sRNA-Detect, we detected 422 putative sRNAs from R. capsulatus RNA-seq data. Based on their sequence similarity to sRNAs in a sRNA collection, consisting of published putative sRNAs from 23 additional bacterial species, and RNA databases, the sequences of 124 putative sRNAs were conserved in at least one other bacterial species; and, 19 putative sRNAs were assigned a predicted function. We bioinformatically characterized all putative sRNAs and applied machine learning approaches to calculate the probability of a nucleotide sequence to be a bona fide sRNA. The resulting quantitative model was able to correctly classify 95.2% of sequences in a validation set. We found that putative cis-targets for antisense and partially overlapping sRNAs were enriched with protein-coding genes involved in primary metabolic processes, photosynthesis, compound binding, and with genes forming part of macromolecular complexes. We performed differential expression analysis to compare the wild type strain to a mutant lacking the response regulator CtrA, an important regulator of gene expression in R. capsulatus, and identified 18 putative sRNAs with differing levels in the two strains. Finally, we validated the existence and expression patterns of four novel sRNAs by Northern blot analysis.

Several Hfq-dependent alterations in physiology of Yersinia enterocolitica O:3 are mediated by derepression of the transcriptional regulator RovM

In bacteria, the RNA chaperone Hfq enables pairing of small regulatory RNAs with their target mRNAs and therefore is a key player of post-transcriptional regulation network. As a global regulator, Hfq is engaged in the adaptation to external environment, regulation of metabolism and bacterial virulence. In this study we used RNA-sequencing and quantitative proteomics (LC-MS/MS) to elucidate the role of this chaperone in the physiology and virulence of Yersinia enterocolitica serotype O:3. This global approach revealed the profound impact of Hfq on gene and protein expression. Furthermore, the role of Hfq in the cell morphology, metabolism, cell wall integrity, resistance to external stresses and pathogenicity was evaluated. Importantly, our results revealed that several alterations typical for the hfq-negative phenotype were due to derepression of the transcriptional factor RovM. The overexpression of RovM caused by the loss of Hfq chaperone resulted in extended growth defect, alterations in the lipid A structure, motility and biofilm formation defects, as well as changes in mannitol utilization. Furthermore, in Y. enterocolitica RovM only in the presence of Hfq affected the abundance of RpoS. Finally, the impact of hfq and rovM mutations on the virulence was assessed in the mouse infection model.

Functional and Structural Analysis of a Highly-Expressed Yersinia pestis Small RNA following Infection of Cultured Macrophages

Non-coding small RNAs (sRNAs) are found in practically all bacterial genomes and play important roles in regulating gene expression to impact bacterial metabolism, growth, and virulence. We performed transcriptomics analysis to identify sRNAs that are differentially expressed in Yersinia pestis that invaded the human macrophage cell line THP-1, compared to pathogens that remained extracellular in the presence of host. Using ultra high-throughput sequencing, we identified 37 novel and 143 previously known sRNAs in Y. pestis. In particular, the sRNA Ysr170 was highly expressed in intracellular Yersinia and exhibited a log2 fold change ~3.6 higher levels compared to extracellular bacteria. We found that knock-down of Ysr170 expression attenuated infection efficiency in cell culture and growth rate in response to different stressors. In addition, we applied selective 2’-hydroxyl acylation analyzed by primer extension (SHAPE) analysis to determine the secondary structure of Ysr170 and observed structural changes resulting from interactions with the aminoglycoside antibiotic gentamycin and the RNA chaperone Hfq. Interestingly, gentamicin stabilized helix 4 of Ysr170, which structurally resembles the native gentamicin 16S ribosomal binding site. Finally, we modeled the tertiary structure of Ysr170 binding to gentamycin using RNA motif modeling. Integration of these experimental and structural methods can provide further insight into the design of small molecules that can inhibit function of sRNAs required for pathogen virulence.

RNA Regulators: Formidable Modulators of Yersinia Virulence

A large repertoire of RNA-based regulatory mechanisms, including a plethora of cis- and trans-acting noncoding RNAs (ncRNAs), sensory RNA elements, regulatory RNA-binding proteins, and RNA-degrading enzymes have been uncovered lately as key players in the regulation of metabolism, stress responses, and virulence of the genus Yersinia. Many of them are strictly controlled in response to fluctuating environmental conditions sensed during the course of the infection, and certain riboregulators have already been shown to be crucial for virulence. Some of them are highly conserved among the family Enterobacteriaceae, while others are genus-, species-, or strain-specific and could contribute to the difference in Yersinia pathogenicity. Importantly, the analysis of Yersinia riboregulators has not only uncovered crucial elements and regulatory mechanisms governing host-pathogen interactions, it also revealed exciting new venues for the design of novel anti-infectives.

Plasmid pPCP1-derived sRNA HmsA promotes biofilm formation of Yersinia pestis

Background

The ability of Yersinia pestis to form a biofilm is an important characteristic in flea transmission of this pathogen. Y. pestis laterally acquired two plasmids (pPCP1and pMT1) and the ability to form biofilms when it evolved from Yersinia pseudotuberculosis. Small regulatory RNAs (sRNAs) are thought to play a crucial role in the processes of biofilm formation and pathogenesis.

Results

A pPCP1-derived sRNA HmsA (also known as sR084) was found to contribute to the enhanced biofilm formation phenotype of Y. pestis. The concentration of c-di-GMP was significantly reduced upon deletion of the hmsA gene in Y. pestis. The abundance of mRNA transcripts determining exopolysaccharide production, crucial for biofilm formation, was measured by primer extension, RT-PCR and lacZ transcriptional fusion assays in the wild-type and hmsA mutant strains. HmsA positively regulated biofilm synthesis-associated genes (hmsHFRS, hmsT and hmsCDE), but had no regulatory effect on the biofilm degradation-associated gene hmsP. Interestingly, the recently identified biofilm activator sRNA, HmsB, was rapidly degraded in the hmsA deletion mutant. Two genes (rovM and rovA) functioning as biofilm regulators were also found to be regulated by HmsA, whose regulatory effects were consistent with the HmsA-mediated biofilm phenotype.

Conclusion

HmsA potentially functions as an activator of biofilm formation in Y. pestis, implying that sRNAs encoded on the laterally acquired plasmids might be involved in the chromosome-based regulatory networks implicated in Y. pestis-specific physiological processes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-016-0793-5) contains supplementary material, which is available to authorized users.

RNA-based mechanisms of virulence control in Enterobacteriaceae

Enteric pathogens of the family Enterobacteriaceae colonize various niches within animals and humans in which they compete with intestinal commensals and are attacked by the host immune system. To survive these hostile environments they possess complex, multilayer regulatory networks that coordinate the control of virulence factors, host-adapted metabolic functions and stress resistance. An important part of these intricate control networks are RNA-based control systems that enable the pathogen to fine-tune its responses. Recent next-generation sequencing approaches revealed a large repertoire of conserved and species-specific riboregulators, including numerous cis- and trans-acting non-coding RNAs, sensory RNA elements (RNA thermometers, riboswitches), regulatory RNA-binding proteins and RNA degrading enzymes which regulate colonization factors, toxins, host defense processes and virulence-relevant physiological and metabolic processes. All of which are important cues for pathogens to sense and respond to fluctuating conditions during the infection. This review covers infection-relevant riboregulators of E. coli, Salmonella, Shigella and Yersinia, highlights their versatile regulatory mechanisms, complex target regulons and functions, and discusses emerging topics and future challenges to fully understand and exploit RNA-based control to combat bacterial infections.

The RNA Chaperone Hfq Is Essential for Virulence and Modulates the Expression of Four Adhesins in Yersinia enterocolitica

In Enterobacteriaceae, the RNA chaperone Hfq mediates the interaction of small RNAs with target mRNAs, thereby modulating transcript stability and translation. This post-transcriptional control helps bacteria adapt quickly to changing environmental conditions. Our previous mutational analysis showed that Hfq is involved in metabolism and stress survival in the enteropathogen Yersinia enterocolitica. In this study we demonstrate that Hfq is essential for virulence in mice and influences production of surface pathogenicity factors, in particular lipopolysaccharide and adhesins mediating interaction with host tissue. Hfq inhibited the production of Ail, the Ail-like protein OmpX and the MyfA pilin post-transcriptionally. In contrast Hfq promoted production of two major autotransporter adhesins YadA and InvA. While protein secretion in vitro was not affected, hfq mutants exhibited decreased protein translocation by the type III secretion system into host cells, consistent with decreased production of YadA and InvA. The influence of Hfq on YadA resulted from a complex interplay of transcriptional, post-transcriptional and likely post-translational effects. Hfq regulated invA by modulating the expression of the transcriptional regulators rovA, phoP and ompR. Therefore, Hfq is a global coordinator of surface virulence determinants in Y. enterocolitica suggesting that it constitutes an attractive target for developing new antimicrobial strategies.

Reannotation of Yersinia pestis Strain 91001 Based on Omics Data

Yersinia pestis is among the most dangerous human pathogens, and systematic research of this pathogen is important in bacterial pathogenomics research. To fully interpret the biological functions, physiological characteristics, and pathogenesis of Y. pestis, a comprehensive annotation of its entire genome is necessary. The emergence of omics-based research has brought new opportunities to better annotate the genome of this pathogen. Here, the complete genome of Y. pestis strain 91001 was reannotated using genomics and proteogenomics data. One hundred and thirty-seven unreliable coding sequences were removed, and 41 homologous genes were relocated with their translational initiation sites, while the functions of seven pseudogenes and 392 hypothetical genes were revised. Moreover, annotations of noncoding RNAs, repeat sequences, and transposable elements have also been incorporated. The reannotated results are freely available at http://tody.bmi.ac.cn.

Genetic Regulation of Yersinia pestis

Y. pestis exhibits dramatically different traits of pathogenicity and transmission, albeit their close genetic relationship with its ancestor-Y. pseudotuberculosis, a self-limiting gastroenteric pathogen. Y. pestis is evolved into a deadly pathogen and transmitted to mammals and/or human beings by infected flea biting or directly contacting with the infected animals. Various kinds of environmental changes are implicated into its complex life cycle and pathogenesis. Dynamic regulation of gene expression is critical for environmental adaptation or survival, primarily reflected by genetic regulation mediated by transcriptional factors and small regulatory RNAs at the transcriptional and posttranscriptional level, respectively. The effects of genetic regulation have been shown to profoundly influence Y. pestis physiology and pathogenesis such as stress resistance, biofilm formation, intracellular survival, and replication. In this chapter, we mainly summarize the progresses on popular methods of genetic regulation and on regulatory patterns and consequences of many key transcriptional and posttranscriptional regulators, with a particular emphasis on how genetic regulation influences the biofilm and virulence of Y. pestis.

Yersinia pestis and Yersinia pseudotuberculosis infection: a regulatory RNA perspective

Yersinia pestis, responsible for causing fulminant plague, has evolved clonally from the enteric pathogen, Y. pseudotuberculosis, which in contrast, causes a relatively benign enteric illness. An ~97% nucleotide identity over 75% of their shared protein coding genes is maintained between these two pathogens, leaving much conjecture regarding the molecular determinants responsible for producing these vastly different disease etiologies, host preferences and transmission routes. One idea is that coordinated production of distinct factors required for host adaptation and virulence in response to specific environmental cues could contribute to the distinct pathogenicity distinguishing these two species. Small non-coding RNAs that direct posttranscriptional regulation have recently been identified as key molecules that may provide such timeous expression of appropriate disease enabling factors. Here the burgeoning field of small non-coding regulatory RNAs in Yersinia pathogenesis is reviewed from the viewpoint of adaptive colonization, virulence and divergent evolution of these pathogens.

Thermal control of virulence factors in bacteria: a hot topic

Pathogenic bacteria sense environmental cues, including the local temperature, to control the production of key virulence factors. Thermal regulation can be achieved at the level of DNA, RNA or protein and although many virulence factors are subject to thermal regulation, the exact mechanisms of control are yet to be elucidated in many instances. Understanding how virulence factors are regulated by temperature presents a significant challenge, as gene expression and protein production are often influenced by complex regulatory networks involving multiple transcription factors in bacteria. Here we highlight some recent insights into thermal regulation of virulence in pathogenic bacteria. We focus on bacteria which cause disease in mammalian hosts, which are at a significantly higher temperature than the outside environment. We outline the mechanisms of thermal regulation and how understanding this fundamental aspect of the biology of bacteria has implications for pathogenesis and human health.

Transcriptomic profiling of Yersinia pseudotuberculosis reveals reprogramming of the Crp regulon by temperature and uncovers Crp as a master regulator of small RNAs

One hallmark of pathogenic yersiniae is their ability to rapidly adjust their life-style and pathogenesis upon host entry. In order to capture the range, magnitude and complexity of the underlying gene control mechanisms we used comparative RNA-seq-based transcriptomic profiling of the enteric pathogen Y. pseudotuberculosis under environmental and infection-relevant conditions. We identified 1151 individual transcription start sites, multiple riboswitch-like RNA elements, and a global set of antisense RNAs and previously unrecognized trans-acting RNAs. Taking advantage of these data, we revealed a temperature-induced and growth phase-dependent reprogramming of a large set of catabolic/energy production genes and uncovered the existence of a thermo-regulated ‘acetate switch’, which appear to prime the bacteria for growth in the digestive tract. To elucidate the regulatory architecture linking nutritional status to virulence we also refined the CRP regulon. We identified a massive remodelling of the CRP-controlled network in response to temperature and discovered CRP as a transcriptional master regulator of numerous conserved and newly identified non-coding RNAs which participate in this process. This finding highlights a novel level of complexity of the regulatory network in which the concerted action of transcriptional regulators and multiple non-coding RNAs under control of CRP adjusts the control of Yersinia fitness and virulence to the requirements of their environmental and virulent life-styles.

Many bacterial pathogens cycle between environmental sources and mammalian hosts. Adaptation to the different natural habitats and host niches is achieved through complex regulatory networks which adjust synthesis of the large repertoire of crucial virulence factors and fitness determinants. To uncover underlying control circuits, we determined the first in-depth single-nucleotide resolution transcriptome of Yersinia. This revealed important novel genetic information, such as global locations of transcriptional start sites, non-coding RNAs, potential riboswitches and provided a set of virulence-relevant expression profiles, which constitute a valuable tool for the research community. The analysis further uncovered a temperature-induced global reprogramming of central metabolic functions, likely to support intestinal colonization of the pathogen. This is accompanied by a major reorganization of the CRP regulon, which involves a multitude of regulatory RNAs. The primary consequence is a fine-tuned, coordinated control of metabolism and virulence through a plethora of environmentally controlled regulatory RNAs allowing rapid adaptation and high flexibility during life-style changes.

HmsB enhances biofilm formation in Yersinia pestis

The hmsHFRS operon is responsible for biosynthesis and translocation of biofilm matrix exopolysaccharide. Yersinia pestis expresses the two sole diguanylate cyclases HmsT and HmsD and the sole phosphodiesterase HmsP, which are specific for biosynthesis and degradation, respectively, of 3′,5′-cyclic diguanosine monophosphate (c-di-GMP), a second messenger promoting exopolysaccharide production. In this work, the phenotypic assays indicates that Y. pestis sRNA HmsB enhances the production of c-di-GMP, exopolysaccharide, and biofilm. Further gene regulation experiments disclose that HmsB stimulates the expression of hmsB, hmsCDE, hmsT, and hmsHFRS but represses that of hmsP. HmsB most likely acts as a major activator of biofilm formation in Y. pestis. This is the first report of regulation of Yersinia biofilm formation by a sRNA. Data presented here will promote us to gain a deeper understanding of the complex regulatory circuits controlling Yersinia biofilm formation.

Coregulation of host-adapted metabolism and virulence by pathogenic yersiniae

Deciphering the principles how pathogenic bacteria adapt their metabolism to a specific host microenvironment is critical for understanding bacterial pathogenesis. The enteric pathogenic Yersinia species Yersinia pseudotuberculosis and Yersinia enterocolitica and the causative agent of plague, Yersinia pestis, are able to survive in a large variety of environmental reservoirs (e.g., soil, plants, insects) as well as warm-blooded animals (e.g., rodents, pigs, humans) with a particular preference for lymphatic tissues. In order to manage rapidly changing environmental conditions and interbacterial competition, Yersinia senses the nutritional composition during the course of an infection by special molecular devices, integrates this information and adapts its metabolism accordingly. In addition, nutrient availability has an impact on expression of virulence genes in response to C-sources, demonstrating a tight link between the pathogenicity of yersiniae and utilization of nutrients. Recent studies revealed that global regulatory factors such as the cAMP receptor protein (Crp) and the carbon storage regulator (Csr) system are part of a large network of transcriptional and posttranscriptional control strategies adjusting metabolic changes and virulence in response to temperature, ion and nutrient availability. Gained knowledge about the specific metabolic requirements and the correlation between metabolic and virulence gene expression that enable efficient host colonization led to the identification of new potential antimicrobial targets.

Discovery of ethanol-responsive small RNAs in Zymomonas mobilis

Zymomonas mobilis is a bacterium that can produce ethanol by fermentation. Due to its unique metabolism and efficient ethanol production, Z. mobilis has attracted special interest for biofuel energy applications; an important area of study is the regulation of those specific metabolic pathways. Small RNAs (sRNAs) have been studied as molecules that function as transcriptional regulators in response to cellular stresses. While sRNAs have been discovered in various organisms by computational prediction and experimental approaches, their discovery in Z. mobilis has not yet been reported. In this study, we have applied transcriptome analysis and computational predictions to facilitate identification and validation of 15 novel sRNAs in Z. mobilis. We furthermore characterized their expression in the context of high and low levels of intracellular ethanol. Here, we report that 3 of the sRNAs (Zms2, Zms4, and Zms6) are differentially expressed under aerobic and anaerobic conditions, when low and high ethanol productions are observed, respectively. Importantly, when we tested the effect of ethanol stress on the expression of sRNAs in Z. mobilis, Zms2, Zms6, and Zms18 showed differential expression under 5% ethanol stress conditions. These data suggest that in this organism regulatory RNAs can be associated with metabolic functions involved in ethanol stress responses.

Mechanisms of post-transcriptional gene regulation in bacterial biofilms

Biofilms are characterized by a dense multicellular community of microorganisms that can be formed by the attachment of bacteria to an inert surface and to each other. The development of biofilm involves the initial attachment of planktonic bacteria to a surface, followed by replication, cell-to-cell adhesion to form microcolonies, maturation, and detachment. Mature biofilms are embedded in a self-produced extracellular polymeric matrix composed primarily of bacterial-derived exopolysaccharides, specialized proteins, adhesins, and occasionally DNA. Because the synthesis and assembly of biofilm matrix components is an exceptionally complex process, the transition between its different phases requires the coordinate expression and simultaneous regulation of many genes by complex genetic networks involving all levels of gene regulation. The finely controlled intracellular level of the chemical second messenger molecule, cyclic-di-GMP is central to the post-transcriptional mechanisms governing the switch between the motile planktonic lifestyle and the sessile biofilm forming state in many bacteria. Several other post-transcriptional regulatory mechanisms are known to dictate biofilm development and assembly and these include RNA-binding proteins, small non-coding RNAs, toxin-antitoxin systems, riboswitches, and RNases. Post-transcriptional regulation is therefore a powerful molecular mechanism employed by bacteria to rapidly adjust to the changing environment and to fine tune gene expression to the developmental needs of the cell. In this review, we discuss post-transcriptional mechanisms that influence the biofilm developmental cycle in a variety of pathogenic bacteria.

Genome-wide analysis of small RNAs expressed by Yersinia pestis identifies a regulator of the Yop-Ysc type III secretion system

Small noncoding RNA (sRNA) molecules are integral components of the regulatory machinery for many bacterial species and are known to posttranscriptionally regulate metabolic and stress-response pathways, quorum sensing, virulence factors, and more. The Yop-Ysc type III secretion system (T3SS) is a critical virulence component for the pathogenic Yersinia species, and the regulation of this system is tightly controlled at each step from transcription to translocation of effectors into host cells. The contribution of sRNAs to the regulation of the T3SS in Yersinia has been largely unstudied, however. Previously, our lab identified a role for the sRNA chaperone protein Hfq in the regulation of components of the T3SS in the gastrointestinal pathogen Yersinia pseudotuberculosis. Here we present data demonstrating a similar requirement for Hfq in the closely related species Yersinia pestis. Through deep sequencing analysis of the Y. pestis sRNA-ome, we found 63 previously unidentified putative sRNAs in this species. We identified a Yersinia-specific sRNA, Ysr141, carried by the T3SS plasmid pCD1 that is required for the production of multiple T3SS proteins. In addition, we show that Ysr141 targets an untranslated region upstream of yopJ to posttranscriptionally activate the synthesis of the YopJ protein. Furthermore, Ysr141 may be an unstable and/or processed sRNA, which could contribute to its function in the regulation of the T3SS. The discovery of an sRNA that influences the synthesis of the T3SS adds an additional layer of regulation to this tightly controlled virulence determinant of Y. pestis.

The RNA chaperone Hfq impacts growth, metabolism and production of virulence factors in Yersinia enterocolitica

To adapt to changes in environmental conditions, bacteria regulate their gene expression at the transcriptional but also at the post-transcriptional level, e.g. by small RNAs (sRNAs) which modulate mRNA stability and translation. The conserved RNA chaperone Hfq mediates the interaction of many sRNAs with their target mRNAs, thereby playing a global role in fine-tuning protein production.

In this study, we investigated the significance of Hfq for the enteropathogen Yersina enterocolitica serotype O:8. Hfq facilitated optimal growth in complex and minimal media. Our comparative protein analysis of parental and hfq-negative strains suggested that Hfq promotes lipid metabolism and transport, cell redox homeostasis, mRNA translation and ATP synthesis, and negatively affects carbon and nitrogen metabolism, transport of siderophore and peptides and tRNA synthesis. Accordingly, biochemical tests indicated that Hfq represses ornithine decarboxylase activity, indole production and utilization of glucose, mannitol, inositol and 1,2-propanediol. Moreover, Hfq repressed production of the siderophore yersiniabactin and its outer membrane receptor FyuA. In contrast, hfq mutants exhibited reduced urease production. Finally, strains lacking hfq were more susceptible to acidic pH and oxidative stress. Unlike previous reports in other Gram-negative bacteria, Hfq was dispensable for type III secretion encoded by the virulence plasmid.

Using a chromosomally encoded FLAG-tagged Hfq, we observed increased production of Hfq-FLAG in late exponential and stationary phases. Overall, Hfq has a profound effect on metabolism, resistance to stress and modulates the production of two virulence factors in Y. enterocolitica, namely urease and yersiniabactin.

The exoribonuclease Polynucleotide Phosphorylase influences the virulence and stress responses of yersiniae and many other pathogens

Microbes are incessantly challenged by both biotic and abiotic stressors threatening their existence. Therefore, bacterial pathogens must possess mechanisms to successfully subvert host immune defenses as well as overcome the stress associated with host-cell encounters. To achieve this, bacterial pathogens typically experience a genetic re-programming whereby anti-host/stress factors become expressed and eventually translated into effector proteins. In that vein, the bacterial host-cell induced stress-response is similar to any other abiotic stress to which bacteria respond by up-regulating specific stress-responsive genes. Following the stress encounter, bacteria must degrade unnecessary stress responsive transcripts through RNA decay mechanisms. The three pathogenic yersiniae (Yersinia pestis, Y. pseudo-tuberculosis, and Y. enterocolitica) are all psychrotropic bacteria capable of growth at 4°C; however, cold growth is dependent on the presence of an exoribonuclease, polynucleotide phosphorylase (PNPase). PNPase has also been implicated as a virulence factor in several notable pathogens including the salmonellae, Helicobacter pylori, and the yersiniae [where it typically influences the type three secretion system (TTSS)]. Further, PNPase has been shown to associate with ribonuclease E (endoribonuclease), RhlB (RNA helicase), and enolase (glycolytic enzyme) in several Gram-negative bacteria forming a large, multi-protein complex known as the RNA degradosome. This review will highlight studies demonstrating the influence of PNPase on the virulence potentials and stress responses of various bacterial pathogens as well as focusing on the degradosome-dependent and -independent roles played by PNPase in yersiniae stress responses.

Identification of novel sRNAs in mycobacterial species

Bacterial small RNAs (sRNAs) are short transcripts that typically do not encode proteins and often act as regulators of gene expression through a variety of mechanisms. Regulatory sRNAs have been identified in many species, including Mycobacterium tuberculosis, the causative agent of tuberculosis. Here, we use a computational algorithm to predict sRNA candidates in the mycobacterial species M. smegmatis and M. bovis BCG and confirmed the expression of many sRNAs using Northern blotting. Thus, we have identified 17 and 23 novel sRNAs in M. smegmatis and M. bovis BCG, respectively. We have also applied a high-throughput technique (Deep-RACE) to map the 5' and 3' ends of many of these sRNAs and identified potential regulators of sRNAs by analysis of existing ChIP-seq datasets. The sRNAs identified in this work likely contribute to the unique biology of mycobacteria.

Determination of sRNA expressions by RNA-seq in Yersinia pestis grown in vitro and during infection

Background

Small non-coding RNAs (sRNAs) facilitate host-microbe interactions. They have a central function in the post-transcriptional regulation during pathogenic lifestyles. Hfq, an RNA-binding protein that many sRNAs act in conjunction with, is required for Y. pestis pathogenesis. However, information on how Yersinia pestis modulates the expression of sRNAs during infection is largely unknown.

Methodology and Principal Findings

We used RNA-seq technology to identify the sRNA candidates expressed from Y. pestis grown invitro and in the infected lungs of mice. A total of 104 sRNAs were found, including 26 previously annotated sRNAs, by searching against the Rfam database with 78 novel sRNA candidates. Approximately 89% (93/104) of these sRNAs from Y. pestis are shared with its ancestor Y. pseudotuberculosis. Ninety-seven percent of these sRNAs (101/104) are shared among more than 80 sequenced genomes of 135 Y. pestis strains. These 78 novel sRNAs include 62 intergenic and 16 antisense sRNAs. Fourteen sRNAs were selected for verification by independent Northern blot analysis. Results showed that nine selected sRNA transcripts were Hfq-dependent. Interestingly, three novel sRNAs were identified as new members of the transcription factor CRP regulon. Semi-quantitative analysis revealed that Y. pestis from the infected lungs induced the expressions of six sRNAs including RyhB1, RyhB2, CyaR/RyeE, 6S RNA, RybB and sR039 and repressed the expressions of four sRNAs, including CsrB, CsrC, 4.5S RNA and sR027.

Conclusions and Significance

This study is the first attempt to subject RNA from Y. pestis-infected samples to direct high-throughput sequencing. Many novel sRNAs were identified and the expression patterns of relevant sRNAs in Y. pestis during invitro growth and invivo infection were revealed. The annotated sRNAs accounted for the most abundant sRNAs either expressed in bacteria grown invitro or differentially expressed in the infected lungs. These findings suggested these sRNAs may have important functions in Y. pestis physiology or pathogenesis.

Small RNA-mediated regulation of host-pathogen interactions

The rise in antimicrobial drug resistance, alongside the failure of conventional research to discover new antibiotics, will inevitably lead to a public health crisis that can drastically curtail our ability to combat infectious disease. Thus, there is a great global health need for development of antimicrobial countermeasures that target novel cell molecules or processes. RNA represents a largely unexploited category of potential targets for antimicrobial design. For decades, control of cellular behavior was thought to be the exclusive purview of protein-based regulators. The recent discovery of small RNAs (sRNAs) as a universal class of powerful RNA-based regulatory biomolecules has the potential to revolutionize our understanding of gene regulation in practically all biological functions. In general, sRNAs regulate gene expression by base-pairing with multiple downstream target mRNAs to prevent translation of mRNA into protein. In this review, we will discuss recent studies that document discovery of bacterial, viral, and human sRNAs and their molecular mechanisms in regulation of pathogen virulence and host immunity. Illuminating the functional roles of sRNAs in virulence and host immunity can provide the fundamental knowledge for development of next-generation antibiotics using sRNAs as novel targets.