Examining the Association between Mitochondrial Genome Variation and Coronary Artery Disease
============================================================================================

* Baiba Vilne
* Aniket Sawant
* Irina Rudaka

## Abstract

**Background** Large-scale genome-wide association studies have identified hundreds of single-nucleotide variants (SNVs) significantly associated with coronary artery disease (CAD). However, collectively, these explain <20% of the heritability.

**Hypothesis** Here, we hypothesize that mitochondrial (MT) SNVs might present one potential source of this “missing heritability”.

**Methods** We analyzed 265 MT-SNVs in ∼500,000 UK Biobank individuals, exploring two different CAD definitions: a more stringent (myocardial infarction and/or revascularization; HARD=20,405), and a more inclusive (also angina and chronic ischemic heart disease; SOFT=34,782).

**Results** In HARD cases, the most significant (P<0.05) associations were for m.295C>T (control region) and m.12612A>G (ND5), found more frequently in cases (OR=1.05), potentially related to reduced cardiorespiratory fitness in response to exercise, as well as for m.12372G>A (ND5) and m.11467A>G (ND4), present more frequently in controls (OR=0.97), previously associated with lower ROS production rate. In SOFT cases, four MT-SNVs survived multiple testing correction (at FDR<5%), all potentially conferring increased CAD risk. Of those, m.11251A>G (ND4) and m.15452C>A (CYB) have previously shown significant associations with body height. In line with this, we observed that CAD cases were slightly less physically active and their average body height was ∼2.00 cm lower compared to controls, both traits known to be related to an increased CAD risk. Gene-based tests identified CO2 associated with HARD/SOFT CAD, whereas ND3 and CYB associated with SOFT cases (P<0.05), dysfunction of which has been related to MT oxidative stress, obesity/T2D (CO2), BMI (ND3) and angina/exercise intolerance (CYB). Finally, we observed that macro-haplogroup I was significantly (P<0.05) more frequent in HARD cases vs. controls (3.35% vs. 3.08%), potentially associated with response to exercise.

**Conclusions** We found only spurious associations between MT genome variation and HARD/SOFT CAD and conclude that more MT-SNV data in even larger study cohorts may be needed to conclusively determine the role of MT-DNA in CAD.

Keywords
*   Coronary Artery Disease
*   Mitochondria
*   Mitochondrial DNA Variants
*   Haplogroups
*   Association
*   Common
*   Rare Variants

## 1. Introduction

Coronary artery disease (CAD) and its major complication myocardial infarction is the most common cardiovascular disease and the main leading cause of morbidity and mortality worldwide. CAD is posing a huge socio-economic burden to the society and health systems [65] and its prevalence is expected to increase in the coming years [29, 71, 87]. CAD is a multi-factorial disease with complex etiology, considered to be driven by both environment/lifestyle and genetic factors [22, 28, 101]. Over the last 14 years, several large-scale genome-wide association studies and their meta-analysis have identified numerous common genetic variants associated with CAD risk [20, 27, 44, 77, 79, 85, 88, 96, 97, 106] and explored their functional consequences [1, 12, 49–51, 63, 76, 89, 102, 109]. However, collectively, these variants explain only a small proportion (∼20%) of the disease heritability [20, 56]. Genetic variations of the mitochondrial (MT) DNA have remained out of focus for a long time and present an underexplored potential source of the “missing heritability” of several complex traits, including CAD [11, 45, 91].

The human MT-DNA is a maternally-inherited, double-stranded, circular, histone-free “chromosome” of 16,596 base pairs (bp). Each mitochondrion contains 2 to 10 copies of MT-DNA and, depending on the tissue energy requirement, each human cell may contain hundreds of mitochondria [4]. MT-DNA encodes 37 genes corresponding to subunits ND1 to 6 (and 4L) of respiratory complex I, catalytic subunits I-III (CO1-3) of cytochrome c oxidase (respiratory complex IV), subunits adenosine triphosphate 6 and 8 (ATP6 and 8) of F1F0 ATPase and cytochrome b of respiratory complex III. The remaining genes encode 2 ribosomal RNAs (16S and 12S rRNAs) and 22 transfer RNAs (tRNAs), used for mitochondrial protein synthesis [11, 92, 104]. All of them are involved in oxidative phosphorylation (OXPHOS), the process by which ATP, the major source of energy, is being synthesized [25, 37, 74].

A toxic by-product of OXPHOS is the production of reactive oxygen species (ROS), unstable compounds which can generate free radicals [58]. Mitochondria are the primary source of endogenous ROS [58]. By antioxidant defense, cells can manage certain levels of free radical production. However, if threshold levels are exceeded, a state of oxidative stress occurs [68], which is known to play a vital role in the pathogenesis of atherosclerosis and CAD [69, 102]. Many of the common CAD risk factors such as age, hypertension, hyperglycemia, high cholesterol levels and smoking are also known to perturb mitochondrial function and increase oxidative stress [11]. In particular, the lifestyle of modern societies is increasingly related to reduced physical activity and calorically excessive diets [21, 70], favoring ROS production and possibly increasing the risk of CAD [61, 82], due to mitochondrial dysfunction, reduced bioenergetic capacity and disrupted redox homeostasis [6, 67].

Although the role for mitochondrial dysfunction in CAD etiology is well-established, the role of the mitochondrial genome in this process has not been extensively investigated [11]. Several forms of cardiovascular disease have been related to the presence of pathogenic mitochondrial genome mutations. However, the vast majority of mitochondrial genetic variation, “natural” single nucleotide variants (SNVs) have not been directly linked to disease pathogenesis [11]. During evolution, a number of such MT-SNVs have accumulated in mitochondrial genomes subdividing the human population into several discrete (geographic region specific) mitochondrial phylogenetic clades or haplogroups [95]. As the mitochondrial genome does not undergo DNA recombination, haplogroups are relatively stable and enable the clustering of individuals based on their shared maternal ancestry [95]. These clusters are often associated with different racial/ethnic groups [11]. Considering that family history and race/ethnicity is known to influence CAD risk, it is reasonable to assume that mitochondrial haplogroups may contribute to this heritable modulator of CAD susceptibility [11].

Increasing evidence suggests that non-pathogenic “natural” MT-SNVs and haplogroups may be associated with alterations in mitochondrial function [11, 25, 103], including subtle differences in OXPHOS capacity and the generation of ROS [30, 84], thus, potentially influencing the onset of and differences in CAD risk of different human populations [11, 80]. For example, MT-SNVs m.8701A>G and m.10398A>G have been demonstrated to lead to lower matrix pH and alterations in intracellular calcium levels [48]. In addition, several other MT-SNVs have been previously identified as significantly associated with HDL cholesterol and triglycerides levels [34], as well as with several metabolic traits including insulin levels and body weight [57]. Mitochondria bearing macro-haplogroup H-defining SNVs have been observed to possess several differences in basal oxygen consumption, maximal respiratory capacity and oxidant production relative to other European haplogroups, especially macro-haplogroup J [72] and African macro-haplogroup L [59]. Similarly, macro-haplogroup T has been associated with CAD and diabetic retinopathy in Middle European Caucasian populations [55]. On the other hand, several haplogroups in both Europe and Asia have been related to increased longevity, likely via reduced oxidant generation, which is a driver of the ageing process [16, 23].

In this study, we hypothesize that mitochondrial genome variation might present a potential source of the “missing heritability” of CAD. To explore this hypothesis, we performed: (1) association analyses of common/low-frequency MT-SNVs (MAF>0.01; N=111) with CAD; (2) gene-based tests to investigate the cumulative impact of all MT-SNVs on the mitochondrial genes in relation with CAD and (3) comparisons of mitochondrial haplogroup frequencies of individuals with CAD. In all cases, we explored two different CAD definitions (as previously used by Nelson et al. [77]): a more stringent (HARD=20,405), considering only myocardial infarction and/or revascularization and a more inclusive (SOFT=34,782), including all HARD CAD cases, as well as also angina and chronic ischemic heart disease vs. controls in a cohort of ca. 500,000 UK Biobank individuals. The complete workflow of this analysis is summarized in Fig.1.

![Figure 1](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/02/10/2022.02.09.22270723/F1.medium.gif)

[Figure 1](http://medrxiv.org/content/early/2022/02/10/2022.02.09.22270723/F1)

Figure 1 
A complete workflow of the analyses performed.

## 2. Materials and Methods

### 2.1. Study population, disease phenotypes and quality filtering

The UK Biobank [13] is a large population-based prospective cohort study from the United Kingdom with genetic and deep phenotypic (∼7,221 phenotypes [http://www.nealelab.is/uk-biobank](http://www.nealelab.is/uk-biobank)) data on ca. 500,000 individuals aged 40 to 69. We downloaded these data (application ID 61684) and used similar CAD case definition, as previously described by Nelson et al. [77] for UK Biobank. HARD CAD cases included individuals with fatal or non-fatal myocardial infarction (MI), percutaneous transluminal coronary angioplasty (PTCA) or coronary artery bypass grafting (CABG). SOFT CAD included individuals meeting the HARD CAD definition as well as those with chronic ischemic heart disease (IHD) and angina (Fig.1). In HESIN hospital episodes’ data and death registry data from diagnosis and operation (primary and secondary causes), MI was defined as hospital admission or cause of death due to ICD9 410-412, ICD10 I21-I24, I25.2; PTCA was defined as hospital admission for PTCA (OPCS-4 K49, K50.1, K75); CABG was defined as hospital admission for CABG (OPCS-4 K40-K46); and angina or chronic IHD was defined as hospital admission or death due to ICD9 413, 414.0, 414.8, 414.9, ICD10 I20, I25.1, I25.5-I25.9. In UK Biobank self-reported data, cases were defined as having ‘vascular/heart problems diagnosed by doctor’ or ‘non-cancer illnesses that self-reported as angina or heart attack’. Self-reported surgery included PTCA, CABG or triple heart bypass. All participants not defined as CAD cases using the SOFT definition were considered as controls in the analysis. For a complete list of definition codes, see Supplementary Table 1. We subsequently performed individual-level filtering (Fig.1) by removing missingness or heterozygosity outliers, participants with self-reported vs. genetically inferred sex mismatches or putative sex chromosome aneuploidy, individuals that were not of European (EUR) ancestry and individuals having withdrawn their consent at the time of analysis. We also identified closely related participants (kinship coefficient >0.088 i.e., first or second-degree relative pairs), preferentially retaining CAD cases or relative with the highest call-rate.

The following individual characteristics were also extracted in order to characterize HARD/SOFT CAD cases vs. controls: age at recruitment (field #21022), sex (field #31), BMI (field #21001), height (field #50), hypertension (fields #4080 and #4079), hypercholesterolemia (self-reported data and ICD9/10) and (self-reported) use of cholesterol lowering drugs, insulin and blood pressure medications (field #6153), type 2 diabetes (T2D, fields #41201, #41202 and #4120, E11), glycemic control, obesity, smoking status (‘ever smoked’: field #20160 and ‘current’ from ‘smoking status’ field #20116), family history of heart disease (fields #20107, #20110, #20111 in any first-degree relative, i.e. father, mother or sibling, respectively). Data at the time of first assessment were obtained and processed to binary (yes / no) values or mean values for fields with continuous data with multiple readings at the time of first assessment. The statistical tests used were Mann-Whitney test for comparison of the continuous measures between HARD vs. SOFT cases vs. controls and the Chi-square test for goodness of fit, when comparing nominal/binary data in the same groups.

### 2.2. Genotype data quality control

In UK Biobank [13], genotyping was performed using Affymetrix UK biobank Axiom (450,000 samples) and Affymetrix UK BiLEVE Axiom (50,000 samples) arrays (Fig.1) and the autosomal genetic data were then imputed to the Haplotype Reference Consortium panel and UK10K4 + 1000 Genomes panel. We downloaded the genotype data for the 265 MT DNA variants for all 500,000 individuals and pre-processed MT DNA data as previously described in ref. [57]. In brief, we first made sure that the reference alleles match the latest MT Cambridge Revised Sequence (rCRS) of the Human MT DNA positions. After setting all potential heterozygotes to missing, further quality control of genotyped individuals included filtering for missingness by individual <0.1 and missingness by SNV<0.1 with PLINK [14]. For common/low-frequency variant association analyses, we also required that the minor allele frequency (MAF)>0.01. An overview of the filtering of MT-SNVs is provided in Fig.1

### 2.3. MT-SNV association analyses

For common and low-frequency (MAF >0.01; n=111) variants, we performed single marker tests to explore their associations with HARD and SOFT CAD (Fig.1) using SNPTEST v2.5.4 with the frequentist test and expected method, as previously described by ref. [9]. We used as covariates the array (UK Biobank vs UK BiLEVE), sex, birth year and the first five principal components of the autosomal genotype data, provided by the UK Biobank, similar to ref. [77] and Benjamini-Hochberg (BH) [8] adjustment for multiple testing was applied to calculate the false discovery rate (FDR). MT-SNV annotations were performed using the manually-curated database, HmtVar ([https://www.hmtvar.uniba.it/](https://www.hmtvar.uniba.it/)).

### 2.4. MT-gene-based association analyses

To also consider the potential effects of rare (MAF ≤ 0.01) variants on CAD risk, we assigned all SNVs to MT genes based on MITOMAP [https://www.mitomap.org/MITOMAP](https://www.mitomap.org/MITOMAP) and used the R software package SKAT (v2.0.1) [62] to perform MT-gene-based (additionally including the whole mitochondrion as our region of interest, MT) association analyses with HARD and SOFT CAD phenotypes (Fig.1), again using as covariates the array (UK Biobank vs UK BiLEVE), sex, birth year and the first five principal components of the autosomal genotype data, provided by the UK Biobank, similar to ref. [77] and obtain resampled residuals (n.Resampling=1000, type.Resampling=“bootstrap”) to compute resampling P value.

### 2.5. Haplogroup assignment

We used the PhyloTree Build 17 [98] as implemented in HaploGrep (v2.2.8) [53] to estimate mitochondrial haplogroups in our dataset. Thereafter, we assigned individuals to one of the major European haplogroups (H, I, J, K, R, T, U, V, W, X), or to a group of ‘others’ [https://www.mitomap.org/foswiki/pub/MITOMAP/WebHome/simple-tree-mitomap-2019.pdf](https://www.mitomap.org/foswiki/pub/MITOMAP/WebHome/simple-tree-mitomap-2019.pdf). Fisher’s exact test [31] was used to calculate the statistical significance of the overlaps between haplogroups and HARD and SOFT CAD phenotypes. Benjamini-Hochberg (BH) [8] adjustment for multiple testing was applied to calculate the false discovery rate (FDR) (Fig.1).

## 3. Results

### 3.1. Characteristics of study subjects

The current study included ca. 500,000 genotyped individuals from the UK Biobank [13] 48,700 with an inclusive CAD phenotype (SOFT) that incorporated self-reported angina or other evidence of chronic coronary heart disease, of which 28,503 had a more stringently defined CAD phenotype (HARD) of myocardial infarction (Fig.1), similarly as in ref. [77]. All participants (N=453,805) not defined as CAD cases using the SOFT definition will be considered as controls in the analysis. After this step of quality control 45,285 individuals were removed, 24,770 cases with the HARD, 42,079 cases with the SOFT CAD phenotype and 415,271 controls remained. Finally, further quality control of genotyped individuals included filtering for missingness by individual <0.1 and missingness by SNV<0.1 was performed, considering the 265 mitochondrial variants present on the UK Biobank or UK BiLEVE arrays (as described in section 2.2. in Methods). As a result, a further set of 54,474 individuals were removed, leaving us with 20,405 cases with the HARD and 34,782 cases with the SOFT phenotype vs. 356,563 controls. Individual characteristics, in terms of common CAD risk factors, of these individuals are summarized in Table 1.

View this table:
[Table 1](http://medrxiv.org/content/early/2022/02/10/2022.02.09.22270723/T1)

Table 1 
Individual characteristics of HARD and SOFT CAD cases vs. controls. \***| Represents statistically significant (of P<0.001) difference between HARD and SOFT CAD cases vs. controls, whereas +++ represents statistically significant (of P<0.001) difference between HARD vs. SOFT CAD cases.

On average, 77/67% vs. 43% of the individuals were men in HARD/SOFT CAD cases vs. controls, respectively and the individuals’ mean age was 63 years in cases vs. 57 years in controls. 49.88%/50.04% vs. 41.52% of the individuals displayed Systolic blood pressure >140 mmHg. Hypercholesterolemia was reported by 51.20%/44.78% vs. 6.00% HARD/SOFT CAD cases vs. controls, respectively, while hypertriglyceridemia and poor glycemic control (defined as serum triglyceride levels >5.6 mmol/dL and HbA1c levels >64 mmol/mol, respectively) was found in 1.55%/1.43% and 3.59%/3.20% HARD/SOFT CAD cases vs. 0.80% (both traits) in controls, respectively. 20.80%, 19.41% and 4.14% had T2D. On average, the BMI of cases was 28.00 kg/m2 vs. 26.51 kg/m2 in controls and 35.17%/36.10% HARD/SOFT CAD cases vs. 22.65% controls displayed obesity (defined as BMI ≥ 30 kg/m2). The average body height of HARD/SOFT CAD cases was 173.92/174.16 cm (male) and 160.37/160.68 cm (female) vs. 175.95 cm and 162.61 cm in male and female controls, respectively. Sufficient physical activity (i.e., meeting the 2017 UK Physical activity guidelines of 150 minutes of moderate activity per week or 75 minutes of vigorous activity) was seen in 51.31%/51.08% and 54.60% of HARD/SOFT CAD cases vs. controls, respectively. 72.23%/69.56% vs. 60.21% of the individuals had a positive smoking history (answered ‘yes’ to ‘ever smoked’) and 13.74%/12.90% vs. 9.27% were current smokers. Finally, 59.14%, 57.35% and 41.41% had positive family history of heart disease. In all cases, the differences between CAD cases and controls were statistically significant (P<0.001, Table 1).

### 3.2. MT-SNV associations with HARD and SOFT CAD phenotypes

After quality control of genotyped individuals (including filtering for missingness by individual <0.1 and missingness by SNV<0.1, as described in section 2.2. in Methods), from the 265 MT-SNVs present on the UK Biobank or UK BiLEVE arrays, 243 remained for further analyses. For the genotyped common and low-frequency MT-SNVs (MAF>0.01; N=111, of those N=39 with MAF>0.05; genotyping rate >0.99) in UK Biobank, we performed single marker association analyses with HARD (N=20,405) and SOFT (N=34,782) CAD phenotypes, adjusting for the array, sex, birth year and first five principal components.

In HARD cases, no MT-SNVs survived multiple testing correction, the most significant (nominal P<0.05) findings (Table 2 and Fig.2) being for m.295C>T (rs41528348, P=0.0118, MAF=0.10, OR=1.05; 95% CI 1.02-1.09, in control region/CR, tagging macro-haplogroup J) and m.12612A>G (rs28359172, P=0.0158, MAF=0.10, OR=1.05; 95% CI 1.02-1.08, synonymous, in ND5 gene, tagging macro-haplogroup J), both more frequent in cases, thus potentially conferring increased CAD risk. In addition, four more MT-SNVs were found more frequently in controls: m.12372G>A (rs2853499, P=0.0059, MAF=0.22, OR=0.97; 95% CI 0.95-0.99, synonymous, in ND5 gene, tagging macro-haplogroup U), m.11467A>G (rs2853493, P=0.0065, MAF=0.22, OR=0.97; 95% CI 0.95-1.00, synonymous, in ND4 gene, tagging macro-haplogroup U), m.15301G>A (rs193302991, P=0.0115, MAF=0.04, OR=0.97; 95% CI 0.92-1.03, synonymous, in CYB gene) and m.7768A>G (rs41534044, P=0.0185, MAF=0.04, OR=0.91; 95% CI 0.86-0.96, synonymous, in CO2 gene). For a complete list of MT-SNV associations with HARD CAD phenotypes, see Supplementary Table 2.

![Figure 2](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/02/10/2022.02.09.22270723/F2.medium.gif)

[Figure 2](http://medrxiv.org/content/early/2022/02/10/2022.02.09.22270723/F2)

Figure 2 
A solar plot of HARD CAD and common and low-frequency (MAF>0.01; N=111) MT SNV associations.

View this table:
[Table 2](http://medrxiv.org/content/early/2022/02/10/2022.02.09.22270723/T2)

Table 2 
HARD CAD and common and low-frequency (MAF>0.01; N=111) MT SNV most significant associations. CR=control region.

In SOFT cases, four MT-SNVs survived multiple testing correction (at FDR<5%; Table 3 and Fig.3), all potentially conferring increased CAD risk: m.10400C>T (rs28358278, P=0.0007, MAF=0.02, OR=1.28; 95% CI 1.21-1.35, non-synonymous/Thr → Ala, in ND3 gene, tagging macro-haplogroup M), m.11251A>G (rs869096886, P=0.0011, MAF=0.20, OR=1.03; 95% CI 1.01-1.05, synonymous, in ND4 gene, tagging macro-haplogroups J and T), and two MT-SNVs in CYB gene - m.15452C>A (rs193302994, P=0.0017, MAF=0.20, OR=1.03; 95% CI 1.01-1.05, non-synonymous/Leu → Ile, tagging macro-haplogroups J and T) and m.15301G>A (rs193302991, P=0.0010, MAF=0.04, OR=1.03; 95% CI 0.99-1.07, synonymous). For a complete list of MT-SNVs associations with SOFT CAD phenotype, see Supplementary Table 3.

![Figure 3](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/02/10/2022.02.09.22270723/F3.medium.gif)

[Figure 3](http://medrxiv.org/content/early/2022/02/10/2022.02.09.22270723/F3)

Figure 3 
A solar plot of SOFT CAD and common and low-frequency (MAF >0.01; N=111) MT SNV associations. CR=control region.

View this table:
[Table 3](http://medrxiv.org/content/early/2022/02/10/2022.02.09.22270723/T3)

Table 3 
SOFT CAD and common and low-frequency (MAF>0.01; N=111) MT SNV most significant associations. CR=control region. *MT-SNVs that survived multiple testing correction (at FDR<5%).

### 3.3. MT gene-based associations with HARD and SOFT CAD phenotypes

We next sought to also consider the potential effects of rare (MAF ≤0.01) MT-SNVs on CAD risk, hence we performed MT-gene-based association analyses with HARD and SOFT CAD phenotypes, additionally including the whole mitochondrion as our region of interest (MT). As a result, we observed that in both, HARD and SOFT cases, CO2 displayed gene-based association at nominal significance (P<0.05), while CYB and ND3 were also associated (nominal P<0.05) with SOFT CAD phenotype (Table 4). When considering the whole mitochondrion (MT), no significant associations with CAD were observed (N=243; P=0.07, Table 4).

View this table:
[Table 4](http://medrxiv.org/content/early/2022/02/10/2022.02.09.22270723/T4)

Table 4 
MT gene-based associations with HARD and SOFT CAD phenotypes. MT=mitochondrion; N (ALL)=the number of MT-SNVs in the region; N (TESTED)=the number of MT-SNVs from the region considered in the gene-based test.

### 3.4. MT-haplogroup associations with HARD and SOFT CAD phenotypes

Different human mitochondrial haplogroups may result in differences in mitochondrial function that may influence susceptibility to CAD. Hence, we estimated all the mitochondrial haplogroups in our dataset (Table 5 and Supplementary Tables 4 and 5).

View this table:
[Table 5](http://medrxiv.org/content/early/2022/02/10/2022.02.09.22270723/T5)

Table 5 
Haplogroup assignment in HARD and SOFT CAD cases vs. controls prior to further assigning individuals to one of the major European haplogroups. * Indicates that a haplogorup survived multiple testing correction (at FDR<5%).

Three haplogroups survived multiple testing correction (at FDR<5%) in both HARD and SOFT cases vs. controls: M45a (0.59% and 0.55 % vs. 0.39%, OR=1.52 and OR=1.42, respectively), G2b1a2 (0.28% and 0.26% vs. 0.16%, OR=1.73 and OR=1.60, respectively) and U2b2 (0.11% and 0.09% vs. 0.04%, OR=2.56 and OR=2.31, respectively). In HARD cases, also haplogroup M57b1 was significantly (at FDR<5%) over-represented in cases vs. controls (0.02% vs. <0.01%, OR=33.06), while haplogroup L2c was significantly (at FDR<5%) under-represented cases vs. controls (0.01% vs.0.07%, OR=0.18) (Table 5 and Supplementary Table 4). In SOFT cases, also haplogroup M3a was significantly (at FDR<5%) over-represented in cases vs. controls (0.21% vs. 0.13%, OR=1.67, Table 5 and Supplementary Table 5).

When further assigning individuals to one of the major European haplogroups (Fig.4). As a result, 43.28%, 3.19%, 10.70%, 8.25%, 0.22%, 9.52%, 13.70%, 2.65%, 2.01%, 1.34% and 5.14% of individuals belonged do the haplogroup H, I, J, K, R, T, U, V, W, X or ‘Others’, respectively. Overall, the frequencies of the major European mitochondrial haplogroups did not differ significantly (at FDR<5%) between CAD patients and control subjects (Fig.4). Only the frequency of haplogroup I was significantly (nominal P<0.05) higher in patients with HARD CAD phenotype vs. controls (3.35% vs. 3.08%, OR=1.09) and the haplogroup R was significantly (nominal P<0.001 and P<0.01) higher in patients with HARD and SOFT CAD phenotype vs. controls (0.26% and 0.23% vs. 0.16%, OR=1.70 and OR=1.49, respectively; Fig.4).

![Figure 4](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/02/10/2022.02.09.22270723/F4.medium.gif)

[Figure 4](http://medrxiv.org/content/early/2022/02/10/2022.02.09.22270723/F4)

Figure 4 
Frequencies (%) of mitochondrial (MT) haplogroups within HARD and SOFT CAD phenotypes vs. controls.

## 4. Discussion

Over the last 14 years, several large-scale genome-wide association studies have found hundreds of single-nucleotide variants (SNVs) significantly associated with CAD, however, these explain <20% of the heritability. In this study, we hypothesize that mitochondrial (MT) SNVs might present one potential source of the “missing heritability”.

We analyzed 265 common/low-frequency (MAF ≥1%) and rare (MAF<1%) MT-SNVs in ∼500,000 UK Biobank individuals, exploring two different CAD definitions, HARD (N=20,405) and SOFT (N=34,782) (Fig.1.), as previously proposed by [77], and using the array, sex, birth year and first five principal components as covariates. Overall, the differences in the prevalence of common risk factors among CAD cases (both HARD and SOFT phenotypes) and controls were statistically significant (P<0.001; Table 1), male gender, older age, hypertension, hypercholesterolemia, obesity, T2D, physical inactivity, shorter body statue, smoking and positive family history demonstrating predominance in CAD patients.

When performing common and low-frequency MT SNVs (MAF>0.01; N=111) association analyses in these individuals, we observed that in HARD cases, no MT-SNVs survived multiple testing correction, the most significant (nominal P<0.05) findings being for m.295C>T, m.12612A>G, m.12372G>A, m.11467A>G, m.15301G>A and m.7768A>G (Table 2 and Fig.2). m.295C>T (rs41528348, P=0.0118, MAF=0.10), is a control region (CR) sequence variant, and the non-coding regions are well known to contain regulatory elements affecting gene expression [24]. m.295C>T has been previously associated with low maximal oxygen uptake (VO2max) in response to a standardized aerobic exercise training program [100]. VO2max is a heritable trait, demonstrating a maternal effect [24]. VO2max is defined by the ability of the cardiorespiratory system to deliver oxygen to the exercising muscles, e.g., the heart rate and maximal cardiac output [7]. High VO2max (i.e., cardiorespiratory fitness) have been associated with a reduction in cardiovascular events and vice versa - low VO2max plays a role in increasing the risk of cardiovascular disease [54, 60, 93, 107]. Interestingly, an earlier study in UK Biobank [108], reported a significant association between m.295C>T and mean corpuscular hemoglobin (MCH) as well as mean corpuscular volume (MCV) (Fig.5). Both traits have been shown to increase with training [10]. In line with these earlier reports, our results show that m.295C>T was more frequent (OR=1.05; 95% CI 1.02-1.09, P=0.0118) in HARD cases, thus, potentially conferring a decreased cardiorespiratory fitness/exercise capacity and increased CAD risk. Of note, m.295C>T is also tagging macro-haplogroup J, and it has been previously shown that VO2max is lower in J than in non-J macro-haplogroup individuals [72] and excellence in endurance performance was less frequent among macro-haplogroup J individuals [52].

![Figure 5](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/02/10/2022.02.09.22270723/F5.medium.gif)

[Figure 5](http://medrxiv.org/content/early/2022/02/10/2022.02.09.22270723/F5)

Figure 5 
Previous findings for HARD and SOFT CAD-MT-SNV associations in UK Biobank by ref. by ref. [108]. MCH: Mean corpuscular hemoglobin; MCV: Mean corpuscular volume; RDW: Red blood cell distribution width; RBC#: Red blood cell count; WBC#: White blood cell count; LYMPH#: Lymphocyte count; MONO%: Monocyte percentage of white cells; MCV: Mean corpuscular volume; PCT: Plateletcrit; Cr: Creatinine; eGFRCr: Estimated glomerular filtration rate creatinine; eGFRCrCy: Estimated glomerular filtration rate creatinine and cystatin C. *Created with [BioRender.com](http://BioRender.com)*.

m.12612A>G (rs28359172, P=0.0158, MAF=0.10) is a synonymous (V92V) sequence variant in the subunit 5 of NADH dehydrogenase (ND5), also tagging macro-haplogroup J and demonstrating significant associations with MCH and MCV (Fig.5) in UK Biobank [108]. In our study, it was also more frequent (OR=1.05; 95% CI 1.02-1.08, P=0.0158) in HARD cases. Therefore, it is reasonable to assume that m.295C>T also may be related to a decreased cardiorespiratory fitness/exercise ability and increased CAD risk. m.12372G>A (rs2853499, P=0.0059, MAF=0.22) is a synonymous (L12L) sequence variant in the ND5 gene, found more often in controls vs. HARD CAD cases (OR=0.97; 95% CI 0.95-0.99) and tagging macro-haplogroup U, members of which have been shown to represented higher brain pH than others [83]. m.11467A>G (rs2853493, P=0.0065, MAF=0.22) represents another synonymous (L236L) sequence variant also tagging macro-haplogroup U and more frequently found in controls vs. HARD CAD cases (OR=0.97; 95% CI 0.95-1.00), however, it is located in the subunit 4 of NADH dehydrogenase (ND4). m.11467A>G (rs2853493) has also been also associated with altered brain pH [83]. Higher brain pH has been shown to correlate with a lower brain activity and conferring protection against psychiatric disorders [41] and attention deficit/hyperactivity [15]. This protective role of macro-haplogroup U could be partially explained by reduced load of harmful reactions [15], as pH is known to play a role in mitochondrial ROS generation [90] and cybrids of macro-haplogroup U have been demonstrated to display lower inner mitochondrial membrane potential and mitochondrial protein synthesis and cell growth capacity [39]. Moreover, macro-haplogroup U has been associated with increased longevity that may be at least partly explained by its lower ROS production rate [17]. In our study, m.12372G>A (rs2853499) was found more frequently in controls, suggesting potential favorable differences in pH and ROS production in these subjects. Interestingly, in addition to associations with MCH and MCV, m.12372G>A displayed significant (P< 1e - 5) associations with eight additional blood cell and kidney-related traits in UK Biobank (Fig.5) [108]. Endurance time during exercise has been related to pre-exercise blood pH and demonstrated to increase with increasing pH [47]. m.7768A>G (rs41534044, P=0.0185, MAF=0.04) represents another synonymous (M61M) sequence variant located in the subunit 2 of mitochondrial encoded cytochrome c oxidase (CO2), more frequently found in controls vs. HARD CAD cases (OR=0.91; 95% CI 0.86-0.96) and displaying significant (P< 1e - 5) associations with several phenotypic traits in UK Biobank [108], including MCH and MCV, as well as the Monocyte percentage of white cells (MONO%) and Plateletcrit (PCT)(Fig.5).

In SOFT cases, four MT-SNVs survived multiple testing correction (at FDR<5%; Table 3 and Fig.3), all potentially conferring increased CAD risk: m.10400C>T, m.11251A>G, m.15452C>A and 15301G>A. m.11251A>G (rs869096886, P=0.0011, MAF=0.20) represents a synonymous sequence variant in ND4 gene and m.15452C>A (rs193302994, P=0.0017, MAF=0.20) is a non-synonymous (Leu → Ile) sequence variant in CYB gene, both were found more frequently in SOFT CAD cases vs. controls (OR=1.03; 95% CI 1.01-1.05). m.11251A>G (rs869096886) and m.15452C>A (rs193302994) are tagging macro-haplogroups J and T, potentially related to a decreased cardiorespiratory fitness/exercise capacity [52, 72] and thus increased CAD risk, as discussed above. Moreover, both MT-SNVs displayed significant (P< 1e - 5) associations with body height in UK Biobank (Fig.5) [108]. In line with this, we observed that the average body height of both male and female CAD cases was ∼2.00 cm lower compared to controls (Table 1). Previous studies have demonstrated that shorter body height is related to an increased CAD risk [78], which could be partly explained by altered lipid profile.

Gene-based tests revealed that in both, HARD and SOFT cases, CO2 displayed gene-based association at nominal significance (P<0.05, Table 4). The CO2 gene encodes for the second subunit of cytochrome c oxidase (COX, complex IV). Dysfunction of COX has been previously associated with mitochondrial oxidative stress, obesity and T2D [43]. CYB and ND3 were also associated (nominal P<0.05) with SOFT CAD phenotype (Table 4). Somatic variations in CYB have been previously related to hypertrophic cardiomyopathy (one of its clinical manifestations being angina) and exercise intolerance [40]. Recently, a large gene-based meta-analysis of mitochondrial genes with several metabolic traits identified ND3 associated with BMI (P<1e-03) [57].

All haplogroups demonstrating significant (at FDR<5%) associations in our study (M45a, G2b1a2, U2b2 with both HARD/SOFT, M57b1 and L2c (under-represented) with HARD and M3a with SOFT CAD phenotypes; Table 5 and Supplementary Table 4) were with a frequency <1%, while other studies have considered only haplogroups with a frequency ≥5% [26]. Low counts in the less common haplogroups may lead to a false-positive result [86]. Although, this should be addressed by performing multiple testing correction, grouping the less frequent haplogroups may be another approach to tackle this [86]. Hence, we also assigned individuals to one of the major European haplogroups (Fig.4) for comparison. As a result, we observed that 43.28% of the individuals belonged to the macro-haplogroup H, 13.70% to the macro-haplogroup U, 10.70% to the macro-haplogroup J and 9.52% to the macro-haplogroup T (Fig.4), in line with earlier reports in other European populations [2]. Overall, the frequencies of the major European mitochondrial haplogroups did not differ significantly (at FDR<5%) between CAD patients and control subjects (Fig.4). Only the frequency of haplogroup I was significantly (nominal P<0.05) higher in patients with HARD CAD phenotype vs. controls (3.35% vs. 3.08%, OR=1.09) and the macro-haplogroup R was significantly (nominal P<0.001 and P<0.01) higher in patients with HARD and SOFT CAD phenotype vs. controls (0.26% and 0.23% vs. 0.16%, OR=1.70 and OR=1.49, respectively) (Fig.4). A recent study [108] identified a significant effect of macro-haplogroup I on the MCH levels in UK Biobank. As previously discussed, this trait is known to increase with training [10], thus, again being potentially related to cardiorespiratory fitness/exercise capacity and thus CAD risk. Of note, however, we were able to assign most samples reliably into haplogroups as the mtDNA haplogroups were deduced from genotyping arrays with limited numbers of (high quality) SNVs being profiled, hence the quality score for haplogroup assignment ranged from 0.50 to 0.86, with a median of 0.68. Hence, we were not able to exclude samples with quality scores for haplogroup assignment <0.8 (as in ref. [36]). Moreover, we performed Fisher’s exact test that did not allow us to adjust for co-variates, hence it is possible that known and unknown potential confounding factors might have influenced these results. Though, in which cases to adjust for which covariates and whether it will increase or decrease the study power and/or bias, is still a matter of intense debate [5, 99].

Several other limitations should also be acknowledged. It is well-known that exceptionally large cohorts are required to reliably associate genetic variations with complex traits 86. The power for detecting causal MT-SNVs and haplogroups has been compared with that in the nuclear genome given equal effect sizes, estimating that the sample size needed for the mitochondrial studies would be roughly the same as that needed for the nuclear genome studies 73. Previous power calculations for ischemic stroke (assuming an additive model) [3] revealed a maximum power of 73% to detect SNVs with OR=1.4 and MAF=0.30, while for SNVs conferring OR=1.20 and MAF=0.20, the study power dropped to 4.6% and further to 0.001% for OR=1.10 and MAF=0.10. This study concluded that “prohibitively large sample sizes” would be needed to achieve sufficient power to detect individual MT DNA variants [3]. In line with this, we observe that in HARD CAD cases, where N=20,405, no MT-SNVs survived multiple testing correction, while when increasing N to 34,782 in SOFT CAD cases, four MT-SNVs survived multiple testing correction (FDR<5%). Hence, even larger sample sizes may be needed to reliably associate MT SNVs and haplogroups with CAD.

In addition to study individuals, also the number of MT-SNVs studied was limited. In UK Biobank [13], genotyping was performed using Affymetrix UK biobank Axiom (450,000 samples) and Affymetrix UK BiLEVE Axiom (50,000 samples) arrays, which included 265 genotyped MT DNA variants. After quality control procedures (Fig.1), 243 MT-SNVs remained for further analyses, 111 of those where common or low-frequency (MAF>0.01) and could be used for single marker association analyses. However, this is clearly not a representative set of MT-SNVs and, as previously recognized, some regions may be not well-covered, such as the hypervariable regions [32, 34]. Clearly, whole genome sequencing or targeted-sequencing of MT-DNA, considering their ability to achieve a deep genome coverage, would allow the identification many more MT-SNV (especially the low-frequency/rare variants; MAF<0.01), also improving the detection of haplogroups and allowing to investigate heteroplasmy, a phenomenon characteristic to MT DNA [18, 92].

Heteroplasmy denotes the coexistence of MT DNA genomes with wild-type inherited SNVs and somatic variants in varying ratios, which are dynamically determined and display patterns of cell and tissue specificity, and may differ even within the same mitochondrion [92], determining the clinical presentation of disease phenotypes [32, 33, 81]. In this study, we were limited to genotype calls from arrays, which are restricted in terms of minor alleles and do not allow to capture heteroplasmy [32, 33]. Moreover, MT DNA content was assessed only in blood cells, whereas previous studies have identified additional six vascular and metabolic tissues relevant to CAD [35, 94]. Therefore, whole genome sequencing/targeted-sequencing of MT-DNA, across several vascular and metabolic tissues relevant to CAD may be necessary [35, 94, 101] in order to characterize the full landscape of mitochondrial genetic variations and their potential contribution to these complex disease phenotypes. Especially, considering that the energy requirements and thus sensitivity to the changes in mitochondrial function differ for different cells and tissues and hence may be important in determining the phenotypic effect of MT SNVs [11].

We also do not consider mitochondrial DNA copy number (MT DNA-CN), representing the number of mitochondria per cell and the number of MT DNA per mitochondrion [6, 19]. Each mitochondrion contains multiple copies of MT DNA and different cells and tissues contain different numbers (up to 7,000) of mitochondria, again displaying patterns of cell and tissue variability [6, 19]. MT DNA-CN is believed to serve as an indirect biomarker that would capture the underlying mitochondrial activity and function, such as energy production capacity and metabolic characteristics, thus possibly playing a causative role in health and disease [19]. Decreased MT DNA-CN has been previously associated with an increased risk of developing cardiovascular disease (CVD) outcomes [6]. More recently, similar analyses in the UK Biobank demonstrated a possible causal role of lower MT DNA-CN on higher CAD risk [67]. In an even larger cohort (of 408,361 individuals from TOPMed and UK Biobank), a decline in MT DNA-CN was observed in elderly individuals (>65 years) and lower MT DNA-CN levels also demonstrated an age-independent associations with hypertension, hyperlipidemia, T2D and obesity, i.e., the well-known CAD risk factors [64]. However, none of these studies compared the MT DNA-CN levels between HARD vs. SOFT CAD phenotypes, which could be a subject of future studies. However, considering that MT DNA-CN varies greatly across cell and tissue types, again profiling of several vascular and metabolic tissues relevant to CAD may be necessary for such investigations [35, 94, 101].

Yet another important aspect not considered here is the nuclear genome, considering the co-evolution of mitochondria and eukaryotic cells [74]. The mitochondrial genome encodes only 37 genes, mainly components of the OXPHOS machinery, whereas the remaining ∼1000-1500 mitochondrial proteins are all encoded by the nuclear genome [42]. The importance of common genetic variation in the nuclear genome regulating MT heteroplasmy and DNA-CN is an active area of research [19, 38, 46, 75]. Moreover, genetic variants in nuclear genes could lead to oxidative disorders or modulate the mitochondrial variants [81] and mild nuclear gene variants could potentially become clinically relevant when combined with an incompatible MT DNA [105]. Additive interactions (epistasis) between mitochondrial variants in MT-ND2 gene and nuclear variants in genes responsible for mitochondrial replication and transcription have been demonstrated to influence the BMI and obesity phenotype [66]. Similarly, our previous investigations have demonstrated the role of nuclear encoded mitochondria imported genes in coordinating the response to hypercholesterolemia and atherosclerotic lesion expansion, and foam cell formation [102]. Hence, further analyses also considering these additional variations will be required.

Finally, similar to SNVs in nuclear genome, even if (mitochondria) genome-wide significant associations with HARD/SOFT CAD phenotypes would be identified, their functional consequences would need to be determined in the CAD relevant tissues [35, 94]. Currently, functional studies for MT-SNVs are not readily available, however, several novel experimental animal models (e.g. mice strains displaying DNA haplogroups similar to those observed in humans) may be available in the near future, allowing to investigate/define the potential causality of the relationship between inherited “natural” non-pathogenic MT-SNVs and potential alterations in mitochondrial function (e.g. oxygen consumption and oxidant production, cellular ATP levels) and their relation to alterations in cardiovascular function and CAD risk [11, 74, 81].

In conclusion, we found only spurious MT-SNV, gene and haplogorup associations with HARD and SOFT CAD phenotypes and conclude that whole genome sequencing/targeted-sequencing of MT-DNA, across several vascular and metabolic tissues relevant to CAD in even larger study cohorts (N>50,000), followed by functional studies in animal models, may be necessary to conclusively determine the role of MT-SNVs, genes and haplogorup in modulating the risk of CAD. Therefore, whole genome sequencing/targeted-sequencing of MT-DNA, across several vascular and metabolic tissues relevant to CAD may be necessary [35, 94, 101] in order to characterize the full landscape of mitochondrial genetic variations and their potential contribution to these complex disease phenotypes. Especially, considering that the energy requirements and thus sensitivity to the changes in mitochondrial function differ for different cells and tissues and hence may be important in determining the phenotypic effect of MT SNVs [11].

## Supporting information

Supplementary Tables [[supplements/270723_file12.xls]](pending:yes)

## Data Availability

All data produced will be submitted to the UK Biobank Resource ([https://www.ukbiobank.ac.uk/](https://www.ukbiobank.ac.uk/)) after successful publishing of this manuscript. 

## Disclosures

None.

## Acknowledgments and Sources of Funding

This research was funded by the Latvian Council of Science, project, “Using Machine Learning to Model the Complex Interplay Between Diet, Genetic Factors and Mitochondria in Coronary Artery Disease”, project No. lzp-2020/2-0338. The present study has been conducted using the UK Biobank Resource (application ID 61684) that is available to researchers.

*   Received February 9, 2022.
*   Revision received February 9, 2022.
*   Accepted February 10, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  [1].Aherrahrou R, Aherrahrou Z, Schunkert H, Erdmann J (2017) Coronary artery disease associated gene Phactr1 modulates severity of vascular calcification in vitro. Biochemical and biophysical research communications 491: 396–402
    
    

2.  [2].Aldámiz-Echevarría T, Resino S, Bellón JM, Jiménez-Sousa MA, Miralles P, Medrano LM, Carrero A, Díez C, Pérez-Latorre L, Fanciulli C, Garcia-Broncano P, Berenguer J (2019) European mitochondrial haplogroups predict liver-related outcomes in patients coinfected with HIV and HCV: a retrospective study. Journal of Translational Medicine 17
    
    

3.  [3].Anderson CD, Biffi A, Nalls MA, Devan WJ, Schwab K, Ayres AM, Valant V, Ross OA, Rost NS, Saxena R, Viswanathan A, Worrall BB, Brott TG, Goldstein JN, Brown D, Broderick JP, Norrving B, Greenberg SM, Silliman SL, Hansen BM, et al. (2013) Common variants within oxidative phosphorylation genes influence risk of ischemic stroke and intracerebral hemorrhage. Stroke 44: 612–619
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToic3Ryb2tlYWhhIjtzOjU6InJlc2lkIjtzOjg6IjQ0LzMvNjEyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDIvMTAvMjAyMi4wMi4wOS4yMjI3MDcyMy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

4.  [4].Anderson S, Bankier AT, Barrell BG, de Bruijn MH, Coulson AR, Drouin J, Eperon IC, Nierlich DP, Roe BA, Sanger F, Schreier PH, Smith AJ, Staden R, Young IG (1981) Sequence and organization of the human mitochondrial genome. Nature 290: 457–465
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/290457a0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7219534&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1981LK04700029&link_type=ISI) 

5.  [5].Aschard H, Vilhjálmsson BJ, Joshi AD, Price AL, Kraft P (2015) Adjusting for heritable covariates can bias effect estimates in genome-wide association studies. American journal of human genetics 96: 329–339
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2014.12.021&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25640676&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

6.  [6].Ashar FN, Zhang Y, Longchamps RJ, Lane J, Moes A, Grove ML, Mychaleckyj JC, Taylor KD, Coresh J, Rotter JI, Boerwinkle E, Pankratz N, Guallar E, Arking DE (2017) Association of Mitochondrial DNA Copy Number with Cardiovascular Disease. JAMA Cardiology 2: 1247
    
    

7.  [7].Bassett DR (2000) Limiting factors for maximum oxygen uptake and determinants of endurance performance. Medicine Science in Sports Exercise: 70
    
    

8.  [8].Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc
    
    

9.  [9].Børte S, Zwart JA, Skogholt AH, Gabrielsen ME, Thomas LF, Fritsche LG, Surakka I, Nielsen JB, Zhou W, Wolford BN, Vigeland MD, Hagen K, Kristoffersen ES, Nyholt DR, Chasman DI, Brumpton BM, Willer CJ, Winsvold BS (2020) Mitochondrial genome-wide association study of migraine – the HUNT Study. Cephalalgia 40: 625–634
    
    

10. [10].Branch JD, Pate RR, Bourque SP, Convertino VA, Durstine JL, Ward DS (1997) Effects of exercise mode on hematologic adaptations to endurance training in adult females. Aviation space and environmental medicine 68: 788–794
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9293346&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

11. [11].Bray AW, Ballinger SW (2017) Mitochondrial DNA mutations and cardiovascular disease. Current Opinion in Cardiology 32: 267–274
    
    

12. [12].Brænne I, Civelek M, Vilne B, Di Narzo A, Johnson AD, Zhao Y, Reiz B, Codoni V, Webb TR, Foroughi Asl H, Hamby SE, Zeng L, Trégouët DA, Hao K, Topol EJ, Schadt EE, Yang X, Samani NJ, Björkegren JLM, Erdmann J, et al. (2015) Prediction of Causal Candidate Genes in Coronary Artery Disease Loci. Arteriosclerosis thrombosis and vascular biology 35: 2207– 2217
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYXR2YmFoYSI7czo1OiJyZXNpZCI7czoxMDoiMzUvMTAvMjIwNyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzEwLzIwMjIuMDIuMDkuMjIyNzA3MjMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

13. [13].Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, Motyer A, Vukcevic D, Delaneau O, O’Connell J, Cortes A, Welsh S, Young A, Effingham M, McVean G, Leslie S, Allen N, Donnelly P, Marchini J (2018) The UK Biobank resource with deep phenotyping and genomic data. Nature 562: 203–209
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-018-0579-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30305743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

14. [14].Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ (2015) Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience 4: 7, original DateCompleted: 20150227
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13742-015-0047-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25722852&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

15. [15].Chang X, Liu Y, Mentch F, Glessner J, Qu H, Nguyen K, Sleiman PMA, Hakonarson H (2020) Mitochondrial DNA haplogroups and risk of attention deficit and hyperactivity disorder in European Americans. Translational psychiatry 10: 370
    
    

16. [16].Chen A, Raule N, Chomyn A, Attardi G (2012) Decreased reactive oxygen species production in cells with mitochondrial haplogroups associated with longevity. PloS one 7: e46473
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0046473&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23144696&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

17. [17].Chen A, Raule N, Chomyn A, Attardi G (2012) Decreased reactive oxygen species production in cells with mitochondrial haplogroups associated with longevity. PloS one 7: e46473
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0046473&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23144696&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

18. [18].Chen R, Aldred MA, Xu W, Zein J, Bazeley P, Comhair SAA, Meyers DA, Bleecker ER, Liu C, Erzurum SC, Hu B, (SARP) NSARP (2021) Comparison of whole genome sequencing and targeted sequencing for mitochondrial DNA. Mitochondrion 58: 303–310
    
    

19. [19].Chong M, Mohammadi-Shemirani P, Perrot N, Nelson W, Morton RW, Narula S, Lali R, Khan I, Khan M, Judge C, Machipisa T, Cawte N, O’Donnell M, Pigeyre M, Akhabir L, Paré G (2021) GWAS and ExWAS of blood Mitochondrial DNA copy number identifies 73 loci and highlights a potential causal role in dementia
    
    

20. [20].The CARDIoGRAMplusC4D Consortium, Deloukas P, Kanoni S, Willenborg C, Farrall M, Assimes TL, Thompson JR, Ingelsson E, Saleheen D, Erdmann J, Goldstein BA, Stirrups K, König IR, Cazier JB, Johansson A, Hall AS, Lee JY, Willer CJ, Chambers JC, Esko T, et al. (2013) Large-scale association analysis identifies new risk loci for coronary artery disease. Nat Genet 45: 25–33
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2480&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23202125&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

21. [21].Cordain L, Eaton SB, Sebastian A, Mann N, Lindeberg S, Watkins BA, O’Keefe JH, Brand-Miller J (2005) Origins and evolution of the Western diet: health implications for the 21st century. The American journal of clinical nutrition 81: 341–354
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpjbiI7czo1OiJyZXNpZCI7czo4OiI4MS8yLzM0MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzEwLzIwMjIuMDIuMDkuMjIyNzA3MjMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

22. [22].Davey Smith G, Ebrahim S, Lewis S, Hansell AL, Palmer LJ, Burton PR (2005) Genetic epidemiology and public health: hope, hype, and future prospects. Lancet London England 366: 1484–1498
    
    

23. [23].De Benedictis G, Rose G, Carrieri G, De Luca M, Falcone E, Passarino G, Bonafe M, Monti D, Baggio G, Bertolini S, Mari D, Mattace R, Franceschi C (1999) Mitochondrial DNA inherited variants are associated with successful aging and longevity in humans. FASEB journal official publication of the Federation of American Societies for Experimental Biology 13: 1532–1536
    
    

24. [24].Dionne FT, Turcotte L, Thibault MC, Boulay MR, Skinner JS, Bouchard C (1993) Mitochondrial DNA sequence polymorphism, VO2max, and response to endurance training. Medicine and science in sports and exercise 25: 766–774
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1249/00005768-199307000-00002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8350696&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1993LM59200001&link_type=ISI) 

25. [25].Earp MA, Brooks-Wilson A, Cook L, Le N (2013) Inherited common variants in mitochondrial DNA and invasive serous epithelial ovarian cancer risk. BMC Res Notes 6: 425
    
    

26. [26].Ebner S, Lang R, Mueller EE, Eder W, Oeller M, Moser A, Koller J, Paulweber B, Mayr JA, Sperl W, Kofler B (2011) Mitochondrial Haplogroups, Control Region Polymorphisms and Malignant Melanoma: A Study in Middle European Caucasians. PLoS ONE 6: e27192
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0027192&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22174736&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

27. [27].Erdmann J, Grosshennig A, Braund PS, König IR, Hengstenberg C, Hall AS, Linsel-Nitschke P, Kathiresan S, Wright B, Trégouët DA, Cambien F, Bruse P, Aherrahrou Z, Wagner AK, Stark K, Schwartz SM, Salomaa V, Elosua R, Melander O, Voight BF, et al. (2009) New susceptibility locus for coronary artery disease on chromosome 3q22.3. Nat Genet 41: 280– 282
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.307&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19198612&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000263640200006&link_type=ISI) 

28. [28].Erdmann J, Kessler T, Munoz Venegas L, Schunkert H (2018) A decade of genome-wide association studies for coronary artery disease: the challenges ahead. Cardiovascular research 114: 1241–1257
    
    

29. [29].Erdmann J, Willenborg C, Nahrstaedt J, Preuss M, König IR, Baumert J, Linsel-Nitschke P, Gieger C, Tennstedt S, Belcredi P, Aherrahrou Z, Klopp N, Loley C, Stark K, Hengstenberg C, Bruse P, Freyer J, Wagner AK, Medack A, Lieb W, et al. (2011) Genome-wide association study identifies a new locus for coronary artery disease on chromosome 10p11.23. Eur Heart J 32: 158–168
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/eurheartj/ehq405&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21088011&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000286215500009&link_type=ISI) 

30. [30].Esposito LA, Melov S, Panov A, Cottrell BA, Wallace DC (1999) Mitochondrial disease in mouse results in increased oxidative stress. Proc Natl Acad Sci U S A 96: 4820–4825
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czo5OiI5Ni85LzQ4MjAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMi8xMC8yMDIyLjAyLjA5LjIyMjcwNzIzLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

31. [31].Fisher RA (1922) On the interpretation of X2 from contingency tables, and the calculation of P. Journal of the Royal Statistical Society 85 (1): 87–94
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/2340521&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000189066700003&link_type=ISI) 

32. [32].Flaquer A, Baumbach C, Kriebel J, Meitinger T, Peters A, Waldenberger M, Grallert H, Strauch K (2014) Mitochondrial genetic variants identified to be associated with BMI in adults. PLoS One 9: e105116
    
    

33. [33].Flaquer A, Baumbach C, Ladwig KH, Kriebel J, Waldenberger M, Grallert H, Baumert J, Meitinger T, Kruse J, Peters A, Emeny R, Strauch K (2015) Mitochondrial genetic variants identified to be associated with posttraumatic stress disorder. Transl Psychiatry 5: e524
    
    

34. [34].Flaquer A, Rospleszcz S, Reischl E, Zeilinger S, Prokisch H, Meitinger T, Meisinger C, Peters A, Waldenberger M, Grallert H, Strauch K (2015) Mitochondrial GWA Analysis of Lipid Profile Identifies Genetic Variants to Be Associated with HDL Cholesterol and Triglyceride Levels. PLoS One 10: e0126294
    
    

35. [35].Franzén O, Ermel R, Cohain A, Akers NK, Di Narzo A, Talukdar HA, Foroughi-Asl H, Giambartolomei C, Fullard JF, Sukhavasi K, Köks S, Gan LM, Giannarelli C, Kovacic JC, Betsholtz C, Losic B, Michoel T, Hao K, Roussos P, Skogsberg J, et al. (2016) Cardiometabolic risk loci share downstream cis- and trans-gene regulation across tissues and diseases. Science New York NY 353: 827–830
    
    

36. [36].Gonalves VF, Giamberardino SN, Crowley JJ, Vawter MP, Saxena R, Bulik CM, Yilmaz Z, Hultman CM, Sklar P, Kennedy JL, Sullivan PF, Knight J (2018) Examining the role of common and rare mitochondrial variants in schizophrenia. PloS one 13: e0191153
    
    

37. [37].Gutierrez J, Ballinger SW, Darley-Usmar VM, Landar A (2006) Free radicals, mitochondria, and oxidized lipids: the emerging role in signal transduction in vascular cells. Circ Res 99: 924–932
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNpcmNyZXNhaGEiO3M6NToicmVzaWQiO3M6ODoiOTkvOS85MjQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMi8xMC8yMDIyLjAyLjA5LjIyMjcwNzIzLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

38. [38].Guyatt AL, Brennan RR, Burrows K, Guthrie PAI, Ascione R, Ring SM, Gaunt TR, Pyle A, Cordell HJ, Lawlor DA, Chinnery PF, Hudson G, Rodriguez S (2019) A genome-wide association study of mitochondrial DNA copy number in two population-based cohorts. Human genomics 13: 6
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s40246-018-0190-2&link_type=DOI) 

39. [39].Gómez-Durán A, Pacheu-Grau D, López-Gallardo E, Díez-Sánchez C, Montoya J, López-Pérez MJ, Ruiz-Pesini E (2010) Unmasking the causes of multifactorial disorders: OXPHOS differences between mitochondrial haplogroups. Human molecular genetics 19: 3343–3353
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddq246&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20566709&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000280704800005&link_type=ISI) 

40. [40].Hagen CM, Aidt FH, Havndrup O, Hedley PL, Jespersgaard C, Jensen M, Kanters JK, Moolman-Smook JC, Møller DV, Bundgaard H, Christiansen M (2013) MT-CYB mutations in hypertrophic cardiomyopathy. Molecular Genetics Genomic Medicine 1: 54–65
    
    

41. [41].Hagihara H, Catts VS, Katayama Y, Shoji H, Takagi T, Huang FL, Nakao A, Mori Y, Huang KP, Ishii S, Graef IA, Nakayama KI, Shannon Weickert C, Miyakawa T (2018) Decreased Brain pH as a Shared Endophenotype of Psychiatric Disorders. Neuropsychopharmacology official publication of the American College of Neuropsychopharmacology 43: 459–468
    
    

42. [42].Hock MB, Kralli A (2009) Transcriptional control of mitochondrial biogenesis and function. Annual review of physiology 71: 177–203
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.physiol.010908.163119&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19575678&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264489600009&link_type=ISI) 

43. [43].Holvoet P, Vanhaverbeke M, Bloch K, Baatsen P, Sinnaeve P, Janssens S (2016) Low MT-CO1 in Monocytes and Microvesicles Is Associated with Outcome in Patients with Coronary Artery Disease. Journal of the American Heart Association 5
    
    

44. [44].Howson JMM, Zhao W, Barnes DR, Ho WK, Young R, Paul DS, Waite LL, Freitag DF, Fauman EB, Salfati EL, Sun BB, Eicher JD, Johnson AD, Sheu WHH, Nielsen SF, Lin WY, Surendran P, Malarstig A, Wilk JB, Tybjaerg-Hansen A, et al. (2017) Fifteen new risk loci for coronary artery disease highlight arterial-wall-specific mechanisms. Nature genetics 49: 1113–1119
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3874&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28530674&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

45. [45].Hudson G, Gomez-Duran A, Wilson IJ, Chinnery PF (2014) Recent mitochondrial DNA mutations increase the risk of developing common late-onset human diseases. PLoS Genet 10: e1004369
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1004369&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24852434&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

46. [46].Hägg S, Jylhävä J, Wang Y, Czene K, Grassmann F (2020) Deciphering the genetic and epidemiological landscape of mitochondrial DNA abundance. Human Genetics 140: 849–861
    
    

47. [47].Jones NL, Sutton JR, Taylor R, Toews CJ (1977) Effect of pH on cardiorespiratory and metabolic responses to exercise. Journal of applied physiology respiratory environmental and exercise physiology 43: 959–964
    
    

48. [48].Kazuno Aa, Munakata K, Nagai T, Shimozono S, Tanaka M, Yoneda M, Kato N, Miyawaki A, Kato T (2006) Identification of mitochondrial DNA polymorphisms that alter mitochondrial matrix pH and intracellular calcium dynamics. PLoS genetics 2: e128
    
    

49. [49].Kessler T, Vilne B, Schunkert H (2016) The impact of genome-wide association studies on the pathophysiology and therapy of cardiovascular disease. EMBO molecular medicine 8: 688–701
    
    

50. [50].Kessler T, Wobst J, Wolf B, Eckhold J, Vilne B, Hollstein R, von Ameln S, Dang TA, Sager HB, Moritz Rumpf P, Aherrahrou R, Kastrati A, Björkegren JLM, Erdmann J, Lusis AJ, Civelek M, Kaiser FJ, Schunkert H (2017) Functional Characterization of the GUCY1A3 Coronary Artery Disease Risk Locus. Circulation 136: 476–489
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjk6IjEzNi81LzQ3NiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzEwLzIwMjIuMDIuMDkuMjIyNzA3MjMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

51. [51].Kessler T, Zhang L, Liu Z, Yin X, Huang Y, Wang Y, Fu Y, Mayr M, Ge Q, Xu Q, Zhu Y, Wang X, Schmidt K, de Wit C, Erdmann J, Schunkert H, Aherrahrou Z, Kong W (2015) ADAMTS-7 inhibits re-endothelialization of injured arteries and promotes vascular remodeling through cleavage of thrombospondin-1. Circulation 131: 1191–1201
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjExOiIxMzEvMTMvMTE5MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzEwLzIwMjIuMDIuMDkuMjIyNzA3MjMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

52. [52].Kiiskilä J, Jokelainen J, Kytövuori L, Mikkola I, Härkönen P, Keinänen-Kiukaanniemi S, Majamaa K (2021) Association of mitochondrial DNA haplogroups J and K with low response in exercise training among Finnish military conscripts. BMC Genomics 22
    
    

53. [53].Kloss-Brandstätter A, Pacher D, Schönherr S, Weissensteiner H, Binna R, Specht G, Kronenberg F (2011) HaploGrep: a fast and reliable algorithm for automatic classification of mitochondrial DNA haplogroups. Hum Mutat 32: 25–32
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/humu.21382&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20960467&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

54. [54].Kodama S (2009) Cardiorespiratory Fitness as a Quantitative Predictor of All-Cause Mortality and Cardiovascular Events in Healthy Men and Women. JAMA 301: 2024
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2009.681&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19454641&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266159600027&link_type=ISI) 

55. [55].Kofler B, Mueller EE, Eder W, Stanger O, Maier R, Weger M, Haas A, Winker R, Schmut O, Paulweber B, Iglseder B, Renner W, Wiesbauer M, Aigner I, Santic D, Zimmermann FA, Mayr JA, Sperl W (2009) Mitochondrial DNA haplogroup T is associated with coronary artery disease and diabetic retinopathy: a case control study. BMC Med Genet 10: 35
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2350-10-35&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19383124&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

56. [56].Kovacic JC (2017) Unraveling the Complex Genetics of Coronary Artery Disease *. Journal of the American College of Cardiology 69: 837–840
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czo0OiJhY2NqIjtzOjU6InJlc2lkIjtzOjg6IjY5LzcvODM3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDIvMTAvMjAyMi4wMi4wOS4yMjI3MDcyMy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

57. [57].Kraja AT, Liu C, Fetterman JL, Graff M, Have CT, Gu C, Yanek LR, Feitosa MF, Arking DE, Chasman DI, Young K, Ligthart S, Hill WD, Weiss S, Luan J, Giulianini F, Li-Gao R, Hartwig FP, Lin SJ, Wang L, et al. (2019) Associations of Mitochondrial and Nuclear Mitochondrial Variants and Genes with Seven Metabolic Traits. American journal of human genetics 104: 112–138
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2018.12.001&link_type=DOI) 

58. [58].Krzywanski DM, Moellering DR, Fetterman JL, Dunham-Snary KJ, Sammy MJ, Ballinger SW (2011) The mitochondrial paradigm for cardiovascular disease susceptibility and cellular function: a complementary concept to Mendelian genetics. Laboratory investigation a journal of technical methods and pathology 91: 1122–1135
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/labinvest.2011.95&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21647091&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

59. [59].Krzywanski DM, Moellering DR, Westbrook DG, Dunham-Snary KJ, Brown J, Bray AW, Feeley KP, Sammy MJ, Smith MR, Schurr TG, Vita JA, Ambalavanan N, Calhoun D, Dell’Italia L, Ballinger SW (2016) Endothelial Cell Bioenergetics and Mitochondrial DNA Damage Differ in Humans Having African or West Eurasian Maternal Ancestry. Circulation Cardiovascular genetics 9: 26–36
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiY2lyY2N2ZyI7czo1OiJyZXNpZCI7czo2OiI5LzEvMjYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMi8xMC8yMDIyLjAyLjA5LjIyMjcwNzIzLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

60. [60].Kurl S, Laukkanen JA, Rauramaa R, Lakka TA, Sivenius J, Salonen JT (2003) Cardiorespiratory Fitness and the Risk for Stroke in Men. Archives of Internal Medicine 163: 1682
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archinte.163.14.1682&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12885683&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000184369900007&link_type=ISI) 

61. [61].Laufs U, Wassmann S, Czech T, Münzel T, Eisenhauer M, Böhm M, Nickenig G (2005) Physical inactivity increases oxidative stress, endothelial dysfunction, and atherosclerosis. Arteriosclerosis thrombosis and vascular biology 25: 809–814
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYXR2YmFoYSI7czo1OiJyZXNpZCI7czo4OiIyNS80LzgwOSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzEwLzIwMjIuMDIuMDkuMjIyNzA3MjMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

62. [62].Lee S, Fuchsberger C, Kim S, Scott L (2016) An efficient resampling method for calibrating single and gene-based rare variant association analysis in case-control studies. Biostatistics Oxford England 17: 1–15
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/biostatistics/kxv033&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26363037&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

63. [63].Lempiäinen H, Brænne I, Michoel T, Tragante V, Vilne B, Webb TR, Kyriakou T, Eichner J, Zeng L, Willenborg C, Franzen O, Ruusalepp A, Goel A, van der Laan SW, Biegert C, Hamby S, Talukdar HA, Foroughi Asl H, CVgenes@target consortium C, Pasterkamp G, et al. (2018) Network analysis of coronary artery disease risk genes elucidates disease mechanisms and druggable targets. Scientific reports 8: 3434
    
    

64. [64].Liu X, Longchamps RJ, Wiggins KL, Raffield LM, Bielak LF, Zhao W, Pitsillides A, Blackwell TW, Yao J, Guo X, Kurniansyah N, Thyagarajan B, Pankratz N, Rich SS, Taylor KD, Peyser PA, Heckbert SR, Seshadri S, Cupples LA, Boerwinkle E, et al. (2021) Association of mitochondrial DNA copy number with cardiometabolic diseases. Cell Genomics 1: 100006
    
    

65. [65].Lopez AD, Mathers CD, Ezzati M, Jamison DT, Murray CJL (2006) Global and regional burden of disease and risk factors, 2001: systematic analysis of population health data. Lancet London England 367: 1747–1757
    
    

66. [66].Ludwig-Slomczynska AH, Seweryn MT, Kapusta P, Pitera E, Handelman SK, Mantaj U, Cyganek K, Gutaj P, Dobrucka Å, Wender-Ozegowska E, Malecki MT, Wolkow PP (2020) Mitochondrial GWAS and association of nuclear - mitochondrial epistasis with BMI in T1DM patients. BMC medical genomics 13: 97
    
    

67. [67].Luo J, Noordam R, Jukema JW, van Dijk KW, Hägg S, Grassmann F, le Cessie S, van Heemst D (2021) Low mitochondrial copy number drives atherogenic cardiovascular disease: evidence from prospective cohort analyses in the UK Biobank combined with Mendelian Randomization
    
    

68. [68].Madamanchi NR, Runge MS (2007) Mitochondrial dysfunction in atherosclerosis. Circ Res 100: 460–473
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNpcmNyZXNhaGEiO3M6NToicmVzaWQiO3M6OToiMTAwLzQvNDYwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDIvMTAvMjAyMi4wMi4wOS4yMjI3MDcyMy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

69. [69].Madamanchi NR, Vendrov A, Runge MS (2005) Oxidative Stress and Vascular Disease. Arteriosclerosis Thrombosis and Vascular Biology 25: 29–38
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYXR2YmFoYSI7czo1OiJyZXNpZCI7czo3OiIyNS8xLzI5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDIvMTAvMjAyMi4wMi4wOS4yMjI3MDcyMy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

70. [70].Manson JE, Skerrett PJ, Greenland P, VanItallie TB (2004) The escalating pandemics of obesity and sedentary lifestyle. A call to action for clinicians. Archives of internal medicine 164: 249–258
    
    

71. [71].Maouche S, Schunkert H (2012) Strategies beyond genome-wide association studies for atherosclerosis. Arterioscler Thromb Vasc Biol 32: 170–181
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYXR2YmFoYSI7czo1OiJyZXNpZCI7czo4OiIzMi8yLzE3MCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzEwLzIwMjIuMDIuMDkuMjIyNzA3MjMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

72. [72].Martínez-Redondo D, Marcuello A, Casajús JA, Ara I, Dahmani Y, Montoya J, Ruiz-Pesini E, López-Pérez MJ, Díez-Sánchez C (2010) Human mitochondrial haplogroup H: the highest VO2max consumer–is it a paradox? Mitochondrion 10: 102–107
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.mito.2009.11.005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19900587&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000274989300003&link_type=ISI) 

73. [73].McRae AF, Byrne EM, Zhao ZZ, Montgomery GW, Visscher PM (2008) Power and SNP tagging in whole mitochondrial genome association studies. Genome Res 18: 911–917
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjg6IjE4LzYvOTExIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDIvMTAvMjAyMi4wMi4wOS4yMjI3MDcyMy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

74. [74].Moreira JD, Gopal DM, Kotton DN, Fetterman JL (2021) Gaining Insight into Mitochondrial Genetic Variation and Downstream Pathophysiology: What Can i(PSCs) Do? Genes 12
    
    

75. [75].Nandakumar P, Tian C, O’Connell J, 23andMe Research Team, Hinds D, Paterson AD, Sondheimer N (2021) Nuclear genome-wide associations with mitochondrial heteroplasmy. Science advances 7
    
    

76. [76].Neiburga K, Vilne B, Bauer S, Bongiovanni D, Ziegler T, Lachmann M, Wengert S, Hawe J, Güldener U, Westerlund A, Li L, Pang S, Yang C, Saar K, Huebner N, Maegdefessel L, Consortium DB, Lange R, Krane M, Schunkert H, et al. (2021) Vascular Tissue Specific miRNA Profiles Reveal Novel Correlations with Risk Factors in Coronary Artery Disease. Biomolecules 11: 1683
    
    

77. [77].Nelson CP, Goel A, Butterworth AS, Kanoni S, Webb TR, Marouli E, Zeng L, Ntalla I, Lai FY, Hopewell JC, Giannakopoulou O, Jiang T, Hamby SE, Di Angelantonio E, Assimes TL, Bottinger EP, Chambers JC, Clarke R, Palmer CNA, Cubbon RM, et al. (2017) Association analyses based on false discovery rate implicate new loci for coronary artery disease. Nature genetics 49: 1385–1391
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3913&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28714975&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

78. [78].Nelson CP, Hamby SE, Saleheen D, Hopewell JC, Zeng L, Assimes TL, Kanoni S, Willenborg C, Burgess S, Amouyel P, Anand S, Blankenberg S, Boehm BO, Clarke RJ, Collins R, Dedoussis G, Farrall M, Franks PW, Groop L, Hall AS, et al. (2015) Genetically determined height and coronary artery disease. The New England journal of medicine 372: 1608–1618
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1404881&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25853659&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

79. [79].Nikpay M, Goel A, Won HH, Hall LM, Willenborg C, Kanoni S, Saleheen D, Kyriakou T, Nelson CP, Hopewell JC, Webb TR, Zeng L, Dehghan A, Alver M, Armasu SM, Auro K, Bjonnes A, Chasman DI, Chen S, Ford I, et al. (2015) A comprehensive 1,000 Genomes-based genome-wide association meta-analysis of coronary artery disease. Nature genetics 47: 1121–1130
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3396&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26343387&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

80. [80].Rea IM, McNerlan SE, Archbold GP, Middleton D, Curran MD, Young IS, Ross OA (2013) Mitochondrial J haplogroup is associated with lower blood pressure and anti-oxidant status: findings in octo/nonagenarians from the BELFAST Study. Age Dordrecht Netherlands 35: 1445–1456
    
    

81. [81].Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, Grody WW, Hegde M, Lyon E, Spector E, Voelkerding K, Rehm HL, ACMGLQAC (2015) Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med 17: 405–424
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/gim.2015.30&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25741868&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

82. [82].Roberts CK, Vaziri ND, Barnard RJ (2002) Effect of diet and exercise intervention on blood pressure, insulin, oxidative stress, and nitric oxide availability. Circulation 106: 2530– 2532
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjExOiIxMDYvMjAvMjUzMCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAyLzEwLzIwMjIuMDIuMDkuMjIyNzA3MjMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

83. [83].Rollins B, Martin MV, Sequeira PA, Moon EA, Morgan LZ, Watson SJ, Schatzberg A, Akil H, Myers RM, Jones EG, Wallace DC, Bunney WE, Vawter MP (2009) Mitochondrial variants in schizophrenia, bipolar disorder, and major depressive disorder. PloS one 4: e4913
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0004913&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19290059&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

84. [84].Ruiz-Pesini E, Lapeña AC, Díez-Sánchez C, Pérez-Martos A, Montoya J, Alvarez E, Díaz M, Urriés A, Montoro L, López-Pérez MJ, Enríquez JA (2000) Human mtDNA haplogroups associated with high or reduced spermatozoa motility. Am J Hum Genet 67: 682–696
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/303040&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10936107&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000089002300015&link_type=ISI) 

85. [85].Samani NJ, Erdmann J, Hall AS, Hengstenberg C, Mangino M, Mayer B, Dixon RJ, Meitinger T, Braund P, Wichmann HE, Barrett JH, König IR, Stevens SE, Szymczak S, Tregouet DA, Iles MM, Pahlke F, Pollard H, Lieb W, Cambien F, et al. (2007) Genomewide association analysis of coronary artery disease. N Engl J Med 357: 443–453
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa072366&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17634449&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248436300005&link_type=ISI) 

86. [86].Samuels DC, Carothers AD, Horton R, Chinnery PF (2006) The power to detect disease associations with mitochondrial DNA haplogroups. Am J Hum Genet 78: 713–720
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/502682&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16532401&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000236755900017&link_type=ISI) 

87. [87].Sayols-Baixeras S, Lluís-Ganella C, Lucas G, Elosua R (2014) Pathogenesis of coronary artery disease: focus on genetic risk factors and identification of genetic variants. Appl Clin Genet 7: 15–32
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24520200&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

88. [88].Schunkert H, König IR, Kathiresan S, Reilly MP, Assimes TL, Holm H, Preuss M, Stewart AFR, Barbalic M, Gieger C, Absher D, Aherrahrou Z, Allayee H, Altshuler D, Anand SS, Andersen K, Anderson JL, Ardissino D, Ball SG, Balmforth AJ, et al. (2011) Large-scale association analysis identifies 13 new susceptibility loci for coronary artery disease. Nat Genet 43: 333–338
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.784&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21378990&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

89. [89].Schunkert H, von Scheidt M, Kessler T, Stiller B, Zeng L, Vilne B (2018) Genetics of coronary artery disease in the light of genome-wide association studies. Clinical research in cardiology official journal of the German Cardiac Society 107: 2–9
    
    

90. [90].Selivanov VA, Zeak JA, Roca J, Cascante M, Trucco M, Votyakova TV (2008) The role of external and matrix pH in mitochondrial reactive oxygen species generation. The Journal of biological chemistry 283: 29292–29300
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamJjIjtzOjU6InJlc2lkIjtzOjEyOiIyODMvNDMvMjkyOTIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMi8xMC8yMDIyLjAyLjA5LjIyMjcwNzIzLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

91. [91].Sobenin IA, Sazonova MA, Postnov AY, Bobryshev YV, Orekhov AN (2013) Changes of mitochondria in atherosclerosis: possible determinant in the pathogenesis of the disease. Atherosclerosis 227: 283–288
    
    

92. [92].Stefano GB, Kream RM (2016) Mitochondrial DNA heteroplasmy in human health and disease. Biomedical Reports 4: 259–262
    
    

93. [93].Sui X, LaMonte MJ, Blair SN (2007) Cardiorespiratory Fitness as a Predictor of Nonfatal Cardiovascular Events in Asymptomatic Women and Men. American Journal of Epidemiology 165: 1413–1423
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwm031&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17406007&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000247240700011&link_type=ISI) 

94. [94].Talukdar HA, Asl HF, Jain RK, Ermel R, Ruusalepp A, Franzén O, Kidd BA, Readhead B, Giannarelli C, Kovacic JC, et al. (2016) Cross-Tissue Regulatory Gene Networks in Coronary Artery Disease. Cell Systems 2: 196–208
    
    

95. [95].Torroni A, Petrozzi M, D’Urbano L, Sellitto D, Zeviani M, Carrara F, Carducci C, Leuzzi V, Carelli V, Barboni P, Negri AD, Scozzari R (1997) Haplotype and phylogenetic analyses suggest that one European-specific mtDNA background plays a role in the expression of Leber hereditary optic neuropathy by increasing the penetrance of the primary mutations 11778 and 14484. Am J Hum Genet 60: 1107–1121
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9150158&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997WX76000011&link_type=ISI) 

96. [96].Trégouët DA, König IR, Erdmann J, Munteanu A, Braund PS, Hall AS, Grosshennig A, Linsel-Nitschke P, Perret C, DeSuremain M, Meitinger T, Wright BJ, Preuss M, Balmforth AJ, Ball SG, Meisinger C, Germain C, Evans A, Arveiler D, Luc G, et al. (2009) Genome-wide haplotype association study identifies the SLC22A3-LPAL2-LPA gene cluster as a risk locus for coronary artery disease. Nature genetics 41: 283–285
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.314&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19198611&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000263640200007&link_type=ISI) 

97. [97].van der Harst P, Verweij N (2018) Identification of 64 Novel Genetic Loci Provides an Expanded View on the Genetic Architecture of Coronary Artery Disease. Circulation research 122: 433–443
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNpcmNyZXNhaGEiO3M6NToicmVzaWQiO3M6OToiMTIyLzMvNDMzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDIvMTAvMjAyMi4wMi4wOS4yMjI3MDcyMy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

98. [98].van Oven M (2015) PhyloTree Build 17: Growing the human mitochondrial DNA tree. Forensic Science International Genetics Supplement Series 5: e392–e394
    
    

99. [99].Vansteelandt S, Goetgeluk S, Lutz S, Waldman I, Lyon H, Schadt EE, Weiss ST, Lange C (2009) On the adjustment for covariates in genetic association analysis: a novel, simple principle to infer direct causal effects. Genetic epidemiology 33: 394–405
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/gepi.20393&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19219893&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000267605700004&link_type=ISI) 

100.[100].Vellers HL, Verhein KC, Burkholder AB, Lee J, Kim Y, Lightfoot JT, Shi M, Weinberg CR, Sarzynski MA, Bouchard C, Kleeberger SR (2020) Association between Mitochondrial DNA Sequence Variants and V·O2 max Trainability. Medicine Science in Sports Exercise 52: 2303–2309
    
    

101.[101].Vilne B, Schunkert H (2018) Integrating Genes Affecting Coronary Artery Disease in Functional Networks by Multi-OMICs Approach. Frontiers in cardiovascular medicine 5: 89
    
    

102.[102].Vilne B, Skogsberg J, Foroughi Asl H, Talukdar HA, Kessler T, Björkegren JLM, Schunkert H (2017) Network analysis reveals a causal role of mitochondrial gene activity in atherosclerotic lesion formation. Atherosclerosis 267: 39–48
    
    

103.[103].Wallace DC (2015) Mitochondrial DNA variation in human radiation and disease. Cell 163: 33–38
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2015.08.067&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26406369&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

104.[104].Wallace DC (2018) Mitochondrial genetic medicine. Nature Genetics 50: 1642–1649
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0264-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30374071&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 

105.[105].Wallace DC, Chalkia D (2013) Mitochondrial DNA genetics and the heteroplasmy conundrum in evolution and disease. Cold Spring Harb Perspect Biol 5: a021220
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImNzaHBlcnNwZWN0IjtzOjU6InJlc2lkIjtzOjEyOiI1LzExL2EwMjEyMjAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMi8xMC8yMDIyLjAyLjA5LjIyMjcwNzIzLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

106.[106].Webb TR, Erdmann J, Stirrups KE, Stitziel NO, Masca NGD, Jansen H, Kanoni S, Nelson CP, Ferrario PG, König IR, Eicher JD, Johnson AD, Hamby SE, Betsholtz C, Ruusalepp A, Franzén O, Schadt EE, Björkegren JLM, Weeke PE, Auer PL, et al. (2017) Systematic Evaluation of Pleiotropy Identifies 6 Further Loci Associated with Coronary Artery Disease. Journal of the American College of Cardiology 69: 823–836
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czo0OiJhY2NqIjtzOjU6InJlc2lkIjtzOjg6IjY5LzcvODIzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDIvMTAvMjAyMi4wMi4wOS4yMjI3MDcyMy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

107.[107].Wessel TR (2004) Relationship of Physical Fitness vs Body Mass Index with Coronary Artery Disease and Cardiovascular Events in Women. JAMA 292: 1179
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.292.10.1179&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15353530&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000223723800017&link_type=ISI) 

108.[108].Yonova-Doing E, Calabrese C, Gomez-Duran A, Schon K, Wei W, Karthikeyan S, Chinnery PF, Howson JMM (2021) An atlas of mitochondrial DNA genotype-phenotype associations in the UK Biobank. Nature genetics 53: 982–993
    
    

109.[109].Zhao Z, Yang Y, Zeng Y, He M (2016) A microfluidic ExoSearch chip for multiplexed exosome detection towards blood-based ovarian cancer diagnosis. Lab on a chip 16: 489– 496
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1039/C5LC01117E&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26645590&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F02%2F10%2F2022.02.09.22270723.atom)