This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Does TCGA data only use host reads to make pre-processed data (CNV, Mutations status, RNA-seq)?

I have a question regarding data available on TCGA, including copy number variation, mutation status and RNA-seq transcriptome profiling data. The following are the bioinformatics pipelines used for pre-processing.

  • Copy Number Variation: ASCAT3 pipeline
  • Mutation status: MuTect2 pipeline
  • RNA-seq: STAR pipeline

My question is, are these pre-processed data (not raw reads) processed from only host genomics/transcriptomics reads? If so, I am curious what kind of filter is used to filter out non-host reads (i.e. viral reads) from BAM files.

Your replies will be valuable in helping me understand the bioinformatics pipelines of data pre-processing. Thank you very much.

rna-seq copy-number-variation mutation-status tcga
• 1,183 views • link
updated by Zhenyu Zhang Zhenyu Zhang • written by yahn • yahn

1 answer

There are viral reads in sequencing runs but they don't align to the human reference. The alignment is the filter of which you speak. The reads are unmapped in the BAM files.

https://pmc.ncbi.nlm.nih.gov/articles/PMC5828528/

There are endogenous retroviruses (ERVs) in the human genome that are basically viral fossils but I don't think most viral sequences would align with them.

Some reference genomes include decoy sequences to intercept viral reads with significant homology to human genes such as the Epstein-Barr virus (EBV).

• 0 views • link
written by Jeremy Leipzig Jeremy Leipzig

Log in to answer this question.

More posts like this