The Utility of Shallow RNA-Seq for Documenting Differential Gene Expression in Genes with High and Low Levels of Expression

Abstract

The sequencing depth necessary for documenting differential gene expression using RNA-Seq has been little explored outside of model systems. In particular, the depth required to analyze large-scale patterns of differential transcription factor expression is not known. The goal of the present study is to explore the effectiveness of shallow (relatively low read depth) RNA-Seq. We focus on two tissues in the honey bee: the sting gland and the digestive tract. The sting gland is an experimentally well-understood tissue that we use to benchmark the utility of this approach. We use the digestive tract to test the results obtained with the sting gland, and to conduct RNA-Seq between tissue types. Using a list of experimentally verified genes conferring tissue-specific functions in the sting gland, we show that relatively little read depth is necessary to identify them. We argue that this result should be broadly applicable, since genes important for tissue-specific functions often have robust expression patterns, and because we obtained similar results in our analysis of the digestive tract. Furthermore, we demonstrate that the differential expression of transcription factors, which are transcribed at low levels compared to other genes, can nevertheless often be determined using shallow RNA-Seq. Overall, we find over 150 differentially expressed transcription factors in our tissues at a read depth of only 12 million. This work shows the utility of low-depth sequencing for identifying genes important for tissue-specific functions. It also verifies the often-held belief that transcription factors show low levels of expression, while demonstrating that, in spite of this, they are frequently amenable to shallow RNA-Seq. Our findings should be of benefit to researchers using RNA-Seq in many different biological systems.

Virus and dsRNA-triggered transcriptional responses reveal key components of honey bee antiviral defense

“…Each pool was paired-end sequenced (2 × 100 nt) on a HiSeq. 2500 using a TruSeq SBS sequencing kit version 4, yielding ~12 million reads per sample, corresponding to at least 9.7 fold coverage (Supplementary Table S2 ), which is in the range of coverage reported in other honey bee transcriptome studies 24 , 25 , 122 . Sequence data was deposited into the NCBI Sequence Read Archive under accession number SRP101337 and is linked with NCBI BioProject #PRJNA377749.…”

An integrated RNAseq-1H NMR metabolomics approach to understand soybean primary metabolism regulation in response to Rhizoctonia foliar blight disease

“…Barcoding of samples permits pooling of multiple replicates allowing for the use of shallow-end sequencing to study a biological system. Due to the diminishment of the number of reads per sample, shallow-end sequencing can effectively be used to determine the fluctuations occurring in differentially expressed (DE) transcripts across a larger array of conditions than typically studied [ 29 , 30 ]. It has been suggested that analyzing the top 2500 transcripts can in fact cover over 80% of the biological information in an Arabidopsis transcriptomic project [ 30 ], and as little as 100,000 reads per samples are required for accurate prediction of mRNA fluctuations in human studies [ 31 ] suggesting that shallow-end RNAseq can provide substantial information.…”

Pathogen-Free PTI Induction with a Plant Conditioner

“…the degree of their really expression. From the RPM values we obtained the individual expression levels of the enzyme sequences reacting to the treatments applied at a given time, that is, the degree of their expression [106].…”