Question: relationship between coverage and expression level in RNAseq
1
gravatar for shl198
4.3 years ago by
shl198350
United States
shl198350 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.3 years ago by Gary450 • written 4.3 years ago by shl198350
2
gravatar for Diego D.
4.3 years ago by
Diego D.40
Chile
Diego D.40 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.3 years ago by Diego D.40
2
gravatar for Gary
4.3 years ago by
Gary450
Taiwan/Taichung/China Medical University Hospital
Gary450 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.3 years ago by Gary450
Please log in to add an answer.

Help
Access

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