Elite athletes’ genetic predisposition for altered risk of complex metabolic traits

Abstract

Background

Genetic variants may predispose humans to elevated risk of common metabolic morbidities such as obesity and Type 2 Diabetes (T2D). Some of these variants have also been shown to influence elite athletic performance and the response to exercise training. We compared the genotype distribution of five genetic Single Nucleotide Polymorphisms (SNPs) known to be associated with obesity and obesity co-morbidities (IGF2BP2 rs4402960, LPL rs320, LPL rs328, KCJN rs5219, and MTHFR rs1801133) between athletes (all male, n = 461; endurance athletes n = 254, sprint/power athletes n = 207), and controls (all male, n = 544) in Polish and Russian samples. We also examined the association between these SNPs and the athletes’ competition level (‘elite’ and ‘national’ level). Genotypes were analysed by Single-Base Extension and Real-Time PCR. Multinomial logistic regression analyses were conducted to assess the association between genotypes and athletic status/competition level.

Results

IGF2BP2 rs4402960 and LPL rs320 were significantly associated with athletic status; sprint/power athletes were twice more likely to have the IGF2BP2 rs4402960 risk (T) allele compared to endurance athletes (OR = 2.11, 95% CI = 1.03-4.30, P <0.041), and non-athletic controls were significantly less likely to have the T allele compared to sprint/power athletes (OR = 0.62, 95% CI =0.43-0.89, P <0.0009). The control group was significantly more likely to have the LPL rs320 risk (G) allele compared to endurance athletes (OR = 1.26, 95% CI = 1.05-1.52, P <0.013). Hence, endurance athletes were the “protected” group being significantly (p < 0.05) less likely to have the risk allele compared to sprint/power athletes (IGF2BP2 rs4402960) and significantly (p < 0.05) less likely to have the risk allele compared to controls (LPL rs320). The other 3 SNPs did not show significant differences between the study groups.

Conclusions

Male endurance athletes are less likely to have the metabolic risk alleles of IGF2BP2 rs4402960 and LPL rs320, compared to sprint/power athletes and controls, respectively. These results suggest that some SNPs across the human genome have a dual effect and may predispose endurance athletes to reduced risk of developing metabolic morbidities, whereas sprint/power athletes might be predisposed to elevated risk.

Background

Complex metabolic diseases such as Obesity and Type 2 Diabetes (2TD) and physical activity levels have long been recognised as being closely-related. For instance, it has been shown that elite athletes or former elite athletes tend to have longer life expectancies, and lower risks of complex metabolic diseases such as obesity and T2D, than matched sedentary controls [1-3]. Genetic factors seem to play a role in elite athlete development, on one hand [4,5], and the predisposing for complex metabolic diseases, on the other hand [6]. Recently, we [6] and others [7,8] hypothesised that genetic Single Nucleotide Polymorphisms (SNPs), including SNPs identified in Genome Wide Association Studies (GWAS) that have been associated with increased risk for complex metabolic diseases, would also be candidates to influence athletic performance/physical activity levels.

The A/T polymorphism (rs9939609) in the fat mass and obesity associated (FTO) gene, was discovered in two separate GWAS [9,10], and is an example of specific variant associated with obesity, 2TD, and physical activity levels. Recent meta-analysis combining data from adults and children, and an adolescent population (overall 54 studies of n = 218,166 and n = 19,268, respectively) have shown that physically active people with the FTO risk allele are 30% less likely to be obese compared to their inactive counterparts [11]. Visfatin, a recently discovered adipokine that contributes to glucose and obesity-related conditions, is another gene that potentially influences both exercise-related phenotypes and complex metabolic diseases. rs4730153 within the Visfatin was associated with aerobic exercise training-induced changes in glucose and obesity-related phenotypes [12]. The peroxisome proliferator-activated receptor gamma coactivator1α (PPARGC1A) Gly482Ser SNP was also associated with increased risk of obesity and type 2 diabetes [13] on one hand, and with elite athletic performance [14-17], on the other hand.

The outcomes of the abovementioned studies assist with understanding the genomic link between complex metabolic diseases and athletic performance; however the widely accepted hypothesis is that there are likely to be many other uncovered variants with dual effects. In that sense, elite athletes represent the end point of the human physical activity continuum with a “rare” and distinguished phenotype, and hence are an excellent model to study.

Potential obesity and T2D-related genetic variants that may influence athletic performance as well are located in the IGF2BP2, LPL, KCJN, and MTHFR genes. IGF2BP2 rs4402960 G > T variant is associated with predisposition to T2D and obesity. GWAS studies have indicated that the risk allele for T2D and obesity is the T allele. Animal model and human studies implicate this variant with reduced beta-cell function, insulin secretion and sensitivity and with raised fasting glucose levels [18-20]. Importantly, recent studies suggest a potential role for IGF2BP2 in skeletal muscle cell proliferation and differentiation [21]. LPL rs320 and rs328 SNPs have been associated with plasma lipids levels, through the protein’s role in the uptake of Free Fatty Acids (FFA) from the plasma to tissues, including muscle cells [22-24]. Thus, it has been hypothesised that these SNPs may alter the availability of FFA to muscle cells and to the utilization of fat by muscles. The obesity risk allele/genotype for both rs320 and rs328 are G allele and the GG genotype. KCNJ11 is an ATP-sensitive K+ (KATP) channel, which couples cell metabolism with membrane excitability in various cell types, including muscle cells. The protein’s known function is mainly related to diabetes phenotypes [25]. However, it was also found to be association with impaired exercise stress response in several models. The E23K SNP at codon 23 of the KCNJ11 gene (rs5219) results in substitution of glutamic acid to lysine, and may cause modest reductions in ATP sensitivity, which could influence muscle response to exercise. The metabolic risk allele/genotype in rs5219 is T/TT. MTHFR is a key enzyme in one carbon cycle. MTHFR C677T SNP results in elevated plasma homocysteine, which has been linked to reduced mobility and muscle functioning in the elderly (women) and has been associated with T2D. The risk allele/genotype in rs1801133 is T/TT [26,27].

Therefore, we studied the association between these five genetic variants associated with both obesity and obesity co-morbidities (IGF2BP2 rs4402960, LPL rs320, LPL rs328, KCJN rs5219, and MTHFR rs1801133) and elite athletic status in a relatively-large cohort (n = 929, from Poland and Russia) of sprint/power and endurance athletes. We also examined the association between these variants and athletic status according to the athletes’ level of competition (‘elite’ and ‘national’ level). We hypothesised that the obesity and/or co-morbidities risk allele/genotype in each of these variants would be underrepresented in elite athletes compared to controls.

Methods

The study was approved by the Pomeranian Medical University Ethics Committee, Poland, and the Ural State University of Physical Culture, Russia, and written informed consent was obtained from each participant. The study complied with the guidelines set out in the Declaration of Helsinki and the ethics policy of the Szczecin University [28].

Participants

A total of 929 male participants from Russia (n = 281) and Poland (n = 648) were involved in the study. The Russian participants included 177 athletes (mean age = 26.3, SD = 10.3) and 104 unrelated sedentary controls (mean age = 31.2, SD = 10.4). The Polish participants were 208 athletes (mean age = 28.6, SD = 6.2) and 440 unrelated sedentary controls (mean age = 22.4, SD = 2.5). All athletes were ranked in the top 10 nationally in their sport discipline and grouped as being either ‘elite-level’ or ‘national-level’ based on their best personal performance. Those in the elite group had participated in international competitions such as World and European Championships, and/or Olympic Games, whereas those in the national-level group had participated in national competitions only. Athletes were further classified as endurance (events requiring predominantly aerobic energy production including long distance and duration events or sprint/power athletes (events requiring predominantly anaerobic energy production).

Russian athletes

This group included 70 (elite n = 10; 14%) endurance athletes and 107 (elite n = 44; 41%) sprint/power athletes. Athletes in classified as endurance athletes included cross country skiers (n = 38), marathon runners (n = 2), rowers (n = 13), 1500-5000 m runners (n = 4), 5000/10 000 m long distance skaters (n = 6), 3000 m steeple chase (n = 1), 1500 m swimmers (n = 2) and walkers (n = 4). Athletes classified as sprint/power athletes came from sports including ski-cross/ alpine skiing (n = 2), discus throw (n = 1), Greco-roman wrestling (n = 10), pole vault (n = 1), 200 m sprint (n = 2), 100 m sprint (n = 10), 400 m sprint (n = 1), 500 m short distance skating (n = 23), 50/100 m short distance swimming (n = 8), and power lifters (n = 49).

Polish athletes

This group included 108 (elite n = 65; 60%) endurance athletes and 100 (elite n = 64; 64%) sprint/power athletes. Athletes classified as endurance athletes included canoeists (n = 10), cross country skiers (n = 2), cyclists (n = 14), 10 00 m/marathon runners (n = 25), rowers (n = 41), 1500 m swimmers (n = 10), and triathletes (n = 6). The sprint/power athletes included 100/200 m runners (n = 34), archers (n = 4), weight/ power lifters (n = 42), high jumpers (n = 1), javelin throwers (n = 1), long jumpers (n = 4), vaulters (n = 3), shooters (n = 1), shot putters (n = 5), 500 m short-distance skaters (n = 1) sky jumpers (n = 2) and 50/100 m swimmers (n = 2).

Genotyping

In the Polish cohort, Genomic DNA was isolated from buccal epithelium using GenElute Mammalian Genomic DNA Miniprep Kit (Sigma, Hamburg, Germany) according to the manufacturer’s instructions. In the Russian cohort, Genomic DNA was isolated from buccal epithelium or peripheral blood, during the years 2011-2013, using the Diatom™ DNA Prep kit (Cat. # D 1025, IsoGene Lab Ltd, Russia). Genotyping was performed as previously described [29]. In the Russian cohort, genotyping of four SNPs (IGFBP2 rs4402960, KCJN11 rs5219, LPL rs320 and rs328) was performed by Single-Base Extension (SBE) method. The sequence surrounding each SNP was obtained from the Genome Reference Consortium Human genome build 37 assembly from the Ensembl Project (www.ensembl.org). The Primer3web software v. 4.0.0 (http://bioinfo.ut.ee/primer3) was used for designing the PCR primers. Genotyping of the MTHFR rs1801133 polymorphism was performed by using a TaqMan® SNP Genotyping Assay with a StepOne™ Real-Time PCR System (Applied Biosystems, Foster city, CA, USA). The assay ID was C___1202883_20. The results were analysed by using TaqMan® Genotyper Software (Applied Biosystem). K562 DNA High Molecular Weight from Promega Corp. (Cat # DD2011, Madison, WI, USA) served as a positive control sample. Genetic profile of the K562 DNA was as follow: IGFBP2 rs4402960 – G/G, KCJN11 rs5219 – C/T, LPL rs320 – G/G, LPL rs328 – G/C, and MTHFR rs1801133 – G/G.

In the Polish cohort, all samples were genotyped in duplicate using allelic discrimination assays with Taqman® probes (Applied Biosystems, Carlsbad, California, USA) on a CFX96 Touch™ Real-Time PCR Detection System (Bio-Rad, Hercules, California, USA). To discriminate rs4402960, rs320, rs328, rs5219, and rs1801133 alleles, TaqMan® Pre-Designed SNP Genotyping Assays were used (assay IDs: C___2165199_10, C___1843003_10, C____901792_1_, C__11654065_10, C___1202883_20 respectively), including appropriate primers and fluorescently labeled (FAM and VIC) MGB™ probes to detect the alleles.

Genotyping reliability across two laboratories

As previously described [29] genotyping was performed in duplicate in the same Laboratory for accuracy. Two independent investigators have called the genotyping score in each laboratory-100% of the genotypes could be called. For the purpose of results reliability across two laboratories in two different countries (Russia and Poland), different DNA samples (one for each SNP, positive or negative controls) were shipped from Russia to Poland and were genotyped by TaqMan assays. The results of the genotyping were in 100% agreement across the two laboratories.

Statistical analysis

Chi squared tests were used to test for the presence of Hardy-Weinberg equilibrium (HWE). HWE was tested separately for each SNP. Genotype frequencies were compared according to athletic status (i.e. controls, endurance, or sprint/power athlete) using Fisher’s exact test. Multinomial logistic regression analyses were conducted to assess the association between genotype and athletic status/competition level. Nationality was adjusted for in the first stage of analysis as there were nationality distribution differences in each athletic status groups and the control group. The homozygous non-risk allele genotype was chosen as the reference genotype for each analysis, with comparisons made to the heterozygous genotype and the homozygous risk allele genotype (co-dominant models). Additional comparisons were made to assess the dominant and recessive models, as described in our work [30]. Significance between these planned comparisons was accepted when p ≤ 0.05. Odds ratios with 95% confidence intervals were also calculated for estimation of the risk effect.

Results

Genotype frequencies distribution for IGFBP2 rs4402960, LPL rs320, LPL rs328, KCJN rs5219, and MTHFR rs1801133 amongst all participants is presented in Tables 1, 2, 3, 4, and 5. In the pooled cohort of Russian and Polish controls, genotype distributions for each of the five SNPs was in agreement with HWE (p-value > 0.05). In the Polish cohort LPL rs320 deviated from HWE (P = 0.026), however LPL rs320 was in agreement with HWE in the Russian cohort (p = 0.7) (Table 2). The analyses for all the SNPs was performed on the pooled cohort, hence the HWE deviation in the Polish cohort had no effect on the results.

Table 1 Genotype frequency distributions for IGFBP2 rs4402960

Table 2 Genotype frequency distributions for LPL rs320

Table 3 Genotype distributions for LPL rs328

Table 4 Genotype frequency distributions for KCJN rs5219

Table 5 Genotype frequency distributions for MTHFR rs1801133

IGF2BP2 rs4402960 was significantly associated with athletic status (Table 6). The control participants were less likely than sprint/power athletes to have the TT (increased risk) genotype compared to GG genotype (OR: 0.62 [0.43-0.89]; p = 0.009), and TT and GT combined (OR: 0.59 [0.42- 0.84]; p = 0.003). The sprint/power athletes were more likely than endurance to be TT compared to GG (OR: 2.11 [1.01-3.95]; p = 0.041), and GT and TT combined (OR: 2.00 [1.01- 3.95]; p = 0.045).

Table 6 Ratios of genotype distributions according to athlete type for IGFBP2 rs4402960

LPL rs320 was also significantly associated with athletic status. Table 7 shows that the control group is more likely than the endurance athletes to have the GT genotype compared to TT genotype (OR: 1.26 [1.05-1.52]; p = 0.013). Controls are also more likely than endurance to have the risk-related GG&GT genotypes compared to the TT genotype (OR: 1.22 [1.03-1.46]; p = 0.024).

Table 7 Ratios of genotype distributions according to athlete type for LPL rs320

There were no differences between the studied groups and the control group across LPL rs328, KCJN rs5219, and MTHFR rs1801133 genotypes (Tables 8, 9 and 10). Furthermore, no significantly greater/lesser odds ratios were observed for any of the genotypes in either competition level.

Table 8 Ratios of genotype distributions according to athlete type for LPL rs328

Table 9 Ratios of genotype distributions according to athlete type for KCJN rs5219

Table 10 Ratios of genotype distributions according to athlete type for MTHFR rs1801133

Finally, Tables 1, 2, 3, 4, and 5 show the percentage of genotypes present in elite-level and national-level athletes according to nationality and athletic status. No significant genotype differences were observed between elite-level and national-level athletes in all SNP and across nationalities (all p > 0.05).

Discussion

We studied the association between five obesity and co-morbidities-related genetic variants (IGF2BP2 rs4402960, LPL rs320, LPL rs328, KCJN rs5219, and MTHFR rs1801133) and athletic status in a well-defined (athletic level, ethnicity, gender) athletic population. We found a significant association between IGF2BP2 rs4402960 and LPL rs320 and athletic status; endurance athletes are less likely to have the metabolic risk IGF2BP2 T and LPL rs320 G alleles compared with sprint/power athletes and controls, respectively. These results suggest that male endurance athletes might be genetically predisposed toward a reduced risk of developing metabolic morbidities, compared with sprint/power athletes and the general population.

Previous studies have demonstrated that genetic variants associated with predisposition to obesity are also associated with responsiveness to exercise training [31-36]. Only a handful of variants, however, were replicated in multiple cohorts mainly due to variability in exercise training level, different ethnicity, gender, age, and cohorts with different metabolic states. To overcome some of the past studies challenges, including variability in physical activity status, different ethnicity and gender, we recruited a relatively- large cohort of Caucasians athletes with a well-defined athletic phenotype.

IGF2BP2, also referred to as IMP2, belongs to a mRNA-binding protein family involved in the development and stimulation of insulin action. The IGF binding protein family plays a role in modulation of IGF2 translation in a tissue-specific and developmental manner [37,38]. Several GWAS have found that carriers of the minor alleles in SNPs rs1470579 and rs4402960 have moderately increased risk for T2D. This association was confirmed across different ethnicities and populations [37-46]. Furthermore, a recent meta-analysis of 48 independent studies confirmed this association in European, East Asian and South Asian populations [47].

The intron 2 G > T substitution in the IGF2BP2 rs4402960 is particularly interesting and has attracted the most attention in obesity and T2D studies. The SNP is located in the second, large IGF2BP2 intron; thus, it is not yet clear how it generates its effect, whether directly through regulatory effects or indirectly through other genes. However, in the context of T2D, animal model and human studies implicate a role for this variant in beta-cell function, insulin secretion and sensitivity, and with elevated fasting glucose levels [18-20]. Importantly, recent studies suggest a potential role for IGF2BP2 protein in skeletal muscle cell proliferation and differentiation [21]. In the present study we have demonstrated that endurance athletes are less likely to have the metabolic risk alleles of IGF2BP2 compared to sprint/power athletes who are twice as much likely to have the metabolic risk allele (homozygote) compared to endurance athletes.

An additional finding in the present study is that endurance athletes are less likely to have the metabolic risk, G allele, of LPL rs320, compared with controls. LPL plays a pivotal role in lipid metabolism by hydrolysing triglyceride -rich lipoproteins. Dysfunction of LPL protein increased the susceptibility for developing several common diseases, including atherosclerosis and obesity [22,47-50]. LPL rs320 or HindIII (intron 8) is a common variant in the LPL gene that has been associated with plasma lipid profile [22,24,51-54]. Although a large number of variants have been identified in the LPL gene, rs320 is of particular interest because of its common occurrence in many populations. Due to LPL rs320′s location within an intron, it was not initially considered functional but rather in linkage disequilibrium with a putative functional variant, such as LPL rs328. However, recent findings suggests that the LPL rs320 may be functional by altering the binding of a transcription factor and impacting LPL expression [49]. We found that sedentary controls are more likely to have the risk variant compared with endurance athletes and thus, might in more risk to develop elevated blood lipids and Cardio Vascular Disease [55].

A possible explanation to the underrepresentation of metabolic diseases risk alleles in endurance athletes arising from studies that evaluated the overall risk of athletes for metabolic and cardiovascular disease. Guo et al., [56] have shown that professional strength-oriented athletes at the heaviest-weight-class are at a significant increased risk for cardiometabolic disease compared with those at all other weight categories. Similarly, Urho et al., [57] found that, compared with controls, strength/power-sports athletes had a higher risk for high body mass index (BMI), whereas former endurance athletes had the lowest odds ratios for T2D and ischemic heart disease. These studies reinforce our hypothesis that endurance athletes would be at lower risk for complex metabolic diseases compared to sprint/power athletes, and controls, and genetics might be, at least partly, behind these differences.

Conclusions

In conclusion, we found a significant association between IGF2BP2 and LPL SNPs and athletic status in males: endurance athletes are less likely to have the metabolic risk alleles of IGF2BP2 rs4402960 and LPL rs320, compared to sprint/power athletes and controls. These results suggest that some SNPs across the human genome have dual effect and may predispose endurance athletes to reduced risk of developing metabolic morbidities, whereas sprint/power athletes might be predisposed to elevated risk. These results need to be confirmed in athlete cohorts with different geographical backgrounds. Future studies should also measure obesity-related intermediate phenotypes, such as fasting blood glucose levels and plasma lipids that could lend support for the associations.