Hi,
I have a question about genome alignment. I am working with RNA-Seq dataset to study the impact of Liquid Culture in response to virus of different doses in Human. I was exploring what could be the good strategy or best in practice method for genome mapping.
Maybe;
Identity the reads mapping to virus and exclude it from the analysis. After this, extract the unaligned reads and map against hg38 genome and perform quantification to count the genes followed by downstream analysis in either edgeR or Deseq2.
Additionally, I was thinking about the below scenarios:
HISAT or STAR aligner > then quantify using RSEM or FeatureCounts.(OR)
Bowtie/Bowtie2 and store the unmapped reads as fastq files, then use these unmapped fastq file in HISAT OR STAR to align against hg38> Quantify(OR)
BBSplit. Use this?Thank you very much for the help.
Toufiq
Maybe; Identity the reads mapping to virus and exclude it from the analysis. After this, extract the unaligned reads and map against hg38 genome
Not like this, As ATPoint said, you want to make a combined virus + human genome, and align to that. You do not want to ever align to a partial reference; you do not want the aligner to force reads to align to one thing when they would align better to something else. A combined reference will make sure that all the reads end up aligned to where they really belong.
RSEM is smarter about read counting than FeatureCounts, so all things being equal, use that. STAR's transcriptome output was designed to be compatible with RSEM, so I'd use that for aligning.
I feel you could use salmon and put the viral genome as decoy.
Rajendra KC Thank you for suggestions. Yes, I have this pipeline too. I have to first explore how to put the viral genome as decoy in the human genome. I have not done this earlier.
Thank you for the suggestions ATpoint I will try this way. Have you come across any example/resource on how to build combined index (viral genome into hg38) using STAR? This will be helpful to me in building.
swbarnes2 and ATpoint thank you very much. I will indeed try as suggested. I was earlier using the STAR hg38 index built by the collaborator. I will explore how to construct one which combined genome.
swbarnes2, ATpoint, and GenoMax. I figured out how to cat the fasta files together of both Human and Virus. The question now arises is about cat on gtf files. There's no gtf file for the viral organism I am working on. There are only fasta files for this virus. How to deal with this? Thank you.
Log in to answer this question.
Are the viral transcripts ending up in the final dataset? If so you may want to see if there is any correlation with the initial dosing. So while you could simply split and remove viral reads, doing what ATPoint suggest may be the way to go.