Article Text

PDF

Mitochondrial haplogroups associated with elite Japanese athlete status
  1. Eri Mikami1,2,
  2. Noriyuki Fuku2,
  3. Hideyuki Takahashi3,
  4. Nao Ohiwa3,
  5. Robert A Scott4,
  6. Yannis P Pitsiladis4,
  7. Mitsuru Higuchi5,
  8. Takashi Kawahara3,
  9. Masashi Tanaka2
  1. 1Graduate School of Sport Sciences, Waseda University, Saitama, Japan
  2. 2Department of Genomics for Longevity and Health, Tokyo Metropolitan Institute of Gerontology, Tokyo, Japan
  3. 3Japan Institute of Sports Sciences, Tokyo, Japan
  4. 4Integrative and Systems Biology, Faculty of Biomedical and Life Sciences (FBLS),University of Glasgow, Glasgow, UK
  5. 5Waseda University, Saitama, Japan
  1. Correspondence to Dr Noriyuki Fuku, Department of Genomics for Longevity and Health, Tokyo Metropolitan Institute of Gerontology, 35-2 Sakae-cho, Itabashi-ku, Tokyo 173-0015, Japan; nfuku{at}tmig.or.jp

Abstract

Purpose It has been hypothesised that certain mitochondrial haplogroups, which are defined by the presence of a characteristic cluster of tightly linked mitochondrial DNA polymorphisms, would be associated with elite Japanese athlete status. To examine this hypothesis, the frequencies of mitochondrial haplogroups found in elite Japanese athletes were compared with those in the general Japanese population.

Methods Subjects comprised 139 Olympic athletes (79 endurance/middle-power athletes (EMA), 60 sprint/power athletes (SPA)) and 672 controls (CON). Two mitochondrial DNA fragments containing the hypervariable sequence I (m16024–m16383) of the major non-coding region and the polymorphic site at m.5178C>A within the NADH dehydrogenase subunit 2 gene were sequenced, and subjects were classified into 12 major mitochondrial haplogroups (ie, F, B, A, N9a, N9b, M7a, M7b, M*, G2, G1, D5 or D4). The mitochondrial haplogroup frequency differences among EMA, SPA and CON were then examined.

Results EMA showed an excess of haplogroup G1 (OR 2.52, 95% CI 1.05 to 6.02, p=0.032), with 8.9% compared with 3.7% in CON, whereas SPA displayed a greater proportion of haplogroup F (OR 2.79, 95% CI 1.28 to 6.07, p=0.007), with 15.0% compared with 6.0% in CON.

Conclusions The results suggest that mitochondrial haplogroups G1 and F are associated with elite EMA and SPA status in Japanese athletes, respectively.

Statistics from Altmetric.com

Elite athletic performance is a complex, multifactorial phenotype. A number of familial studies have assessed the relative contribution of genetic and environmental factors to physical performance-related traits and estimated that there is a significant genetic component to phenotypes such as maximal oxygen uptake (V.O2max)1 2 and muscle strength.3 4 Therefore, it is commonly believed that genetic factors are likely determinants of elite athlete status, which may require a synergy of advantageous physical performance phenotypes. Over 200 genes, in both the nuclear and mitochondrial genomes, have recently been suggested to have an effect on physical performance and health-related fitness.5

Mitochondria are essential to all higher organisms for sustaining life, and are extremely important in energy metabolism, providing 36 molecules of ATP per glucose molecule in contrast to the two ATP molecules produced by glycolysis. Although most DNA is packaged in chromosomes within the nucleus, mitochondria also possess their own circular DNA: mitochondrial DNA (mtDNA). The 16 569-bp human mtDNA contains 13 genes for mitochondrial oxidative phosphorylation (OXPHOS), as well as two ribosomal RNA and 22 transfer RNA genes that are necessary for protein synthesis within mitochondria.6 Unlike nuclear DNA, mtDNA is inherited maternally. Interestingly, in familial studies, aerobic capacity has been found to have stronger maternal inheritance than paternal.7 8 It has also been reported that mtDNA polymorphisms can influence the interindividual variations in aerobic capacity and its trainability.9 10 Furthermore, patients with mutations in mtDNA commonly present with exercise intolerance, muscle weakness and increased production of lactate.11 It is clear, therefore, that mtDNA polymorphisms represent a promising candidate to contain variants influencing physical performance.

The matrilineal inheritance of mtDNA and linear accumulation of polymorphisms has allowed the construction of detailed mtDNA phylogenies.12 These phylogenies display the variation and diversity of human mtDNA and allow haplogroup identification through the analysis of a small number of haplogroup-specific polymorphisms. This matrilineal pattern of descent means that individual haplotypes share linked complexes of polymorphisms common to all sequences in a haplogroup. We previously reported that elite Kenyan endurance athletes differed in their mitochondrial haplogroup distribution compared with the general Kenyan population,13 although this was not the case in elite Ethiopian endurance athletes.14 Mitochondrial haplogroups have previously been associated with the elite endurance athlete status15 16 and V.O2max17 in Europeans. Mitochondrial haplogroup distributions display geographical diversity, and certain haplogroups can be regionally specific. Indeed, almost all mitochondrial haplogroups in Africans and Europeans are not present in Asians, including Japanese.

There are highly diverse mtDNA sequences in Africans because of a time depth for mtDNA lineages of approximately 200 000 years.18 Although mtDNA lineages in east Asians are shorter than in Africans (approximately 30 000–55 000 years), they also have a relatively long history. Japanese people are descendents of immigrants from the northern route, who had adapted to cold climates and/or famine, and those from the southern route.19 Therefore, acquired mtDNA variations may have permitted the founding populations of modern Japanese individuals to adapt to extreme environmental conditions such as cold climates and/or famine, and consequently influencing modern human phenotypes related to altered bioenergetics or mitochondrial function.20 Indeed, we previously reported that certain mitochondrial haplogroups were associated with metabolic disorders closely related to mitochondrial function, such as obesity,21 type 2 diabetes mellitus,22,,24 metabolic syndrome22 25 and myocardial infarction26 in Japanese individuals. Therefore, here we hypothesise that mitochondrial genetic variation in Japanese people may also associate with exercise capacity, because certain mitochondrial haplogroups are related not only to mitochondrial dysfunction, but also to enhanced mitochondrial function. Given elite athletes are a select population with greater exercise capacity, it is of interest to study the genetic determinants of their superior exercise performance. In the present study, we assessed the mitochondrial haplogroup frequencies of elite Japanese athletes compared with the general Japanese population in order to determine the influence of mitochondrial haplogroups on elite Japanese athlete status.

Methods

Subjects

The subjects consisted of 141 elite Japanese athletes (110 men and 31 women). All athletes had represented Japan at the Olympic Games. A total of 141 athletes was classified as endurance/middle-power athletes (EMA, 81) and sprint/power athletes (SPA, 60) based on the criteria of Yang et al.27 The EMA group included 13 endurance runners competing in events of 800 m or more, 10 sailing athletes, seven swimmers competing in events of 200 m or more, seven rowers, five long-distance cyclists, seven canoeists, nine volleyball players, six basketball players, six hockey players, four soccer players, four water polo players, two boxers and one modern pentathlete. The SPA comprised 18 track and field athletes (seven sprinters competing in events of ≤400 m, six jumpers and five throwers), nine swimmers competing in events of 100 m or less, seven gymnasts, seven competitive fencers, six divers, five wrestlers, four weightlifters, two short-distance track cyclists and two judo athletes. Mitochondrial haplogroup frequencies determined in 672 general Japanese populations (387 men and 285 women) from our Human Mitochondrial Genome Single Nucleotide Polymorphism Database (http://mtsnp.tmig.or.jp/mtsnp/index_e.shtml) represented the control group (CON).19 Written consent was obtained from each subject, and the study was approved by the Ethics Committees of the Japan Institute of Sports Sciences and Tokyo Metropolitan Institute of Gerontology.

Data collection and analysis

Total DNA was isolated from venous blood by the use of QIAamp DNA Blood Maxi Kit (QIAGEN, Hilden, Germany). Two fragments of mtDNA including the hypervariable sequence I (m.16024–m.16383) and nucleotide position m.5178 were amplified by PCR as detailed in table 1. The forward primer was a 38-mer oligonucleotide, consisting of an 18-base sequence of a universal forward sequencing primer (−21M13: 5′–TGTAAAACGACGGCCAGT–3′) connected on its 5′ side to the 3′ side of a 20-base light-strand-specific sequence (table 1). PCR amplifications were carried out in a final reaction volume of 10 μl, containing 20 ng of total DNA, 0.25 units of TaKaRa Ex TaqHS (Takara, Shiga, Japan), 1 μl of 10×Ex Taq buffer, 0.2 mM of each dNTP and 0.5 μM of each primer. The PCR conditions used were an initial denaturation step at 94°C for 5 min, followed by 40 cycles of denaturation at 94°C for 15 s, annealing at 60°C for 15 s, and extension at 72°C for 1 min, with a final extension of 10 min at 72°C. Following this, sequence reactions were performed by use of the PCR template, −21M13 forward primer, and a BigDye Terminator Cycle Sequencing Kit version 3.1 (Applied Biosystems, Foster City, California, USA). Sequences were analysed with an automated DNA sequencer Applied Biosystems 3130xl Genetic Analyzer (Applied Biosystems). In some cases, the light-strand sequences could not be determined on stretches of C due to the m.16189T>C transition. In such cases, the heavy-strand (reverse) sequences were also determined. For this purpose, the primer H81 (table 1) was used. Each of the mtDNA sequences was compared with the original6 and the revised28 Cambridge reference sequences, and polymorphisms were confirmed visually by the use of DNA Sequencing Software Sequencher version 4.2.2 (Gene Codes Corporation, Ann Arbor, Michigan, USA).

Table 1

Primers used for analysis

Haplogroup classification

Subjects were classified into 12 major mitochondrial haplogroups (ie, F, B, A, N9a, N9b, M7a, M7b, M*, G2, G1, D5 or D4) on the basis of the presence of hypervariable sequence I polymorphisms and several protein-coding-region polymorphisms (table 2).19 Data from two of 141 individuals were deemed to be unusable because of incomplete haplogroup classification.

Table 2

Polymorphisms used to determine mitochondrial haplogroups

The mitochondrial haplogroups were determined by the combination of both the main polymorphism and one of the secondary polymorphisms in the mitochondrial DNA. Four protein-coding-region polymorphisms (m.15851A>G [Cytb], m.15860A>G [Cytb], m.15874A>G [Cytb] and m.5178C>A [ND2]) were used to determine haplogroups B, G1, D4 and D5, respectively (Cytb, Cytochrome b; ND2, NADH dehydrogenase subunit 2).

Statistical analysis

Mitochondrial haplogroup frequencies between EMA and CON, and between SPA and CON were compared by a χ2 test by the use of JMP version 8 (SAS Institute Japan, Tokyo, Japan). For haplogroup analysis, we compared each haplogroup versus the sum of all other haplogroups. The p value, OR and 95% CI were calculated. A p value of less than 0.05 was considered statistically significant.

Results

The distribution of mitochondrial haplogroups in EMA and SPA compared with CON is shown in figure 1. When the frequency of each haplogroup versus the sum of all others was compared between EMA or SPA and CON, it was found that EMA displayed an excess of haplogroup G1 (OR 2.52, 95% CI 1.05 to 6.02, p=0.032; table 3), with 8.9% in EMA compared with 3.7% in CON. In addition, the frequency of haplogroup B in EMA tended to be lower than in CON (OR 0.45, 95% CI 0.18 to 1.14, p=0.084; table 3), although this difference was not statistically significant. On the other hand, SPA showed a greater proportion of haplogroup F (OR 2.79, 95% CI 1.28 to 6.07, p=0.007; table 3), with 15.0% in SPA compared with 6.0% in CON. The overall frequency of mitochondrial haplogroups F and G1 in our control samples was 7.0% and 3.5% (394 and 198 of 5651); 6.7% and 3.0% in Gunma Prefecture (95 and 43 of 1418) (unpublished data), 7.6% and 3.8% in Gifu Prefecture (208 and 104 of 2748) (unpublished data) and 6.1% and 3.4% in Tokyo (91 and 51 of 1493) (unpublished data). Performance-associated haplogroups F and G1 thus have a wide geographical distribution in Japan.

Figure 1

Mitochondrial haplogroup distribution in endurance athletes, sprint/power athletes and controls. Significant differences from controls (grey bars) are highlighted by asterisks (*p<0.05, **p<0.01). Endurance athletes displayed an excess of the haplogroup G1 compared with controls. Sprint/power athletes displayed an excess of haplogroup F compared with controls.

Table 3

Mitochondrial haplogroup distribution of subject groups

Discussion

We found that the frequencies of mitochondrial haplogroups differed significantly between elite Japanese athletes and the general Japanese population. EMA displayed an excess of haplogroup G1, whereas SPA showed an excess of haplogroup F. This is the first study to report that Asian-specific mitochondrial haplogroups G1 and F are associated with elite athlete status. We recently reported that mitochondrial haplogroups were associated with elite Kenyan athlete status,13 but not with elite Ethiopian14 and Jamaican/African-American athlete status (Deason et al, 2010, manuscript in preparation). Elite Kenyan athletes displayed an excess of haplogroup L0 and a dearth of haplogroup L3 compared with the general Kenyan population. Niemi et al16 also found differences in mitochondrial haplogroup frequencies between elite Finnish sprint and endurance athletes; no endurance athletes belonged to mitochondrial haplogroup K and subhaplogroup J2. Furthermore, Castro et al15 reported that haplogroup T was negatively associated with elite Spanish endurance athlete status. These associations with elite athlete status may suggest that these haplogroups contain polymorphisms that influence some aspect of exercise performance or its trainability.

The greater diversity of mtDNA sequences among African population is due to the fact that the most recent ancestor of modern humans (so called ‘mitochondrial Eve’) originated from Africa.29 30 As such, each of the African haplogroups (L0–L3) has deep genetic roots. Haplogroup L3 is proposed to be the ancestor of all non-African populations. European haplogroups (H, I, J, K, S, T, U, V, W, etc) belong to macrohaplogroup N,31 whereas Asian haplogroups belong to both macrohaplogroups N and M (haplogroups A, B, F and N9 to macrohaploup N; and haplogroups M7a, M7b, M8, D and G to macrohaploup M).19 Both macrohaplogroups N and M have a common root with haplogroup L3.32 Africans, Europeans and Asians thus have different phylogenetic structures. Therefore, the mtDNA polymorphisms of these ethnic groups and their effects on exercise performance must be analysed according to the haplogroup frequencies in each population.

Mitochondrial haplogroup G1 is characterised by three polymorphisms: m.8200T>C, m.15323G>A and m.15497G>A, in the protein-coding region of the mtDNA.19 32 Two of these, namely, m.15323G>A and m.15497G>A, are non-synonymous substitutions. These polymorphisms cause the Ala193Thr and Gly251Ser replacements in Cytb, which is a subunit of the complex III. Previously, the m.15498G>A mutation causing the Gly251Asp replacement in the Cytb was reported in a patient with mitochondrial myopathy and exercise intolerance.33 Interestingly, in middle-aged individuals, we reported an association of obesity-related phenotypes with m.15497G>A transition (Cytb: Gly251Ser), which is a polymorphism characterising haplogroup G1.34 It is reasonable to speculate that this Gly251Ser replacement is accompanied by functional alterations of Cytb.

We previously hypothesised two possibilities to explain the association between the obesity-related phenotype and the m.15497G>A transition characterising mitochondrial haplogroup G1: (1) increased efficiency of mitochondrial energy conservation at the cytochrome bc1 complex resulting in decreased energy consumption; or (2) inhibiting reduction in ubiquinone at the Qo site (one of the ubiquinone-binding sites of the complex III) resulting in a reduced β-oxidation of fatty acid, which leads to fat accumulation.34 In the present study, haplogroup G1 was associated with elite Japanese EMA status. The main function of mitochondria is to produce ATP by OXPHOS; and while the uncoupling of mitochondrial OXPHOS generates heat, it concomitantly reduces the production of ATP.35 Conversely, more tightly coupled OXPHOS would be expected to decrease heat production and result in higher efficiency of ATP production. This improved efficiency of ATP production could explain, at least partly, the association between haplogroup G1 and endurance/middle-power performance reported in the present study. However, this energy conservation by mitochondrial OXPHOS may predispose to obesity in sedentary individuals later on in life; a phenomenon commonly referred to as the ‘thrifty’ genotype and/or ‘thrifty’ phenotype.20 These hypotheses require further investigation.

In the present study, we found that the frequency of mitochondrial haplogroup F was significantly higher in SPA than in CON. This haplogroup is defined by four polymorphisms: m.3970C>T, m.6392T>C, m.10310G>A and m.13928G>C, in the coding region of the mtDNA.19 32 One of the polymorphisms, namely m.13928G>C, causes Ser531Thr replacement in the ND5, which is a subunit of complex I. Sprint performance relies more on anaerobic glycolysis than OXPHOS. Certain mtDNA polymorphisms and/or mitochondrial haplogroups may influence the regulation of ATP production not only by the OXPHOS system in the mitochondria but also by the glycolytic pathway in the cytosol. We reported that the peak cytosolic calcium levels after histamine stimulation were higher in the cybrids with mitochondrial macrohaplogroup N than in those with mitochondrial macrohaplogroup M.36 This result suggests that the cytosolic calcium response may be enhanced in the cybrids with macrohaplogroup N. Sprint/power performance-related mitochondrial haplogroup F is a component of macrohaplogroup N. Therefore, haplogroup F may be related to the calcium dynamics in the cell including skeletal muscle. Calcium regulates glycogen breakdown in the skeletal muscle. During activation of contraction in skeletal muscle, calcium is released from the sarcoplasmic reticulum, and also activates phosphorylase kinase. The physiological significance of this calcium activation process is that muscle contraction is triggered by a transient increase in the level of cytosolic calcium through its release from intracellular reservoirs, which may also include mitochondrial calcium, by nerve impulses. The rate of glycogen breakdown is linked to the rate of muscle contraction, an important regulatory link because glycogen breakdown in muscle provides fuel for glycolysis, which, in turn, generates the ATP required for muscle contraction. Mitochondrial haplogroup F, which is associated with SPA status in the present study, may thus influence the rate of glycolytic ATP production and/or the rate of muscle contraction. These hypotheses require further investigation by functional analysis of this haplogroup.

In conclusion, we found significant associations both between mitochondrial haplogroup G1 and elite EMA status and between mitochondrial haplogroup F and elite SPA status in Japanese athletes. These associations may implicate these mitochondrial haplogroups in determining elite athlete status in Japanese individuals (as previously suggested in other populations). It should be noted that the sample size of the present study was relatively small but in line with other similar studies.15 16 Nevertheless, our previous findings suggest physiological plausibility.34 Further investigation will require more detailed analysis of mtDNA.

What is already known on this topic

Mitochondrial haplogroups, which are a set of mitochondrial DNA polymorphisms, appear to influence physical performance. Indeed, previous studies have reported associations between certain mitochondrial haplogroups and elite endurance athlete status or V.O2max in African and European, but not in Asian individuals.

What this study adds

In the present study, Asian-specific mitochondrial haplogroup G1 was found to associate with elite Japanese EMA status, whereas haplogroup F was related to elite Japanese SPA status. Mitochondria may therefore regulate not only aerobic metabolism but also anaerobic metabolism.

References

View Abstract

Footnotes

  • Funding This work was supported in part by a Grant-in-Aid for Young Scientists (A-21680050 and B-18700541 to NF), a Grant-in-Aid for Exploratory Research (20650113 to NF), a Grant-in-Aid for Scientific Research (A-15200051 to MT), and a Grant-in-Aid for the Global COE (Sport Sciences for the Promotion of Active Life to Waseda University) from the Ministry of Education, Culture, Sports, Science and Technology; by a research project in Japan Institute of Sports Sciences (to TK); and by grants for scientific research from The Uehara Memorial Foundation (to NF) and from The Takeda Science Foundation (to MT).

  • Competing interests None.

  • Ethics approval This study was conducted with the approval of the Japan Institute of Sports Sciences and Tokyo Metropolitan Institute of Gerontology.

  • Patient consent Obtained.

  • Provenance and peer review Not commissioned; externally peer reviewed.

Request permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.