Search

Background: Mitochondria play an important role in cellular metabolism, and their dysfunction is postulated to be involved in metabolic disturbances. Mitochondrial DNA is present in multiple copies per cell. The quantification of mitochondrial DNA copy number (mtDNA‐CN) might be used to assess mitochondrial dysfunction. Objectives: We aimed to investigate the cross‐sectional association of mtDNA‐CN with type 2 diabetes and the potential mediating role of metabolic syndrome. Methods: We examined 4812 patients from the German Chronic Kidney Disease (GCKD) study and 9364 individuals from the Cooperative Health Research in South Tyrol (CHRIS) study. MtDNA‐CN was measured in whole blood using a plasmid‐normalized qPCR‐based assay. Results: In both studies, mtDNA‐CN showed a significant correlation with most metabolic syndrome parameters: mtDNA‐CN decreased with increasing number of metabolic syndrome components. Furthermore, individuals with low mtDNA‐CN had significantly higher odds of metabolic syndrome (OR = 1.025; 95% CI = 1.011–1.039, P = 3.19 × 10−4, for each decrease of 10 mtDNA copies) and type 2 diabetes (OR = 1.027; 95% CI = 1.012–1.041; P = 2.84 × 10−4) in a model adjusted for age, sex, smoking and kidney function in the meta‐analysis of both studies. Mediation analysis revealed that the association of mtDNA‐CN with type 2 diabetes was mainly mediated by waist circumference in the GCKD study (66%) and by several metabolic syndrome parameters, especially body mass index and triglycerides, in the CHRIS study (41%). Conclusions: Our data show an inverse association of mtDNA‐CN with higher risk of metabolic syndrome and type 2 diabetes. A major part of the total effect of mtDNA‐CN on type 2 diabetes is mediated by obesity parameters. [ABSTRACT FROM AUTHOR]

Rationale & Objective: Afamin is a vitamin E-binding glycoprotein primarily expressed in liver and kidney. This study investigated whether serum afamin concentrations are associated with kidney function and incident kidney failure., Study Design: Prospective cohort study with 6.5 years follow-up., Setting & Participants: 5,041 Caucasian patients enrolled in the German Chronic Kidney Disease (GCKD) study with measured afamin concentrations and either an estimated glomerular filtration rate (eGFR) of 30-60 mL/min/1.73m 2 or an eGFR >60 mL/min/1.73m 2 with a urinary albumin to creatinine ratio (UACR) >300 mg/g at study entry., Exposure: Serum afamin concentrations (mg/L)., Outcome: Incident kidney failure (initiation of kidney replacement therapy or kidney-related death)., Analytical Approach: Generalized linear regression and quantile regression models fit to investigate the association of afamin concentrations with eGFR and UACR. Adjusted Cox regression analysis to examine the association of afamin concentrations with incident failure., Results: The mean (±SD) afamin concentration at study entry was 73.2±17.6 mg/L. Higher afamin concentrations were associated with better kidney function with a 2.60 ml/min/1.73m 2 higher eGFR (95% CI 2.30-2.89) and a 5.97 mg/g lower UACR (95% CI 3.04-8.90) for each 10 mg/L higher level of afamin concentration in adjusted analysis. During follow-up, each 10 mg/L higher level of afamin concentration was associated with a 14% lower risk of kidney failure (HR=0.86, 95%CI: 0.81 to 0.92, P<0.001)., Limitations: Residual confounding. Potential limited generalizability to non-Caucasian populations and people with mild stages of CKD or no CKD., Conclusions: Higher serum afamin concentrations appear to be associated with a higher eGFR, less albuminuria, and a lower risk for future kidney failure in patients with CKD., (Copyright © 2024. Published by Elsevier Inc.)

Background: 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.73 m2 and an overt proteinuria. Cardiovascular, cerebrovascular 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 2947 participants (56.5%), of which 680 participants (13.0%) died. 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 presumed diabetic kidney disease and men were more prone to experiencing events., Conclusion: This comprehensive explorative tool to visualize adverse events (https://www.gckd.org/studienhintergrund/previous-study-results/event-analysis/), their combination, mortality and multimorbidity among persons with CKD may serve as a valuable resourec for patient care, identification of high-risk groups, health services and public health policy planning., (© The Author(s) 2024. Published by Oxford University Press on behalf of the ERA.)

Background: Repetitive genome regions, such as variable number of tandem repeats (VNTR) or short tandem repeats (STR), are major constituents of the uncharted dark genome and evade conventional sequencing approaches. The protein-coding LPA kringle IV type-2 (KIV-2) VNTR (5.6 kb per unit, 1-40 units per allele) is a medically highly relevant example with a particularly intricate structure, multiple haplotypes, intragenic homologies, and an intra-VNTR STR. It is the primary regulator of plasma lipoprotein(a) [Lp(a)] concentrations, an important cardiovascular risk factor. Lp(a) concentrations vary widely between individuals and ancestries. Multiple variants and functional haplotypes in the LPA gene and especially in the KIV-2 VNTR strongly contribute to this variance., Methods: We evaluated the performance of amplicon-based nanopore sequencing with unique molecular identifiers (UMI-ONT-Seq) for SNP detection, haplotype mapping, VNTR unit consensus sequence generation, and copy number estimation via coverage-corrected haplotypes quantification in the KIV-2 VNTR. We used 15 human samples and low-level mixtures (0.5 to 5%) of KIV-2 plasmids as a validation set. We then applied UMI-ONT-Seq to extract KIV-2 VNTR haplotypes in 48 multi-ancestry 1000 Genome samples and analyzed at scale a poorly characterized STR within the KIV-2 VNTR., Results: UMI-ONT-Seq detected KIV-2 SNPs down to 1% variant level with high sensitivity, specificity, and precision (0.977 ± 0.018; 1.000 ± 0.0005; 0.993 ± 0.02) and accurately retrieved the full-length haplotype of each VNTR unit. Human variant levels were highly correlated with next-generation sequencing (R 2  = 0.983) without bias across the whole variant level range. Six reads per UMI produced sequences of each KIV-2 unit with Q40 quality. The KIV-2 repeat number determined by coverage-corrected unique haplotype counting was in close agreement with droplet digital PCR (ddPCR), with 70% of the samples falling even within the narrow confidence interval of ddPCR. We then analyzed 62,679 intra-KIV-2 STR sequences and explored KIV-2 SNP haplotype patterns across five ancestries., Conclusions: UMI-ONT-Seq accurately retrieves the SNP haplotype and precisely quantifies the VNTR copy number of each repeat unit of the complex KIV-2 VNTR region across multiple ancestries. This study utilizes the KIV-2 VNTR, presenting a novel and potent tool for comprehensive characterization of medically relevant complex genome regions at scale., (© 2024. The Author(s).)

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., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (© 2024. The Author(s).)

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.73m 2 or an eGFR > 60 mL/min/1.73m 2 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 eGFR creatinine ). 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., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

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.

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., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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 Author(s) 2022. Published by Oxford University Press on behalf of the Guarantors of Brain.)

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., Competing Interests: Conflict of interest F. K. has served on the advisory boards and has received lecture fees from Novartis, Amgen, and Kaneka. S. C., C. L., and L. F. have received lecture fees from Novartis. The other authors disclose no conflicts of interest., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

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., (© The Author(s) 2021. Published by Oxford University Press.)

Background: Sequencing quality has improved over the last decade for long-reads, allowing for more accurate detection of somatic low-frequency variants. In this study, we used mixtures of mitochondrial samples with different haplogroups (i.e., a specific set of mitochondrial variants) to investigate the applicability of nanopore sequencing for low-frequency single nucleotide variant detection. Methods: We investigated the impact of base-calling, alignment/mapping, quality control steps, and variant calling by comparing the results to a previously derived short-read gold standard generated on the Illumina NextSeq. For nanopore sequencing, six mixtures of four different haplotypes were prepared, allowing us to reliably check for expected variants at the predefined 5%, 2%, and 1% mixture levels. We used two different versions of Guppy for base-calling, two aligners (i.e., Minimap2 and Ngmlr), and three variant callers (i.e., Mutserve2, Freebayes, and Nanopanel2) to compare low-frequency variants. We used F 1 score measurements to assess the performance of variant calling. Results: We observed a mean read length of 11 kb and a mean overall read quality of 15. Ngmlr showed not only higher F 1 scores but also higher allele frequencies (AF) of false-positive calls across the mixtures (mean F 1 score = 0.83; false-positive allele frequencies < 0.17) compared to Minimap2 (mean F 1 score = 0.82; false-positive AF < 0.06). Mutserve2 had the highest F 1 scores (5% level: F 1 score >0.99, 2% level: F 1 score >0.54, and 1% level: F 1 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., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Lüth, Schaake, Grünewald, May, Trinh and Weissensteiner.)

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.73m 2 or an eGFR >60 ml/min/1.73m 2 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., (© 2021 The Authors. Journal of Internal Medicine published by John Wiley & Sons Ltd on behalf of Association for Publication of The Journal of Internal Medicine.)

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., (Copyright © 2021 Elsevier Inc. All rights reserved.)

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.

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.

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.

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., (© 2021 Weissensteiner et al.; Published by Cold Spring Harbor Laboratory Press.)

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.

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., (Copyright © 2020 International Society of Nephrology. Published by Elsevier Inc. All rights reserved.)

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

A breakthrough article published in PNAS by Luo et al. challenges a central dogma in biology which states that the mitochondrial DNA (mtDNA) in humans is inherited exclusively from the mother. We re-analyzed original FASTQ files and results reported by Luo et al. to investigate methodological issues (e.g. nuclear mitochondrial DNA or NUMTs, DNA rearrangements) that could lead to biological misinterpretations. A comprehensive analysis of their data reveals several methodological and analytical issues that must be carefully addressed before challenging the current paradigm. We first show that the probability of the findings described by the authors is extremely small (most likely below 10 -37 ). 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. Applying the principle of reductio ad absurdum, we demonstrate that the nuclear factor invoked by the authors to explain the phenomenon would need to be extraordinarily complex and precise to preclude linear accumulation of mtDNA lineages across generations, which would make the appearance of mixed haplotypes a much more frequent event in the population. 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., Competing Interests: Declaration of Competing Interest The authors declare that they have no conflicts of interest., (Copyright © 2020 Elsevier B.V. All rights reserved.)

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

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., (Copyright © 2019 International Society of Nephrology. Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2019 Coassin et al. Published by The American Society for Biochemistry and Molecular Biology, Inc.)

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., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

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., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Telomere length is considered as a biological marker for aging. It is expected that telomeres shorten with age and with conditions associated with oxidative stress and inflammation. Both are present in patients with chronic kidney disease (CKD) who have a very high cardiovascular risk. We investigated whether CKD duration is associated with relative telomere length (RTL) in 4802 patients from the German Chronic Kidney Disease (GCKD) study. We measured RTL in each sample in quadruplicates using a quantitative polymerase chain reaction (qPCR). We observed a U-shaped association of RTL with CKD duration: the longest RTL was found in those 339 patients who reported the shortest disease duration (<6 months) and shorter RTL in 2108 patients with duration between 6 months and less than 5 years. Most importantly, those 2331 patients who reported a CKD duration of 5 years and more had significantly longer RTL compared to those with intermediate CKD duration (6 months to less than 5 years): mean 0.954, 95%CI 0.946-0.961 versus 0.937, 95%CI 0.929-0.944, p=0.002). Due to the cross-sectional nature of the study these surprising results have to be considered with caution and as hypothesis-generating. Whether the longer RTL in patients with long-lasting disease is caused by an activation of telomerase to counteract the shortening of RTL due to oxidative stress and inflammation or whether they are caused by a survival bias needs to be investigated in longitudinal studies. Our data are in support of a higher plasticity of shortening and elongations of RTL as until recently anticipated., (Copyright © 2015 Elsevier Inc. All rights reserved.)

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., (Copyright © 2015 Elsevier Ireland Ltd. All rights reserved.)

Background: Oral squamous cell carcinoma (OSCC) is mainly caused by smoking and alcohol abuse and shows a five-year survival rate of ~50%. We aimed to explore the variation of somatic mitochondrial DNA (mtDNA) mutations in primary oral tumors, recurrences and metastases., Methods: We performed an in-depth validation of mtDNA next-generation sequencing (NGS) on an Illumina HiSeq 2500 platform for its application to cancer tissues, with the goal to detect low-level heteroplasmies and to avoid artifacts. Therefore we genotyped the mitochondrial genome (16.6 kb) from 85 tissue samples (tumors, recurrences, resection edges, metastases and blood) collected from 28 prospectively recruited OSCC patients applying both Sanger sequencing and high-coverage NGS (~35,000 reads per base)., Results: We observed a strong correlation between Sanger sequencing and NGS in estimating the mixture ratio of heteroplasmies (r = 0.99; p<0.001). Non-synonymous heteroplasmic variants were enriched among cancerous tissues. The proportions of somatic and inherited variants in a given gene region were strongly correlated (r = 0.85; p<0.001). Half of the patients shared mutations between benign and cancerous tissue samples. Low level heteroplasmies (<10%) were more frequent in benign samples compared to tumor samples, where heteroplasmies >10% were predominant. Four out of six patients who developed a local tumor recurrence showed mutations in the recurrence that had also been observed in the primary tumor. Three out of five patients, who had tumor metastases in the lymph nodes of their necks, shared mtDNA mutations between primary tumors and lymph node metastases. The percentage of mutation heteroplasmy increased from the primary tumor to lymph node metastases., Conclusions: We conclude that Sanger sequencing is valid for heteroplasmy quantification for heteroplasmies ≥10% and that NGS is capable of reliably detecting and quantifying heteroplasmies down to the 1%-level. The finding of shared mutations between primary tumors, recurrences and metastasis indicates a clonal origin of malignant cells in oral cancer.

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.

Objective: To evaluate the association between a family history of cardiovascular and metabolic diseases or risk factors and the presence of peripheral artery disease (PAD)., Methods: Participants were recruited within one PAD case-control study (CAVASIC, n = 481) and two population-based studies (KORA-F3, n = 3118; KORA-F4, n = 1325). In the KORA studies, an ankle-brachial-index <0.9 and/or symptomatic claudication was defined as PAD. For myocardial infarction, stroke, diabetes mellitus, hypertension, PAD, and obesity, family risk scores (FamRS) were calculated taking into account the number and age of diseased parents and siblings and regressed on prevalent PAD or ankle-brachial-index., Results: A significant association with PAD was found for family history of myocardial infarction and hypertension in a combined analysis of all studies and for family history of PAD in the case-control study. A combined family history score was derived from FamRS myocardial infarction, stroke, diabetes mellitus, hypertension and PAD. A positive family history of at least two and/or a strong positive family history of at least one of these diseases was found to be associated with increased risk of prevalent PAD in all three studies, even after adjustment for classical risk factors (OR = 1.93, 95%CI = [1.53-2.44], p = 2.6 × 10(-8))., Conclusion: The presence of a positive family history for cardiovascular and metabolic diseases or risk factors was shown to be associated with the presence of PAD. A FamRS calculation tool that considers age and family size can guide a physician to perform extended examinations to prevent complications and progression of PAD., (Copyright © 2014 Elsevier Ireland Ltd. All rights reserved.)

Background: Myanmar is the largest country in mainland Southeast Asia with a population of 55 million people subdivided into more than 100 ethnic groups. Ruled by changing kingdoms and dynasties and lying on the trade route between India and China, Myanmar was influenced by numerous cultures. Since its independence from British occupation, tensions between the ruling Bamar and ethnic minorities increased., Results: Our aim was to search for genetic footprints of Myanmar's geographic, historic and sociocultural characteristics and to contribute to the picture of human colonization by describing and dating of new mitochondrial DNA (mtDNA) haplogroups. Therefore, we sequenced the mtDNA control region of 327 unrelated donors and the complete mitochondrial genome of 44 selected individuals according to highest quality standards., Conclusion: Phylogenetic analyses of the entire mtDNA genomes uncovered eight new haplogroups and three unclassified basal M-lineages. The multi-ethnic population and the complex history of Myanmar were reflected in its mtDNA heterogeneity. Population genetic analyses of Burmese control region sequences combined with population data from neighboring countries revealed that the Myanmar haplogroup distribution showed a typical Southeast Asian pattern, but also Northeast Asian and Indian influences. The population structure of the extraordinarily diverse Bamar differed from that of the Karen people who displayed signs of genetic isolation. Migration analyses indicated a considerable genetic exchange with an overall positive migration balance from Myanmar to neighboring countries. Age estimates of the newly described haplogroups point to the existence of evolutionary windows where climatic and cultural changes gave rise to mitochondrial haplogroup diversification in Asia.