Influence of RNA extraction methods and library selection schemes on RNA-seq data.

BACKGROUND: Gene expression analysis by RNA sequencing is now widely used in a number of applications surveying the whole transcriptomes of cells and tissues. The recent introduction of ribosomal RNA depletion protocols, such as RiboZero, has extended the view of the polyadenylated transcriptome to the poly(A)- fraction of the RNA. However, substantial amounts of intronic transcriptional activity has been reported in RiboZero protocols, raising issues regarding their potential nuclear origin and the impact on the actual sequence depth in exonic regions. RESULTS: Using HEK293 human cells as source material, we assessed here the impact of the two commonly used RNA extraction methods and of the library construction protocols (rRNA depletion versus mRNA) on 1) the relative abundance of intronic reads and 2) on the estimation of gene expression values. We benchmarked the rRNA depletion-based sequencing with a specific analysis of the cytoplasmic and nuclear transcriptome fractions, suggesting that the large majority of the intronic reads correspond to unprocessed nuclear transcripts rather than to independent transcriptional units. We show that Qiagen or TRIzol extraction methods retain differentially nuclear RNA species, and that consequently, rRNA depletion-based RNA sequencing protocols are particularly sensitive to the extraction methods. CONCLUSIONS: We could show that the combination of Trizol-based RNA extraction with rRNA depletion sequencing protocols led to the largest fraction of intronic reads, after the sequencing of the nuclear transcriptome. We discuss here the impact of the various strategies on gene expression and alternative splicing estimation measures. Further, we propose guidelines and a double selection strategy for minimizing the expression biases, without loss of information.

The BTB and CNC homology 1 (BACH1) target genes are involved in the oxidative stress response and in control of the cell cycle.

The regulation of gene expression in response to environmental signals and metabolic imbalances is a key step in maintaining cellular homeostasis. BTB and CNC homology 1 (BACH1) is a heme-binding transcription factor repressing the transcription from a subset of MAF recognition elements at low intracellular heme levels. Upon heme binding, BACH1 is released from the MAF recognition elements, resulting in increased expression of antioxidant response genes. To systematically address the gene regulatory networks involving BACH1, we combined chromatin immunoprecipitation sequencing analysis of BACH1 target genes in HEK 293 cells with knockdown of BACH1 using three independent types of small interfering RNAs followed by transcriptome profiling using microarrays. The 59 BACH1 target genes identified by chromatin immunoprecipitation sequencing were found highly enriched in genes showing expression changes after BACH1 knockdown, demonstrating the impact of BACH1 repression on transcription. In addition to known and new BACH1 targets involved in heme degradation (HMOX1, FTL, FTH1, ME1, and SLC48A1) and redox regulation (GCLC, GCLM, and SLC7A11), we also discovered BACH1 target genes affecting cell cycle and apoptosis pathways (ITPR2, CALM1, SQSTM1, TFE3, EWSR1, CDK6, BCL2L11, and MAFG) as well as subcellular transport processes (CLSTN1, PSAP, MAPT, and vault RNA). The newly identified impact of BACH1 on genes involved in neurodegenerative processes and proliferation provides an interesting basis for future dissection of BACH1-mediated gene repression in neurodegeneration and virus-induced cancerogenesis.

Prediction of alternative isoforms from exon expression levels in RNA-Seq experiments.

Alternative splicing, polyadenylation of pre-messenger RNA molecules and differential promoter usage can produce a variety of transcript isoforms whose respective expression levels are regulated in time and space, thus contributing specific biological functions. However, the repertoire of mammalian alternative transcripts and their regulation are still poorly understood. Second-generation sequencing is now opening unprecedented routes to address the analysis of entire transcriptomes. Here, we developed methods that allow the prediction and quantification of alternative isoforms derived solely from exon expression levels in RNA-Seq data. These are based on an explicit statistical model and enable the prediction of alternative isoforms within or between conditions using any known gene annotation, as well as the relative quantification of known transcript structures. Applying these methods to a human RNA-Seq dataset, we validated a significant fraction of the predictions by RT-PCR. Data further showed that these predictions correlated well with information originating from junction reads. A direct comparison with exon arrays indicated improved performances of RNA-Seq over microarrays in the prediction of skipped exons. Altogether, the set of methods presented here comprehensively addresses multiple aspects of alternative isoform analysis. The software is available as an open-source R-package called Solas at http://cmb.molgen.mpg.de/2ndGenerationSequencing/Solas/.

Gene expression variation in Down's syndrome mice allows prioritization of candidate genes.

BACKGROUND: Down's syndrome (DS), or trisomy 21, is a complex developmental disorder that exhibits many clinical signs that vary in occurrence and severity among patients. The molecular mechanisms responsible for DS have thus far remained elusive. We argue here that normal variation in gene expression in the population contributes to the heterogeneous clinical picture of DS, and we estimated the amplitude of this variation in 50 mouse orthologs of chromosome 21 genes in brain regions of Ts65Dn (a mouse model of DS). We analyzed the RNAs of eight Ts65Dn and eight euploid mice by real-time polymerase chain reaction. RESULTS: In pooled RNAs, we confirmed that trisomic/euploid gene expression ratios were close to 1.5. However, we observed that inter-individual gene expression levels spanned a broad range of values. We identified three categories of genes: genes with expression levels consistently higher in Ts65Dn than in euploids (9, 17, and 7 genes in cerebellum, cortex, and midbrain, respectively); genes whose expression levels partially overlap between the two groups (10, 9, and 14 genes); and genes with intermingled expression, which cannot be used to differentiate trisomics from euploids (12, 5 and 9 genes). Of the genes in the first category, App, Cbr1, and Mrps6 exhibited tight regulation in the three tissues and are therefore attractive candidates for further research. CONCLUSION: This is the first analysis addressing inter-individual gene expression levels as a function of trisomy. We propose a strategy allowing discrimination between candidates for the constant features of DS and those genes that may contribute to the partially penetrant signs of DS.

Ataxin-2 interacts with the DEAD/H-box RNA helicase DDX6 and interferes with P-bodies and stress granules.

Tight control of translation is fundamental for eukaryotic cells, and deregulation of proteins implicated contributes to numerous human diseases. The neurodegenerative disorder spinocerebellar ataxia type 2 is caused by a trinucleotide expansion in the SCA2 gene encoding a lengthened polyglutamine stretch in the gene product ataxin-2, which seems to be implicated in cellular RNA-processing pathways and translational regulation. Here, we substantiate a function of ataxin-2 in such pathways by demonstrating that ataxin-2 interacts with the DEAD/H-box RNA helicase DDX6, a component of P-bodies and stress granules, representing cellular structures of mRNA triage. We discovered that altered ataxin-2 levels interfere with the assembly of stress granules and cellular P-body structures. Moreover, ataxin-2 regulates the intracellular concentration of its interaction partner, the poly(A)-binding protein, another stress granule component and a key factor for translational control. Thus, our data imply that the cellular ataxin-2 concentration is important for the assembly of stress granules and P-bodies, which are main compartments for regulating and controlling mRNA degradation, stability, and translation.

Transcript level alterations reflect gene dosage effects across multiple tissues in a mouse model of down syndrome.

Human trisomy 21, which results in Down syndrome (DS), is one of the most complicated congenital genetic anomalies compatible with life, yet little is known about the molecular basis of DS. It is generally accepted that chromosome 21 (Chr21) transcripts are overexpressed by about 50% in cells with an extra copy of this chromosome. However, this assumption is difficult to test in humans due to limited access to tissues, and direct support for this idea is available for only a few Chr21 genes or in a limited number of tissues. The Ts65Dn mouse is widely used as a model for studies of DS because it is at dosage imbalance for the orthologs of about half of the 284 Chr21 genes. Ts65Dn mice have several features that directly parallel developmental anomalies of DS. Here we compared the expression of 136 mouse orthologs of Chr21 genes in nine tissues of the trisomic and euploid mice. Nearly all of the 77 genes which are at dosage imbalance in Ts65Dn showed increased transcript levels in the tested tissues, providing direct support for a simple model of increased transcription proportional to the gene copy number. However, several genes escaped this rule, suggesting that they may be controlled by additional tissue-specific regulatory mechanisms revealed in the trisomic situation.