Imperial College London
Portrait Picture

DrDerekHuntley

Faculty of Natural Sciences, Department of Life Sciences

Senior Teaching Fellow
 
 
 
//

Contact

 

+44 (0)20 7594 7478d.huntley

 
 
//

Location

 

239Sir Ernst Chain BuildingSouth Kensington Campus

//

Summary

 

Publications

Citation

BibTex format

@article{Maude:2019:10.3389/fcell.2019.00201,
author = {Maude, H and Davidson, M and Charitakis, N and Diaz, L and Bowers, W and Gradovich, E and Andrew, T and Huntley, D},
doi = {10.3389/fcell.2019.00201},
journal = {Frontiers in Cell and Developmental Biology},
title = {NUMT confounding biases mitochondrial heteroplasmy calls in favor of the reference allele},
url = {http://dx.doi.org/10.3389/fcell.2019.00201},
volume = {7},
year = {2019}
}
Download

RIS format (EndNote, RefMan)

TY  - JOUR
AB  - Homology between mitochondrial DNA (mtDNA) and nuclear DNA of mitochondrial origin (nuMTs) causes confounding when aligning short sequence reads to the reference human genome, as the true sequence origin cannot be determined. Using a systematic in silico approach, we here report the impact of all potential mitochondrial variants on alignment accuracy and variant calling. A total of 49,707 possible mutations were introduced across the 16,569bp reference mitochondrial genome (16,569 x 3 alternative alleles), one variant at-at-time. The resulting in silico fragmentation and alignment to the entire reference genome (GRCh38) revealed preferential mapping of mutated mitochondrial fragments to nuclear loci, as variants increased loci similarity to nuMTs, for a total of 807, 362 and 41 variants at 333, 144 and 27 positions when using 100bp, 150bp and 300bp single end fragments. We subsequently modelled these affected variants at 50% heteroplasmy and carried out variant calling, observing bias in the reported allele frequencies in favor of the reference allele. Four variants (chrM:6023A, chrM:4456T, chrM:5147A and chrM:7521A) including a possible hypertension factor, chrM:4456T, caused 100% loss of coverage at the mutated position (with all 100bp single-end fragments aligning to homologous, nuclear positions instead of chrM), rendering these variants undetectable when aligning to the entire reference genome. Furthermore, four mitochondrial variants reported to be pathogenic were found to cause significant loss of coverage and select Haplogroup-defining SNPs were shown to exacerbate the loss of coverage caused by surrounding variants. Increased fragment length and use of paired-end reads both improved alignment accuracy.
AU  - Maude,H
AU  - Davidson,M
AU  - Charitakis,N
AU  - Diaz,L
AU  - Bowers,W
AU  - Gradovich,E
AU  - Andrew,T
AU  - Huntley,D
DO  - 10.3389/fcell.2019.00201
PY  - 2019///
SN  - 2296-634X
TI  - NUMT confounding biases mitochondrial heteroplasmy calls in favor of the reference allele
T2  - Frontiers in Cell and Developmental Biology
UR  - http://dx.doi.org/10.3389/fcell.2019.00201
UR  - http://hdl.handle.net/10044/1/73167
VL  - 7
ER  -
Download