Transcriptomic profiling of the oyster pathogen Vibrio splendidus opens a window on the evolutionary dynamics of the small RNA repertoire in the Vibrio genus

High-Resolution Small RNAs Landscape Provides Insights into Alkane Adaptation in the Marine Alkane-Degrader Alcanivorax dieselolei B-5

Alkanes are widespread in the ocean, and Alcanivorax is one of the most ubiquitous alkane-degrading bacteria in the marine ecosystem. Small RNAs (sRNAs) are usually at the heart of regulatory pathways, but sRNA-mediated alkane metabolic adaptability still remains largely unknown due to the difficulties of identification. Here, differential RNA sequencing (dRNA-seq) modified with a size selection (~50-nt to 500-nt) strategy was used to generate high-resolution sRNAs profiling in the model species Alcanivorax dieselolei B-5 under alkane (n-hexadecane) and non-alkane (acetate) conditions. As a result, we identified 549 sRNA candidates at single-nucleotide resolution of 5'-ends, 63.4% of which are with transcription start sites (TSSs), and 36.6% of which are with processing sites (PSSs) at the 5'-ends. These sRNAs originate from almost any location in the genome, regardless of intragenic (65.8%), antisense (20.6%) and intergenic (6.2%) regions, and RNase E may function in the maturation of sRNAs. Most sRNAs locally distribute across the 15 reference genomes of Alcanivorax, and only 7.5% of sRNAs are broadly conserved in this genus. Expression responses to the alkane of several core conserved sRNAs, including 6S RNA, M1 RNA and tmRNA, indicate that they may participate in alkane metabolisms and result in more actively global transcription, RNA processing and stresses mitigation. Two novel CsrA-related sRNAs are identified, which may be involved in the translational activation of alkane metabolism-related genes by sequestering the global repressor CsrA. The relationships of sRNAs with the characterized genes of alkane sensing (ompS), chemotaxis (mcp, cheR, cheW2), transporting (ompT1, ompT2, ompT3) and hydroxylation (alkB1, alkB2, almA) were created based on the genome-wide predicted sRNA-mRNA interactions. Overall, the sRNA landscape lays the ground for uncovering cryptic regulations in critical marine bacterium, among which both the core and species-specific sRNAs are implicated in the alkane adaptive metabolisms.

[Regulatory role of small RNA srn821978 in mutacin IV expression in Streptococcus mutans]

OBJECTIVE

To analyze the role of small RNA srn821798 in posttranscriptional regulation of mutacin IV expression in Streptococcus mutans.

METHODS

The potential target genes of srn821978 were predicted using RNAhybrid, RNAPredator and IntaRNA. We collected 10 Streptococcus mutans (S.muans) strains with high expression of mutacin IV and another 10 S.muans strains that did not express mutacin IV screened by inhibition zone test, and the expression levels of srn821798 and the candidate target genes in these strains were detected by qPCR. Using synthesized mimics and inhibitors of srn821798, we constructed S.muans strains with high or low srn821798 expression via electroporation based on the standard strain of S.muans UA159, and analyzed the expression levels of srn821798 and its candidate target genes in these strains. We also examined the binding ability of srn821798 to its target gene sepM using electrophoresis and a dual- luciferase reporter system.

RESULTS

The expression levels of the candidate target genes of srn821798 including sepM, comD, comE, nlmA and nlmB were significantly higher while the expression level of srn821798 was significantly lower in clinical S.muans strains with high expression of mutacin IV than in those without mutacin IV expression (P < 0.05). Although the expression levels of the candidate target genes in strains with up- regulated or down- regulated srn821798 expression did not differ significantly from those in the standard strain, the expression level of sepM showed a trend of differential distribution, and srn821798 was predicted to have a strong binding ability to sepM action site.

CONCLUSION

srn821798 may play a regulatory role in the expression of mutacin IV in S.muans, but the underlying mechanism remains to be explored.

The Pseudoalteromonas multipartite genome: distribution and expression of pangene categories, and a hypothesis for the origin and evolution of the chromid

Bacterial genomes typically consist of one large chromosome, but can also include secondary replicons. These so-called multipartite genomes are scattered on the bacterial tree of life with the majority of cases belonging to Proteobacteria. Within the class gamma-proteobacteria, multipartite genomes are restricted to the two families Vibrionaceae and Pseudoalteromonadaceae. Whereas the genome of vibrios is well studied, information on the Pseudoalteromonadaceae genome is much scarcer. We have studied Pseudoalteromonadaceae with respect to the origin of the chromid, how pangene categories are distributed, how genes are expressed relative to their genomic location, and identified chromid hallmark genes. We calculated the Pseudoalteromonadaceae pangenome based on 25 complete genomes and found that core/softcore are significantly overrepresented in late replicating sectors of the chromid, regardless of how the chromid is replicated. On the chromosome, core/softcore and shell/cloud genes are only weakly overrepresented at the chromosomal replication origin and termination sequences, respectively. Gene expression is trending downwards with increasing distance from the chromosomal oriC, whereas the chromidal expression pattern is more complex. Moreover, we identified 78 chromid hallmark genes, and BLASTp searches suggest that the majority of them were acquired from the ancestral gene pool of Alteromonadales. Finally, our data strongly suggest that the chromid originates from a plasmid that was acquired in a relatively recent event. In summary, this study extends our knowledge on multipartite genomes, and helps us understand how and why secondary replicons are acquired, why they are maintained, and how they are shaped by evolution.

A Small Regulatory RNA Generated from the malK 5' Untranslated Region Targets Gluconeogenesis in Vibrio Species

Vsr217 is a small RNA from Vibrio tasmaniensis LGP32, a pathogen associated with mortality events affecting juvenile oysters. The vsr217 gene is located within the 5' untranslated region (UTR) of malK, encoding the ATPase component of the maltose importer, and is conserved within the genus Vibrio. In the presence of maltose, vsr217 is regulated by MalT, the positive regulator of the maltose regulon. vsr217 is required in cis for the full expression of malK. In addition, Vsr217 acts in trans to downregulate the expression of fbp encoding fructose-1,6-bisphosphatase, an enzyme involved in gluconeogenesis. Thus, in the presence of maltose, the induction of Vsr217 is expected to promote glycolysis by negatively regulating the expression of a key enzyme of gluconeogenesis. IMPORTANCE Juvenile pacific oysters have been subject in recent years to summer mortality episodes with deep economic consequences. The pathogen Vibrio tasmaniensis has been associated with such mortality events. For bacterial pathogens, survival within the host requires profound metabolic adaptations according to available resources. All kinds of regulatory elements, including noncoding RNAs, orchestrate this response. Oysters are rich in glycogen, a precursor of maltose, and we previously reported that V. tasmaniensis maltose-regulated genes are strongly induced during oyster infection. Here, we report the dual mechanism by which a small regulatory RNA, generated from the 5' untranslated region of a gene belonging to the maltose regulon, acts both in cis and trans. In cis, it stimulates growth on maltose, and in trans, it downregulates the expression of a gene associated with gluconeogenesis, thus coordinating maltose utilization with central carbon metabolism.

Biosecurity implications of drifting marine plastic debris: Current knowledge and future research

The introduction and spread of marine non-indigenous species (NIS) and pathogens into new habitats are a major threat to biodiversity, ecosystem services, human health, and can have substantial economic consequences. Shipping is considered the main vector for marine biological invasions; less well understood is the increased spread of marine NIS and pathogens rafting on marine plastic debris (MPD). Despite an increasing research interest and recent progress in characterizing the plastisphere, this manuscript highlights critical knowledge gaps and research priorities towards a better understanding of the biosecurity implications of MPD. We advocate for future research to (i) investigate plastisphere community succession and the factors influencing NIS propagules and pathogens recruitment through robust experimental investigations; (ii) combine microscopy and molecular approaches to effectively assess the presence of specific taxa; (iii) include additional genetic markers to thoroughly characterize the biodiversity associated with MPD and explore the presence of specific marine pests.

Vibrionaceae core, shell and cloud genes are non-randomly distributed on Chr 1: An hypothesis that links the genomic location of genes with their intracellular placement

Background

The genome of Vibrionaceae bacteria, which consists of two circular chromosomes, is replicated in a highly ordered fashion. In fast-growing bacteria, multifork replication results in higher gene copy numbers and increased expression of genes located close to the origin of replication of Chr 1 (ori1). This is believed to be a growth optimization strategy to satisfy the high demand of essential growth factors during fast growth. The relationship between ori1-proximate growth-related genes and gene expression during fast growth has been investigated by many researchers. However, it remains unclear which other gene categories that are present close to ori1 and if expression of all ori1-proximate genes is increased during fast growth, or if expression is selectively elevated for certain gene categories.

Results

We calculated the pangenome of all complete genomes from the Vibrionaceae family and mapped the four pangene categories, core, softcore, shell and cloud, to their chromosomal positions. This revealed that core and softcore genes were found heavily biased towards ori1, while shell genes were overrepresented at the opposite part of Chr 1 (i.e., close to ter1). RNA-seq of Aliivibrio salmonicida and Vibrio natriegens showed global gene expression patterns that consistently correlated with chromosomal distance to ori1. Despite a biased gene distribution pattern, all pangene categories contributed to a skewed expression pattern at fast-growing conditions, whereas at slow-growing conditions, softcore, shell and cloud genes were responsible for elevated expression.

Conclusion

The pangene categories were non-randomly organized on Chr 1, with an overrepresentation of core and softcore genes around ori1, and overrepresentation of shell and cloud genes around ter1. Furthermore, we mapped our gene distribution data on to the intracellular positioning of chromatin described for V. cholerae, and found that core/softcore and shell/cloud genes appear enriched at two spatially separated intracellular regions. Based on these observations, we hypothesize that there is a link between the genomic location of genes and their cellular placement.

The ilvGMEDA Operon Is Regulated by Transcription Attenuation in Vibrio alginolyticus ZJ-T

Bacteria synthesize amino acids according to their availability in the environment or, in the case of pathogens, within the host. We explored the regulation of the biosynthesis of branched-chain amino acids (BCAAs) (l-leucine, l-valine, and l-isoleucine) in Vibrio alginolyticus, a marine fish and shellfish pathogen and an emerging opportunistic human pathogen. In this species, the ilvGMEDA operon encodes the main pathway for biosynthesis of BCAAs. Its upstream regulatory region shows no sequence similarity to the corresponding region in Escherichia coli or other Enterobacteriaceae, and yet we show that this operon is regulated by transcription attenuation. The translation of a BCAA-rich peptide encoded upstream of the structural genes provides an adaptive response similar to the E. coli canonical model. This study of a nonmodel Gram-negative organism highlights the mechanistic conservation of transcription attenuation despite the absence of primary sequence conservation.IMPORTANCE This study analyzes the regulation of the biosynthesis of branched-chain amino acids (leucine, valine, and isoleucine) in Vibrio alginolyticus, a marine bacterium that is pathogenic to fish and humans. The results highlight the conservation of the main regulatory mechanism with that of the enterobacterium Escherichia coli, suggesting that such a mechanism appeared early during the evolution of Gram-negative bacteria, allowing adaptation to a wide range of environments.

csrB Gene Duplication Drives the Evolution of Redundant Regulatory Pathways Controlling Expression of the Major Toxic Secreted Metalloproteases in Vibrio tasmaniensis LGP32

The conserved CsrB sRNAs are an example of sibling sRNAs, i.e., sRNAs which are present in multiple copies in genomes. This report illustrates how new copies arise through gene duplication events and highlights two evolutionary advantages of having such multiple copies: differential regulation of the multiple copies allows integration of different input signals into the regulatory network of which they are parts, and the high redundancy that they provide confers a strong robustness to the system.

CsrBs are bacterial highly conserved and multiple-copy noncoding small RNAs (sRNAs) that play major roles in cell physiology and virulence. In the Vibrio genus, they are known to be regulated by the two-component system VarS/VarA. They modulate the well-characterized quorum sensing pathway controlling virulence and luminescence in Vibrio cholerae and Vibrio harveyi, respectively. Remarkably, Vibrio tasmaniensis LGP32, an oyster pathogen that belongs to the Splendidus clade, was found to have four copies of csrB, named csrB1-4, compared to two to three copies in other Vibrio species. Here, we show that the extra csrB4 copy results from a csrB3 gene duplication, a characteristic of the Splendidus clade. Interestingly, csrB genes are regulated in different ways in V. tasmaniensis, with csrB1 expression being independent of the VarS/VarA system. We found that a complex regulatory network involving CsrBs, quorum sensing, and the stationary-phase sigma factor σS redundantly but differentially controls the production of two secreted metalloproteases, Vsm and PrtV, the former being a major determinant of the V. tasmaniensis extracellular product toxicity. In particular, we identified a novel VarS/VarA-dependent but CsrB-independent pathway that controls positively both Vsm production and PrtV production as well as rpoS expression. Altogether, our data show that a csrB gene duplication event in V. tasmaniensis supported the evolution of the regulatory network controlling the expression of major toxic secreted metalloproteases, thereby increasing redundancy and enabling the integration of additional input signals.

IMPORTANCE The conserved CsrB sRNAs are an example of sibling sRNAs, i.e., sRNAs which are present in multiple copies in genomes. This report illustrates how new copies arise through gene duplication events and highlights two evolutionary advantages of having such multiple copies: differential regulation of the multiple copies allows integration of different input signals into the regulatory network of which they are parts, and the high redundancy that they provide confers a strong robustness to the system.

Defining the transcription landscape of the Gram-negative marine bacterium Vibrio harveyi

Vibrio harveyi is a Gram-negative pathogenic bacterium ubiquitously present in natural aquatic systems. Although environmental adaptability in V. harveyi may be enabled by profound reprogramming of gene expression previously observed during responses to starvation, suboptimal temperatures and other stress factors, the key characteristics of V. harveyi transcripts and operons, such as their boundaries and size as well as location of small RNA genes, remain largely unknown. To reveal the main features of the V. harveyi transcriptome, total RNA of this organism was analyzed by differential RNA sequencing (dRNA-seq). Analysis of the dRNA-seq data made it possible to define the primary transcriptome of V. harveyi along with cis-acting regulatory elements (riboswitches and leader sequences) and new genes. The latter encode a number of putative polypeptides and new phylogenetically conserved antisense RNAs potentially involved in the post-transcriptional control of gene expression.

Identification of a Novel Small RNA srvg23535 in Vibrio alginolyticus ZJ-T and Its Characterization With Phenotype MicroArray Technology

Small non-coding RNAs (sRNAs) are important modulators of gene expression and are involved in the pathogenesis and survival of prokaryotes. However, few studies have been conducted with Vibrio alginolyticus, which limits our ability to probe the regulation of virulence and environmental adaptation by sRNAs in this opportunistic pathogen. In this study, the sRNA candidate srvg23535 was identified in V. alginolyticus ZJ-T. The precise transcript end, secondary structure, and sequence conservation were determined. A srvg23535 null mutant was constructed and characterized by using Phenotype MicroArray (PM) technology. In silico target prediction was conducted by IntaRNA and TargetRNA2. Subsequently, a 107 nt transcript was validated with a sigma70 promoter at the 5' end and a Rho-independent terminator at the 3' end. The sRNA srvg23535 had four stem-loop structures and was conserved among Vibrio harveyi, Vibrio parahaemolyticus, and Vibrio splendidus. Deletion of srvg23535 in V. alginolyticus ZJ-T led to a weaker utilization of D-mannose, D-melibiose, lactulose, and inosine as carbon sources but stronger utilization of L-cysteine as nitrogen source. Moreover, the srvg2353 mutant showed stronger resistance to osmotic stress but weaker resistance to pH stress. Additionally, a total of 22 common targets were identified and several were related to the observed phenotype of the mutant. This study indicated that the novel sRNA, srvg23535, is conserved and restricted to Vibrio spp., affecting the utilization of several carbon and nitrogen sources and the response to osmotic and pH stress. These results extend our understanding of sRNA regulation in V. alginolyticus and provide a significant resource for the further study of the precise target mRNAs of srvg23535, which may provide targets for antibacterial therapeutic or attenuated vaccines against Vibrio spp.

The conserved regulatory RNA RsaE down-regulates the arginine degradation pathway in Staphylococcus aureus

RsaE is a regulatory RNA highly conserved amongst Firmicutes that lowers the amount of mRNAs associated with the TCA cycle and folate metabolism. A search for new RsaE targets in Staphylococcus aureus revealed that in addition to previously described substrates, RsaE down-regulates several genes associated with arginine catabolism. In particular, RsaE targets the arginase rocF mRNA via direct interactions involving G-rich motifs. Two duplicated C-rich motifs of RsaE can independently downregulate rocF expression. The faster growth rate of ΔrsaE compared to its parental strain in media containing amino acids as sole carbon source points to an underlying role for RsaE in amino acid catabolism. Collectively, the data support a model in which RsaE acts as a global regulator of functions associated with metabolic adaptation.

Expansion of the SOS regulon of Vibrio cholerae through extensive transcriptome analysis and experimental validation

BACKGROUND

The SOS response is an almost ubiquitous response of cells to genotoxic stresses. The full complement of genes in the SOS regulon for Vibrio species has only been addressed through bioinformatic analyses predicting LexA binding box consensus and in vitro validation. Here, we perform whole transcriptome sequencing from Vibrio cholerae treated with mitomycin C as an SOS inducer to characterize the SOS regulon and other pathways affected by this treatment.

RESULTS

Comprehensive transcriptional profiling allowed us to define the full landscape of promoters and transcripts active in V. cholerae. We performed extensive transcription start site (TSS) mapping as well as detection/quantification of the coding and non-coding RNA (ncRNA) repertoire in strain N16961. To improve TSS detection, we developed a new technique to treat RNA extracted from cells grown in various conditions. This allowed for identification of 3078 TSSs with an average 5'UTR of 116 nucleotides, and peak distribution between 16 and 64 nucleotides; as well as 629 ncRNAs. Mitomycin C treatment induced transcription of 737 genes and 28 ncRNAs at least 2 fold, while it repressed 231 genes and 17 ncRNAs. Data analysis revealed that in addition to the core genes known to integrate the SOS regulon, several metabolic pathways were induced. This study allowed for expansion of the Vibrio SOS regulon, as twelve genes (ubiEJB, tatABC, smpA, cep, VC0091, VC1190, VC1369-1370) were found to be co-induced with their adjacent canonical SOS regulon gene(s), through transcriptional read-through. Characterization of UV and mitomycin C susceptibility for mutants of these newly identified SOS regulon genes and other highly induced genes and ncRNAs confirmed their role in DNA damage rescue and protection.

CONCLUSIONS

We show that genotoxic stress induces a pervasive transcriptional response, affecting almost 20% of the V. cholerae genes. We also demonstrate that the SOS regulon is larger than previously known, and its syntenic organization is conserved among Vibrio species. Furthermore, this specific co-localization is found in other γ-proteobacteria for genes recN-smpA and rmuC-tatABC, suggesting SOS regulon conservation in this phylum. Finally, we comment on the limitations of widespread NGS approaches for identification of all RNA species in bacteria.

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.

Assessment of Bona Fide sRNAs in Staphylococcus aureus

Bacterial regulatory RNAs have been extensively studied for over a decade, and are progressively being integrated into the complex genetic regulatory network. Transcriptomic arrays, recent deep-sequencing data and bioinformatics suggest that bacterial genomes produce hundreds of regulatory RNAs. However, while some have been authenticated, the existence of the others varies according to strains and growth conditions, and their detection fluctuates with the methodologies used for data acquisition and interpretation. For example, several small RNA (sRNA) candidates are now known to be parts of UTR transcripts. Accurate annotation of regulatory RNAs is a complex task essential for molecular and functional studies. We defined bona fide sRNAs as those that (i) likely act in trans and (ii) are not expressed from the opposite strand of a coding gene. Using published data and our own RNA-seq data, we reviewed hundreds of Staphylococcus aureus putative regulatory RNAs using the DETR'PROK computational pipeline and visual inspection of expression data, addressing the question of which transcriptional signals correspond to sRNAs. We conclude that the model strain HG003, a NCTC8325 derivative commonly used for S. aureus genetic regulation studies, has only about 50 bona fide sRNAs, indicating that these RNAs are less numerous than commonly stated. Among them, about half are associated to the S. aureus sp. core genome and a quarter are possibly expressed in other Staphylococci. We hypothesize on their features and regulation using bioinformatic approaches.

Carrageenan catabolism is encoded by a complex regulon in marine heterotrophic bacteria

Macroalgae contribute substantially to primary production in coastal ecosystems. Their biomass, mainly consisting of polysaccharides, is cycled into the environment by marine heterotrophic bacteria using largely uncharacterized mechanisms. Here we describe the complete catabolic pathway for carrageenans, major cell wall polysaccharides of red macroalgae, in the marine heterotrophic bacterium Zobellia galactanivorans. Carrageenan catabolism relies on a multifaceted carrageenan-induced regulon, including a non-canonical polysaccharide utilization locus (PUL) and genes distal to the PUL, including a susCD-like pair. The carrageenan utilization system is well conserved in marine Bacteroidetes but modified in other phyla of marine heterotrophic bacteria. The core system is completed by additional functions that might be assumed by non-orthologous genes in different species. This complex genetic structure may be the result of multiple evolutionary events including gene duplications and horizontal gene transfers. These results allow for an extension on the definition of bacterial PUL-mediated polysaccharide digestion.

Carrageenans, major cell wall polysaccharides of red macroalgae, are metabolised by marine heterotrophic bacteria through unclear mechanisms. Here, the authors identify an unusual polysaccharide-utilization locus encoding carrageenan catabolism in a marine bacterium, and characterise the complete pathway.

Unraveling the evolution and coevolution of small regulatory RNAs and coding genes in Listeria

BACKGROUND

Small regulatory RNAs (sRNAs) are widely found in bacteria and play key roles in many important physiological and adaptation processes. Studying their evolution and screening for events of coevolution with other genomic features is a powerful way to better understand their origin and assess a common functional or adaptive relationship between them. However, evolution and coevolution of sRNAs with coding genes have been sparsely investigated in bacterial pathogens.

RESULTS

We designed a robust and generic phylogenomics approach that detects correlated evolution between sRNAs and protein-coding genes using their observed and inferred patterns of presence-absence in a set of annotated genomes. We applied this approach on 79 complete genomes of the Listeria genus and identified fifty-two accessory sRNAs, of which most were present in the Listeria common ancestor and lost during Listeria evolution. We detected significant coevolution between 23 sRNA and 52 coding genes and inferred the Listeria sRNA-coding genes coevolution network. We characterized a main hub of 12 sRNAs that coevolved with genes encoding cell wall proteins and virulence factors. Among them, an sRNA specific to L. monocytogenes species, rli133, coevolved with genes involved either in pathogenicity or in interaction with host cells, possibly acting as a direct negative post-transcriptional regulation.

CONCLUSIONS

Our approach allowed the identification of candidate sRNAs potentially involved in pathogenicity and host interaction, consistent with recent findings on known pathogenicity actors. We highlight four sRNAs coevolving with seven internalin genes, some of which being important virulence factors in Listeria.

Impact of bacterial sRNAs in stress responses

Bacterial life is harsh and involves numerous environmental and internal challenges that are perceived as stresses. Consequently, adequate responses to survive, cope with, and counteract stress conditions have evolved. In the last few decades, a class of small, non-coding RNAs (sRNAs) has been shown to be involved as key players in stress responses. This review will discuss — primarily from an enterobacterial perspective — selected stress response pathways that involve antisense-type sRNAs. These include themes of how bacteria deal with severe envelope stress, threats of DNA damage, problems with poisoning due to toxic sugar intermediates, issues of iron homeostasis, and nutrient limitation/starvation. The examples discussed highlight how stress relief can be achieved, and how sRNAs act mechanistically in regulatory circuits. For some cases, we will propose scenarios that may suggest why contributions from post-transcriptional control by sRNAs, rather than transcriptional control alone, appear to be a beneficial and universally selected feature.

A Genome-Wide Prediction and Identification of Intergenic Small RNAs by Comparative Analysis in Mesorhizobium huakuii 7653R

In bacteria, small non-coding RNAs (sRNAs) are critical regulators of cellular adaptation to changes in metabolism, physiology, or the external environment. In the last decade, more than 2000 of sRNA families have been reported in the Rfam database and have been shown to exert various regulatory functions in bacterial transcription and translation. However, little is known about sRNAs and their functions in Mesorhizobium. Here, we predicted putative sRNAs in the intergenic regions (IGRs) of M. huakuii 7653R by genome-wide comparisons with four related Mesorhizobial strains. The expression and transcribed regions of candidate sRNAs were analyzed using a set of high-throughput RNA deep sequencing data. In all, 39 candidate sRNAs were found, with 5 located in the symbiotic megaplasmids and 34 in the chromosome of M. huakuii 7653R. Of these, 24 were annotated as functional sRNAs in the Rfam database and 15 were recognized as putative novel sRNAs. The expression of nine selected sRNAs was confirmed by Northern blotting, and most of the nine selected sRNAs were highly expressed in 28 dpi nodules and under symbiosis-mimicking conditions. For those putative novel sRNAs, functional categorizations of their target genes were performed by analyzing the enriched GO terms. In addition, MH_s15 was shown to be an abundant and conserved sRNA.

Divergence of functional effects among bacterial sRNA paralogs

BACKGROUND

Non-coding small RNAs (sRNAs) regulate a variety of important biological processes across all life domains, including bacteria. However, little is known about the functional evolution of sRNAs in bacteria, which might occur via changes in sRNA structure and/or stability or changes in interactions between sRNAs and their associated regulatory networks, including target mRNAs. The sRNA Pxr functions as a developmental gatekeeper in the model cooperative bacterium Myxococcus xanthus. Specifically, Pxr prevents the initiation of fruiting body development when nutrients are abundant. Previous work has shown that Pxr appears to have a recent origin within a sub-clade of the myxobacteria, which allowed us to infer the most recent common ancestor of pxr and examine the divergence of Pxr since its origin.

RESULTS

To test for inter-specific divergence in functional effects, extant pxr homologs from several species and their inferred ancestor were introduced into an M. xanthus deletion mutant lacking pxr. Both the inferred ancestral pxr and all extant alleles from species containing only one copy of pxr were found to control development in M. xanthus in a qualitatively similar manner to the native M. xanthus allele. However, multiple paralogs present in Cystobacter species exhibited divergent effects, with two paralogs controlling M. xanthus development but two others failing to do so. These differences may have occurred through changes in gene expression caused by apparent structural differences in the sRNA variants encoded by these paralogs.

CONCLUSIONS

Taken together, our results suggest that Pxr plays a common fundamental role in developmental gene regulation across diverse species of myxobacteria but also that the functional effects of some Pxr variants may be evolving in some lineages.

Construction of a fur null mutant and RNA-sequencing provide deeper global understanding of the Aliivibrio salmonicida Fur regulon

BACKGROUND

The ferric uptake regulator (Fur) is a transcription factor and the main regulator of iron acquisition in prokaryotes. When bound to ferric iron, Fur recognizes its DNA binding site and generally executes its function by repressing transcription of its target genes. Due to its importance in virulence, the Fur regulon is well studied for several model bacteria. In our previous work, we used computational predictions and microarray to gain insights into Fur-regulation in Aliivibrio salmonicida, and have identified a number of genes and operons that appear to be under direct control of Fur. To provide a more accurate and deeper global understanding of the biological role of Fur we have now generated an A. salmonicida fur knock-out strain and used RNA-sequencing to compare gene expression between the wild-type and fur null mutant strains.

RESULTS

An A. salmonicida fur null mutant strain was constructed. Biological assays demonstrate that deletion of fur results in loss of fitness, with reduced growth rates, and reduced abilities to withstand low-iron conditions, and oxidative stress. When comparing expression levels in the wild-type and the fur null mutant we retrieved 296 differentially expressed genes distributed among 18 of 21 functional classes of genes. A gene cluster encoding biosynthesis of the siderophore bisucaberin represented the highest up-regulated genes in the fur null mutant. Other highly up-regulated genes all encode proteins important for iron acquisition. Potential targets for the RyhB sRNA was predicted from the list of down-regulated genes, and significant complementarities were found between RyhB and mRNAs of the fur, sodB, cysN and VSAL_I0422 genes. Other sRNAs with potential functions in iron homeostasis were identified.

CONCLUSION

The present work provides by far the most comprehensive and deepest understanding of the Fur regulon in A. salmonicida to date. Our data also contribute to a better understanding of how Fur plays a key role in iron homeostasis in bacteria in general, and help to show how Fur orchestrates iron uptake when iron levels are extremely low.

The malS-5′UTR regulates hisG, a key gene in the histidine biosynthetic pathway in Salmonella enterica serovar Typhi

Bacterial noncoding RNAs (ncRNA) regulate diverse cellular processes, including virulence and environmental fitness. The malS 5′ untranslated region (named malS-5′UTR) was identified as a regulatory ncRNA that increases the invasive capacity of Salmonella enterica serovar Typhi. An IntaRNA search suggested base pairing between malS-5′UTR and hisG mRNA, a key gene in the histidine biosynthetic pathway. Overexpression of malS-5′UTR markedly reduced bacterial growth in minimal medium without histidine. Overexpression of malS-5′UTR increased mRNA from his operon genes, independently of the bax gene, and decreased HisG protein in Salmonella Typhi. RNA structure analysis showed base pairing of the malS-5′UTR RNA with the hisG mRNA across the ribosome binding site. Thus, we propose that malS-5′UTR inhibited hisG translation, probably by base pairing to the Shine–Dalgarno sequence.

Accelerating Discovery and Functional Analysis of Small RNAs with New Technologies

Over the past decade, bacterial small RNAs (sRNAs) have gone from a biological curiosity to being recognized as a major class of regulatory molecules. High-throughput methods for sampling the transcriptional output of bacterial cells demonstrate that sRNAs are universal features of bacterial transcriptomes, are plentiful, and appear to vary extensively over evolutionary time. With ever more bacteria coming under study, the question becomes how can we accelerate the discovery and functional characterization of sRNAs in diverse organisms. New technologies built on high-throughput sequencing are emerging that can rapidly provide global insight into the numbers and functions of sRNAs in bacteria of interest, providing information that can shape hypotheses and guide research. In this review, we describe recent developments in transcriptomics (RNA-seq) and functional genomics that we expect to help us develop an integrated, systems-level view of sRNA biology in bacteria.

The emergence of Vibrio pathogens in Europe: ecology, evolution, and pathogenesis (Paris, 11-12th March 2015)

Global change has caused a worldwide increase in reports of Vibrio-associated diseases with ecosystem-wide impacts on humans and marine animals. In Europe, higher prevalence of human infections followed regional climatic trends with outbreaks occurring during episodes of unusually warm weather. Similar patterns were also observed in Vibrio-associated diseases affecting marine organisms such as fish, bivalves and corals. Basic knowledge is still lacking on the ecology and evolutionary biology of these bacteria as well as on their virulence mechanisms. Current limitations in experimental systems to study infection and the lack of diagnostic tools still prevent a better understanding of Vibrio emergence. A major challenge is to foster cooperation between fundamental and applied research in order to investigate the consequences of pathogen emergence in natural Vibrio populations and answer federative questions that meet societal needs. Here we report the proceedings of the first European workshop dedicated to these specific goals of the Vibrio research community by connecting current knowledge to societal issues related to ocean health and food security.

5′-UTR of malS increases the invasive capacity of Salmonella enterica serovar Typhi by influencing the expression of bax

ABSTRACT 

Aim: An RNA-seq analysis recently identified a 236-nucleotide transcript upstream from malS in Salmonella enterica serovar Typhi. Here, we investigated its molecular characteristics and function. Materials & methods: RACE and northern blotting were used to determine the molecular characteristics of the sequence, and mutagenesis, microarray, immunoblotting and an invasion assay were used to investigate the functions of the transcript. Results: The transcript was identified as the malS 5′-untranslated region (UTR), which could influence the expression of the flagellar and SPI-1 genes and the invasion of HeLa cells by S. Typhi. Deletion of bax increased the expression of the invasion genes and the invasive capacity of S. Typhi, whereas the expression of the malS 5′-UTR reduced the expression of bax. Conclusion: The malS 5′-UTR reduces the expression of bax and increases the invasive capacity of S. Typhi.

How do base-pairing small RNAs evolve?

The increasing numbers of characterized base-pairing small RNAs (sRNAs) and the identification of these regulators in a broad range of bacteria are allowing comparisons between species and explorations of sRNA evolution. In this review, we describe some examples of trans-encoded base-pairing sRNAs that are species-specific and others that are more broadly distributed. We also describe examples of sRNA orthologs where different features are conserved. These examples provide the background for a discussion of mechanisms of sRNA evolution and selective pressures on the sRNAs and their mRNA target(s).

BactPepDB: a database of predicted peptides from a exhaustive survey of complete prokaryote genomes

With the recent progress in complete genome sequencing, mining the increasing amount of genomic information available should in theory provide the means to discover new classes of peptides. However, annotation pipelines often do not consider small reading frames likely to be expressed. BactPepDB, available online at http://bactpepdb.rpbs.univ-paris-diderot.fr, is a database that aims at providing an exhaustive re-annotation of all complete prokaryotic genomes—chromosomal and plasmid DNA—available in RefSeq for coding sequences ranging between 10 and 80 amino acids. The identified peptides are classified as (i) previously identified in RefSeq, (ii) entity-overlapping (intragenic) or intergenic, and (iii) potential pseudogenes—intergenic sequences corresponding to a portion of a previously annotated larger gene. Additional information is related to homologs within order, predicted signal sequence, transmembrane segments, disulfide bonds, secondary structure, and the existence of a related 3D structure in the Protein Databank. As a result, BactPepDB provides insights about candidate peptides, and provides information about their conservation, together with some of their expected biological/structural features. The BactPepDB interface allows to search for candidate peptides in the database, or to search for peptides similar to a query, according to the multiple properties predicted or related to genomic localization.

Database URL:http://www.yeastgenome.org/

Primary transcriptome map of the hyperthermophilic archaeon Thermococcus kodakarensis

Background

Prokaryotes have relatively small genomes, densely-packed with protein-encoding sequences. RNA sequencing has, however, revealed surprisingly complex transcriptomes and here we report the transcripts present in the model hyperthermophilic Archaeon, Thermococcus kodakarensis, under different physiological conditions.

Results

Sequencing cDNA libraries, generated from RNA isolated from cells under different growth and metabolic conditions has identified >2,700 sites of transcription initiation, established a genome-wide map of transcripts, and consensus sequences for transcription initiation and post-transcription regulatory elements. The primary transcription start sites (TSS) upstream of 1,254 annotated genes, plus 644 primary TSS and their promoters within genes, are identified. Most mRNAs have a 5'-untranslated region (5'-UTR) 10 to 50 nt long (median = 16 nt), but ~20% have 5'-UTRs from 50 to 300 nt long and ~14% are leaderless. Approximately 50% of mRNAs contain a consensus ribosome binding sequence. The results identify TSS for 1,018 antisense transcripts, most with sequences complementary to either the 5'- or 3'-region of a sense mRNA, and confirm the presence of transcripts from all three CRISPR loci, the RNase P and 7S RNAs, all tRNAs and rRNAs and 69 predicted snoRNAs. Two putative riboswitch RNAs were present in growing but not in stationary phase cells. The procedure used is designed to identify TSS but, assuming that the number of cDNA reads correlates with transcript abundance, the results also provide a semi-quantitative documentation of the differences in T. kodakarensis genome expression under different growth conditions and confirm previous observations of substrate-dependent specific gene expression. Many previously unanticipated small RNAs have been identified, some with relative low GC contents (≤50%) and sequences that do not fold readily into base-paired secondary structures, contrary to the classical expectations for non-coding RNAs in a hyperthermophile.

Conclusion

The results identify >2,700 TSS, including almost all of the primary sites of transcription initiation upstream of annotated genes, plus many secondary sites, sites within genes and sites resulting in antisense transcripts. The T. kodakarensis genome is small (~2.1 Mbp) and tightly packed with protein-encoding genes, but the transcriptomes established also contain many non-coding RNAs and predict extensive RNA-based regulation in this model Archaeon.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-684) contains supplementary material, which is available to authorized users.

Evolutionary patterns of Escherichia coli small RNAs and their regulatory interactions

Most bacterial small RNAs (sRNAs) are post-transcriptional regulators of gene expression, exerting their regulatory function by base-pairing with their target mRNAs. While it has become evident that sRNAs play central regulatory roles in the cell, little is known about their evolution and the evolution of their regulatory interactions. Here we used the prokaryotic phylogenetic tree to reconstruct the evolutionary history of Escherichia coli sRNAs and their binding sites on target mRNAs. We discovered that sRNAs currently present in E. coli mainly accumulated inside the Enterobacteriales order, succeeding the appearance of other types of noncoding RNAs and concurrently with the evolution of a variant of the Hfq protein exhibiting a longer C-terminal region. Our analysis of the evolutionary ages of sRNA-mRNA interactions revealed that while all sRNAs were evolutionarily older than most of their known binding sites on mRNA targets, for quite a few sRNAs there was at least one binding site that coappeared with or preceded them. It is conceivable that the establishment of these first interactions forced selective pressure on the sRNAs, after which additional targets were acquired by fitting a binding site to the active region of the sRNA. This conjecture is supported by the appearance of many binding sites on target mRNAs only after the sRNA gain, despite the prior presence of the target gene in ancestral genomes. Our results suggest a selective mechanism that maintained the sRNAs across the phylogenetic tree, and shed light on the evolution of E. coli post-transcriptional regulatory network.

The highly conserved bacterial RNase YbeY is essential in Vibrio cholerae, playing a critical role in virulence, stress regulation, and RNA processing

YbeY, a highly conserved protein, is an RNase in E. coli and plays key roles in both processing of the critical 3′ end of 16 S rRNA and in 70 S ribosome quality control under stress. These central roles account for YbeY's inclusion in the postulated minimal bacterial genome. However, YbeY is not essential in E. coli although loss of ybeY severely sensitizes it to multiple physiological stresses. Here, we show that YbeY is an essential endoribonuclease in Vibrio cholerae and is crucial for virulence, stress regulation, RNA processing and ribosome quality control, and is part of a core set of RNases essential in most representative pathogens. To understand its function, we analyzed the rRNA and ribosome profiles of a V. cholerae strain partially depleted for YbeY and other RNase mutants associated with 16 S rRNA processing; our results demonstrate that YbeY is also crucial for 16 S rRNA 3′ end maturation in V. cholerae and that its depletion impedes subunit assembly into 70 S ribosomes. YbeY's importance to V. cholerae pathogenesis was demonstrated by the complete loss of mice colonization and biofilm formation, reduced cholera toxin production, and altered expression levels of virulence-associated small RNAs of a V. cholerae strain partially depleted for YbeY. Notably, the ybeY genes of several distantly related pathogens can fully complement an E. coli ΔybeY strain under various stress conditions, demonstrating the high conservation of YbeY's activity in stress regulation. Taken together, this work provides the first comprehensive exploration of YbeY's physiological role in a human pathogen, showing its conserved function across species in essential cellular processes.

Bacteria adapt and survive unfavorable environments by quickly changing their gene expression and physiology, for example as pathogens do during infection of host cells. Gene expression is often determined by RNA turnover, a balance between transcription and RNA decay carried out by multiple RNases. The recently identified RNase YbeY was shown in E. coli to participate in rRNA maturation and 70 S ribosome quality control, however YbeY's roles in other organisms and the extent of functional conservation is unknown. Here, we show that YbeY is an essential RNase in the pathogen Vibrio cholerae, critical for cell fitness and general stress tolerance. We demonstrate that YbeY is crucial for 16 S rRNA 3′ end maturation, assembly of functional 70 S ribosomes and ribosome quality control. Moreover, YbeY regulates virulence-associated small RNAs and its depletion leads to an overall reduction in pathogenesis, exemplified by significantly decreased biofilm formation, mouse colonization and cholera toxin production. We also show that YbeY belongs to a minimal core set of RNases essential in most representative pathogens. The multifaceted roles of YbeY in several essential cellular processes and its highly conserved function across bacterial species, suggest that YbeY could be an attractive new antimicrobial target.

Small RNAs in the Vibrionaceae: an ocean still to be explored

In bacteria, the discovery of noncoding small RNAs (sRNAs) as modulators of gene expression in response to environmental signals has brought new insights into bacterial gene regulation, including control of pathogenicity. The Vibrionaceae constitute a family of marine bacteria of which many are responsible for infections affecting not only humans, such as Vibrio cholerae but also fish and marine invertebrates, representing the major cause of mortality in farmed marine species. They are able to colonize many habitats, existing as planktonic forms, in biofilms or associated with various hosts. This high adaptability is linked to their capacity to generate genetic diversity, in part through lateral gene transfer, but also by varying gene expression control. In the recent years, several major studies have illustrated the importance of small regulatory sRNAs in the Vibrionaceae for the control of pathogenicity and adaptation to environment and nutrient sources such as chitin, especially in V. cholerae and Vibrio harveyi. The existence of a complex regulatory network controlled by quorum sensing has been demonstrated in which sRNAs play central roles. This review covers major advances made in the discovery and elucidation of functions of Vibrionaceae sRNAs within the last 10 years.

Comprehensive discovery and characterization of small RNAs in Corynebacterium glutamicum ATCC 13032

Background

Recent discoveries on bacterial transcriptomes gave evidence that small RNAs (sRNAs) have important regulatory roles in prokaryotic cells. Modern high-throughput sequencing approaches (RNA-Seq) enable the most detailed view on transcriptomes offering an unmatched comprehensiveness and single-base resolution. Whole transcriptome data obtained by RNA-Seq can be used to detect and characterize all transcript species, including small RNAs. Here, we describe an RNA-Seq approach for comprehensive detection and characterization of small RNAs from Corynebacterium glutamicum, an actinobacterium of high industrial relevance and model organism for medically important Corynebacterianeae, such as C. diphtheriae and Mycobacterium tuberculosis.

Results

In our RNA-Seq approach, total RNA from C. glutamicum ATCC 13032 was prepared from cultures grown in minimal medium at exponential growth or challenged by physical (heat shock, cold shock) or by chemical stresses (diamide, H₂O₂, NaCl) at this time point. Total RNA samples were pooled and sequencing libraries were prepared from the isolated small RNA fraction. High throughput short read sequencing and mapping yielded over 800 sRNA genes. By determining their 5′- and 3′-ends and inspection of their locations, these potential sRNA genes were classified into UTRs of mRNAs (316), cis-antisense sRNAs (543), and trans-encoded sRNAs (262). For 77 of trans-encoded sRNAs significant sequence and secondary structure conservation was found by a computational approach using a whole genome alignment with the closely related species C. efficiens YS-314 and C. diphtheriae NCTC 13129. Three selected trans-encoded sRNAs were characterized by Northern blot analysis and stress-specific transcript patterns were found.

Conclusions

The study showed comparable numbers of sRNAs known from genome-wide surveys in other bacteria. In detail, our results give deep insight into the comprehensive equipment of sRNAs in C. glutamicum and provide a sound basis for further studies concerning the functions of these sRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-714) contains supplementary material, which is available to authorized users.

Comparative analysis of non-coding RNAs in the antibiotic-producing Streptomyces bacteria

Background

Non-coding RNAs (ncRNAs) are key regulatory elements that control a wide range of cellular processes in all bacteria in which they have been studied. Taking advantage of recent technological innovations, we set out to fully explore the ncRNA potential of the multicellular, antibiotic-producing Streptomyces bacteria.

Results

Using a comparative RNA sequencing analysis of three divergent model streptomycetes (S. coelicolor, S. avermitilis and S. venezuelae), we discovered hundreds of novel cis-antisense RNAs and intergenic small RNAs (sRNAs). We identified a ubiquitous antisense RNA species that arose from the overlapping transcription of convergently-oriented genes; we termed these RNA species ‘cutoRNAs’, for convergent untranslated overlapping RNAs. Conservation between different classes of ncRNAs varied greatly, with sRNAs being more conserved than antisense RNAs. Many species-specific ncRNAs, including many distinct cutoRNA pairs, were located within antibiotic biosynthetic clusters, including the actinorhodin, undecylprodigiosin, and coelimycin clusters of S. coelicolor, the chloramphenicol cluster of S. venezuelae, and the avermectin cluster of S. avermitilis.

Conclusions

These findings indicate that ncRNAs, including a novel class of antisense RNA, may exert a previously unrecognized level of regulatory control over antibiotic production in these bacteria. Collectively, this work has dramatically expanded the ncRNA repertoire of three Streptomyces species and has established a critical foundation from which to investigate ncRNA function in this medically and industrially important bacterial genus.

RNA at 92 °C: the non-coding transcriptome of the hyperthermophilic archaeon Pyrococcus abyssi

The non-coding transcriptome of the hyperthermophilic archaeon Pyrococcus abyssi is investigated using the RNA-seq technology. A dedicated computational pipeline analyzes RNA-seq reads and prior genome annotation to identify small RNAs, untranslated regions of mRNAs, and cis-encoded antisense transcripts. Unlike other archaea, such as Sulfolobus and Halobacteriales, P. abyssi produces few leaderless mRNA transcripts. Antisense transcription is widespread (215 transcripts) and targets protein-coding genes that appear to evolve more rapidly than average genes. We identify at least three novel H/ACA-like guide RNAs among the newly characterized non-coding RNAs. Long 5′ UTRs in mRNAs of ribosomal proteins and amino-acid biosynthesis genes strongly suggest the presence of cis-regulatory leaders in these mRNAs. We selected a high-interest subset of non-coding RNAs based on their strong promoters, high GC-content, phylogenetic conservation, or abundance. Some of the novel small RNAs and long 5′ UTRs display high GC contents, suggesting unknown structural RNA functions. However, we were surprised to observe that most of the high-interest RNAs are AU-rich, which suggests an absence of stable secondary structure in the high-temperature environment of P. abyssi. Yet, these transcripts display other hallmarks of functionality, such as high expression or high conservation, which leads us to consider possible RNA functions that do not require extensive secondary structure.

Detection of non-coding RNA in bacteria and archaea using the DETR'PROK Galaxy pipeline

RNA-seq experiments are now routinely used for the large scale sequencing of transcripts. In bacteria or archaea, such deep sequencing experiments typically produce 10-50 million fragments that cover most of the genome, including intergenic regions. In this context, the precise delineation of the non-coding elements is challenging. Non-coding elements include untranslated regions (UTRs) of mRNAs, independent small RNA genes (sRNAs) and transcripts produced from the antisense strand of genes (asRNA). Here we present a computational pipeline (DETR'PROK: detection of ncRNAs in prokaryotes) based on the Galaxy framework that takes as input a mapping of deep sequencing reads and performs successive steps of clustering, comparison with existing annotation and identification of transcribed non-coding fragments classified into putative 5' UTRs, sRNAs and asRNAs. We provide a step-by-step description of the protocol using real-life example data sets from Vibrio splendidus and Escherichia coli.