- Research article
- Open Access
Cardiac performance, biomarkers and gene expression studies in previously sedentary men participating in half-marathon training
BMC Sports Science, Medicine and Rehabilitation volume 6, Article number: 6 (2014)
The mechanisms through which exercise reduces cardiovascular disease are not fully understood. We used echocardiograms, cardiac biomarkers and gene expression to investigate cardiovascular effects associated with exercise training.
Nineteen sedentary men (22–37 years) completed a 17-week half-marathon training program. Serial measurements of resting heart rate, blood pressure, maximum oxygen consumption, lipids, C-reactive protein, cardiac troponin T, echocardiograms and blood for gene expression were obtained from baseline to peak training. Controls included 22 sedentary men who did not exercise.
Among the training group, VO2 max increased from 37.1 to 42.0 ml/kg/min (p < 0.001). Significant changes were seen in left ventricular wall thickness and mass, stroke volume, resting heart rate and blood pressure (p < 0.001). The control group demonstrated no significant changes. Expression profiling in the training group identified 10 significantly over-expressed and 53 significantly under-expressed loci involved in inflammatory pathways. Dividing the training group into high and low responders based on percent change in VO2 max identified loci that differentiated these two groups at baseline and after training.
Intensive exercise training leads to significant increase in cardiac and hemodynamic performance, and significant changes in expression of genes involved in immune and inflammatory response.
Long-term exercise training has been associated with weight loss, improvements in lipid profile and blood pressure, and cardiac remodeling [1–3]. Despite this, the cardiovascular response to high-intensity training in previously untrained individuals is not well characterized. Importantly, even with extreme marathon and half-marathon (HM) exercise training, individuals may experience varying levels of cardiovascular and physiological changes, and these differences may have a biological basis. For example, HM and marathon training have been associated with increases in oxidative stress and myocardial damage [4–9]. HM and marathon training also induces the release of troponin I and T, markers of myocardial injury; elevations in cardiac troponin T (cTnT) are reported in up to 39-56% of athletes immediately after performance . Although the mechanism and implications are not well defined, these biomarkers usually return to normal within 24 to 48 hours, and the magnitude of biomarker peak is inversely related to the level of training prior to the event . The purpose of our study was to assess time-related cardiovascular performance and function and corresponding gene expression changes in participants undergoing HM training. Our goal is to begin to identify a genomic basis for response to exercise by correlating specific changes in gene expression with the physiologic changes following HM training.
Apparently healthy, untrained and sedentary males (n = 47), between 18 and 45 years of age, participated in the study. Candidates were required to have engaged in less than 2 hours of aerobic training per week prior to enrollment, have an average maximum oxygen consumption (VO2 max) less than 50 ml/kg/min, and not previously completed a HM. Twenty-four of the subjects, who were enrolled in a local HM Training Program, were recruited into the training group (TRAIN), and twenty-three male subjects of similar age were recruited into the control group (CON). The CON group did not participate in the training program or undergo any lifestyle modification over the 16 week duration of the study. All subjects volunteered for the study and signed informed consents. TRAIN subjects were compensated for meals, travel time, and training program fees. Exclusion criteria included smoking, known metabolic or cardiovascular conditions (diabetes or dyslipidemia), chronic inflammatory diseases, or chronic use of non-steroidal anti-inflammatory drugs. The protocol was approved by the University of Miami Institutional Review Board for human subjects research.
Initial screening and evaluation included a complete medical and training history, baseline VO2 max test, resting heart rate, blood pressure, height, weight, body composition assessment, serum biomarkers, and complete echocardiographic evaluation. Body composition was assessed through dual-energy X-ray absorptiometry scanning. Blood for the lipid panels and serum biomarkers was drawn via the antecubital vein using the LeukoLock Protection system and stored at −80°C for analyses with the post-training RNA samples. Lipid panels were done on fasted blood samples. TRAIN and CON subjects had heart rate, blood pressure, height, weight, and serum biomarkers repeated at 6, 11, and 16 weeks after the baseline evaluation. Echocardiograms, VO2 max, and body fat evaluation were recorded at baseline and 16 weeks after the baseline evaluation.
Of the 24 subjects enrolled in the TRAIN group, 5 subjects were unable to complete the training due to injury or time constraints and withdrew prior to the peak training data (week 16) collection. Three subjects were able to complete the peak training data collection, but were unable to complete the HM due to injury. Thus, a total of 19 subjects completed the program with peak training data collected. Of these 19 subjects, 16 subjects completed the 2007 Miami HM and had blood drawn within 30 minutes after the race (see Figure 1). Of the CON group, one individual was unable to continue the study as a result of scheduling difficulty, leaving a total of 22 controls that completed all data collection.
Two-dimensional echocardiograms with pulsed-Doppler imaging were acquired on each subject using a commercial system (Philips HD11). Standard measurements were performed according to the American Society of Echocardiography guidelines . Interventricular septal wall thickness, LV internal diameter, and LV posterior wall thickness were measured in systole and diastole in the parasternal long axis view. The LV end-diastolic volume, LV end-systolic volume, stroke volume, and ejection fraction were calculated via modified bi-planer calculation based on the apical 2 and 4 chamber view. LV mass was calculated using the area length method. With the pulsed-Doppler trans-mitral filling velocities (peak early [E’] and late [atrial-A’]), diastolic function was assessed and an E/A ratio was determined. One certified echocardiogram technician was responsible for acquiring all images for each participant. The measurements and calculations were independently confirmed by an experienced echocardiographer. Both were blinded to each subject’s group status.
Blood was collected in EDTA and standard red top tubes (serum) and processed immediately. Tubes were immediately centrifuged, and the serum was separated and stored in a -80° freezer. Serum was analyzed for cTnT , CRP , and lipids (HDL, LDL, and TG) . Blood for RNA extraction was collected at both baseline and peak training (week 16) using the Ambion LeukoLOCK system that captured and preserved the blood mononuclear cell layer.
Clinical and biomarker data were analyzed using SPSS for Windows 15.0 software (SPSS, Inc, Chicago, IL). Frequency and descriptive statistics were calculated to check all relevant characteristics of the data. Independent samples t-tests were used to compare baseline variable differences between the TRAIN and CON groups. Paired samples t-tests were used to compare changes from baseline to peak training or baseline to post-HM.
Microarray analysis was performed on blood collected from each subject at baseline and at 16 weeks. Total RNA was extracted and analyzed using the Affymetrix GeneChip Human Gene ST array, which has 28,869 unique genes (or transcripts) which are captured by ~26 probes per gene. Data were normalized using Robust Multi-Array (RMA) Normalization. Significance Analysis of Microarray (SAM) Software  was used to determine statistically significant genes that met predetermined False Discovery Rate (FDR) thresholds. SAM uses non-parametric statistics and gene specific t-tests to identify statistically significant genes and to measure the strength of relationships between gene expression and a specific response variable. SAM also uses permutation tests to estimate the FDR after accounting for multiple testing corrections . Genes with FDR p-values < 0.15 and fold-change greater than 1.5 or less than −1.5 were considered significant. Pathway analysis was conducted using MetaCore™ (GeneGo, Inc) software.
Two separate gene expression analyses were performed. The first compared differences in gene expression profiles between the CON and TRAIN subjects before and after the HM training. The second analysis focused only on the TRAIN subjects, by dividing them into high and low responders based on the percent change of their relative VO2 max (rVO2 max) over the course of the training program. Relative VO2 max takes into account and individual’s body weight and is defined as milliliters of oxygen per kilogram of bodyweight per minute. The mean percent change in rVO2 max was calculated and individuals who fell in the upper quartile (upper 25%) of percent change in rVO2 max were labeled as high-responders (10 individuals) and those who fell below the upper quartile were considered low-responders (9 individuals). Gene expression differences between the two groups were studied at baseline and after the 16 weeks of training.
Of the 16 subjects who completed the HM, the average time to completion was 2:22:21 (H:Min:Sec), SD = 26:33. The fastest time was 1:39:00 and the slowest time was 3:11:58. No runner required medical attention after the race. There were no significant differences between the TRAIN and CON groups on any baseline characteristics other than triglycerides, which were significantly higher (p = 0.03) in the TRAIN group (Table 1).
Comparisons of basic physiological variables were made from baseline to the peak training time-point (17 weeks) in both groups (Table 1). VO2 max increased significantly (p = 0.00004) among those in the TRAIN group during the course of the training program. In addition, the following were significantly decreased: (1) resting heart rate (p = 0.002), (2) percent body fat (p = 0.044), and (3) systolic (p = 0.0001) and diastolic (p = 0.007) blood pressure. Average body weight (p = 0.08) and HDL (p = 0.055) decreased slightly (p = 0.08), but no other measured biomarkers (LDL, Tg, CRP, cTnT) showed significant improvements. There were no significant changes for any variables among the CON subjects (Table 1). Table 2 shows comparison of echocardiographic changes between TRAIN and CON groups.
In comparing the gene expression profiles between TRAIN and CON individuals before training, only 4 genes (CMKLR1, FOS, IFIT1, and PDK4) showed significant changes (FC > 1.5, FDR <0.05) (See Table 3). All had increased levels of gene expression in the TRAIN group compared to the CON group. No genes were significantly decreased.
When comparing profiles from blood collected in TRAIN and CON individuals after the HM training program, ten genes (CMKLR1, FCRL6, GPR56, HSPA6, IFNG, IL1B, KLRC3, MIR15A) showed significantly elevated expression (FC > 1.5, FDR < .05) (Table 3) and 53 genes (including AREG,DEFA3) were found to be significantly decreased (FC < −1.5, FDR < .05) in the TRAIN group (Table 4). According to pathway analysis, the biological pathways that were significantly up-regulated in the TRAIN individuals after training (p <0.001) were those involved in inflammatory and immune response, including migration inhibitory factor (MIF) in innate immunity response (p = .00014), bacterial infections in normal airways (p = .0002), apoptosis and survival Nitric Oxide synthesis and signaling (p = .0002) and bacterial infections in CF airways (p = .0002) (data not shown). One biological pathway, translation insulin regulation of translation (p = .00008), showed significantly decreased expression in TRAIN individuals after HM training (data not shown).
The TRAIN group was categorized into 10 high and 9 low responders based on percent change in their rVO2 max. Within the TRAIN group alone, 16 genes showed trending evidence of decreased expression (FC > −1.5, FDR < .06) in high vs. low responders, and one gene showed increased gene expression (FC > 1.5, FDR < .06) in the high compared to low responders group (Table 5). Biological processes that were highlighted included those involved in regulation of DNA transcription, translation and other metabolic processes.
This is the first prospective study evaluating the combination of changes in cardiac structure and function, peripheral blood cell gene expression, inflammatory markers, and potential for myocardial injury associated with training for, and performance of, a HM in previously sedentary subjects. The decreases in blood pressure, resting heart rate, and body fat percentage, and the increases in VO2 max, observed in our TRAIN group over the course of the program are consistent with previously well-documented responses to exercise training [15, 16] and validate the training effect in our HM participants. More interesting is the effect of HM training on biomarkers related to myocardial injury, cardiovascular structural changes by echocardiogram, serum lipids, gene expression and serum inflammatory markers.
In our study, LV mass and wall thickness, and LV end-diastolic and stroke volumes increased significantly after 16 weeks of training. The E to A ratio also decreased throughout training, consistent with a decrease in diastolic relaxation commonly seen in hypertrophied ventricles. These echocardiographic changes have been well-established and recently described to occur in endurance athletes in as little as 90 days [17–19]. This phenomenon, known as the “athlete’s heart,” is thought to be due to a long-sustained high cardiac output resulting in a volume-loaded heart that can lead to an increase in LV cavity, LV mass, and wall thickness [20–22].
Both TRAIN and CON groups had the same inclusion criteria, which enabled us to maximize the comparison after training. When comparing gene expression profiles in blood collected before training in TRAIN versus CON men, 4 genes showed significant changes in gene expression (FC = (1.5-1.7) FDR < .05) suggesting that the two groups were fairly similar before training. Through pathway analysis, the biological pathways that were significantly up regulated in the TRAIN individuals after training were those primarily involved in the body’s inflammatory and immune response. In looking at the specific genes involved of the 53 genes that showed significantly decreased expression in TRAIN individuals after training, Amphiregulin (AREG) was most significantly decreased (FC = −2.42). The protein coded by AREG is a member of the epidermal growth factor family. AREG has both growth factor and inhibitory properties (NCBI, PubMed Gene ID:374). Changes in AREG expression in human neutrophils have been related to exercise in previous studies. For example, following 30 minutes of aerobic exercise, AREG was shown to increase by 3.4 fold . Although the change is in the opposite direction of the change found in the current study this could be due to the fact that Radom-Aizik et al.  examined gene expression changes following acute exercise, whereas our focus is on gene expression changes due to long term, consistent exercise.
Another gene significantly down-regulated in the TRAIN individuals; defensin, alpha 3, neutrophil-specific (DEFA3) also has a role in both the immune and inflammatory response (NCBI, PubMed Gene ID: 1668). A previous study linked increased expression of DEFA3 at both the mRNA and protein level to inflammatory diseases . It is well established that regular exercise helps reduce the amount of total body inflammation and therefore helps prevent the induction of many chronic inflammatory diseases including atherosclerosis. A decrease in total body inflammation is usually associated with a decrease in C-reactive protein, which was not observed in our study. This data suggests that perhaps gene expression changes associated with a decrease in inflammation occur prior to changes in cardiac biomarkers or that subjects with normal CRP at baseline are less likely to show a measurable response to exercise associated with changes in gene expression. Perhaps with longer training, or elevated baseline levels of CPR, similar changes in the cardiac biomarkers would also be seen.
When comparing the TRAIN high responders and low responders, several gene expression differences were noted. Two genes that showed increased expression in high responders before training included FGFBP2 (FC = 1.76) and CD160 (FC = 1.55). Both are involved in immune response (NCBI, PubMED Gene ID: 11126).
Of the genes that showed decreased expression in high responders after exercise of particular note is LRRFIP1 (FC = −1.57). LRRFIP1 has been identified as having a role as both a regulator of toll-like receptor pathway signaling and transcriptional repression. Its expression has been previously shown to be induced by nicotine and proinflammatory cytokines .
In summary, in previously sedentary individuals, intensive training for a HM over 17 weeks leads to increased cardiac and hemodynamic performance that is not associated with longitudinal changes in known cardiac biomarkers. Elevated cTnT occurred in only 1 of 16 subjects who completed the HM. Training is also associated with changes in blood gene expression that reflects the anabolic state of the body, as well as changes in the body’s immune and inflammatory response that could be more sensitive precursors to changes in serum biomarkers. Whereas some genes involved in the immune/inflammatory response showed increased expression, others showed decreased expression implying that exercise brings about changes that both suppress aspects of the immune system while stimulating others. Genes with changes in gene expression as response to HM-training include AREG, DEFA3, CD160, GPR56, KLRC3, and LRRFIP1. The results presented here are from a single sample and do not withstand conservative corrections for multiple comparisons. The study’s sample is limited to men and therefore cannot necessarily be generalized to the broader population which could include females However, they generate interesting hypotheses about a possible role of genes involved in immunity and inflammation in responding to intensive exercise training, which bear examination in other data sets.
Bernadet P: Benefits of physical activity in the prevention of cardiovascular diseases. J Cardiovasc Pharmacol. 1995, 25 (Suppl 1): S3-S8.
Giannuzzi P, Mezzani A, 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, Working Group on Cardiac Rehabilitation and Exercise Physiology. European Society of Cardiology: Physical activity for primary and secondary prevention. Position paper of the Working Group on Cardiac Rehabilitation and Exercise Physiology of the European Society of Cardiology. Eur J Cardiovasc Prev Rehabil. 2003, 10: 319-327. 10.1097/01.hjr.0000086303.28200.50.
Physical activity and cardiovascular health. NIH consensus development panel on physical activity and cardiovascular health. JAMA. 1996, 276: 241-246. http://jama.jamanetwork.com/article.aspx?articleid=405583.
Siegel AJ, Lewandrowski EL, Chun KY, Sholar MB, Fischman AJ, Lewandrowski KB: Changes in cardiac markers including B-natriuretic peptide in runners after the Boston marathon. Am J Cardiol. 2001, 88: 920-923. 10.1016/S0002-9149(01)01910-5.
Briviba K, Watzl B, Nickel K, Kulling S, Bos K, Haertel S, Rechkemmer G, Bub A: A half-marathon and a marathon run induce oxidative DNA damage, reduce antioxidant capacity to protect DNA against damage and modify immune function in hobby runners. Redox Rep. 2005, 10: 325-331. 10.1179/135100005X83716.
Neilan TG, Januzzi JL, Lee-Lewandrowski E, Ton-Nu TT, Yoerger DM, Jassal DS, Lewandrowski KB, Siegel AJ, Marshall JE, Douglas PS, Lawlor D, Picard MH, Wood MJ: Myocardial injury and ventricular dysfunction related to training levels among nonelite participants in the Boston marathon. Circulation. 2006, 114: 2325-2333. 10.1161/CIRCULATIONAHA.106.647461.
Neilan TG, Yoerger DM, Douglas PS, Marshall JE, Halpern EF, Lawlor D, Picard MH, Wood MJ: Persistent and reversible cardiac dysfunction among amateur marathon runners. Eur Heart J. 2006, 27: 1079-1084.
Fortescue EB, Shin AY, Greenes DS, Mannix RC, Agarwal S, Feldman BJ, Shah MI, Rifai N, Landzberg MJ, Newburger JW, Almond CS: Cardiac troponin increases among runners in the Boston Marathon. Ann Emerg Med. 2007, 49: 137-143. 10.1016/j.annemergmed.2006.09.024. 143.e1
Bouchard C, Sarzynski MA, Rice TK, Kraus WE, Church TS, Sung YJ, Rao DC, Rankinen T: Genomic predictors of the maximal O uptake response to standardized exercise training programs. J Appl Physiol. 2011, 110: 1160-1170. 10.1152/japplphysiol.00973.2010.
Shave R, George KP, Atkinson G, Hart E, Middleton N, Whyte G, Gaze D, Collinson PO: Exercise-induced cardiac troponin T release: a meta-analysis. Med Sci Sports Exerc. 2007, 39: 2099-2106. 10.1249/mss.0b013e318153ff78.
Schiller NB, Shah PM, Crawford M, DeMaria A, Devereux R, Feigenbaum H, Gutgesell H, Reichek N, Sahn D, Schnittger I: Recommendations for quantitation of the left ventricle by two-dimensional echocardiography. American society of echocardiography committee on standards, subcommittee on quantitation of two-dimensional echocardiograms. J Am Soc Echocardiogr. 1989, 2: 358-367.
Alpert JS, Thygesen K, Antman E, Bassand JP: Myocardial infarction redefined–a consensus document of the joint European society of cardiology/American college of cardiology committee for the redefinition of myocardial infarction. J Am Coll Cardiol. 2000, 36: 959-969. 10.1016/S0735-1097(00)00804-4.
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98: 5116-5121. 10.1073/pnas.091062498.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. Ann Stat. 1995, 57: 289-300.
Kelley GA, Kelley KS: Efficacy of aerobic exercise on coronary heart disease risk factors. Prev Cardiol. 2008, 11: 71-75. 10.1111/j.1751-7141.2008.08037.x.
Kaplan NM: Non-drug treatment of hypertension. Ann Intern Med. 1985, 102: 359-373. 10.7326/0003-4819-102-3-359.
Morganroth J, Maron BJ, Henry WL, Epstein SE: Comparative left ventricular dimensions in trained athletes. Ann Intern Med. 1975, 82: 521-524. 10.7326/0003-4819-82-4-521.
Baggish AL, Wang F, Weiner RB, Elinoff JM, Tournoux F, Boland A, Picard MH, Hutter AM, Wood MJ: Training-specific changes in cardiac structure and function: a prospective and longitudinal assessment of competitive athletes. J Appl Physiol. 2008, 104: 1121-1128. 10.1152/japplphysiol.01170.2007.
Maron BJ: Structural features of the athlete heart as defined by echocardiography. J Am Coll Cardiol. 1986, 7: 190-203. 10.1016/S0735-1097(86)80282-0.
Khamis RY, Mayet J: Echocardiographic assessment of left ventricular hypertrophy in elite athletes. Heart. 2008, 94: 1254-1255.
Naylor LH, George K, O’Driscoll G, Green DJ: The athlete’s heart: a contemporary appraisal of the ‘Morganroth hypothesis’. Sports Med. 2008, 38: 69-90. 10.2165/00007256-200838010-00006.
Mihl C, Dassen WR, Kuipers H: Cardiac remodelling: concentric versus eccentric hypertrophy in strength and endurance athletes. Neth Heart J. 2008, 16: 129-133. 10.1007/BF03086131.
Radom-Aizik S, Kaminski N, Hayek S, Halkin H, Cooper DM, Ben-Dov I: Effects of exercise training on quadriceps muscle gene expression in chronic obstructive pulmonary disease. J Appl Physiol. 2007, 102: 1976-1984. 10.1152/japplphysiol.00577.2006.
Froy O, Sthoeger ZM: Defensins in systemic lupus erythematosus. Ann N Y Acad Sci. 2009, 1173: 365-369. 10.1111/j.1749-6632.2009.04622.x.
Arakawa R, Bagashev A, Song L, Maurer K, Sullivan KE: Characterization of LRRFIP1. Biochem Cell Biol. 2010, 88: 899-906. 10.1139/o10-014.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/2052-1847/6/6/prepub
The authors wish to thank Dr. Martin Bilsker, Director of the Echocardiography lab in Jackson Memorial Hospital, Miami, FL, for reading the echocardiograms in this study. Also special thanks to Diego Duque and Philips for the use of their Philips HD 11 Digital Ultrasound Machine. Thanks to the laboratory and analytical staff of the Hussman Institute for Human Genomics for their support in particular Dr. Demet Tekin. Dr. Myerburg is supported in part by a grant from the Fondation Leducq, Paris, France (Trans-Atlantic Network of Excellence “Alliance Against Sudden Cardiac Death”). Research conducted in this paper was funded in part by funds from Florida’s Office of Tourism, Trade, and Economic Development.
The authors declare that they have no competing interests.
DDV: study design, RNA extraction, gene expression and pathway analysis, manuscript preparation. GC: echocardiogram data collection, cardiac biomarker analysis, manuscript preparation. MS: ascertainment, data collection. RJM: study design, cardiac biomarker analysis, manuscript editing. KJ: ascertainment, data collection. LN: microarray and gene expression analysis. AP: data collection, study design, editing. DS: study design, gene expression analysis, manuscript editing. PJGC: study design, manuscript editing. ER statistical analysis, gene expression and pathway analysis, manuscript editing. All authors read and approved the final manuscript.
Danica D Vance, Gordon L Chen contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Vance, D.D., Chen, G.L., Stoutenberg, M. et al. Cardiac performance, biomarkers and gene expression studies in previously sedentary men participating in half-marathon training. BMC Sports Sci Med Rehabil 6, 6 (2014). https://doi.org/10.1186/2052-1847-6-6
- Gene expression