Genetic variations related to maternal whole blood mitochondrial DNA copy number: a genome-wide and candidate gene study

Nuclear and mitochondrial genetic variants associated with mitochondrial DNA copy number

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

The Association Between Polymorphisms in Cell-Cycle Genes and Mitochondrial DNA Copy Number in Coke Oven Workers

The mitochondrial DNA (mtDNA) copy number is a vital component in maintaining normal mitochondrial function. It is affected by environmental and occupational exposures, as well as polymorphisms in nuclear genes. Nonetheless, the specific roles of polymorphisms in cell-cycle genes and mtDNA copy number are still unknown. This study enrolled a sample of 544 coke oven workers and 238 non-exposed controls so as to assess the effect of exposure of coke oven emissions (COEs) and polymorphisms in cell-cycle genes on the mtDNA copy number. We found that the mtDNA copy number in the exposed group (0.60 ± 0.29) was significantly lower than that in the control group (1.03 ± 0.31) (t =18.931, P < 0.001). The analysis of covariance showed that both the rs1801270 (CA+CC) and the rs1059234 (CT+CC) in p21 gene were associated with lower mtDNA copy number in the exposed group (P = 0.001). Generalized linear models indicated COEs-exposure (β = -0.432, P < 0.001) and rs1059234 (CT+CC) in p21 gene (β = -0.060, P = 0.024) were the factors in mtDNA copy number reduction. In conclusion, this study suggests that the decrease of the mtDNA copy number is associated with COEs-exposure and the rs1059234 (CT+CC) in the p21 gene.

Influence of environmental factors and genetic variation on mitochondrial DNA copy number

Mitochondrial DNA copy number (mtDNA CN) has been shown to be highly heritable and associated with traits of interest in humans. However, studies are lacking in the literature for livestock species such as beef cattle. In this study, 2,371 individuals from a crossbred beef population comprising the Germplasm Evaluation program from the U.S. Meat Animal Research Center had samples of blood, leucocyte, or semen collected for low-pass sequencing (LPS) that resulted in both nuclear DNA (nuDNA) and mitochondrial DNA (mtDNA) sequence reads. Mitochondrial DNA CN was estimated based on the ratio of mtDNA to nuDNA coverages. Genetic parameters for mtDNA CN were estimated from an animal model based on a genomic relationship matrix (~87K SNP from the nuDNA). Different models were used to test the effects of tissue, sex, age at sample collection, heterosis, and breed composition. Maternal effects, assessed by fitting a maternal additive component and by fitting eleven SNP on the mtDNA, were also obtained. As previously reported, mtDNA haplotypes were used to classify individuals into Taurine haplogroups (T1, T2, T3/T4, and T5). Estimates of heritability when fitting fixed effects in addition to the intercept were moderate, ranging from 0.11 to 0.31 depending on the model. From a model ignoring contemporary group, semen samples had the lowest mtDNA CN, as expected, followed by blood and leucocyte samples (P ≤ 0.001). The effect of sex and the linear and quadratic effects of age were significant (P ≤ 0.02) depending on the model. When significant, females had greater mtDNA CN than males. The effects of heterosis and maternal heterosis were not significant (P ≥ 0.47). The estimates of maternal and mtDNA heritability were near zero (≤0.03). Most of the samples (98%) were classified as haplogroup T3. Variation was observed in the mtDNA within Taurine haplogroups, which enabled the identification of 24 haplotypes. These results suggest that mtDNA CN is under nuclear genetic control and would respond favorably to selection.

Genome-wide analysis of mitochondrial DNA copy number reveals loci implicated in nucleotide metabolism, platelet activation, and megakaryocyte proliferation

Mitochondrial DNA copy number (mtDNA-CN) measured from blood specimens is a minimally invasive marker of mitochondrial function that exhibits both inter-individual and intercellular variation. To identify genes involved in regulating mitochondrial function, we performed a genome-wide association study (GWAS) in 465,809 White individuals from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and the UK Biobank (UKB). We identified 133 SNPs with statistically significant, independent effects associated with mtDNA-CN across 100 loci. A combination of fine-mapping, variant annotation, and co-localization analyses was used to prioritize genes within each of the 133 independent sites. Putative causal genes were enriched for known mitochondrial DNA depletion syndromes (p = 3.09 × 10-15) and the gene ontology (GO) terms for mtDNA metabolism (p = 1.43 × 10-8) and mtDNA replication (p = 1.2 × 10-7). A clustering approach leveraged pleiotropy between mtDNA-CN associated SNPs and 41 mtDNA-CN associated phenotypes to identify functional domains, revealing three distinct groups, including platelet activation, megakaryocyte proliferation, and mtDNA metabolism. Finally, using mitochondrial SNPs, we establish causal relationships between mitochondrial function and a variety of blood cell-related traits, kidney function, liver function and overall (p = 0.044) and non-cancer mortality (p = 6.56 × 10-4).

Genome-wide association study of mitochondrial copy number

Mitochondrial DNA copy number (mtDNAcn) variation has been associated with increased risk of several human diseases in epidemiological studies. The quantification of mtDNAcn performed with real-time PCR is currently considered the de facto standard among several techniques. However, the heterogeneity of the laboratory methods (DNA extraction, storage, processing) used could give rise to results that are difficult to compare and reproduce across different studies. Several lines of evidence suggest that mtDNAcn is influenced by nuclear and mitochondrial genetic variability, however this relation is largely unexplored. The aim of this work was to elucidate the genetic basis of mtDNAcn variation. We performed a genome-wide association study (GWAS) of mtDNAcn in 6836 subjects from the ESTHER prospective cohort, and included, as replication set, the summary statistics of a GWAS that used 295 150 participants from the UK Biobank. We observed two novel associations with mtDNAcn variation on chromosome 19 (rs117176661), and 12 (rs7136238) that reached statistical significance at the genome-wide level. A polygenic score that we called mitoscore including all known single nucleotide polymorphisms explained 1.11% of the variation of mtDNAcn (p = 5.93 × 10-7). In conclusion, we performed a GWAS on mtDNAcn, adding to the evidence of the genetic background of this trait.

Deciphering the genetic and epidemiological landscape of mitochondrial DNA abundance

Mitochondrial (MT) dysfunction is a hallmark of aging and has been associated with most aging-related diseases as well as immunological processes. However, little is known about aging, lifestyle and genetic factors influencing mitochondrial DNA (mtDNA) abundance. In this study, mtDNA abundance was estimated from the weighted intensities of probes mapping to the MT genome in 295,150 participants from the UK Biobank. We found that the abundance of mtDNA was significantly elevated in women compared to men, was negatively correlated with advanced age, higher smoking exposure, greater body-mass index, higher frailty index as well as elevated red and white blood cell count and lower mortality. In addition, several biochemistry markers in blood-related to cholesterol metabolism, ion homeostasis and kidney function were found to be significantly associated with mtDNA abundance. By performing a genome-wide association study, we identified 50 independent regions genome-wide significantly associated with mtDNA abundance which harbour multiple genes involved in the immune system, cancer as well as mitochondrial function. Using mixed effects models, we estimated the SNP-heritability of mtDNA abundance to be around 8%. To investigate the consequence of altered mtDNA abundance, we performed a phenome-wide association study and found that mtDNA abundance is involved in risk for leukaemia, hematologic diseases as well as hypertension. Thus, estimating mtDNA abundance from genotyping arrays has the potential to provide novel insights into age- and disease-relevant processes, particularly those related to immunity and established mitochondrial functions.

A case for measuring both cellular and cell-free mitochondrial DNA as a disease biomarker in human blood

Circulating mitochondrial DNA (mtDNA), widely studied as a disease biomarker, comprises of mtDNA located within mitochondria, indicative of mitochondrial function, and cell-free (cf) mtDNA linked to inflammation. The purpose of this study was to determine the ranges of, and relationship between, cellular and cf mtDNA in human blood. Whole blood from 23 controls (HC) and 20 patients with diabetes was separated into peripheral blood mononuclear cells (PBMCs), plasma, and serum. Total DNA was isolated and mtDNA copy numbers were determined using absolute quantification. Cellular mtDNA content in PBMCs was higher than in peripheral blood and a surprisingly high level of cf mtDNA was present in serum and plasma of HC, with no direct relationship between cellular and cf mtDNA content within individuals. Diabetes patients had similar levels of cellular mtDNA compared to healthy participants but a significantly higher cf mtDNA content. Furthermore, only in patients with diabetes, we observed a correlation between whole blood and plasma mtDNA levels, indicating that the relationship between cellular and cf mtDNA content is affected by disease status. In conclusion, when evaluating mtDNA in human blood as a biomarker of mitochondrial dysfunction, it is important to measure both cellular and cf mtDNA.

A genome-wide association study of mitochondrial DNA copy number in two population-based cohorts

BACKGROUND

Mitochondrial DNA copy number (mtDNA CN) exhibits interindividual and intercellular variation, but few genome-wide association studies (GWAS) of directly assayed mtDNA CN exist. We undertook a GWAS of qPCR-assayed mtDNA CN in the Avon Longitudinal Study of Parents and Children (ALSPAC) and the UK Blood Service (UKBS) cohort. After validating and harmonising data, 5461 ALSPAC mothers (16-43 years at mtDNA CN assay) and 1338 UKBS females (17-69 years) were included in a meta-analysis. Sensitivity analyses restricted to females with white cell-extracted DNA and adjusted for estimated or assayed cell proportions. Associations were also explored in ALSPAC children and UKBS males.

RESULTS

A neutrophil-associated locus approached genome-wide significance (rs709591 [MED24], β (change in SD units of mtDNA CN per allele) [SE] - 0.084 [0.016], p = 1.54e-07) in the main meta-analysis of adult females. This association was concordant in magnitude and direction in UKBS males and ALSPAC neonates. SNPs in and around ABHD8 were associated with mtDNA CN in ALSPAC neonates (rs10424198, β [SE] 0.262 [0.034], p = 1.40e-14), but not other study groups. In a meta-analysis of unrelated individuals (N = 11,253), we replicated a published association in TFAM (β [SE] 0.046 [0.017], p = 0.006), with an effect size much smaller than that observed in the replication analysis of a previous in silico GWAS.

CONCLUSIONS

In a hypothesis-generating GWAS, we confirm an association between TFAM and mtDNA CN and present putative loci requiring replication in much larger samples. We discuss the limitations of our work, in terms of measurement error and cellular heterogeneity, and highlight the need for larger studies to better understand nuclear genomic control of mtDNA copy number.

Mutation in Cytochrome B gene causes debility and adverse effects on health of sheep

Cytochrome B is the mitochondrial protein, which functions as part of the electron transport chain and is the main subunit of transmembrane cytochrome bc1 and b6f complexes affecting energy metabolism through oxidative phosphorylation. The present study was conducted to study the effect of mutation of Cytochrome B gene on the health condition of sheep, which the first report of association of mitochondrial gene with disease traits in livestock species. Non-synonymous substitutions (F33 L and D171N) and Indel mutations were observed for Cytochrome B gene, leading to a truncated protein, where anemia, malfunctioning of most of the vital organs as liver, kidney and mineral status was observed and debility with exercise intolerance and cardiomyopathy in extreme cases were depicted. These findings were confirmed by bioinformatics analysis, haematological and biochemical data analysis, and other phenotypical physiological data pertaining to different vital organs. The molecular mechanism of cytochrome B mutation was that the mutant variant interferes with the site of heme binding (iron containing) domain and calcium binding essential for electron transport chain. Mutation at amino acid site 33 is located within transmembrane helix A, a hydrophobic environment at the Qi site and close to heme binding domain, and mutation effects these domain and diseases occur. Thermodynamic stability was also observed to decrease in mutant variant. Sheep Cytochrome B being genetically more similar to the human, it may be used as a model for studying human diseases related to cytochrome B defects. Future prospect of the study includes the therapeutic application of recombinant protein, gene therapy and marker-assisted selection of disease-resistant livestock.

Genetic variants in nuclear DNA along with environmental factors modify mitochondrial DNA copy number: a population-based exome-wide association study

Background

Mitochondrial DNA (mtDNA) copy number has been found associated with multiple diseases, including cancers, diabetes and so on. Both environmental and genetic factors could affect the copy number of mtDNA. However, limited study was available about the relationship between genetic variants and mtDNA copy number. What’s more, most of previous studies considered only environmental or genetic factors. Therefore, it’s necessary to explore the genetic effects on mtDNA copy number with the consideration of PM_2.5 exposure and smoking.

Results

A multi-center population-based study was performed with 301 subjects from Zhuhai, Wuhan and Tianjin. Personal 24-h PM_2.5 exposure levels, smoking and mtDNA copy number were evaluated. The Illumina Human Exome BeadChip, which contained 241,305 single nucleotide variants, was used for genotyping. The association analysis was conducted in each city and meta-analysis was adopted to combine the overall effect among three cities. Seven SNPs showed significant association with mtDNA copy number with P value less than 1.00E-04 after meta-analysis. The following joint analysis of our identified SNPs showed a significant allele-dosage association between the number of variants and mtDNA copy number (P = 5.02 × 10^− 17). Further, 11 genes were identified associated with mtDNA copy number using gene-based analysis with a P value less than 0.01.

Conclusion

This study was the first attempt to evaluate the genetic effects on mtDNA copy number with the consideration of personal PM_2.5 exposure level. Our findings could provide more evidences that genetic variants played important roles in modulating the copy number of mtDNA.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5142-7) contains supplementary material, which is available to authorized users.

New insights into mitogenomic phylogeny and copy number in eight indigenous sheep populations based on the ATP synthase and cytochrome c oxidase genes

The origins and phylogeny of different sheep breeds has been widely studied using polymorphisms within the mitochondrial hypervariable region. However, little is known about the mitochondrial DNA (mtDNA) content and phylogeny based on mtDNA protein-coding genes. In this study, we assessed the phylogeny and copy number of the mtDNA in eight indigenous (population size, n=184) and three introduced (n=66) sheep breeds in China based on five mitochondrial coding genes (COX1, COX2, ATP8, ATP6 and COX3). The mean haplotype and nucleotide diversities were 0.944 and 0.00322, respectively. We identified a correlation between the lineages distribution and the genetic distance, whereby Valley-type Tibetan sheep had a closer genetic relationship with introduced breeds (Dorper, Poll Dorset and Suffolk) than with other indigenous breeds. Similarly, the Median-joining profile of haplotypes revealed the distribution of clusters according to genetic differences. Moreover, copy number analysis based on the five mitochondrial coding genes was affected by the genetic distance combining with genetic phylogeny; we also identified obvious non-synonymous mutations in ATP6 between the different levels of copy number expressions. These results imply that differences in mitogenomic compositions resulting from geographical separation lead to differences in mitochondrial function.