Principal-component analysis for assessment of population stratification in mitochondrial medical genetics

Maternal genetic architecture of Guizhou's Hmong-Mien populations via whole mitogenome

The Hmong-Mien (HM) language family, majorly distributed across southern China and Southeast Asia, has remained underexplored in population genetics, particularly concerning whole mitogenome studies. In this study, we sequenced the whole mitogenomes of 261 individuals from Guizhou Hmong-Mien-speaking populations (HM-G), comprising Miao, Yao, and She individuals. The haplogroup distribution was dominated by southern East Asian haplogroups (B, M7, and F). The neutrality test revealed significantly negative values, and mismatch distribution analyses showed a pronounced unimodal distribution, indicating high genetic diversity and recent population expansion in populations. For a comprehensive understanding of the matrilineal genetic background of the HM-speaking population, we merged whole mitogenome data from 83 populations worldwide. Our findings showed that the HM-G exhibited relatively close genetic distances to HM- and Tai-Kadai-speaking populations from East and Southeast Asia. These results provide crucial insights into the genetic structure and evolutionary history of HM-G.

Mitochondrial-derived microproteins: from discovery to function

Given the uniqueness of the mitochondria, and the fact that they have their own genome, mitochondrial-derived microproteins (MDPs) are similar to, but different from, nuclear-encoded microproteins. The discovery of an increasing number of microproteins from this organelle and the importance of mitochondria to cellular and organismal health make it a priority to study this novel class of proteins in search of possible therapeutic targets and cures. In this review, we discuss the history of MDP discovery, describe the function of each MDP, and conclude with future goals and techniques to help discover more MDPs.

Whole mitochondrial genome analysis of Tai-Kadai-speaking populations in Southwest China

As a single matrilineal gene, human mitochondrial DNA plays a very important role in the study of population genetics. The whole mitogenome sequences of 287 individuals of the Tai-Kadai-speaking population in Guizhou were obtained. It was discovered that there were 82, 104, and 94 haplotypes in 83 Bouyei individuals, 107 Dong individuals, and 97 Sui individuals, respectively; and the haplotype diversity in Bouyei, Dong, and Sui groups was 1.000 ± 0.02, 0.9993 ± 0.0015, and 0.999 ± 0.002, respectively. The result of neutrality tests of the Tai-Kadai-speaking population in Guizhou showed significant negative values, and the analysis of mismatch distribution showed an obvious unimodal distribution. The results implied that Guizhou Tai-Kadai-speaking populations had high genetic diversities and may have experienced recent population expansion. In addition, the primary haplogroups of studied populations were M*, F, B, D, and R*, implying that they may origin from Southern China. The matrilineal genetic structure of the Tai-Kadai-speaking populations in Guizhou was analyzed by merging the mitogenome data of 79 worldwide populations as reference data. The results showed that there were close relationships between studied populations and other Tai-Kadai as well as some Austronesian populations in East and Southeast Asia. Overall, the mitogenome data generated in this study will provide important data for the study of genetic structure of Tai-Kadai speaking populations.

Mitochondria-derived peptides in aging and healthspan

The mechanisms that explain mitochondrial dysfunction in aging and healthspan continue to be studied, but one element has been unexplored: microproteins. Small open reading frames in circular mitochondria DNA can encode multiple microproteins, called mitochondria-derived peptides (MDPs). Currently, eight MDPs have been published: humanin, MOTS-c, and SHLPs 1–6. This Review describes recent advances in microprotein discovery with a focus on MDPs. It discusses what is currently known about MDPs in aging and how this new understanding could add to the way we understand age-related diseases including type 2 diabetes, cancer, and neurodegenerative diseases at the genomic, proteomic, and drug-development levels.

Mitochondrial genome-wide analysis of nuclear DNA methylation quantitative trait loci

Mitochondria have a complex communication network with the surrounding cell and can alter nuclear DNA methylation (DNAm). Variation in the mitochondrial DNA (mtDNA) has also been linked to differential DNAm. Genome-wide association studies have identified numerous DNAm quantitative trait loci, but these studies have not examined the mitochondrial genome. Herein, we quantified nuclear DNAm from blood and conducted a mitochondrial genome-wide association study of DNAm, with an additional emphasis on sex- and prediabetes-specific heterogeneity. We used the Young Finns Study (n = 926) with sequenced mtDNA genotypes as a discovery sample and sought replication in the Ludwigshafen Risk and Cardiovascular Health study (n = 2317). We identified numerous significant associations in the discovery phase (P < 10−9), but they were not replicated when accounting for multiple testing. In total, 27 associations were nominally replicated with a P < 0.05. The replication analysis presented no evidence of sex- or prediabetes-specific heterogeneity. The 27 associations were included in a joint meta-analysis of the two cohorts, and 19 DNAm sites associated with mtDNA variants, while four other sites showed haplogroup associations. An expression quantitative trait methylation analysis was performed for the identified DNAm sites, pinpointing two statistically significant associations. This study provides evidence of a mitochondrial genetic control of nuclear DNAm with little evidence found for sex- and prediabetes-specific effects. The lack of a comparable mtDNA data set for replication is a limitation in our study and further studies are needed to validate our results.

Human genetic analyses of organelles highlight the nucleus in age-related trait heritability

Most age-related human diseases are accompanied by a decline in cellular organelle integrity, including impaired lysosomal proteostasis and defective mitochondrial oxidative phosphorylation. An open question, however, is the degree to which inherited variation in or near genes encoding each organelle contributes to age-related disease pathogenesis. Here, we evaluate if genetic loci encoding organelle proteomes confer greater-than-expected age-related disease risk. As mitochondrial dysfunction is a 'hallmark' of aging, we begin by assessing nuclear and mitochondrial DNA loci near genes encoding the mitochondrial proteome and surprisingly observe a lack of enrichment across 24 age-related traits. Within nine other organelles, we find no enrichment with one exception: the nucleus, where enrichment emanates from nuclear transcription factors. In agreement, we find that genes encoding several organelles tend to be 'haplosufficient,' while we observe strong purifying selection against heterozygous protein-truncating variants impacting the nucleus. Our work identifies common variation near transcription factors as having outsize influence on age-related trait risk, motivating future efforts to determine if and how this inherited variation then contributes to observed age-related organelle deterioration.

Examining the effect of mitochondrial DNA variants on blood pressure in two Finnish cohorts

High blood pressure (BP) is a major risk factor for many noncommunicable diseases. The effect of mitochondrial DNA single-nucleotide polymorphisms (mtSNPs) on BP is less known than that of nuclear SNPs. We investigated the mitochondrial genetic determinants of systolic, diastolic, and mean arterial BP. MtSNPs were determined from peripheral blood by sequencing or with genome-wide association study SNP arrays in two independent Finnish cohorts, the Young Finns Study and the Finnish Cardiovascular Study, respectively. In total, over 4200 individuals were included. The effects of individual common mtSNPs, with an additional focus on sex-specificity, and aggregates of rare mtSNPs grouped by mitochondrial genes were evaluated by meta-analysis of linear regression and a sequence kernel association test, respectively. We accounted for the predicted pathogenicity of the rare variants within protein-encoding and the tRNA regions. In the meta-analysis of 87 common mtSNPs, we did not observe significant associations with any of the BP traits. Sex-specific and rare-variant analyses did not pinpoint any significant associations either. Our results are in agreement with several previous studies suggesting that mtDNA variation does not have a significant role in the regulation of BP. Future studies might need to reconsider the mechanisms thought to link mtDNA with hypertension.

MitoIMP: A Computational Framework for Imputation of Missing Data in Low-Coverage Human Mitochondrial Genome

The incompleteness of partial human mitochondrial genome sequences makes it difficult to perform relevant comparisons among multiple resources. To deal with this issue, we propose a computational framework for deducing missing nucleotides in the human mitochondrial genome. We applied it to worldwide mitochondrial haplogroup lineages and assessed its performance. Our approach can deduce the missing nucleotides with a precision of 0.99 or higher in most human mitochondrial DNA lineages. Furthermore, although low-coverage mitochondrial genome sequences often lead to a blurred relationship in the multidimensional scaling analysis, our approach can correct this positional arrangement according to the corresponding mitochondrial DNA lineages. Therefore, our framework will provide a practical solution to compensate for the lack of genome coverage in partial and fragmented human mitochondrial genome sequences. In this study, we developed an open-source computer program, MitoIMP, implementing our imputation procedure. MitoIMP is freely available from https://github.com/omics-tools/mitoimp.

Worldwide human mitochondrial haplogroup distribution from urban sewage

Community level genetic information can be essential to direct health measures and study demographic tendencies but is subject to considerable ethical and legal challenges. These concerns become less pronounced when analyzing urban sewage samples, which are ab ovo anonymous by their pooled nature. We were able to detect traces of the human mitochondrial DNA (mtDNA) in urban sewage samples and to estimate the distribution of human mtDNA haplogroups. An expectation maximization approach was used to determine mtDNA haplogroup mixture proportions for samples collected at each different geographic location. Our results show reasonable agreement with both previous studies of ancient evolution or migration and current US census data; and are also readily reproducible and highly robust. Our approach presents a promising alternative for sample collection in studies focusing on the ethnic and genetic composition of populations or diseases associated with different mtDNA haplogroups and genotypes.

Complimentary Methods for Multivariate Genome-Wide Association Study Identify New Susceptibility Genes for Blood Cell Traits

Genome-wide association studies (GWAS) have found hundreds of novel loci associated with full blood count (FBC) phenotypes. However, most of these studies were performed in a single phenotype framework without putting into consideration the clinical relatedness among traits. In this work, in addition to the standard univariate GWAS, we also use two different multivariate methods to perform the first multiple traits GWAS of FBC traits in ∼7000 individuals from the Ugandan General Population Cohort (GPC). We started by performing the standard univariate GWAS approach. We then performed our first multivariate method, in this approach, we tested for marker associations with 15 FBC traits simultaneously in a multivariate mixed model implemented in GEMMA while accounting for the relatedness of individuals and pedigree structures, as well as population substructure. In this analysis, we provide a framework for the combination of multiple phenotypes in multivariate GWAS analysis and show evidence of multi-collinearity whenever the correlation between traits exceeds the correlation coefficient threshold of r ² >=0.75. This approach identifies two known and one novel loci. In the second multivariate method, we applied principal component analysis (PCA) to the same 15 correlated FBC traits. We then tested for marker associations with each PC in univariate linear mixed models implemented in GEMMA. We show that the FBC composite phenotype as assessed by each PC expresses information that is not completely encapsulated by the individual FBC traits, as this approach identifies three known and five novel loci that were not identified using both the standard univariate and multivariate GWAS methods. Across both multivariate methods, we identified six novel loci. As a proof of concept, both multivariate methods also identified known loci, HBB and ITFG3. The two multivariate methods show that multivariate genotype-phenotype methods increase power and identify novel genotype-phenotype associations not found with the standard univariate GWAS in the same dataset.

Discovery of mitochondrial DNA variants associated with genome-wide blood cell gene expression: a population-based mtDNA sequencing study

The effect of mitochondrial DNA (mtDNA) variation on peripheral blood transcriptomics in health and disease is not fully known. Sex-specific mitochondrially controlled gene expression patterns have been shown in Drosophila melanogaster but in humans, evidence is lacking. Functional variation in mtDNA may also have a role in the development of type 2 diabetes and its precursor state, i.e. prediabetes. We examined the associations between mitochondrial single-nucleotide polymorphisms (mtSNPs) and peripheral blood transcriptomics with a focus on sex- and prediabetes-specific effects. The genome-wide blood cell expression data of 19 637 probes, 199 deep-sequenced mtSNPs and nine haplogroups of 955 individuals from a population-based Young Finns Study cohort were used. Significant associations were identified with linear regression and analysis of covariance. The effects of sex and prediabetes on the associations between gene expression and mtSNPs were studied using random-effect meta-analysis. Our analysis identified 53 significant expression probe-mtSNP associations after Bonferroni correction, involving 7 genes and 31 mtSNPs. Eight probe-mtSNP signals remained independent after conditional analysis. In addition, five genes showed differential expression between haplogroups. The meta-analysis did not show any significant differences in linear model effect sizes between males and females but identified the association between the OASL gene and mtSNP C16294T to show prediabetes-specific effects. This study pinpoints new independent mtSNPs associated with peripheral blood transcriptomics and replicates six previously reported associations, providing further evidence of the mitochondrial genetic control of blood cell gene expression. In addition, we present evidence that prediabetes might lead to perturbations in mitochondrial control.

Comparing the Utility of Mitochondrial and Nuclear DNA to Adjust for Genetic Ancestry in Association Studies

Mitochondrial genome-wide association studies identify mitochondrial single nucleotide polymorphisms (mtSNPs) that associate with disease or disease-related phenotypes. Most mitochondrial and nuclear genome-wide association studies adjust for genetic ancestry by including principal components derived from nuclear DNA, but not from mitochondrial DNA, as covariates in statistical regression analyses. Furthermore, there is no standard when controlling for genetic ancestry during mitochondrial and nuclear genetic interaction association scans, especially across ethnicities with substantial mitochondrial genetic heterogeneity. The purpose of this study is to (1) compare the degree of ethnic variation captured by principal components calculated from microarray-defined nuclear and mitochondrial DNA and (2) assess the utility of mitochondrial principal components for association studies. Analytic techniques used in this study include a principal component analysis for genetic ancestry, decision-tree classification for self-reported ethnicity, and linear regression for association tests. Data from the Health and Retirement Study, which includes self-reported White, Black, and Hispanic Americans, was used for all analyses. We report that (1) mitochondrial principal component analysis (PCA) captures ethnic variation to a similar or slightly greater degree than nuclear PCA in Blacks and Hispanics, (2) nuclear and mitochondrial DNA classify self-reported ethnicity to a high degree but with a similar level of error, and 3) mitochondrial principal components can be used as covariates to adjust for population stratification in association studies with complex traits, as demonstrated by our analysis of height-a phenotype with a high heritability. Overall, genetic association studies might reveal true and robust mtSNP associations when including mitochondrial principal components as regression covariates.

Complex spatio-temporal distribution and genomic ancestry of mitochondrial DNA haplogroups in 24,216 Danes

Mitochondrial DNA (mtDNA) haplogroups (hgs) are evolutionarily conserved sets of mtDNA SNP-haplotypes with characteristic geographical distribution. Associations of hgs with disease and physiological characteristics have been reported, but have frequently not been reproducible. Using 418 mtDNA SNPs on the PsychChip (Illumina), we assessed the spatio-temporal distribution of mtDNA hgs in Denmark from DNA isolated from 24,642 geographically un-biased dried blood spots (DBS), collected from 1981 to 2005 through the Danish National Neonatal Screening program. ADMIXTURE was used to establish the genomic ancestry of all samples using a reference of 100K+ autosomal SNPs in 2,248 individuals from nine populations. Median-joining analysis determined that the hgs were highly variable, despite being typically Northern European in origin, suggesting multiple founder events. Furthermore, considerable heterogeneity and variation in nuclear genomic ancestry was observed. Thus, individuals with hg H exhibited 95%, and U hgs 38.2% - 92.5%, Danish ancestry. Significant clines between geographical regions and rural and metropolitan populations were found. Over 25 years, macro-hg L increased from 0.2% to 1.2% (p = 1.1*E-10), and M from 1% to 2.4% (p = 3.7*E-8). Hg U increased among the R macro-hg from 14.1% to 16.5% (p = 1.9*E-3). Genomic ancestry, geographical skewedness, and sub-hg distribution suggested that the L, M and U increases are due to immigration. The complex spatio-temporal dynamics and genomic ancestry of mtDNA in the Danish population reflect repeated migratory events and, in later years, net immigration. Such complexity may explain the often contradictory and population-specific reports of mito-genomic association with disease.

Schizophrenia-associated mt-DNA SNPs exhibit highly variable haplogroup affiliation and nuclear ancestry: Bi-genomic dependence raises major concerns for link to disease

Mitochondria play a significant role in human diseases. However, disease associations with mitochondrial DNA (mtDNA) SNPs have proven difficult to replicate. An analysis of eight schizophrenia-associated mtDNA SNPs, in 23,743 Danes without a psychiatric diagnosis and 2,538 schizophrenia patients, revealed marked inter-allelic differences in mitochondrial haplogroup affiliation and nuclear ancestry. This bi-genomic dependence could entail population stratification. Only two mitochondrial SNPs, m.15043A and m.15218G, were significantly associated with schizophrenia. However, these associations disappeared when corrected for haplogroup affiliation and nuclear ancestry. The extensive bi-genomic dependence documented here is a major concern when interpreting historic, as well as designing future, mtDNA association studies.

Meta-analysis identifies mitochondrial DNA sequence variants associated with walking speed

Declines in walking speed are associated with a variety of poor health outcomes including disability, comorbidity, and mortality. While genetic factors are putative contributors to variability in walking, few genetic loci have been identified for this trait. We examined the role of mitochondrial genomic variation on walking speed by sequencing the entire mitochondrial DNA (mtDNA). Data were meta-analyzed from 1758 Lifestyle Interventions and Independence for Elders (LIFE) Study and replication data from 730 Health, Aging, and Body Composition (HABC) Study participants with baseline walking speed information. Participants were 69+ years old of diverse racial backgrounds (African, European, and other race/ethnic groups) and had a wide range of mean walking speeds [4-6 m (0.78-1.09 m/s) and 400 m (0.83-1.24 m/s)]. Meta-analysis across studies and racial groups showed that m.12705C>T, ND5 variant was significantly associated (p < 0.0001) with walking speed at both short and long distances. Replication and meta-analysis also identified statistically significant walking speed associations (p < 0.0001) between the m.5460.G>A, ND2 and m.309C>CT, HV2 variants at short and long distances, respectively. All results remained statistically significant after multiple comparisons adjustment for 499 mtDNA variants. The m.12705C>T variant can be traced to the beginnings of human global migration and that cells carrying this variant display altered tRNA expression. Significant pooled effects related to stopping during the long-distance walk test were observed across OXPHOS complexes I (p = 0.0017) and III (p = 0.0048). These results suggest that mtDNA-encoded variants are associated with differences in walking speed among older adults, potentially identifying those at risk of developing mobility impairments.

Integrated Analysis of Genetic Ancestry and Genomic Alterations across Cancers

Disparities in cancer care have been a long-standing challenge. We estimated the genetic ancestry of The Cancer Genome Atlas patients, and performed a pan-cancer analysis on the influence of genetic ancestry on genomic alterations. Compared with European Americans, African Americans (AA) with breast, head and neck, and endometrial cancers exhibit a higher level of chromosomal instability, while a lower level of chromosomal instability was observed in AAs with kidney cancers. The frequencies of TP53 mutations and amplification of CCNE1 were increased in AAs in the cancer types showing higher levels of chromosomal instability. We observed lower frequencies of genomic alterations affecting genes in the PI3K pathway in AA patients across cancers. Our result provides insight into genomic contribution to cancer disparities.

Mitochondrial DNA sequence data reveals association of haplogroup U with psychosis in bipolar disorder

Converging genetic, postmortem gene-expression, cellular, and neuroimaging data implicate mitochondrial dysfunction in bipolar disorder. This study was conducted to investigate whether mitochondrial DNA (mtDNA) haplogroups and single nucleotide variants (SNVs) are associated with sub-phenotypes of bipolar disorder. MtDNA from 224 patients with Bipolar I disorder (BPI) was sequenced, and association of sequence variations with 3 sub-phenotypes (psychosis, rapid cycling, and adolescent illness onset) was evaluated. Gene-level tests were performed to evaluate overall burden of minor alleles for each phenotype. The haplogroup U was associated with a higher risk of psychosis. Secondary analyses of SNVs provided nominal evidence for association of psychosis with variants in the tRNA, ND4 and ND5 genes. The association of psychosis with ND4 (gene that encodes NADH dehydrogenase 4) was further supported by gene-level analysis. Preliminary analysis of mtDNA sequence data suggests a higher risk of psychosis with the U haplogroup and variation in the ND4 gene implicated in electron transport chain energy regulation. Further investigation of the functional consequences of this mtDNA variation is encouraged.

Mitonuclear linkage disequilibrium in human populations

There is extensive evidence from model systems that disrupting associations between co-adapted mitochondrial and nuclear genotypes can lead to deleterious and even lethal consequences. While it is tempting to extrapolate from these observations and make inferences about the human-health effects of altering mitonuclear associations, the importance of such associations may vary greatly among species, depending on population genetics, demographic history and other factors. Remarkably, despite the extensive study of human population genetics, the statistical associations between nuclear and mitochondrial alleles remain largely uninvestigated. We analysed published population genomic data to test for signatures of historical selection to maintain mitonuclear associations, particularly those involving nuclear genes that encode mitochondrial-localized proteins (N-mt genes). We found that significant mitonuclear linkage disequilibrium (LD) exists throughout the human genome, but these associations were generally weak, which is consistent with the paucity of population genetic structure in humans. Although mitonuclear LD varied among genomic regions (with especially high levels on the X chromosome), N-mt genes were statistically indistinguishable from background levels, suggesting that selection on mitonuclear epistasis has not preferentially maintained associations involving this set of loci at a species-wide level. We discuss these findings in the context of the ongoing debate over mitochondrial replacement therapy.

Ulcerative Colitis Is Under Dual (Mitochondrial and Nuclear) Genetic Control

BACKGROUND

Cellular oxidative stress and genetic susceptibility have been implicated in the multifactorial etiology of ulcerative colitis (UC). The nuclear genome association with UC has been intensely investigated, but the role of the mitochondrial DNA (mtDNA) has received far less attention and may account for part of the missing heritability. This study is a comprehensive analysis of the mtDNA contribution to UC susceptibility.

METHODS

The association of mitochondrial single-nucleotide polymorphisms (mtSNPs) and haplogroups with UC was tested in 488 cases and 833 controls of European ancestry from the NIDDK IBD Genetics Consortium Ulcerative Colitis Genome-Wide Association Study available through dbGaP and from the Illumina Genotype Control Database (studies 64 and 65).

RESULTS

No evidence of population stratification could be detected using 218 ancestry informative markers for European Americans. Seven of the 58 tested mtSNPs were nominally associated with UC, and A10550G in MT-ND4L withstands the Bonferroni correction (P = 1.29E-06, odds ratio [ORG] [95% confidence interval (CI)] = 4.80 [2.54-9.05], 10550G allele: 8.1% of patients and 1.9% of controls). A10550G remains equally associated after conditional analyses on the 11 UC genome-wide association studies (GWAS) top SNPs (6.35E-07 < Pcond < 4.58E-06), which suggests that it constitutes an independent risk factor from nuclear-encoded susceptibility loci. We detected additive (but not multiplicative) epistatic interactions between A10550G and all 11 top GWAS hits. Subhaplogroup K1 (P = 0.021, OR [95% CI] = 1.71 [1.08-2.69]) increased the risk for UC, whereas the U5b lineage conferred protection (P = 0.016, OR [95% CI] = 0.34 [0.14-0.82]).

CONCLUSIONS

These results suggest that UC has a dual mitochondrial and nuclear genetic control that warrants further replication in independent data sets and reinforces its etiopathogenic complexity.

Mitochondrial DNA sequence variation in multiple sclerosis

OBJECTIVE

To assess the influence of common mitochondrial DNA (mtDNA) sequence variation on multiple sclerosis (MS) risk in cases and controls part of an international consortium.

METHODS

We analyzed 115 high-quality mtDNA variants and common haplogroups from a previously published genome-wide association study among 7,391 cases from the International Multiple Sclerosis Genetics Consortium and 14,568 controls from the Wellcome Trust Case Control Consortium 2 project from 7 countries. Significant single nucleotide polymorphism and haplogroup associations were replicated in 3,720 cases and 879 controls from the University of California, San Francisco.

RESULTS

An elevated risk of MS was detected among haplogroup JT carriers from 7 pooled clinic sites (odds ratio [OR] = 1.15, 95% confidence interval [CI] = 1.07-1.24, p = 0.0002) included in the discovery study. The increased risk of MS was observed for both haplogroup T (OR = 1.17, 95% CI = 1.06-1.29, p = 0.002) and haplogroup J carriers (OR = 1.11, 95% CI = 1.01-1.22, p = 0.03). These haplogroup associations with MS were not replicated in the independent sample set. An elevated risk of primary progressive (PP) MS was detected for haplogroup J participants from 3 European discovery populations (OR = 1.49, 95% CI = 1.10-2.01, p = 0.009). This elevated risk was borderline significant in the US replication population (OR = 1.43, 95% CI = 0.99-2.08, p = 0.058) and remained significant in pooled analysis of discovery and replication studies (OR = 1.43, 95% CI = 1.14-1.81, p = 0.002). No common individual mtDNA variants were associated with MS risk.

CONCLUSIONS

Identification and validation of mitochondrial genetic variants associated with MS and PPMS may lead to new targets for treatment and diagnostic tests for identifying potential responders to interventions that target mitochondria.

Mitochondrial DNA Sequence Variation Associated With Peripheral Nerve Function in the Elderly

BACKGROUND

Mitochondrial dysfunction is a prominent hallmark of many sensory neuropathies. The purpose of this study was to assess the influence of mitochondrial DNA sequence variation on peripheral nerve function in the population-based Health, Aging, and Body Composition Study.

METHODS

We investigated the role of common mitochondrial DNA variation (n = 1,580) and complete mitochondrial DNA sequences (n = 138) on peroneal motor nerve conduction velocity and amplitude, average vibration detection threshold, and monofilament sensitivity.

RESULTS

Nominal associations among common mitochondrial DNA variants and haplogroups were identified but were not statistically significant after adjustment for multiple comparisons. Sequence-based approaches were used to identify aggregate variant associations across the 16S rRNA (weighted-sum, p = 2E-05 and variable threshold, p = 9E-06) for nerve conduction velocity. Several of these rare 16S variants occurred at or near sites with earlier disease associations and are also in close proximity to the peptidyl transferase center, which is the catalytic center of the 16S rRNA CONCLUSIONS: These results suggest that sequence variation related to mitochondrial protein synthesis/assembly is associated with peripheral nerve function and may provide insight into targets for intervention or new clinical strategies to preserve nerve function in late life.

The co-occurrence of mtDNA mutations on different oxidative phosphorylation subunits, not detected by haplogroup analysis, affects human longevity and is population specific

To re-examine the correlation between mtDNA variability and longevity, we examined mtDNAs from samples obtained from over 2200 ultranonagenarians (and an equal number of controls) collected within the framework of the GEHA EU project. The samples were categorized by high-resolution classification, while about 1300 mtDNA molecules (650 ultranonagenarians and an equal number of controls) were completely sequenced. Sequences, unlike standard haplogroup analysis, made possible to evaluate for the first time the cumulative effects of specific, concomitant mtDNA mutations, including those that per se have a low, or very low, impact. In particular, the analysis of the mutations occurring in different OXPHOS complex showed a complex scenario with a different mutation burden in 90+ subjects with respect to controls. These findings suggested that mutations in subunits of the OXPHOS complex I had a beneficial effect on longevity, while the simultaneous presence of mutations in complex I and III (which also occurs in J subhaplogroups involved in LHON) and in complex I and V seemed to be detrimental, likely explaining previous contradictory results. On the whole, our study, which goes beyond haplogroup analysis, suggests that mitochondrial DNA variation does affect human longevity, but its effect is heavily influenced by the interaction between mutations concomitantly occurring on different mtDNA genes.

Recent mitochondrial DNA mutations increase the risk of developing common late-onset human diseases

Mitochondrial DNA (mtDNA) is highly polymorphic at the population level, and specific mtDNA variants affect mitochondrial function. With emerging evidence that mitochondrial mechanisms are central to common human diseases, it is plausible that mtDNA variants contribute to the “missing heritability” of several complex traits. Given the central role of mtDNA genes in oxidative phosphorylation, the same genetic variants would be expected to alter the risk of developing several different disorders, but this has not been shown to date. Here we studied 38,638 individuals with 11 major diseases, and 17,483 healthy controls. Imputing missing variants from 7,729 complete mitochondrial genomes, we captured 40.41% of European mtDNA variation. We show that mtDNA variants modifying the risk of developing one disease also modify the risk of developing other diseases, thus providing independent replication of a disease association in different case and control cohorts. High-risk alleles were more common than protective alleles, indicating that mtDNA is not at equilibrium in the human population, and that recent mutations interact with nuclear loci to modify the risk of developing multiple common diseases.

There is a growing body of evidence indicating that mitochondrial dysfunction, a result of genetic variation in the mitochondrial genome, is a critical component in the aetiology of a number of complex traits. Here, we take advantage of recent technical and methodological advances to examine the role of common mitochondrial DNA variants in several complex diseases. By examining over 50,000 individuals, from 11 different diseases we show that mitochondrial DNA variants can both increase or decrease an individual's risk of disease, replicating and expanding upon several previously reported studies. Moreover, by analysing several large disease groups in tandem, we are able to show a commonality of association, with the same mitochondrial DNA variants associated with several distinct disease phenotypes. These shared genetic associations implicate a shared underlying functional effect, likely changing cellular energy, which manifests as distinct phenotypes. Our study confirms the important role that mitochondrial DNA variation plays on complex traits and additionally supports the utility of a GWAS-based approach for analysing mitochondrial genetics.

Genetic variation of oxidative phosphorylation genes in stroke and Alzheimer's disease

Previous research implicates alterations in oxidative phosphorylation (OXPHOS) in the development of Alzheimer's disease (AD). We sought to test whether genetic variants within OXPHOS genes increase the risk of AD. We first used gene-set enrichment analysis to identify associations, and then applied a previously replicated stroke genetic risk score to determine if OXPHOS genetic overlap exists between stroke and AD. Gene-set enrichment analysis identified associations between variation in OXPHOS genes and AD versus control status (p = 0.012). Conversion from cognitively normal controls to mild cognitive impairment was also associated with the OXPHOS gene-set (p = 0.045). Subset analyses demonstrated association for complex I genes (p < 0.05), but not for complexes II-V. Among neuroimaging measures, hippocampal volume and entorhinal cortex thickness were associated with OXPHOS genes (all p < 0.025). The stroke genetic risk score demonstrated association with clinical status, baseline and longitudinal imaging measures (p < 0.05). OXPHOS genetic variation influences clinical status and neuroimaging intermediates of AD. OXPHOS genetic variants associated with stroke are also linked to AD progression. Further studies are needed to explore functional consequences of these OXPHOS variants.

Mitochondrial DNA sequence associations with dementia and amyloid-β in elderly African Americans

Mitochondrial dysfunction occurs early in the course of several neurodegenerative diseases, and is potentially related to increased oxidative damage and amyloid-β (Aβ) formation in Alzheimer's disease. The goals of this study were to assess mtDNA sequence associations with dementia risk, 10-year cognitive change, and markers of oxidative stress and Aβ among 1089 African-Americans in the population-based Health, Aging, and Body Composition Study. Participants were free of dementia at baseline, and incidence was determined in 187 (18%) cases over 10 to 12 follow-up years. Haplogroup L1 participants were at increased risk for developing dementia (odds ratio = 1.88, 95% confidence interval = 1.23-2.88, p = 0.004), lower plasma Aβ42 levels (p = 0.03), and greater 10-year decline on the Digit Symbol Substitution Test (p = 0.04) when compared with common haplogroup L3. The p.V193I, ND2 substitution was associated with significantly higher Aβ42 levels (p = 0.0012), and this association was present in haplogroup L3 (p = 0.018) but not L1 (p = 0.90) participants. All associations were independent of potential confounders, including APOEε4 status and nuclear genetic ancestry. Identification of mtDNA sequence variation associated with dementia risk and cognitive decline may contribute to the development of new treatment targets and diagnostic tests that identify responders to interventions targeting mitochondria.

Pathway-based approaches for sequencing-based genome-wide association studies

For analyzing complex trait association with sequencing data, most current studies test aggregated effects of variants in a gene or genomic region. Although gene-based tests have insufficient power even for moderately sized samples, pathway-based analyses combine information across multiple genes in biological pathways and may offer additional insight. However, most existing pathway association methods are originally designed for genome-wide association studies, and are not comprehensively evaluated for sequencing data. Moreover, region-based rare variant association methods, although potentially applicable to pathway-based analysis by extending their region definition to gene sets, have never been rigorously tested. In the context of exome-based studies, we use simulated and real datasets to evaluate pathway-based association tests. Our simulation strategy adopts a genome-wide genetic model that distributes total genetic effects hierarchically into pathways, genes, and individual variants, allowing the evaluation of pathway-based methods with realistic quantifiable assumptions on the underlying genetic architectures. The results show that, although no single pathway-based association method offers superior performance in all simulated scenarios, a modification of Gene Set Enrichment Analysis approach using statistics from single-marker tests without gene-level collapsing (weighted Kolmogrov-Smirnov [WKS]-Variant method) is consistently powerful. Interestingly, directly applying rare variant association tests (e.g., sequence kernel association test) to pathway analysis offers a similar power, but its results are sensitive to assumptions of genetic architecture. We applied pathway association analysis to an exome-sequencing data of the chronic obstructive pulmonary disease, and found that the WKS-Variant method confirms associated genes previously published.

Common variants within oxidative phosphorylation genes influence risk of ischemic stroke and intracerebral hemorrhage

BACKGROUND AND PURPOSE

Previous studies demonstrated association between mitochondrial DNA variants and ischemic stroke (IS). We investigated whether variants within a larger set of oxidative phosphorylation (OXPHOS) genes encoded by both autosomal and mitochondrial DNA were associated with risk of IS and, based on our results, extended our investigation to intracerebral hemorrhage (ICH).

METHODS

This association study used a discovery cohort of 1643 individuals, a validation cohort of 2432 individuals for IS, and an extension cohort of 1476 individuals for ICH. Gene-set enrichment analysis was performed on all structural OXPHOS genes, as well as genes contributing to individual respiratory complexes. Gene-sets passing gene-set enrichment analysis were tested by constructing genetic scores using common variants residing within each gene. Associations between each variant and IS that emerged in the discovery cohort were examined in validation and extension cohorts.

RESULTS

IS was associated with genetic risk scores in OXPHOS as a whole (odds ratio [OR], 1.17; P=0.008) and complex I (OR, 1.06; P=0.050). Among IS subtypes, small vessel stroke showed association with OXPHOS (OR, 1.16; P=0.007), complex I (OR, 1.13; P=0.027), and complex IV (OR, 1.14; P=0.018). To further explore this small vessel association, we extended our analysis to ICH, revealing association between deep hemispheric ICH and complex IV (OR, 1.08; P=0.008).

CONCLUSIONS

This pathway analysis demonstrates association between common genetic variants within OXPHOS genes and stroke. The associations for small vessel stroke and deep ICH suggest that genetic variation in OXPHOS influences small vessel pathobiology. Further studies are needed to identify culprit genetic variants and assess their functional consequences.

The utility of mitochondrial and y chromosome phylogenetic data to improve correction for population stratification

Genome-wide association (GWA) studies have become a standard approach for discovering and validating genomic polymorphisms putatively associated with phenotypes of interest. Accounting for population structure in GWA studies is critical to attain unbiased parameter measurements and control Type I error. One common approach to accounting for population structure is to include several principal components derived from the entire autosomal dataset, which reflects population structure signal. However, knowing which components to include is subjective and generally not conclusive. We examined how phylogenetic signal from mitochondrial DNA (mtDNA) and chromosome Y (chr:Y) markers is concordant with principal component data based on autosomal markers to determine whether mtDNA and chr:Y phylogenetic data can help guide principal component selection. Using HAPMAP and other original data from individuals of multiple ancestries, we examined the relationships of mtDNA and chr:Y phylogenetic signal with the autosomal PCA using best subset logistic regression. We show that while the two approaches agree at times, this is independent of the component order and not completely represented in the Eigen values. Additionally, we use simulations to demonstrate that our approach leads to a slightly reduced Type I error rate compared to the standard approach. This approach provides preliminary evidence to support the theoretical concept that mtDNA and chr:Y data can be informative in locating the PCs that are most associated with evolutionary history of populations that are being studied, although the utility of such information will depend on the specific situation.

Population genetic structure and origins of Native Hawaiians in the multiethnic cohort study

The population genetic structure of Native Hawaiians has yet to be comprehensively studied, and the ancestral origins of Polynesians remain in question. In this study, we utilized high-resolution genome-wide SNP data and mitochondrial genomes of 148 and 160 Native Hawaiians, respectively, to characterize their population structure of the nuclear and mitochondrial genomes, ancestral origins, and population expansion. Native Hawaiians, who self-reported full Native Hawaiian heritage, demonstrated 78% Native Hawaiian, 11.5% European, and 7.8% Asian ancestry with 99% belonging to the B4 mitochondrial haplogroup. The estimated proportions of Native Hawaiian ancestry for those who reported mixed ancestry (i.e. 75% and 50% Native Hawaiian heritage) were found to be consistent with their self-reported heritage. A significant proportion of Melanesian ancestry (mean = 32%) was estimated in 100% self-reported Native Hawaiians in an ADMIXTURE analysis of Asian, Melanesian, and Native Hawaiian populations of K = 2, where K denotes the number of ancestral populations. This notable proportion of Melanesian admixture supports the “Slow-Boat” model of migration of ancestral Polynesian populations from East Asia to the Pacific Islands. In addition, approximately 1,300 years ago a single, strong expansion of the Native Hawaiian population was estimated. By providing important insight into the underlying population structure of Native Hawaiians, this study lays the foundation for future genetic association studies of this U.S. minority population.

Evaluating mitochondrial DNA variation in autism spectrum disorders

Despite the increasing speculation that oxidative stress and abnormal energy metabolism may play a role in Autism Spectrum Disorders (ASD), and the observation that patients with mitochondrial defects have symptoms consistent with ASD, there are no comprehensive published studies examining the role of mitochondrial variation in autism. Therefore, we have sought to comprehensively examine the role of mitochondrial DNA (mtDNA) variation with regard to ASD risk, employing a multi-phase approach. In phase 1 of our experiment, we examined 132 mtDNA single-nucleotide polymorphisms (SNPs) genotyped as part of our genome-wide association studies of ASD. In phase 2 we genotyped the major European mitochondrial haplogroup-defining variants within an expanded set of autism probands and controls. Finally in phase 3, we resequenced the entire mtDNA in a subset of our Caucasian samples (∼400 proband-father pairs). In each phase we tested whether mitochondrial variation showed evidence of association to ASD. Despite a thorough interrogation of mtDNA variation, we found no evidence to suggest a major role for mtDNA variation in ASD susceptibility. Accordingly, while there may be attractive biological hints suggesting the role of mitochondria in ASD our data indicate that mtDNA variation is not a major contributing factor to the development of ASD.

Mitochondrial DNA sequence variation is associated with free-living activity energy expenditure in the elderly

The decline in activity energy expenditure underlies a range of age-associated pathological conditions, neuromuscular and neurological impairments, disability, and mortality. The majority (90%) of the energy needs of the human body are met by mitochondrial oxidative phosphorylation (OXPHOS). OXPHOS is dependent on the coordinated expression and interaction of genes encoded in the nuclear and mitochondrial genomes. We examined the role of mitochondrial genomic variation in free-living activity energy expenditure (AEE) and physical activity levels (PAL) by sequencing the entire (~16.5 kilobases) mtDNA from 138 Health, Aging, and Body Composition Study participants. Among the common mtDNA variants, the hypervariable region 2 m.185G>A variant was significantly associated with AEE (p=0.001) and PAL (p=0.0005) after adjustment for multiple comparisons. Several unique nonsynonymous variants were identified in the extremes of AEE with some occurring at highly conserved sites predicted to affect protein structure and function. Of interest is the p.T194M, CytB substitution in the lower extreme of AEE occurring at a residue in the Qi site of complex III. Among participants with low activity levels, the burden of singleton variants was 30% higher across the entire mtDNA and OXPHOS complex I when compared to those having moderate to high activity levels. A significant pooled variant association across the hypervariable 2 region was observed for AEE and PAL. These results suggest that mtDNA variation is associated with free-living AEE in older persons and may generate new hypotheses by which specific mtDNA complexes, genes, and variants may contribute to the maintenance of activity levels in late life.

Mitochondrial DNA sequence variation associated with dementia and cognitive function in the elderly

Mitochondrial dysfunction is a prominent hallmark of Alzheimer's disease (AD). Mitochondrial DNA (mtDNA) damage may be a major cause of abnormal reactive oxidative species production in AD or increased neuronal susceptibility to oxidative injury during aging. The purpose of this study was to assess the influence of mtDNA sequence variation on clinically significant cognitive impairment and dementia risk in the population-based Health, Aging, and Body Composition (Health ABC) Study. We first investigated the role of common mtDNA haplogroups and individual variants on dementia risk and 8-year change on the Modified Mini-Mental State Examination (3MS) and Digit Symbol Substitution Test (DSST) among 1,631 participants of European genetic ancestry. Participants were free of dementia at baseline and incidence was determined in 273 cases from hospital and medication records over 10-12 follow-up years. Participants from haplogroup T had a statistically significant increased risk of developing dementia (OR = 1.86, 95% CI = 1.23, 2.82, p = 0.0008) and haplogroup J participants experienced a statistically significant 8-year decline in 3MS (β = -0.14, 95% CI = -0.27, -0.03, p = 0.0006), both compared with common haplogroup H. The m.15244A>G, p.G166G, CytB variant was associated with a significant decline in DSST score (β = -0.58, 95% CI -0.89, -0.28, p = 0.00019) and the m.14178T>C, p.I166V, ND6 variant was associated with a significant decline in 3MS score (β = -0.87, 95% CI -1.31, -3.86, p = 0.00012). Finally, we sequenced the complete ~16.5 kb mtDNA from 135 Health ABC participants and identified several highly conserved and potentially functional nonsynonymous variants unique to 22 dementia cases and aggregate sequence variation across the hypervariable 2-3 regions that influences 3MS and DSST scores.

Mitochondrial DNA sequence variation and risk of pancreatic cancer

Although the mitochondrial genome exhibits high mutation rates, common mitochondrial DNA (mtDNA) variation has not been consistently associated with pancreatic cancer. Here, we comprehensively examined mitochondrial genomic variation by sequencing the mtDNA of participants (cases = 286, controls = 283) in a San Francisco Bay Area pancreatic cancer case-control study. Five common variants were associated with pancreatic cancer at nominal statistical significance (P < 0.05) with the strongest finding for mt5460g in the ND2 gene [OR = 3.9; 95% confidence interval (CI), 1.5-10; P = 0.004] which encodes an A331T substitution. Haplogroup K was nominally associated with reduced pancreatic cancer risk (OR = 0.32; 95% CI, 0.13-0.76; P = 0.01) when compared with the most common haplogroup, H. A total of 19 haplogroup-specific rare variants yielded nominal statistically significant associations (P < 0.05) with pancreatic cancer risk, with the majority observed in genes involved in oxidative phosphorylation. Weighted-sum statistics were used to identify an aggregate effect of variants in the 22 mitochondrial tRNAs on pancreatic cancer risk (P = 0.02). While the burden of singleton variants in the HV2 and 12S RNA regions was three times higher among European haplogroup N cases than controls, the prevalence of singleton variants in ND4 and ND5 was two to three times higher among African haplogroup L cases than in controls. Together, the results of this study provide evidence that aggregated common and rare variants and the accumulation of singleton variants are important contributors to pancreatic cancer risk.

Mitochondrial DNA variation in human metabolic rate and energy expenditure

The role of climate in driving selection of mtDNA as Homo sapiens migrated out of Africa into Eurasia remains controversial. We evaluated the role of mtDNA variation in resting metabolic rate (RMR) and total energy expenditure (TEE) among 294 older, community-dwelling African and European American adults from the Health, Aging and Body Composition Study. Common African haplogroups L0, L2 and L3 had significantly lower RMRs than European haplogroups H, JT and UK with haplogroup L1 RMR being intermediate to these groups. This study links mitochondrial haplogroups with ancestry-associated differences in metabolic rate and energy expenditure.

HaploGrep: a fast and reliable algorithm for automatic classification of mitochondrial DNA haplogroups

An ongoing source of controversy in mitochondrial DNA (mtDNA) research is based on the detection of numerous errors in mtDNA profiles that led to erroneous conclusions and false disease associations. Most of these controversies could be avoided if the samples' haplogroup status would be taken into consideration. Knowing the mtDNA haplogroup affiliation is a critical prerequisite for studying mechanisms of human evolution and discovering genes involved in complex diseases, and validating phylogenetic consistency using haplogroup classification is an important step in quality control. However, despite the availability of Phylotree, a regularly updated classification tree of global mtDNA variation, the process of haplogroup classification is still time-consuming and error-prone, as researchers have to manually compare the polymorphisms found in a population sample to those summarized in Phylotree, polymorphism by polymorphism, sample by sample. We present HaploGrep, a fast, reliable and straight-forward algorithm implemented in a Web application to determine the haplogroup affiliation of thousands of mtDNA profiles genotyped for the entire mtDNA or any part of it. HaploGrep uses the latest version of Phylotree and offers an all-in-one solution for quality assessment of mtDNA profiles in clinical genetics, population genetics and forensics. HaploGrep can be accessed freely at http://haplogrep.uibk.ac.at.

Common mitochondrial sequence variants in ischemic stroke

OBJECTIVE

Rare mitochondrial mutations cause neurologic disease, including ischemic stroke and MRI white matter changes. We investigated whether common mitochondrial genetic variants influence risk of sporadic ischemic stroke and, in patients with stroke, the volume of white matter hyperintensity (WMHV).

METHODS

In this multicenter, mitochondrial genome-wide association study (GWAS), 2284 ischemic stroke cases and 1728 controls from the International Stroke Genetics Consortium were genotyped for 64 mitochondrial single nucleotide polymorphisms (SNPs). Imputation resulted in 144 SNPs, which were tested in each cohort and in meta-analysis for ischemic stroke association. A genetic score of all mitochondrial variants was also tested in association with ischemic stroke.

RESULTS

No individual SNP reached adjusted significance in meta-analysis. A genetic score comprised of the summation of contributions from individual variants across the mitochondrial genome showed association with ischemic stroke in meta-analysis (odds ratio [OR] = 1.13, p < 0.0001) with minimal heterogeneity (I(2) = 0.00). This ischemic stroke score was robust to permutation, and was also associated with WMHV in 792 nested case individuals with ischemic stroke (p = 0.037).

INTERPRETATION

In this mitochondrial GWAS of ischemic stroke, a genetic score comprised of the sum of all common variants in the mitochondrial genome showed association with ischemic stroke. In an independent analysis of a related trait, this same score correlated with WMHV in stroke cases. Despite this aggregate association, no individual variant reached significance. Substantially larger studies will be required to identify precise sequence variants influencing cerebrovascular disease.