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

Work in recent years has led to the recognition of the importance of small regulatory RNAs (sRNAs) in bacterial regulation networks. New high-throughput sequencing technologies are paving the way to the exploration of an expanding sRNA world in nonmodel bacteria. In the Vibrio genus, compared to the enterobacteriaceae, still a limited number of sRNAs have been characterized, mostly in Vibrio cholerae, where they have been shown to be important for virulence, as well as in Vibrio harveyi. In addition, genome-wide approaches in V. cholerae have led to the discovery of hundreds of potential new sRNAs. Vibrio splendidus is an oyster pathogen that has been recently associated with massive mortality episodes in the French oyster growing industry. Here, we report the first RNA-seq study in a Vibrio outside of the V. cholerae species. We have uncovered hundreds of candidate regulatory RNAs, be it cis-regulatory elements, antisense RNAs, and trans-encoded sRNAs. Conservation studies showed the majority of them to be specific to V. splendidus. However, several novel sRNAs, previously unidentified, are also present in V. cholerae. Finally, we identified 28 trans sRNAs that are conserved in all the Vibrio genus species for which a complete genome sequence is available, possibly forming a Vibrio “sRNA core.”     

Expression of the ISPpu9 transposase of Pseudomonas putida KT2440 is regulated by two small RNAs and the secondary structure of the mRNA 5′-untranslated region

Insertion sequences (ISs) are mobile genetic elements that only carry the information required for their own transposition. Pseudomonas putida KT2440, a model bacterium, has seven copies of an IS called ISPpu9 inserted into repetitive extragenic palindromic sequences. This work shows that the gene for ISPpu9 transposase, tnp, is regulated by two small RNAs (sRNAs) named Asr9 and Ssr9, which are encoded upstream and downstream of tnp, respectively. The tnp mRNA has a long 5′-untranslated region (5′-UTR) that can fold into a secondary structure that likely includes the ribosome-binding site (RBS). Mutations weakening this structure increased tnp mRNA translation. Asr9, an antisense sRNA complementary to the 5′-UTR, was shown to be very stable. Eliminating Asr9 considerably reduced tnp mRNA translation, suggesting that it helps to unfold this secondary structure, exposing the RBS. Ectopic overproduction of Asr9 increased the transposition frequency of a new ISPpu9 entering the cell by conjugation, suggesting improved tnp expression. Ssr9 has significant complementarity to Asr9 and annealed to it in vitro forming an RNA duplex; this would sequester it and possibly facilitate its degradation. Thus, the antisense Asr9 sRNA likely facilitates tnp expression, improving transposition, while Ssr9 might counteract Asr9, keeping tnp expression low.     

The Arabidopsis RNA-Directed DNA Methylation Argonautes Functionally Diverge Based on Their Expression and Interaction with Target Loci

Argonaute (AGO) effectors of RNA silencing bind small RNA (sRNA) molecules and mediate mRNA cleavage, translational repression, or epigenetic DNA modification. In many organisms, these targeting mechanisms are devolved to different products of AGO multigene families. To investigate the basis of AGO functional diversification, we characterized three closely related Arabidopsis thaliana AGOs (AGO4, AGO6, and AGO9) implicated in RNA-directed DNA methylation. All three AGOs bound 5′ adenosine 24-nucleotide sRNAs, but each exhibited different preferences for sRNAs from different heterochromatin-associated loci. This difference was reduced when AGO6 and AGO9 were expressed from the AGO4 promoter, indicating that the functional diversification was partially due to differential expression of the corresponding genes. However, the AGO4-directed pattern of sRNA accumulation and DNA methylation was not fully recapitulated with AGO6 or AGO9 expressed from the AGO4 promoter. Here, we show that sRNA length and 5′ nucleotide do not account for the observed functional diversification of these AGOs. Instead, the selectivity of sRNA binding is determined by the coincident expression of the AGO and sRNA-generating loci, and epigenetic modification is influenced by interactions between the AGO protein and the different target loci. These findings highlight the importance of tissue specificity and AGO-associated proteins in influencing epigenetic modifications.     

Combined in silico and experimental identification of the Pyrococcus abyssi H/ACA sRNAs and their target sites in ribosomal RNAs

How far do H/ACA sRNPs contribute to rRNA pseudouridylation in Archaea was still an open question. Hence here, by computational search in three Pyrococcus genomes, we identified seven H/ACA sRNAs and predicted their target sites in rRNAs. In parallel, we experimentally identified 17 Ψ residues in P. abyssi rRNAs. By in vitro reconstitution of H/ACA sRNPs, we assigned 15 out of the 17 Ψ residues to the 7 identified H/ACA sRNAs: one H/ACA motif can guide up to three distinct pseudouridylations. Interestingly, by using a 23S rRNA fragment as the substrate, one of the two remaining Ψ residues could be formed in vitro by the aCBF5/aNOP10/aGAR1 complex without guide sRNA. Our results shed light on structural constraints in archaeal H/ACA sRNPs: the length of helix H2 is of 5 or 6 bps, the distance between the ANA motif and the targeted U residue is of 14 or 15 nts, and the stability of the interaction formed by the substrate rRNA and the 3′-guide sequence is more important than that formed with the 5′-guide sequence. Surprisingly, we showed that a sRNA–rRNA interaction with the targeted uridine in a single-stranded 5′-UNN-3′ trinucleotide instead of the canonical 5′-UN-3′ dinucleotide is functional.     

RNA-Seq of the Nucleolus Reveals Abundant SNORD44-Derived Small RNAs

Small non-coding RNAs represent RNA species that are not translated to proteins, but which have diverse and broad functional activities in physiological and pathophysiological states. The knowledge of these small RNAs is rapidly expanding in part through the use of massive parallel (deep) sequencing efforts. We present here the first deep sequencing of small RNomes in subcellular compartments with particular emphasis on small RNAs (sRNA) associated with the nucleolus. The vast majority of the cellular, cytoplasmic and nuclear sRNAs were identified as miRNAs. In contrast, the nucleolar sRNAs had a unique size distribution consisting of 19–20 and 25 nt RNAs, which were predominantly composed of small snoRNA-derived box C/D RNAs (termed as sdRNA). Sequences from 47 sdRNAs were identified, which mapped to both 5′ and 3′ ends of the snoRNAs, and retained conserved box C or D motifs. SdRNA reads mapping to SNORD44 comprised 74% of all nucleolar sdRNAs, and were confirmed by Northern blotting as comprising both 20 and 25 nt RNAs. A novel 120 nt SNORD44 form was also identified. The expression of the SNORD44 sdRNA and 120 nt form was independent of Dicer/Drosha–mediated processing pathways but was dependent on the box C/D snoRNP proteins/sno-ribonucleoproteins fibrillarin and NOP58. The 120 nt SNORD44-derived RNA bound to fibrillarin suggesting that C/D sno-ribonucleoproteins are involved in regulating the stability or processing of SNORD44. This study reveals sRNA cell-compartment specific expression and the distinctive unique composition of the nucleolar sRNAs.     

Natural Selection Shapes Maintenance of Orthologous sRNAs in Divergent Host-Restricted Bacterial Genomes

Historically it has been difficult to study the evolution of bacterial small RNAs (sRNAs) across distantly related species. For example, identifying homologs of sRNAs is often difficult in genomes that have undergone multiple structural rearrangements. Also, some types of regulatory sRNAs evolve at rapid rates. The high degree of genomic synteny among divergent host-restricted bacterial lineages, including intracellular symbionts, is conducive to sRNA maintenance and homolog identification. In turn, symbiont genomes can provide us with novel insights into sRNA evolution. Here, we examine the sRNA expression profile of the obligate symbiont of psyllids, Carsonella ruddii, which has one of the smallest cellular genomes described. Using RNA-seq, we identified 36 and 32 antisense sRNAs (asRNAs) expressed by Carsonella from the psyllids Bactericera cockerelli (Carsonella-BC) and Diaphorina citri (Carsonella-DC), respectively. The majority of these asRNAs were associated with genes that are involved in essential amino acid biosynthetic pathways. Eleven of the asRNAs were conserved in both Carsonella lineages and the majority were maintained by selection. Notably, five of the corresponding coding sequences are also the targets of conserved asRNAs in a distantly related insect symbiont, Buchnera. We detected differential expression of two asRNAs for genes involved in arginine and leucine biosynthesis occurring between two distinct Carsonella-BC life stages. Using asRNAs identified in Carsonella, Buchnera, and Profftella which are all endosymbionts, and Escherichia coli, we determined that regions upstream of these asRNAs encode unique conserved patterns of AT/GC richness, GC skew, and sequence motifs which may be involved in asRNA regulation.     

Hfq assists small RNAs in binding to the coding sequence of ompD mRNA and in rearranging its structure

The bacterial protein Hfq participates in the regulation of translation by small noncoding RNAs (sRNAs). Several mechanisms have been proposed to explain the role of Hfq in the regulation by sRNAs binding to the 5′-untranslated mRNA regions. However, it remains unknown how Hfq affects those sRNAs that target the coding sequence. Here, the contribution of Hfq to the annealing of three sRNAs, RybB, SdsR, and MicC, to the coding sequence of Salmonella ompD mRNA was investigated. Hfq bound to ompD mRNA with tight, subnanomolar affinity. Moreover, Hfq strongly accelerated the rates of annealing of RybB and MicC sRNAs to this mRNA, and it also had a small effect on the annealing of SdsR. The experiments using truncated RNAs revealed that the contributions of Hfq to the annealing of each sRNA were individually adjusted depending on the structures of interacting RNAs. In agreement with that, the mRNA structure probing revealed different structural contexts of each sRNA binding site. Additionally, the annealing of RybB and MicC sRNAs induced specific conformational changes in ompD mRNA consistent with local unfolding of mRNA secondary structure. Finally, the mutation analysis showed that the long AU-rich sequence in the 5′-untranslated mRNA region served as an Hfq binding site essential for the annealing of sRNAs to the coding sequence. Overall, the data showed that the functional specificity of Hfq in the annealing of each sRNA to the ompD mRNA coding sequence was determined by the sequence and structure of the interacting RNAs.     

SRD: a Staphylococcus regulatory RNA database

An overflow of regulatory RNAs (sRNAs) was identified in a wide range of bacteria. We designed and implemented a new resource for the hundreds of sRNAs identified in Staphylococci, with primary focus on the human pathogen Staphylococcus aureus. The “Staphylococcal Regulatory RNA Database” (SRD, http://srd.genouest.org/) compiled all published data in a single interface including genetic locations, sequences and other features. SRD proposes novel and simplified identifiers for Staphylococcal regulatory RNAs (srn) based on the sRNA''s genetic location in S. aureus strain N315 which served as a reference. From a set of 894 sequences and after an in-depth cleaning, SRD provides a list of 575 srn exempt of redundant sequences. For each sRNA, their experimental support(s) is provided, allowing the user to individually assess their validity and significance. RNA-seq analysis performed on strains N315, NCTC8325, and Newman allowed us to provide further details, upgrade the initial annotation, and identified 159 RNA-seq independent transcribed sRNAs. The lists of 575 and 159 sRNAs sequences were used to predict the number and location of srns in 18 S. aureus strains and 10 other Staphylococci. A comparison of the srn contents within 32 Staphylococcal genomes revealed a poor conservation between species. In addition, sRNA structure predictions obtained with MFold are accessible. A BLAST server and the intaRNA program, which is dedicated to target prediction, were implemented. SRD is the first sRNA database centered on a genus; it is a user-friendly and scalable device with the possibility to submit new sequences that should spread in the literature.     

Optimization of phage λ promoter strength for synthetic small regulatory RNA-based metabolic engineering

Synthetic small regulatory RNAs (sRNAs) are gene-silencing tools that can be used to tune gene expression in prokaryotes. A recent study by our group proposed rational design principles, introduced a regulatory system that may be used to implement synthetic sRNAs, and showed their utility in metabolic engineering. The regulatory system employed the strong phage λ P_R promoter to tightly control synthetic sRNA production. Here, we fine-tuned the strength of the P_R promoter via mutagenesis in order to optimize the level of synthetic sRNAs while maintaining the ability of the promoter to be regulated by CI proteins. Five mutant promoters of different strengths, ranging from 24 to 87% of that of the wild-type P_R promoter, were identified and confirmed to be repressed by CI proteins. A mutated promoter with only 40% of the original strength still produced enough synthetic sRNAs to inhibit the translation of the target mRNA to ~10% of the original level. As a practical application, we tested our promoters as drivers for a synthetic anti-murE sRNA, which was used to adjust the production of cadaverine. As the promoter strength decreased, the cadaverine titer first increased and then dropped. A mutated promoter with 39% of the original strength achieved the improved cadaverine titer of 2.15 g/L. The mutant promoters developed in this study should prove useful for tuning the expression levels of synthetic sRNAs for metabolic engineering.     

Deficiency of the tRNATyr:Ψ35-synthase aPus7 in Archaea of the Sulfolobales order might be rescued by the H/ACA sRNA-guided machinery

Up to now, Ψ formation in tRNAs was found to be catalysed by stand-alone enzymes. By computational analysis of archaeal genomes we detected putative H/ACA sRNAs, in four Sulfolobales species and in Aeropyrum pernix, that might guide Ψ35 formation in pre-tRNA^Tyr(GUA). This modification is achieved by Pus7p in eukarya. The validity of the computational predictions was verified by in vitro reconstitution of H/ACA sRNPs using the identified Sulfolobus solfataricus H/ACA sRNA. Comparison of Pus7-like enzymes encoded by archaeal genomes revealed amino acid substitutions in motifs IIIa and II in Sulfolobales and A. pernix Pus7-like enzymes. These conserved RNA:Ψ-synthase- motifs are essential for catalysis. As expected, the recombinant Pyrococcus abyssi aPus7 was fully active and acted at positions 35 and 13 and other positions in tRNAs, while the recombinant S. solfataricus aPus7 was only found to have a poor activity at position 13. We showed that the presence of an A residue 3′ to the target U residue is required for P. abyssi aPus7 activity, and that this is not the case for the reconstituted S. solfataricus H/ACA sRNP. In agreement with the possible formation of Ψ35 in tRNA^Tyr(GUA) by aPus7 in P. abyssi and by an H/ACA sRNP in S. solfataricus, the A36G mutation in the P. abyssi tRNA^Tyr(GUA) abolished Ψ35 formation when using P. abyssi extract, whereas the A36G substitution in the S. solfataricus pre-tRNA^Tyr did not affect Ψ35 formation in this RNA when using an S. solfataricus extract.