Question: Chip-Seq For Mitochondrial Genome
0
gravatar for kanwarjag
7.2 years ago by
kanwarjag1.0k
United States
kanwarjag1.0k wrote:

I have a chipseq data where I am looking for mitochondrial genes differential binding. The samples are for human cell lines. My question is if I align reads with hg19 genome which as such will have mitochondrial genome will that give me a correct representation or I have to specifically take mitochondrial genome as reference for using any aligner

Thanks

ADD COMMENTlink modified 7.2 years ago by Sukhdeep Singh10k • written 7.2 years ago by kanwarjag1.0k
0
gravatar for Sukhdeep Singh
7.2 years ago by
Netherlands
Sukhdeep Singh10k wrote:

EIther way, you will get the information.

Map all reads and filter the Bed/Bam file for chrM or just subset the reference fasta index before mapping for chrM and then map and proceed with downstream analysis. The former will take more time, depending on the file size.

ADD COMMENTlink written 7.2 years ago by Sukhdeep Singh10k

mitochondrial DNA is circular, does some precautions/procedures need to be taken to deal with that? I assume almost all aligners take linear genome as reference?

ADD REPLYlink written 5.8 years ago by epigene490

Oh, I never thought about that, have a no idea. People at BWA/Bowtie support could answer.
There are papers for circular genome mappings though.
http://www.dcs.kcl.ac.uk/staff/csi/publications/ITAB_2010.pdf

ADD REPLYlink written 5.8 years ago by Sukhdeep Singh10k

Hmm, being discussed here, 2 strategies
1) Either add the beginning of the genome/chromosome to the end for the length of the reads, so you don't loose anything.
2) With Bowtie, use  -M 1 so that the read is assigned randomly. (quite not clear to me)
3) Use NovoAlign

http://www.researchgate.net/post/NGS_sequencing_issue_When_the_regerence_genome_is_circular

ADD REPLYlink written 5.8 years ago by Sukhdeep Singh10k
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: 1115 users visited in the last hour