A quantitative comparison of sRNA-based and protein-based gene regulation

Small, non-coding RNAs (sRNAs) play important roles as genetic regulators in prokaryotes. sRNAs act post-transcriptionally via complementary pairing with target mRNAs to regulate protein expression. We use a quantitative approach to compare and contrast sRNAs with conventional transcription factors (TFs) to better understand the advantages of each form of regulation. In particular, we calculate the steady-state behavior, noise properties, frequency-dependent gain (amplification), and dynamical response to large input signals of both forms of regulation. While the mean steady-state behavior of sRNA-regulated proteins exhibits a distinctive tunable threshold-linear behavior, our analysis shows that transcriptional bursting leads to significantly higher intrinsic noise in sRNA-based regulation than in TF-based regulation in a large range of expression levels and limits the ability of sRNAs to perform quantitative signaling. Nonetheless, we find that sRNAs are better than TFs at filtering noise in input signals. Additionally, we find that sRNAs allow cells to respond rapidly to large changes in input signals. These features suggest a niche for sRNAs in allowing cells to transition quickly yet reliably between distinct states. This functional niche is consistent with the widespread appearance of sRNAs in stress-response and quasi-developmental networks in prokaryotes.Comment: 26 pages, 8 figures; accepted for publication in Molecular Systems Biolog

An improved protocol for small RNA library construction using High Definition adapters

Next generation sequencing of small RNA (sRNA) libraries is widely used for studying sRNAs in various biological systems. However, cDNA libraries of sRNAs are biased for molecules that are ligated to adapters more or less efficiently than other molecules. One approach to reduce this ligation bias is to use a pool of adapters instead of a single adapter sequence, which allows many sRNAs to be ligated efficiently. We previously developed High Definition (HD) adapters for the Illumina sequencing platform, which contain degenerate nucleotides at the ligating ends of the adapters. However, the current commercial kits produced a large amount of 5’ adapter – 3’ adapter ligation product without the cDNA insert when HD adapters were used to replace the kit adapters. Here, we report a protocol to generate sRNA libraries using HD adapters with dramatically reduced adapter-adapter product. This protocol was developed from the procedure invented by Vaidyanathan et al. The libraries can be completed within two days and can be used for various biological and clinical samples. As examples for using this protocol, we constructed sRNA libraries using total RNA extracted from cultured mammalian cells and plant leaf tissue. The PCR products contained a very small amount of adapter-adapter product. Bioinformatic analysis of the sequencing data revealed sRNAs with diverse sequences and many different miRNA families

An improved method for surface immobilisation of RNA: application to small Non-Coding RNA - mRNA pairing

Characterisation of RNA and its intermolecular interactions is increasing in importance as the inventory of known RNA functions continues to expand. RNA-RNA interactions are central to post-transcriptional gene regulation mechanisms in bacteria, and the interactions of bacterial small non-coding RNAs (sRNAs) with their mRNA targets are the subject of much current research. The technology of surface plasmon resonance (SPR) is an attractive approach to studying these interactions since it is highly sensitive, and allows interaction measurements to be recorded in real-time. Whilst a number of approaches exist to label RNAs for surface-immobilisation, the method documented here is simple, quick, efficient, and utilises the high-affinity streptavidin-biotin interaction. Specifically, we ligate a biotinylated nucleotide to the 3' end of RNA using T4 RNA ligase. Although this is a previously recognised approach, we have optimised the method by our discovery that the incorporation of four or more adenine nucleotides at the 3' end of the RNA (a poly-A-tail) is required in order to achieve high ligation efficiencies. We use this method within the context of investigating small non-coding RNA (sRNA)-mRNA interactions through the application of surface technologies, including quantitative SPR assays. We first focus on validating the method using the recently characterised Escherichia coli sRNA-mRNA pair, MicA-ompA, specifically demonstrating that the addition of the poly-A-tail to either RNA does not affect its subsequent binding interactions with partner molecules. We then apply this method to investigate the novel interactions of a Vibrio cholerae Qrr sRNA with partner mRNAs, hapR and vca0939; RNA-RNA pairings that are important in mediating pathogenic virulence. The calculated binding parameters allow insights to be drawn regarding sRNA-mRNA interaction mechanisms

Changes in Aphid Host Plant Diet Influence the Small-RNA Expression Profiles of Its Obligate Nutritional Symbiont, Buchnera.

Plants are a difficult food resource to use, and herbivorous insects have evolved a variety of mechanisms that allow them to fully exploit this poor nutritional resource. One such mechanism is the maintenance of bacterial symbionts that aid in host plant feeding and development. The majority of these intracellular symbionts have highly eroded genomes that lack many key regulatory genes; consequently, it is unclear if these symbionts can respond to changes in the insect's diet to facilitate host plant use. There is emerging evidence that symbionts with highly eroded genomes express small RNAs (sRNAs), some of which potentially regulate gene expression. In this study, we sought to determine if the reduced genome of the nutritional symbiont (Buchnera) in the pea aphid responds to changes in the aphid's host plant diet. Using transcriptome sequencing (RNA-seq), Buchnera sRNA expression profiles were characterized within two Buchnera life stages when pea aphids fed on either alfalfa or fava bean. Overall, this study demonstrates that Buchnera sRNA expression changes not only with life stage but also with changes in aphid host plant diet. Of the 321 sRNAs characterized in this study, 47% were previously identified and 22% showed evidence of conservation in two or more Buchnera taxa. Functionally, 13 differentially expressed sRNAs were predicted to target genes related to pathways involved in essential amino acid biosynthesis. Overall, results from this study reveal that host plant diet influences the expression of conserved and lineage-specific sRNAs in Buchnera and that these sRNAs display distinct host plant-specific expression profiles among biological replicates.IMPORTANCE In general, the genomes of intracellular bacterial symbionts are reduced compared to those of free-living relatives and lack many key regulatory genes. Many of these reduced genomes belong to obligate mutualists of insects that feed on a diet that is deficient in essential nutrients, such as essential amino acids. It is unclear if these symbionts respond with their host to changes in insect diet, because of their reduced regulatory capacity. Emerging evidence suggests that these symbionts express small RNAs (sRNAs) that regulate gene expression at the posttranscriptional level. Therefore, in this study, we sought to determine if the reduced genome of the nutritional symbiont Buchnera in the pea aphid responds to changes in the aphid's host plant diet. This study demonstrates for the first time that Buchnera sRNAs, some conserved in two or more Buchnera lineages, are differentially expressed when aphids feed on different plant species and potentially target genes within essential amino acid biosynthesis pathways

On the Species Specificity of Acceptor RNA and Attachment Enzymes

One of the steps in protein biosynthesis appears to be the attachment of each amino acid to a specific acceptor (SRNA) molecule. According to the adaptor hypothesis, each SRNA molecule would then fit to a specific complementary base sequence on a linear RNA template, specifying the sequence of amino acids in the resultant protein [1,2]. An adaptor molecule thus could have two specificities: one recognizing the correct amino acid and activating enzyme; the other, the proper position on the template. The correctness of the amino-acid sequence therefore would depend upon the precision and constancy of the adaptors. However, the structures of the enzymes and adaptors are presumably under the genetic control of the organism and might be subject to heritable modifications. It is therefore conceivable that one or both ends of an adaptor might change sufficiently to cause occasional errors and, in the long run, an alteration of the genetic code might evolve. This notion, prompted by genetic observations [3] which suggested that mutation of a bacterium might modify its translation of genetic information, lead to the present comparison of the specificities of the acceptor RNA and activating enzymes of different organisms. Several differences in specificity have been reported previously. Berg et al. [4] demonstrated that SRNA from Escherichia coli contains two distinguishable acceptors for methionine. An enzyme prepared from yeast could attach methionine to one of these, while the enzyme from E. coli could attach to both. Webster found, in pig liver, a difference between the nuclear and cytoplasmic attachment enzymes for alanine. Rendi and Ochoa [6] noted that, for leucine, the enzymes in yeast and in E. coli could attach only to their homologous SRNA. Furthermore, in the case of leucine, rat liver enzyme and SRNA were interchangeable with those from E. coli. The observations presented below show that whether the enzymes and/or acceptors from two organisms are interchangeable depends upon not only the organisms in question but also the particular amino aci

Deep Sequencing Analysis of RNAs from Citrus Plants Grown in a Citrus Sudden Death-Affected Area Reveals Diverse Known and Putative Novel Viruses.

Citrus sudden death (CSD) has caused the death of approximately four million orange trees in a very important citrus region in Brazil. Although its etiology is still not completely clear, symptoms and distribution of affected plants indicate a viral disease. In a search for viruses associated with CSD, we have performed a comparative high-throughput sequencing analysis of the transcriptome and small RNAs from CSD-symptomatic and -asymptomatic plants using the Illumina platform. The data revealed mixed infections that included Citrus tristeza virus (CTV) as the most predominant virus, followed by the Citrus sudden death-associated virus (CSDaV), Citrus endogenous pararetrovirus (CitPRV) and two putative novel viruses tentatively named Citrus jingmen-like virus (CJLV), and Citrus virga-like virus (CVLV). The deep sequencing analyses were sensitive enough to differentiate two genotypes of both viruses previously associated with CSD-affected plants: CTV and CSDaV. Our data also showed a putative association of the CSD-symptomatic plants with a specific CSDaV genotype and a likely association with CitPRV as well, whereas the two putative novel viruses showed to be more associated with CSD-asymptomatic plants. This is the first high-throughput sequencing-based study of the viral sequences present in CSD-affected citrus plants, and generated valuable information for further CSD studies

Time Series Data Mining Algorithms for Identifying Short RNA in Arabidopsis thaliana

The class of molecules called short RNAs (sRNAs) are known to play a key role in gene regulation. Th are typically sequences of nucleotides between 21-25 nucleotides in length. They are known to play a key role in gene regulation. The identification, clustering and classification of sRNA has recently become the focus of much research activity. The basic problem involves detecting regions of interest on the chromosome where the pattern of candidate matches is somehow unusual. Currently, there are no published algorithms for detecting regions of interest, and the unpublished methods that we are aware of involve bespoke rule based systems designed for a specific organism. Work in this very new field has understandably focused on the outcomes rather than the methods used to obtain the results. In this paper we propose two generic approaches that place the specific biological problem in the wider context of time series data mining problems. Both methods are based on treating the occurrences on a chromosome, or “hit count” data, as a time series, then running a sliding window along a chromosome and measuring unusualness. This formulation means we can treat finding unusual areas of candidate RNA activity as a variety of time series anomaly detection problem. The first set of approaches is model based. We specify a null hypothesis distribution for not being a sRNA, then estimate the p-values along the chromosome. The second approach is instance based. We identify some typical shapes from known sRNA, then use dynamic time warping and fourier trans-form based distance to measure how closely the candidate series matches. We demonstrate that these methods can find known sRNA on Arabidopsis thaliana chromosomes and illustrate the benefits of the added information provided by these algorithms

Global analysis of the sugarcane microtranscriptome reveals a unique composition of small RNAs associated with axillary bud outgrowth

Axillary bud outgrowth determines shoot architecture and is under the control of endogenous hormones and a fine-tuned gene-expression network, which probably includes small RNAs (sRNAs). Although it is well known that sRNAs act broadly in plant development, our understanding about their roles in vegetative bud outgrowth remains limited. Moreover, the expression profiles of microRNAs (miRNAs) and their targets within axillary buds are largely unknown. Here, we employed sRNA next-generation sequencing as well as computational and gene-expression analysis to identify and quantify sRNAs and their targets in vegetative axillary buds of the biofuel crop sugarcane (Saccharum spp.). Computational analysis allowed the identification of 26 conserved miRNA families and two putative novel miRNAs, as well as a number of trans-acting small interfering RNAs. sRNAs associated with transposable elements and protein-encoding genes were similarly represented in both inactive and developing bud libraries. Conversely, sequencing and quantitative reverse transcription-PCR results revealed that specific miRNAs were differentially expressed in developing buds, and some correlated negatively with the expression of their targets at specific stages of axillary bud development. For instance, the expression patterns of miR159 and its target GAMYB suggested that they may play roles in regulating abscisic acid-signalling pathways during sugarcane bud outgrowth. Our work reveals, for the first time, differences in the composition and expression profiles of diverse sRNAs and targets between inactive and developing vegetative buds that, together with the endogenous balance of specific hormones, may be important in regulating axillary bud outgrowth

Discovery Of Ethanol-Responsive Small Rnas In Zymomonas Mobilis

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