|
|
||||||||
1Pennington Biomedical Research Center, Human Genomics Laboratory, Louisiana State University System, Baton Rouge, Louisiana; 2Merikoski Rehabilitation and Research Center, Oulu; 3Division of Cardiology, Department of Medicine, University of Oulu, Oulu; and 4Lapland Central Hospital, Rovaniemi, Finland
Submitted 11 November 2005 ; accepted in final form 19 February 2006
| ABSTRACT |
|---|
|
|
|---|
genotype; exercise training; cardiovascular autonomic function
Regular endurance training shifts the cardiac autonomic balance toward vagal dominance (3, 7, 13, 17, 19, 47). Cardiac vagal activity increases markedly even after 2 wk of regular training (30, 50). Previous studies have also shown that postexercise HR recovery improves after endurance training in sedentary healthy males (46) and in patients with coronary artery disease (18). However, a substantial heritable component may also be involved in the regulation of HR behavior in response to training (40, 41), and this may partly explain the large interindividual variation in HR recovery.
Muscarinic acetylcholine receptors play a fundamental role in cardiac function via vagally mediated regulation of the autonomic nervous system. The human heart expresses predominantly muscarinic acetylcholine receptor subtype M2 (CHRM2) (6, 8). Fisher et al. (14) showed that vagally induced bradycardiac responses were totally abolished in CHRM2-deficient mice in vivo, suggesting the exclusive role of CHRM2 in HR regulation (14). Therefore, we tested the hypothesis that, among healthy individuals, the CHRM2 gene polymorphisms are associated with 1-min HR recovery after peak exercise in the sedentary state and after endurance training.
| METHODS |
|---|
|
|
|---|
The subjects were recruited by newspaper ads, which attracted 355 replies. All smokers, subjects with body mass index (BMI) > 32 kg/m2, subjects who did regular physical training more than two times per week, and subjects with diabetes mellitus, asthma, or cardiovascular disorders were excluded. We invited 108 subjects to our laboratory (Dept. of Exercise and Medical Physiology, Merikoski Rehabilitation and Research Centre, Oulu, Finland) for more specific assessment. The number of subjects was selected on the basis of a priori power analysis to give 80% power to the responsiveness of maximal aerobic power [peak oxygen consumption (
O2 peak)] after 2 wk of endurance training. According the previous studies (16, 30, 31, 50), 8% improvement in
O2 peak can be expected after short-term endurance training with SD of approximately ±7%, giving the estimate of needing at least 10 subjects for both the intervention and the control group. The subjects were randomized into a training group (n = 90) and a control group (n = 18). The unbalanced study design was used to particularly assess the contribution of genotypes to actual training effect. The average heterozygosity and allele frequencies of DNA sequence variation of interest are available from the National Center for Biotechnology Information (NCBI) database. These frequencies were used to estimate the hypothetical number of subjects of each genotype group at the present population. During the study, 13 subjects dropped out (10 from the training group and 3 from the control group) because of personal or health-related problems. Finally, 80 subjects (36 men and 44 women) performed the endurance training program. The control group consisted of 15 subjects (9 men and 6 women). The Ethical Committee of the Northern Ostrobothnia Hospital District, Oulu, Finland, approved the protocol.
Experimental Design
Sequence of tests. On their first laboratory day, the subjects completed a health status questionnaire, gave written informed consent, and were assessed for BMI. Resting electrocardiogram (12-lead ECG) was recorded to confirm their cardiac health status, and overnight R-R intervals were recorded to evaluate autonomic regulation by the HR variability method. On the second day, blood pressure was recorded and maximal exercise testing performed. Use of alcohol or strenuous physical activity was not allowed during the test days and on the two preceding days.
The endurance training period was 2 wk, including five consecutive sessions per week (MondayFriday). The short intervention period was chosen on the basis of the previous studies, which shows that a well-controlled, short-term training intervention increased significantly
O2 peak (16, 20, 30, 31, 45, 50) and cardiac vagal activity (30, 50). At the end of the training intervention, all measures were repeated 48 h after the last training session. The control group was tested similarly to the training group, and they were asked to maintain their habitual physical activity level during the 2-wk period. Finally, blood samples were collected 2 days after the last exercise testing in the laboratory of Internal Medicine at the University Hospital of Oulu.
Assessment of
O2 peak, maximal HR, and HR recovery.
The subjects performed a graded maximal exercise test on an 839E Monark cycle ergometer (Stockholm, Sweden). The test was started at 25 W, and work rate was increased by 25 W every 2 min until exhaustion. Pedaling was started at 50 rpm and increased by 5 rpm up to 90 rpm, to achieve maximal effort. The subjects were encouraged to continue cycling until they could no longer maintain the required pace, at which time the test was terminated. After the termination of the test, the subjects were asked to remain seated on the bicycle, and HR recording was continued for 6 min. They were not allowed to talk or move during the first minute after the test, which was followed by a 5-min cool-down period at a work rate of 25 W. HR recovery was determined as maximal HR (mean of 10 s) HR at 1 min after exercise (mean of 5 s). Ventilation, gas exchange (M909 Ergospirometer, Medikro, Kuopio, Finland), and HR (Cardiolife TEC-7721K, Nihon Kohden, Tokyo) were monitored continuously during the protocol. The highest oxygen uptake value measured during the test (1-min collection) was taken as the
O2 peak. All subjects fulfilled the criteria for
O2 peak given in the literature (i.e., respiratory exchange ratio >1.1 or maximal HR within ±10 beats of the age-appropriate reference value) (21). HR was recorded during the exercise test and recovery with a Polar R-R Recorder (Polar Electro, Kempele, Finland) and saved in a computer for further analysis with the HEARTS software (Heart Signal, Kempele, Finland).
Assessment of resting blood pressure. Resting systolic (SBP) and diastolic blood pressures (DBP) were assessed by using an electronic sphygmomanometer (Omron M4, Omron Healthcare). The subjects lay in the supine position in a quiet room for at least 10 min before the blood pressure measurements. The blood pressure measurements were performed in a supine position at the same time of the day before and after the training program.
Assessment of HR variability. The R-R intervals were recorded overnight (from midnight to 6 AM) with a Polar R-R Recorder at an accuracy of 1 ms and saved in a computer for further analysis of HR variability with the Hearts software. All R-R intervals were edited by visual inspection on the basis of ECG portions to exclude all undesirable beats, which accounted for <2% in every subject's recording. The details of this analysis and the filtering technique have been described previously (22). The subjects were asked to go to bed before midnight and to stay in bed until 6 AM on the R-R interval recording days.
The mean HR and the SD of all R-R intervals were used as time-domain measures of HR variability. An autoregressive model was used to estimate the power spectrum densities of R-R interval variability. Low-frequency power (LF, 0.040.15 Hz) and high-frequency power (HF, 0.150.4 Hz) were calculated from the segments of the 512 R-R interval overnight recording (36). A logarithmic transformation to the natural base was performed on both spectral components of HR variability.
Exercise training. The endurance training program consisted of cycling on an 818E Monark cycle ergometer for 40 min (Stockholm, Sweden). Each exercise session consisted of a 5-min warm-up period (cycling at 50-W and 75-W resistance for women and men, respectively), followed by 30 min of cycling at a resistance that elicited a HR of 7080% of maximal HR, and ended with a 5-min cool-down period (cycling at 50-W and 75-W resistance for women and men, respectively). Exercise intensity was closely monitored using a Polar Electro HR monitor A1, and the mean HR of each training session was recorded. Endurance exercise training was planned on the basis of the recommendations of American College of Sports Medicine (2). A professional instructor supervised each training session.
Determination of genotypes. EDTA blood samples were obtained from all subjects, and DNA was isolated by using simple salting-out procedure. Six single-nucleotide polymorphisms (SNPs) were selected from the NCBI dbSNP database: rs2278098 in intron 3 (1780 bp downstream from exon 3), rs2061174 in intron 4 (25,573 from exon 4), rs324640 and rs324650 in intron 5 (10,571 and 5,906 bp, respectively, upstream of exon 6), rs8191992 in the 3'-untranslated region of exon 6, and rs1378650 located 152 bp downstream of exon 6. These SNPs were chosen on the basis of their clinical relevance (48) and their capturing role of common genetic variation at the CHRM2 locus. The SNPs were genotyped by using template-directed dye-terminator incorporation with the fluorescence polarization detection (FP-TDI) method. A DNA sequence containing the SNP was amplified with PCR, and after cleaning with shrimp alkaline phosphatase and exonuclease I, the PCR product was used as a template in the AcycloPrime-FP reaction. The SNP detection primers were designed so as to locate their 3'-end immediately upstream of the polymorphic (SNP) site. The SNP detection PCR reaction utilized a mutant thermostable polymerase and a pair of AcycloTerminators labeled with R110 and TAMRA (AcycloPrime-FP SNP kit, PerkinElmer Life Sciences, Boston, MA), representing the possible alleles for the SNP of interest. Allele detection was done by measuring the changes in fluorescence polarization after excitation of the samples by plane-polarized light by using a Victor2 FP Plate Reader (PerkinElmer Life Sciences) on a 384-well plate format. Allele calling was done by using the SNPscorer genotyping software (PerkinElmer Life Sciences). Haplotypes were constructed using the SHEsis software package (38).
Statistical Analyses
Chi-square test was used to verify whether the observed genotype frequencies were in Hardy-Weinberg equilibrium, and pairwise linkage disequilibrium (LD) between the SNPs was assessed using the ldmax program available in the GOLD software package (1). The normal Gaussian distribution of the data was verified by the Kolmogorov-Smirnov goodness-of-fit test. The difference in change in
O2 peak after training within the training and control group was analyzed by a two-factor ANOVA for repeated measures with time and interventions followed by post hoc analysis (Student's paired t-test). The associations between the CHRM2 gene polymorphisms, the end-point phenotypes, and the associations with the haplotypes were tested with an analysis of covariance by using the general linear model procedure of the SAS software package. The phenotypes were adjusted for age, sex, and BMI. The training response phenotypes were adjusted in addition to the baseline value of the phenotype. Pearson's correlation coefficients were calculated to study the associations between the HR phenotypes, HR variability phenotypes, physical performance phenotypes, and resting blood pressure before and after the training program. Pearson's correlation analysis was also used to test the relationships between HR recovery and the other measured phenotypes separately within the three genotype groups (rs324640 in Intron 5; C/C, C/T, T/T). Values are given as means ± SD. The SAS statistical software package (SAS Institute) was used for the analyses.
| RESULTS |
|---|
|
|
|---|
|
|
All the subjects in the training group maintained the training frequency of 10 times during 2 wk. The intensity of training was 74 ± 3% of maximal HR, and the mean duration of exercise was 40 ± 3 min/session. The change in
O2 peak was +8.1 ± 5.7% (P < 0.001) for the training group and 0.5 ± 7.0% (P = 0.941) for the control group. There were no differences in the change of
O2 peak between the genotypes in the CHRM2 SNPs (e.g., for rs324640 SNP the improvement of
O2 peak was 7.5 ± 7.2% for C/C, 8.2 ± 7.2% for C/T, and 8.0 ± 5.7% for T/T, P = 0.924).
Association Between CHRM2 SNPs and HR Recovery
Maximal HR was not significantly associated with the CHRM2 SNPs. However, HR recovery differed significantly between the rs324640 genotypes at baseline (C/C, 33 ± 10; C/T, 33 ± 7; and T/T, 40 ± 11 beats/min, P = 0.008, Fig. 1A, Table 3). The association was even stronger (P = 0.001) after the 2-wk endurance training period: the C/C homozygotes showed 6 and 12 beats/min lower HR recovery than the C/T heterozygotes and the T/T homozygotes, respectively (Fig. 1B). The change in HR recovery after the training program also differed between the rs324640 genotypes (C/C, 3 ± 7; C/T, 2 ± 7; and T/T, 2 ± 8 beats/min, P = 0.038). Moreover, rs8191992 SNP was associated with HR recovery: the A/A homozygotes and the A/T heterozygotes showed lower HR recovery values than the T/T homozygotes both at baseline (P = 0.025, Table 4) and after the endurance training program (P = 0.005, Table 4). The SNPs rs2278098, rs2061174, and rs1378650 were not associated with HR recovery.
|
|
|
Overnight HR was not significantly associated with CHRM2 SNPs. However, the LF-to-HF ratio showed significantly higher values in the rs8191992 A/A homozygotes than in the other genotypes both before (P = 0.036) and after (P = 0.046) endurance training (Table 4). The rs324640 C/C homozygotes showed a similar trend (P = 0.064 and P = 0.063, respectively), but the difference was not statistically significant (Table 3).
Association Between CHRM2 SNPs and Blood Pressure
Resting DBP was significantly higher in the rs324640 C/C homozygotes than in the other genotypes both before (P = 0.049) and after (P = 0.045) the training program (Table 3). Similar trends were also observed for SBP, but the differences did not reach statistical significance.
Correlation Between SBP and HR Recovery Between rs324640 Genotypes
There were no statistically significant correlations between HR recovery, HR variability, and the other hemodynamic phenotypes before (e.g., HR recovery vs. HF power, r = 0.016; and HR recovery vs. LF-to-HF ratio, r = 0.064) or after the training program (e.g., HR recovery vs. HF power, r = 0.033; and HR recovery vs. LF-to-HF ratio, r = 0.135). However, the relationship between resting SBP and HR recovery varied considerably between the rs324640 genotypes, especially in the trained state (Fig. 2): the common T/T homozygotes (Fig. 2F) showed a strong positive correlation [r = 0.517, 95% confidence interval (CI) from 0.20 to 0.74], whereas the C/C homozygotes (Fig. 2B) showed a negative correlation (r = 0.353, CI from 0.72 to 0.17) and the C/T heterozygotes (Fig. 2D) an intermediate correlation (r = 0.087, CI from 0.22 to 0.38).
|
The CHRM2 haplotypes were constructed by using the "best" option of the SHEsis software haplotyping function. The haplotype analyses confirmed the associations detected with individual SNPs. The haplotype consisting of the minor alleles of the rs324640 and rs8191992 variants (CA) was associated with lower HR recovery values in the trained state (P = 0.007). The minor allele haplotype also showed a stronger association with the LF-to-HF ratio both at baseline (P = 0.004) and after the training period (P = 0.007) compared with the association detected in rs8191992 alone. The common allele (TT) haplotype of the rs324640 and rs8191992 variants was associated with higher HR recovery both before (P = 0.003) and after the training program (P = 0.003). Furthermore, the TT haplotype showed lower values of resting DBP at the baseline (P = 0.037) and after the training (P = 0.011).
| DISCUSSION |
|---|
|
|
|---|
Physiological Background of HR Recovery
CHRM2 plays a fundamental role in cardiac autonomic regulation (14). Activation of cardiac vagal efferents leads to release of acetylcholine, which acts on cardiac CHRM2 to decrease HR (6, 8). Several studies have shown that the main physiological mechanism underlying postexercise cardiodeceleration is vagal reactivation (4, 10, 24, 25, 35). A recent study by Smith and colleagues (43) showed that mongrel dogs vulnerable to ventricular fibrillation had reduced HR recovery compared with dogs resistant to malignant arrhythmias. The differences in HR recovery and cardiac vagal activity were eliminated by administration of atropine, which confirmed the dominant role of vagal activity in the control of HR after exercise (43). Furthermore, acetylcholine receptor resistance in rats resulted in reduced cardiac muscarinic receptor function, leading to cardiovagal insufficiency (34). As a novel finding of the present study, we found that genetic variation at the CHRM2 locus is associated with interindividual variability in the modulation of HR after maximal exercise.
The role of HR recovery as a clean index of cardiac vagal function has been questioned because the sympathetic withdrawal is also involved (37). In the present study, HR recovery did not associate with HF power of R-R intervals, which is the most commonly used index of cardiac vagal outflow. Analysis of HF power from ambulatory R-R interval recordings includes also some methodological problems, e.g., the saturation of HF power during the very high vagal outflow (26). This may be one reason for the low association between HR recovery and ambulatory measured HF power. Second, the CHRM2 gene polymorphisms were slightly but significantly associated with LF-to-HF ratio, known as an index of sympathovagal balance. Taken together, it is possible that both branches of autonomic regulation contribute to the magnitude of HR recovery after maximal exercise.
Association Between CHRM2 SNPs and Central Nervous System
The central nervous system plays an important role in cardiovascular regulation (15). Locally released acetylcholine acts through CHRM2 within the ventrolateral medulla in the brain, increasing the activity of sympathetic preganglionic neurons and consequently elevating blood pressure and HR (15, 27, 28). It was previously observed that the same SNPs of the CHRM2 gene as reported in the present study showed a strong association with alcohol dependence, depressive disorder, and various biopotential waves, some of which are also associated with cardiovascular dysfunctions (48). Our subjects were healthy but were not assessed for alcohol dependence or depressive disorders. However, we were able to confirm the association between the CHRM2 polymorphisms and early cardiovascular dysfunction in various autonomic markers by, for instance, verifying slow HR recovery in untrained and trained states and an elevated LF-to-HF ratio in long-term R-R interval recording.
Effects of Training on HR Recovery
Endurance training has been suggested to protect the heart against harmful cardiac events by increasing vagal tone (5). Cardiac vagal activity has been shown to increase after 2 wk of regular endurance training (30, 50). Previous studies on sedentary healthy men and patients with severe coronary artery disease have reported endurance training induced improvement in cardiac autonomic regulation, quantified as improved postexercise HR recovery (18, 46). In the present study, the minor allele homozygotes of the rs324640 and rs8191992 variants showed less optimal HR recovery than the common allele homozygotes, both in the sedentary state and especially after the 2-wk training program. Furthermore, the determinant of sympathovagal balance, i.e., LF-to-HF ratio, showed an association with the minor allele (CA) haplotype of the rs324640 and rs8191992 variants, confirming the contribution of the CHRM2 gene locus to the autonomic regulation of the heart.
Role of Acetylcholine Receptors in Blood Pressure Regulation
Altered autonomic modulation of HR has been shown to be associated with elevated blood pressure (23) and early stage of hypertension among normotensive men (42). Systemic administration of acetylcholine and other muscarinic agonists lowers blood pressure, presumably because of vasodilatation caused by the activation of muscarinic acetylcholine receptors located on vascular endothelial cells (14). Fisher et al. (14) reported recently that the blood pressure response in CHRM3-deficient mice was greatly attenuated after the administration of muscarinic agonist compared with wild-type control mice. Furthermore, they observed that the CHRM2-mediated decrease in HR was the major factor responsible for the pronounced reduction in blood pressure after the administration of muscarinic agonist, suggesting a dominant role of CHRM receptors in the control of blood pressure level (14). In our healthy normotensive series, we observed that rs324640 C/C homozygotes also tended to have higher resting blood pressure values before and after endurance training. Furthermore, the relationship between resting blood pressure and postexercise HR recovery varied markedly across the rs324640 genotypes, especially after endurance training. Taken together, these results suggest that the CHRM2 gene locus may contribute to interindividual variation in blood pressure levels. However, this hypothesis should be confirmed in larger populations.
Limitations
The present study was limited by the fact that a relatively small number of subjects were examined, which may emphasize the preliminary nature of the results. Therefore, the results must obviously be confirmed in other populations. Second, to date, there is no functional evidence of the associated CHRM2 SNPs, although the role of CHRM2 in the regulation of human HR is evident at the receptor level.
Implications
The importance of HR recovery after both maximal (10, 25) and submaximal exercise (11) as a predictor of mortality has been well established in large population-based cohorts. The independent predictive value is evident both in healthy subjects (25) and in different patient groups (10, 32, 33, 44). Indeed, the American Heart Association has proposed that HR recovery could be used to assess the risk for developing clinical coronary disease among asymptomatic subjects (29). In this respect, the findings of the present study may have clinical relevance and contribute to our understanding of the genetic background of HR recovery after exercise. In summary, DNA sequence variation in the CHRM2 gene locus is a potential modifier of postexercise HR recovery in healthy individuals. Whether the same findings apply to other ethnic groups and to cardiac patients remains to be confirmed in future studies.
| GRANTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
P. Palatini Exercise testing in asymptomatic subjects: from diagnostic test to prognostic tool? Eur. Heart J., August 1, 2008; 29(15): 1803 - 1806. [Full Text] [PDF] |
||||
![]() |
D. M. Dick, F. Aliev, J. C. Wang, R. A. Grucza, M. Schuckit, S. Kuperman, J. Kramer, A. Hinrichs, S. Bertelsen, J. P. Budde, et al. Using Dimensional Models of Externalizing Psychopathology to Aid in Gene Identification Arch Gen Psychiatry, March 1, 2008; 65(3): 310 - 318. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. A. Dewland, A. S. Androne, F. A. Lee, R. J. Lampert, and S. D. Katz Effect of acetylcholinesterase inhibition with pyridostigmine on cardiac parasympathetic function in sedentary adults and trained athletes Am J Physiol Heart Circ Physiol, July 1, 2007; 293(1): H86 - H92. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Ingelsson, M. G. Larson, R. S. Vasan, C. J. O'Donnell, X. Yin, J. N. Hirschhorn, C. Newton-Cheh, J. A. Drake, S. L. Musone, N. L. Heard-Costa, et al. Heritability, Linkage, and Genetic Associations of Exercise Treadmill Test Responses Circulation, June 12, 2007; 115(23): 2917 - 2924. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |