Question: Detection of mtDNA
0
gravatar for gianvitop95
12 weeks ago by
gianvitop950 wrote:

Hi Community, I have a question about mtDNA. If I extract the whole cellular genome and then I sequence it with Illumina or other NGS tecniques, can I use a bioinformatic tool to extract from all the data output (such as file FASTA) just the mitocondrial DNA? Which tools should I use? Thank you in advance.

ADD COMMENTlink modified 12 weeks ago by Jeremy Leipzig18k • written 12 weeks ago by gianvitop950
2

Hi, your question is ambiguous. Please explain if it is just hypothetical or if you are dealing with real data and which.

What do you mean by:

all the data output (such as file FASTA) just the mitocondrial DNA

ADD REPLYlink written 12 weeks ago by Michael Dondrup45k

Hi, I'm sorry for the ambiguity. It is an hypothetical question. I was wondering if it is possible extract mtDNA sequences from a WGS method, like using Illumina. The output would be the sequences of the whole cellular genome (nuclear and mitocondrial) but what if I wanna extract just the mitocondrial data ? Is it possible? Thank you for the help.

ADD REPLYlink written 12 weeks ago by gianvitop950
2
gravatar for ATpoint
12 weeks ago by
ATpoint13k
Germany
ATpoint13k wrote:

Yes, technically this is possible from WGS data, but you'll not get much output. To put it in numbers, an average WGS file, like 30-50x with 500-1000mio fragments sequenced gives you a few hundred thousand reads mapping to chrM. The question would be what you want to do with it. If you want high coverage chrM reads, download any ATAC-seq experiment, as inherent in the library preparation method (especially with the original protocol from 2013), you'll have plenty of chrM fragments in your library. Given you have an indexed BAM file, simply run for extraction:

samtools view -bo chrM.bam total.bam chrM

Of course the chrM must be in the reference genome fasta file that you align against. What exactly do you want to do?

ADD COMMENTlink modified 12 weeks ago • written 12 weeks ago by ATpoint13k

Why would this be when mitochondria outnumber autosomal chromosomes in each cell by several thousand X?

ADD REPLYlink written 12 weeks ago by Jeremy Leipzig18k
1

You mean in WGS? The kits to extract genomic DNA typically do not capture the small circular plasmid-like mitochondrial genome. Edit: This should raise the question if the alignments to chrM are indeed mtDNA or some false-positives originating from mitochondrial homologues in the genome.

ADD REPLYlink modified 12 weeks ago • written 12 weeks ago by ATpoint13k
1
gravatar for WouterDeCoster
12 weeks ago by
Belgium
WouterDeCoster36k wrote:

What if you would just align your reads, and then go on with those aligned to the mitochondrial genome?

ADD COMMENTlink written 12 weeks ago by WouterDeCoster36k
0
gravatar for h.mon
12 weeks ago by
h.mon23k
Brazil
h.mon23k wrote:

There are many programs that aim at fishing mitochondrial reads (particularly from shotgun whole genome sequencing) and assembling them into a single contig, e.g. MITObim and Norgal. This is complicated by the fact there are "nuclear mitochondrial DNA segment", or NUMTs, spread over eukaryotic genomes - the WikiPedia page provides lots of references.

ADD COMMENTlink modified 12 weeks ago • written 12 weeks ago by h.mon23k
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: 1607 users visited in the last hour