On the utility of RNA sample pooling to optimize cost and statistical power in RNA sequencing experiments

Background: In gene expression studies, RNA sample pooling is sometimes considered because of budget constraints or lack of sufficient input material. Using microarray technology, RNA sample pooling strategies have been reported to optimize both the cost of data generation as well as the statistical power for differential gene expression (DGE) analysis. For RNA sequencing, with its different quantitative output in terms of counts and tunable dynamic range, the adequacy and empirical validation of RNA sample pooling strategies have not yet been evaluated. In this study, we comprehensively assessed the utility of pooling strategies in RNA-seq experiments using empirical and simulated RNA-seq datasets. Result: The data generating model in pooled experiments is defined mathematically to evaluate the mean and variability of gene expression estimates. The model is further used to examine the trade-off between the statistical power of testing for DGE and the data generating costs. Empirical assessment of pooling strategies is done through analysis of RNA-seq datasets under various pooling and non-pooling experimental settings. Simulation study is also used to rank experimental scenarios with respect to the rate of false and true discoveries in DGE analysis. The results demonstrate that pooling strategies in RNA-seq studies can be both cost-effective and powerful when the number of pools, pool size and sequencing depth are optimally defined. Conclusion: For high within-group gene expression variability, small RNA sample pools are effective to reduce the variability and compensate for the loss of the number of replicates. Unlike the typical cost-saving strategies, such as reducing sequencing depth or number of RNA samples (replicates), an adequate pooling strategy is effective in maintaining the power of testing DGE for genes with low to medium abundance levels, along with a substantial reduction of the total cost of the experiment. In general, pooling RNA samples or pooling RNA samples in conjunction with moderate reduction of the sequencing depth can be good options to optimize the cost and maintain the power

A Novel Gonadotropin-Releasing Hormone 1 (Gnrh1) Enhancer-Derived Noncoding RNA Regulates Gnrh1 Gene Expression in GnRH Neuronal Cell Models.

Gonadotropin-releasing hormone (GnRH), a neuropeptide released from a small population of neurons in the hypothalamus, is the central mediator of the hypothalamic-pituitary-gonadal axis, and is required for normal reproductive development and function. Evolutionarily conserved regulatory elements in the mouse, rat, and human Gnrh1 gene include three enhancers and the proximal promoter, which confer Gnrh1 gene expression specifically in GnRH neurons. In immortalized mouse hypothalamic GnRH (GT1-7) neurons, which show pulsatile GnRH release in culture, RNA sequencing and RT-qPCR revealed that expression of a novel long noncoding RNA at Gnrh1 enhancer 1 correlates with high levels of GnRH mRNA expression. In GT1-7 neurons, which contain a transgene carrying 3 kb of the rat Gnrh1 regulatory region, both the mouse and rat Gnrh1 enhancer-derived noncoding RNAs (GnRH-E1 RNAs) are expressed. We investigated the characteristics and function of the endogenous mouse GnRH-E1 RNA. Strand-specific RT-PCR analysis of GnRH-E1 RNA in GT1-7 cells revealed GnRH-E1 RNAs that are transcribed in the sense and antisense directions from distinct 5' start sites, are 3' polyadenylated, and are over 2 kb in length. These RNAs are localized in the nucleus and have a half-life of over 8 hours. In GT1-7 neurons, siRNA knockdown of mouse GnRH-E1 RNA resulted in a significant decrease in the expression of the Gnrh1 primary transcript and Gnrh1 mRNA. Over-expression of either the sense or antisense mouse GnRH-E1 RNA in immature, migratory GnRH (GN11) neurons, which do not express either GnRH-E1 RNA or GnRH mRNA, induced the transcriptional activity of co-transfected rat Gnrh1 gene regulatory elements, where the induction requires the presence of the rat Gnrh1 promoter. Together, these data indicate that GnRH-E1 RNA is an inducer of Gnrh1 gene expression. GnRH-E1 RNA may play an important role in the development and maturation of GnRH neurons

RNA-seq transcriptional profiling of peripheral blood leukocytes from cattle infected with Mycobacterium bovis

Bovine tuberculosis, caused by infection with Mycobacterium bovis, is a major endemic disease affecting cattle populations worldwide, despite the implementation of stringent surveillance and control programs in many countries. The development of high-throughput functional genomics technologies, including gene expression microarrays and RNA-sequencing (RNA-seq), has enabled detailed analysis of the host transcriptome to M. bovis infection, particularly at the macrophage and peripheral blood level. In the present study, we have analyzed the peripheral blood leukocyte (PBL) transcriptome of eight natural M. bovis-infected and eight age- and sex-matched non-infected control Holstein-Friesian animals using RNA-seq. In addition, we compared gene expression profiles generated using RNA-seq with those previously generated using the high-density Affymetrix(®) GeneChip(®) Bovine Genome Array platform from the same PBL-extracted RNA. A total of 3,250 differentially expressed (DE) annotated genes were detected in the M. bovis-infected samples relative to the controls (adjusted P-value ≤0.05), with the number of genes displaying decreased relative expression (1,671) exceeding those with increased relative expression (1,579). Ingenuity(®) Systems Pathway Analysis (IPA) of all DE genes revealed enrichment for genes with immune function. Notably, transcriptional suppression was observed among several of the top-ranking canonical pathways including Leukocyte Extravasation Signaling. Comparative platform analysis demonstrated that RNA-seq detected a larger number of annotated DE genes (3,250) relative to the microarray (1,398), of which 917 genes were common to both technologies and displayed the same direction of expression. Finally, we show that RNA-seq had an increased dynamic range compared to the microarray for estimating differential gene expression

Transcriptional Regulation of Gene Expression in \u3cem\u3eTetrahymena thermophila\u3c/em\u3e

The only well-characterized study of gene expression in Tetrahymena thermophila (1) demonstrates that the temperature dependent expression of the Ser H3 gene is regulated at the level of mRNA stability. A run-on transcription assay was developed to determine if regulation of RNA stability was a major mechanism regulating gene expression in Tetrahymena or if transcriptional regulation dominates. The relative transcriptional activities of 14 Tetrahymena genes were determined in different physiological/developmental states (growing, starved and conjugating) in which many of the genes showed striking differences in RNA abundance. In every case except Ser H3, changes in transcription accompanied changes in RNA abundance. Thus differential transcription, not differential RNA degradation, is the major mechanism regulating RNA abundance in Tetrahymena

Broad application of a simple and affordable protocol for isolating plant RNA

BACKGROUND: Standard molecular biological methods involve the analysis of gene expression in living organisms under diverse environmental and developmental conditions. One of the most direct approaches to quantify gene expression is the isolation of RNA. Most techniques used to quantify gene expression require the isolation of RNA, usually from a large number of samples. While most published protocols, including those for commercial reagents, are either labour intensive, use hazardous chemicals and/or are costly, a previously published protocol for RNA isolation in Arabidopsis thaliana yields high amounts of good quality RNA in a simple, safe and inexpensive manner. FINDINGS: We have tested this protocol in tomato and wheat leaves, as well as in Arabidopsis leaves, and compared the resulting RNA to that obtained using a commercial phenol-based reagent. Our results demonstrate that this protocol is applicable to other plant species, including monocots, and offers yield and purity at least comparable to those provided by commercial phenol-based reagents. CONCLUSIONS: Here, we show that this previously published RNA isolation protocol can be easily extended to other plant species without further modification. Due to its simplicity and the use of inexpensive reagents, this protocol is accessible and affordable and can be easily implemented to work on different plant species in laboratories worldwide

CCR5AS lncRNA variation differentially regulates CCR5, influencing HIV disease outcome.

Multiple genome-wide studies have identified associations between outcome of human immunodeficiency virus (HIV) infection and polymorphisms in and around the gene encoding the HIV co-receptor CCR5, but the functional basis for the strongest of these associations, rs1015164A/G, is unknown. We found that rs1015164 marks variation in an activating transcription factor 1 binding site that controls expression of the antisense long noncoding RNA (lncRNA) CCR5AS. Knockdown or enhancement of CCR5AS expression resulted in a corresponding change in CCR5 expression on CD4+ T cells. CCR5AS interfered with interactions between the RNA-binding protein Raly and the CCR5 3' untranslated region, protecting CCR5 messenger RNA from Raly-mediated degradation. Reduction in CCR5 expression through inhibition of CCR5AS diminished infection of CD4+ T cells with CCR5-tropic HIV in vitro. These data represent a rare determination of the functional importance of a genome-wide disease association where expression of a lncRNA affects HIV infection and disease progression