Impact of exercise training on endothelial transcriptional profiles in healthy swine: a genome-wide microarray analysis

Interaction of genetic background and exercise training intensity on endothelial function in mouse aorta

The purpose of this study was to characterize the genetic contribution to endothelial adaptation to exercise training. Vasoreactivity was assessed in aortas from four inbred mouse strains (129S1, B6, NON, and SJL) after 4 weeks of moderate intensity continuous exercise training (MOD), high intensity interval training (HIT) or in sedentary controls (SED). Intrinsic variations in endothelium-dependent vasorelaxation (EDR) to acetylcholine (ACh) as well as vasocontractile responses were observed across SED groups. For responses to exercise training, there was a significant interaction between mouse strain and training intensity on EDR. Exercise training had no effect on EDR in aortas from 129S1 and B6 mice. In NON, EDR was improved in aortas from MOD and HIT compared with respective SED, accompanied by diminished responses to PE in those groups. Interestingly, EDR was impaired in aorta from SJL HIT compared with SED. The transcriptional activation of endothelial genes was also influenced by the interaction between mouse strain and training intensity. The number of genes altered by HIT was greater than MOD, and there was little overlap between genes altered by HIT and MOD. HIT was associated with gene pathways for inflammatory responses. NON MOD genes showed enrichment for vessel growth pathways. These findings indicate that exercise training has non-uniform effects on endothelial function and transcriptional activation of endothelial genes depending on the interaction between genetic background and training intensity.

Vascular Adaptation to Exercise in Humans: Role of Hemodynamic Stimuli

On the 400th anniversary of Harvey's Lumleian lectures, this review focuses on "hemodynamic" forces associated with the movement of blood through arteries in humans and the functional and structural adaptations that result from repeated episodic exposure to such stimuli. The late 20th century discovery that endothelial cells modify arterial tone via paracrine transduction provoked studies exploring the direct mechanical effects of blood flow and pressure on vascular function and adaptation in vivo. In this review, we address the impact of distinct hemodynamic signals that occur in response to exercise, the interrelationships between these signals, the nature of the adaptive responses that manifest under different physiological conditions, and the implications for human health. Exercise modifies blood flow, luminal shear stress, arterial pressure, and tangential wall stress, all of which can transduce changes in arterial function, diameter, and wall thickness. There are important clinical implications of the adaptation that occurs as a consequence of repeated hemodynamic stimulation associated with exercise training in humans, including impacts on atherosclerotic risk in conduit arteries, the control of blood pressure in resistance vessels, oxygen delivery and diffusion, and microvascular health. Exercise training studies have demonstrated that direct hemodynamic impacts on the health of the artery wall contribute to the well-established decrease in cardiovascular risk attributed to physical activity.

Transcriptomic effects of metformin in skeletal muscle arteries of obese insulin-resistant rats

We examined the effects of metformin, a commonly used antidiabetic drug, on gene expression in multiple arteries. Specifically, transcriptional profiles of feed arteries and second branch order arterioles in the soleus, gastrocnemius, and diaphragm muscles as well as aortic endothelial scrapes were examined from obese insulin-resistant Otsuka Long-Evans Tokushima Fatty rats treated with ( n = 9) or without ( n = 10) metformin from 20 to 32 weeks of age. Metformin-treated rats exhibited a reduction in body weight, adiposity, and HbA1c ( P < 0.05). The greatest number of differentially expressed genes (FDR < 15%) between those treated with and without metformin was found in the red gastrocnemius 2a arterioles (93 genes), followed by the diaphragm 2a arterioles (62 genes), and soleus 2a arterioles (15 genes). We also found that two genes were differentially expressed in aortic endothelial cells (LETMD1 and HMGCS2, both downregulated), one gene in the gastrocnemius feed artery (BLNK, downregulated), and no genes in the soleus and diaphragm feed arteries and white gastrocnemius 2a arterioles. No single gene was altered by metformin across all vessels examined. This study provides evidence that metformin treatment produces distinct gene expression effects throughout the arterial tree in a rat model of obesity and insulin resistance. Genes whose expression was modulated with metformin do not appear to have a clear connection with its known mechanisms of action. These findings support the notion that vascular gene regulation in response to oral pharmacological therapy, such as metformin, is vessel specific. Impact statement This study provides evidence that metformin treatment produces artery-specific gene expression effects. The genes whose expression was modulated with metformin do not appear to have a clear connection with its known mechanisms of action.

Regular Exercise Reduces Endothelial Cortical Stiffness in Western Diet-Fed Female Mice

We recently showed that Western diet-induced obesity and insulin resistance promotes endothelial cortical stiffness in young female mice. Herein, we tested the hypothesis that regular aerobic exercise would attenuate the development of endothelial and whole artery stiffness in female Western diet-fed mice. Four-week-old C57BL/6 mice were randomized into sedentary (ie, caged confined, n=6) or regular exercise (ie, access to running wheels, n=7) conditions for 16 weeks. Exercise training improved glucose tolerance in the absence of changes in body weight and body composition. Compared with sedentary mice, exercise-trained mice exhibited reduced endothelial cortical stiffness in aortic explants (sedentary 11.9±1.7 kPa versus exercise 5.5±1.0 kPa; P<0.05), as assessed by atomic force microscopy. This effect of exercise was not accompanied by changes in aortic pulse wave velocity (P>0.05), an in vivo measure of aortic stiffness. In comparison, exercise reduced femoral artery stiffness in isolated pressurized arteries and led to an increase in femoral internal artery diameter and wall cross-sectional area (P<0.05), indicative of outward hypertrophic remodeling. These effects of exercise were associated with an increase in femoral artery elastin content and increased number of fenestrae in the internal elastic lamina (P<0.05). Collectively, these data demonstrate for the first time that the aortic endothelium is highly plastic and, thus, amenable to reductions in stiffness with regular aerobic exercise in the absence of changes in in vivo whole aortic stiffness. Comparatively, the same level of exercise caused destiffening effects in peripheral muscular arteries, such as the femoral artery, that perfuse the working limbs.

Exercise training causes differential changes in gene expression in diaphragm arteries and 2A arterioles of obese rats

We employed next-generation, transcriptome-wide RNA sequencing (RNA-Seq) technology to assess the effects of two different exercise training protocols on transcriptional profiles in diaphragm second-order arterioles (D2a) and in the diaphragm feed artery (DFA) from Otsuka Long Evans Tokushima Fatty (OLETF) rats. Arterioles were isolated from the diaphragm of OLETF rats that underwent an endurance exercise training program (EX; n = 13), interval sprint training program (SPRINT; n = 14), or remained sedentary (Sed; n = 12). Our hypothesis was that exercise training would have similar effects on gene expression in the diaphragm and soleus muscle arterioles because diaphragm blood flow increases during exercise to a similar extent as in soleus. Results reveal that several canonical pathways that were significantly altered by exercise in limb skeletal muscles were not among the pathways significantly changed in the diaphragm arterioles including actin cytoskeleton signaling, role of NFAT in regulation of immune response, protein kinase A signaling, and protein ubiquitination pathway. EX training altered the expression of a smaller number of genes than did SPRINT in the DFA but induced a larger number of genes with altered expression in the D2a than did SPRINT. In fact, FDR differential expression analysis (FDR, 10%) indicated that only two genes exhibited altered expression in D2a of SPRINT rats. Very few of the genes that exhibited altered expression in the DFA or D2a were also altered in limb muscle arterioles. Finally, results indicate that the 2a arterioles of soleus muscle (S2a) from endurance-trained animals and the DFA of SPRINT animals exhibited the largest number of genes with altered expression.

Exercise-induced differential changes in gene expression among arterioles of skeletal muscles of obese rats

Using next-generation, transcriptome-wide RNA sequencing (RNA-Seq) technology we assessed the effects of exercise training on transcriptional profiles in skeletal muscle arterioles isolated from the soleus and gastrocnemius muscles of Otsuka Long Evans Tokushima Fatty (OLETF) rats that underwent an endurance exercise training program (EX; n = 13), interval sprint training program (SPRINT; n = 14), or remained sedentary (Sed; n = 12). We hypothesized that the greatest effects of exercise would be in the gastrocnemius arterioles. Results show that EX caused the largest number of changes in gene expression in the soleus and white gastrocnemius 2a arterioles with little to no changes in the feed arteries. In contrast, SPRINT caused substantial changes in gene expression in the feed arteries. IPA canonical pathway analysis revealed 18 pathways with significant changes in gene expression when analyzed across vessels and revealed that EX induces increased expression of the following genes in all arterioles examined: Shc1, desert hedgehog protein (Dhh), adenylate cyclase 4 (Adcy4), G protein binding protein, alpha (Gnat1), and Bcl2l1 and decreased expression of ubiquitin D (Ubd) and cAMP response element modulator (Crem). EX increased expression of endothelin converting enzyme (Ece1), Hsp90b, Fkbp5, and Cdcl4b in four of five arterioles. SPRINT had effects on expression of Crem, Dhh, Bcl2l1, and Ubd that were similar to EX. SPRINT also increased expression of Nfkbia, Hspa5, Tubb 2a and Tubb 2b, and Fkbp5 in all five arterioles and increased expression of Gnat1 in all but the soleus second-order arterioles. Many contractile and/or structural protein genes were increased by SPRINT in the gastrocnemius feed artery, but the same genes exhibited decreased expression in red gastrocnemius arterioles. We conclude that training-induced changes in arteriolar gene expression patterns differ by muscle fiber type composition and along the arteriolar tree.

Transcriptomic-Wide Discovery of Direct and Indirect HuR RNA Targets in Activated CD4+ T Cells

Due to poor correlation between steady state mRNA levels and protein product, purely transcriptomic profiling methods may miss genes posttranscriptionally regulated by RNA binding proteins (RBPs) and microRNAs (miRNAs). RNA immunoprecipitation (RIP) methods developed to identify in vivo targets of RBPs have greatly elucidated those mRNAs which may be regulated via transcript stability and translation. The RBP HuR (ELAVL1) and family members are major stabilizers of mRNA. Many labs have identified HuR mRNA targets; however, many of these analyses have been performed in cell lines and oftentimes are not independent biological replicates. Little is known about how HuR target mRNAs behave in conditional knock-out models. In the present work, we performed HuR RIP-Seq and RNA-Seq to investigate HuR direct and indirect targets using a novel conditional knock-out model of HuR genetic ablation during CD4+ T activation and Th2 differentiation. Using independent biological replicates, we generated a high coverage RIP-Seq data set (>160 million reads) that was analyzed using bioinformatics methods specifically designed to find direct mRNA targets in RIP-Seq data. Simultaneously, another set of independent biological replicates were sequenced by RNA-Seq (>425 million reads) to identify indirect HuR targets. These direct and indirect targets were combined to determine canonical pathways in CD4+ T cell activation and differentiation for which HuR plays an important role. We show that HuR may regulate genes in multiple canonical pathways involved in T cell activation especially the CD28 family signaling pathway. These data provide insights into potential HuR-regulated genes during T cell activation and immune mechanisms.

Identification of genes whose expression is altered by obesity throughout the arterial tree

We used next-generation RNA sequencing (RNA-Seq) technology on the whole transcriptome to identify genes whose expression is consistently affected by obesity across multiple arteries. Specifically, we examined transcriptional profiles of the iliac artery as well as the feed artery, first, second, and third branch order arterioles in the soleus, gastrocnemius, and diaphragm muscles from obese Otsuka Long-Evans Tokushima Fatty (OLETF) and lean Long-Evans Tokushima Otsuka (LETO) rats. Within the gastrocnemius and soleus muscles, the number of genes differentially expressed with obesity tended to increase with increasing branch order arteriole number (i.e., decreasing size of the artery). This trend was opposite in the diaphragm. We found a total of 15 genes that were consistently upregulated with obesity (MIS18A, CTRB1, FAM151B, FOLR2, PXMP4, OAS1B, SREBF2, KLRA17, SLC25A44, SNX10, SLFN3, MEF2BNB, IRF7, RAD23A, LGALS3BP) and five genes that were consistently downregulated with obesity (C2, GOLGA7, RIN3, PCP4, CYP2E1). A small fraction (∼9%) of the genes affected by obesity was modulated across all arteries examined. In conclusion, the present study identifies a select number of genes (i.e., 20 genes) whose expression is consistently altered throughout the arterial network in response to obesity and provides further insight into the heterogeneous vascular effects of obesity. Although there is no known direct function of the majority of 20 genes related to vascular health, the obesity-associated upregulation of SREBF2, LGALS3BP, IRF7, and FOLR2 across all arteries is suggestive of an unfavorable vascular phenotypic alteration with obesity. These data may serve as an important resource for identifying novel therapeutic targets against obesity-related vascular complications.

Transcriptome-wide RNA sequencing analysis of rat skeletal muscle feed arteries. I. Impact of obesity

We employed next-generation RNA sequencing (RNA-Seq) technology to determine the influence of obesity on global gene expression in skeletal muscle feed arteries. Transcriptional profiles of the gastrocnemius and soleus muscle feed arteries (GFA and SFA, respectively) and aortic endothelial cell-enriched samples from obese Otsuka Long-Evans Tokushima Fatty (OLETF) and lean Long-Evans Tokushima Otsuka (LETO) rats were examined. Obesity produced 282 upregulated and 133 downregulated genes in SFA and 163 upregulated and 77 downregulated genes in GFA [false discovery rate (FDR) < 10%] with an overlap of 93 genes between the arteries. In LETO rats, there were 89 upregulated and 114 downregulated genes in the GFA compared with the SFA. There were 244 upregulated and 275 downregulated genes in OLETF rats (FDR < 10%) in the GFA compared with the SFA, with an overlap of 76 differentially expressed genes common to both LETO and OLETF rats in both the GFA and SFA. A total of 396 transcripts were found to be differentially expressed between LETO and OLETF in aortic endothelial cell-enriched samples. Overall, we found 1) the existence of heterogeneity in the transcriptional profile of the SFA and GFA within healthy LETO rats, 2) that this between-vessel heterogeneity was markedly exacerbated in the hyperphagic, obese OLETF rat, and 3) a greater number of genes whose expression was altered by obesity in the SFA compared with the GFA. Also, results indicate that in OLETF rats the GFA takes on a relatively more proatherogenic phenotype compared with the SFA.

Transcriptome-wide RNA sequencing analysis of rat skeletal muscle feed arteries. II. Impact of exercise training in obesity

We employed next-generation RNA sequencing (RNA-Seq) technology to determine the extent to which exercise training alters global gene expression in skeletal muscle feed arteries and aortic endothelial cells of obese Otsuka Long-Evans Tokushima Fatty (OLETF) rats. Transcriptional profiles of the soleus and gastrocnemius muscle feed arteries (SFA and GFA, respectively) and aortic endothelial cell-enriched samples from rats that underwent an endurance exercise training program (EndEx; n = 12) or a interval sprint training program (IST; n = 12) or remained sedentary (Sed; n = 12) were examined. In response to EndEx, there were 39 upregulated (e.g., MANF) and 20 downregulated (e.g., ALOX15) genes in SFA and 1 upregulated (i.e., Wisp2) and 1 downregulated (i.e., Crem) gene in GFA [false discovery rate (FDR) < 10%]. In response to IST, there were 305 upregulated (e.g., MANF, HSPA12B) and 324 downregulated genes in SFA and 101 upregulated and 66 downregulated genes in GFA, with an overlap of 32 genes between arteries. Furthermore, in aortic endothelial cells, there were 183 upregulated (e.g., eNOS, SOD-3) and 141 downregulated (e.g., ATF3, Clec1b, npy, leptin) genes with EndEx and 71 upregulated and 69 downregulated genes with IST, with an overlap of 35 between exercise programs. Expression of only two genes (Tubb2b and Slc9a3r2) was altered (i.e., increased) by exercise in all three arteries. The finding that both EndEx and IST produced greater transcriptional changes in the SFA compared with the GFA is intriguing when considering the fact that treadmill bouts of exercise are associated with greater relative increases in blood flow to the gastrocnemius muscle compared with the soleus muscle.

Exercise training and vascular cell phenotype in a swine model of familial hypercholesterolaemia: conduit arteries and veins

NEW FINDINGS

What is the central question of this study? Does endurance exercise training cause anti-atherogenic effects on the endothelium in a swine model of familial hypercholesterolaemia (FH), and how are these effects distributed across veins, arteries and multiple vascular territories within each system? What is the main finding and its importance? Coronary artery endothelium-dependent vasomotor function was depressed in sedentary FH pigs compared with sedentary control animals, and exercise training did not change vasomotor function within FH. In systemic conduit arteries and veins, few effects of FH on endothelial cell protein expression were noted, including both pro- and anti-atherogenic changes. These findings suggest that exercise training does not produce a consistently improved endothelial cell phenotype in either coronary or systemic conduit vessels in this swine model of FH. Exercise training has emerged as an intervention for the primary and secondary prevention of coronary artery disease, but the mechanisms through which training reduces relative risk are not completely understood. The goal of this study was to investigate the impact of endurance exercise training on vasomotor function and vascular cell phenotype in coronary arteries and systemic conduit arteries and veins against a background of advanced atherosclerosis. We tested the hypothesis that exercise training restores endothelial vasomotor function and produces an anti-atherogenic endothelial and smooth muscle cell phenotype in familial hypercholesterolaemic (FH) swine. The study included 30 FH (15 exercised and 15 sedentary) and 13 non-FH control male castrated swine. The exercise-training intervention consisted of treadmill running 5 days per week for 16-20 weeks. Tissues sampled at sacrifice included vascular rings from the coronary circulation for vasomotor function experiments (dose-dependent bradykinin-induced vasorelaxation) and endothelial cells (ECs) from isolated segments of the thoracic aorta, the carotid, brachial, femoral and renal arteries, as well as each corresponding regionally associated vein, and from the abdominal vena cava, the right coronary and internal mammary arteries. Smooth muscle cells were sampled from the right coronary artery only. Vascular cell phenotype was assessed by immunoblotting for a host of both pro- and anti-atherogenic markers [e.g. endothelial nitric oxide synthase, p67phox, superoxide dismutase 1 (SOD1)]. Coronary artery endothelium-dependent vasomotor function was depressed in sedentary FH pigs compared with sedentary control pigs, and exercise training did not change vasomotor function within FH. In contrast, only scattered effects of FH on EC phenotype were noted across the vasculature, which included both pro- and anti-atherogenic changes in EC protein expression (e.g. increased endothelial nitric oxide synthase in carotid artery ECs, decreased p67phox in brachial artery ECs, but decreased expression of the antioxidant protein SOD1 in thoracic vena cava; all P < 0.05). In thoracic vena cava ECs, this deficit was corrected by exercise training, while no other effects of exercise were observed in conduit vessel EC phenotype. Thus, while exercise training abrogated the adverse effect of hypercholesterolaemia on thoracic vena cava SOD1 expression, it appears that exercise training does not produce a consistently improved EC phenotype in either coronary or systemic conduit vessels in this FH swine model.

Exercise training does not increase muscle FNDC5 protein or mRNA expression in pigs

BACKGROUND

Exercise training elevates circulating irisin and induces the expression of the FNDC5 gene in skeletal muscles of mice. Our objective was to determine whether exercise training also increases FNDC5 protein or mRNA expression in the skeletal muscles of pigs as well as plasma irisin.

METHODS

Castrated male pigs of the Rapacz familial hypercholesterolemic (FHM) strain and normal (Yucatan miniature) pigs were sacrificed after 16-20 weeks of exercise training. Samples of cardiac muscle, deltoid and triceps brachii muscle, subcutaneous and epicardial fat were obtained and FNDC5 mRNA, along with that of 6 other genes, was measured in all tissues of FHM pigs by reverse transcription polymerase chain reaction. FNDC protein in deltoid and triceps brachii was determined by Western blotting in both FHM and normal pigs. Citrate synthase activity was measured in the muscle samples of all pigs as an index of exercise training. Irisin was measured by an ELISA assay.

RESULTS

There was no statistically significant effect of exercise training on FNDC5 gene expression in epicardial or subcutaneous fat, deltoid muscle, triceps brachii muscle or heart muscle. Exercise-training elevated circulating levels of irisin in the FHM pigs and citrate synthase activity in deltoid and triceps brachii muscle. A similar increase in citrate synthase activity was seen in muscle extracts of exercise-trained normal pigs but there was no alteration in circulating irisin.

CONCLUSION

Exercise training in pigs does not increase FNDC5 mRNA or protein in the deltoid or triceps brachii of FHM or normal pigs while increasing circulating irisin only in the FHM pigs. These data indicate that the response to exercise training in normal pigs is not comparable to that seen in mice.

Vascular transcriptional alterations produced by juvenile obesity in Ossabaw swine

We adopted a transcriptome-wide microarray analysis approach to determine the extent to which vascular gene expression is altered as a result of juvenile obesity and identify obesity-responsive mRNAs. We examined transcriptional profiles in the left anterior descending coronary artery (LAD), perivascular fat adjacent to the LAD, and descending thoracic aorta between obese (n = 5) and lean (n = 6) juvenile Ossabaw pigs (age = 22 wk). Obesity was experimentally induced by feeding the animals a high-fat/high-fructose corn syrup/high-cholesterol diet for 16 wk. We found that expression of 189 vascular cell genes in the LAD and expression of 165 genes in the thoracic aorta were altered with juvenile obesity (false discovery rate ≤ 10%) with an overlap of only 28 genes between both arteries. Notably, a number of genes found to be markedly upregulated in the LAD of obese pigs are implicated in atherosclerosis, including ACP5, LYZ, CXCL14, APOE, PLA2G7, LGALS3, SPP1, ITGB2, CYBB, and P2RY12. Furthermore, pathway analysis revealed the induction of proinflammatory and pro-oxidant pathways with obesity primarily in the LAD. Gene expression in the LAD perivascular fat was minimally altered with juvenile obesity. Together, we provide new evidence that obesity produces artery-specific changes in pretranslational regulation with a clear upregulation of proatherogenic genes in the LAD. Our data may offer potential viable drug targets and mechanistic insights regarding the molecular precursors involved in the origins of overnutrition and obesity-associated vascular disease. In particular, our results suggest that the oxidized LDL/LOX-1/NF-κB signaling axis may be involved in the early initiation of a juvenile obesity-induced proatherogenic coronary artery phenotype.

Heterogeneity of endothelial cell phenotype within and amongst conduit vessels of the swine vasculature

The purpose of this study was to investigate the extent of endothelial cell phenotypic heterogeneity throughout the swine vasculature, with a focus on the conduit vessels of the arterial and venous circulations. We tested the hypothesis that atheroprone arteries exhibit higher expression of markers of inflammation and oxidative stress than do veins and atheroresistant arteries. The study sample included tissues from 79 castrated, male swine. Immediately after the animals were killed, endothelial cells were mechanically scraped from isolated segments of the thoracic and abdominal aorta, carotid, brachial, femoral and renal arteries, and the vein regionally associated with each of these vessels, as well as the internal mammary and right coronary arteries. Cells were also taken from two regions of the aortic arch contrasted by atheroprone versus atherosusceptible haemodynamics. Endothelial cell phenotype was assessed by either immunoblotting or quantitative real-time PCR for a host of both pro- and anti-atherogenic markers (e.g. endothelial nitric oxide synthase, p67phox, cyclo-oxygenase-1 and superoxide dismutase 1). Marked heterogeneity across the vasculature was observed in the expression of both pro- and anti-atherogenic markers, at both the protein and transcriptional levels. In particular, the coronary vascular endothelium expressed higher levels of the oxidative stress marker p67phox (P < 0.05 versus other arteries). In addition, differential expression of endothelial nitric oxide synthase and KLF4 was evident between atheroprone and atherosusceptible regions of the aorta, while expression of endothelial nitric oxide synthase, KLF2, KLF4 and cyclo-oxygenase-1 was lower in both areas of the aortic arch compared with the internal mammary artery. Conduit arteries typically expressed higher levels of both pro- and anti-atherogenic markers relative to their associated veins. We show, for the first time, that endothelial cell phenotype is variable within vessels, across six major vascular territories, and between the arterial and venous circulations. Importantly, even straight vessel segments from systemic conduit arteries (e.g. brachial and carotid arteries) exhibited regional phenotypic heterogeneity; a finding not expected on the basis of local haemodynamic forces alone.