cDNA Library Generation for the Analysis of Small RNAs by High-Throughput Sequencing

The RNome of a cell is highly diverse and consists besides messenger RNAs (mRNAs), transfer RNAs (tRNAs), and ribosomal RNAs (rRNAs) also of other small and long transcript entities without apparent coding potential. This class of molecules, commonly referred to as non-protein-coding RNAs (ncRNAs), is involved in regulating numerous biological processes and thought to contribute to cellular complexity. Therefore, much effort is put into their identification and further functional characterization. Here we provide a cost-effective and reliable method for cDNA library construction of small RNAs in the size range of 20-500 residues. The effectiveness of the described method is demonstrated by the analysis of ribosome-associated small RNAs in the eukaryotic model organism Trypanosoma brucei.

Length suppression in histone messenger RNA 3′-end maturation: Processing defects of insertion mutant premessenger RNAs can be compensated by insertions into the U7 small nuclear RNA

Efficient 3′-end processing of cell cycle-regulated mammalian histone premessenger RNAs (pre-mRNAs) requires an upstream stem–loop and a histone downstream element (HDE) that base pairs with the U7 small ribonuclearprotein. Insertions between these elements have two effects: the site of cleavage moves in concert with the HDE and processing efficiency declines. We used Xenopus oocytes to ask whether compensatory length insertions in the human U7 RNA could restore the fidelity and efficiency of processing of mouse histone insertion pre-mRNAs. An insertion of 5 nt into U7 RNA that extends its complementary to the HDE compensated for both defects in processing of a 5-nt insertion substrate; a noncomplementary insertion into U7 did not. Yet, the noncomplementary insertion mutant U7 was shown to be active on insertion substrates further mutated to allow base pairing. Our results suggest that the histone pre-mRNA becomes rigidified upstream of its HDE, allowing the bound U7 small ribonucleoprotein to measure from the HDE to the cleavage site. Such a mechanism may be common to other RNA measuring systems. To our knowledge, this is the first demonstration of length suppression in an RNA processing system.

Sense- and antisense-mediated gene silencing in tobacco is inhibited by the same viral suppressors and is associated with accumulation of small RNAs

Antisense-mediated gene silencing (ASGS) and posttranscriptional gene silencing (PTGS) with sense transgenes markedly reduce the steady-state mRNA levels of endogenous genes similar in transcribed sequence. RNase protection assays established that silencing in tobacco plants transformed with plant-defense-related class I sense and antisense chitinase (CHN) transgenes is at the posttranscriptional level. Infection of tobacco plants with cucumber mosaic virus strain FN and a necrotizing strain of potato virus Y, but not with potato virus X, effectively suppressed PTGS and ASGS of both the transgenes and homologous endogenes. This suggests that ASGS and PTGS share components associated with initiation and maintenance of the silent state. Small, ca. 25-nt RNAs (smRNA) of both polarities were associated with PTGS and ASGS in CHN transformants as reported for PTGS in other transgenic plants and for RNA interference in Drosophila. Similar results were obtained with an antisense class I β-1,3-glucanase transformant showing that viral suppression and smRNAs are a more general feature of ASGS. Several current models hold that diverse signals lead to production of double-stranded RNAs, which are processed to smRNAs that then trigger PTGS. Our results provide direct evidence for mechanistic links between ASGS and PTGS and suggest that ASGS could join a common PTGS pathway at the double-stranded RNA step.

Optimized methods for extracting circulating small RNAs from long-term stored equine samples.

Circulating miRNAs in body fluids, particularly serum, are promising candidates for future routine biomarker profiling in various pathologic conditions in human and veterinary medicine. However, reliable standardized methods for miRNA extraction from equine serum and fresh or archived whole blood are sorely lacking. We systematically compared various miRNA extraction methods from serum and whole blood after short and long-term storage without addition of RNA stabilizing additives prior to freezing. Time of storage at room temperature prior to freezing did not affect miRNA quality in serum. Furthermore, we showed that miRNA of NGS-sufficient quality can be recovered from blood samples after >10 years of storage at -80 °C. This allows retrospective analyses of miRNAs from archived samples.

Structure and dynamics of a gene network model incorporating small RNAs

As advances in molecular biology continue to reveal additional layers of complexity in gene regulation, computational models need to incorporate additional features to explore the implications of new theories and hypotheses. It has recently been suggested that eukaryotic organisms owe their phenotypic complexity and diversity to the exploitation of small RNAs as signalling molecules. Previous models of genetic systems are, for several reasons, inadequate to investigate this theory. In this study, we present an artificial genome model of genetic regulatory networks based upon previous work by Torsten Reil, and demonstrate how this model generates networks with biologically plausible structural and dynamic properties. We also extend the model to explore the implications of incorporating regulation by small RNA molecules in a gene network. We demonstrate how, using these signals, highly connected networks can display dynamics that are more stable than expected given their level of connectivity.

Phylogenetic relationships of Trypanosoma chelodina and Trypanosoma binneyi from Australian tortoises and platypuses inferred from small subunit rRNA analyses

Trypanosome infections are often difficult to detect by conventional microscopy and their pleomorphy often confounds differential diagnosis. Molecular techniques are now being used to diagnose infections and to determine phylogenetic relationships between species. Complete small subunit rRNA gene sequences were determined for isolates of Trypanosoma chelodina from the Brisbane River tortoise (Emydura signata), the saw-shelled tortoise (Elseya latisternum), and the eastern snake-necked tortoise (Chelodina longicollis) from southeast Queensland, Australia. Partial sequence data were also obtained for T. binneyi from a platypus (Ornithorhynchus anatinus) from Tasmania. Phylogenetic relationships between T. chelodina, T. binneyi and other species were examined by maximum parsimony and likelihood methods. The Australian tortoise and platypus trypanosomes did not exhibit any close phylogenetic relationships with those of mammals, reptiles or amphibians, but were closely related to each other, and to fish trypanosomes. This contra-indicates their co-evolution with their vertebrate hosts but does not exclude co-evolution with different groups of invertebrate vectors, notably insects and leeches.

Non-coding RNAs: the architects of eukaryotic complexity

Around 98% of all transcriptional output in humans is noncoding RNA. RNA-mediated gene regulation is widespread in higher eukaryotes and complex genetic phenomena like RNA interference, co-suppression, transgene silencing, imprinting, methylation, and possibly position-effect variegation and transvection, all involve intersecting pathways based on or connected to RNA signaling. I suggest that the central dogma is incomplete, and that intronic and other non-coding RNAs have evolved to comprise a second tier of gene expression in eukaryotes, which enables the integration and networking of complex suites of gene activity. Although proteins are the fundamental effectors of cellular function, the basis of eukaryotic complexity and phenotypic variation may lie primarily in a control architecture composed of a highly parallel system of trans-acting RNAs that relay state information required for the coordination and modulation of gene expression, via chromatin remodeling, RNA-DNA, RNA-RNA and RNA-protein interactions. This system has interesting and perhaps informative analogies with small world networks and dataflow computing.

Groundwater from infiltration galleries used for small public water supply systems: contamination with pesticides and endocrine disruptors

Infiltration galleries are among the oldest known means used for small public water fountains. Owing to its ancestral origin they are usually associated with high quality water. Thirty-one compounds, including pesticides and estrogens from different chemical families, were analysed in waters from infiltration galleries collected in Alto Douro Demarcated Wine region (North of Portugal). A total of twelve compounds were detected in the water samples. Nine of these compounds are described as presenting evidence or potential evidence of interfering with the hormone system of humans and wildlife. Although concentrations of the target analytes were relatively low, many of them below their limit of quantification, four compounds were above quantification limit and two of them even above the legal limit of 0.1 lg/L: dimethoate (30.38 ng/L), folpet (64.35 ng/L), terbuthylazine-desethyl (22.28 to 292.36 ng/L) and terbuthylazine (22.49 to 369.33 ng/L).

Structured RNAs in the ENCODE selected regions of the human genome.

Functional RNA structures play an important role both in the context of noncoding RNA transcripts as well as regulatory elements in mRNAs. Here we present a computational study to detect functional RNA structures within the ENCODE regions of the human genome. Since structural RNAs in general lack characteristic signals in primary sequence, comparative approaches evaluating evolutionary conservation of structures are most promising. We have used three recently introduced programs based on either phylogenetic-stochastic context-free grammar (EvoFold) or energy directed folding (RNAz and AlifoldZ), yielding several thousand candidate structures (corresponding to approximately 2.7% of the ENCODE regions). EvoFold has its highest sensitivity in highly conserved and relatively AU-rich regions, while RNAz favors slightly GC-rich regions, resulting in a relatively small overlap between methods. Comparison with the GENCODE annotation points to functional RNAs in all genomic contexts, with a slightly increased density in 3'-UTRs. While we estimate a significant false discovery rate of approximately 50%-70% many of the predictions can be further substantiated by additional criteria: 248 loci are predicted by both RNAz and EvoFold, and an additional 239 RNAz or EvoFold predictions are supported by the (more stringent) AlifoldZ algorithm. Five hundred seventy RNAz structure predictions fall into regions that show signs of selection pressure also on the sequence level (i.e., conserved elements). More than 700 predictions overlap with noncoding transcripts detected by oligonucleotide tiling arrays. One hundred seventy-five selected candidates were tested by RT-PCR in six tissues, and expression could be verified in 43 cases (24.6%).

Structured RNAs in the ENCODE selected regions of the human genome

Functional RNA structures play an important role both in the context of noncoding RNA transcripts as well as regulatory elements in mRNAs. Here we present a computational study to detect functional RNA structures within the ENCODE regions of the human genome. Since structural RNAs in general lack characteristic signals in primary sequence, comparative approaches evaluating evolutionary conservation of structures are most promising. We have used three recently introduced programs based on either phylogenetic–stochastic context-free grammar (EvoFold) or energy directed folding (RNAz and AlifoldZ), yielding several thousand candidate structures (corresponding to ∼2.7% of the ENCODE regions). EvoFold has its highest sensitivity in highly conserved and relatively AU-rich regions, while RNAz favors slightly GC-rich regions, resulting in a relatively small overlap between methods. Comparison with the GENCODE annotation points to functional RNAs in all genomic contexts, with a slightly increased density in 3′-UTRs. While we estimate a significant false discovery rate of ∼50%–70% many of the predictions can be further substantiated by additional criteria: 248 loci are predicted by both RNAz and EvoFold, and an additional 239 RNAz or EvoFold predictions are supported by the (more stringent) AlifoldZ algorithm. Five hundred seventy RNAz structure predictions fall into regions that show signs of selection pressure also on the sequence level (i.e., conserved elements). More than 700 predictions overlap with noncoding transcripts detected by oligonucleotide tiling arrays. One hundred seventy-five selected candidates were tested by RT-PCR in six tissues, and expression could be verified in 43 cases (24.6%).

GacA-controlled activation of promoters for small RNA genes in Pseudomonas fluorescens.

The Gac/Rsm signal transduction pathway positively regulates secondary metabolism, production of extracellular enzymes, and biocontrol properties of Pseudomonas fluorescens CHA0 via the expression of three noncoding small RNAs, termed RsmX, RsmY, and RsmZ. The architecture and function of the rsmY and rsmZ promoters were studied in vivo. A conserved palindromic upstream activating sequence (UAS) was found to be necessary but not sufficient for rsmY and rsmZ expression and for activation by the response regulator GacA. A poorly conserved linker region located between the UAS and the -10 promoter sequence was also essential for GacA-dependent rsmY and rsmZ expression, suggesting a need for auxiliary transcription factors. One such factor involved in the activation of the rsmZ promoter was identified as the PsrA protein, previously recognized as an activator of the rpoS gene and a repressor of fatty acid degradation. Furthermore, the integration host factor (IHF) protein was found to bind with high affinity to the rsmZ promoter region in vitro, suggesting that DNA bending contributes to the regulated expression of rsmZ. In an rsmXYZ triple mutant, the expression of rsmY and rsmZ was elevated above that found in the wild type. This negative feedback loop appears to involve the translational regulators RsmA and RsmE, whose activity is antagonized by RsmXYZ, and several hypothetical DNA-binding proteins. This highly complex network controls the expression of the three small RNAs in response to cell physiology and cell population densities.

CHARACTERIZATION OF NrsZ, A NOVEL SMALL REGULATORY NON-CODING RNA INVOLVED IN THE ADAPTATION OF PSEUDOMONAS AERUGINOSA PAO1

The opportunistic ubiquitous pathogen Pseudomonas aeruginosa strain PAOl is a versatile Gram-negative bacterium that has the extraordinary capacity to colonize a wide diversity of ecological niches and to cause severe and persistent infections in humans. To ensure an optimal coordination of the genes involved in nutrient utilization, this bacterium uses the NtrB/C and/or the CbrA/B two-component systems, to sense nutrients availability and to regulate in consequence the expression of genes involved in their uptake and catabolism. NtrB/C is specialized in nitrogen utilization, while the CbrA/B system is involved in both carbon and nitrogen utilization and both systems activate their target genes expression in concert with the alternative sigma factor RpoN. Moreover, the NtrB/C and CbrA/B two- component systems regulate the secondary metabolism of the bacterium, such as the production of virulence factors. In addition to the fine-tuning transcriptional regulation, P. aeruginosa can rapidly modulate its metabolism using small non-coding regulatory RNAs (sRNAs), which regulate gene expression at the post-transcriptional level by diverse and sophisticated mechanisms and contribute to the fast physiological adaptability of this bacterium. In our search for novel RpoN-dependent sRNAs modulating the nutritional adaptation of P. aeruginosa PAOl, we discovered NrsZ (Nitrogen regulated sRNA), a novel RpoN-dependent sRNA that is induced under nitrogen starvation by the NtrB/C two-component system. NrsZ has a unique architecture, formed of three similar stem-loop structures (SL I, II and II) separated by variant spacer sequences. Moreover, this sRNA is processed in short individual stem-loop molecules, by internal cleavage involving the endoribonuclease RNAse E. Concerning NrsZ functions in P. aeruginosa PAOl, this sRNA was shown to trigger the swarming motility and the rhamnolipid biosurfactants production. This regulation is due to the NrsZ-mediated activation of rhlA expression, a gene encoding for an enzyme essential for swarming motility and rhamnolipids production. Interestingly, the SL I structure of NrsZ ensures its regulatory function on rhlA expression, suggesting that the similar SLs are the functional units of this modular sRNA. However, the regulatory mechanism of action of NrsZ on rhlA expression activation remains unclear and is currently being investigated. Additionally, the NrsZ regulatory network was investigated by a transcriptome analysis, suggesting that numerous genes involved in both primary and secondary metabolism are regulated by this sRNA. To emphasize the importance of NrsZ, we investigated its conservation in other Pseudomonas species and demonstrated that NrsZ is conserved and expressed under nitrogen limitation in Pseudomonas protegens Pf-5, Pseudomonas putida KT2442, Pseudomonas entomophila L48 and Pseudomonas syringae pv. tomato DC3000, strains having different ecological features, suggesting an important role of NrsZ in the adaptation of Pseudomonads to nitrogen starvation. Interestingly the architecture of the different NrsZ homologs is similarly composed by SL structures and variant spacer sequences. However, the number of SL repetitions is not identical, and one to six SLs were predicted on the different NrsZ homologs. Moreover, NrsZ is processed in short molecules in all the strains, similarly to what was previously observed in P. aeruginosa PAOl, and the heterologous expression of the NrsZ homologs restored rhlA expression, swarming motility and rhamnolipids production in the P. aeruginosa NrsZ mutant. In many aspects, NrsZ is an atypical sRNA in the bacterial panorama. To our knowledge, NrsZ is the first described sRNA induced by the NtrB/C. Moreover, its unique modular architecture and its processing in similar short SL molecules suggest that NrsZ belongs to a novel family of bacterial sRNAs. -- L'agent pathogène opportuniste et ubiquitaire Pseudomonas aeruginosa souche PAOl est une bactérie Gram négative versatile ayant l'extraordinaire capacité de coloniser différentes niches écologiques et de causer des infections sévères et persistantes chez l'être humain. Afin d'assurer une coordination optimale des gènes impliqués dans l'utilisation de différents nutriments, cette bactérie se sert de systèmes à deux composants tel que NtrB/C et CbrA/B afin de détecter la disponibilité des ressources nutritives, puis de réguler en conséquence l'expression des gènes impliqués dans leur importation et leur catabolisme. Le système NtrB/C régule l'utilisation des sources d'azote alors que le système CbrA/B est impliqué à la fois dans l'utilisation des sources de carbone et d'azote. Ces deux systèmes activent l'expression de leurs gènes-cibles de concert avec le facteur sigma alternatif RpoN. En outre, NtrB/C et CbrA/B régulent aussi le métabolisme secondaire, contrôlant notamment la production d'importants facteurs de virulence. En plus de toutes ces régulations génétiques fines ayant lieu au niveau transcriptionnel, P. aeruginosa est aussi capable de moduler son métabolisme en se servant de petits ARNs régulateurs non-codants (ARNncs), qui régulent l'expression génétique à un niveau post- transcriptionnel par divers mécanismes sophistiqués et contribuent à rendre particulièrement rapide l'adaptation physiologique de cette bactérie. Au cours de nos recherches sur de nouveaux ARNncs dépendant du facteur sigma RpoN et impliqués dans l'adaptation nutritionnelle de P. aeruginosa PAOl, nous avons découvert NrsZ (Nitrogen regulated sRNA), un ARNnc induit par la cascade NtrB/C-RpoN en condition de carence en azote. NrsZ a une architecture unique, composée de trois structures en tige- boucle (TB I, II et III) hautement similaires et séparées par des « espaceurs » ayant des séquences variables. De plus, cet ARNnc est clivé en petits fragments correspondant au trois molécules en tige-boucle, par un processus de clivage interne impliquant l'endoribonucléase RNase E. Concernant les fonctions de NrsZ chez P. aeruginosa PAOl, cet ARNnc est capable d'induire la motilité de type « swarming » et la production de biosurfactants, nommés rhamnolipides. Cette régulation est due à l'activation par NrsZ de l'expression de rhlA, un gène essentiel pour la motilité de type swarming et pour la production de rhamnolipides. Étonnamment, la structure TB I est capable d'assurer à elle seule la fonction régulatrice de NrsZ sur l'expression de rhlA, suggérant que ces molécules TBs sont les unités fonctionnelles de cet ARNnc modulaire. Cependant, le mécanisme moléculaire par lequel NrsZ active l'expression de rhlA demeure à ce jour incertain et est actuellement à l'étude. En plus, le réseau de régulations médiées par NrsZ a été étudié par une analyse de transcriptome qui a indiqué que de nombreux gènes impliqués dans le métabolisme primaire ou secondaire seraient régulés par NrsZ. Pour accentuer l'importance de NrsZ, nous avons étudié sa conservation dans d'autres espèces de Pseudomonas. Ainsi, nous avons démontré que NrsZ est conservé et exprimé en situation de carence d'azote par les souches Pseudomonas protegens Pf-5, Pseudomonas putida KT2442, Pseudomonas entomophila L48, Pseudomonas syringae pv. tomato DC3000, quatre espèces ayant des caractéristiques écologiques très différentes, suggérant que NrsZ joue un rôle important dans l'adaptation du genre Pseudomonas envers la carence en azote. Chez toutes les souches étudiées, les différents homologues de NrsZ présentent une architecture similaire faite de TBs conservées et d'espaceurs. Cependant, le nombre de TBs n'est pas identique et peut varier de une à six copies selon la souche. Les différentes versions de NrsZ sont clivées en petites molécules dans ces quatre souches, comme il a été observé chez P. aeruginosa PAOl. De plus, l'expression hétérologue des différentes variantes de NrsZ est capable de restaurer l'expression de rhlA, la motilité swarming et la production de rhamnolipides dans une souche de P. aeruginosa dont nrsZ a été inactivé. Par bien des aspects, NrsZ est un ARNnc atypique dans le monde bactérien. À notre connaissance, NrsZ est le premier ARNnc décrit comme étant régulé par le système NtrB/C. De plus, son unique architecture modulaire et son clivage en petites molécules similaires suggèrent que NrsZ appartient à une nouvelle famille d'ARNncs bactériens.

The small RNA PhrS stimulates synthesis of the Pseudomonas aeruginosa quinolone signal.

Quorum sensing, a cell-to-cell communication system based on small signal molecules, is employed by the human pathogen Pseudomonas aeruginosa to regulate virulence and biofilm development. Moreover, regulation by small trans-encoded RNAs has become a focal issue in studies of virulence gene expression of bacterial pathogens. In this study, we have identified the small RNA PhrS as an activator of PqsR synthesis, one of the key quorum-sensing regulators in P. aeruginosa. Genetic studies revealed a novel mode of regulation by a sRNA, whereby PhrS uses a base-pairing mechanism to activate a short upstream open reading frame to which the pqsR gene is translationally coupled. Expression of phrS requires the oxygen-responsive regulator ANR. Thus, PhrS is the first bacterial sRNA that provides a regulatory link between oxygen availability and quorum sensing, which may impact on oxygen-limited growth in P. aeruginosa biofilms.

Genome-wide search reveals a novel GacA-regulated small RNA in Pseudomonas species.

BACKGROUND: Small RNAs (sRNAs) are widespread among bacteria and have diverse regulatory roles. Most of these sRNAs have been discovered by a combination of computational and experimental methods. In Pseudomonas aeruginosa, a ubiquitous Gram-negative bacterium and opportunistic human pathogen, the GacS/GacA two-component system positively controls the transcription of two sRNAs (RsmY, RsmZ), which are crucial for the expression of genes involved in virulence. In the biocontrol bacterium Pseudomonas fluorescens CHA0, three GacA-controlled sRNAs (RsmX, RsmY, RsmZ) regulate the response to oxidative stress and the expression of extracellular products including biocontrol factors. RsmX, RsmY and RsmZ contain multiple unpaired GGA motifs and control the expression of target mRNAs at the translational level, by sequestration of translational repressor proteins of the RsmA family. RESULTS: A combined computational and experimental approach enabled us to identify 14 intergenic regions encoding sRNAs in P. aeruginosa. Eight of these regions encode newly identified sRNAs. The intergenic region 1698 was found to specify a novel GacA-controlled sRNA termed RgsA. GacA regulation appeared to be indirect. In P. fluorescens CHA0, an RgsA homolog was also expressed under positive GacA control. This 120-nt sRNA contained a single GGA motif and, unlike RsmX, RsmY and RsmZ, was unable to derepress translation of the hcnA gene (involved in the biosynthesis of the biocontrol factor hydrogen cyanide), but contributed to the bacterium's resistance to hydrogen peroxide. In both P. aeruginosa and P. fluorescens the stress sigma factor RpoS was essential for RgsA expression. CONCLUSION: The discovery of an additional sRNA expressed under GacA control in two Pseudomonas species highlights the complexity of this global regulatory system and suggests that the mode of action of GacA control may be more elaborate than previously suspected. Our results also confirm that several GGA motifs are required in an sRNA for sequestration of the RsmA protein.

Small RNA-dependent expression of secondary metabolism is controlled by Krebs cycle function in Pseudomonas fluorescens.

Pseudomonas fluorescens CHA0, an antagonist of phytopathogenic fungi in the rhizosphere of crop plants, elaborates and excretes several secondary metabolites with antibiotic properties. Their synthesis depends on three small RNAs (RsmX, RsmY, and RsmZ), whose expression is positively controlled by the GacS-GacA two-component system at high cell population densities. To find regulatory links between primary and secondary metabolism in P. fluorescens and in the related species Pseudomonas aeruginosa, we searched for null mutations that affected central carbon metabolism as well as the expression of rsmY-gfp and rsmZ-gfp reporter constructs but without slowing down the growth rate in rich media. Mutation in the pycAB genes (for pyruvate carboxylase) led to down-regulation of rsmXYZ and secondary metabolism, whereas mutation in fumA (for a fumarase isoenzyme) resulted in up-regulation of the three small RNAs and secondary metabolism in the absence of detectable nutrient limitation. These effects required the GacS sensor kinase but not the accessory sensors RetS and LadS. An analysis of intracellular metabolites in P. fluorescens revealed a strong positive correlation between small RNA expression and the pools of 2-oxoglutarate, succinate, and fumarate. We conclude that Krebs cycle intermediates (already known to control GacA-dependent virulence factors in P. aeruginosa) exert a critical trigger function in secondary metabolism via the expression of GacA-dependent small RNAs.