Comprehensive association testing of common mitochondrial DNA variation in metabolic disease

Many lines of evidence implicate mitochondria in phenotypic variation: (a) rare mutations in mitochondrial proteins cause metabolic, neurological, and muscular disorders; (b) alterations in oxidative phosphorylation are characteristic of type 2 diabetes, Parkinson disease, Huntington disease, and other diseases; and (c) common missense variants in the mitochondrial genome (mtDNA) have been implicated as having been subject to natural selection for adaptation to cold climates and contributing to "energy deficiency" diseases today. To test the hypothesis that common mtDNA variation influences human physiology and disease, we identified all 144 variants with frequency >1% in Europeans from >900 publicly available European mtDNA sequences and selected 64 tagging single-nucleotide polymorphisms that efficiently capture all common variation (except the hypervariable D-loop). Next, we evaluated the complete set of common mtDNA variants for association with type 2 diabetes in a sample of 3,304 diabetics and 3,304 matched nondiabetic individuals. Association of mtDNA variants with other metabolic traits (body mass index, measures of insulin secretion and action, blood pressure, and cholesterol) was also tested in subsets of this sample. We did not find a significant association of common mtDNA variants with these metabolic phenotypes. Moreover, we failed to identify any physiological effect of alleles that were previously proposed to have been adaptive for energy metabolism in human evolution. More generally, this comprehensive association-testing framework can readily be applied to other diseases for which mitochondrial dysfunction has been implicated.     

The First Modern Human Dispersals across Africa

The emergence of more refined chronologies for climate change and archaeology in prehistoric Africa, and for the evolution of human mitochondrial DNA (mtDNA), now make it feasible to test more sophisticated models of early modern human dispersals suggested by mtDNA distributions. Here we have generated 42 novel whole-mtDNA genomes belonging to haplogroup L0, the most divergent clade in the maternal line of descent, and analysed them alongside the growing database of African lineages belonging to L0’s sister clade, L1’6. We propose that the last common ancestor of modern human mtDNAs (carried by “mitochondrial Eve”) possibly arose in central Africa ~180 ka, at a time of low population size. By ~130 ka two distinct groups of anatomically modern humans co-existed in Africa: broadly, the ancestors of many modern-day Khoe and San populations in the south and a second central/eastern African group that includes the ancestors of most extant worldwide populations. Early modern human dispersals correlate with climate changes, particularly the tropical African “megadroughts” of MIS 5 (marine isotope stage 5, 135–75 ka) which paradoxically may have facilitated expansions in central and eastern Africa, ultimately triggering the dispersal out of Africa of people carrying haplogroup L3 ~60 ka. Two south to east migrations are discernible within haplogroup LO. One, between 120 and 75 ka, represents the first unambiguous long-range modern human dispersal detected by mtDNA and might have allowed the dispersal of several markers of modernity. A second one, within the last 20 ka signalled by L0d, may have been responsible for the spread of southern click-consonant languages to eastern Africa, contrary to the view that these eastern examples constitute relicts of an ancient, much wider distribution.     

Mitochondrial DNA Variant Discovery and Evaluation in Human Cardiomyopathies through Next-Generation Sequencing

Mutations in mitochondrial DNA (mtDNA) may cause maternally-inherited cardiomyopathy and heart failure. In homoplasmy all mtDNA copies contain the mutation. In heteroplasmy there is a mixture of normal and mutant copies of mtDNA. The clinical phenotype of an affected individual depends on the type of genetic defect and the ratios of mutant and normal mtDNA in affected tissues. We aimed at determining the sensitivity of next-generation sequencing compared to Sanger sequencing for mutation detection in patients with mitochondrial cardiomyopathy. We studied 18 patients with mitochondrial cardiomyopathy and two with suspected mitochondrial disease. We “shotgun” sequenced PCR-amplified mtDNA and multiplexed using a single run on Roche''s 454 Genome Sequencer. By mapping to the reference sequence, we obtained 1,300× average coverage per case and identified high-confidence variants. By comparing these to >400 mtDNA substitution variants detected by Sanger, we found 98% concordance in variant detection. Simulation studies showed that >95% of the homoplasmic variants were detected at a minimum sequence coverage of 20× while heteroplasmic variants required >200× coverage. Several Sanger “misses” were detected by 454 sequencing. These included the novel heteroplasmic 7501T>C in tRNA serine 1 in a patient with sudden cardiac death. These results support a potential role of next-generation sequencing in the discovery of novel mtDNA variants with heteroplasmy below the level reliably detected with Sanger sequencing. We hope that this will assist in the identification of mtDNA mutations and key genetic determinants for cardiomyopathy and mitochondrial disease.     

Heteroplasmic mitochondrial DNA variants in cardiovascular diseases

Mitochondria are implicated in the pathogenesis of cardiovascular diseases (CVDs) but the reasons for this are not well understood. Maternally-inherited population variants of mitochondrial DNA (mtDNA) which affect all mtDNA molecules (homoplasmic) are associated with cardiometabolic traits and the risk of developing cardiovascular disease. However, it is not known whether mtDNA mutations only affecting a proportion of mtDNA molecules (heteroplasmic) also play a role. To address this question, we performed a high-depth (~1000-fold) mtDNA sequencing of blood DNA in 1,399 individuals with hypertension (HTN), 1,946 with ischemic heart disease (IHD), 2,146 with ischemic stroke (IS), and 723 healthy controls. We show that the per individual burden of heteroplasmic single nucleotide variants (mtSNVs) increases with age. The age-effect was stronger for low-level heteroplasmies (heteroplasmic fraction, HF, 5–10%), likely reflecting acquired somatic events based on trinucleotide mutational signatures. After correcting for age and other confounders, intermediate heteroplasmies (HF 10–95%) were more common in hypertension, particularly involving non-synonymous variants altering the amino acid sequence of essential respiratory chain proteins. These findings raise the possibility that heteroplasmic mtSNVs play a role in the pathophysiology of hypertension.     

Allelic Spectra of Risk SNPs Are Different for Environment/Lifestyle Dependent versus Independent Diseases

Genome-wide association studies (GWAS) have generated sufficient data to assess the role of selection in shaping allelic diversity of disease-associated SNPs. Negative selection against disease risk variants is expected to reduce their frequencies making them overrepresented in the group of minor (<50%) alleles. Indeed, we found that the overall proportion of risk alleles was higher among alleles with frequency <50% (minor alleles) compared to that in the group of major alleles. We hypothesized that negative selection may have different effects on environment (or lifestyle)-dependent versus environment (or lifestyle)-independent diseases. We used an environment/lifestyle index (ELI) to assess influence of environmental/lifestyle factors on disease etiology. ELI was defined as the number of publications mentioning “environment” or “lifestyle” AND disease per 1,000 disease-mentioning publications. We found that the frequency distributions of the risk alleles for the diseases with strong environmental/lifestyle components follow the distribution expected under a selectively neutral model, while frequency distributions of the risk alleles for the diseases with weak environmental/lifestyle influences is shifted to the lower values indicating effects of negative selection. We hypothesized that previously selectively neutral variants become risk alleles when environment changes. The hypothesis of ancestrally neutral, currently disadvantageous risk-associated alleles predicts that the distribution of risk alleles for the environment/lifestyle dependent diseases will follow a neutral model since natural selection has not had enough time to influence allele frequencies. The results of our analysis suggest that prediction of SNP functionality based on the level of evolutionary conservation may not be useful for SNPs associated with environment/lifestyle dependent diseases.     

Bioenergetics in human evolution and disease: implications for the origins of biological complexity and the missing genetic variation of common diseases

Two major inconsistencies exist in the current neo-Darwinian evolutionary theory that random chromosomal mutations acted on by natural selection generate new species. First, natural selection does not require the evolution of ever increasing complexity, yet this is the hallmark of biology. Second, human chromosomal DNA sequence variation is predominantly either neutral or deleterious and is insufficient to provide the variation required for speciation or for predilection to common diseases. Complexity is explained by the continuous flow of energy through the biosphere that drives the accumulation of nucleic acids and information. Information then encodes complex forms. In animals, energy flow is primarily mediated by mitochondria whose maternally inherited mitochondrial DNA (mtDNA) codes for key genes for energy metabolism. In mammals, the mtDNA has a very high mutation rate, but the deleterious mutations are removed by an ovarian selection system. Hence, new mutations that subtly alter energy metabolism are continuously introduced into the species, permitting adaptation to regional differences in energy environments. Therefore, the most phenotypically significant gene variants arise in the mtDNA, are regional, and permit animals to occupy peripheral energy environments where rarer nuclear DNA (nDNA) variants can accumulate, leading to speciation. The neutralist–selectionist debate is then a consequence of mammals having two different evolutionary strategies: a fast mtDNA strategy for intra-specific radiation and a slow nDNA strategy for speciation. Furthermore, the missing genetic variation for common human diseases is primarily mtDNA variation plus regional nDNA variants, both of which have been missed by large, inter-population association studies.     

Mitochondrial haplogroup T is negatively associated with the status of elite endurance athlete

Mitochondrial function is absolutely necessary to supply the energy required for muscles, and germ line mutations in mitochondrial genes have been related with impaired cardiac function and exercise intolerance. In addition, alleles at several polymorphic sites in mtDNA define nine common haplogroups, and some of these haplogroups have been implicated in the risk of developing several diseases. In this study, we analysed the association between mtHaplogroups and the capacity to reach the status of elite endurance athlete. DNA was obtained from blood leukocytes of 95 Spanish elite endurance athletes and 250 healthy male population controls. We analysed eight mitochondrial polymorphisms and the frequencies were statistically compared between elite athletes and controls. Haplogroup T, specifically defined by 13368A, was significantly less frequent among elite endurance athletes (p =0.012, Fisher's exact test). Our study suggests that allele 13368A and mitochondrial haplogroup T might be a marker negatively associated with the status of elite endurance athlete. This mitochondrial variant could be related with a lower capacity to respond to endurance training, through unknown mechanisms involving a less efficient mitochondrial workload.     

Yeast Cells Expressing the Human Mitochondrial DNA Polymerase Reveal Correlations between Polymerase Fidelity and Human Disease Progression

Mutations in the human mitochondrial polymerase (polymerase-γ (Pol-γ)) are associated with various mitochondrial disorders, including mitochondrial DNA (mtDNA) depletion syndrome, Alpers syndrome, and progressive external opthamalplegia. To correlate biochemically quantifiable defects resulting from point mutations in Pol-γ with their physiological consequences, we created “humanized” yeast, replacing the yeast mtDNA polymerase (MIP1) with human Pol-γ. Despite differences in the replication and repair mechanism, we show that the human polymerase efficiently complements the yeast mip1 knockouts, suggesting common fundamental mechanisms of replication and conserved interactions between the human polymerase and other components of the replisome. We also examined the effects of four disease-related point mutations (S305R, H932Y, Y951N, and Y955C) and an exonuclease-deficient mutant (D198A/E200A). In haploid cells, each mutant results in rapid mtDNA depletion, increased mutation frequency, and mitochondrial dysfunction. Mutation frequencies measured in vivo equal those measured with purified enzyme in vitro. In heterozygous diploid cells, wild-type Pol-γ suppresses mutation-associated growth defects, but continuous growth eventually leads to aerobic respiration defects, reduced mtDNA content, and depolarized mitochondrial membranes. The severity of the Pol-γ mutant phenotype in heterozygous diploid humanized yeast correlates with the approximate age of disease onset and the severity of symptoms observed in humans.     

Disruption of Mitochondrial DNA Replication in Drosophila Increases Mitochondrial Fast Axonal Transport In Vivo

Mutations in mitochondrial DNA polymerase (pol γ) cause several progressive human diseases including Parkinson''s disease, Alper''s syndrome, and progressive external ophthalmoplegia. At the cellular level, disruption of pol γ leads to depletion of mtDNA, disrupts the mitochondrial respiratory chain, and increases susceptibility to oxidative stress. Although recent studies have intensified focus on the role of mtDNA in neuronal diseases, the changes that take place in mitochondrial biogenesis and mitochondrial axonal transport when mtDNA replication is disrupted are unknown. Using high-speed confocal microscopy, electron microscopy and biochemical approaches, we report that mutations in pol γ deplete mtDNA levels and lead to an increase in mitochondrial density in Drosophila proximal nerves and muscles, without a noticeable increase in mitochondrial fragmentation. Furthermore, there is a rise in flux of bidirectional mitochondrial axonal transport, albeit with slower kinesin-based anterograde transport. In contrast, flux of synaptic vesicle precursors was modestly decreased in pol γ−α mutants. Our data indicate that disruption of mtDNA replication does not hinder mitochondrial biogenesis, increases mitochondrial axonal transport, and raises the question of whether high levels of circulating mtDNA-deficient mitochondria are beneficial or deleterious in mtDNA diseases.     

Mitochondrial regulation of epigenetics and its role in human diseases

Most pathogenic mitochondrial DNA (mtDNA) mutations induce defects in mitochondrial oxidative phosphorylation (OXPHOS). However, phenotypic effects of these mutations show a large degree of variation depending on the tissue affected. These differences are difficult to reconcile with OXPHOS as the sole pathogenic factor suggesting that additional mechanisms contribute to lack of genotype and clinical phenotype correlationship. An increasing number of studies have identified a possible effect on the epigenetic landscape of the nuclear genome as a consequence of mitochondrial dysfunction. In particular, these studies demonstrate reversible or irreversible changes in genomic DNA methylation profiles of the nuclear genome. Here we review how mitochondria damage checkpoint (mitocheckpoint) induces epigenetic changes in the nucleus. Persistent pathogenic mutations in mtDNA may also lead to epigenetic changes causing genomic instability in the nuclear genome. We propose that “mitocheckpoint” mediated epigenetic and genetic changes may play key roles in phenotypic variation related to mitochondrial diseases or host of human diseases in which mitochondrial defect plays a primary role.     

Derived SNP alleles are used more frequently than ancestral alleles as risk-associated variants in common human diseases

Evolutionary aspects of the genetic architecture of common human diseases remain enigmatic. The results of more than 200 genome-wide association studies published to date were compiled in a catalog (). We used cataloged data to determine whether derived (mutant) alleles are associated with higher risk of human disease more frequently than ancestral alleles. We placed all allelic variants into ten categories of population frequency (0%-100%) in 10% increments. We then analyzed the relationship between allelic frequency, evolutionary status of the polymorphic site (ancestral versus derived), and disease risk status (risk versus protection). Given the same population frequency, derived alleles are more likely to be risk associated than ancestral alleles, as are rarer alleles. The common interpretation of this association is that negative selection prevents fixation of the risk variants. However, disease stratification as early or late onset suggests that weak selection against risk-associated alleles is unlikely a major factor shaping genetic architecture of common diseases. Our results clearly suggest that the duration of existence of an allele in a population is more important. Alleles existing longer tend to show weaker linkage disequilibrium with neighboring alleles, including the causal alleles, and are less likely to tag a SNP-disease association.     

Survival of the Curviest: Noise-Driven Selection for Synergistic Epistasis

A major goal of human genetics is to elucidate the genetic architecture of human disease, with the goal of fueling improvements in diagnosis and the understanding of disease pathogenesis. The degree to which epistasis, or non-additive effects of risk alleles at different loci, accounts for common disease traits is hotly debated, in part because the conditions under which epistasis evolves are not well understood. Using both theory and evolutionary simulation, we show that the occurrence of common diseases (i.e. unfit phenotypes with frequencies on the order of 1%) can, under the right circumstances, be expected to be driven primarily by synergistic epistatic interactions. Conditions that are necessary, collectively, for this outcome include a strongly non-linear phenotypic landscape, strong (but not too strong) selection against the disease phenotype, and “noise” in the genotype-phenotype map that is both environmental (extrinsic, time-correlated) and developmental (intrinsic, uncorrelated) and, in both cases, neither too little nor too great. These results suggest ways in which geneticists might identify, a priori, those disease traits for which an “epistatic explanation” should be sought, and in the process better focus ongoing searches for risk alleles.     

Association of Mitochondrial Genetic Variation with Carotid Atherosclerosis

In human pathology, several diseases are associated with somatic mutations in the mitochondrial genome (mtDNA). Even though mitochondrial dysfunction leads to increased oxidative stress, the role of mitochondrial mutations in atherosclerosis has not received much attention so far. In this study we analyzed the association of mitochondrial genetic variation with the severity of carotid atherosclerosis, as assessed by carotid intima-media thickness (cIMT) and the presence of coronary heart disease (CHD) in 190 subjects from Moscow, Russia, a population with high CHD occurrence. cIMT was measured by high-resolution B-mode ultrasonography and mtDNA heteroplasmies by a pyrosequencing-based method. We found that heteroplasmies for several mutations in the mtDNA in leukocytes, including C3256T, T3336C, G12315A, G13513A, G14459A, G14846A, and G15059A mutations, were significantly (p<0.001) associated with both the severity of carotid atherosclerosis and the presence of CHD. These findings indicate that somatic mitochondrial mutations have a role in the development of atherosclerosis.     

Somatic mtDNA Mutation Spectra in the Aging Human Putamen

The accumulation of heteroplasmic mitochondrial DNA (mtDNA) deletions and single nucleotide variants (SNVs) is a well-accepted facet of the biology of aging, yet comprehensive mutation spectra have not been described. To address this, we have used next generation sequencing of mtDNA-enriched libraries (Mito-Seq) to investigate mtDNA mutation spectra of putamen from young and aged donors. Frequencies of the “common” deletion and other “major arc” deletions were significantly increased in the aged cohort with the fold increase in the frequency of the common deletion exceeding that of major arc deletions. SNVs also increased with age with the highest rate of accumulation in the non-coding control region which contains elements necessary for translation and replication. Examination of predicted amino acid changes revealed a skew towards pathogenic SNVs in the coding region driven by mutation bias. Levels of the pathogenic m.3243A>G tRNA mutation were also found to increase with age. Novel multimeric tandem duplications that resemble murine control region multimers and yeast ρ⁻ mtDNAs, were identified in both young and aged specimens. Clonal ∼50 bp deletions in the control region were found at high frequencies in aged specimens. Our results reveal the complex manner in which the mitochondrial genome alters with age and provides a foundation for studies of other tissues and disease states.     

Mitochondrial nt3010G-nt3970C haplotype is implicated in high-altitude adaptation of Tibetans

Tibetans are well adapted to living and thriving in high-altitude environments. Mitochondria are central links to oxygen consumption, and variations in mitochondrial DNA (mtDNA) could play a role in high-altitude adaptation. Alleles at several polymorphic sites in mtDNA define common haplotypes, or haplogroups, including some that have been implicated in the risk of developing certain diseases. However, few reports have determined whether relationships exist between haplogroups and high-altitude adaptation in the Tibetan population. The D4 haplogroup is a major haplogroup of the Han Chinese. In the present study, genotypes of 12 polymorphisms were determined in members of a Tibetan population (n = 72), low altitude-Han (la-Han, n = 144), and high altitude-Han (ha-Han, n = 227) populations using polymerase chain reaction-restriction fragment length polymorphism and polymerase chain reaction-ligase detection reaction assays. The mitochondrial haplogroup D4 was negatively associated with high-altitude adaptation in Tibetans (P = 0.001 vs. la-Han, OR = 0.166, 95% CI = 0.048-0.567; P = 0.009 vs. ha-Han OR = 0.232, 95% CI = 0.069-0.778). The frequency of the nt3010G-nt3970C haplotype was significantly higher in Tibetans than in la-Han (P = 0.000) and ha-Han (P = 0.001) subjects. Findings in the present study suggest that unique mitochondrial variations determine a genetic background that is associated with high-altitude adaptation in the Tibetan population.     

Assessing Mitochondrial DNA Variation and Copy Number in Lymphocytes of ~2,000 Sardinians Using Tailored Sequencing Analysis Tools

DNA sequencing identifies common and rare genetic variants for association studies, but studies typically focus on variants in nuclear DNA and ignore the mitochondrial genome. In fact, analyzing variants in mitochondrial DNA (mtDNA) sequences presents special problems, which we resolve here with a general solution for the analysis of mtDNA in next-generation sequencing studies. The new program package comprises 1) an algorithm designed to identify mtDNA variants (i.e., homoplasmies and heteroplasmies), incorporating sequencing error rates at each base in a likelihood calculation and allowing allele fractions at a variant site to differ across individuals; and 2) an estimation of mtDNA copy number in a cell directly from whole-genome sequencing data. We also apply the methods to DNA sequence from lymphocytes of ~2,000 SardiNIA Project participants. As expected, mothers and offspring share all homoplasmies but a lesser proportion of heteroplasmies. Both homoplasmies and heteroplasmies show 5-fold higher transition/transversion ratios than variants in nuclear DNA. Also, heteroplasmy increases with age, though on average only ~1 heteroplasmy reaches the 4% level between ages 20 and 90. In addition, we find that mtDNA copy number averages ~110 copies/lymphocyte and is ~54% heritable, implying substantial genetic regulation of the level of mtDNA. Copy numbers also decrease modestly but significantly with age, and females on average have significantly more copies than males. The mtDNA copy numbers are significantly associated with waist circumference (p-value = 0.0031) and waist-hip ratio (p-value = 2.4×10^-5), but not with body mass index, indicating an association with central fat distribution. To our knowledge, this is the largest population analysis to date of mtDNA dynamics, revealing the age-imposed increase in heteroplasmy, the relatively high heritability of copy number, and the association of copy number with metabolic traits.     

Mitochondrial DNA Genetics and the Heteroplasmy Conundrum in Evolution and Disease

The unorthodox genetics of the mtDNA is providing new perspectives on the etiology of the common “complex” diseases. The maternally inherited mtDNA codes for essential energy genes, is present in thousands of copies per cell, and has a very high mutation rate. New mtDNA mutations arise among thousands of other mtDNAs. The mechanisms by which these “heteroplasmic” mtDNA mutations come to predominate in the female germline and somatic tissues is poorly understood, but essential for understanding the clinical variability of a range of diseases. Maternal inheritance and heteroplasmy also pose major challengers for the diagnosis and prevention of mtDNA disease.     

Genetic Dissection of the Drosophila melanogaster Female Head Transcriptome Reveals Widespread Allelic Heterogeneity

Modern genetic mapping is plagued by the “missing heritability” problem, which refers to the discordance between the estimated heritabilities of quantitative traits and the variance accounted for by mapped causative variants. One major potential explanation for the missing heritability is allelic heterogeneity, in which there are multiple causative variants at each causative gene with only a fraction having been identified. The majority of genome-wide association studies (GWAS) implicitly assume that a single SNP can explain all the variance for a causative locus. However, if allelic heterogeneity is prevalent, a substantial amount of genetic variance will remain unexplained. In this paper, we take a haplotype-based mapping approach and quantify the number of alleles segregating at each locus using a large set of 7922 eQTL contributing to regulatory variation in the Drosophila melanogaster female head. Not only does this study provide a comprehensive eQTL map for a major community genetic resource, the Drosophila Synthetic Population Resource, but it also provides a direct test of the allelic heterogeneity hypothesis. We find that 95% of cis-eQTLs and 78% of trans-eQTLs are due to multiple alleles, demonstrating that allelic heterogeneity is widespread in Drosophila eQTL. Allelic heterogeneity likely contributes significantly to the missing heritability problem common in GWAS studies.     

Genotype Imputation Reference Panel Selection Using Maximal Phylogenetic Diversity

The recent dramatic cost reduction of next-generation sequencing technology enables investigators to assess most variants in the human genome to identify risk variants for complex diseases. However, sequencing large samples remains very expensive. For a study sample with existing genotype data, such as array data from genome-wide association studies, a cost-effective approach is to sequence a subset of the study sample and then to impute the rest of the study sample, using the sequenced subset as a reference panel. The use of such an internal reference panel identifies population-specific variants and avoids the problem of a substantial mismatch in ancestry background between the study population and the reference population. To efficiently select an internal panel, we introduce an idea of phylogenetic diversity from mathematical phylogenetics and comparative genomics. We propose the “most diverse reference panel”, defined as the subset with the maximal “phylogenetic diversity”, thereby incorporating individuals that span a diverse range of genotypes within the sample. Using data both from simulations and from the 1000 Genomes Project, we show that the most diverse reference panel can substantially improve the imputation accuracy compared to randomly selected reference panels, especially for the imputation of rare variants. The improvement in imputation accuracy holds across different marker densities, reference panel sizes, and lengths for the imputed segments. We thus propose a novel strategy for planning sequencing studies on samples with existing genotype data.