Selection-, age-, and exercise-dependence of skeletal muscle gene expression patterns in a rat model of metabolic fitness.

Intrinsic aerobic exercise capacity can influence many complex traits including obesity and aging. To study this connection we established two rat lines by divergent selection of untrained aerobic capacity. After 32 generations the high capacity runners (HCR) and low capacity runners (LCR) differed in endurance running distance and body fat, blood glucose, other health indicators, and natural life span. To understand the interplay among genetic differences, chronological age, and acute exercise we performed microarray-based gene expression analyses in skeletal muscle with a 2×2×2 design to simultaneously compare HCR and LCR, old and young animals, and rest and exhaustion. Transcripts for mitochondrial function are expressed higher in HCRs than LCRs at both rest and exhaustion and for both age groups. Expression of cell adhesion and extracellular matrix genes tend to decrease with age. This and other age effects are more prominent in LCRs than HCRs, suggesting that HCRs have a slower aging process and this may be partly due to their better metabolic health. Strenuous exercise mainly affects transcription regulation and cellular response. The effects of any one factor often depend on the other two. For example, there are ∼140 and ∼110 line-exercise "interacting" genes for old and young animals, respectively. Many genes highlighted in our study are consistent with prior reports, but many others are novel. The gene- and pathway-level statistics for the main effects, either overall or stratified, and for all possible interactions, represent a rich reference dataset for understanding the interdependence among lines, aging, and exercise.

Increased glucose metabolism and glycerolipid formation by fatty acids and GPR40 receptor signaling underlies the fatty acid potentiation of insulin secretion.

Acute fatty acid (FA) exposure potentiates glucose-stimulated insulin secretion in ß cells through metabolic and receptor-mediated effects. We assessed the effect of fatty acids on the dynamics of the metabolome in INS-1 cells following exposure to [U-(13)C]glucose to assess flux through metabolic pathways. Metabolite profiling showed a fatty acid-induced increase in long chain acyl-CoAs that were rapidly esterified with glucose-derived glycerol-3-phosphate to form lysophosphatidic acid, mono- and diacylglycerols, and other glycerolipids, some implicated in augmenting insulin secretion. Glucose utilization and glycolytic flux increased, along with a reduction in the NADH/NAD(+) ratio, presumably by an increase in conversion of dihydroxyacetone phosphate to glycerol-3-phosphate. The fatty acid-induced increase in glycolysis also resulted in increases in tricarboxylic cycle flux and oxygen consumption. Inhibition of fatty acid activation of FFAR1/GPR40 by an antagonist decreased glycerolipid formation, attenuated fatty acid increases in glucose oxidation, and increased mitochondrial FA flux, as evidenced by increased acylcarnitine levels. Conversely, FFAR1/GPR40 activation in the presence of low FA increased flux into glycerolipids and enhanced glucose oxidation. These results suggest that, by remodeling glucose and lipid metabolism, fatty acid significantly increases the formation of both lipid- and TCA cycle-derived intermediates that augment insulin secretion, increasing our understanding of mechanisms underlying ß cell insulin secretion.

Sex modulates the diet-induced changes to the plasma lipidome in a rat model of cardiorespiratory fitness.

OBJECTIVE: Individuals with higher intrinsic cardiorespiratory fitness (CRF) experience decreased rates of cardiometabolic disease and mortality, and high CRF is associated with increased utilization of fatty acids (FAs) for energy. Studies suggest a complex relationship between CRF, diet, and sex with health outcomes, but this interaction is understudied. We hypothesized that FA utilization differences by fitness and sex could be detected in the plasma metabolome when rats or humans were fed a high carbohydrate (HC) or high fat (HF) diet. METHODS: Male and female rats selectively bred for low (LCR) and high (HCR) CRF were fed a chow diet or a sucrose-free HF (45 % fat) or HC (10 % fat) diet. Plasma samples were collected at days 0, 3, and 14. Human plasma data was collected from male and female participants who were randomized into a HC or HF diet for 21 days. Samples were analyzed using liquid chromatography-mass spectrometry and regression statistics were used to quantify the effect of diet, CRF, and sex on the lipidome. RESULTS: In rats, the baseline lipidome is more significantly influenced by sex than by CRF, especially as elevated diglycerides, triglycerides, phosphatidylcholines, and lysophosphatidylcholines in males. A dynamic response to diet was observed 3 days after diet, but after 14 days of either diet, the lipidome was modulated by sex with a larger effect size than by diet. Data from the human study also suggests a sex-dependent response to diet with opposite directionality of affect compared to rats, highlighting species-dependent responses to dietary intervention.

Intrinsic cardiorespiratory fitness modulates clinical and molecular response to caloric restriction.

OBJECTIVE: Caloric restriction (CR) is one extrinsic intervention that can improve metabolic health, and it shares many phenotypical parallels with intrinsic high cardiorespiratory fitness (CRF), including reduced adiposity, increased cardiometabolic health, and increased longevity. CRF is a highly heritable trait in humans and has been established in a genetic rat model selectively bred for high (HCR) and low (LCR) CRF, in which the HCR live longer and have reduced body weight compared to LCR. This study addresses whether the inherited high CRF phenotype occurs through similar mechanisms by which CR promotes health and longevity. METHODS: We compared HCR and LCR male rats fed ad libitum (AL) or calorically restricted (CR) for multiple physiological, metabolic, and molecular traits, including running capacity at 2, 8, and 12 months; per-hour metabolic cage activity over daily cycles at 6 and 12 months; and plasma lipidomics, liver and muscle transcriptomics, and body composition after 12 months of treatment. RESULTS: LCR-CR developed a physiological profile that mirrors the high-CRF phenotype in HCR-AL, including reduced adiposity and increased insulin sensitivity. HCR show higher spontaneous activity than LCR. Temporal modeling of hourly energy expenditure (EE) dynamics during the day, adjusted for body weight and hourly activity levels, suggest that CR has an EE-suppressing effect, and high-CRF has an EE-enhancing effect. Pathway analysis of gene transcripts indicates that HCR and LCR both show a response to CR that is similar in the muscle and different in the liver. CONCLUSIONS: CR provides LCR a health-associated positive effect on physiological parameters that strongly resemble HCR. Analysis of whole-body EE and transcriptomics suggests that HCR and LCR show line-dependent responses to CR that may be accreditable to difference in genetic makeup. The results do not preclude the possibility that CRF and CR pathways may converge.

Skeletal muscle gene expression in response to resistance exercise: sex specific regulation.

BACKGROUND: The molecular mechanisms underlying the sex differences in human muscle morphology and function remain to be elucidated. The sex differences in the skeletal muscle transcriptome in both the resting state and following anabolic stimuli, such as resistance exercise (RE), might provide insight to the contributors of sexual dimorphism of muscle phenotypes. We used microarrays to profile the transcriptome of the biceps brachii of young men and women who underwent an acute unilateral RE session following 12 weeks of progressive training. Bilateral muscle biopsies were obtained either at an early (4 h post-exercise) or late recovery (24 h post-exercise) time point. Muscle transcription profiles were compared in the resting state between men (n = 6) and women (n = 8), and in response to acute RE in trained exercised vs. untrained non-exercised control muscle for each sex and time point separately (4 h post-exercise, n = 3 males, n = 4 females; 24 h post-exercise, n = 3 males, n = 4 females). A logistic regression-based method (LRpath), following Bayesian moderated t-statistic (IMBT), was used to test gene functional groups and biological pathways enriched with differentially expressed genes. RESULTS: This investigation identified extensive sex differences present in the muscle transcriptome at baseline and following acute RE. In the resting state, female muscle had a greater transcript abundance of genes involved in fatty acid oxidation and gene transcription/translation processes. After strenuous RE at the same relative intensity, the time course of the transcriptional modulation was sex-dependent. Males experienced prolonged changes while females exhibited a rapid restoration. Most of the biological processes involved in the RE-induced transcriptional regulation were observed in both males and females, but sex specificity was suggested for several signaling pathways including activation of notch signaling and TGF-beta signaling in females. Sex differences in skeletal muscle transcriptional regulation might implicate a mechanism behind disproportional muscle growth in males as compared with female counterparts after RE training at the same relative intensity. CONCLUSIONS: Sex differences exist in skeletal muscle gene transcription both at rest and following acute RE, suggesting that sex is a significant modifier of the transcriptional regulation in skeletal muscle. The findings from the present study provide insight into the molecular mechanisms for sex differences in muscle phenotypes and for muscle transcriptional regulation associated with training adaptations to resistance exercise.

Extended regions of suspected mis-assembly in the rat reference genome.

We performed whole-genome sequencing for eight inbred rat strains commonly used in genetic mapping studies. They are the founders of the NIH heterogeneous stock (HS) outbred colony. We provide their sequences and variant calls to the rat genomics community. When analyzing the variant calls we identified regions with unusually high levels of heterozygosity. These regions are consistent across the eight inbred strains, including Brown Norway, which is the basis of the rat reference genome. These regions show higher read depths than other regions in the genome and contain higher rates of apparent tri-allelic variant sites. The evidence suggests that these regions may correspond to duplicated segments that were incorrectly overlaid as a single segment in the reference genome. We provide masks for these regions of suspected mis-assembly as a resource for the community to flag potentially false interpretations of mapping or functional results.

Expression of macrophage genes within skeletal muscle correlates inversely with adiposity and insulin resistance in humans.

Local inflammation in obese adipose tissue has been shown to contribute to insulin resistance; however, the role of macrophage infiltration within skeletal muscle is still debatable. This study aimed to evaluate the association of skeletal muscle macrophage gene expression with adiposity levels and insulin sensitivity in obese patients. Twenty-two nondiabetic obese patients and 23 healthy lean controls were included. Obese patients underwent a 3-month weight loss intervention. Macrophage gene expression in skeletal muscle (quantitative real-time polymerase chain reaction), body composition (dual-energy X-ray absorptiometry), and insulin sensitivity (homeostatic model assessment (HOMA) and oral glucose tolerance test) were compared between groups and their associations were analyzed. To validate skeletal muscle findings, we repeated the analyses with macrophage gene expression in adipose tissue. Expression levels of macrophage genes (CD68, CD11b, CD206, CD16, CD40, and CD163) were lower in skeletal muscle tissue of obese versus lean participants. Macrophage gene expression was also found to be inversely associated with adiposity, fasting insulin, and HOMA (r = -0.4 ∼ -0.6, p < 0.05), as well as positively associated with insulin sensitivity (r = 0.4 ∼ 0.8, p < 0.05). On the other hand, adipose tissue macrophage gene expression showed higher levels in obese versus lean participants, presenting a positive association with adiposity levels. Macrophage gene expression, in both skeletal and adipose tissue samples, was only minimally affected by the weight loss intervention. In contrast with the established positive relationship between adiposity and macrophage gene expression, an unexpected inverse correlation between these 2 variables was observed in skeletal muscle tissue. Additionally, muscle macrophage gene expression was inversely correlated with insulin resistance.

Islet microvasculature in islet hyperplasia and failure in a model of type 2 diabetes.

Gene expression profiling of islets from pre-diabetic male Zucker diabetic fatty (ZDF) rats showed increased expression of hypoxia-related genes, prompting investigation of the vascular integrity of the islets. The islet microvasculature was increased approximately twofold in young male ZDF rats by both morphometric analysis and quantifying mRNA levels of endothelial markers. ZDF rats at 12 weeks of age showed a significant reduction in the number of endothelial cells, which was prevented by pretreatment with pioglitazone. Light and electron microscopy of normoglycemic 7-week-old ZDF rats showed thickened endothelial cells with loss of endothelial fenestrations. By 12 weeks of age, there was disruption of the endothelium and intra-islet hemorrhage. Islets from 7- and 12-week-old ZDF rats showed an approximate three- and twofold increase in vascular endothelial growth factor (VEGF)-A mRNA and VEGF protein secretion, respectively, compared with lean controls. Thrombospondin-1 mRNA increased in 7- and 12-week-old rats by 2- and 10-fold, respectively, and was reduced by 50% in 12-week-old rats pretreated with pioglitazone. Islets from young male control rats induced migration of endothelial cells in a collagen matrix only after pretreatment with matrix metalloproteinase (MMP)-9. Islets from 7-week-old ZDF rats showed a fivefold increase in migration score compared with wild-type controls, even without MMP-9 treatment. Islets from 15-week-old ZDF rats did not induce migration; rather, they caused a significant rounding up of the duct-derived cells, suggesting a toxic effect. These data suggest that in the ZDF rat model of type 2 diabetes, an inability of the islet to maintain vascular integrity may contribute to beta-cell failure.

Genetic analysis of a rat model of aerobic capacity and metabolic fitness.

Aerobic capacity is a strong predictor of all-cause mortality and can influence many complex traits. To explore the biological basis underlying this connection, we developed via artificial selection two rat lines that diverge for intrinsic (i.e. inborn) aerobic capacity and differ in risk for complex disease traits. Here we conduct the first in-depth pedigree and molecular genetic analysis of these lines, the high capacity runners (HCR) and low capacity runners (LCR). Our results show that both HCR and LCR lines maintain considerable narrow-sense heritability (h(2)) for the running capacity phenotype over 28 generations (h(2) = 0.47 ± 0.02 and 0.43 ± 0.02, respectively). To minimize inbreeding, the lines were maintained by rotational mating. Pedigree records predict that the inbreeding coefficient increases at a rate of <1% per generation, ~37-38% slower than expected for random mating. Genome-wide 10K SNP genotype data for generations 5, 14, and 26 demonstrate substantial genomic evolution: between-line differentiation increased progressively, while within-line diversity deceased. Genome-wide average heterozygosity decreased at a rate of <1% per generation, consistent with pedigree-based predictions and confirming the effectiveness of rotational breeding. Linkage disequilibrium index r(2) decreases to 0.3 at ~3 Mb, suggesting that the resolution for mapping quantitative trait loci (QTL) can be as high as 2-3 cM. To establish a test population for QTL mapping, we conducted an HCR-LCR intercross. Running capacity of the F1 population (n=176) was intermediate of the HCR and LCR parentals (28 pairs); and the F2 population (n=645) showed a wider range of phenotypic distribution. Importantly, heritability in the F0-F2 pedigree remained high (h(2)~0.6). These results suggest that the HCR-LCR lines can serve as a valuable system for studying genomic evolution, and a powerful resource for mapping QTL for a host of characters relevant to human health.

Islet hypertrophy following pancreatic disruption of Smad4 signaling.

To investigate the role of transforming growth factor (TGF)-beta family signaling in the adult pancreas, a transgenic mouse (E-dnSmad4) was created that expresses a dominant-negative Smad4 protein driven by a fragment of the elastase promoter. Although E-dnSmad4 mice have normal growth, pancreas weight, and pancreatic exocrine and ductal histology, beginning at 4-6 wk of age, E-dnSmad4 mice show an age-dependent increase in the size of islets. In parallel, an expanded population of replicating cells expressing the E-dnSmad4 transgene is found in the stroma between the enlarged islets and pancreatic ducts. Despite the marked enlargement, E-dnSmad4 islets contain normal ratios and spatial organization of endocrine cell subtypes and have normal glucose homeostasis. Replication of cells derived from primary duct cultures of wild-type mice, but not E-dnSmad4 mice, was inhibited by the addition of TGF-beta family proteins, demonstrating a cell-autonomous effect of the transgene. These data show that, in the adult pancreas, TGF-beta family signaling plays a role in islet size by regulating the growth of a pluripotent progenitor cell residing in the periductal stroma of the pancreas.