Open Access

A genetic predisposition score for muscular endophenotypes predicts the increase in aerobic power after training: the CAREGENE study

  • Tom Thomaes1Email author,
  • Martine Thomis2,
  • Steven Onkelinx1,
  • Robert Fagard3,
  • Gert Matthijs4,
  • Roselien Buys1,
  • Dirk Schepers1,
  • Véronique Cornelissen1 and
  • Luc Vanhees1
BMC Genetics201112:84

DOI: 10.1186/1471-2156-12-84

Received: 15 March 2011

Accepted: 3 October 2011

Published: 3 October 2011



It is widely accepted that genetic variability might explain a large part of the observed heterogeneity in aerobic capacity and its response to training. Significant associations between polymorphisms of different genes with muscular strength, anaerobic phenotypes and body composition have been reported. Muscular endophenotypes are positively correlated with aerobic capacity, therefore, we tested the association of polymorphisms in twelve muscular related genes on aerobic capacity and its response to endurance training.


935 Coronary artery disease patients (CAD) who performed an incremental exercise test until exhaustion at baseline and after three months of training were included. Polymorphisms of the genes were detected using the invader assay. Genotype-phenotype association analyses were performed using ANCOVA. Different models for a genetic predisposition score (GPS) were constructed based on literature and own data and were related to baseline and response VO2 scores.


Carriers of the minor allele in the R23K polymorphism of the glucocorticoid receptor gene (GR) and the ciliary neurotrophic factor gene (CNTF) had a significantly higher increase in peakVO2 after training (p < 0.05). Carriers of the minor allele (C34T) in the adenosine monophosphate deaminase (AMPD1) gene had a significantly lower relative increase (p < 0.05) in peakVO2. GPS of data driven models were significantly associated with the increase in peakVO2 after training.


In CAD patients, suggestive associations were found in the GR, CNTF and the AMPD1 gene with an improved change in aerobic capacity after three months of training. Additionally data driven models with a genetic predisposition score (GPS) showed a significant predictive value for the increase in peakVO2.


Cardiac Rehabilitation Ischemic Heart Disease Exercise Test Physical capacity genetic association aerobic power


In patients with ischemic heart disease aerobic power and skeletal muscle strength are often reduced in comparison with age and sex-matched healthy subjects [1]. Regular physical activity improves aerobic power and skeletal muscle strength and is accompanied with an increase in survival in these patients. Moreover, aerobic power has been shown to be an independent prognostic factor of survival in patients with heart disease [24]. Therefore, participation in a cardiovascular exercise training program is strongly recommended in coronary artery disease (CAD) patients [57]. Recently, we reported an average increase in aerobic power of 26% ± 22% (Standard Deviation; SD) after a training period of three months in 1909 CAD patients [8]. A considerable interindividual variation could be observed with regard to aerobic power and its response to physical training. We identified 12 independent determinants, which accounted for 21% of the total variation in the aerobic power response to three months of training [8]. However, a large part of the individual variation remained unexplained. Evidence from healthy families and twin studies showed that a large part (50 to 66%) of the interindividual variation in aerobic power and its training response can be explained by genetic factors [911]. Likewise, multivariate genetic analysis on muscle strength phenotypes indicated that additive genetic factors could explain 66-78% of the variance in maximal elbow flexor torque [12]. Dynamic resistance training programs induce changes in muscle fibre size, fibre type and muscle metabolism [1316] which can lead to a greater oxygen uptake by the muscle and a higher peak oxygen uptake (peakVO2) [17]. Conley et al. [18] already pointed out that muscle volume is an important determinant in the response of peakVO2 to training. Therefore by increasing the amount of muscle mass used in the various exercises, independently of the mode of (resistance) training, the haemodynamic responses become more dynamic (aerobic) in nature. This was confirmed by Longhurst et al. [19] who attributed the slight increase of peakVO2 in weight lifters to the dynamic components in the exercises. In a subset of 80 CAD-patients we observed a significant increase (10-15%) in muscle force after three months of cardiac rehabilitation [20]. In agreement with Kostka et al. [21] we reported a correlation of r = 0.60 between baseline peakVO2 and isometric and isokinetic muscle strength. Further, a correlation of r = 0.32 between change in peakVO2 and change in isometric muscle force could be observed [20]. We therefore hypothesized that genes encoding for the different muscular subsystems, e.g. muscular structure, metabolism, cytokines, growth or differentiation factors, neurotrophic factors and hormones might be biologically plausible candidate genes that could be involved in aerobic power variability and its response to training. However, so far, research has mainly focussed on the associations between polymorphisms of these various genes and muscle strength, muscle function or their response to training and not on the association with aerobic power [22, 23].

We therefore aimed to identify whether previously described polymorphisms in these genes account for differences in aerobic power or its response to training in a cohort of 935 cardiac patients included in the CAREGENE study (CArdiac REhabilitation and GENetics of Exercise performance and training effect). The significance of a genetic predisposition score (GPS) for 7 literature based replicated strength-related polymorphisms was tested in predicting baseline aerobic power. Additionally, data driven models for GPS were tested in predicting baseline aerobic power and the response of peakVO2 to training.


Study population

A detailed description of the CAREGENE study design, eligibility criteria and methodology are presented elsewhere [24]. In brief, 935 biologically unrelated Caucasian CAD patients (76 women; mean age 56 ± 0.3 yrs) who had performed a graded cycle ergometer test until exhaustion, before and after three months of physical training, were included. The study protocol was approved by the Ethical Committee of the Faculty of Medicine of the Catholic University of Leuven and written informed consent was obtained from each participant.

Exercise training

Patients completed a three months ambulatory supervised exercise training program (90 minutes/session and involving cycling, running, arm ergometry, rowing, predominantly isotonic calisthenics and relaxation). Training frequency averaged 2.27 ± 0.02 times/week and training intensity was 79.7 ± 0.35% of the maximal heart rate (HRmax). The latter was calculated as: (training heart rate/peak heart rate) * 100, where the mean training heart rate of the last 3 exercise sessions and peak heart rate of the exercise test after training were used.

Exercise testing

A detailed description of the exercise test methodology has been reported previously [24]. Before and after the three months training program patients completed a maximal graded exercise test on a cycle ergometer with respiratory gas analysis. The initial workload of 20 W was increased until exhaustion by 30 W every 3 min and after a protocol change in the year 2000 the initial workload of 20 W was increased by 20 W per min. Heart rate was measured throughout the test with a 12-lead ECG and blood pressure was measured every 2 minutes. Respiratory data were measured through breath-by-breath analysis. PeakVO2 was defined as the highest 15 seconds average of VO2 obtained at the end of the test and was expressed as mL.min-1. The percent of predicted peakVO2 was calculated as peakVO2 divided by maximal predicted VO2, using the values reported by Wasserman et al. [25].

Genotype determinations

21 single nucleotide polymorphisms (SNPs) in 12 muscular related genes were genotyped (Table 1). Selection of genes and SNPs was based on a literature search at the beginning of the study. SNPs had to be associated with muscular phenotypes or aerobic phenotypes. DNA was extracted from white blood cells using the 'salting-out' method and the Invader TM assay (Third Wave Technologies) was used for genotyping [26]. The Invader assay combines structure-specific cleavage enzymes and a universal fluorescent resonance energy transfer (FRET) system. When the gene-specific probes bind the target, these enzymes will cleave. This mechanism warrants the specificity for distinguishing between alleles, whereas the FRET system generates an amplified readout [27]. In every experiment synthetic target oligonucleotides were used as controls. Third Wave Technologies designed and provided all reaction components. Genotyping was performed in a 96-well format. By combining Probe mix (211 μl), Cleavase enzyme (79 μl) FRET mix (317 μl) and MAP buffer (26 μl) the reaction mixture was prepared. Of this mixture 6 mL was added into a 96-well plate. Six micro litres of no target blank, synthetic target oligonucleotides, or genomic DNA samples (50 ng/μl) were added. After short centrifugation and incubation at 63°C for four hours, fluorescent intensities were measured using a fluorescence microtiter plate reader. By calculating the ratios of the net wild type and net mutant (or minor allele) signals [28], genotypes were determined. The analysis was repeated once on those samples where no genotype could be obtained during initial testing. Samples were excluded from the analysis when genotyping failed twice (this failure rate varies between 1.4% and 5.1%, data not shown).
Table 1

Overview of the investigated genes and polymorphisms within muscular structural or functional subsystems




Genotype frequency (n patients)


adenosine monophosphate deaminase (AMPD1)


CC (652)

CT (239)

TT (24)

Muscular structure

Alpha-actinin 3 (ACTN3)


TT (303)

TC (444)

CC (176)



GG (468)

GA (386)

AA (66)


Myosine light chain kinase


GG (589)

AG (289)

AA (38)




GG (749)

GT (157)

TT (9)


Interleukin 6 (IL-6)


GG (332)

GC (411)

CC (176)


Interleukin 15 receptor


CC (238)

CA (468)

AA (206)


alpha (IL-15Ra)


TT (223)

GT (480)

GG (211)



GG (227)

GA (459)

AA (203)

Growth or differentiation factors

Insulin-like growth factor 2 (IGF-II)


GG (448)

AG (388)

AA (80)


Activin-type II receptor B (ACVR2B)


TT (283)

CT (468)

CC (162)


Follistatin (FST)


CC (605)

CT (278)

TT (30)



CC (606)

CT (278)

TT (35)



TT (601)

CT (275)

CC (36)


v-akt murine thymoma viral oncogene homolog 1 (AKT1)


GG (452)

GT (377)

TT (76)

Neurotrophic factors

Ciliary Neurotrophic Factor (CNTF)


GG (664)

GA (226)

AA (21)


Ciliary Neurotrophic


CC (580)

CT (300)

TT (41)


Factor Receptor (CNTFR)


CC (695)

CT (187)

TT (12)


Glucocorticoid Receptor


GG (857)

GA (54)

AA (1)



AA (843)

AG (70)

GG (2)



CC (399)

GC (403)

GG (120)

Genotype frequencies (N) in the CAREGENE study included in brackets

Statistical methods

Data were analyzed using SAS statistical software ® version 9.1 for Windows (SAS Institute Inc, Cary, NC, USA). Data are reported as means ± SE or as number of patients with percentage for dichotomous variables. To test whether the observed genotype frequencies were in Hardy-Weinberg equilibrium a χ2- test with one degree of freedom was used. Distributions were checked for normality with the Shapiro-Wilk statistic. Comparisons between the exercise test at baseline and after training were made by paired Student's t-test; comparisons across genotypes by AN(C)OVA, followed by Fisher's protected least significance difference (LSD) method if significant. Where appropriate, categorical data were tested by χ2 or by Fisher's exact test. ANCOVA was used with age, gender, stature, body mass and baseline peakVO2 as covariates to test potential relationships between gene variation and the increase in peakVO2 after training. Gene polymorphisms with a minor allele frequency of less than 2% were analysed as carriers vs. non-carriers of the minor allele.

Given the large set of gene variants and 3 phenotypes under study, multiple testing issues probably induce false positive interpretation of association findings. To partially overcome the multiple testing problem and to test for the additive effect of cumulative (response in) VO2 uptake- alleles, a genetic predisposition score (GPS) was calculated and its predictive value tested. Based on recent literature, seven polymorphisms were selected for this 'increasing allele' approach for baseline peakVO2 [2931]. Each of these polymorphisms has been associated with muscle strength or function in at least three independent replicated studies with evidence for the same increasing allele in all studies. The selected polymorphisms were Alpha-actinin 3 (ACTN3 R577X), Myosine light chain kinase (MLCK C49T), Insulin-like growth factor 2 (IGF-II ApaI), Glucocorticoid receptor (GR, R23K), Interleukin 15 receptor alpha (IL15Ra BstnI), Ciliary Neurotrophic Factor Receptor (CNTFR C-1703T and C174T). Additional file 1 shows which alleles were chosen as the 'increasing alleles'. An additive effect was hypothesized and equal weights were given for each increasing allele, because no well-defined effect sizes are known for the different SNPs and weighting of increasing alleles might only have limited effect [32]. The genetic predisposition scores was calculated for each individual by adding the strength-increasing alleles (0, 1 or 2) from all seven variants and was analyzed in a GLM approach with age, gender, stature and body mass as covariates. Additionally, based on our own data, backward regression analysis was applied to detect three subsets of SNPs to be associated with the three phenotypes; baseline peakVO2, the simple difference (ml/min) in peakVO2 (ml/min) and the relative difference (%) in peakVO2 after training. Based on these subsets of variants we created three models to test the predictive value of its GPS. Finally, when the predictive value of the GPS was significant in the model, patients were categorized in the upper 10% peakVO2 (High group) and lower 10% peakVO2 (Low group) using sex-specific cut-off scores. Similarly, high and low responder groups were defined. Different logistic regression models were used to test the power of the GPS to predict 'high group' status, and to produce the receiver operating characteristic curve (ROC curve) and the area under the receiver operating characteristic curve (AUC). Predicted probability was calculated for a hypothetical average patient. Adjusted odds ratios were calculated to estimate the effect of the GPS score as the odds per increasing allele to belong to the high peakVO2 group or to the high responder group. All statistical tests were performed two-sided at a significance level of 5%.


Descriptive characteristics of the overall study cohort (n = 935) are presented in Table 2. At baseline, overall aerobic power averaged 1716 ± 16 mL.min-1 or 77.9 ± 0.6% (± SE) of the predicted normal value (men: 1766 ± 16 ml/min, women: 1142 ± 31 ml/min, p < 0.0001). Three months of physical training resulted in a significant average increase of the aerobic power by 24.2 ± 0.6% (p < 0.001) (men: 24.1 ± 17.4%, women: 25.1 ± 21.1, p = 0.47), ranging from a decrease of 33.6% to an increase of 111.1%. Baseline peakVO2 and its change after training follow the normal Gaussian distribution (p > 0.05). At peak exercise, respiratory gas exchange ratio averaged 1.14 ± 0.003 at baseline and 1.13 ± 0.002 after training (p > 0.05). Peak ventilatory equivalent for oxygen was 38 ± 0.2 and 37 ± 0.2 respectively (p > 0.05). None of the patient characteristics and training characteristics were significantly different amongst the different polymorphisms. χ2 tests revealed no deviations from the Hardy-Weinberg equilibrium except for G-174C of the IL-6 gene (p = 0.017).
Table 2

Clinical characteristics for biologically unrelated Caucasian CAD patients (n = 935) in the CAREGENE study


Overall Cohort


N (%)


76 (8)

Age (years)

56 ± 0.3

Body mass index (kg.m-2)

25.8 ± 0.1

History of diabetes

49 (5)

History of hypertension

251 (27)

Current smoking

45 (5)

Past smoking

681 (73)

Complaints of angina pectoris in daily life

41 (4)

dyspnoea in daily life

149 (16)


630 (67)


252 (27)


333 (36)


377 (40)


470 (50)

Angina pectoris

23 (2)

AMI, acute myocardial infarction; CBS, coronary bypass surgery; PCI, percutaneous coronary Intervention. Some patients had more than one pathology. Data are presented as means ± SE for continuous variables and as numbers (percentage) for dichotomous variables.

Single polymorphism associations

Significant associations with the simple difference (ml/min) and/or relative difference (%) in aerobic capacity after three months of training could be observed for 3 out of the 21 investigated polymorphisms which were genotyped (Table 3 and Additional file 2, Table S1).
Table 3

Genotype-phenotype association analysis for muscular subsystem gene polymorphisms and baseline aerobic power and changes after training in the CAREGENE study












N (%)





(%) b





652 (71)

1711 ± 15

p = 0.40

396 ± 10

p = 0.11

24.8 ± 0.6

p = 0.04



263 (29)

1734 ± 23


367 ± 16


22.4 ± 0.99





664 (73)

1724 ± 15

p = 0.43

369 ± 9

p = 0.001

23.1 ± 0.6

p = 0.002



226 (25)

1701 ± 25


416 ± 16


25.1 ± 1.0



21 (2)

1801 ± 81


519 ± 52


33.9 ± 3.3





857 (94)

1715 ± 13

p = 0.08

383 ± 9

p = 0.02

23.8 ± 0.6

p = 0.04



55 (6)

1805 ± 50


464 ± 34*


28.5 ± 2.2*


AMPD1, Adenosine monophosphate deaminase; CNTF, Ciliary neurotrophic factor; GR, Glucocorticoid receptor

a. Mean ± SE, corrected for gender, age, height and weight;

b. Mean ± SE, corrected for gender, age, height, weight and baseline peakVO2

In the R23K polymorphisms of the GR gene, carriers of the minor allele had a significantly larger increase in peakVO2, both in simple difference (ml/min) (p < 0.05) and in relative (p < 0.05) term, compared to the homozygous genotype groups. Further, homozygous CC genotype group members for the C34T polymorphism of the AMPD1 gene showed a significantly larger relative increase in peakVO2 compared to the carriers of the minor allele (p < 0.05), whereas the increase in peakVO2 was significantly larger in the AA genotype group of the G-6A polymorphism of the CNTF gene compared to the homozygous GG group (p < 0.01).

Genetic predisposition score

A theoretical maximal GPS of 14 increasing alleles could occur for seven selected polymorphisms based on recent literature. 848 patients had complete data for all of these polymorphisms. Figure 1 shows the distribution of GPS within our patient group with four patients exhibiting only one increasing allele up to six patients having nine increasing alleles (max GPS = 9). To avoid the possibility of false positive results by the smaller groups at the two tails, we combined groups 1 and 2 and groups 8 and 9. ANCOVA analysis showed no overall significant differences between the patient groups categorized by GPS scores and baseline peakVO2 (Figure 1) or simple and relative difference in peakVO2 after training. However, based on linear regression (b = 13.5, p = 0.06) analysis, we observed a trend that higher numbers of increasing alleles (higher GPS) result in a higher baseline peakVO2.
Figure 1

Genetic predisposition score for muscular subsystems gene polymorphisms (literature based) and baseline peakVO 2 in the CAREGENE study. Left Y-axis: Number of patients in each increasing alleles group (bar graph). Right Y-axis: Baseline peakVO2 for each increasing alleles group (square dots ± SE) corrected for age, gender, height and body mass. X-axis: GPS - Number of increasing alleles. Regression line for baseline peakVO2

Backward regression analyses revealed three subsets of genes to be associated with baseline peakVO2 (ACTN3, Q523R and GR, R23K), simple difference in peakVO2 (ACVR2B, rs2268757; CNTF, G-6A and GR, R23K) and relative difference in peakVO2 (AMPD1, C34T; ACVR2B, rs2268757; CNTF, G-6A and GR, R23K). Significant differences were found in the models for simple (p < 0.01) and relative difference (p < 0.001) in peakVO2 (Figure 2 and 3).
Figure 2

Genetic predisposition score for 3 polymorphisms associated with the change (ml/min) in peakVO 2 in the CAREGENE study. Left Y-axis: Number of patients in each increasing alleles group (bar graph). Right Y-axis: Change in peakVO2 (ml/min) for each increasing alleles group (square dots ± SE) corrected for age, gender, height, body mass, and baseline peakVO2. X-axis: GPS - Number of increasing alleles. Regression line for the change in peakVO2 (ml/min)
Figure 3

Genetic predisposition score for 4 polymorphisms associated with relative change (%) in peakVO 2 in the CAREGENE study. Left Y-axis: Number of patients in each increasing alleles group (bar graph). Right Y-axis: Relative change in peakVO2 (%) for each increasing alleles group (square dots ± SE) corrected for age, gender, height, body mass, and baseline peakVO2. X-axis: GPS - Number of increasing alleles. Regression line for the relative change in peakVO2

High vs. Low aerobic capacity groups

Additional analyses of the high vs. low responder groups were performed for the GPS of the simple and relative difference after training. In the model for simple difference in peakVO2, GPS alone was statistically significant to predict a high increase in peakVO2 (AUC: 0.62; 95%CI: 0.54-0.69, p < 0.01), (Odds ratio: 1.62; 95%CI: 1.13-2.33) (Figure 4). In the model for relative increase in peakVO2, GPS alone was also statistically significant to predict a high relative increase in peakVO2 (AUC: 0.63; 95%CI: 0.56-0.71, p < 0.01), (Odds ratio: 1.68; 95%CI: 1.22-2.31) (Figure 5). Both models were comparable to a model with age and weight to predict the simple and relative difference in peakVO2 (AUC of 0.69 and 0.60 respectively). When GPS was added to the model with age and weight, the predictive value increased with 2% (simple difference) and 5% (relative difference). However this increase in predictive value was not significantly different from the model with GPS alone or the model with age and weight. Height and gender had no significant contribution to the model and were omitted. A patient with a high GPS has a higher probability to end up in the group with the 10% highest increase in peakVO2 values after training and vice versa. The GPS had however no predictive value for the baseline peakVO2 for both the literature based model as the own-data model.
Figure 4

Overlay of three different models to predict high vs. low responder in peakVO 2 change (ml/min) after training. Model GPS (AUC: 0.62; 95% CI: 0.54 - 0.69). Model age and weight (AUC: 0.69; 95% CI: 0.61 - 0.77) Model GPS, age and weight (AUC: 0.71; 95% CI: 0.63 - 0.78)
Figure 5

Overlay of three different models to predict high vs. low responder in relative peakVO 2 change (%) after training. Model GPS (AUC: 0.63; 95% CI: 0.56 - 0.71). Model age and weight (AUC: 0.60; 95% CI: 0.52 - 0.69) Model GPS, age and weight (AUC: 0.65; 95% CI: 0.57 - 0.73)


Since dynamic resistance training can induce changes in muscle size, fiber type and muscle metabolism [1316] and can lead to a greater oxygen uptake, we hypothesized that genes coding for these muscular subsystems might also be involved in aerobic power or its change after training through the different pathways of muscular endophenotypes (i.e. structure, growth factors, hormones, ...). We found suggestive associations in polymorphisms of the GR gene, AMPD1 gene and CNTF gene with aerobic capacity or its change after training in CAD patients.

Single polymorphism associations


The mutation in codon 23 results in a change of arginine to lysine. The mutation in R23K polymorphism, which is a part of the ER22/23EK polymorphism complex (as often described in the literature), has been associated with a decreased glucocorticoid sensitivity which can lead to decreased catabolism and less muscle atrophy [3335]. Van Rossum et al. [36] showed that male ER22/23EK-carriers had a higher muscle mass and strength. In concordance with these results we found that carriers of the minor allele of the R23K polymorphism had a significantly higher increase in aerobic power after three months of training. Although no significant differences in peakVO2 were apparent at baseline, a tendency could be observed that carriers of the minor allele had a higher baseline peakVO2 (p = 0.08).


The C34T mutation results in a premature stop codon of protein synthesis which is the main cause of AMPD deficiency. Fishbein et al. [37] proposed that a deficiency of AMPD causes muscular weakness or cramping after exercise. Rico-Sanz et al. [38] showed that individuals homozygous for the T allele in the C34T polymorphism of the AMPD1 gene have diminished exercise capacity and cardiorespiratory responses to exercise in the sedentary state and that the training response of ventilatory phenotypes during maximal exercise is more limited in TT. Due to the small number of individuals in the TT homozygous group, we performed a combined CT + TT vs. CC analysis. We showed that T allele carriers had a significantly lower relative increase in peakVO2 after three months of training compared to individuals homozygous for the C allele (p < 0.05) which agrees with earlier findings that the C34T mutation can lead to a decreased exercise capacity and a lower response to exercise (Rico-Sanz 2003).


Circulating CNTF has trophic effects on neuronal and muscular tissues. Earlier De Mars et al. [39] and Roth et al. [40] hypothesized that, with regard to the CNTF G-6A polymorphism, carriers of the minor allele would have lower muscle strength than wild-type carriers due to lower levels of circulating CNTF. However, both research groups were unable to confirm their hypothesis. Subsequently, Roth et al. [40] showed that individuals heterozygous for the G-6A polymorphism possess significantly greater muscle strength and muscle quality at relatively fast contraction speeds compared to homozygous wild-type individuals. Therefore it is plausible that improved muscle strength and quality in these individuals may also result in a more pronounced increase in aerobic power. We showed that A allele carriers of the G-6A polymorphism had a larger increase in aerobic capacity after three months of training (p < 0.01). The homozygous minor allele group had almost 11% additional increase in peakVO2 after three months of training compared with the homozygous wild-type group. To our knowledge, aerobic capacity has never been studied in the G-6A polymorphisms of the CNTF gene, so it might be that individuals homozygous for the A allele who exhibit higher muscle strength and quality at fast contraction speeds [40] also have a higher aerobic capacity increase, measured on a bicycle treadmill. Unfortunately we do not have muscle strength measurements in this patient group so this hypothesis remains to be confirmed in additional studies.

Genetic predisposition score

To partially overcome the multiple testing problem and to test for the additive effect of cumulative (response in) VO2 uptake- alleles, an 'increasing allele analysis' is performed, as was proposed by studies of Ruiz et al. [29], Santiago et al. [30] and Williams et al. [31]. In calculating a genetic predisposition score (GPS) every polymorphism is considered to have an equal influence on the phenotype, although weighting approaches have also been taken [32]. Based on this hypothesis, every additional occurrence of a 'beneficial' allele for muscular strength would result in an increased peakVO2. We constructed an 'increasing allele analysis' for baseline aerobic capacity with 7 polymorphisms, based on the literature. Although no overall significant differences between different allele groups were found, based on linear regression (b = 13.5, p = 0.06) analysis, we observed a trend that a higher GPS results in a higher baseline peakVO2. An increase of peakVO2 by 13.5 ml/min per allele equals approximately 1% increase in peakVO2 with every additional increasing allele.

Based on backward regression analysis of the complete group of polymorphisms, we identified three significant subsets of polymorphisms associated with the three phenotypes; baseline peakVO2, simple and relative difference in peakVO2. The model for relative change showed remarkable resemblance with the single polymorphism association analysis (only rs2268757 of the ACVR2B gene was not significantly associated in the single polymorphism associations). The GPS models for simple and relative difference in peakVO2 showed significant associations with the change in peakVO2 after training. The GPS for the simple difference showed that patients with a GPS of 0 had an increase of 335.3 ± 18.0 ml/min whereas patients with a GPS of 2 or more would increase with approximately 420 ml/min. The GPS for the relative change showed similar results. GPS of 1 or less had an increase of 16.7% whereas patients with GPS of 4 or more increased their peakVO2 with approximately 25%. The findings in both models were confirmed with an analysis of the 10% highest responders versus the 10% lowest responders for both simple difference (ml/min) and relative difference (%) in peakVO2 after training. The odds ratio estimate for both models showed that with every additional increasing allele, the chance of belonging to the high responder peakVO2 group is 62% and 68% higher for simple and relative response respectively.

In conclusion, although we found some suggestive associations between SNPs with aerobic capacity or its improvement after three months of training, one should carefully interpret the results. The CAREGENE study population is not a random sample of the total population, as they are all Caucasian CAD patients who had been hospitalized before baseline testing (± 8 weeks post-operation). These CAD patients who probably were deconditioned by their disease have proceeded through a spontaneous recovery which might also differ individually. Generalization is therefore limited. It might be that our individuals are also genetically predisposed for CAD and that this patient group might have distorted allele frequencies for the genes under study. However allele frequencies did not differ to the reported CEPH frequencies in the HapMap database. One could speculate that if the identified gene variants are not specifically related to the CAD-status of the patients, and if underlying pathways that drive the adaptations do not differ between CAD and healthy individuals, that these might also play a role in exercise-induced changes in aerobic capacity in healthy individuals.

Literature search at the beginning of the study led to a selection of genes and polymorphisms in which associations were found with muscle structure or muscle function. We are aware that the field is evolving quickly and 'new' candidate gene variants have become available in the mean time. We hypothesized that if patients had better muscular performance by mutations of their genes they might also have a better aerobic capacity or the ability to improve more beneficially. This could be generated by a better local oxygen extraction by the muscles, which leads to a higher peakVO2. In a substudy on 80 patients we found a correlation of r = 0.60 between baseline quadriceps isometric and isokinetic muscle force and baseline peakVO2 and a correlation of r = 0.32 between change in peakVO2 and change in isometric muscle force [20]. Although these phenotypic correlations have not been decomposed into genetic versus environmental covariation between the endurance and strength phenotypes, they are however indicative for at least partial pleiotrophic gene action. For more closely related phenotypes e.g. walking speed and lower leg power, substantial genetic correlations have been reported [41]. Although we describe some pathways within our associations, the lack of an objective measurement of muscle force or function makes it hard to confirm these hypotheses. In addition one could also question if these genes for muscular structure and function are 'good endophenotypes' for improved peakVO2 as some studies have indicated that the 'increasing allele' for endurance is not necessarily the same as for maximal strength. With regard to genetic influence on aerobic capacity and/or its improvement after training, more research is definitively needed with additional muscular testing.

Furthermore we set the significance level at p < 0.05, however when applying Bonferroni correction for multiple testing none of our results would reach significance. This study therefore also included a genotypic predisposition score analysis, in which additive effects on seven polymorphisms were studied, as well as a data-driven GPS analysis was performed. Although meaningful phenotypic correlations between muscle strength and endurance performance in this specific population do suggest underlying pleiotrophic gene actions, the contribution of each specific genetic variant is expected to be small.


Suggestive associations were found for changes in aerobic capacity and single polymorphisms of the CNTF, AMPD1 and GR gene and a genetic predisposition score showed a significant predictive value for this GPS to estimate an individuals' potential to belong to the higher peakVO2 responder group, or inversely to predict the risk of belonging to the low responder peakVO2 group. Although no strong associations were found, it is reasonable to assume that mutations in genes coding for muscular structure and function can lead to changes in aerobic capacity.



This study was supported by grants from the Fund for Scientific Research -Flanders 'Fonds voor Wetenschappelijk Onderzoek - Vlaanderen', Belgium (F.W.O. grant G.0624.08 and G.0124.02) and from the Research Council of the University of Leuven 'Onderzoeksraad K.U.Leuven', Belgium (grant OT/07/064 and OT/01/046).

Authors’ Affiliations

Cardiovascular Rehabilitation Unit, Department of Rehabilitation Sciences, Katholieke Universiteit Leuven
Research Centre for Exercise and Health, Department of Biomedical Kinesiology, Katholieke Universiteit Leuven
Hypertension and Cardiovascular Rehabilitation Unit, Department of Cardiovascular Diseases, UZ Leuven
Department of Human Genetics, Centre for Human Genetics of the University Hospitals, UZ Leuven


  1. Ghroubi S, Chaari M, Elleuch H, Massmoudi K, Abdenadher M, Trabelssi I, Akrout M, Feki H, Frikha I, Dammak J, Kammoun S, Zouari N, Elleuch MH: The isokinetic assessment of peripheral muscle function in patients with coronary artery disease: correlations with cardiorespiratory capacity. Ann Readapt Med Phys. 2007, 50 (5): 295-301.View ArticlePubMed
  2. Kavanagh T, Mertens DJ, Hamm LF, Beyene J, Kennedy J, Corey P, Shephard RJ: Prediction of long-term prognosis in 12 169 men referred for cardiac rehabilitation. Circ. 2002, 106: 666-671. 10.1161/01.CIR.0000024413.15949.ED.View Article
  3. Myers J, Prakash M, Froelicher V, Do D, Partington S, Atwood JE: Exercise capacity and mortality among men referred for exercise testing. N Engl J Med. 2002, 346: 793-801. 10.1056/NEJMoa011858.View ArticlePubMed
  4. Vanhees L, Fagard R, Thijs L, Staessen J, Amery A: Prognostic significance of peak exercise capacity in patients with coronary artery disease. J Am Coll Cardio. 1994, 23: 358-363. 10.1016/0735-1097(94)90420-0.View Article
  5. Borjesson M, Assanelli D, Carre F, Dugmore D, Panhuyzen-Goedkoop NM, Seiler C, Senden J, Solberg EE: ESC Study Group of Sports Cardiology: recommendations for participation in leisure-time physical activity and competitive sports for patients with ischaemic heart disease. Eur J Cardiovasc Prev Rehabil. 2006, 13: 137-149. 10.1097/01.hjr.0000199494.46708.5a.View ArticlePubMed
  6. Giannuzzi P, Saner H, Bjornstad H, Fioretti P, Mendes M, Cohen-Solal A, Dugmore L, Hambrecht R, Hellemans I, McGee H, Perk J, Vanhees L, Veress G: Secondary prevention through cardiac rehabilitation: position paper of the Working Group on Cardiac Rehabilitation and Exercise Physiology of the European Society of Cardiology. Eur Heart J. 2003, 24: 1273-1278. 10.1016/S0195-668X(03)00198-2.View ArticlePubMed
  7. Mezzani A, Giannuzzi P: Physical activity for cardiovascular disease prevention. Ital Heart J. 2003, 4: 739-744.PubMed
  8. Vanhees L, Stevens A, Schepers D, Defoor J, Rademakers F, Fagard R: Determinants of the effects of physical training and of the complications requiring resuscitation during exercise in patients with cardiovascular disease. Eur J Cardiovasc Prev Rehabil. 2004, 11: 304-312. 10.1097/01.hjr.0000136458.44614.a2.View ArticlePubMed
  9. Bouchard C, Daw EW, Rice T, Perusse L, Gagnon J, Province MA, Leon AS, Rao DC, Skinner JS, Wilmore JH: Familial resemblance for VO2max in the sedentary state: the HERITAGE family study. Med Sci Sports Exerc. 1998, 30: 252-258.View ArticlePubMed
  10. Fagard R, Bielen E, Amery A: Heritability of aerobic power and anaerobic energy generation during exercise. J Appl Physiol. 1991, 70: 357-362.PubMed
  11. Prud'homme D, Bouchard C, Leblanc C, Landry F, Fontaine E: Sensitivity of maximal aerobic power to training is genotype-dependent. Med Sci Sports Exerc. 1984, 16: 489-493. 10.1249/00005768-198410000-00012.View ArticlePubMed
  12. Thomis MA, Van Leemputte M, Maes HH, Blimkie CJR, Claessens AL, Marchal G, Willems E, Vlietinck RF, Beunen GP: Multivariate genetic analysis of maximal isometric muscle force at different elbow angles. J Appl Physiol. 1997, 82: 959-967.PubMed
  13. Blazevich AJ, Cannevan D, Coleman DR, Horne S: Influence of concentric and eccentric resistance training on architectural adaptation in human quadriceps muscles. J Appl Physiol. 2007, 103: 1565-1575. 10.1152/japplphysiol.00578.2007.View ArticlePubMed
  14. Blazevich AJ, Gill ND, Bronks R, Newton R: Training-Specific Muscle Architecture Adaptation after 5-wk Training in Athletes. Med Sci Sports Exerc. 2003, 35: 2013-2022. 10.1249/01.MSS.0000099092.83611.20.View ArticlePubMed
  15. Ogawa T, Spina RJ, Martin WH, Kohrt WM, Schechtman KB, Holloszy JO, Ehsani AA: Effects of aging, sex and physical training on cardiovascular responses to exercise. Circ. 1992, 86: 494-503.View Article
  16. Spina RJ, Ogawa T, Kohrt WM, Martin WH, Holloszy JO, Ehsani AA: Differences in cardiovascular adaptations to endurance exercise training between older men and women. J Appl Physiol. 1993, 75: 849-55.PubMed
  17. Cornelissen VA, Fagard RH: Effect of resistance training on resting blood pressure: a meta-analysis of randomized controlled trials. J Hypertens. 2005, 23: 251-259. 10.1097/00004872-200502000-00003.View ArticlePubMed
  18. Conley KE, Esselman PC, Jubrias SA, Cress ME, Inglin B, Mogadam C, Schoene RB: Ageing, muscle properties and maximal O(2) uptake rate in humans. J Physiol. 2000, 526: 211-217. 10.1111/j.1469-7793.2000.00211.x.PubMed CentralView ArticlePubMed
  19. Longhurst JC, Stebbins CL: The isometric athlete. Cardiol Clin. 1992, 10: 281-294.PubMed
  20. Cornelissen VA, Defoor JGM, Stevens A, Schepers D, Hespel P, Decramer M, Mortelmans L, Dobbels F, Vanhaecke J, Fagard R, Vanhees L: Effect of Creatine Supplementation as a Potential Adjuvant Therapy to Exercise Training in Cardiac Patients: a randomized controlled trial. Clin Rehabil. 2010, 24: 988-999. 10.1177/0269215510367995.View ArticlePubMed
  21. Kostka T, Rahmani A, Berthouze SE, Lacour J-R, Bonnefoy M: Quadriceps muscle function in relation to habitual physical activity and VO2max in men and women aged more than 65 Years. J Gerontol A Biol Sci Med Sci. 2000, 55: B481-488. 10.1093/gerona/55.10.B481.View ArticlePubMed
  22. Bray MS, Hagberg JM, Pérusse L, Rankinen T, Roth SM, Wolfarth B, Bouchard C: The human gene map for performance and health-related fitness phenotypes: the 2006-2007 update. Med Sci Sports Exerc. 2009, 41: 35-73.View ArticlePubMed
  23. Stewart CEH, Rittweger J: Adaptive processes in skeletal muscle: Molecular regulators and genetic influences. J Musculoskelet Neuronal Interact. 2006, 6: 73-86.PubMed
  24. Defoor J, Martens K, Zielinska D, Matthijs G, Van Nerum H, Schepers D, Fagard R, Vanhees L: The CAREGENE study: polymorphisms of the beta1-adrenoceptor gene and aerobic power in coronary artery disease. Eur Heart J. 2006, 27: 808-816. 10.1093/eurheartj/ehi737.View ArticlePubMed
  25. Wasserman K, Hansen J, Sue D, Whipp B: Physiology of exercise. Principles of Exercise Testing and Interpretation. Edited by: Whipp B. 1999, Philadelphia: Lea & Febiger, 10-61.
  26. Miller SA, Dykes DD, Polesky HF: A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988, 16: 1215-10.1093/nar/16.3.1215.PubMed CentralView ArticlePubMed
  27. De AM, Lyamichev VI, Eis PS, Iszczyszyn W, Kwiatkowski RW, Law SM, Olson MC, Rasmussen EB: Invader technology for DNA and RNA analysis: principles and applications. Expert Rev Mol Diagn. 2002, 2: 487-496. 10.1586/14737159.2.5.487.View Article
  28. Hessner MJ, Budish MA, Friedman KD: Genotyping of factor V G1691A (Leiden) without the use of PCR by invasive cleavage of oligonucleotide probes. Clin Chem. 2000, 46: 1051-1056.PubMed
  29. Ruiz JR, Gómez-Gallego F, Santiago C, González-Freire M, Verde Z, Foster C, Lucia A: Is there an optimum endurance polygenic profile?. J Physiol. 2009, 587: 1527-1534. 10.1113/jphysiol.2008.166645.PubMed CentralView ArticlePubMed
  30. Santiago C, Ruiz JR, Muniesa CA, González-Freire M, Gómez-Gallego F, Lucia A: Does the polygenic profile determine the potential for becoming a world-class athlete? Insights from the sport of rowing. Scand J Med Sci Sports. 2009, 20: e188-e194.View ArticlePubMed
  31. Williams AG, Folland JP: Similarity of polygenic profiles limits the potential for elite human physical performance. J Physiol. 2008, 586: 113-121.PubMed CentralView ArticlePubMed
  32. Janssens AC, Moonesinghe R, Yang Q, Steyerberg EW, van Duijn CM, Khoury MJ: The impact of genotype frequencies on the clinical validity of genomic profiling for predicting common chronic diseases. Genet Med. 2007, 9: 528-35. 10.1097/GIM.0b013e31812eece0.View ArticlePubMed
  33. Czerwinski SM, Hickson RC: Glucocorticoid receptor activation during exercise in muscle. J Appl Physiol. 1990, 68: 1615-1620.PubMed
  34. Menconi M, Fareed M, O'Neal P, Poylin V, Wei W, Hasselgren PO: Role of glucocorticoids in the molecular regulation of muscle wasting. Crit Care Med. 2007, 35 (Suppl 9): S602-608.View ArticlePubMed
  35. Peeters GM, van Schoor NM, van Rossum EF, Visser M, Lips P: The relationship between cortisol, muscle mass and muscle strength in older persons and the role of genetic variations in the glucocorticoid receptor. Clin Endocrinol (Oxf). 2008, 69: 673-82. 10.1111/j.1365-2265.2008.03212.x.View Article
  36. Van Rossum EF, Voorhoeve PG, te Velde SJ, Koper JW, Delemarre-van de Waal HA, Kemper HC, Lamberts SW: The ER22/23EK polymorphism in the glucocorticoid receptor gene is associated with a beneficial body composition and muscle strength in young edults. J Clin Endocrinol Metab. 2004, 89: 4004-4009. 10.1210/jc.2003-031422.View ArticlePubMed
  37. Fishbein WN, Armbrustmacher VW, Griffin JL: Myoadenylate deaminase deficiency: a new disease of muscle. Science. 1978, 200: 545-548. 10.1126/science.644316.View ArticlePubMed
  38. Rico-Sanz J, Rankinen T, Joanisse DR, Leon AS, Skinner JS, Wilmore JH, Rao DC, Bouchard C: Associations between cardiorespiratory responses to exercise and the C34T AMPD1 gene polymorphism in the HERITAGE Family study. Physiol Genomics. 2003, 14: 161-166.View ArticlePubMed
  39. De Mars G, Windelinckx A, Beunen G, Delecluse C, Lefevre J, Thomis MA: Polymorphisms in the CNTF and CNTF receptor genes are associated with muscle strength in men and women. J Appl Physiol. 2007, 102: 1824-1831. 10.1152/japplphysiol.00692.2006.View ArticlePubMed
  40. Roth SM, Schrager MA, Ferrell RE, Riechman SE, Metter EJ, Lynch NA, Lindle RS, Hurley BF: CNTF genotype is associated with muscular strength and quality in humans across the adult age span. J Appl Physiol. 2001, 90: 1205-1210.PubMed
  41. Tiainen K, Pajala S, Sipilä S, Kaprio J, Koskenvuo M, Alén M, Heikkinen E, Tolvanen A, Rantanen T: Genetic effects in common on maximal walking speed and muscle performance in older women. Scand J Med Sci Sports. 2007, 17: 274-280.PubMed


© Thomaes et al; licensee BioMed Central Ltd. 2011

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.