Contamination detection in sequencing studies using the mitochondrial phylogeny

mtDNA-Server 2: advancing mitochondrial DNA analysis through highly parallelized data processing and interactive analytics

Over the past decade, mtDNA-Server established itself as one of the most widely used variant calling web-services for human mitochondrial genomes. The service accepts sequencing data in BAM format and returns an annotated variant analysis report for both homoplasmic and heteroplasmic variants. In this work we present mtDNA-Server 2, which includes several new features highly requested by the community. Most importantly, it includes (a) the integration of a novel variant calling mode that accurately call insertions, deletions and single nucleotide variants at once, (b) the integration of additional quality control and input validation modules, (c) a method to estimate the required coverage to minimize false positives and (d) an interactive analytics dashboard. Furthermore, we migrated the complete analysis workflow to the Nextflow workflow manager for improved parallelization, reproducibility and local execution. Recognizing the importance of insertions and deletions as well as offering novel quality control, validation and reporting features, mtDNA-Server 2 provides researchers and clinicians a new state-of-the-art analysis platform for interpreting mitochondrial genomes. mtDNA-Server 2 is available via mitoverse, our analysis platform that offers a centralized place for mtDNA analysis in the cloud. The web-service, source code and its documentation are freely accessible at https://mitoverse.i-med.ac.at.

Late Neolithic collective burial reveals admixture dynamics during the third millennium BCE and the shaping of the European genome

The third millennium BCE was a pivotal period of profound cultural and genomic transformations in Europe associated with migrations from the Pontic-Caspian steppe, which shaped the ancestry patterns in the present-day European genome. We performed a high-resolution whole-genome analysis including haplotype phasing of seven individuals of a collective burial from ~2500 cal BCE and of a Bell Beaker individual from ~2300 cal BCE in the Paris Basin in France. The collective burial revealed the arrival in real time of steppe ancestry in France. We reconstructed the genome of an unsampled individual through its relatives' genomes, enabling us to shed light on the early-stage admixture patterns, dynamics, and propagation of steppe ancestry in Late Neolithic Europe. We identified two major Neolithic/steppe-related ancestry admixture pulses around 3000/2900 BCE and 2600 BCE. These pulses suggest different population expansion dynamics with striking links to the Corded Ware and Bell Beaker cultural complexes.

mtDNA analysis using Mitopore

Mitochondrial DNA (mtDNA) analysis is crucial for the diagnosis of mitochondrial disorders, forensic investigations, and basic research. Existing pipelines are complex, expensive, and require specialized personnel. In many cases, including the diagnosis of detrimental single nucleotide variants (SNVs), mtDNA analysis is still carried out using Sanger sequencing. Here, we developed a simple workflow and a publicly available webserver named Mitopore that allows the detection of mtDNA SNVs, indels, and haplogroups. To simplify mtDNA analysis, we tailored our workflow to process noisy long-read sequencing data for mtDNA analysis, focusing on sequence alignment and parameter optimization. We implemented Mitopore with eliBQ (eliminate bad quality reads), an innovative quality enhancement that permits the increase of per-base quality of over 20% for low-quality data. The whole Mitopore workflow and webserver were validated using patient-derived and induced pluripotent stem cells harboring mtDNA mutations. Mitopore streamlines mtDNA analysis as an easy-to-use fast, reliable, and cost-effective analysis method for both long- and short-read sequencing data. This significantly enhances the accessibility of mtDNA analysis and reduces the cost per sample, contributing to the progress of mtDNA-related research and diagnosis.

High number of mitochondrial DNA alterations in postmortem brain tissue of patients with schizophrenia compared to healthy controls

Previous studies have shown mitochondrial dysfunction in schizophrenia (SZ) patients, which may be caused by mitochondrial DNA (mtDNA) alterations. However, there are few studies in SZ that have analyzed mtDNA in brain samples by next-generation sequencing (NGS). To address this gap, we used mtDNA-targeted NGS and qPCR to characterize mtDNA alterations in brain samples from patients with SZ (n = 40) and healthy controls (HC) (n = 40). 35 % of SZ patients showed mtDNA alterations, a significantly higher prevalence compared to 10 % of HC. Specifically, SZ patients had a significantly higher frequency of deletions (35 vs. 5 in HC), with a mean number of deletions of 3.8 in SZ vs. 1.0 in HC. Likely pathogenic missense variants were also significantly more frequent in patients with SZ than in HC (10 vs. three HC), encompassing 14 variants in patients and three in HC. The pathogenic tRNA variant m.3243A>G was identified in one SZ patient with a high heteroplasmy level of 32.2 %. While no significant differences in mtDNA copy number (mtDNA-CN) were observed between SZ and HC, antipsychotic users had significantly higher mtDNA-CN than non-users. These findings suggest a potential role for mtDNA alterations in the pathophysiology of SZ that require further validation and functional studies.

MitoH3: Mitochondrial Haplogroup and Homoplasmic/Heteroplasmic Variant Calling Pipeline for Alzheimer's Disease Sequencing Project

Background

Mitochondrial DNA (mtDNA) is a double-stranded circular DNA and has multiple copies in each cell. Excess heteroplasmy, the coexistence of distinct variants in copies of mtDNA within a cell, may lead to mitochondrial impairments. Accurate determination of heteroplasmy in whole-genome sequencing (WGS) data has posed a significant challenge because mitochondria carrying heteroplasmic variants cannot be distinguished during library preparation. Moreover, sequencing errors, contamination, and nuclear mtDNA segments can reduce the accuracy of heteroplasmic variant calling.

Objective

To efficiently and accurately call mtDNA homoplasmic and heteroplasmic variants from the large-scale WGS data generated from the Alzheimer's Disease Sequencing Project (ADSP), and test their association with Alzheimer's disease (AD).

Methods

In this study, we present MitoH3-a comprehensive computational pipeline for calling mtDNA homoplasmic and heteroplasmic variants and inferring haplogroups in the ADSP WGS data. We first applied MitoH3 to 45 technical replicates from 6 subjects to define a threshold for detecting heteroplasmic variants. Then using the threshold of 5% ≤variant allele fraction≤95%, we further applied MitoH3 to call heteroplasmic variants from a total of 16,113 DNA samples with 6,742 samples from cognitively normal controls and 6,183 from AD cases.

Results

This pipeline is available through the Singularity container engine. For 4,311 heteroplasmic variants identified from 16,113 samples, no significant variant count difference was observed between AD cases and controls.

Conclusions

Our streamlined pipeline, MitoH3, enables computationally efficient and accurate analysis of a large number of samples.

CD4+ T cell mitochondrial genotype in Multiple Sclerosis: a cross-sectional and longitudinal analysis

Multiple Sclerosis (MS) is a chronic autoimmune demyelinating disease of the central nervous system (CNS), with a largely unknown etiology, where mitochondrial dysfunction likely contributes to neuroaxonal loss and brain atrophy. Mirroring the CNS, peripheral immune cells from patients with MS, particularly CD4+ T cells, show inappropriate mitochondrial phenotypes and/or oxidative phosphorylation (OxPhos) insufficiency, with a still unknown contribution of mitochondrial DNA (mtDNA). We hypothesized that mitochondrial genotype in CD4+ T cells might influence MS disease activity and progression. Thus, we performed a retrospective cross-sectional and longitudinal study on patients with a recent diagnosis of either Clinically Isolated Syndrome (CIS) or Relapsing-Remitting MS (RRMS) at two timepoints: 6 months (VIS1) and 36 months (VIS2) after disease onset. Our primary outcomes were the differences in mtDNA extracted from CD4+ T cells between: (I) patients with CIS/RRMS (PwMS) at VIS1 and age- and sex-matched healthy controls (HC), in the cross-sectional analysis, and (II) different diagnostic evolutions in PwMS from VIS1 to VIS2, in the longitudinal analysis. We successfully performed mtDNA whole genome sequencing (mean coverage: 2055.77 reads/base pair) in 183 samples (61 triplets). Nonetheless, mitochondrial genotype was not associated with a diagnosis of CIS/RRMS, nor with longitudinal diagnostic evolution.

Mitochondrial DNA integrity and metabolome profile are preserved in the human induced pluripotent stem cell reference line KOLF2.1J

Quality control of human induced pluripotent stem cells (iPSCs) is critical to ensure reproducibility of research. Recently, KOLF2.1J was characterized and published as a male iPSC reference line to study neurological disorders. Emerging evidence suggests potential negative effects of mtDNA mutations, but its integrity was not analyzed in the original publication. To assess mtDNA integrity, we conducted a targeted mtDNA analysis followed by untargeted metabolomics analysis. We found that KOLF2.1J mtDNA integrity was intact at the time of publication and is still preserved in the commercially distributed cell line. In addition, the basal KOLF2.1J metabolome profile was similar to that of the two commercially available iPSC lines IMR90 and iPSC12, but clearly distinct from an in-house-generated ERCC6R683X/R683X iPSC line modeling Cockayne syndrome. Conclusively, we validate KOLF2.1J as a reference iPSC line, and encourage scientists to conduct mtDNA analysis and unbiased metabolomics whenever feasible.

Genome sequences of 36,000- to 37,000-year-old modern humans at Buran-Kaya III in Crimea

Populations genetically related to present-day Europeans first appeared in Europe at some point after 38,000-40,000 years ago, following a cold period of severe climatic disruption. These new migrants would eventually replace the pre-existing modern human ancestries in Europe, but initial interactions between these groups are unclear due to the lack of genomic evidence from the earliest periods of the migration. Here we describe the genomes of two 36,000-37,000-year-old individuals from Buran-Kaya III in Crimea as belonging to this newer migration. Both genomes share the highest similarity to Gravettian-associated individuals found several thousand years later in southwestern Europe. These genomes also revealed that the population turnover in Europe after 40,000 years ago was accompanied by admixture with pre-existing modern human populations. European ancestry before 40,000 years ago persisted not only at Buran-Kaya III but is also found in later Gravettian-associated populations of western Europe and Mesolithic Caucasus populations.

North and East African mitochondrial genetic variation needs further characterization towards precision medicine

INTRODUCTION

Mitochondria are maternally inherited cell organelles with their own genome, and perform various functions in eukaryotic cells such as energy production and cellular homeostasis. Due to their inheritance and manifold biological roles in health and disease, mitochondrial genetics serves a dual purpose of tracing the history as well as disease susceptibility of human populations across the globe. This work requires a comprehensive catalogue of commonly observed genetic variations in the mitochondrial DNAs for all regions throughout the world. So far, however, certain regions, such as North and East Africa have been understudied.

OBJECTIVES

To address this shortcoming, we have created the most comprehensive quality-controlled North and East African mitochondrial data set to date and use it for characterizing mitochondrial genetic variation in this region.

METHODS

We compiled 11 published cohorts with novel data for mitochondrial genomes from 159 Sudanese individuals. We combined these 641 mitochondrial sequences with sequences from the 1000 Genomes (n = 2504) and the Human Genome Diversity Project (n = 828) and used the tool haplocheck for extensive quality control and detection of in-sample contamination, as well as Nanopore long read sequencing for haplogroup validation of 18 samples.

RESULTS

Using a subset of high-coverage mitochondrial sequences, we predict 15 potentially novel haplogroups in North and East African subjects and observe likely phylogenetic deviations from the established PhyloTree reference for haplogroups L0a1 and L2a1.

CONCLUSION

Our findings demonstrate common hitherto unexplored variants in mitochondrial genomes of North and East Africa that lead to novel phylogenetic relationships between haplogroups present in these regions. These observations call for further in-depth population genetic studies in that region to enable the prospective use of mitochondrial genetic variation for precision medicine.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans

Mitochondrial DNA (mtDNA) is a maternally inherited, high-copy-number genome required for oxidative phosphorylation1. Heteroplasmy refers to the presence of a mixture of mtDNA alleles in an individual and has been associated with disease and ageing. Mechanisms underlying common variation in human heteroplasmy, and the influence of the nuclear genome on this variation, remain insufficiently explored. Here we quantify mtDNA copy number (mtCN) and heteroplasmy using blood-derived whole-genome sequences from 274,832 individuals and perform genome-wide association studies to identify associated nuclear loci. Following blood cell composition correction, we find that mtCN declines linearly with age and is associated with variants at 92 nuclear loci. We observe that nearly everyone harbours heteroplasmic mtDNA variants obeying two principles: (1) heteroplasmic single nucleotide variants tend to arise somatically and accumulate sharply after the age of 70 years, whereas (2) heteroplasmic indels are maternally inherited as mixtures with relative levels associated with 42 nuclear loci involved in mtDNA replication, maintenance and novel pathways. These loci may act by conferring a replicative advantage to certain mtDNA alleles. As an illustrative example, we identify a length variant carried by more than 50% of humans at position chrM:302 within a G-quadruplex previously proposed to mediate mtDNA transcription/replication switching2,3. We find that this variant exerts cis-acting genetic control over mtDNA abundance and is itself associated in-trans with nuclear loci encoding machinery for this regulatory switch. Our study suggests that common variation in the nuclear genome can shape variation in mtCN and heteroplasmy dynamics across the human population.

Mitochondrial DNA in Human Diversity and Health: From the Golden Age to the Omics Era

Mitochondrial DNA (mtDNA) is a small fraction of our hereditary material. However, this molecule has had an overwhelming presence in scientific research for decades until the arrival of high-throughput studies. Several appealing properties justify the application of mtDNA to understand how human populations are-from a genetic perspective-and how individuals exhibit phenotypes of biomedical importance. Here, I review the basics of mitochondrial studies with a focus on the dawn of the field, analysis methods and the connection between two sides of mitochondrial genetics: anthropological and biomedical. The particularities of mtDNA, with respect to inheritance pattern, evolutionary rate and dependence on the nuclear genome, explain the challenges of associating mtDNA composition and diseases. Finally, I consider the relevance of this single locus in the context of omics research. The present work may serve as a tribute to a tool that has provided important insights into the past and present of humankind.

Next-generation sequencing reveals mitogenome diversity in plasma extracellular vesicles from colorectal cancer patients

BACKGROUND

Recent reports have demonstrated that the entire mitochondrial genome can be secreted in extracellular vesicles (EVs), but the biological attributes of this cell-free mitochondrial DNA (mtDNA) remain insufficiently understood. We used next-generation sequencing to compare plasma EV-derived mtDNA to that of whole blood (WB), peripheral blood mononuclear cells (PBMCs), and formalin-fixed paraffin-embedded (FFPE) tumor tissue from eight rectal cancer patients and WB and fresh-frozen (FF) tumor tissue from eight colon cancer patients.

METHODS

Total DNA was isolated before the mtDNA was enriched by PCR with either two primer sets generating two long products or multiple primer sets (for the FFPE tumors), prior to the sequencing. mtDNA diversity was assessed as the total variant number, level of heteroplasmy (mutant mtDNA copies mixed with wild-type copies), variant distribution within the protein-coding genes, and the predicted functional effect of the variants in the different sample types. Differences between groups were compared by paired Student's t-test or ANOVA with Dunnett's multiple comparison tests when comparing matched samples from patients. Mann-Whitney U test was used when comparing differences between the cancer types and patient groups. Pearson correlation analysis was performed.

RESULTS

In both cancer types, EV mtDNA presented twice as many variants and had significantly more low-level heteroplasmy than WB mtDNA. The EV mtDNA variants were clustered in the coding regions, and the proportion of EV mtDNA variants that were missense mutations (i.e., estimated to moderately affect the mitochondrial protein function) was significantly higher than in WB and tumor tissues. Nonsense mutations (i.e., estimated to highly affect the mitochondrial protein function) were only observed in the tumor tissues and EVs.

CONCLUSION

Taken together, plasma EV mtDNA in CRC patients exhibits a high degree of diversity.

TRIAL REGISTRATION

ClinicalTrials.gov: NCT01816607 . Registered 22 March 2013.

Population analysis of complete mitogenomes for 334 samples from El Salvador

The use of mitochondrial DNA (mtDNA) in the field of forensic genetics is widely spread mainly due to its advantages when identifying highly degraded samples. In this sense, massive parallel sequencing has made the analysis of the whole mitogenome more accessible, noticeably increasing the informativeness of mtDNA haplotypes. The civil war (1980-1992) in El Salvador caused many deaths and disappearances (including children) all across the country and the economic and social instability after the war forced many people to emigration. For this reason, different organizations have collected DNA samples from relatives with the aim of identifying missing people. Thus, we present a dataset containing 334 complete mitogenomes from the Salvadoran general population. To the best of our knowledge, this is the first publication of a nationwide forensic-quality complete mitogenome database of any Latin American country. We found 293 different haplotypes, with a random match probability of 0.0041 and 26.6 mean pairwise differences, which is similar to other Latin American populations, and which represent a marked improvement from the values obtained with just control region sequences. These haplotypes belong to 54 different haplogroups, being 91% of them of Native American origin. Over a third (35.9%) of the individuals carried at least a heteroplasmic site (excluding length heteroplasmies). Ultimately, the present database aims to represent mtDNA haplotype diversity in the general Salvadoran populations as a basis for the identification of people that disappeared during or after the civil war.

HaploCart: Human mtDNA haplogroup classification using a pangenomic reference graph human mtDNA haplogroup inference

Current mitochondrial DNA (mtDNA) haplogroup classification tools map reads to a single reference genome and perform inference based on the detected mutations to this reference. This approach biases haplogroup assignments towards the reference and prohibits accurate calculations of the uncertainty in assignment. We present HaploCart, a probabilistic mtDNA haplogroup classifier which uses a pangenomic reference graph framework together with principles of Bayesian inference. We demonstrate that our approach significantly outperforms available tools by being more robust to lower coverage or incomplete consensus sequences and producing phylogenetically-aware confidence scores that are unbiased towards any haplogroup. HaploCart is available both as a command-line tool and through a user-friendly web interface. The C++ program accepts as input consensus FASTA, FASTQ, or GAM files, and outputs a text file with the haplogroup assignments of the samples along with the level of confidence in the assignments. Our work considerably reduces the amount of data required to obtain a confident mitochondrial haplogroup assignment.

Nanopore sequencing identifies a higher frequency and expanded spectrum of mitochondrial DNA deletion mutations in human aging

Mitochondrial DNA (mtDNA) deletion mutations cause many human diseases and are linked to age-induced mitochondrial dysfunction. Mapping the mutation spectrum and quantifying mtDNA deletion mutation frequency is challenging with next-generation sequencing methods. We hypothesized that long-read sequencing of human mtDNA across the lifespan would detect a broader spectrum of mtDNA rearrangements and provide a more accurate measurement of their frequency. We employed nanopore Cas9-targeted sequencing (nCATS) to map and quantitate mtDNA deletion mutations and develop analyses that are fit-for-purpose. We analyzed total DNA from vastus lateralis muscle in 15 males ranging from 20 to 81 years of age and substantia nigra from three 20-year-old and three 79-year-old men. We found that mtDNA deletion mutations detected by nCATS increased exponentially with age and mapped to a wider region of the mitochondrial genome than previously reported. Using simulated data, we observed that large deletions are often reported as chimeric alignments. To address this, we developed two algorithms for deletion identification which yield consistent deletion mapping and identify both previously reported and novel mtDNA deletion breakpoints. The identified mtDNA deletion frequency measured by nCATS correlates strongly with chronological age and predicts the deletion frequency as measured by digital PCR approaches. In substantia nigra, we observed a similar frequency of age-related mtDNA deletions to those observed in muscle samples, but noted a distinct spectrum of deletion breakpoints. NCATS-mtDNA sequencing allows the identification of mtDNA deletions on a single-molecule level, characterizing the strong relationship between mtDNA deletion frequency and chronological aging.

Nuclear genetic control of mtDNA copy number and heteroplasmy in humans

Human mitochondria contain a high copy number, maternally transmitted genome (mtDNA) that encodes 13 proteins required for oxidative phosphorylation. Heteroplasmy arises when multiple mtDNA variants co-exist in an individual and can exhibit complex dynamics in disease and in aging. As all proteins involved in mtDNA replication and maintenance are nuclear-encoded, heteroplasmy levels can, in principle, be under nuclear genetic control, however this has never been shown in humans. Here, we develop algorithms to quantify mtDNA copy number (mtCN) and heteroplasmy levels using blood-derived whole genome sequences from 274,832 individuals of diverse ancestry and perform GWAS to identify nuclear loci controlling these traits. After careful correction for blood cell composition, we observe that mtCN declines linearly with age and is associated with 92 independent nuclear genetic loci. We find that nearly every individual carries heteroplasmic variants that obey two key patterns: (1) heteroplasmic single nucleotide variants are somatic mutations that accumulate sharply after age 70, while (2) heteroplasmic indels are maternally transmitted as mtDNA mixtures with resulting levels influenced by 42 independent nuclear loci involved in mtDNA replication, maintenance, and novel pathways. These nuclear loci do not appear to act by mtDNA mutagenesis, but rather, likely act by conferring a replicative advantage to specific mtDNA molecules. As an illustrative example, the most common heteroplasmy we identify is a length variant carried by >50% of humans at position m.302 within a G-quadruplex known to serve as a replication switch. We find that this heteroplasmic variant exerts cis -acting genetic control over mtDNA abundance and is itself under trans -acting genetic control of nuclear loci encoding protein components of this regulatory switch. Our study showcases how nuclear haplotype can privilege the replication of specific mtDNA molecules to shape mtCN and heteroplasmy dynamics in the human population.

Community-engaged ancient DNA project reveals diverse origins of 18th-century African descendants in Charleston, South Carolina

In this study, we present the results of community-engaged ancient DNA research initiated after the remains of 36 African-descended individuals dating to the late 18th century were unearthed in the port city of Charleston, South Carolina. The Gullah Society of Charleston, along with other Charleston community members, initiated a collaborative genomic study of these ancestors of presumed enslaved status, in an effort to visibilize their histories. We generated 18 low-coverage genomes and 31 uniparental haplotypes to assess their genetic origins and interrelatedness. Our results indicate that they have predominantly West and West-Central African genomic ancestry, with one individual exhibiting some genomic affiliation with populations in the Americas. Most were assessed as genetic males, and no autosomal kin were identified among them. Overall, this study expands our understanding of the colonial histories of African descendant populations in the US South.

Sık Kullanılan Bazı Hücre Hatları için Kalite Kontrol: Mikoplazma Kontaminasyon Tespiti, Sitokrom B ve Sitokrom Oksidaz Alt Birim I Genlerinin DNA Dizi Analizlerinin Gerçekleştirilmesi

Amaç: Laboratuvarlarda sık kullanılan serviks epitelyal karsinom (HeLa), insan periferal kan promiyelösitik lösemi (HL-60), fare C3/bağ dokusu (L929), Madin Darby köpek böbrek (MDCK), fare nöroblastom (Neuro-2a) gibi bazı hücre hatlarının mikoplazma kontaminasyon kontrollerinin yapılması, kimlik doğrulamalarının gerçekleştirilmesi ve klonalitelerinin belirlenmesidir.Yöntem: Bu çalışmada üç farklı türe ait beş hücre hattı kullanılmıştır. Çalışılan tüm hatların Bisbenzimid (Hoechst 33258) ile deoksiribonükleik asit (DNA) floresan işaretlemesi yapılarak mikoplazma kontaminasyonu kontrolleri gerçekleştirilmiştir. Hücre hatlarından DNA izolasyonları yapılmış, elde edilen DNA örneklerinden sitokrom B (CYTB) geninin bölgesel amplifikasyonu için L14816 ve H15173 primerleri; sitokrom oksidaz alt birim I (COI) geni için ise LCO 1490 and HCO 2198 primerleri kullanılmıştır. İlgili amplifikasyonların DNA dizi analiz sonuçları, biyoinformatik araçlar kullanılarak referans dizilerle karşılaştırmalı olarak değerlendirilmiştir.Bulgular: Çalışmada ilgili hücrelerin, Bisbenzimid (Hoechst 33258) ile üretici firmanın protokollerine göre belirlenen konsantrasyon ve sürede yapılan boyama sonucunda mikoplazma kontaminasyonuna rastlanılmamıştır. Ayrıca CYTB gen bölgesi için veritabanında yer alan referans dizi ile yapılan karşılaştırma sonucu HL-60 için %97; "HeLa, L929, MDCK, Neuro-2a” hücre hatları için ise %98 oranında benzerlik bulunmuştur. COI gen bölgesi için ise bu benzerlik oranları “HeLa, HL-60, L929, MDCK ve Neuro-2a” hücre hatları için sırasıyla %95, %99, %96, %96 ve %98 olarak bulunmuştur.Sonuç: Bu bağlamda, çalışmadan elde edilen Bisbenzimid (Hoechst 33258) işaretleme ve DNA dizi analiz sonuçları, pek çok araştırmada kullanılan bu hücre hatlarının kalitesi konusunda kabul edilebilir bir belirteç ve güven sağlamıştır.   

Digging into the admixture strata of current-day Canary Islanders based on mitogenomes

The conquest of the Canary Islands by Europeans began at the beginning of the 15th century and culminated in 1496 with the surrender of the aborigines. The collapse of the aboriginal population during the conquest and the arrival of settlers caused a drastic change in the demographic composition of the archipelago. To shed light on this historical process, we analyzed 896 mitogenomes of current inhabitants from the seven main islands. Our findings confirm the continuity of aboriginal maternal contributions and the persistence of their genetic footprints in the current population, even at higher levels (>60% on average) than previously evidenced. Moreover, the age estimates for most autochthonous founder lineages support a first aboriginal arrival to the islands at the beginning of the first millennium. We also revealed for the first time that the main recognizable genetic influences from Europe are from Portuguese and Galicians.

Genomic evidence for ancient human migration routes along South America's Atlantic coast

An increasing body of archaeological and genomic evidence has hinted at a complex settlement process of the Americas by humans. This is especially true for South America, where unexpected ancestral signals have raised perplexing scenarios for the early migrations into different regions of the continent. Here, we present ancient human genomes from the archaeologically rich Northeast Brazil and compare them to ancient and present-day genomic data. We find a distinct relationship between ancient genomes from Northeast Brazil, Lagoa Santa, Uruguay and Panama, representing evidence for ancient migration routes along South America's Atlantic coast. To further add to the existing complexity, we also detect greater Denisovan than Neanderthal ancestry in ancient Uruguay and Panama individuals. Moreover, we find a strong Australasian signal in an ancient genome from Panama. This work sheds light on the deep demographic history of eastern South America and presents a starting point for future fine-scale investigations on the regional level.

Dissecting the multifaceted contribution of the mitochondrial genome to autism spectrum disorder

Autism spectrum disorder (ASD) is a clinically heterogeneous class of neurodevelopmental conditions with a strong, albeit complex, genetic basis. The genetic architecture of ASD includes different genetic models, from monogenic transmission at one end, to polygenic risk given by thousands of common variants with small effects at the other end. The mitochondrial DNA (mtDNA) was also proposed as a genetic modifier for ASD, mostly focusing on maternal mtDNA, since the paternal mitogenome is not transmitted to offspring. We extensively studied the potential contribution of mtDNA in ASD pathogenesis and risk through deep next generation sequencing and quantitative PCR in a cohort of 98 families. While the maternally-inherited mtDNA did not seem to predispose to ASD, neither for haplogroups nor for the presence of pathogenic mutations, an unexpected influence of paternal mtDNA, apparently centered on haplogroup U, came from the Italian families extrapolated from the test cohort (n = 74) when compared to the control population. However, this result was not replicated in an independent Italian cohort of 127 families and it is likely due to the elevated paternal age at time of conception. In addition, ASD probands showed a reduced mtDNA content when compared to their unaffected siblings. Multivariable regression analyses indicated that variants with 15%-5% heteroplasmy in probands are associated to a greater severity of ASD based on ADOS-2 criteria, whereas paternal super-haplogroups H and JT were associated with milder phenotypes. In conclusion, our results suggest that the mtDNA impacts on ASD, significantly modifying the phenotypic expression in the Italian population. The unexpected finding of protection induced by paternal mitogenome in term of severity may derive from a role of mtDNA in influencing the accumulation of nuclear de novo mutations or epigenetic alterations in fathers' germinal cells, affecting the neurodevelopment in the offspring. This result remains preliminary and needs further confirmation in independent cohorts of larger size. If confirmed, it potentially opens a different perspective on how paternal non-inherited mtDNA may predispose or modulate other complex diseases.

The kringle IV type 2 domain variant 4925G>A causes the elusive association signal of the LPA pentanucleotide repeat

Lipoprotein(a) [Lp(a)] concentrations are regulated by the LPA gene mainly via the large kringle IV-type 2 (KIV-2) copy number variation and multiple causal variants. Early studies suggested an effect of long pentanucleotide repeat (PNR) alleles (10 and 11 repeats, PNR10 and PNR11) in the LPA promoter on gene transcription and found an association with lower Lp(a). Subsequent in vitro studies showed no effects on mRNA transcription, but the association with strongly decreased Lp(a) remained consistent. We investigated the isolated and combined effect of PNR10, PNR11, and the frequent splice site variant KIV-2 4925G>A on Lp(a) concentrations in the Cooperative Health Research in the Region of Augsburg F4 study by multiple quantile regression in single-SNP and joint models. Data on Lp(a), apolipoprotein(a) Western blot isoforms, and variant genotypes were available for 2,858 individuals. We found a considerable linkage disequilibrium between KIV-2 4925G>A and the alleles PNR10 and PNR11. In single-variant analysis adjusted for age, sex, and the shorter apo(a) isoform, we determined that both PNR alleles were associated with a highly significant Lp(a) decrease (PNR10: β = -14.43 mg/dl, 95% CI: -15.84, -13.02, P = 3.33e-84; PNR11: β = -17.21 mg/dl, 95% CI: -20.19, -14.23, P = 4.01e-29). However, a joint model, adjusting the PNR alleles additionally for 4925G>A, abolished the effect on Lp(a) (PNR10: β = +0.44 mg/dl, 95% CI: -1.73, 2.60, P = 0.69; PNR11: β = -1.52 mg/dl, 95% CI: -6.05, 3.00, P = 0.51). Collectively, we conclude that the previously reported Lp(a) decrease observed in pentanucleotide alleles PNR10 or PNR11 carriers results from a linkage disequilibrium with the frequent splicing mutation KIV-2 4925G>A.

The role of mitochondrial genome abundance in Alzheimer's disease

Mitochondrial dysfunction is an early and prominent feature of Alzheimer's disease (AD), with impaired energy metabolism preceding the onset of clinical symptoms. Here we propose an update to the mitochondrial dysfunction hypothesis of AD based on recent results examining the role of mitochondrial genome abundance in AD. In a large post mortem study, we show that lower brain mitochondrial genome abundance is associated with a greater odds of AD neuropathological change and worse cognitive performance. We hypothesize that lower mitochondrial genome abundance impairs mitochondrial function by reducing mitochondrial bioenergetics, thereby impacting neuronal and glial cell function. However, it remains to be determined if mitochondrial dysfunction causes, mediates, or is a by-product of AD pathogenesis. Additional support for this hypothesis will be generated by linking peripheral blood mitochondrial genome abundance to AD and establishing clinical trials of compounds that upregulate total mitochondrial genome abundance or boost mitochondrial mass.

A method for multiplexed full-length single-molecule sequencing of the human mitochondrial genome

Methods to reconstruct the mitochondrial DNA (mtDNA) sequence using short-read sequencing come with an inherent bias due to amplification and mapping. They can fail to determine the phase of variants, to capture multiple deletions and to cover the mitochondrial genome evenly. Here we describe a method to target, multiplex and sequence at high coverage full-length human mitochondrial genomes as native single-molecules, utilizing the RNA-guided DNA endonuclease Cas9. Combining Cas9 induced breaks, that define the mtDNA beginning and end of the sequencing reads, as barcodes, we achieve high demultiplexing specificity and delineation of the full-length of the mtDNA, regardless of the structural variant pattern. The long-read sequencing data is analysed with a pipeline where our custom-developed software, baldur, efficiently detects single nucleotide heteroplasmy to below 1%, physically determines phase and can accurately disentangle complex deletions. Our workflow is a tool for studying mtDNA variation and will accelerate mitochondrial research.

Accurate analysis of mitochondrial DNA is important for mitochondrial disease clinical research and diagnostics. Here, authors present a method using Cas9 cleavage, nanopore sequencing and a custom pipeline to identify pathogenic variants, deletions and accurately quantify heteroplasmy to below 1%.

A bioinformatics pipeline for estimating mitochondrial DNA copy number and heteroplasmy levels from whole genome sequencing data

Mitochondrial diseases are a heterogeneous group of disorders that can be caused by mutations in the nuclear or mitochondrial genome. Mitochondrial DNA (mtDNA) variants may exist in a state of heteroplasmy, where a percentage of DNA molecules harbor a variant, or homoplasmy, where all DNA molecules have the same variant. The relative quantity of mtDNA in a cell, or copy number (mtDNA-CN), is associated with mitochondrial function, human disease, and mortality. To facilitate accurate identification of heteroplasmy and quantify mtDNA-CN, we built a bioinformatics pipeline that takes whole genome sequencing data and outputs mitochondrial variants, and mtDNA-CN. We incorporate variant annotations to facilitate determination of variant significance. Our pipeline yields uniform coverage by remapping to a circularized chrM and by recovering reads falsely mapped to nuclear-encoded mitochondrial sequences. Notably, we construct a consensus chrM sequence for each sample and recall heteroplasmy against the sample's unique mitochondrial genome. We observe an approximately 3-fold increased association with age for heteroplasmic variants in non-homopolymer regions and, are better able to capture genetic variation in the D-loop of chrM compared to existing software. Our bioinformatics pipeline more accurately captures features of mitochondrial genetics than existing pipelines that are important in understanding how mitochondrial dysfunction contributes to disease.

Benchmarking Low-Frequency Variant Calling With Long-Read Data on Mitochondrial DNA

Background: Sequencing quality has improved over the last decade for long-reads, allowing for more accurate detection of somatic low-frequency variants. In this study, we used mixtures of mitochondrial samples with different haplogroups (i.e., a specific set of mitochondrial variants) to investigate the applicability of nanopore sequencing for low-frequency single nucleotide variant detection.

Methods: We investigated the impact of base-calling, alignment/mapping, quality control steps, and variant calling by comparing the results to a previously derived short-read gold standard generated on the Illumina NextSeq. For nanopore sequencing, six mixtures of four different haplotypes were prepared, allowing us to reliably check for expected variants at the predefined 5%, 2%, and 1% mixture levels. We used two different versions of Guppy for base-calling, two aligners (i.e., Minimap2 and Ngmlr), and three variant callers (i.e., Mutserve2, Freebayes, and Nanopanel2) to compare low-frequency variants. We used F₁ score measurements to assess the performance of variant calling.

Results: We observed a mean read length of 11 kb and a mean overall read quality of 15. Ngmlr showed not only higher F₁ scores but also higher allele frequencies (AF) of false-positive calls across the mixtures (mean F₁ score = 0.83; false-positive allele frequencies < 0.17) compared to Minimap2 (mean F₁ score = 0.82; false-positive AF < 0.06). Mutserve2 had the highest F₁ scores (5% level: F₁ score >0.99, 2% level: F₁ score >0.54, and 1% level: F₁ score >0.70) across all callers and mixture levels.

Conclusion: We here present the benchmarking for low-frequency variant calling with nanopore sequencing by identifying current limitations.

Benchmarking the Effectiveness and Accuracy of Multiple Mitochondrial DNA Variant Callers: Practical Implications for Clinical Application

Mitochondrial DNA (mtDNA) mutations contribute to human disease across a range of severity, from rare, highly penetrant mutations causal for monogenic disorders to mutations with milder contributions to phenotypes. mtDNA variation can exist in all copies of mtDNA or in a percentage of mtDNA copies and can be detected with levels as low as 1%. The large number of copies of mtDNA and the possibility of multiple alternative alleles at the same DNA nucleotide position make the task of identifying allelic variation in mtDNA very challenging. In recent years, specialized variant calling algorithms have been developed that are tailored to identify mtDNA variation from whole-genome sequencing (WGS) data. However, very few studies have systematically evaluated and compared these methods for the detection of both homoplasmy and heteroplasmy. A publicly available synthetic gold standard dataset was used to assess four mtDNA variant callers (Mutserve, mitoCaller, MitoSeek, and MToolBox), and the commonly used Genome Analysis Toolkit “best practices” pipeline, which is included in most current WGS pipelines. We also used WGS data from 126 trios and calculated the percentage of maternally inherited variants as a metric of calling accuracy, especially for homoplasmic variants. We additionally compared multiple pathogenicity prediction resources for mtDNA variants. Although the accuracy of homoplasmic variant detection was high for the majority of the callers with high concordance across callers, we found a very low concordance rate between mtDNA variant callers for heteroplasmic variants ranging from 2.8% to 3.6%, for heteroplasmy thresholds of 5% and 1%. Overall, Mutserve showed the best performance using the synthetic benchmark dataset. The analysis of mtDNA pathogenicity resources also showed low concordance in prediction results. We have shown that while homoplasmic variant calling is consistent between callers, there remains a significant discrepancy in heteroplasmic variant calling. We found that resources like population frequency databases and pathogenicity predictors are now available for variant annotation but still need refinement and improvement. With its peculiarities, the mitochondria require special considerations, and we advocate that caution needs to be taken when analyzing mtDNA data from WGS data.

The Value of Whole-Genome Sequencing for Mitochondrial DNA Population Studies: Strategies and Criteria for Extracting High-Quality Mitogenome Haplotypes

Whole-genome sequencing (WGS) data present a readily available resource for mitochondrial genome (mitogenome) haplotypes that can be utilized for genetics research including population studies. However, the reconstruction of the mitogenome is complicated by nuclear mitochondrial DNA (mtDNA) segments (NUMTs) that co-align with the mtDNA sequences and mimic authentic heteroplasmy. Two minimum variant detection thresholds, 5% and 10%, were assessed for the ability to produce authentic mitogenome haplotypes from a previously generated WGS dataset. Variants associated with NUMTs were detected in the mtDNA alignments for 91 of 917 (~8%) Swedish samples when the 5% frequency threshold was applied. The 413 observed NUMT variants were predominantly detected in two regions (nps 12,612–13,105 and 16,390–16,527), which were consistent with previously documented NUMTs. The number of NUMT variants was reduced by ~97% (400) using a 10% frequency threshold. Furthermore, the 5% frequency data were inconsistent with a platinum-quality mitogenome dataset with respect to observed heteroplasmy. These analyses illustrate that a 10% variant detection threshold may be necessary to ensure the generation of reliable mitogenome haplotypes from WGS data resources.

Mitochondrial DNA variation across 56,434 individuals in gnomAD

Genomic databases of allele frequency are extremely helpful for evaluating clinical variants of unknown significance; however, until now, databases such as the Genome Aggregation Database (gnomAD) have focused on nuclear DNA and have ignored the mitochondrial genome (mtDNA). Here we present a pipeline to call mtDNA variants that addresses three technical challenges: (i) detecting homoplasmic and heteroplasmic variants, present respectively in all or a fraction of mtDNA molecules, (ii) circular mtDNA genome, and (iii) misalignment of nuclear sequences of mitochondrial origin (NUMTs). We observed that mtDNA copy number per cell varied across gnomAD cohorts and influenced the fraction of NUMT-derived false-positive variant calls, which can account for the majority of putative heteroplasmies. To avoid false positives, we excluded contaminated samples, cell lines, and samples prone to NUMT misalignment due to few mtDNA copies. Furthermore, we report variants with heteroplasmy greater than 10%. We applied this pipeline to 56,434 whole genome sequences in the gnomAD v3.1 database that includes individuals of European (58%), African (25%), Latino (10%), and Asian (5%) ancestry. Our gnomAD v3.1 release contains population frequencies for 10,850 unique mtDNA variants at more than half of all mtDNA bases. We report frequencies within each nuclear ancestral population and mitochondrial haplogroup. Homoplasmic variants account for most variant calls (98%) and unique variants (85%). We observed that 1/250 individuals carry a pathogenic mtDNA variant with heteroplasmy above 10%. These mtDNA population allele frequencies are freely accessible and will aid in diagnostic interpretation and research studies.

mitoDataclean: A machine learning approach for the accurate identification of cross-contamination-derived tumor mitochondrial DNA mutations

Next-generation sequencing (NGS) of mitochondrial DNA (mtDNA) has widespread applications in aging and cancer studies. However, cross-contamination of mtDNA constitutes a major concern. Previous methods for the detection of mtDNA contamination mainly focus on haplogroup-level phylogeny, but neglect haplotype-level differences, leading to limited sensitivity and accuracy. In this study, we present mitoDataclean, a random-forest-based machine learning package for accurate identification of cross-contamination, evaluation of contamination levels and detection of contamination-derived variants in mtDNA NGS data. Comprehensive optimization of mitoDataclean revealed that training simulation with mixtures of small haplogroup distance and low polymorphic difference was critical for optimal modeling. Compared with existing methods, mitoDataclean exhibited significantly improved sensitivity and accuracy for the detection of sample contamination in simulated data. In addition, mitoDataclean achieved area under the curve values of 0.91 and 0.97 for discerning genuine and contamination-derived mtDNA variants in a simulated Western dataset and private sequencing contamination data, respectively, suggesting that this tool may be applicable for different populations and samples with different sources of contamination. Finally, mitoDataclean was further evaluated in several private and public datasets and showed a robust ability for contamination detection. Altogether, our study demonstrates that mitoDataclean may be used for accurate detection of contaminated samples and contamination-derived variants in mtDNA NGS data. This article is protected by copyright. All rights reserved.

Identifying Somatic Mitochondrial DNA Mutations

Mitochondria are cellular organelles that play an essential role in eukaryotes, producing the energy needed for a cell to survive. Beyond the ~3.2 Gb of nuclear genomic DNA, each human cell has hundreds of mitochondria which carry one or a few copies of the 16.5 kb circular mitochondria DNA (mtDNA). Despite its small size, the circular genome encodes 37 genes, including 13 proteins that generate respiratory chain complexes together with other proteins of nuclear origin. Similar to nuclear genome, mtDNA in cancer cells frequently harbor somatically acquired alterations. Whole-genome or whole-exome sequencing of the tumor and its matched normal tissues (frequently blood or adjacent non-tumor tissues) enables sensitive and efficient detection of somatic mtDNA mutations. Because each cancer cell commonly carries hundreds to thousands of mtDNA copies, detection of mtDNA mutations is dependent on the heteroplasmic level of each mutation. Here, we describe strategies to accurately identify somatic mtDNA mutations in cancer genome studies.

Estimation of intrafamilial DNA contamination in family trio genome sequencing using deviation from Mendelian inheritance

With the increasing number of sequencing projects involving families, quality control tools optimized for family genome sequencing are needed. However, accurately quantifying contamination in a DNA mixture is particularly difficult when genetically related family members are the sources. We developed TrioMix, a maximum likelihood estimation (MLE) framework based on Mendel's law of inheritance, to quantify DNA mixture between family members in genome sequencing data of parent-offspring trios. TrioMix can accurately deconvolute any intrafamilial DNA contamination, including parent-offspring, sibling-sibling, parent-parent, and even multiple familial sources. In addition, TrioMix can be applied to detect genomic abnormalities that deviate from Mendelian inheritance patterns, such as uniparental disomy (UPD) and chimerism. A genome-wide depth and variant allele frequency plot generated by TrioMix facilitates tracing the origin of Mendelian inheritance deviations. We showed that TrioMix could accurately deconvolute genomes in both simulated and real data sets.

From Forensics to Clinical Research: Expanding the Variant Calling Pipeline for the Precision ID mtDNA Whole Genome Panel

Despite a multitude of methods for the sample preparation, sequencing, and data analysis of mitochondrial DNA (mtDNA), the demand for innovation remains, particularly in comparison with nuclear DNA (nDNA) research. The Applied Biosystems™ Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific, USA) is an innovative library preparation kit suitable for degraded samples and low DNA input. However, its bioinformatic processing occurs in the enterprise Ion Torrent Suite™ Software (TSS), yielding BAM files aligned to an unorthodox version of the revised Cambridge Reference Sequence (rCRS), with a heteroplasmy threshold level of 10%. Here, we present an alternative customizable pipeline, the PrecisionCallerPipeline (PCP), for processing samples with the correct rCRS output after Ion Torrent sequencing with the Precision ID library kit. Using 18 samples (3 original samples and 15 mixtures) derived from the 1000 Genomes Project, we achieved overall improved performance metrics in comparison with the proprietary TSS, with optimal performance at a 2.5% heteroplasmy threshold. We further validated our findings with 50 samples from an ongoing independent cohort of stroke patients, with PCP finding 98.31% of TSS's variants (TSS found 57.92% of PCP's variants), with a significant correlation between the variant levels of variants found with both pipelines.

Whole mitochondrial genomes assembled from thermally altered forensic bones and teeth

The recovery and analysis of genetic material obtained from thermally altered human bones and teeth are increasingly important to forensic investigations, especially in cases where soft-tissue identification is no longer possible. Although little is known about how these fire-related processes affect DNA degradation over time, next-generation sequencing technology in combination with traditional osteobiographical applications may provide us clues to these questions. In this study, we compare whole mitochondrial genome data generated using two different DNA extraction methods from 27 thermally altered samples obtained from fire victims (Maricopa County, Arizona) . DNA extracts were converted to double-stranded DNA libraries and enriched for whole mitochondrial DNA (mtDNA) using synthetic biotinylated RNA baits, then sequenced on an Illumina MiSeq. We processed the mitochondrial data using an in-house computational pipeline (MitoPipe1.0) composed of ancient DNA and modern genomics applications, then compared the resulting information across the two extraction types and five burn categories. Our analysis shows that DNA fragmentation increases with temperature, but that the acute insult from fire combined with the lack of water is insufficient to produce 5' and 3' terminal deamination characteristic of ancient DNA. Our data also suggest an acute and significant point of DNA degradation between 350 °C and 550 °C, and that the likelihood of generating high quality mtDNA haplogroup calls decreases significantly at temperatures > 550 °C. This research is part of a concerted effort to understand how fire affects our ability to generate genetic profiles suitable for forensic identification purposes.

An in-depth analysis of the mitochondrial phylogenetic landscape of Cambodia

Cambodia harbours a variety of human aboriginal populations that have scarcely been studied in terms of genetic diversity of entire mitochondrial genomes. Here we present the matrilineal gene pool of 299 Cambodian refugees from three different ethnic groups (Cham, Khmer, and Khmer Loeu) deriving from 16 Cambodian districts. After establishing a DNA-saving high-throughput strategy for mitochondrial whole-genome Sanger sequencing, a HaploGrep based workflow was used for quality control, haplogroup classification and phylogenetic reconstruction. The application of diverse phylogenetic algorithms revealed an exciting picture of the genetic diversity of Cambodia, especially in relation to populations from Southeast Asia and from the whole world. A total of 224 unique haplotypes were identified, which were mostly classified under haplogroups B5a1, F1a1, or categorized as newly defined basal haplogroups or basal sub-branches of R, N and M clades. The presence of autochthonous maternal lineages could be confirmed as reported in previous studies. The exceptional homogeneity observed between and within the three investigated Cambodian ethnic groups indicates genetic isolation of the whole population. Between ethnicities, genetic barriers were not detected. The mtDNA data presented here increases the phylogenetic resolution in Cambodia significantly, thereby highlighting the need for an update of the current human mtDNA phylogeny.