Skip to main content

Sexual dimorphism in myocardial acylcarnitine and triglyceride metabolism

Abstract

Background

Cardiovascular disease is the leading cause of death among diabetic patients. Importantly, recent data highlight the apparent sexual dimorphism in the pathogenesis of cardiovascular disease in diabetics with respect to both frequency- and age-related risk factors. The disposition to cardiovascular disease among diabetic patients has been attributed, at least in part, to excess lipid supply to the heart culminating in lipotoxicity of the heart and downstream derangements. A confounding factor in obese animal models of diabetes is that increased peripheral lipid availability to the heart can induce cardio-metabolic remodeling independent of the underlying pathophysiology of diabetes, thus masking the diabetic phenotype. To that end, we hypothesized that the use of non-obese diabetic (NOD) animal models will reveal metabolic signatures of diabetes in a sex-specific manner.

Methods

To test this hypothesis, male and female NOD Goto-Kakizaki (GK) rats were used to assess the expression profile of 84 genes involved in lipid metabolism. In parallel, targeted lipidomics analysis was performed to characterize sex differences in homeostasis of non-esterified fatty acids (NEFA), acylcarnitines (AC), and triglycerides (TG).

Results

Our analysis revealed significant sex differences in the expression of a broad range of genes involved in transport, activation, and utilization of lipids. Furthermore, NOD male rats exhibited enhanced oxidative metabolism and accumulation of TG, whereas female NOD rats exhibited reduced TG content coupled with accumulation of AC species. Multi-dimensional statistical analysis identified saturated AC16:0, AC18:0, and AC20:0 as dominant metabolites in mediating sex differences in AC metabolism. Confocal microscopy of rat cardiomyocytes exposed to AC14:0, AC16:0, and AC18:0 confirmed induction of ROS with AC18:0 being more potent followed by AC14:0.

Conclusion

Overall, we demonstrate sex differences in myocardial AC and TG metabolism with implications for therapy and diagnosis of diabetic cardiovascular disease.

Background

Heart disease is the leading cause of death among diabetic patients, independent of macro- and micro- vascular diseases [13]. While diabetes increases the risk of cardiovascular disease in both males and females, recent data suggests that heart disease is fivefold more common in diabetic women compared with twofold in men [4]. For example, in the setting of myocardial infarction, diabetic women have a greater progression to heart failure and poorer outcome compared with diabetic males [5]. This “sexual dimorphism” is in contrast to the established notion of a “female advantage” in that premenopausal women in general are better protected from most common forms of cardiovascular disease compared to men [4, 69]. Thus, T2DM nullifies the female advantage in the prevalence of cardiovascular disease [4, 10].

Multiple lines of evidence suggest that the predisposition to heart failure among diabetic patients is a consequence of severe alterations in myocardial substrate metabolism [11, 12]. In particular, insulin resistance shifts the balance of substrate utilization such that the diabetic heart relies almost exclusively on fatty acids (FAs) for its energy needs [13]. The increased reliance on FA oxidation competes with glucose oxidation further exacerbating deficiencies in insulin signaling [14]. Excess lipid supply beyond the oxidative capacity of the heart results in accumulation of triglycerides (TG) and lipid intermediates, leading to lipotoxicity of the heart and downstream derangements, including generation of reactive oxygen species (ROS), inflammation, and cell death [1520].

On the one hand, the observed sexual dimorphism may be explained by the higher systemic availability of lipids [21, 22] resulting in increased myocardial FA metabolism. Indeed, it is well established that females have higher distribution of visceral fat with higher levels of circulating FA and TGs [22]. As an omnivore, the heart will necessarily utilize circulating substrates in proportion to their availability in blood. To that end, in vivo metabolic imaging studies in diabetic humans have demonstrated that, in fact, female hearts exhibit higher flux of FAs compared to age-matched males [23, 24], attributed to higher systemic availability of FA. On the other hand, however, there may be intrinsic sex differences in cardiac lipid metabolism beyond increased peripheral supply of lipids to the heart that may account for the sexual dimorphism. One confounding factor in obese and dietary animal models of diabetes is that increased peripheral lipid availability to the heart can induce cardio-metabolic remodeling independent of the underlying pathophysiology of diabetes, thus masking the diabetic phenotype. In this regard, it is worth noting that, in fact, 20 % of the population is non-obese diabetic (NOD) [25, 26]. This suggests that studies utilizing NOD humans and/or non-obese animal models of diabetes may unmask mechanisms relevant to pathophysiology of diabetes.

In this work, we sought to determine whether there are intrinsic sex differences in myocardial FA metabolism in the setting of diabetes, independent of systemic effects seen in obesity or dietary animal models of diabetes. We previously demonstrated that NOD Goto-Kakizaki (GK) rats exhibit a genomic disposition to increased FA metabolism [27]. We hypothesized that utilizing both sexes of this model will reveal metabolic signatures of diabetes independent of obesity, in a sex-specific manner. Quantitative expression profiling of genes involved in FA, acylcarnitines (AC), and TG metabolism as well as targeted lipidomic profiling on cardiac tissue extracts was performed to characterize genomic-metabolic variations between the groups. Finally, multi-dimensional statistical analysis was performed to identify and confirm signatures of sex differences in metabolism that can induce ROS.

Methods

Animal model and study design

All experiments were conducted according to the protocol approved by the animal experiment committee at Washington University School of Medicine at Saint Louis (IACUC Animal Welfare Assurance # A-3381-01). NOD and Wistar rats (Charles River, USA) were maintained under standard housing conditions and fed with commercial pelleted chow (Harlan, USA) until ready for the study at age of 24 weeks. A balanced factorial design at two levels was used to investigate sex and phenotype differences and possible interactions between these two factors. The following four animal groups were tested: (1) Wistar male, (2) Wistar female, (3) GK male, and (4) GK female with N = 4–6 per group. Blood samples were taken via the tail vein from anesthetized animals (overnight fast, typically 10 h) in each group and immediately centrifuged. The serum fractions were assayed for TG, non-esterified fatty acids (NEFA), glucose, and insulin levels. All animals were sacrificed at the end of the study, and their hearts were snap-frozen for subsequent analyses. Lipidomic analysis was performed on individual tissue extracts (N = 6 per group), while gene expression studies was performed on representative samples (N = 4 per group), randomly chosen from each group. FA and TG were extracted from 100 mg of pulverized cardiac tissue samples using a chloroform/methanol mixture. Commercial kits for FA (Wako Diagnostics, VA) and TG analyses (Infinity Triglycerides Reagent, Thermo Scientific, USA) were used for both plasma and tissue samples.

RNA isolation and real-time RT-PCR

Quantitative real-time RT-PCR (qPCR) was performed as described previously [27]. The RNeasy lipid tissue extraction kit (Qiagen, MD) was used to extract total RNA from pulverized heart tissue. RNA concentration, purity, and integrity were determined by spectrophotometric analysis and Agilent Bio-analyzer. qPCR was conducted in a 384-well plate array format for 84 genes involved in FA metabolism (Cat No: PARN- 0007E, SA Biosciences, MD). From the panel of housekeeping genes, lactate dehydrogenase A (ldhA) was chosen for normalizing the data as it displayed no variation among test groups. For assessing ROS response in AC-treated cardiomyocytes, PrimeTime® primers for Sod1, Sod2, and Cat were sourced from IDT DNA technologies (Coralville, IA).

Mass spectrometric detection of NEFAs, ACs, and TGs

Mass spectrometric analysis was performed at the Diabetic Cardiovascular Disease Center (DCDC) Lipidomics Core Facility using established methodology. Briefly, a modified Bligh-Dyer extraction [28] was used to extract NEFAs and TGs, and a protein precipitation method was used to extract ACs from the heart homogenates in the presence of internal standards (d4-NEFA16:0, TG17:1-17:1-17:1, and d3-AC16:0). To improve the sensitivity of MS, NEFAs were derivatized with dimethylaminopropylamine, dimethylaminopyridine, and 1-ethyl-3-[3-dimethylaminopropyl] carbodiimide hydrochloride into amides, and ACs were derivatized with methanol and acetyl chloride into methyl esters. Sample analysis was performed with a Shimadzu 10A HPLC system coupled to a TSQ Quantum Ultra triple quadrupole mass spectrometer operated in SRM mode under ESI(+).

Echocardiography

Echocardiography was performed using non-invasive ultrasound imaging with the Vevo2100 Ultrasound System (Visual Sonics Inc., Toronto, ON, Canada) at ages 23–24 weeks as described previously [29]. Briefly, rats were anesthetized with 1 % Isoflurane and secured to the imaging platform. Complete two-dimensional, M-mode, and Doppler examinations using a 21-MHz transducer were performed to quantify left ventricular structure as well as diastolic and systolic function. To account for differences in body weight, all of the structural dimensions were indexed to body weight.

Confocal microscopy to assess induction of ROS by ACs

Immortalized neonatal rat cardiomyocytes were obtained from Applied Biological Materials, Inc. Canada (Catalog No. T0015) and maintained in PriGrow III media (ABM Good, Canada). The method for microscopic determination of ROS has been described previously [30, 31] and was used here with minor modifications appropriate for this study. Briefly, cells were seeded at 1 × 105 per 35-mm glass-bottomed μdish (Catalog no. 81158 iBidi, Germany) and incubated overnight. The cells were serum starved for 8 h and then treated with 10 μM CM-H2DCFDA (Molecular Probes, C6827), a chloromethyl (CM) derivative of H2DCFDA, in PBS at 37 °C for 30 min. The treatment was followed by two quick washes using ice cold PBS. Based on previous work [30], cells were incubated with 25 μM of each of the following ACs for 3 h at 37 °C: palmitoyl-l-carnitine (AC16:0) (P1645, Sigma-Aldrich), stearoyl-l-carnitine (AC18:0) (50494231, Fisher Scientific), and myristoyl-l-carnitine (AC14:0) (41570, Crystal Chem, IL). The 3-h incubation period was selected to remain within the stability range of the fluorophore per manufacturer’s recommendation. Cells were then examined using a Nikon A1Rsi + laser scanning confocal microscope at ×40 magnification using laser excitation at 488 and emission 543 nm at the Bakewell Neuroimaging Laboratory at Washington University. ROIs were defined using NIS-Elements (Nikon, USA) on 30–35 images by filtering background and saturated intensities as well as filtering for objects/debris that were less than 5 μm in size as rat cardiomyocytes are at least 10 μm in size [32]. The sum intensity of ROIs for each image from the test groups were compared with the control by t test.

Data and statistical analysis

All data sets from the two level factorial designs were first examined to verify equal variance among test groups using residual analysis, and if necessary, the data were log2 transformed. The data were analyzed using a two-way ANOVA with interaction model (sex, phenotype, sex*phenotype). ANOVA was followed with group contrast comparisons using Fisher’s least squares difference to determine significance. Student’s t test was used to determine significant statistical differences in fluorescence intensities obtained from our microscopic analysis of genes Sod1, Sod2, and Cat between control and AC-treated groups.

Sex differences in the distribution of lipid species

We also sought to assess whether the fraction (as percent) of species enriched/depleted in each lipid category (i.e., FA, TG, and AC) was significantly different within and between sexes. This analysis ignored the statistical findings of the ANOVA for individual species; rather, the number of species enriched or depleted in the respective category was counted. The significance of the overall fold change in the number of species enriched or depleted (GK to Wistar) profile was determined by comparing to a binomial distribution with the assumption that if the variable had no effect, then the number of species enriched or depleted would be random with a binomial distribution with mean = 0.5.

Sex differences in metabolic pathways

Genes were categorized into the metabolic functional groups. Because the interpretation of the ANOVA interaction term for the genes was more nuanced than for the metabolites (i.e., many interactions were due to differences in magnitude rather than just sign), a rank-sum test, which considers both sign and magnitude, was used to determine significant differences in metabolic pathways between the sexes. The GK/Wistar fold change sign and magnitude was used as the test parameter and sex was the categorical variable. The data were analyzed for low-expressing genes using a criterion that the Ct count must not be greater than 2 standard deviations from the average of all genes. Therefore, any gene with a Ct count >28 was considered non-expressing.

Principal component analysis

Principal component analysis (PCA) was performed on the gene and metabolite data separately using the correlation matrix of the standardized data. Only those genes or metabolites that had at least one significant contrast comparison in the ANOVA analysis were included in the PCA. As in above, low-expressing genes were omitted from the analysis. Radar plots were used to present the PCA results. Partek Genomics Suite (St. Louis, MO) was used for ANOVA, STATA (College Station, TX) for the binomial distribution tests, and JMP Genomics for the PCA. The data are presented as mean ± standard error of the mean (SEM) and a P value ≤0.05 was considered significant. While the statistics may have been performed on transformed data, the results are graphed in their original metric as a convenience to the reader.

Results

Baseline characteristics

We found no significant differences in serum FA level among the four groups. Serum TG levels were significantly elevated in NOD rats with males and females exhibiting 1.38- and 2.2-fold increase over corresponding Wistar rats, respectively (Table 1). Overall, these baseline characteristics are in agreement with previous reports on GK rats [33, 34].

Table 1 Baseline measurements

Echocardiograph: anatomical and functional parameters

Weight-indexed measures of LVPW, LVID, and LVM were significantly higher in NOD rats compared to corresponding Wistar control rats (Table 2). Functionally, NOD rats exhibited a significant reduction in parameters of early component mitral annular diastolic motion (E′) and LV filling velocity in the atrial component (A) compared to Wistars, as well as both phenotypic and sex differences in early left ventricular filling velocity (E) (Table 2).

Table 2 Echocardiograph measurements

Sex differences in expression profile of genes involved in lipid metabolism

Figure 1a summarizes the ANOVA results for the differential expression of genes categorized by their metabolic functionality. The number of genes modulated between NOD and Wistar by sex is shown for each functional group. Overall, the pathway analysis indicates a significant increase in the expression of genes involved in beta-oxidation in NOD male rats relative to female rats, concomitant with reduced expression of genes involved FA activation and AC metabolism in NOD female rats. Other significant sex differences in pathway expression include intracellular FA trafficking and triacylglycerol metabolism. To characterize which genes dominate the variance among the groups, we performed PCA to reduce the dimensionality of the data.

Fig. 1
figure 1

Sex differences in expression profile of genes involved in lipid metabolism. a Number of genes modulated between NOD and Wistar rats by sex from ANOVA model. A section sign (§) indicates that pathways differ between sexes by the rank-sum test (P ≤ 0.05). b 3D score plot of first three principal components (PCs) for genes from principal component analysis. c Shows the radar plot of loading factors of PCs for gene analysis. The purple line and circles represent PC1, the brown line and circles represent PC2, and green line and circles represent PC3. While absolute loading values are used to plot, the closed and open circles represent (+) and (–) loading values, respectively. P < 0.05 is the significance level for all tests

Figure 1b depicts the segregation of the four test groups by PCA. The first three principal components (PCs) accounted for 85 % of the variability. Figure 1c is a radar plot of absolute gene loading values greater than 0.5 for each PCs. The genes were sorted by decreasing order of the absolute magnitude of their primary component. The PCs are associated with general fold change patterns that can be discerned from examination of Additional file 1: Table S1. PC1 accounts for 57.4 % of the variation in the data and is associated with genes that exhibit sex differences and a significant interaction term in the ANOVA model. Those genes with a positive PC1 value (PC1pos) are downregulated in females compared to males, while those with a negative value (PC1neg) are upregulated in females. Further, the interaction is far more prevalent in genes associated with PC1neg in that the sex difference fold changes are higher in Wistar than in the NOD phenotype. Generally, genes associated with PC2(+) are elevated in NOD rats, while those with PC2(−) are lower in NOD rats. PC3 has too few genes to ascribe a pattern. Overall, the loading values can be used to correlate the contribution of a gene/pathway to the overall phenotype of sex differences in metabolism.

Sex differences in NEFA homeostasis

The enrichment profile of NEFAs obtained by LC-MS/MS is summarized in Fig. 2a (Additional file 2: Figure S1 depicts male and female data separately). Our results reveal that of the NEFAs profiled, every measured NEFA species is significantly enriched in male NOD rats compared to male Wistar rats in agreement with overall tissue NEFA concentration summarized in Table 1. In females, however, only NEFAs 20:1 and 22:1 are significantly depleted in NOD, while NEFAs 20:3 and 22:4 are significantly enriched in NODs. Based on the binomial distribution analysis, there is a significant sex difference in the distribution of NEFAs with male NOD exhibiting significant enrichment of lipid species (Fig. 2b). The data in Fig. 2c suggests that the concentration of both long and very long NEFAs is higher in NOD compared to Wistar rats, but the differences only reach statistical significance in males. Analysis of NEFA profile based on degree of saturation reveals that saturated, mono-, and poly-unsaturated fatty acids accumulate to a higher extent in NOD compared with Wistar rats, but the differences are only significant in males (Fig. 2d). Furthermore, NOD rats overall exhibit significant accumulation of mono- and poly-unsaturated fatty acids compared to Wistar rats.

Fig. 2
figure 2

Sex differences in NEFA homeostasis. a Profile of NEFA in cardiac tissue, as measured by LC-MS/MS, presented as fold change between NOD and Wistar rats by sex. An asterisk (*) indicates that the fold change is statistically significant. A section sign (§) indicates that the ANOVA interaction term is significant, and therefore, the fold changes in males and females also are. A dagger (†) indicates that in the ANOVA model, the overall sex difference is statistically significant. b Distribution of enriched and depleted of NEFA moities. The asterisk (*) indicates that the distribution is significantly different than a binomial distribution with mean = 0.5. The dagger (†) indicates that the distribution in males and females differ significantly. c, d The distribution of FAs based on their chain-length and degree of saturation, respectively. Data is presented as mean ± SEM (N = 6). An asterisk (*) indicates that the difference in adjacent bars is statistically significant. A section sign (§) indicates that the ANOVA interaction term is significant. P < 0.05 is the significance level for all tests

Sex differences in TG homeostasis

Analysis of the individual TG species levels by LC-MS/MS indicates that all measured TG species were depleted in NOD females compared to Wistar females (Fig. 3a) (Additional file 3: Figure S2 depicts male and female data separately). Moreover, the binomial distribution analysis suggests that the number of species depleted in females was statistically significant, compared to males (see “Methods” section) (Fig. 3b). Consistent with these data, the spectrophotometric analysis of total TG in cardiac tissue suggests that female NOD rats exhibited significantly reduced TG concentrations than the Wistar female, while male NOD rats exhibited a higher average TG content. Although the difference in males is not significant, the significant interaction suggests that the differences in total TG accumulation between NOD and Wistar rats differs significantly in males and females (Fig. 3c).

Fig. 3
figure 3

Sex differences in TG homeostasis. a Profile of TG species in cardiac tissue, as measured by LC-MS/MS, presented as fold change between NOD and Wistar rats by sex. An asterisk (*) indicates that the fold change is statistically significant. A section sign (§) indicates that the ANOVA interaction term is significant, and therefore, the fold changes in males and females are also significant. A dagger (†) indicates that in the ANOVA model, the overall sex difference is statistically significant. b Distribution of enriched and depleted of TG moities. The asterisk (*) indicates that the distribution is significantly different than a binomial distribution with mean = 0.5. The dagger (†) indicates that the distribution in males and females differ significantly. c The total amount of TG present in the cardiac tissue as measured by spectrophotometric analysis, mean ± SEM (N = 6). P < 0.05 is the significance level for all tests

Sex differences in AC homeostasis

The LC-MS/MS metabolic profile of ACs revealed significant sex differences in the ACs profiled (Fig. 4a) (Additional file 4: Figure S3 depicts male and female data separately). Moreover, our data suggests sex differences in abundance of ACs with female rats exhibiting significantly higher abundance of AC (Fig. 4b) and species (Fig. 4c) compared to male rats. Since tissue ACs are elevated in female rats, we sought to determine whether differences in serum TG or NEFA levels could account for the higher AC levels in tissue. To test this possibility, we normalized the AC concentration to the sum of the serum TG and NEFA concentrations for each rat and ran the ANOVA model on this modified data. Our data (not shown) indicate sex differences in AC levels independent of serum TG and NEFA. Thus, the elevated serum levels cannot sufficiently account for the increased AC levels in female.

Fig. 4
figure 4

Sex differences in AC homeostasis. a Profile of AC in cardiac tissue, as measured by LC-MS/MS, presented as fold change between NOD and Wistar rats by sex. An asterisk (*) indicates that the fold change is statistically significant. A dagger (†) indicates that in the ANOVA model, the overall sex difference is statistically significant. b Net ACs level in the tissue, identified by summing up individual AC species. The dagger (†) indicates that in the ANOVA model the overall sex difference is statistically significant. c Distribution of enriched and depleted of AC moities. The asterisk (*) indicates that the distribution is significantly different than a binomial distribution with mean = 0.5. The dagger (†) indicates that the distribution in males and females differ significantly. d Pooled distribution of ACs based on their chain length. e Pooled distribution of ACs based degree of saturation. Data in d and e are presented as mean ± SEM (N = 6). The asterisk (*) indicates that comparison between adjacent bars is significant. The section sign (§) indicates that the ANOVA interaction term is significant. A dagger (†) indicates that in the ANOVA model, the overall sex difference is statistically significant

The AC concentration data was binned by medium, long, and very long chains, and the results plotted in Fig. 4d. The results indicate that AC concentrations are lower in males for all chain lengths and that for the medium chain lengths, the concentrations are statistically lower in NOD females than in Wistar females. These observations are consistent with the information in Fig. 4a, b. It is also observed from Fig. 4a that, for the short chain AC, 4/4 are elevated in females in Wistar rats, while only the longest of the short chain AC is statistically elevated in NOD females compared to males. In the long chain AC, however, 5/10 are statistically higher in female compared to male NOD rats, while only 1/10 is statistically higher in Wistar rats. Finally, partitioning of AC by degree of saturation suggests that female rats exhibit higher levels of saturated and mono-unsaturated AC (Fig. 4e).

Integrated lipidomics analysis of NEFA, AC, and TG

Figure 5a reveals the segregation of the lipidomics data among the four test groups based on PCA analysis. Figure 5b summarizes the results of the PCA and is a radar plot of absolute gene loading values greater than 0.5 for each PC. The general fold change pattern associated with each PC can be discerned from examination of Additional file 5: Table S2. PC1 accounts for 52.4 % of the variation in the data and is associated almost exclusively with TG that are elevated in female Wistar over male Wistars. There are two AC associated with PC1, and both are elevated in females irrespective of phenotype. There are also three NEFAs associated with PC1 that exhibit interactions both with respect to female vs male by phenotype and GK vs Wistar by sex.

Fig. 5
figure 5

Integrated lipidomics analysis of NEFA, AC, and TG and induction of ROS by ACs. a 3D score plot of first three principal components (PCs) for metabolites obtained from PCA. b Radar plot of loading factors of principal components for lipidomic analysis. The purple line and circles represent PC1, the brown line and circles represent PC2, the green line and circles represent PC3, and the yellow line and circle represents PC4. All loadings were positive and hence are depicted using closed circles. c Depicted are representative confocal microscopy images of cardiomyocytes incubated with ethanol (control) and 25 μM of AC14:0, AC16:0, and AC18:0 (bar = 50 μM). d The sum intensities of each image from the test groups were compared with untreated cells. Averages of the image intensities, mean + SEM are shown. e Relative expression levels of genes Sod1, Sod2, and Cat in cardiomyocytes compared to the untreated control cells. Average CT of b-actin and LdhA was used to normalize the data. An asterisk (*) indicates that the comparison identified by the bar is significant at P < 0.05

PC2 accounts for 20.7 % of the variation and is associated primarily with NEFA. PC2 species also have significant interaction terms in their ANOVA models. NEFA are elevated in females among Wistar rats but are depleted in females among NOD rats. Among males, NEFA are elevated in NOD rats, while generally the same among females. PC3 is composed primarily of ACs, which are elevated in females irrespective of phenotype. One AC, 22:4, is also elevated in NOD over Wistar overall. PC4 is also exclusively AC of 20 carbon length but differs from PC3 in that the interaction term is significant in the ANOVA model. AC 20:0 is elevated in females among NOD rats but is unchanged with sex in Wistar, and it is depleted in NOD among male rats.

Of the FA examined, PCA identified FA 18:1, 20:3, and 22:4 as having the most influence in dictating sex differences in FA metabolism, whereas AC16:0 and AC18:0 in PC3 as well as AC20:0 in PC4 were identified as contributing to sex differences in AC metabolism (Fig. 5b). To test the feasibility of the former set of ACs to induce ROS, cardiomyocytes were exposed to 25 μM of each AC14:0, AC16:0, and AC18:0 ACs. Quantitative analysis of confocal microscopy images (Fig. 5c) suggests that the relative efficiency of ROS induction was AC18:0 > AC14:0 > AC16:0 (Fig. 5d). Accordingly, the expression levels of Sod1, Sod2, and Cat genes encoding for key enzymes involved in ROS scavenging were significantly higher in AC-treated cardiomyocytes (Fig. 5e). In agreement with quantitative confocal microscopy images, the magnitude of differential expression was highest in AC18:0-treated cardiomyocytes compared to untreated controls (Fig. 5e).

Discussion

The overarching objective of this work was to identify sex differences in lipid metabolism in a NOD animal model of T2DM. This effort was partly motivated by the lack of adequate animal models to study sex differences in T2DM [35]. Since gene expression alone may be insufficient to address the question of sex differences in metabolism, we pursed a genomic-metabolomic approach to characterize sex differences in metabolism where changes in metabolite species are the ultimate end-product of gene and protein expression differences. Of note, it is well established that women generally have a higher percentage of body fat than men [36], resulting in increased levels of circulating FA and TGs. Since increased peripheral lipid availability to the heart, as is the case in obesity, can induce cardio-metabolic remodeling independent of the underlying pathophysiology of diabetes, we hypothesized that the use of NOD animal models of T2DM will reveal metabolic signatures of diabetes, in a sex-specific manner.

Unlike obese animal models of T2DM, NOD rats did not display significant differences in circulating NEFA levels although they displayed slightly elevated TG, consistent with previous work on NOD rats [3740]. In cardiac tissue, NOD rats exhibited elevated steady state levels of NEFAs, which was particularly higher in male NOD rats. Expression analysis of gene arrays involved in FA metabolism revealed marked differences in mitochondrial beta-oxidation, most notably preference towards FA oxidation in male NOD rats in contrast to down-regulation of gene involved in FA oxidation in female rats (Fig. 1). In parallel, we observed significant sex differences in the expression of genes involved in FA activation and AC metabolism. Moreover, multi-dimensional PCA analysis identified significant differences in the expression levels of a broad set of genes regulating multiple pathways including energy sensing, FA catabolism, acylcarnitine transport and activation, and that of TG biosynthesis and degradation. To assess downstream implications of the differences in gene expression, we performed lipidomic profiling of NEFA, TG, and AC in cardiac tissues.

Our lipidomics data suggests sex differences in TG homeostasis. NOD female rats exhibited significantly lower total TG concentrations than female Wistar rats, and virtually all the TG species assayed were lower in NOD females compared with age-matched Wistar rats. While TGs have long been proposed to adversely affect cardiac function and physiology [16, 18], recent data suggests that TG may also play a cardio-protective role by absorbing excess lipid supply to the heart [41, 42]. The observation that female NOD hearts display significantly lower levels of TG may suggest that FAs were channeled away from TG synthesis or that TGs were converted into other lipid products such as DAG and ceramides.

In addition to sex differences in TG metabolism, our lipidomic analysis also revealed sex differences in AC metabolism with female rats exhibiting higher concentration of ACs than male rats. Long-chain ACs have long ago been documented as arrhythmogenic affecting myocyte electrophysiology [4345] and influencing recovery following ischemic injury [46, 47]. However, recent mechanistic [48, 49] and metabolomic [5053] studies using diabetic animal models have identified intra-mitochondrial disturbances where lipid overload resulted in the accumulation of AC and incomplete FA oxidation [54]. The role of AC in T2DM has been highlighted by recent reports suggesting that serum levels of ACs are altered in pre-diabetic conditions [54, 55] and by reports identifying serum AC species indicative of enhanced risk to T2DM and insulin resistance [31, 5659]. Importantly, accumulation of ACs has been shown to interfere with insulin signaling [48, 6063], activate pro-inflammatory signaling pathways [30], induce cellular uncoupling [64], and generate cytoplasmic ROS [30] in murine monocyte/macrophages.

In the current work, PC3 in Fig. 5c identified AC16:0 and AC18:0 as having the most variability in dictating sex differences in AC metabolism. To test the feasibility of AC to induce ROS in cardiomyocytes, rat cardiomyocytes were exposed to 25 μM of AC16:0, AC18:0, and AC14:0 (the latter based on literature precedent [30]). Our results suggest that AC18 induced over threefold higher production of ROS than untreated cardiomyocytes, followed by AC14:0. In parallel to increased ROS as a response to AC stimuli, expression analysis of genes with established roles in ROS scavenging confirmed activation of ROS mitigating pathways. These observations demonstrate, for the first time, that ACs can induce production of ROS in cardiomyocytes. Further studies are required to fully characterize the specificity and selectivity of ACs in inducing cardiac ROS, insulin resistance, and downstream derangements in cardiomyocyte health and function in a sex-specific manner.

Conclusions

In summary, these observations represent the first such report of sex differences in AC and TG metabolism as a potential metabolic signature of diabetes, which may contribute to downstream derangements in cardiac function in a sex-specific manner. Our analysis revealed sex differences in the expression of a broad range of genes involved in the metabolism of FAs, ACs, and TGs. Furthermore, integrated lipidomics analysis indicated that whereas NOD male rats exhibit enhanced oxidative metabolism activity and accumulation of TG species, female NOD rats tend to exhibit enhanced accumulation of AC species and reduced TG compared to male rats. Finally, we have demonstrated that ACs can induce production of ROS in cardiomyocytes, providing a new dimension for the observed sexual dimorphism in diabetic cardiovascular disease. Some of the potential implications include the following: from a therapeutic point, therapies that minimize accumulation of ACs and potentially other by-products of beta-oxidation may be of interest in treating females. From an imaging standpoint, our data suggests that currently available metabolic tracers, in particular for positron emission tomography (PET), may be unable to discern intrinsic (as opposed to systemic) sex differences in lipid metabolism, suggesting that new panel of metabolic PET tracers may be needed to fully interrogate sex differences in cardiac metabolism and downstream effects. Thus, overall, our data highlight broad implications for therapy and diagnosis of T2DM.

References

  1. Fang ZY, Prins JB, Marwick TH. Diabetic cardiomyopathy: evidence, mechanisms, and therapeutic implications. Endocr Rev. 2004;25:543–67.

    Article  CAS  PubMed  Google Scholar 

  2. Stamler J, Vaccaro O, Neaton JD, Wentworth D. Diabetes, other risk factors, and 12-yr cardiovascular mortality for men screened in the Multiple Risk Factor Intervention Trial. Diabetes Care. 1993;16:434–44.

    Article  CAS  PubMed  Google Scholar 

  3. Wilson PW. Diabetes mellitus and coronary heart disease. Endocrinol Metab Clin North Am. 2001;30:857–81.

    Article  CAS  PubMed  Google Scholar 

  4. Ren J, Ceylan-Isik AF. Diabetic cardiomyopathy: do women differ from men? Endocrine. 2004;25:73–83.

    Article  CAS  PubMed  Google Scholar 

  5. Savage MP, Krolewski AS, Kenien GG, Lebeis MP, Christlieb AR, Lewis SM. Acute myocardial infarction in diabetes mellitus and significance of congestive heart failure as a prognostic factor. Am J Cardiol. 1988;62:665–9.

    Article  CAS  PubMed  Google Scholar 

  6. Skafar DF, Xu R, Morales J, Ram J, Sowers JR. Clinical review 91: female sex hormones and cardiovascular disease in women. J Clin Endocrinol Metab. 1997;82:3913–8.

    CAS  PubMed  Google Scholar 

  7. Ren J. Cardiac health and diabetes mellitus in women: problems and prospects. Minerva Cardioangiol. 2006;54:289–309.

    CAS  PubMed  Google Scholar 

  8. Kaseta JR, Skafar DF, Ram JL, Jacober SJ, Sowers JR. Cardiovascular disease in the diabetic woman. J Clin Endocrinol Metab. 1999;84:1835–8.

    Article  CAS  PubMed  Google Scholar 

  9. Collins P, Stevenson JC, Mosca L. Spotlight on gender. Cardiovasc Res. 2002;53:535–7.

    Article  CAS  PubMed  Google Scholar 

  10. Sowers JR. Diabetes mellitus and cardiovascular disease in women. Arch Intern Med. 1998;158:617–21.

    Article  CAS  PubMed  Google Scholar 

  11. Lopaschuk GD, Belke DD, Gamble J, Itoi T, Schonekess BO. Regulation of fatty acid oxidation in the mammalian heart in health and disease. Biochim Biophys Acta. 1994;1213:263–76.

    Article  CAS  PubMed  Google Scholar 

  12. Lopaschuk GD. Metabolic abnormalities in the diabetic heart. Heart Fail Rev. 2002;7:149–59.

    Article  CAS  PubMed  Google Scholar 

  13. Saddik M, Lopaschuk GD. Triacylglycerol turnover in isolated working hearts of acutely diabetic rats. Can J Physiol Pharmacol. 1994;72:1110–9.

    Article  CAS  PubMed  Google Scholar 

  14. Randle PJ, Garland PB, Hales CN, Newsholme EA. The glucose fatty-acid cycle. Its role in insulin sensitivity and the metabolic disturbances of diabetes mellitus. Lancet. 1963;1:785–9.

    Article  CAS  PubMed  Google Scholar 

  15. Drosatos K, Schulze PC. Cardiac lipotoxicity: molecular pathways and therapeutic implications. Curr Heart Fail Rep. 2013;10:109–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Wende AR, Symons JD, Abel ED. Mechanisms of lipotoxicity in the cardiovascular system. Curr Hypertens Rep. 2012;14:517–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Goldberg IJ, Trent CM, Schulze PC. Lipid metabolism and toxicity in the heart. Cell Metab. 2012;15:805–12.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Wende AR, Abel ED. Lipotoxicity in the heart. Biochim Biophys Acta. 2010;1801:311–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Borradaile NM, Schaffer JE. Lipotoxicity in the heart. Curr Hypertens Rep. 2005;7:412–7.

    Article  CAS  PubMed  Google Scholar 

  20. Schaffer JE. Lipotoxicity: when tissues overeat. Curr Opin Lipidol. 2003;14:281–7.

    Article  CAS  PubMed  Google Scholar 

  21. Kotani K, Tokunaga K, Fujioka S, Kobatake T, Keno Y, Yoshida S, Shimomura I, Tarui S, Matsuzawa Y. Sexual dimorphism of age-related-changes in whole-body fat distribution in the obese. Int J Obes (Lond). 1994;18:207–12.

    CAS  Google Scholar 

  22. Williams CM. Lipid metabolism in women. Proc Nutr Soc. 2004;63:153–60.

    Article  CAS  PubMed  Google Scholar 

  23. Lyons MR, Peterson LR, McGill JB, Herrero P, Coggan AR, Saeed IM, Recklein C, Schechtman KB, Gropler RJ. Impact of sex on the heart’s metabolic and functional responses to diabetic therapies. Am J Physiol Heart Circ Physiol. 2013;305:H1584–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Iozzo P, Lautamaki R, Borra R, Lehto HR, Bucci M, Viljanen A, Parkka J, Lepomaki V, Maggio R, Parkkola R. Contribution of glucose tolerance and gender to cardiac adiposity. J Clin Endocrinol Metab. 2009;94:4472–82.

    Article  CAS  PubMed  Google Scholar 

  25. Dalton M, Cameron AJ, Zimmet PZ, Shaw JE, Jolley D, Dunstan DW, et al. Waist circumference, waist-hip ratio and body mass index and their correlation with cardiovascular disease risk factors in Australian adults. J Intern Med. 2003;254(6):555–63.

    Article  CAS  PubMed  Google Scholar 

  26. Garancini MP, Calori G, Ruotolo G, Manara E, Izzo A, Ebbli E, et al. Prevalence of NIDDM and impaired glucose-tolerance in Italy—an OGTT-based population study. Diabetologia. 1995;38(3):306–13.

    Article  CAS  PubMed  Google Scholar 

  27. Devanathan S, Nemanich ST, Kovacs A, Fettig N, Gropler RJ, Shoghi KI. Genomic and metabolic disposition of non-obese type 2 diabetic rats to increased myocardial fatty acid metabolism. PloS One. 2013. In Press.

  28. Bligh EG, Dyer WJ. A rapid method of total lipid extraction and purification. Can J Biochem Physiol. 1959;37:911–7.

    Article  CAS  PubMed  Google Scholar 

  29. Shoghi KI, Gropler RJ, Sharp T, Herrero P, Fettig N, Su Y, Mitra MS, Kovacs A, Finck BN, Welch MJ. Time course of alterations in myocardial glucose utilization in the Zucker diabetic fatty rat with correlation to gene expression of glucose transporters: a small-animal PET investigation. J Nucl Med. 2008;49:1320–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Rutkowsky JM, Knotts TA, Ono-Moore KD, McCoin CS, Huang S, Schneider D, Singh S, Adams SH, Hwang DH. Acylcarnitines activate proinflammatory signaling pathways. Am J Physiol Endocrinol Metab. 2014;306:E1378–1387.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Aguer C, McCoin CS, Knotts TA, Thrush AB, Ono-Moore K, McPherson R, Dent R, Hwang DH, Adams SH, Harper ME. Acylcarnitines: potential implications for skeletal muscle insulin resistance. FASEB J. 2015;29:336–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Egorova MV, Afanas'ev SA, Popov SV. A simple method for isolation of cardiomyocytes from adult rat heart. Bull Exp Biol Med. 2005;140:370–3.

    Article  CAS  PubMed  Google Scholar 

  33. Almon RR, DuBois DC, Lai W, Xue B, Nie J, Jusko WJ. Gene expression analysis of hepatic roles in cause and development of diabetes in Goto-Kakizaki rats. J Endocrinol. 2009;200(3):331–46.

    Article  CAS  PubMed  Google Scholar 

  34. Chandler MP, Morgan EE, McElfresh TA, Kung TA, Rennison JH, Hoit BD, Young ME. Heart failure progression is accelerated following myocardial infarction in type 2 diabetic rats. Am J Physiol Heart Circ Physiol. 2007;293:H1609–1616.

    Article  CAS  PubMed  Google Scholar 

  35. Franconi F, Seghieri G, Canu S, Straface E, Campesi I, Malorni W. Are the available experimental models of type 2 diabetes appropriate for a gender perspective? Pharmacol Res. 2008;57:6–18.

    Article  PubMed  Google Scholar 

  36. Power ML, Schulkin J. Sex differences in fat storage, fat metabolism, and the health risks from obesity: possible evolutionary origins. Br J Nutr. 2008;99:931–40.

    Article  CAS  PubMed  Google Scholar 

  37. O'Rourke CM, Davis JA, Saltiel AR, Cornicelli JA. Metabolic effects of troglitazone in the Goto-Kakizaki rat, a non-obese and normolipidemic rodent model of non-insulin-dependent diabetes mellitus. Metabolism. 1997;46:192–8.

    Article  PubMed  Google Scholar 

  38. Desrois M, Sidell RJ, Gauguier D, Davey CL, Radda GK, Clarke K. Gender differences in hypertrophy, insulin resistance and ischemic injury in the aging type 2 diabetic rat heart. J Mol Cell Cardiol. 2004;37:547–55.

    Article  CAS  PubMed  Google Scholar 

  39. Nie J, Xue B, Sukumaran S, Jusko WJ, DuBois DC, Almon RR. Differential muscle gene expression as a function of disease progression in Goto-Kakizaki diabetic rats. Mol Cell Endocrinol. 2011;338(1-2):10–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Ao Y, Toy N, Song MK, Go VLW, Yang H. Altered glucose and insulin responses to brain medullary thyrotropin-releasing hormone (TRH)-induced autonomic activation in type 2 diabetic Goto-Kakizaki rats. Endocrinology. 2005;146(12):5425–32.

    Article  CAS  PubMed  Google Scholar 

  41. Holland WL, Knotts TA, Chavez JA, Wang LP, Hoehn KL, Summers SA. Lipid mediators of insulin resistance. Nutr Rev. 2007;65:S39–46.

    Article  PubMed  Google Scholar 

  42. Listenberger LL, Han XL, Lewis SE, Cases S, Farese RV, Ory DS, Schaffer JE. Triglyceride accumulation protects against fatty acid-induced lipotoxicity. Proc Natl Acad Sci U S A. 2003;100:3077–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Ferro F, Ouille A, Tran TA, Fontanaud P, Bois P, Babuty D, Labarthe F, Le Guennec JY. Long-chain acylcarnitines regulate the hERG channel. PLoS One. 2012;7:e41686.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Wu J, Corr PB. Influence of long-chain acylcarnitines on voltage-dependent calcium current in adult ventricular myocytes. Am J Physiol. 1992;263:H410–417.

    CAS  PubMed  Google Scholar 

  45. Yamada KA, Kanter EM, Newatia A. Long-chain acylcarnitine induces Ca2+ efflux from the sarcoplasmic reticulum. J Cardiovasc Pharmacol. 2000;36:14–21.

    Article  CAS  PubMed  Google Scholar 

  46. Corr PB, Heathers GP, Yamada KA. Mechanisms contributing to the arrhythmogenic influences of alpha 1-adrenergic stimulation in the ischemic heart. Am J Med. 1989;87:19S–25S.

    Article  CAS  PubMed  Google Scholar 

  47. Ford DA, Han XL, Horner CC, Gross RW. Accumulation of unsaturated acylcarnitine molecular species during acute myocardial ischemia: metabolic compartmentalization of products of fatty acyl chain elongation in the acylcarnitine pool. Biochemistry. 1996;35:7903–9.

    Article  CAS  PubMed  Google Scholar 

  48. Koves TR, Ussher JR, Noland RC, Slentz D, Mosedale M, Ilkayeva O, Bain J, Stevens R, Dyck JR, Newgard CB, et al. Mitochondrial overload and incomplete fatty acid oxidation contribute to skeletal muscle insulin resistance. Cell Metab. 2008;7:45–56.

    Article  CAS  PubMed  Google Scholar 

  49. Petersen KF, Dufour S, Morino K, Yoo PS, Cline GW, Shulman GI. Reversal of muscle insulin resistance by weight reduction in young, lean, insulin-resistant offspring of parents with type 2 diabetes. Proc Natl Acad Sci U S A. 2012;109:8236–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Huffman KM, Shah SH, Stevens RD, Bain JR, Muehlbauer M, Slentz CA, Tanner CJ, Kuchibhatla M, Houmard JA, Newgard CB, Kraus WE. Relationships between circulating metabolic intermediates and insulin action in overweight to obese, inactive men and women. Diabetes Care. 2009;32:1678–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Huffman KM, Slentz CA, Bateman LA, Thompson D, Muehlbauer MJ, Bain JR, Stevens RD, Wenner BR, Kraus VB, Newgard CB, Kraus WE. Exercise-induced changes in metabolic intermediates, hormones, and inflammatory markers associated with improvements in insulin sensitivity. Diabetes Care. 2011;34:174–6.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Redman LM, Huffman KM, Landerman LR, Pieper CF, Bain JR, Muehlbauer MJ, Stevens RD, Wenner BR, Kraus VB, Newgard CB. Effect of caloric restriction with and without exercise on metabolic intermediates in nonobese men and women. J Clin Endocrinol Metab. 2011;96:E312–321.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Shah SH, Hauser ER, Bain JR, Muehlbauer MJ, Haynes C, Stevens RD, Wenner BR, Dowdy ZE, Granger CB, Ginsburg GS, et al. High heritability of metabolomic profiles in families burdened with premature cardiovascular disease. Mol Syst Biol. 2009;5:258.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Schooneman MG, Vaz FM, Houten SM, Soeters MR. Acylcarnitines: reflecting or inflicting insulin resistance? Diabetes. 2013;62:1–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Mai M, Tonjes A, Kovacs P, Stumvoll M, Fiedler GM, Leichtle AB. Serum levels of acylcarnitines are altered in prediabetic conditions. PLoS One. 2013;8:e82459.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Floegel A, Stefan N, Yu Z, Muhlenbruch K, Drogan D, Joost HG, Fritsche A, Haring HU, Hrabe de Angelis M, Peters A, et al. Identification of serum metabolites associated with risk of type 2 diabetes using a targeted metabolomic approach. Diabetes. 2013;62:639–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Bene J, Marton M, Mohas M, Bagosi Z, Bujtor Z, Oroszlan T, Gasztonyi B, Wittmann I, Melegh B. Similarities in serum acylcarnitine patterns in type 1 and type 2 diabetes mellitus and in metabolic syndrome. Ann Nutr Metab. 2013;62:80–5.

    Article  CAS  PubMed  Google Scholar 

  58. Sampey BP, Freemerman AJ, Zhang J, Kuan PF, Galanko JA, O'Connell TM, Ilkayeva OR, Muehlbauer MJ, Stevens RD, Newgard CB, et al. Metabolomic profiling reveals mitochondrial-derived lipid biomarkers that drive obesity-associated inflammation. PLoS One. 2012;7:e38812.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Mihalik SJ, Goodpaster BH, Kelley DE, Chace DH, Vockley J, Toledo FG, DeLany JP. Increased levels of plasma acylcarnitines in obesity and type 2 diabetes and identification of a marker of glucolipotoxicity. Obesity (Silver Spring). 2010;18:1695–700.

    Article  CAS  Google Scholar 

  60. Godin DV, Lopaschuk GD, McNeill JH. Subcellular myocardial abnormalities in experimental diabetes: role of long-chain acylcarnitines. Can J Cardiol. 1986;2:222–9.

    CAS  PubMed  Google Scholar 

  61. Rodriguez-Gutierrez R, Lavalle-Gonzalez FJ, Martinez-Garza LE, Landeros-Olvera E, Lopez-Alvarenga JC, Torres-Sepulveda MR, Gonzalez-Gonzalez JG, Mancillas-Adame LG, Salazar-Gonzalez B, Villarreal-Perez JZ. Impact of an exercise program on acylcarnitines in obesity: a prospective controlled study. J Int Soc Sports Nutr. 2012;9:22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Muoio DM, Neufer PD. Lipid-induced mitochondrial stress and insulin action in muscle. Cell Metab. 2012;15:595–605.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Su X, Han X, Mancuso DJ, Abendschein DR, Gross RW. Accumulation of long-chain acylcarnitine and 3-hydroxy acylcarnitine molecular species in diabetic myocardium: identification of alterations in mitochondrial fatty acid processing in diabetic myocardium by shotgun lipidomics. Biochemistry. 2005;44:5234–45.

    Article  CAS  PubMed  Google Scholar 

  64. Yamada KA, McHowat J, Yan GX, Donahue K, Peirick J, Kleber AG, Corr PB. Cellular uncoupling induced by accumulation of long-chain acylcarnitine during ischemia. Circ Res. 1994;74:83–95.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank the Preclinical PET/CT imaging facility at Washington University School of Medicine for the technical support in facilitating the studies. The authors thank Dr. Daniel Ory and the staff of the Metabolomics Core facility at the Diabetic Cardiovascular Disease Center (DCDC) for the lipidomics profiling, and Dr. Attila Kovacks at the Mouse Cardiovascular Phenotyping Core for echocardiogram analysis.

Funding

This work was funded by NIH/NIDDK grant R01DK085298 and NIH/NIBIB grant R01EB012284. The work utilized core-services provided by the Diabetes Research Center (DRC) grant P30DK020579, the Nutrition Obesity Research Center (NORC) grant P30DK056341, and the Center for Cardiovascular Research (CCR) supported by NIH Shared Instrumentation Grant S10RR026949.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Kooresh I. Shoghi.

Additional information

Competing interest

The authors declare that they have no competing interests.

Authors’ contributions

SD performed in vitro experiments, interpreted the data, and drafted the manuscript. TDW analyzed and interpreted the data and drafted the manuscript. NF performed the in vivo experiments. RJG interpreted and drafted the manuscript. SN analyzed the data. KIS conceived the study, analyzed and interpreted the data, and wrote the manuscript. All authors read and approved the final manuscript.

Additional files

Additional file 1:

Table S1. Summary statistics, ANOVA and principal component analysis results for genes. (XLSX 31 kb)

Additional file 2:

Figure S1. Concentrations of non-esterified fatty acid moieties in NOD and Wistar rats by sex. (PPTX 163 kb)

Additional file 3:

Figure S2. Concentrations of triglyceride moieties in NOD and Wistar rats by sex. (PPTX 226 kb)

Additional file 4:

Figure S3. Concentrations of acylcarnitine moieties in NOD and Wistar rats by sex. (PPTX 142 kb)

Additional file 5:

Table S2. Summary statistics and ANOVA and principal component analysis results for metabolites. (XLSX 43 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Devanathan, S., Whitehead, T.D., Fettig, N. et al. Sexual dimorphism in myocardial acylcarnitine and triglyceride metabolism. Biol Sex Differ 7, 25 (2016). https://doi.org/10.1186/s13293-016-0077-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13293-016-0077-7

Keywords