A Potential Endurance Algorithm Prediction in the Field of Sports Performance

Abstract

Sport performance is influenced by several factors, including genetic susceptibility. In the past years, specific single nucleotide polymorphisms have been associated to sport performance; however, these effects should be considered in multivariable prediction systems since they are related to a polygenic inheritance. The aim of this study was to design a genetic endurance prediction score (GES) of endurance performance and analyze its association with anthropometric, nutritional and sport efficiency variables in a cross-sectional study within fifteen male cyclists. A statistically significant positive relationship between GES and the VO₂ maximum (P = 0.033), VO₂ VT1 (P = 0.049) and VO₂ VT2 (P < 0.001) was observed. Moreover, additional remarkable associations between genotype and the anthropometric, nutritional and sport performance variables, were achieved. In addition, an interesting link between the habit of consuming caffeinated beverages and the GES was observed. The outcomes of the present study indicate a potential use of this genetic prediction algorithm in the sports’ field, which may facilitate the finding of genetically talented athletes, improve their training and food habits, as well as help in the improvement of physical conditions of amateurs.

Introduction

Athletic performance can be influenced by several factors. These can include extrinsic factors such as the quality of training, the dietary habits, the technology used or the weather conditions (Sacha and Quinn, 2011). Besides, intrinsic factors such as individual genetic profile also play an important role (Peveler and Green, 2010). In fact, the study of the genetic influence on sports performance has become a leading field of research.

Single nucleotide polymorphisms (SNPs) are among the genetic variants implicated in the phenotypic differences that can influence individual physical abilities (Santos et al., 2016). To date, various links between several genetic variants and favorable phenotypes for certain sports have been established (Santos et al., 2016). This suggests that the presence of specific genotypes can predispose an individual to competitive advantages in a particular sport (Zilberman-Schapira et al., 2012).

For example, backed by a meta-analysis, the II genotype of the angiotensin converting enzyme ACE insertion/deletion (I/D) rs4340 has been significantly associated with endurance’s discipline as compared with the ID and DD genotypes (Ma et al., 2013). Specifically, it is suggested that the I allele implies higher blood flow and sugar utilization by muscles, that facilitates endurance performance (Woods et al., 2000). Additionally, considering reported difficulties in genotyping ACE I/D, the ACE rs4343 A and G alleles are admitted as equivalents to the ACE rs4340 I and D alleles, respectively, and considered an alternative method for genotyping of the ACE I/D polymorphism (Glenn et al., 2009).

Moreover, the different genotypes of the peroxisome proliferator activated receptor alpha PPARA rs4253778 have also been associated with athlete’s performance. Concretely, in a meta-analysis carried out by Lopez-Leon et al. (2016) the GG genotype and G allele were significantly more common in endurance athletes as compared to controls.

Besides, the α-actinin 3 (ACTN3) is a sarcomeric actin-binding protein specifically expressed in fast twitch myofibers of the skeletal muscle, required for explosive muscle contraction (Wilson et al., 2012). But the gene encoding this protein also seems to be associated with endurance capacities. In particular, TT rs1815739 carriers were found to be more common in endurance athletes as compared with sedentary individuals in a sample of 395 Israeli (Eynon et al., 2009). They were also described to obtain better results in an endurance test (Pasqua et al., 2016), to exhibit a higher proportion of endurance-associated type I myofibers and to prefer to skate long- than short-distance races (Ahmetov et al., 2011).

Another gene with a potential association with endurance capacities is the Aquaporin 1 AQP1 rs1049305 where C allele carriers were faster in running performances than carriers of the GG genotype (Martinez et al., 2009; Rivera et al., 2011; Saunders et al., 2015).

Moreover, the CC genotype of the peroxisome proliferator activated receptor gamma, coactivator 1 alpha PPARGC1A rs8192678 has been described to be associated with high values of aerobic performance (Akhmetov et al., 2007; Stefan et al., 2007; Eynon et al., 2010b).

And finally, other SNPs such as the beta-3 adrenergic receptor ADRB3 rs4994 (Santiago et al., 2011), the GA-binding protein transcription factor, subunit beta 1 GABPB1 rs12594956 (Eynon et al., 2010a, 2013), the collagen type V alpha 1 chain COL5A1 rs12722 (Brown et al., 2011; Posthumus et al., 2011) and the hemochromatosis HFE rs1799945 (Grealy et al., 2015) have also been associated with better endurance performance, although the scientific evidence support for these associations is still scarce.

Thus, it is clear that there is a genetic influence, but there is still a weak scientific evidence for most of the reported associations. Besides, most studies that associate genetics with physical capacities have focused on individual genes. However, as a polygenic inheritance, various genes can make a contribution to the overall outcome. For these reasons, we consider that it is necessary to create more complex prediction algorithms, including different genetic factors together. Therefore, in the present study, we have developed a predictive algorithm of endurance performance including 11 genes. In this report, we present the associations of the predictive algorithm and each of the SNPs with anthropometric, nutritional and sport performance variables, in a group of 15 semi-professional cyclists.

Materials and Methods

Subjects

Fifteen healthy male members of the Spanish Cycling Federation, body mass index (BMI) 22.3 ± 2.5 aged 40.7 ± 7.0, with at least 5 years of experience in national-level competitions were recruited by the Sports Medicine University Center (Complutense University of Madrid, Madrid, Spain) to participate in the present study. All participants were non-smokers. This research was conducted according to the guidelines laid down in the Declaration of Helsinki and all procedures involving human subjects were approved by the Research Ethics Committee of the IMDEA Food Foundation (PI-0031). Written informed consent to participate in the study was obtained from all subjects.

Design

This was a cross-sectional clinical study where volunteers attended two different centers to complete the study:

• 1)

  The Sports Medicine University Center to carry out a maximal incremental treadmill test.

• 2)

  The Research Institute on Food and Health Sciences “IMDEA Food” (Madrid, Spain) for anthropometric measurements, body composition analysis, dietary records, DNA collection and genotyping.

Methodology

Cardiopulmonary Exercise Test

All participants carried out an incremental (30 w/min) exercise testing on a cycle ergometer (Cardgirus, Barcelona, Spain) after 2 h of fasting. During the test, the heart rate (HR) was measured using a 12-lead wireless electrocardiograph (Norav, Wiesbaden, Germany).

Oxygen uptake (VO₂), carbon dioxide output (VCO₂) and minute ventilation (VE) were assessed using the analyzer Jaeger Oxycon-Pro (Hoechberg, Germany). The respiratory exchange ratio (RER) was calculated as VCO₂/VO₂, while the VO₂ pulse as VO₂/HR (Bergh et al., 2000). Every 2 min cyclists had to estimate their feelings of exertion and pain using the Borg Rating of Perceived Exertion (RPE) Scale (Borg et al., 1985).

All participants achieved maximal exercise criteria: VO₂ plateau (considered as VO₂ max) which was estimated as a VO₂ increase lower than 150 mL/min for two consecutive periods (Petot et al., 2012), RPE higher than 16 (Bergh et al., 2000), RER above 1.10 and a HR upper 90% of the theoretical maximum HR (Howley et al., 1995). The maximum HR was calculated as [(208.75−(0.73 × age)] (Tanaka et al., 2001).

The aerobic ventilatory threshold (VT1) was estimated using the criteria of the ventilatory equivalent for VO₂ (VE/VO₂ ratio), corresponding to the rupture of the linearity in the increment of VE. Finally, the anaerobic ventilatory threshold (VT2) was calculated by the increase of the VCO₂ equivalent principles (VE/VCO₂), as the second rupture of the linearity in the increment of VE (Beaver et al., 1986).

Anthropometry and Lifestyle Parameters

Anthropometric measurements were determined while subjects were wearing light clothing and no shoes. Height was assessed to the nearest 0.1 cm using a stadiometer (Leicester-Biological Medical Technology SL, Barcelona). Body weight, fat mass and muscle mass percentages were evaluated using a BF511 Body Composition Monitor (BF511- OMRON Healthcare UK, LT, Kyoto, Japan). Brachial, contracted arm, waist, hip and leg circumferences were measured with an inextensible tape (KaWe Kirchner & Wilhelm GmbH, Asperg, Germany; range 0–150 cm, 1 mm of precision). A caliper (Holtain Ltd., Crymych, United Kingdom; 10 g/mm² constant pressure; range 0–39 mm and 0.1 mm of precision) was used for biceps, triceps, subscapular, abdominal, supraspinal, front thigh and medial calf skinfolds determinations. Moreover, the diameters of the femur and humerus were also assessed using a small bone caliper (Nonio sliding Bicondyleo, Holtain Ltd., United Kingdom). Systolic and diastolic blood pressures were evaluated using an automatic digital blood pressure monitor Model M3 (OMRON Healthcare UK, LT, Kyoto, Japan) in the right arm, with the patient seated and relaxed. Measurements were taken three times after a 5-min resting period, following World Health Organization (WHO) criteria (Whitworth and Chalmers, 2004). Finally, with the different anthropometric data, the BMI was calculated as the body weight divided by the squared height (kg/m²) and somatotype values (endomorphic, mesomorphic, and ectomorphic values) accordingly to Heath-Carter method (Carter and Heath, 1990).

The food habits of each participant were recorded using a validated 3-day dietary food record and a food frequency questionnaire (Aguirre-Jaime et al., 2008). Subsequently, the composition of the different dietary records was analyzed using the DIAL software (2.16 version Alce Ingeniería, Madrid, Spain). For the calculation of the Healthy Eating Index score (Guenther et al., 2013), the DIAL program gives different values ranging from 0 to 100 considering the daily servings of cereals, vegetables, fruits, dairy products, and meat; the percentage of energy provided by total and saturated fats; the amount of cholesterol and sodium per day and the number of different foods consumed. The final score is classified into five categories: an “excellent diet” (>80 points), a “very good diet” (71–80 points), a “good diet” (61–70 points), an “acceptable diet” (51–60 points), or an “inadequate diet” (0–50 points).

Genotyping of the Population

A sample of 500 μl of peripheral capillary blood of each volunteer was drawn for DNA extraction. To perform the subsequent analysis of the samples, genomic DNA was extracted from the cellular fraction collected by the Genomic QIAamp DNA Blood Kit Mini Kit (QIAGEN, Spain). The samples were genotyped with TaqMan Assays by the high-performance QuantStudio Real-Time PCR (Applied Biosystem, United States).

Statistical Analysis

Data were analyzed using the R Statistical Software Version 3.4.1¹. The description of the qualitative data was made in the form of absolute frequencies and percentages and the quantitative data by mean and standard deviation. The Mann–Whitney U test was used to check for significant differences in the continuous variables (not always normally distributed) for the different genotypes. The Spearman correlation coefficient was used for the association between the algorithm and the other variables. The Bonferroni correction was also applied to control against type-I errors for multiple tests. All the statistical tests were two-tailed. Statistical significance was assumed when P < 0.05.

Selection of SNPs and Design of the GES

An exhaustive literature review of the scientific databases (Pubmed, Medline, Web of Sciences) was carried out to identify all studies that analyzed the relationship between one or more SNP and sports performance. A selection of 11 SNPs was made considering the European frequencies of each SNP according to Ensembl database, the scientific evidence of each association and its availability for TaqMan SNP Genotyping Assay. Among them, nine SNPs were associated with endurance capacities and two SNPs were mainly related to power abilities.

Once the SNPs selection was completed, an algorithm to predict endurance capacities based on the nine SNPs associated with endurance performance was established (Table 1). Depending on the scientific evidence, each SNP was given a different weight to the total GES. Therefore, the two SNPs with the highest scientific evidence based on meta-analysis, ACE and PPARα, were given the highest normalized weight of the GES, contributing 22.2% each of the total value of the algorithm. On the other hand, those SNPs whose evidence was based on at least three studies with positive and conclusive results (ACTN3, AQP1, and PPARGC1A) were given a weight of 11.1% each, of the total value of the algorithm. Besides, to the remaining 4 SNPs, with less scientific evidence and some contradictory results, each one was given a weight of 5.6% to the final GES. Finally, based on the literature, the three possible genotypes of each SNP were classified into three categories: “favorable genotypes” (GES weight × 1), which were those associated with better endurance performance; “intermediate genotypes” (GES weight × 0.5), which were those with a neutral effect on endurance abilities, and “unfavorable genotypes” (GES weight × 0), which were those with a negative effect. For example, an individual with all favorable SNPs, would have a GES value of 100: 22.2 × 1 (ACE) + 22.2 × 1 (PPARα) + 11.1 × 1 (ACTN3) + 11.1 × 1 (AQP1) + 11.1 × 1 (PPARGC1A) + 5.6 × 1 (ADRB3) + 5.6 × 1 (GABPB1) + 5.6 × 1 (COL5A1) + 5.6 × 1 (HFE) = 100. However, an individual with unfavorable ACE and PPARα SNPs and all the other SNPs favorable, would have a score of 56: 22.2 × 0 (ACE) + 22.2 × 0 (PPARα) + 11.1 × 1 (ACTN3) + 11.1 × 1 (AQP1) + 11.1 × 1 (PPARGC1A) + 5.6 × 1 (ADRB3) + 5.6 × 1 (GABPB1) + 5.6 × 1 (COL5A1) + 5.6 × 1 (HFE) = 56.

The selected power-related genotypes were:

• •

  Hypoxia inducible factor 1 alpha subunit HIF1A rs11549465: CC (81.1%) unfavorable, CT (17.7%) intermediate, TT (1.2%) favorable.

• •

  Muscle-specific creatine kinase CKM rs8111989: TT (48.3%) unfavorable, CT (43.3%) intermediate, CC (8.3%) favorable.

Results

Descriptive Analysis

Mean values of anthropometric, body composition, somatotype and blood pressure of all cyclists are shown in Table 2. As expected, the mean BMI and fat mass percentage values were in the range of “normal weight” [18.5–24.9 BMI, 12–20% fat mass percentage according to the Spanish Society for the Study of Obesity (Salas-Salvado et al., 2007)]. The three somatotype components were around the moderate 3–5 rate (Carter and Heath, 1990).

Moreover, the results of the analysis of the 3-day dietary records can be observed in Table 3. Here, it can be appreciated how the energy requirements of the cyclists were slightly higher than the energy intakes, these lasts with a high variability between the participants. Moreover, the mean Healthy Eating Index resulted in a “good diet” which, as explained in the methodology, is considered an intermediate value (61–70 points). This table also shows the average of servings per group of food consumed according to the Validated Food Frequency Questionnaire, where we can observe that vegetables and fruits are the most consumed.

Finally, Table 4 describes the mean results of the cardiopulmonary exercise test carried out at the Sports Medicine University Center.

Association Analyses

The GES Correlates With the Cardiovascular Exercise Test

Once established the prediction algorithm (Table 1), it was related to the results obtained from the functional cardiovascular exercise test (Table 2). The association analysis of the GES and the variables obtained in the cardiovascular exercise test, revealed a statistically significant positive relationship between the GES results and the VO₂ max (P = 0.033, Figure 1), VO₂ in VT1 (P = 0.049, ρ = 0.516) and VO₂ in VT2 (P < 0.001, ρ = 0.813) values.

FIGURE 1

As expected, SNPs included in the GES were also individually associated with specific parameters related to individual performance. The analysis of the association of the different genotypes individually and the cardiopulmonary exercise test showed that GG genotypes for AQP1 showed lower levels of VO₂ in VT1 (32.1 ± 4.6 mL/kg/min vs 42.3 ± 7.4 mL/kg/min) and fewer values of% VO₂ in VT1 in relation to VO₂ max (58.2 ± 7.0% vs 70.1 ± 7.7%) than CC + CG (P = 0.020 and P = 0.030, respectively) (Figures 2A,B). On the other hand, GG genotypes for the PPARα presented higher values of VO₂ in VT2 than CC + CG (52.6 ± 4.6 mL/kg/min vs 42.8 ± 5.9 mL/kg/min, P = 0.030), as shown in Figure 2C.

FIGURE 2

Regarding the GABPB1, AA genotypes had significantly lower maximal HR (Figure 2D) and minor HR in VT1 and VT2 than AC genotypes (167.8 ± 16.5 bpm vs 185.9 ± 4.5 bpm, 120.8 ± 11.0 bpm vs 144.1 ± 7.8 bpm and 149.0 ± 13.3 bpm vs 167.6 ± 6.3 bpm; P < 0.001, P-adjusted for Bonferroni = 0.040 and P = 0.001, respectively). There were no CC genotypes in the studied sample. Moreover, a statistically significant difference was found between AA genotypes for the ACE gene and the AG + GG with respect to the maximum VE (200 ± 16 L/min vs 162 ± 19 L/min, P = 0.040, Figure 2E).

HIF1A Genetic Variant Is Associated With Somatotype

When we studied the relationship between the different SNPs and the results of the somatotype, we found a statistically significant association between HIF1A genotypes and the mesomorphic component. Precisely, among the CC cyclists, 8.33% presented a low mesomorphic value, 91.7% a moderate value and 0% a high value; while amongst the CT individuals, 33.33% presented a moderate result and the rest 66.7% a high value (P = 0.029). There were no TT individuals in the studied sample.

Genes Associated With Dietary Records

When analyzing the association of the different genotypes of each SNP and the results of the dietary variables, we found different interesting associations.

• •

  Subjects GG for the AQP1 presented a lower intake of carbohydrates [% total caloric value (TCV)] than CG + CC individuals (34.2 ± 11.4% vs 46.9 ± 8.5%, P = 0.040), as shown in Figure 3A.

  FIGURE 3

  

• •

  Participants TT for the CKM had a higher consumption of sugars (% TCV) than CC + CT individuals (25.9 ± 9.2% vs 16.5 ± 3.9%, P = 0.010, Figure 3B).

• •

  Cyclists AA for the ACE presented a higher mean glycemic index consumption than AG + GG (58.9 ± 1.1 vs 46.5 ± 11.7, P = 0.010, Figure 3C).

On the other hand, with regard to the food frequency questionnaire, and inverse correlation between the number of coffee and tea rations consumed per day and the result of the GES was observed (P = 0.004), so that the higher the consumption the lower the GES (Figure 4).

FIGURE 4

Discussion

The present research provides new information regarding the link between genetics and sport performance from different angles.

The most important result was the potential validation of an algorithm prediction of genetic susceptibility to endurance abilities. At present, cardiovascular exercise test is considered the gold standard assessment of endurance performance (Hausen et al., 2018). We observed that those subjects with a higher value in the GES, presented significantly better results in the cardiovascular exercise test according to VO₂ in VT1, VO₂ in VT2 and most importantly, VO₂ max which is internationally considered the key measure of aerobic capacity (Hausen et al., 2018). Actually, VO₂ max obtained in this kind of test represents the highest capacity of oxygen consumption during a maximal exercise (Shete et al., 2014). Additionally, VT1 is the point at which the aerobic metabolism is no longer the major energy source and the anaerobic metabolism begins to be used. The exercise intensity at which VT1 occurs is considered the highest submaximal level tolerated by an individual for long time periods (Herdy et al., 2016). Hence, we hypothesized that higher values of VO₂ at VT1 may predict higher VO₂ max, which leads to better endurance conditions. Concerning VT2, it is the point where the demand of oxygen by muscles exceeds the mitochondrial supplies and consequently, energy production begins predominantly anaerobic. It begins when lactate starts to accumulate in blood which is usually at 60–70% VO₂ max (Albouaini et al., 2007). The VO₂ at VT2 is considered a submaximal index of exercise capacity and endurance performance, so that the higher VO₂ at VT2, the better endurance results (Coyle et al., 1988; Herdy et al., 2016). These achievements may complement the association found between a genetic score based on a GWAS study and the improvement of the VO₂ max in a sedentary population carried out by previous researchers (Bouchard et al., 2011).

When we analyzed the association between the different genotypes individually and the cardiovascular test results, we also found some interesting correlations. In this regard, individuals presenting the favorable genotypes for endurance sports according to AQP1 (Martinez et al., 2009; Saunders et al., 2015), showed better results in the cardiovascular exercise test (higher levels of VO₂ in VT1 and % of VO₂ in VT1 in relation to VO₂ max), than GG homozygous. The physiological explanation for this association may lie in the AQP1 encoded membrane protein role. The aquaporin 1 water channel is implicated in the transport of water, maintaining the osmotic balance between the blood and the cells (Frigeri et al., 2004). During prolonged exercise, body temperature is regulated by controlling the distribution of body fluid; water passes from the intracellular to the extracellular spaces and evaporates by sweating. An inadequate loss of sweat during exercise, especially in hot climates, will negatively affect athletic performance (Kenny, 2014). Thus, we hypothesize that CC genotype is associated with an efficient state of cellular hydration and body temperature regulation which leads to a better endurance performance.

Moreover, favorable genotypes for endurance performance of the PPARα (Lopez-Leon et al., 2016) presented significantly higher values of VO₂ in VT2 than CC + CG. This gene regulates the expression of other genes involved in the metabolism of energy, fats and sugars in the skeletal muscle among other tissues (Duval et al., 2004). Given its involvement in these processes, it is hypothesized that it is activated during endurance exercise (Lopez-Leon et al., 2016). Another suggested explanation for the association of PPARα genotypes with endurance performance has to do with the type of fibers in the skeletal muscle. In this sense, in a cohort of 786 Russian athletes it was observed that GG homozygous presented significantly higher percentage of type I muscle fibers than the other genotypes (Ahmetov et al., 2006). These skeletal muscle fibers are classified into two types: type I or slow contraction and type II or rapid contraction fibers. Type I fibers have greater resistance to fatigue and predominate in resistance athletes, while the type II ones are adapted to strong and explosive muscle contractions and predominate in athletes who perform power sports (Cartee et al., 2016).

Regarding the GABPB1, the individuals presenting the AA genotype showed a lower maximal HR and lower HR in VT1 and VT2 than AC genotypes. This result is in accordance with other studies that also consider this genotype favorable for endurance sports (Eynon et al., 2013) as it has been described that endurance athletes present lower HR in maximum aerobic traits (Zaniqueli et al., 2014). The GABPB1 gene encodes the GA-binding protein transcription factor, which is implicated in the regulation of the mitochondrial function generating ATP energy (Dinkova-Kostova and Abramov, 2015) what might explain the implication of this gene in the individual endurance capacity.

Additionally, AA genotypes for the ACE gene presented higher maximum VE than the AG + GG genotypes. The AA genotype of the rs4343 is considered equivalent to the II genotype of the ACE I/D rs4340 (Glenn et al., 2009) which at the same time is the most studied favorable genotype for endurance performance (Ma et al., 2013). Maximal VE has been reported to be directly correlated with VO₂ max (Keramidas et al., 2010; Malekmohammad et al., 2012). Accordingly, it is suggested that muscles of endurance athletes require higher values of VO₂ so that their VE during exercise is higher.

Hypoxia-inducible factor-1 (HIF1) regulates oxygen homeostasis in mammalian cells and in particular, it seems to have a role during high intensity exercise, helping the skeletal muscle to adapt to low oxygen concentrations (Freyssenet, 2007). According to the association between somatotype variables and the different genotypes, individuals carrying the T allele of the HIF1A SNP presented moderate or high mesomorphic values, while none of the CC homozygous had a high mesomorphic component. Precisely, allele T is associated with power-oriented athletes (Drozdovska et al., 2013). This makes sense as high mesomorphic individuals are characterized by high skeletal muscle mass, needed for power anaerobic exercises (Gutnik et al., 2015). Besides, as cycling is more an endurance-oriented sport than a power sport, it also makes sense that the mesomorph was the component with the lowest value among the sample.

When we analyzed the association of the different genotypes and the dietary records, we also found diverse significant associations. The cyclists with a favorable genotype for endurance sports regarding the AQP1 presented a higher intake of carbohydrates which we hypothesize that would probably be due to maintain the glycogen stores needed for long distance exercise (Alghannam et al., 2018). Similarly, individuals with an unfavorable genotype for power sports according to CKM consumed higher amounts of sugars (Chen et al., 2017). Commonly, it is considered that genotypes unfavorable for power are favorable for endurance sports and vice versa, which might explain this association. To our knowledge, this is the first time that AQP1 and CKM genotypes have been associated with dietary intake.

Moreover, individuals presenting a favorable genotype for endurance sports according to ACE gene presented a higher mean glycemic index intake than the other genotypes. Although it seems clear that carbohydrate consumption is needed to maintain glycogen stores, whether these macronutrients are preferable to be complex or with a high glycemic index appears to be controversial, as there are studies that point out that moderate glycemic index diets improves exercise performance (Durkalec-Michalski et al., 2017) while others support the low glycemic index (Durkalec-Michalski et al., 2018). A possible explanation of this controversy may lie in the fact that the unfavorable ACE genotype for endurance performance has also been linked to lower glucose tolerance (Schuler et al., 2017). Thus, we hypothesize that ACE AA individuals might take advantage in endurance sports by consuming a moderate glycemic index diet as they can metabolize glucose in a better way. However, we are aware that more studies in this field are needed to confirm this hypothesis.

Finally, an interesting association between the habit of consuming caffeinated beverages and the result of the GES was observed. A possible explanation for this inverse correlation might be that athletes less genetically predisposed to endurance sports where energy is a key factor, are more likely to use ergogenic aids such as caffeine. However, we have to consider that the food frequency questionnaire used gives overage data of the last year, but we do not have data on dietary intakes in specific time frames such as competitions.

Conclusion

The outcomes of the present study confirm a positive relationship between an endurance prediction algorithm and the results of a cardiopulmonary exercise test. Moreover, AQP1, PPARα, GABPB1, and ACE genes were individually related with endurance performance. Besides, HIF1A showed an association with the somatotype and AQP1, CKM, and ACE genes were associated with the athletes’ dietary intake. In addition, an inverse association between the habit of consuming caffeinated beverages and the GES was observed.

This information may facilitate the design of larger studies implicated in the prediction of sports capacities, which may facilitate the finding of genetically talented athletes, improve their training and dietary habits, as well as help in the improvement of physical conditions of amateur athletes.

References

• 1

  Aguirre-JaimeA.Cabrera de LeonA.Dominguez CoelloS.Borges AlamoC.Carrillo FernandezL.Gavilan BatistaJ. C.et al (2008). Validation of a food intake frequency questionnaire adapted for the study and monitoring of the adult population of the Canary Islands, Spain.Rev. Esp. Salud Publ.82509–518. 10.1590/s1135-57272008000500006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  AhmetovI. I.DruzhevskayaA. M.LyubaevaE. V.PopovD. V.VinogradovaO. L.WilliamsA. G. (2011). The dependence of preferred competitive racing distance on muscle fibre type composition and ACTN3 genotype in speed skaters.Exp. Physiol.961302–1310. 10.1113/expphysiol.2011.060293

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  AhmetovI. I.MozhayskayaI. A.FlavellD. M.AstratenkovaI. V.KomkovaA. I.LyubaevaE. V.et al (2006). PPARalpha gene variation and physical performance in Russian athletes.Eur. J. Appl. Physiol.97103–108. 10.1007/s00421-006-0154-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  AkhmetovI. I.PopovD. V.MozhaiskaiaI. A.MissinaS. S.AstratenkovaI. V.VinogradovaO. L.et al (2007). Association of regulatory genes polymorphisms with aerobic and anaerobic performance of athletes.Ross Fiziol Zh Im I M Sechenova93837–843.

  • Google Scholar

• 5

  AlbouainiK.EgredM.AlahmarA.WrightD. J. (2007). Cardiopulmonary exercise testing and its application.Postgrad. Med. J.83675–682. 10.1136/hrt.2007.121558

  • CrossRef
  • Google Scholar

• 6

  AlghannamA. F.GonzalezJ. T.BettsJ. A. (2018). Restoration of muscle glycogen and functional capacity: role of post-exercise carbohydrate and protein co-ingestion.Nutrients10:253. 10.3390/nu10020253

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  BeaverW. L.WassermanK.WhippB. J. (1986). A new method for detecting anaerobic threshold by gas exchange.J. Appl. Physiol.602020–2027. 10.1152/jappl.1986.60.6.2020

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8

  BerghU.EkblomB.AstrandP. O. (2000). Maximal oxygen uptake “classical” versus “contemporary” viewpoints.Med. Sci. Sports Exerc.3285–88. 10.1097/00005768-200001000-00013

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  BorgG.LjunggrenG.CeciR. (1985). The increase of perceived exertion, aches and pain in the legs, heart rate and blood lactate during exercise on a bicycle ergometer.Eur. J. Appl. Physiol. Occup. Physiol.54343–349. 10.1007/bf02337176

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  BouchardC.SarzynskiM. A.RiceT. K.KrausW. E.ChurchT. S.SungY. J.et al (2011). Genomic predictors of the maximal O(2) uptake response to standardized exercise training programs.J. Appl. Physiol.1101160–1170. 10.1152/japplphysiol.00973.2010

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  BrownJ. C.MillerC. J.PosthumusM.SchwellnusM. P.CollinsM. (2011). The COL5A1 gene, ultra-marathon running performance, and range of motion.Int. J. Sports Physiol. Perform.6485–496. 10.1123/ijspp.6.4.485

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12

  CarteeG. D.AriasE. B.YuC. S.PatakyM. W. (2016). Novel single skeletal muscle fiber analysis reveals a fiber type-selective effect of acute exercise on glucose uptake.Am. J. Physiol. Endocrinol. Metab.311E818–E824. 10.1152/ajpendo.00289.2016

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  CarterJ.HeathB. (1990). Somatotyping: Development And Implications.Cambridge: Cambridge University Press.

  • Google Scholar

• 14

  ChenC.SunY.LiangH.YuD.HuS. (2017). A meta-analysis of the association of CKM gene rs8111989 polymorphism with sport performance.Biol. Sport34323–330. 10.5114/biolsport.2017.69819

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  CoyleE. F.CogganA. R.HopperM. K.WaltersT. J. (1988). Determinants of endurance in well-trained cyclists.J. Appl. Physiol.642622–2630. 10.1152/jappl.1988.64.6.2622

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  Dinkova-KostovaA. T.AbramovA. Y. (2015). The emerging role of Nrf2 in mitochondrial function.Free Radic. Biol. Med.88179–188. 10.1016/j.freeradbiomed.2015.04.036

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  DrozdovskaS. B.DosenkoV. E.AhmetovI. I.IlyinV. N. (2013). The association of gene polymorphisms with athlete status in ukrainians.Biol. Sport30163–167. 10.5604/20831862.1059168

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  Durkalec-MichalskiK.ZawiejaE. E.ZawiejaB. E.JurkowskaD.BuchowskiM. S.JeszkaJ. (2018). Effects of low versus moderate glycemic index diets on aerobic capacity in endurance runners: three-week randomized controlled crossover trial.Nutrients10:370. 10.3390/nu10030370

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  Durkalec-MichalskiK.ZawiejaE. E.ZawiejaB. E.PodgorskiT.JurkowskaD.JeszkaJ. (2017). Influence of low versus moderate glycemic index of diet on substrate oxidation and energy expenditure during incremental exercise in endurance athletes: a randomized counterbalanced cross-over trial.Int. J. Food Sci. Nutr.69741–752. 10.1080/09637486.2017.1411891

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  DuvalC.FruchartJ. C.StaelsB. (2004). PPAR alpha, fibrates, lipid metabolism and inflammation.Arch. Mal. Coeur. Vaiss97665–672.

  • Google Scholar

• 21

  EynonN.AlvesA. J.SagivM.YaminC.SagivM.MeckelY. (2010a). Interaction between SNPs in the NRF2 gene and elite endurance performance.Physiol. Genomics4178–81. 10.1152/physiolgenomics.00199.2009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  EynonN.MeckelY.SagivM.YaminC.AmirR.SagivM.et al (2010b). Do PPARGC1A and PPARalpha polymorphisms influence sprint or endurance phenotypes?Scand. J. Med. Sci. Sports20e145–e150. 10.1111/j.1600-0838.2009.00930.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  EynonN.DuarteJ. A.OliveiraJ.SagivM.YaminC.MeckelY.et al (2009). ACTN3 R577X polymorphism and Israeli top-level athletes.Int. J. Sports Med.30695–698. 10.1055/s-0029-1220731

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  EynonN.RuizJ. R.BishopD. J.SantiagoC.Gomez-GallegoF.LuciaA.et al (2013). The rs12594956 polymorphism in the NRF-2 gene is associated with top-level Spanish athlete’s performance status.J. Sci. Med. Sport16135–139. 10.1016/j.jsams.2012.05.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 25

  FreyssenetD. (2007). Energy sensing and regulation of gene expression in skeletal muscle.J. Appl. Physiol.102529–540. 10.1152/japplphysiol.01126.2005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  FrigeriA.NicchiaG. P.BalenaR.NicoB.SveltoM. (2004). Aquaporins in skeletal muscle: reassessment of the functional role of aquaporin-4.FASEB J.18905–907. 10.1096/fj.03-0987fje

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  GlennK. L.DuZ. Q.EisenmannJ. C.RothschildM. F. (2009). An alternative method for genotyping of the ACE I/D polymorphism.Mol. Biol. Rep.361305–1310. 10.1007/s11033-008-9313-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  GrealyR.HerruerJ.SmithC. L.HillerD.HaselerL. J.GriffithsL. R. (2015). Evaluation of a 7-gene genetic profile for athletic endurance phenotype in ironman championship triathletes.PLoS One10:e0145171. 10.1371/journal.pone.0145171

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  GuentherP. M.CasavaleK. O.ReedyJ.KirkpatrickS. I.HizaH. A.KuczynskiK. J.et al (2013). Update of the healthy eating index: HEI-2010.J. Acad. Nutr. Diet.113569–580. 10.1016/j.jand.2012.12.016

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  GutnikB.ZuozaA.ZuozieneI.AlekrinskisA.NashD.ScherbinaS. (2015). Body physique and dominant somatotype in elite and low-profile athletes with different specializations.Medicina51247–252. 10.1016/j.medici.2015.07.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  HausenM.SoaresP. P.AraujoM. P.EstevesD.JulioH.TauilR.et al (2018). Specificity elicits higher maximal and submaximal cardiorespiratory responses during a new taekwondo aerobic test.Int. J. Sports Physiol. Perform.5667–71. 10.1123/ijspp.2017-0846

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  HerdyA. H.RittL. E.SteinR.AraujoC. G.MilaniM.MenegheloR. S.et al (2016). Cardiopulmonary exercise test: background, applicability and interpretation.Arq. Bras. Cardiol.107467–481. 10.5935/abc.20160171

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 33

  HowleyE. T.BassettD. R.Jr.WelchH. G. (1995). Criteria for maximal oxygen uptake: review and commentary.Med. Sci. Sports Exerc.271292–1301.

  • Google Scholar

• 34

  KennyG. P. (2014). Muscle temperature and sweating during exercise: a new link?Acta Physiol.21211–13. 10.1111/apha.12335

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  KeramidasM. E.DebevecT.AmonM.KounalakisS. N.SimunicB.MekjavicI. B. (2010). Respiratory muscle endurance training: effect on normoxic and hypoxic exercise performance.Eur. J. Appl. Physiol.108759–769. 10.1007/s00421-009-1286-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  Lopez-LeonS.TuvbladC.ForeroD. A. (2016). Sports genetics: the PPARA gene and athletes’ high ability in endurance sports. a systematic review and meta-analysis.Biol. Sport333–6. 10.5604/20831862.1180170

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  MaF.YangY.LiX.ZhouF.GaoC.LiM.et al (2013). The association of sport performance with ACE and ACTN3 genetic polymorphisms: a systematic review and meta-analysis.PLoS One8:e54685. 10.1371/journal.pone.0054685

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  MalekmohammadM.Ahmadi-NejadM.AdimiP.JamaatiH. R.MarashianS. M. (2012). Evaluation of maximum O₂ consumption: using ergo-spirometry in severe heart failure.Acta Med. Iran.50619–623.

  • Google Scholar

• 39

  MartinezJ. L.CarrionA.FlorianM. E.MartinJ. A.Lopez-TaylorJ. R.FaheyT. D.et al (2009). Aquaporin-1 gene DNA variation predicts performance in Hispanic marathon runners.Med. Sport.13251–255. 10.2478/v10036-009-0039-9

  • CrossRef
  • Google Scholar

• 40

  PasquaL. A.BuenoS.ArtioliG. G.LanchaA. H.Jr.MatsudaM.MarqueziniM. V.et al (2016). Influence of ACTN3 R577X polymorphism on ventilatory thresholds related to endurance performance.J. Sports Sci.34163–170. 10.1080/02640414.2015.1040823

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  PetotH.MeillandR.Le MoyecL.Mille-HamardL.BillatV. L. (2012). A new incremental test for VO(2)max accurate measurement by increasing VO(2)max plateau duration, allowing the investigation of its limiting factors.Eur. J. Appl. Physiol.1122267–2276. 10.1007/s00421-011-2196-5

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  PevelerW. W.GreenM. (2010). The effect of extrinsic factors on simulated 20-km time trial performance.J. Strength Cond. Res.243265–3269. 10.1519/JSC.0b013e3181b2c817

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  PosthumusM.SchwellnusM. P.CollinsM. (2011). The COL5A1 gene: a novel marker of endurance running performance.Med. Sci. Sports Exerc.43584–589. 10.1249/MSS.0b013e3181f34f4d

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  RiveraM. A.MartinezJ. L.CarrionA.FaheyT. D. (2011). AQP-1 association with body fluid loss in 10-km runners.Int. J. Sports Med.32229–233. 10.1055/s-0030-1268489

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  SachaJ. J.QuinnJ. M. (2011). The environment, the airway, and the athlete.Ann. Allergy Asthma Immunol.10681–87. 10.1016/j.anai.2010.06.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  Salas-SalvadoJ.RubioM. A.BarbanyM.MorenoB. (2007). SEEDO 2007 Consensus for the evaluation of overweight and obesity and the establishment of therapeutic intervention criteria.Med. Clin.128184–196.

  • Google Scholar

• 47

  SantiagoC.RuizJ. R.BuxensA.ArtiedaM.ArtetaD.Gonzalez-FreireM.et al (2011). Trp64Arg polymorphism in ADRB3 gene is associated with elite endurance performance.Br. J. Sports Med.45147–149. 10.1136/bjsm.2009.061366

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  SantosC. G.Pimentel-CoelhoP. M.BudowleB.de Moura-NetoR. S.Dornelas-RibeiroM.PompeuF. A.et al (2016). The heritable path of human physical performance: from single polymorphisms to the “next generation”.Scand. J. Med. Sci. Sports26600–612. 10.1111/sms.12503

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 49

  SaundersC. J.PosthumusM.O’ConnellK.SeptemberA. V.CollinsM. (2015). A variant within the AQP1 3’-untranslated region is associated with running performance, but not weight changes, during an Ironman Triathlon.J. Sports Sci.331342–1348. 10.1080/02640414.2014.989535

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  SchulerR.OsterhoffM. A.FrahnowT.MohligM.SprangerJ.StefanovskiD.et al (2017). Dietary fat intake modulates effects of a frequent ace gene variant on glucose tolerance with association to type 2 diabetes.Sci. Rep.7:9234. 10.1038/s41598-017-08300-7

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 51

  SheteA. N.ButeS. S.DeshmukhP. R. (2014). A Study of VO2 Max and body fat percentage in female athletes.J. Clin. Diagn. Res.8BC01–BC03. 10.7860/JCDR/2014/10896.5329

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 52

  StefanN.ThamerC.StaigerH.MachicaoF.MachannJ.SchickF.et al (2007). Genetic variations in PPARD and PPARGC1A determine mitochondrial function and change in aerobic physical fitness and insulin sensitivity during lifestyle intervention.J. Clin. Endocrinol. Metab.921827–1833. 10.1210/jc.2006-1785

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  TanakaH.MonahanK. D.SealsD. R. (2001). Age-predicted maximal heart rate revisited.J. Am. Coll. Cardiol.37153–156. 10.1016/s0735-1097(00)01054-8

  • CrossRef
  • Google Scholar

• 54

  WhitworthJ. A.ChalmersJ. (2004). World health organisation-international society of hypertension (WHO/ISH) hypertension guidelines.Clin. Exp. Hypertens.26747–752. 10.1081/ceh-200032152

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  WilsonJ. M.LoennekeJ. P.JoE.WilsonG. J.ZourdosM. C.KimJ. S. (2012). The effects of endurance, strength, and power training on muscle fiber type shifting.J. Strength Cond. Res.261724–1729. 10.1519/JSC.0b013e318234eb6f

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  WoodsD. R.HumphriesS. E.MontgomeryH. E. (2000). The ACE I/D polymorphism and human physical performance.Trends Endocrinol. Metab.11416–420. 10.1016/s1043-2760(00)00310-6

  • CrossRef
  • Google Scholar

• 57

  ZaniqueliD.MorraE. A.DantasE. M.BaldoM. P.CarlettiL.PerezA. J.et al (2014). Heart rate at 4 s after the onset of exercise in endurance-trained men.Can. J. Physiol. Pharmacol.92476–480. 10.1139/cjpp-2014-0042

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  Zilberman-SchapiraG.ChenJ.GersteinM. (2012). On sports and genes.Recent Pat. DNA Gene Seq.6180–188. 10.2174/187221512802717367

  • Pubmed Abstract
  • CrossRef
  • Google Scholar