Question: Detection of mtDNA
0
gravatar for gianvitop95
8 months 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 8 months ago by Jeremy Leipzig18k • written 8 months 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 8 months ago by Michael Dondrup46k

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 8 months ago by gianvitop950
2
gravatar for ATpoint
8 months ago by
ATpoint19k
Germany
ATpoint19k 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 8 months ago • written 8 months ago by ATpoint19k

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

ADD REPLYlink written 8 months 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 8 months ago • written 8 months ago by ATpoint19k
1
gravatar for WouterDeCoster
8 months ago by
Belgium
WouterDeCoster40k wrote:

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

ADD COMMENTlink written 8 months ago by WouterDeCoster40k
0
gravatar for h.mon
8 months ago by
h.mon26k
Brazil
h.mon26k 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 8 months ago • written 8 months ago by h.mon26k
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: 537 users visited in the last hour