Search

Objective Adipose tissue contributes to adverse outcomes in chronic kidney disease (CKD), but there is uncertainty regarding the prognostic relevance of different adiposity measures. We analyzed the associations of neck circumference (NC), waist circumference (WC), and body mass index (BMI) with clinical outcomes in patients with mild to severe CKD. Methods The German Chronic Kidney Disease (GCKD) study is a prospective cohort study, which enrolled Caucasian adults with mild to severe CKD, defined as estimated glomerular filtration rate (eGFR): 30–60 mL/min/1.73 m2, or >60 mL/min/1.73 m2 in the presence of overt proteinuria. Associations of NC, WC and BMI with all-cause death, major cardiovascular events (MACE: a composite of non-fatal stroke, non-fatal myocardial infarction, peripheral artery disease intervention, and cardiovascular death), kidney failure (a composite of dialysis or transplantation) were analyzed using multivariable Cox proportional hazards regression models adjusted for confounders and the Akaike information criteria (AIC) were calculated. Models included sex interactions with adiposity measures. Results A total of 4537 participants (59% male) were included in the analysis. During a 6.5-year follow-up, 339 participants died, 510 experienced MACE, and 341 developed kidney failure. In fully adjusted models, NC was associated with all-cause death in women (HR 1.080 per cm; 95% CI 1.009–1.155), but not in men. Irrespective of sex, WC was associated with all-cause death (HR 1.014 per cm; 95% CI 1.005–1.038). NC and WC showed no association with MACE or kidney failure. BMI was not associated with any of the analyzed outcomes. Models of all-cause death including WC offered the best (lowest) AIC. Conclusion In Caucasian patients with mild to severe CKD, higher NC (in women) and WC were significantly associated with increased risk of death from any cause, but BMI was not.

The detection of contamination is key in all types of sequencing studies, especially for mitochondrial (mtDNA) studies. Contamination can lead to genotype misclassification and ultimately false-positive associations. In mtDNA studies, contamination is usually detected by applying computational expensive approaches using nuclear DNA (nDNA) reads. Here, we describe a novel approach that relies on the unique properties of the small mitochondrial genome to detect contamination fast and accurately.

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.

Biallelic mutations in PINK1 and PRKN cause recessively inherited Parkinson’s disease (PD). Though some studies suggest that PINK1/PRKN monoallelic mutations may not contribute to risk, deep phenotyping assessment showed that PINK1 or PRKN monoallelic pathogenic variants were at a significantly higher rate in PD compared to controls. Given the established role of PINK1 and Parkin in regulating mitochondrial dynamics, we explored mitochondrial DNA (mtDNA) integrity and inflammation as potential disease modifiers in carriers of mutations in these genes. MtDNA integrity, global gene expression and serum cytokine levels were investigated in a large collection of biallelic (n=84) and monoallelic (n=170) carriers of PINK1/PRKN mutations, iPD patients (n=67) and controls (n=90). Affected and unaffected PINK1/PRKN monoallelic mutation carriers can be distinguished by heteroplasmic mtDNA variant load (AUC=0.83, CI:0.74-0.93). Biallelic PINK1/PRKN mutation carriers harbor more heteroplasmic mtDNA variants in blood (p=0.0006, Z=3.63) compared to monoallelic mutation carriers. This enrichment was confirmed in iPSC-derived and postmortem midbrain neurons from biallelic PRKN-PD patients. Lastly, the heteroplasmic mtDNA variant load was found to correlate with IL6 levels in PINK1/PRKN mutation carriers (r=0.57, p=0.0074). PINK1/PRKN mutations predispose individuals to mtDNA variant accumulation in a dose- and disease-dependent manner. MtDNA variant load over time is a potential marker of disease manifestation in PINK1/PRKN mutation carriers.

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.

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.We aimed to investigate the association between apoA-IV and all-cause mortality and cardiovascular outcomes in the German Chronic Kidney Disease study.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.73mMean 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).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.

Heart-type fatty acid binding protein (H-FABP) is a biomarker that has been shown to provide long-term prognostic information in patients with coronary artery disease independently of high-sensitivity troponin T (hs-TNT). We examined the independent associations of H-FABP with cardiovascular outcomes in patients with chronic kidney disease (CKD).Prospective cohort study.4,951 patients enrolled in the German Chronic Kidney Disease (GCKD) study with an estimated glomerular filtration rate of 30-60 mL/min/1.73 mSerum levels of H-FABP and hs-TNT were measured at study entry.Noncardiovascular (non-CV) death, CV death, combined major adverse CV events (MACE), and hospitalization for congestive heart failure (CHF).Hazard ratios (HRs) for associations of H-FABP and hs-TNT with outcomes were estimated using Cox regression analyses adjusted for established risk factors.During a maximum follow-up of 6.5 years, 579 non-CV deaths, 190 CV deaths, 522 MACE, and 381 CHF hospitalizations were observed. In Cox regression analyses adjusted for established risk factors, H-FABP was associated with all 4 outcomes, albeit with lower HRs than those found for hs-TNT. After further adjustment for hs-TNT levels, H-FABP was found to be associated with non-CV death (HR, 1.57 [95% CI, 1.14-2.18]) and MACE (HR, 1.40 [95% CI, 1.02-1.92]) but with neither CV death (HR, 1.64 [95% CI, 0.90-2.99]) nor CHF hospitalizations (HR, 1.02 [95% CI, 0.70-1.49]).Single-point measurements of H-FABP and hs-TNT. Uncertain generalizability to non-European populations.In this large cohort of patients with CKD, H-FABP was associated with non-CV death and MACE, even after adjustment for hs-TNT. Whether measurement of H-FABP improves cardiovascular disease risk prediction in these patients warrants further studies.

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.

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.

Simple Summary In this study, we evaluated the prognostic value of the positron emission tomography–computed tomography (PET-CT) imaging technique in patients with newly diagnosed oral squamous cell carcinomas. PET-CT is routinely used to detect and quantify metabolically active tissues such as tumors. By using receiver operating characteristic (ROC) analysis, we successfully determined an optimal cut-off value for patient stratification in order to predict clinical outcome in this population. Furthermore, other clinical variables and their impact on clinical outcome as well as PET-CT values were evaluated. We show that, based on the determined optimal cut-off value, PET-CT is a reliable and independent predictor for clinical outcome, even in a fully adjusted model. Finally, we analyzed mitochondrial DNA to evaluate if potentially deleterious mutations might be a potential cause of metabolic changes, leading to differences in PET-CT values and consequently, clinical outcome. Abstract 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.

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.

Rationale & Objective Mechanisms underlying the variable course of disease progression in patients with chronic kidney disease (CKD) are incompletely understood. The aim of this study was to identify novel biomarkers of adverse kidney outcomes and overall mortality, which may offer insights into pathophysiologic mechanisms. Study Design Metabolome-wide association study. Setting & Participants 5,087 patients with CKD enrolled in the observational German Chronic Kidney Disease Study. Exposures Measurements of 1,487 metabolites in urine. Outcomes End points of interest were time to kidney failure (KF), a combined end point of KF and acute kidney injury (KF+AKI), and overall mortality. Analytical Approach Statistical analysis was based on a discovery-replication design (ratio 2:1) and multivariable-adjusted Cox regression models. Results After a median follow-up of 4 years, 362 patients died, 241 experienced KF, and 382 experienced KF+AKI. Overall, we identified 55 urine metabolites whose levels were significantly associated with adverse kidney outcomes and/or mortality. Higher levels of C-glycosyltryptophan were consistently associated with all 3 main end points (hazard ratios of 1.43 [95% CI, 1.27-1.61] for KF, 1.40 [95% CI, 1.27-1.55] for KF+AKI, and 1.47 [95% CI, 1.33-1.63] for death). Metabolites belonging to the phosphatidylcholine pathway showed significant enrichment. Members of this pathway contributed to the improvement of the prediction performance for KF observed when multiple metabolites were added to the well-established Kidney Failure Risk Equation. Limitations Findings among patients of European ancestry with CKD may not be generalizable to the general population. Conclusions Our comprehensive screen of the association between urine metabolite levels and adverse kidney outcomes and mortality identifies metabolites that predict KF and represents a valuable resource for future studies of biomarkers of CKD progression.

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% &gt, 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% &gt, 20 variant in mtDNA protein-coding regions relative to those with no mutations (HR = 4.6, (95%CI = 1.3&ndash, 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.

Within-species contamination is a major issue in sequencing studies, especially for mitochondrial studies. Contamination can be detected by analysing the nuclear genome or by inspecting the heteroplasmic sites in the mitochondrial genome. Existing methods using the nuclear genome are computationally expensive, and no suitable tool for detecting contamination in large-scale mitochondrial datasets is available. Here we present haplocheck, a tool that requires only the mitochondrial genome to detect contamination in both mitochondrial and whole-genome sequencing studies. Haplocheck is able to distinguish between contaminated and real heteroplasmic sites using the mitochondrial phylogeny. By applying haplocheck to the 1000 Genomes Project data, we show (1) high concordance in contamination estimates between mitochondrial and nuclear DNA and (2) quantify the impact of mitochondrial copy numbers on the mitochondrial based contamination results. Haplocheck complements leading nuclear DNA based contamination tools, and can therefore be used as a proxy tool in nuclear genome studies.Haplocheck is available both as a command-line tool athttps://github.com/genepi/haplocheckand as a cloud web-service producing interactive reports that facilitates the navigation through the phylogeny of contaminated samples.

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 breakthrough article published in PNAS by Luo et al. (2018) challenges a central dogma in biology which states that the mitochondrial DNA (mtDNA) is inherited exclusively from the mother. By sequencing the mitogenomes of several members of three independent families, the authors inferred an unprecedented pattern of biparental inheritance that requires the participation of an autosomal nuclear factor in the molecular process. However, a comprehensive analysis of their data reveals a number of issues that must be carefully addressed before challenging the current paradigm. Unfortunately, the methods section lacks any description of sample management, validation of their results in independent laboratories was deficient, and the reported findings have been observed at a frequency at complete variance with established evidence. Moreover, the remarkably high (and unusually homogeneous) levels of heteroplasmy reported can be readily detected using classical techniques for DNA sequencing. By reassessing the raw sequencing data with an alternative computational pipeline, we report strong correlation to the NextGENe results provided by the authors on a per sample base. However, the sequencing replicates from the same donors show aberrations in the variants detected that need further investigation to exclude contributions from other sources or methodological artifacts. Finally, applying the principle of reductio ad absurdum, we demonstrate that the nuclear factor invoked by the authors would need to be extraordinarily complex and precise in order to preclude linear accumulation of mtDNA lineages across generations. We discuss alternate scenarios that explain findings of the same nature as reported by Luo et al., in the context of in-vitro fertilization and therapeutic mtDNA replacement ooplasmic transplantation.

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.

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 (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 timeconsuming. 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 The project was supported by the Austrian Cancer Society/Tirol and by the Tyrolean Science Fund (Tiroler Wissenschaftsfonds) SI

Background Chronic kidney disease (CKD) affects 10–15% of the general population and affected individuals are at an increased risk for cardiovascular disease (CVD). Since telomere length is considered to be involved in biological aging, we tested whether relative telomere length (RTL) might be a marker for these two diseases. Methods The German Chronic Kidney Disease (GCKD) study is an ongoing prospective cohort study including patients with CKD of moderate severity. RTL was measured by qPCR in 4955 out of 5217 GCKD patients at baseline. Results RTL was distributed in the cohort with a mean ± SD of 0.95 ± 0.19. CVD was present in 1266 patients. Each decrease of RTL by 0.1 unit was associated with a higher probability for prevalent CVD: OR = 1.06, 95% CI 1.02–1.11, p = 0.007 (adjusted for age, sex, eGFR, BMI, ln-CRP, smoking, hypertension, diabetes, and lipids). Similar findings were observed for history of specific CVD entities, such as coronary artery disease (OR = 1.05, p = 0.025), myocardial infarction (OR = 1.08, p = 0.013) and percutaneous transluminal coronary angioplasty (OR = 1.06, p = 0.032). The strongest associations were found for interventions at the carotid arteries (OR = 1.25, p = 0.001) as well as aortic aneurysms (OR = 1.22, p = 0.001). Conclusions In the presence of CKD there is a significant association between shorter RTL and CVD manifestations. RTL appears to be a marker reflecting changes in homeostasis associated with CKD that may contribute to the excess CVD risk.