Contamination detection in sequencing studies using the mitochondrial phylogeny.

Within-species contamination is a major issue in sequencing studies, especially for mitochondrial studies. Contamination can be detected by analyzing the nuclear genome or by inspecting polymorphic sites in the mitochondrial genome (mtDNA). Existing methods using the nuclear genome are computationally expensive, and no appropriate tool for detecting sample contamination in large-scale mtDNA data sets is available. Here we present haplocheck, a tool that requires only the mtDNA to detect contamination in both targeted mitochondrial and whole-genome sequencing studies. Our in silico simulations and amplicon mixture experiments indicate that haplocheck detects mtDNA contamination accurately and is independent of the phylogenetic distance within a sample mixture. By applying haplocheck to The 1000 Genomes Project Consortium data, we further evaluate the application of haplocheck as a fast proxy tool for nDNA-based contamination detection using the mtDNA and identify the mitochondrial copy number within a mixture as a critical component for the overall accuracy. The haplocheck tool is available both as a command-line tool and as a cloud web service producing interactive reports that facilitates the navigation through the phylogeny of contaminated samples.

Analyzing Low-Level mtDNA Heteroplasmy-Pitfalls and Challenges from Bench to Benchmarking.

Massive parallel sequencing technologies are promising a highly sensitive detection of low-level mutations, especially in mitochondrial DNA (mtDNA) studies. However, processes from DNA extraction and library construction to bioinformatic analysis include several varying tasks. Further, there is no validated recommendation for the comprehensive procedure. In this study, we examined potential pitfalls on the sequencing results based on two-person mtDNA mixtures. Therefore, we compared three DNA polymerases, six different variant callers in five mixtures between 50% and 0.5% variant allele frequencies generated with two different amplification protocols. In total, 48 samples were sequenced on Illumina MiSeq. Low-level variant calling at the 1% variant level and below was performed by comparing trimming and PCR duplicate removal as well as six different variant callers. The results indicate that sensitivity, specificity, and precision highly depend on the investigated polymerase but also vary based on the analysis tools. Our data highlight the advantage of prior standardization and validation of the individual laboratory setup with a DNA mixture model. Finally, we provide an artificial heteroplasmy benchmark dataset that can help improve somatic variant callers or pipelines, which may be of great interest for research related to cancer and aging.

Mitochondrial DNA copy number is associated with mortality and infections in a large cohort of patients with chronic kidney disease.

Damage of mitochondrial DNA (mtDNA) with reduction in copy number has been proposed as a biomarker for mitochondrial dysfunction and oxidative stress. Chronic kidney disease (CKD) is associated with increased mortality and risk of cardiovascular disease, but the underlying mechanisms remain incompletely understood. Here we investigated the prognostic role of mtDNA copy number for cause-specific mortality in 4812 patients from the German Chronic Kidney Disease study, an ongoing prospective observational national cohort study of patients with CKD stage G3 and A1-3 or G1-2 with overt proteinuria (A3) at enrollment. MtDNA was quantified in whole blood using a plasmid-normalized PCR-based assay. At baseline, 1235 patients had prevalent cardiovascular disease. These patients had a significantly lower mtDNA copy number than patients without cardiovascular disease (fully-adjusted model: odds ratio 1.03, 95% confidence interval [CI] 1.01-1.05 per 10 mtDNA copies decrease). After four years of follow-up, we observed a significant inverse association between mtDNA copy number and all-cause mortality, adjusted for kidney function and cardiovascular disease risk factors (hazard ratio 1.37, 95% CI 1.09-1.73 for quartile 1 compared to quartiles 2-4). When grouped by causes of death, estimates pointed in the same direction for all causes but in a fully-adjusted model decreased copy numbers were significantly lower only in infection-related death (hazard ratio 1.82, 95% CI 1.08-3.08). A similar association was observed for hospitalizations due to infections in 644 patients (hazard ratio 1.19, 95% CI 1.00-1.42 in the fully-adjusted model). Thus, our data support a role of mitochondrial dysfunction in increased cardiovascular disease and mortality risks as well as susceptibility to infections in patients with CKD.

Rapid coastal spread of First Americans: novel insights from South America's Southern Cone mitochondrial genomes.

It is now widely agreed that the Native American founders originated from a Beringian source population ~15-18 thousand years ago (kya) and rapidly populated all of the New World, probably mainly following the Pacific coastal route. However, details about the migration into the Americas and the routes pursued on the continent still remain unresolved, despite numerous genetic, archaeological, and linguistic investigations. To examine the pioneering peopling phase of the South American continent, we screened literature and mtDNA databases and identified two novel mitochondrial DNA (mtDNA) clades, here named D1g and D1j, within the pan-American haplogroup D1. They both show overall rare occurrences but local high frequencies, and are essentially restricted to populations from the Southern Cone of South America (Chile and Argentina). We selected and completely sequenced 43 D1g and D1j mtDNA genomes applying highest quality standards. Molecular and phylogeographic analyses revealed extensive variation within each of the two clades and possibly distinct dispersal patterns. Their age estimates agree with the dating of the earliest archaeological sites in South America and indicate that the Paleo-Indian spread along the entire longitude of the American double continent might have taken even <2000 yr. This study confirms that major sampling and sequencing efforts are mandatory for uncovering all of the most basal variation in the Native American mtDNA haplogroups and for clarification of Paleo-Indian migrations, by targeting, if possible, both the general mixed population of national states and autochthonous Native American groups, especially in South America.

Forensic and phylogeographic characterisation of mtDNA lineages from Somalia.

The African mitochondrial (mt) phylogeny is coarsely resolved but the majority of population data generated so far is limited to the analysis of the first hypervariable segment (HVS-1) of the control region (CR). Therefore, this study aimed on the investigation of the entire CR of 190 unrelated Somali individuals to enrich the severely underrepresented African mtDNA pool. The majority (60.5 %) of the haplotypes were of sub-Saharan origin with L0a1d, L2a1h and L3f being the most frequently observed haplogroups. This is in sharp contrast to previous data reported from the Y-chromosome, where only about 5 % of the observed haplogroups were of sub-Saharan provenance. We compared the genetic distances based on population pairwise F (st) values between 11 published East, Central and North African as well as western Asian populations and the Somali sequences and displayed them in a multi-dimensional scaling plot. Genetic proximity evidenced by clustering roughly reflected the relative geographic location of the populations. The sequences will be included in the EMPOP database ( www.empop.org ) under accession number EMP00397 upon publication (Parson and Dür Forensic Sci Int Genet 1:88-92, 2007).

Accumulation of mutations over the entire mitochondrial genome of breast cancer cells obtained by tissue microdissection.

The occurrence of heteroplasmy and mixtures is technically challenging for the analysis of mitochondrial DNA. More than that, observed mutations need to be carefully interpreted in the light of the phylogeny as mitochondrial DNA is a uniparental marker reflecting human evolution. Earlier attempts to explain the role of mtDNA in cancerous tissues led to substantial confusion in medical genetics mainly due to the presentation of low sequence data quality and misinterpretation of mutations representing a particular haplogroup background rather than being cancer-specific. The focus of this study is to characterize the extent and level of mutations in breast cancer samples obtained by tissue microdissection by application of an evaluated full mtDNA genome sequencing protocol. We amplified and sequenced the complete mitochondrial genomes of microdissected breast cancer cells of 15 patients and compared the results to those obtained from paired non-cancerous breast tissue derived from the same patients. We observed differences in the heteroplasmic states of substitutions between cancerous and normal cells, one of which was affecting a position that has been previously reported in lung cancer and another one that has been identified in 16 epithelial ovarian tumors, possibly indicating functional relevance. In the coding region, we found full transitions in two cancerous mitochondrial genomes and 12 heteroplasmic substitutions as compared to the non-cancerous breast cells. We identified somatic mutations over the entire mtDNA of human breast cancer cells potentially impairing the mitochondrial OXPHOS system.

Implications of Standardized Uptake Values of Oral Squamous Cell Carcinoma in PET-CT on Prognosis, Tumor Characteristics and Mitochondrial DNA Heteroplasmy.

Under aerobic conditions, some cancers switch to glycolysis to cover their energy requirements. Taking advantage of this process, functional imaging techniques such as PET-CT can be used to detect and assess tumorous tissues. The aim of this study was to investigate standardized uptake values and mitochondrial DNA mutations in oral squamous cell carcinoma. A cohort of 57 patients underwent 18[F]FDG-PET-CT and standardized uptake values were collected. In 15 patients, data on mitochondrial DNA mutations of the tumor were available. Kaplan-Meier curves were calculated, and correlation analyses as well as univariate Cox proportional hazard models were performed. Using ROC analysis to determine a statistical threshold for SUVmax in PET investigations, a cut-off value was determined at 9.765 MB/mL. Survival analysis for SUVmax in these groups showed a Hazard Ratio of 4 (95% CI 1.7-9) in the high SUVmax group with 5-year survival rates of 23.5% (p = 0.00042). For SUVmax and clinicopathological tumor features, significant correlations were found. A tendency towards higher mtDNA heteroplasmy levels in high SUVmax groups could be observed. We were able to confirm the prognostic value of SUVmax in OSCC, showing higher survival rates at lower SUVmax levels. Correlations between SUVmax and distinct tumor characteristics were highly significant, providing evidence that SUVmax may act as a reliable diagnostic parameter. Correlation analysis of mtDNA mutations suggests an influence on metabolic activity in OSCC.

OXPHOS remodeling in high-grade prostate cancer involves mtDNA mutations and increased succinate oxidation.

Rewiring of energy metabolism and adaptation of mitochondria are considered to impact on prostate cancer development and progression. Here, we report on mitochondrial respiration, DNA mutations and gene expression in paired benign/malignant human prostate tissue samples. Results reveal reduced respiratory capacities with NADH-pathway substrates glutamate and malate in malignant tissue and a significant metabolic shift towards higher succinate oxidation, particularly in high-grade tumors. The load of potentially deleterious mitochondrial-DNA mutations is higher in tumors and associated with unfavorable risk factors. High levels of potentially deleterious mutations in mitochondrial Complex I-encoding genes are associated with a 70% reduction in NADH-pathway capacity and compensation by increased succinate-pathway capacity. Structural analyses of these mutations reveal amino acid alterations leading to potentially deleterious effects on Complex I, supporting a causal relationship. A metagene signature extracted from the transcriptome of tumor samples exhibiting a severe mitochondrial phenotype enables identification of tumors with shorter survival times.

Profiling of Mitochondrial DNA Heteroplasmy in a Prospective Oral Squamous Cell Carcinoma Study.

While a shift in energy metabolism is essential to cancers, the knowledge about the involvement of the mitochondrial genome in tumorigenesis and progression in oral squamous cell carcinoma (OSCC) is still very limited. In this study, we evaluated 37 OSCC tumors and the corresponding benign mucosa tissue pairs by deep sequencing of the complete mitochondrial DNA (mtDNA). After extensive quality control, we identified 287 variants, 137 in tumor and 150 in benign samples exceeding the 1% threshold. Variant heteroplasmy levels were significantly increased in cancer compared to benign tissues (p = 0.0002). Furthermore, pairwise high heteroplasmy frequency difference variants (∆HF% > 20) with potential functional impact were increased in the cancer tissues (p = 0.024). Fourteen mutations were identified in the protein-coding region, out of which thirteen were detected in cancer and only one in benign tissue. After eight years of follow-up, the risk of mortality was higher for patients who harbored at least one ∆HF% > 20 variant in mtDNA protein-coding regions relative to those with no mutations (HR = 4.6, (95%CI = 1.3-17); p = 0.019 in primary tumor carriers). Haplogroup affiliation showed an impact on survival time, which however needs confirmation in a larger study. In conclusion, we observed a significantly higher accumulation of somatic mutations in the cancer tissues associated with a worse prognosis.

Sequencing strategy for the whole mitochondrial genome resulting in high quality sequences.

BACKGROUND: It has been demonstrated that a reliable and fail-safe sequencing strategy is mandatory for high-quality analysis of mitochondrial (mt) DNA, as the sequencing and base-calling process is prone to error. Here, we present a high quality, reliable and easy handling manual procedure for the sequencing of full mt genomes that is also appropriate for laboratories where fully automated processes are not available. RESULTS: We amplified whole mitochondrial genomes as two overlapping PCR-fragments comprising each about 8500 bases in length. We developed a set of 96 primers that can be applied to a (manual) 96 well-based technology, which resulted in at least double strand sequence coverage of the entire coding region (codR). CONCLUSION: This elaborated sequencing strategy is straightforward and allows for an unambiguous sequence analysis and interpretation including sometimes challenging phenomena such as point and length heteroplasmy that are relevant for the investigation of forensic and clinical samples.