Resolving intra-repeat variation in medically relevant VNTRs from short-read sequencing data using the cardiovascular risk gene LPA as a model.

BACKGROUND: Variable number tandem repeats (VNTRs) are highly polymorphic DNA regions harboring many potentially disease-causing variants. However, VNTRs often appear unresolved ("dark") in variation databases due to their repetitive nature. One particularly complex and medically relevant VNTR is the KIV-2 VNTR located in the cardiovascular disease gene LPA which encompasses up to 70% of the coding sequence. RESULTS: Using the highly complex LPA gene as a model, we develop a computational approach to resolve intra-repeat variation in VNTRs from largely available short-read sequencing data. We apply the approach to six protein-coding VNTRs in 2504 samples from the 1000 Genomes Project and developed an optimized method for the LPA KIV-2 VNTR that discriminates the confounding KIV-2 subtypes upfront. This results in an F1-score improvement of up to 2.1-fold compared to previously published strategies. Finally, we analyze the LPA VNTR in > 199,000 UK Biobank samples, detecting > 700 KIV-2 mutations. This approach successfully reveals new strong Lp(a)-lowering effects for KIV-2 variants, with protective effect against coronary artery disease, and also validated previous findings based on tagging SNPs. CONCLUSIONS: Our approach paves the way for reliable variant detection in VNTRs at scale and we show that it is transferable to other dark regions, which will help unlock medical information hidden in VNTRs.

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.

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.

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.

Nuclear and mitochondrial genetic variants associated with mitochondrial DNA copy number.

Mitochondrial DNA copy number (mtDNA-CN) is a biomarker for mitochondrial dysfunction associated with several diseases. Previous genome-wide association studies (GWAS) have been performed to unravel underlying mechanisms of mtDNA-CN regulation. However, the identified gene regions explain only a small fraction of mtDNA-CN variability. Most of this data has been estimated from microarrays based on various pipelines. In the present study we aimed to (1) identify genetic loci for qPCR-measured mtDNA-CN from three studies (16,130 participants) using GWAS, (2) identify potential systematic differences between our qPCR derived mtDNA-CN measurements compared to the published microarray intensity-based estimates, and (3) disentangle the nuclear from mitochondrial regulation of the mtDNA-CN phenotype. We identified two genome-wide significant autosomal loci associated with qPCR-measured mtDNA-CN: at HBS1L (rs4895440, p = 3.39 × 10-13) and GSDMA (rs56030650, p = 4.85 × 10-08) genes. Moreover, 113/115 of the previously published SNPs identified by microarray-based analyses were significantly equivalent with our findings. In our study, the mitochondrial genome itself contributed only marginally to mtDNA-CN regulation as we only detected a single rare mitochondrial variant associated with mtDNA-CN. Furthermore, we incorporated mitochondrial haplogroups into our analyses to explore their potential impact on mtDNA-CN. However, our findings indicate that they do not exert any significant influence on our results.

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.

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.

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.

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.

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 F1 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 F1 scores but also higher allele frequencies (AF) of false-positive calls across the mixtures (mean F1 score = 0.83; false-positive allele frequencies < 0.17) compared to Minimap2 (mean F1 score = 0.82; false-positive AF < 0.06). Mutserve2 had the highest F1 scores (5% level: F1 score >0.99, 2% level: F1 score >0.54, and 1% level: F1 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.