Question: relationship between coverage and expression level in RNAseq
gravatar for shl198
4.9 years ago by
United States
shl198360 wrote:

I wonder what is the relationship between the depth and gene expression level in RNAseq?

For example, let's say I have RNAseq with 10M reads for one sample. If I use the same machine and generate 20M reads for the same sample, then except for the new detected genes, will the expression profile for genes expressed in both case be similar? In other words, if I plot the density VS log2(count) for them, will they have the similar shape?

ADD COMMENTlink modified 4.9 years ago by Gary480 • written 4.9 years ago by shl198360
gravatar for Diego D.
4.9 years ago by
Diego D.50
Diego D.50 wrote:

Yes, if you use the same method to normalize the counts, you will get similar MA plots (in theory). At greater coverage, maybe you will see more dots, because there will be more genes detected.

The reads generated during the sequencing should be randomly distributed, except if the process runs into a problem. After normalization, you should have similar values for genes that were detected in both cases.

More coverage is useful to detect low expressed genes.

You will see significant differences in MA plot if you change the number of biological replicates. With more replicates, there will be more differentially expressed genes.

ADD COMMENTlink written 4.9 years ago by Diego D.50
gravatar for Gary
4.9 years ago by
Taiwan/Taichung/China Medical University Hospital
Gary480 wrote:

Two related papers for your information.

Wang, Y., Ghaffari, N., Johnson, C.D., Braga-Neto, U.M., Wang, H., Chen, R. & Zhou, H. (2011). Evaluation of the coverage and depth of transcriptome by RNA-Seq in chickens. BMC Bioinformatics 12, S5.

Sims, D., Sudbery, I., Ilott, N.E., Heger, A. & Ponting, C.P. (2014). Sequencing depth and coverage: key considerations in genomic analyses. Nat. Rev. Genet. 15, 121-132.

ADD COMMENTlink written 4.9 years ago by Gary480
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1505 users visited in the last hour