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.

Haplogrep 3 - an interactive haplogroup classification and analysis platform.

Over the last decade, Haplogrep has become a standard tool for haplogroup classification in the field of human mitochondrial DNA and is widely used by medical, forensic, and evolutionary researchers. Haplogrep scales well for thousands of samples, supports many file formats and provides an intuitive graphical web interface. Nevertheless, the currently available version has limitations when applying it to large biobank-scale data. In this paper, we present a major upgrade to the software by adding (a) haplogroup summary statistics and variant annotations from various publicly available genome databases, (b) an interface to connect new phylogenetic trees, (c) a new state-of-the-art web framework managing large scale data, (d) algorithmic adaptions to improve FASTA classification using BWA-specific alignment rules and (e) a pre-classification quality control step for VCF samples. These improvements will give researchers the opportunity to classify thousands of samples as usual but providing additional ways to investigate the dataset directly in the browser. The web service and its documentation can be accessed freely without any registration at https://haplogrep.i-med.ac.at.

First mitochondrial genome-wide association study with metabolomics.

In the era of personalized medicine with more and more patient-specific targeted therapies being used, we need reliable, dynamic, faster and sensitive biomarkers both to track the causes of disease and to develop and evolve therapies during the course of treatment. Metabolomics recently has shown substantial evidence to support its emerging role in disease diagnosis and prognosis. Aside from biomarkers and development of therapies, it is also an important goal to understand the involvement of mitochondrial DNA (mtDNA) in metabolic regulation, aging and disease development. Somatic mutations of the mitochondrial genome are also heavily implicated in age-related disease and aging. The general hypothesis is that an alteration in the concentration of metabolite profiles (possibly conveyed by lifestyle and environmental factors) influences the increase of mutation rate in the mtDNA and thereby contributes to a range of pathophysiological alterations observed in complex diseases. We performed an inverted mitochondrial genome-wide association analysis between mitochondrial nucleotide variants (mtSNVs) and concentration of metabolites. We used 151 metabolites and the whole sequenced mitochondrial genome from 2718 individuals to identify the genetic variants associated with metabolite profiles. Because of the high coverage, next-generation sequencing-based analysis of the mitochondrial genome allows for an accurate detection of mitochondrial heteroplasmy and for the identification of variants associated with the metabolome. The strongest association was found for mt715G > A located in the MT-12SrRNA with the metabolite ratio of C2/C10:1 (P-value = 6.82*10-09, ß = 0.909). The second most significant mtSNV was found for mt3714A > G located in the MT-ND1 with the metabolite ratio of phosphatidylcholine (PC) ae C42:5/PC ae C44:5 (P-value = 1.02*10-08, ß = 3.631). A large number of significant metabolite ratios were observed involving PC aa C36:6 and the variant mt10689G > A, located in the MT-ND4L gene. These results show an important interconnection between mitochondria and metabolite concentrations. Considering that some of the significant metabolites found in this study have been previously related to complex diseases, such as neurological disorders and metabolic conditions, these associations found here might play a crucial role for further investigations of such complex diseases. Understanding the mechanisms that control human health and disease, in particular, the role of genetic predispositions and their interaction with environmental factors is a prerequisite for the development of safe and efficient therapies for complex disorders.

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.

Apolipoprotein A-IV concentrations and cancer in a large cohort of chronic kidney disease patients: results from the GCKD study.

BACKGROUND: Chronic kidney disease (CKD) is highly connected to inflammation and oxidative stress. Both favour the development of cancer in CKD patients. Serum apolipoprotein A-IV (apoA-IV) concentrations are influenced by kidney function and are an early marker of kidney impairment. Besides others, it has antioxidant and anti-inflammatory properties. Proteomic studies and small case-control studies identified low apoA-IV as a biomarker for various forms of cancer; however, prospective studies are lacking. We therefore investigated whether serum apoA-IV is associated with cancer in the German Chronic Kidney Disease (GCKD) study. METHODS: These analyses include 5039 Caucasian patients from the prospective GCKD cohort study followed for 6.5 years. Main inclusion criteria were an eGFR of 30-60 mL/min/1.73m2 or an eGFR > 60 mL/min/1.73m2 in the presence of overt proteinuria. RESULTS: Mean apoA-IV concentrations of the entire cohort were 28.9 ± 9.8 mg/dL (median 27.6 mg/dL). 615 patients had a history of cancer before the enrolment into the study. ApoA-IV concentrations above the median were associated with a lower odds for a history of cancer (OR = 0.79, p = 0.02 when adjusted age, sex, smoking, diabetes, BMI, albuminuria, statin intake, and eGFRcreatinine). During follow-up 368 patients developed an incident cancer event and those with apoA-IV above the median had a lower risk (HR = 0.72, 95%CI 0.57-0.90, P = 0.004). Finally, 62 patients died from such an incident cancer event and each 10 mg/dL higher apoA-IV concentrations were associated with a lower risk for fatal cancer (HR = 0.62, 95%CI 0.44-0.88, P = 0.007). CONCLUSIONS: Our data indicate an association of high apoA-IV concentrations with reduced frequencies of a history of cancer as well as incident fatal and non-fatal cancer events in a large cohort of patients with CKD.

Interactive exploration of adverse events and multimorbidity in CKD.

BACKGROUND AND HYPOTHESIS: Persons with chronic kidney disease (CKD) are at increased risk of adverse events, early mortality, and multimorbidity. A detailed overview of adverse event types and rates from a large CKD cohort under regular nephrological care is missing. We generated an interactive tool to enable exploration of adverse events and their combinations in the prospective, observational German CKD (GCKD) study. METHODS: The GCKD study enrolled 5217 participants under regular nephrological care with an estimated glomerular filtration rate of 30-60 or >60 mL/min/1.73m2 and an overt proteinuria. Cardio-, cerebro- and peripheral vascular, kidney, infection, and cancer events, as well as deaths were adjudicated following a standard operation procedure. We summarized these time-to-event data points for exploration in interactive graphs within an R shiny app. Multivariable adjusted Cox models for time to first event were fitted. Cumulative incidence functions, Kaplan-Meier curves and intersection plots were used to display main adverse events and their combinations by sex and CKD etiology. RESULTS: Over a median of 6.5 years, 10 271 events occurred in total and 680 participants (13.0%) died while 2947 participants (56.5%) experienced any event. The new publicly available interactive platform enables readers to scrutinize adverse events and their combinations as well as mortality trends as a gateway to better understand multimorbidity in CKD: incident rates per 1000 patient-years varied by event type, CKD etiology, and baseline characteristics. Incidence rates for the most frequent events and their recurrence were 113.6 (cardiovascular), 75.0 (kidney), and 66.0 (infection). Participants with diabetic kidney disease and men were more prone to experiencing events. CONCLUSION: This comprehensive explorative tool to visualize adverse events (https://gckd.diz.uk-erlangen.de/), their combination, mortality, and multimorbidity among persons with CKD may manifest as a valuable resource for patient care, identification of high-risk groups, health services, and public health policy planning.

Mitochondrial DNA heteroplasmy distinguishes disease manifestation in PINK1/PRKN-linked Parkinson's disease.

Biallelic mutations in PINK1/PRKN cause recessive Parkinson's disease. Given the established role of PINK1/Parkin in regulating mitochondrial dynamics, we explored mitochondrial DNA integrity and inflammation as disease modifiers in carriers of mutations in these genes. Mitochondrial DNA integrity was investigated in a large collection of biallelic (n = 84) and monoallelic (n = 170) carriers of PINK1/PRKN mutations, idiopathic Parkinson's disease patients (n = 67) and controls (n = 90). In addition, we studied global gene expression and serum cytokine levels in a subset. Affected and unaffected PINK1/PRKN monoallelic mutation carriers can be distinguished by heteroplasmic mitochondrial DNA variant load (area under the curve = 0.83, CI 0.74-0.93). Biallelic PINK1/PRKN mutation carriers harbour more heteroplasmic mitochondrial DNA variants in blood (P = 0.0006, Z = 3.63) compared to monoallelic mutation carriers. This enrichment was confirmed in induced pluripotent stem cell-derived (controls, n = 3; biallelic PRKN mutation carriers, n = 4) and post-mortem (control, n = 1; biallelic PRKN mutation carrier, n = 1) midbrain neurons. Last, the heteroplasmic mitochondrial DNA variant load correlated with IL6 levels in PINK1/PRKN mutation carriers (r = 0.57, P = 0.0074). PINK1/PRKN mutations predispose individuals to mitochondrial DNA variant accumulation in a dose- and disease-dependent manner.

Systemic Evidence for Mitochondrial Dysfunction in Age-Related Macular Degeneration as Revealed by mtDNA Copy Number Measurements in Peripheral Blood.

Mitochondrial dysfunction is a common occurrence in the aging process and is observed in diseases such as age-related macular degeneration (AMD). Increased levels of reactive oxygen species lead to damaged mitochondrial DNA (mtDNA), resulting in dysfunctional mitochondria, and, consequently, mtDNA causes further harm in the retinal tissue. However, it is unclear whether the effects are locally restricted to the high-energy-demanding retinal pigment epithelium or are also systematically present. Therefore, we measured mtDNA copy number (mtDNA-CN) in peripheral blood using a qPCR approach with plasmid normalization in elderly participants with and without AMD from the AugUR study (n = 2262). We found significantly lower mtDNA-CN in the blood of participants with early (n = 453) and late (n = 170) AMD compared to AMD-free participants (n = 1630). In regression analyses, we found lower mtDNA-CN to be associated with late AMD when compared with AMD-free participants. Each reduction of mtDNA-CN by one standard deviation increased the risk for late AMD by 24%. This association was most pronounced in geographic atrophy (OR = 1.76, 95% CI 1.19-2.60, p = 0.004), which has limited treatment options. These findings provide new insights into the relationship between mtDNA-CN in blood and AMD, suggesting that it may serve as a more accessible biomarker than mtDNA-CN in the retina.

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.

Apolipoprotein A-IV concentrations and clinical outcomes in a large chronic kidney disease cohort: Results from the GCKD study.

BACKGROUND: Chronic kidney disease (CKD) represents a chronic proinflammatory state and is associated with very high cardiovascular risk. Apolipoprotein A-IV (apoA-IV) has antiatherogenic, antioxidative, anti-inflammatory and antithrombotic properties and levels increase significantly during the course of CKD. OBJECTIVES: We aimed to investigate the association between apoA-IV and all-cause mortality and cardiovascular outcomes in the German Chronic Kidney Disease study. METHODS: This was a prospective cohort study including 5141 Caucasian patients with available apoA-IV measurements and CKD. The majority of the patients had an estimated glomerular filtration rate (eGFR) of 30-60 ml/min/1.73m2 or an eGFR >60 ml/min/1.73m2 in the presence of overt proteinuria. Median follow-up was 6.5 years. The association of apoA-IV with comorbidities at baseline and endpoints during follow-up was modelled adjusting for major confounders. RESULTS: Mean apoA-IV concentrations of the entire cohort were 28.9 ± 9.8 mg/dl. Patients in the highest apoA-IV quartile had the lowest high-sensitivity C-reactive protein values despite the highest prevalence of diabetes, albuminuria and the lowest eGFR. Each 10 mg/dl higher apoA-IV translated into lower odds of prevalent cardiovascular disease (1289 cases, odds ratio = 0.80, 95% confidence interval [CI] 0.72-0.86, p = 0.0000003). During follow-up, each 10 mg/dl higher apoA-IV was significantly associated with a lower risk for all-cause mortality (600 cases, hazard ratio [HR] = 0.81, 95% CI 0.73-0.89, p = 0.00004), incident major adverse cardiovascular events (506 cases, HR = 0.88, 95% CI 0.79-0.99, p = 0.03) and death or hospitalizations due to heart failure (346 cases, HR = 0.84, 95% CI 0.73-0.96, p = 0.01). CONCLUSIONS: These data support a link between elevated apoA-IV concentrations and reduced inflammation in moderate CKD. ApoA-IV appears to be an independent risk marker for reduced all-cause mortality, cardiovascular events and heart failure in a large cohort of patients with CKD.

CovidPhy: A tool for phylogeographic analysis of SARS-CoV-2 variation.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the pathogen responsible for the coronavirus disease 2019 (COVID-19) pandemic. SARS-CoV-2 genomes have been sequenced massively and worldwide and are now available in different public genome repositories. There is much interest in generating bioinformatic tools capable to analyze and interpret SARS-CoV-2 variation. We have designed CovidPhy (http://covidphy.eu), a web interface that can process SARS-CoV-2 genome sequences in plain fasta text format or provided through identity codes from the Global Initiative on Sharing Avian Influenza Data (GISAID) or GenBank. CovidPhy aggregates information available on the large GISAID database (>1.49 M genomes). Sequences are first aligned against the reference sequence and the interface provides different sources of information, including automatic classification of genomes into a pre-computed phylogeny and phylogeographic information, haplogroup/lineage frequencies, and sequencing variation, indicating also if the genome contains known variants of concern (VOC). Additionally, CovidPhy allows searching for variants and haplotypes introduced by the user and includes a list of genomes that are good candidates for being responsible for large outbreaks worldwide, most likely mediated by important superspreading events, indicating their possible geographic epicenters and their relative impact as recorded in the GISAID database.

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.

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.

Results from the German Chronic Kidney Disease (GCKD) study support association of relative telomere length with mortality in a large cohort of patients with moderate chronic kidney disease.

Telomere length is known to be inversely associated with aging and has been proposed as a marker for aging-related diseases. Telomere attrition can be accelerated by oxidative stress and inflammation, both commonly present in patients with chronic kidney disease. Here, we investigated whether relative telomere length is associated with mortality in a large cohort of patients with chronic kidney disease stage G3 and A1-3 or G1-2 with overt proteinuria (A3) at enrollment. Relative telomere length was quantified in peripheral blood by a quantitative PCR method in 4,955 patients from the GCKD study, an ongoing prospective observational cohort. Complete four-year follow-up was available from 4,926 patients in whom we recorded 354 deaths. Relative telomere length was a strong and independent predictor of all-cause mortality. Each decrease of 0.1 relative telomere length unit was highly associated with a 14% increased risk of death (hazard ratio1.14 [95% confidence interval 1.06-1.22]) in a model adjusted for age, sex, baseline eGFR, urine albumin/creatinine ratio, diabetes mellitus, prevalent cardiovascular disease, LDL-cholesterol, HDL-cholesterol, smoking, body mass index, systolic and diastolic blood pressure, C-reactive protein and serum albumin. This translated to a 75% higher risk for those in the lowest compared to the highest quartile of relative telomere length. The association was mainly driven by 117 cardiovascular deaths (1.20 [1.05-1.35]) as well as 67 deaths due to infections (1.27 [1.07-1.50]). Thus, our findings support an association of shorter telomere length with all-cause mortality, cardiovascular mortality and death due to infections in patients with moderate chronic kidney disease.

A comprehensive map of single-base polymorphisms in the hypervariable LPA kringle IV type 2 copy number variation region.

Lipoprotein (a) [Lp(a)] concentrations are among the strongest genetic risk factors for cardiovascular disease and present pronounced interethnic and interindividual differences. Approximately 90% of Lp(a) variance is controlled by the LPA gene, which contains a 5.6-kb-large copy number variation [kringle IV type 2 (KIV-2) repeat] that generates >40 protein isoforms. Variants within the KIV-2 region are not called in common sequencing projects, leaving up to 70% of the LPA coding region currently unaddressed. To completely assess the variability in LPA, we developed a sequencing strategy for this region and report here the first map of genetic variation in the KIV-2 region, a comprehensively evaluated ultradeep sequencing protocol, and an easy-to-use variant analysis pipeline. We sequenced 123 Central-European individuals and reanalyzed public data of 2,504 individuals from 26 populations. We found 14 different loss-of-function and splice-site mutations, as well as >100, partially even common, missense variants. Some coding variants were frequent in one population but absent in others. This provides novel candidates to explain the large ethnic and individual differences in Lp(a) concentrations. Importantly, our approach and pipeline are also applicable to other similar copy number variable regions, allowing access to regions that are not captured by common genome sequencing.

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.

mtDNA-Server: next-generation sequencing data analysis of human mitochondrial DNA in the cloud.

Next generation sequencing (NGS) allows investigating mitochondrial DNA (mtDNA) characteristics such as heteroplasmy (i.e. intra-individual sequence variation) to a higher level of detail. While several pipelines for analyzing heteroplasmies exist, issues in usability, accuracy of results and interpreting final data limit their usage. Here we present mtDNA-Server, a scalable web server for the analysis of mtDNA studies of any size with a special focus on usability as well as reliable identification and quantification of heteroplasmic variants. The mtDNA-Server workflow includes parallel read alignment, heteroplasmy detection, artefact or contamination identification, variant annotation as well as several quality control metrics, often neglected in current mtDNA NGS studies. All computational steps are parallelized with Hadoop MapReduce and executed graphically with Cloudgene. We validated the underlying heteroplasmy and contamination detection model by generating four artificial sample mix-ups on two different NGS devices. Our evaluation data shows that mtDNA-Server detects heteroplasmies and artificial recombinations down to the 1% level with perfect specificity and outperforms existing approaches regarding sensitivity. mtDNA-Server is currently able to analyze the 1000G Phase 3 data (n = 2,504) in less than 5 h and is freely accessible at https://mtdna-server.uibk.ac.at.

HaploGrep 2: mitochondrial haplogroup classification in the era of high-throughput sequencing.

Mitochondrial DNA (mtDNA) profiles can be classified into phylogenetic clusters (haplogroups), which is of great relevance for evolutionary, forensic and medical genetics. With the extensive growth of the underlying phylogenetic tree summarizing the published mtDNA sequences, the manual process of haplogroup classification would be too time-consuming. The previously published classification tool HaploGrep provided an automatic way to address this issue. Here, we present the completely updated version HaploGrep 2 offering several advanced features, including a generic rule-based system for immediate quality control (QC). This allows detecting artificial recombinants and missing variants as well as annotating rare and phantom mutations. Furthermore, the handling of high-throughput data in form of VCF files is now directly supported. For data output, several graphical reports are generated in real time, such as a multiple sequence alignment format, a VCF format and extended haplogroup QC reports, all viewable directly within the application. In addition, HaploGrep 2 generates a publication-ready phylogenetic tree of all input samples encoded relative to the revised Cambridge Reference Sequence. Finally, new distance measures and optimizations of the algorithm increase accuracy and speed-up the application. HaploGrep 2 can be accessed freely and without any registration at http://haplogrep.uibk.ac.at.

Origin and dynamics of admixture in Brazilians and its effect on the pattern of deleterious mutations.

While South Americans are underrepresented in human genomic diversity studies, Brazil has been a classical model for population genetics studies on admixture. We present the results of the EPIGEN Brazil Initiative, the most comprehensive up-to-date genomic analysis of any Latin-American population. A population-based genome-wide analysis of 6,487 individuals was performed in the context of worldwide genomic diversity to elucidate how ancestry, kinship, and inbreeding interact in three populations with different histories from the Northeast (African ancestry: 50%), Southeast, and South (both with European ancestry >70%) of Brazil. We showed that ancestry-positive assortative mating permeated Brazilian history. We traced European ancestry in the Southeast/South to a wider European/Middle Eastern region with respect to the Northeast, where ancestry seems restricted to Iberia. By developing an approximate Bayesian computation framework, we infer more recent European immigration to the Southeast/South than to the Northeast. Also, the observed low Native-American ancestry (6-8%) was mostly introduced in different regions of Brazil soon after the European Conquest. We broadened our understanding of the African diaspora, the major destination of which was Brazil, by revealing that Brazilians display two within-Africa ancestry components: one associated with non-Bantu/western Africans (more evident in the Northeast and African Americans) and one associated with Bantu/eastern Africans (more present in the Southeast/South). Furthermore, the whole-genome analysis of 30 individuals (42-fold deep coverage) shows that continental admixture rather than local post-Columbian history is the main and complex determinant of the individual amount of deleterious genotypes.