Sex-dimorphism in Cardiac Nutrigenomics: effect of Trans fat and/or Monosodium Glutamate consumption
BMC Genomics volume 12, Article number: 555 (2011)
A paucity of information on biological sex-specific differences in cardiac gene expression in response to diet has prompted this present nutrigenomics investigation.
Sexual dimorphism exists in the physiological and transcriptional response to diet, particularly in response to high-fat feeding. Consumption of Trans-fatty acids (TFA) has been linked to substantially increased risk of heart disease, in which sexual dimorphism is apparent, with males suffering a higher disease rate. Impairment of the cardiovascular system has been noted in animals exposed to Monosodium Glutamate (MSG) during the neonatal period, and sexual dimorphism in the growth axis of MSG-treated animals has previously been noted. Processed foods may contain both TFA and MSG.
We examined physiological differences and changes in gene expression in response to TFA and/or MSG consumption compared to a control diet, in male and female C57BL/6J mice.
Heart and % body weight increases were greater in TFA-fed mice, who also exhibited dyslipidemia (P < 0.05). Hearts from MSG-fed females weighed less than males (P < 0.05). 2-factor ANOVA indicated that the TFA diet induced over twice as many cardiac differentially expressed genes (DEGs) in males compared to females (P < 0.001); and 4 times as many male DEGs were downregulated including Gata4, Mef2d and Srebf2. Enrichment of functional Gene Ontology (GO) categories were related to transcription, phosphorylation and anatomic structure (P < 0.01). A number of genes were upregulated in males and downregulated in females, including pro-apoptotic histone deacetylase-2 (HDAC2). Sexual dimorphism was also observed in cardiac transcription from MSG-fed animals, with both sexes upregulating approximately 100 DEGs exhibiting sex-specific differences in GO categories. A comparison of cardiac gene expression between all diet combinations together identified a subset of 111 DEGs significant only in males, 64 DEGs significant in females only, and 74 transcripts identified as differentially expressed in response to dietary manipulation in both sexes.
Our model identified major changes in the cardiac transcriptional profile of TFA and/or MSG-fed mice compared to controls, which was reflected by significant differences in the physiological profile within the 4 diet groups. Identification of sexual dimorphism in cardiac transcription may provide the basis for sex-specific medicine in the future.
There is growing evidence that sexual dimorphism exists in the physiological response to diet and other environmental factors , including differences in insulin and leptin sensitivity [2, 3], particularly in response to high fat diets . The heart requires high levels of ATP in order to maintain contractility, and this essential energy requirement is usually met by the oxidation of fatty acids (FAs), which are the most important source of energy for the adult heart. FAs have also been shown to be involved in signal transduction and the regulation of gene transcription ; however excessive FA uptake may cause detrimental lipotoxicity, resulting in deregulated cardiac function . Several key studies have recently shown that considerable sex-dimorphism exists in both rodents and humans ; and have established sex-specific pathways in cardiac transcription in response to both pressure overload  and dilated cardiomyopathy .However there is a relative scarcity of information concerning sex-specific effects of diet on cardiac gene expression.
Dietary intake of added fats and oils has increased by over 60% during the past 35 years , and use of hydrogenated oils and shortenings in order to prolong shelf-life has resulted in an increase in trans fatty acid (TFA) consumption, which now accounts for between 1.7 - 8% of the global dietary fat intake (0.6 - 3% of total energy intake) . TFA consumption has been linked to the increased risk of obesity , inflammation , coronary heart disease (CHD)  and type 2 diabetes . Unlike saturated fat, food nutrition labeling does not set a recommended intake limit for Trans-fat, however one meta-analysis study found that a 2% increase in energy intake from TFA was associated with a 23% increase in the incidence of CHD . Concern over the deleterious effects of dietary Trans- fat has prompted several countries to regulate the industrial production ; and retrospective meta-analysis of observational controlled dietary trials and prospective cohort studies suggests that replacement of Trans-fat with nonhydrogenated oils would result in significant reduction of the incidence of CHD . Gender differences in the prevalence of CHD are well documented , and many factors appear to be responsible including lipid profile  and genetic susceptibility .
Neonatal administration of high doses of the food flavor enhancer Monosodium Glutamate (MSG) has also been shown to result in obesity  together with stunted heart growth and hypoplasticity in rodents . Furthermore, the cardiovascular system of these animals have been shown to respond deficiently to certain challenges, including an attenuated blood pressure response to the systemic injection of nitric oxide synthase inhibitors, angiotensin II and other vasoactive compounds . Injections of high doses of MSG (4 g/Kg body weight) cause ablation of the Arcuate Nucleus via glutamate-induced degeneration of those areas of the immature brain which are insufficiently protected by a mature blood-brain barrier, including the hypothalamic area . This suggested that the hypothalamus may play an important role in the growth of the heart. It is now apparent that maternal administration of MSG can penetrate the placental barrier and distribute to the embryonic tissues of the fetus. Oral administration of MSG and 3H-labeled Glutamate (3H- Glu) as a tracer to pregnant mice resulted in marked elevations of 3H- Glu in the placenta and in fetal brain, liver and kidney . MSG consumption has increased globally in recent years, with recent estimations of the current average daily intake believed to be up to 10 g/day . Despite the widespread consumption of this common food flavor enhancer, which often occurs in processed and packaged foods in combination with Trans-fat, very few studies have addressed the effect of dietary Trans-fat and/or MSG consumption on cardiac gene expression. Furthermore, although sexual dimorphism at the growth hormone axis has been demonstrated in MSG-treated rodents , little is known about sex-specific gene expression in response to exposure to MSG.
The aim of this work was to establish the effect of dietary TFA and/or MSG on cardiac gene expression in an in vivo animal model, and to examine sex-specific differences as well as commonalities of gene expression in response to diet. The amount of oral MSG used in this study reflects current consumption levels [26, 28] and is 30 - 40 times less that the level previously reported to induce neuronal damage when injected neonatally [21, 23]. Since it has previously been ascertained that MSG excitotoxicity occurs only when the blood brain barrier is vulnerable, for example neonatally [24, 25], and because developmental programming of cardiac gene expression may be affected by maternal nutritional balance , we bred our study animals from females previously established on these diets for 3 weeks prior to mating, in accordance with our previous studies . Exposure to these diets occurred throughout the study, in order to mimic as closely as possible the situation that occurs in humans. The results of our microarray analysis revealed considerable sexual dimorphism in response to the dietary interventions that we tested, with males appearing to exhibit more differential gene expression compared to females.
Consumption of Trans-fat and/or MSG results in sex-specific changes in physiological characteristics
Figure 1A shows our experimental design and microarray analysis. There were no significant differences in food and water intake between the 4 different diet groups or between the sexes (data not shown). Averaged daily MSG intake was 91.21 ± 4.63 and 99.04 ± 2.5 mg/Kg body weight in males and females respectively. Trans-fat (TFA) feeding resulted in elevated weight gain (6-32 weeks) and greater heart weights in both sexes (Table 1, P < 0.001), and females weighed less than males (P ≤ 0.05). Additionally, male mice in the MSG and Trans-fat diet groups had significantly greater weight gain and heart : body weight ratios compared to females (P < 0.001). Analysis of significant sex-specific differences in response to diet indicated that females in the MSG and the MSG+TFA diet group had smaller hearts and lower fasting HDL-C levels compared to males. Females in the TFA diet group also had higher fasting triglyceride and HOMA-IR values compared to female control diet mice. Conversely TFA feeding elevated male T-CHOL and HDL-C levels. There was no apparent difference in the histological appearance of cardiac tissues from the four diet groups (online Additional File 1).
Sexual dimorphism in cardiac gene expression
We used Affymetrix Mouse Gene 1.0 ST expression arrays to determine differences in global cardiac gene/Expressed Sequence Tags (ESTs) in response to diet, in both males and females. Using a False Discovery Rate with a significance level set at 0.05, we identified two subsets of genes/ESTs that were differentially expressed between males (n = 1803) and females (n = 1985) at various stringency levels (absolute fold change > 2.0 > 2.5, > 3.0, adjusted probability value <0.05, Figure 1B and Additional File 2 Table S2a and S2b). To better understand the relevance of these sex-dimorphic genes we used a program based on Gene Ontology (GO) categories of biological processes, molecular function and cellular component. Additional File 3, Table S3 lists the main functional GO categories of transcripts that were differentially expressed between the sexes with a stringency of ± 1.5-fold change in expression (P < 0.01). DEGs were further characterized as being relatively up-regulated either in males or in females. Increased transcription of male genes in the "Biological processes" category included genes involved in cellular macromolecule metabolism, anatomical structure development and blood vessel development. These genes included Notch4, interstitial matrix gene Col3a1 (type 3 collagen alpha 1), fibulin 1 (Fbln1) and integrin αV (Itgav). The predominant male-specific KEGG pathways included the regulation of actin cytoskeleton and antigen processing pathways (Additional File 3, Table S3). Conversely female-specific transcription were enriched in the categories of cellular and metabolic processes, including mitochondrial carbonic anhydrase 5b (Car5B), glutathione S-transferase 4 (Gstm4) and magnesium transporter 1 (Magt1). One KEGG-pathway for Cysteine and methionine metabolism was also upregulated.
Chromosomal localization of differentially expressed cardiac genes
To investigate a possible non-random chromosomal distribution of the DEGs, we used the Fisher exact test to compare the frequencies of significant genes on each chromosome against the total number of genes interrogated on each chromosome. Chromosomes 8,11 and Y were enriched for male-biased genes whereas chromosomes 3,4,12 and X contained more female-biased genes (Figure 2, P < 0.05). Notably we found four genes encoded on the Y chromosome which were male-specific (Sry, Ddx3y, Eif2s3y, and Uty), and 86 female-biased X-linked genes including Ace2, dystrophin, and several transcription factors such as Tceal1, Tceal8 and Eif2s3x.
Trans-fat feeding induces sex-specific alterations in cardiac gene expression
In order to identify the subset of genes and ESTs with a fold difference of ± 1.5 or more within male or female comparisons for mice fed the TFA diet versus control, we used the results from a 2-way ANOVA model described in the methods section. We identified 361 DEGs (148 upregulated and 213 downregulated) that were responsive to diet in males; and approximately half as many diet-responsive DEGs in females: 133 (80 upregulated and 53 downregulated; P <0.01, Figure 3A, Additional File 4, Table S4). Specifically, the Trans-fat diet significantly downregulated the expression of many male cardiac genes involved in the regulation of biological processes, transcription and phosphorylation. Conversely females upregulated the expression of genes with functional ontologies relating to binding, nitrogen compound metabolism and ketone metabolism (Figure 3B-E, P < 0.01). Significant male DEGs downregulated by TFA consumption included key transcription modulator sterol regulatory element binding factor 2 (Srebf2, 3.5-fold, P < 0.01), transcription factors Mef2d (myocyte enhancer factor 2D, 2-fold) and Gata4 (2-fold); vasoconstrictive alpha-2b adrenergic receptor (2.8-fold); major cardiac Na+ channel Scn5a (2-fold), heat shock transcription factor 4 (Hsf4: 9.8-fold), calcium-regulating enzyme phospholipase C δ (1.7-fold), and circadian clock gene Per1 (period homolog 1: 2.6-fold). Signaling molecules calcium/calmodulin-dependent protein kinase kinase 1α (Camkk1) and mitogen-acrivated protein kinase kinase 4 (Map2k4) expression were also halved in male TFA mice. Conversely pro-hypertrophic Hdac2 (histone deacetylase 2), which represses Mef2 activity, was upregulated in males and downregulated in females. Other genes with this pattern of sex-dimorphic expression included plastin 3 (Pls3, 2.3-fold) and tachykinin 3 (Tacr3: 2.4-fold). In both sexes, the TFA diet led to increased expression of fatty acid oxidizing 3-hydroxy-3-methylglutaryl CoA synthase (Hmgcs2) and fructose bisphosphatase 2 (Fbp 2); however expression of TNF superfamily member LIGHT/Tnfsf14 was increased 7-fold in male Trans-fat hearts compared to control, whereas it was only increase 2-fold in female Trans-fat hearts. Diet-induced changes in female hearts were less significant, with only 84 DEGs significantly up-regulated, including progestin adipoQ receptor family member V (Paqr5: 2.6-fold) and RAS guanyl releasing protein 2 (Rasgrp2, 2.2-fold). Enrichment analysis of functional GO pathways within this subset revealed that females exhibited a four-fold less down-regulation of cardiac genes with classified biological function than males. We next investigated the relationship between these genes using Ingenuity Pathways Analysis (IPA). A network of 50 DEGs relating to cardiovascular disease and lipid metabolism among others was identified as being significantly affected by the TFA diet in cardiac tissue from male mice (Figure 4A), whereas a network of 25 functionally related genes involved in female lipid metabolism and cardiovascular development is shown in Figure 4B (P < 0.01).
Monosodium Glutamate consumption alters cardiac gene expression with sex-specific differences
A different pattern of gene expression was induced by MSG feeding, which again showed significant sex-specific differences in cardiac gene expression. 2-way ANOVA analysis revealed a subset of 153 male cardiac DEGs (98 upregulated and 55 downregulated) which were responsive to dietary MSG compared to control diet, and 174 female DEGs (117 upregulated and 57 downregulated). (Figure 1C and Additional File 5, Table S5, P < 0.01). Unlike Trans-fat, MSG upregulated more cardiac DEGs and down-regulated less. Genes/ESTs upregulated in males but not females in response to MSG included Ing4 (inhibitor of growth family member 4), glutathione S-transferase 4 and pancreatic lipase (Pnlip). Transcripts enriched for female-biased expression included myosin 7a (Myo7a), Pdzd3 (PDZ domain containing 3) and aldehyde dehydrogenase gene Aldh1a3. Functional analysis of Gene Ontologies of these subsets showed increased expression of male cellular and metabolic processes, whilst decreasing expression of phosphorylation-related DEGs (Additional File 6, Table S6, P < 0.01). Conversely females upregulated more genes/ESTs with developmental and anatomical function ontologies. Genes/ESTs that were upregulated in both sexes included ionotropic glutamate receptor Grik3 (2.5-fold in males and 1.7-fold in females); nebulin which regulates actin filament length (2.3-fold and 1.8-fold in males and females respectively).
MSG alters Trans-fat induced cardiac gene expression
In order to address the question as to whether the addition of MSG to a Trans-fat diet could affect cardiac gene expression, we next identified a subset of 338 DEGs which showed significant sex-dimorphism in response to these two diets. 197 DEGs (117 upregulated and 81 downregulated) in males showed significant response to the inclusion of MSG in the Trans fat diet, and 141 (60 upregulated, 81 downregulated) in females (Figure 1C and Additional File 7 Table S7, P < 0.01). Of note in both sexes, expression of lipogenic stearoyl-coA desaturase 4 (Scd4) was upregulated in TFA+MSG diet group hearts compared to TFA alone (1.6- and 1.4-fold respectively, P < 0.01) and was accompanied by an increase in Myosin H1, magnesium transporter Slc41a1, and apoptosis-inducing BCL2-interacting killer (Bik). Downregulated DEGs included fibrulin, sclerostin and dystrophin, together with fatty acid-oxidizing enzymes phospholipase A2 (Pla2g15) and aldehyde dehydrogenase (Aldh1a3). Significant differentially expressed cardiac transcripts in TFA + MSG-fed versus TFA-fed mice were examined for enrichment of GO and KEGG classifications (Additional File 8, Table S8). In males, major functional GO categories included metabolism, biosynthesis and transcription, whereas enrichment of female DEGs upregulated by MSG feeding again fell into the categories of developmental and anatomical structural function.
Cellular positional mapping reveals sex-specific differences in diet-deregulated gene expression
A final analysis of significant diet-induced deregulation of cardiac gene expression between all diet combinations co-analyzed yielded a cohort of 97 mapable genes which were deregulated by all dietary manipulations only in males (Figure 5 and Additional File 9, Table S9), and half that number of genes which were only deregulated in females (Figure 5 and Additional File 10, Table S10). Interestingly in this male subset, expression of genes previously demonstrated to be upregulated by TFA feeding (Figure 5A) were downregulated in animals consuming the MSG-enriched diets (Figure 5B and 5C, P < 0.01). Male-specific genes in this subset included stress-related proteins Dnaja1 and Pdia4, together with phosphodiesterase 6a and histone H1e . Females displayed a slightly divergent pattern (Figure 5D-F). Additionally 74 transcripts were identified as being significantly differentially expressed in both sexes in response to all dietary manipulations. Figure 6 shows a heatmap of a selection based on functional relevance, together with a number of X-linked and Y-linked diet-induced DEGs, P < 0.01 with a stringency of ± 1.5-fold change in expression. Eighteen genes were randomly chosen for further analyzed by QRT PCR based on biological relevance. Pearson correlation coefficients between the microarray analysis and QRT PCR were calculated. Ratios of expressions between the diet comparisons calculated from the microarray data set correlated well with the ratio calculated from the real-time PCR data (Figure 7, r = 0.718, P < 0.0001). A complete list of these genes and PCR primers is given in Additional File 11, Table S11.
Our results indicate that dietary manipulation prenatally, and over the first 8 months of life can markedly influence cardiac gene expression in a sex-specific manner. Sexual dimorphism in gene expression exists even in the healthy heart , and chromosome 4 has been shown to have more female-biased differential cardiac gene expression , whereas male-biased enrichment was found in chromosomes 3 and 14. In addition to the X and Y chromosomes which segregated as anticipated, we found female-biased chromosomal enrichment on mouse chromosomes 3, 4 and 12; and male-biased expression on chromosomes 8 and 11. Interestingly, a locus on murine chromosome 4 has been associated with markers of the Metabolic Syndrome, including body weight, insulin, leptin and triglyceride .
The main outcome from our study was that Trans-fat induced more than twice as many DEGs in males compared to females, and downregulated four times as many genes, including a number that have been implicated in cellular pathways leading to hypertrophy such as Gata4 and Mef2d. Expression of Gata4, previously identified as being of importance in stress responsiveness , was halved in males in response to the TFA diet. Gata4 regulates transcription of a large number of cardiac genes via binding to the consensus 5'-WGATAR-3' sequence, and is implicated in cardiac muscle apoptosis, where under conditions of oxidative stress Gata4 downregulation is preceded by MEK/ERK-dependent phosphorylation . Cardiac hypertrophy and CHD are mediated by transcription factor integration of upstream stress signals, and Mef2d-deficient mice have an attenuated hypertrophic response to cardiac stress . Mef2d activity is repressed by the histone deacetylase Hdac2, which we found to be increased in male TFA mice. Histone deacetylases are also regulated by calmodulin and calmodulin kinase, and are among the negative regulators of cardiac hypertrophy which attenuates Mef2 transcriptional activity . In the present study, we found a reduction of expression of Camkk1 and an increase in expression of Camk2b, suggesting that Trans-fat may deregulate both the calmodulin kinase/Mef2/Hdac pathway and the PKC/MAPK/GATA4 pathways of cardiac gene transcription.
Trans-fat feeding deregulated other stress-associated genes in males including phospholipase C δ1, the main cardiac isoform implicated in oxidative stress-induced REDOX signaling , and LIGHT/Tnfsf14, demonstrated to be upregulated in clinical and experimental heart failure . Additionally, key circadian clock gene Per1 was decreased by more than half in the hearts of male TFA-fed mice. Downregulation of Per1 by glucose has been shown to be accompanied by changes in the expression of genes involved in lipid metabolism, transcription and cell cycle , such as sterol regulatory element binding factor (Srebf2), hydroxy-3-methylglutaryl CoA synthase (Hmgcs) and stearoyl-coA desaturase (Scd), all of which were differentially regulated by Trans-fat in our study. Whilst it is known that Trans-fat consumption is linked to cardiovascular disease , the mechanism behind this is incompletely understood. Elaidic acid, the main TFA in partially hydrogenated vegetable shortening, incorporates into cardiac tissues ; and it has been suggested that TFAs are oxidized at a slower rate than cis-fatty acids , resulting in the accumulation of reactive oxidation intermediates , which could lead to oxidative stress and the activation of the Gata4, Mefd2 and Hdac-associated transcriptional pathways. Taken together, our results suggest that Trans-fat feeding results in significant impairment of transcriptional pathways relating to cellular stress, particularly in males. It has recently been shown that male C57BL/6J mice are more vulnerable than females to the impact of a high-fat diet in terms of ensuing weight gain and deleterious metabolic peturbations ; and it is therefore possible that the sex-differences in Trans fat -induced cardiac gene expression seen in our study may have contributed to the resultant increase in male weight gain and heart: body ratio compared to females.
The MSG diet resulted in a different pattern of sexual dimorphism in cardiac transcription, with fairly equal numbers of DEGs being up- and down-regulated. However there were a number of genes which were overexpressed in males but not in females, and vice versa. The functional categories of overexpressed genes differed too, with females showing more enrichment of genes associated with developmental processes for example, whereas male upregulated transcripts were enriched for cellular and metabolic function. Neonatal injections of MSG result in smaller hearts with less DNA synthesizing activity ; and hypotension has been observed in female, but not male rats treated neonatally with MSG . In our study, female hearts from MSG-diet group mice weighed significantly less than male hearts, however an examination of genes deregulated in our MSG mice did not provide many clues as to the mechanism responsible for this apparent difference.
We have previously shown that MSG alters Trans fat -induced hepatic and adipose tissue gene expression with the result that several key transcripts involved in ß-oxidation of lipids were downregulated in TFA+MSG treated mice ; and in our present study we found that the combination of TFA+MSG decreased the expression of several fatty acid-oxidizing enzymes such as phospholipase A2 and aldehyde dehydrogenase,whilst increasing the expression of lipogenic cardiac-specific stearoyl-coA desaturase 4 (Scd4), previously shown to be regulated by leptin and dietary factors .
When we compared cardiac gene transcription between all dietary combinations co-analyzed, we again found that more significant DEGs were induced in males compared to females. Our analysis identified a subset of over 100 transcripts which significantly changed only in male mice in response to diet, and approximately half this number of transcripts differentially expressed in females only, and a similar number of shared DEGs were common to both sexes. This suggests that males may have a higher cardiac transcriptional response to dietary manipulation, and indeed greater transcriptional response to exercise has previously been demonstrated in males but not females . Several theories have been put forward as to this disparity, the cardioprotective effect of estrogen being one possibility , especially in light of epidemiological data indicating a lower incidence of female CHD, at least until menopause, after which this dissimilarity declines. However, this is unlikely to account for all the sex-specific differences observed, since after correcting for lifestyle factors, the discrepancy in CHD rates was reduced .
By using microarrays to analyze diet-induced cardiac gene expression patterns, we have identified genes which are differentially regulated and may play important roles in the development of CHD associated with Trans-fat consumption. However there are some issues to consider when interpreting the results from the present study. For the microarray experiments, we employed a pooling design in order to reduce biological variation, and to lower the costs associated with sample size and chip availability . In order to address a possible loss of sensitivity and an increase in the incidence of false positives, we employed stringent criteria for the identification of the diet-regulated DEGs. One limitation of this study is the fact that these experiments were only performed on C57Bl/6J mice of both sexes. It would be advantageous but beyond the scope of the present study, to ascertain whether similar results could be obtained using a different genetic background, or possibly a higher animal model.
We have found substantial evidence for sexual dimorphism in response to dietary manipulation, with males apparently showing the greater cardiac transcriptional response. In particular, we found that Trans-fat feeding in male mice resulted in deregulation of the expression of genes involved in cellular stress pathways. Excessive TFA consumption is believed to be a risk factor in human heart disease, the incidence of which appears to be higher in male patients. Our present observations may provoke further studies which could result in sex-specific medical intervention in the future.
Animals and diets
Our study animals were bred from C57BL/6J mice obtained from The Jackson Laboratory (Maine, USA). Female breeders were housed/caged and fed a standard chow diet until 6 weeks of age whereupon they were placed on one of 4 different dietary regimens for a period of 3 weeks prior to mating as described previously . The four diet regimens used in this study were:  Standard Chow (Control diet) with ad lib drinking water.  Standard Chow, with ad lib drinking water containing 0.64 g/L Monosodium Glutamate (MSG diet). Trans-fat diet of 20% (w/w) Partially Hydrogenated Vegetable Shortening (Test Diet #5C4M containing 8.68% w/w Trans fatty acids; Test Diet Purina, USA), with ad lib drinking water (TFA diet). Trans-fat Diet #5C4M together with ad lib drinking water containing 0.64 g/L Monosodium Glutamate (TFA+MSG diet). See Table 2 for diet composition. Following mating, the 4 groups of dams were maintained on their respective diets throughout the gestation and nursing period. Male and female offspring (n = 20 per diet group) used in the following experiments were weighed, weaned onto the same diets and maintained on their respective dietary regimens until they reached 32 weeks of age. Average body weight was assessed at 6 and 32 weeks of age; water intake and food consumption was assessed at 6 and again at 28-30 weeks. The breeding and care of the animals were in accordance with the protocols approved by the Animal Care and Use Committee of the King Faisal Specialist Hospital & Research Centre.
Measurement of fasting serum cholesterols, triglyceride, glucose, insulin and HOMA-IR levels
Serum Triglyceride (TG), T-CHOL and HDL-C concentrations were measured in overnight fasted 32-week old mice using the Reflovet Plus instrument (Roche, F. Hoffmann-La Roche Ltd, Basel, Switzerland) according to the manufacturer's instructions. Overnight fasting blood glucose was measured using the Ascensia Contour glucometer (Bayer HealthCare, IN, USA). Fasting Serum insulin was measured using the ultrasensitive mouse insulin ELISA kit from Mercodia (Uppsala, Sweden), as described previously . Homeostatic Model Assessment Index (HOMA-IR) values, a measure of insulin resistance, were calculated according to the established formula: (fasting serum insulin μIU/ml) * (fasting serum glucose mM)/22.5 .
Animals were euthanized at 32 weeks of age by xylazine/ketamine intramuscular injection, and the hearts were rapidly removed and rinsed in saline solution. After weighing, the hearts were snap-frozen for RNA extraction. Total RNA was prepared from snap-frozen cardiac tissues from 32-week-old male and female mice (n = 8 per diet group) using Qiagen RNeasy Kit (Quiagen, CA, USA) according to the manufacturer's instructions and stored at -80°C. The integrity of total RNA was measured using a 2100 Bioanalyzer instrument and an RNA 6000 Nano LabChip assay (Agilent Technologies, CA, USA). RNA concentrations were determined by absorption at 260-nm wavelength with an ND-1000 spectrometer (Nanodrop Technologies, DE, USA).
Gene expression analysis
Gene expression in these 64 samples was analyzed using 16 GeneChip (R) Mouse Gene 1.0 ST arrays representing 28,853 genes. We used 2 chips per diet group, and applied pooled RNA from four mice per chip (see Figure 1 for design of the microarray experiments). Targets were prepared and microarrays were processed as described in the Affymetrix GeneChip Whole Transcript Expression Analysis manual using commercially available Affymetrix GeneChip WT cDNA Synthesis Kit, WT cDNA Amplification Kit, and WT Terminal Labeling Kit as per manufacturer's instructions. Briefly, approximately 200 ng of total RNA was used to synthesize double-stranded DNA with random hexamers tagged with a T7 promoter sequence. The cDNA was used as a template for in vitro transcription. In the second cycle cDNA synthesis, random primers were used in reverse transcription to convert the cRNA into single-stranded DNA, which was fragmented, labeled, and hybridized to the array for 16 hours using the Fluidics 450 station. Arrays were scanned using the Affymetrix 3000 7G scanner and GeneChip Operating Software version 1.4 to produce.CEL intensity files. This software also provided summary reports by which array QA metrics were evaluated including average background, average signal, and 3'/5' expression ratios for spike-in controls, β-actin, and GAPDH. Microarray data was deposited at the MIAME compliant NCBI gene expression hybridization array data repository (GEO: http://ncbi.nlm.nih.gov/geo) under accession #GSE22881.
Morphological and biochemical data were analyzed using one-way analysis of variance (ANOVA) followed by Tukey's test. Differences were considered significant at P < 0.05. For identification of differentially expressed genes (hereafter termed DEGs), microarray analysis was performed using the Partek genomic suite software version 6.3 (Partek, MO, USA) as previously described . Probe set data were summarized and background adjusted using the GC-Robust Multi-Array (GCRMA) algorithm . This method is used to reduce discrepancies in hybridization patterns that might result from variables in target amplification, hybridization conditions, staining or probe array lot. All data were normalized using non-linear transformation termed Quantile Normalization to improve the linearity, normality & homoscedasticity of the data, which was further filtered to remove noise and extreme expression values. A 2-way analysis of variance (ANOVA) was performed to detect significant differences among the diet groups and sexes. Probability values were adjusted for multiple testing using the False Discovery Rate (FDR) method .This method uses a controlled FDR while adjusting for testing simultaneously across multiple subgroup comparisons of sex and diet groups. Contrasts were included in the model based on the comparison of interest. All further sublists were created using genes that passed the FDR adjusted ANOVA p-value as well as fold change criteria. See Figure 1 for the design and statistical analysis of the microarray experiments. All resultant genes and expressed sequence tags (ESTs) meeting the criteria for differential expression were ascribed genome-wide significance using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) annotation . Testing for chromosomal enrichment of the DEGs was performed using the Fisher exact test. For detection of Gene Ontology (http://www.geneontology.org) categories and Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.kegg.com) pathways and pathway analysis, with significant overrepresentation of transcripts in diet groups compared to control. Heat maps of the gene/ESTs expression data were also generated using Partek. Networks of biologically related genes were created using Ingenuity Pathways Analysis (IPA: http://www.ingenuity.com). Sex-specific pathways identified from microarray analysis in the comparison of gene expression in Trans-fat fed mice compared to control were used to build literature-based male and female-specific pathways. In the two networks, genes or gene products are represented as colored nodes, and biological relationships between nodes are depicted as lines connecting the individual nodes. P ≤ 0.01 was considered significant, with a stringency of ± 1.5-fold difference in expression levels.
Confirmation of microarray results was performed by quantitative QPCR as previously described . Total RNA (2 μg) was reverse transcribed using Promega Reverse Transcription system (Promega, WI, USA). Subsequently, QPCR reactions were performed in triplicate using SYBR-Green 1 master mix (Applied Biosystems) and 10 ng cDNA as template. No template and no reverse transcriptase controls were included and products were analysed by gel electrophoresis. Real time RT-PCR values for each target gene were calculated as a ratio of target gene expression level to the 18-S ribosomal expression level in the same specimen. Statistical significance was assessed using a two-tailed t test assuming unequal variance of the biological replicates. Intron-spanning gene specific primers were designed using Primer 3 software and sequences are provided in Additional File 11. DEGs that differed significantly (P < 0.01) in their regulation between the diet groups' microarray analysis were selected, based on their biological relevance, and validated with the same samples by QRT- PCR analysis. Ratios of expressions between the diet comparisons were calculated from the microarray data set and Pearson correlation analysis between the QRT-PCR and microarray data were calculated.
Kyoto Encyclopedia of Genes and Genome
Trans Fatty Acid
Differentially Expressed Genes
Expressed Sequence Tags
Gabory A, Attig L, Junien C: Sexual dimorphism in environmental epigenetic programming. Mol Cell Endocrinol. 2009, 304 (1-2): 8-18. 10.1016/j.mce.2009.02.015.
Gómez-Pérez Y, Amengual-Cladera E, Català-Niell A, Thomàs-Moyà E, Gianotti M, Proenza AM, Lladó I: Gender dimorphism in high-fat-diet-induced insulin resistance in skeletal muscle of aged rats. Cell Physiol Biochem. 2008, 22 (5-6): 539-48. 10.1159/000185538.
Guerra B, Fuentes T, Delgado-Guerra S, Guadalupe-Grau A, Olmedillas H, Santana A, Ponce-Gonzalez JG, Dorado C, Calbet JA: Gender Dimorphism in Skeletal Muscle Leptin Receptors, Serum Leptin and Insulin Sensitivity. PLoS ONE. 2008, 3 (10): e3466-10.1371/journal.pone.0003466.
Zammaretti F, Panzica G, Eva C: Sex-dependent regulation of hypothalamic neuropeptide Y-Y1 receptor gene expression in moderate/high fat, high-energy diet-fed mice. J Physiol. 2007, 583 (Pt 2): 445-54.
Clarke SD: The multi-dimensional regulation of gene expression by fatty acids: polyunsaturated fats as nutrient sensors. Curr Opin Lipidol. 2004, 15: 13-18. 10.1097/00041433-200402000-00004.
Park TS, Yamashita H, Blaner WS, Goldberg IJ: Lipids in the heart: a source of fuel and a source of toxins. Curr Opin Lipidol. 2007, 18: 277-282. 10.1097/MOL.0b013e32814a57db.
Isensee J, Witt H, Pregla R, Hetzer R, Regitz-Zagrosek V, Noppinger PR: Sexually dimorphic gene expression in the heart of mice and men. J Mol Med. 2008, 86 (1): 61-74. 10.1007/s00109-007-0240-z.
Witt H, Schubert C, Jaekel J, Fliegner D, Penkalla A, Tiemann K, Stypmann J, Roepcke S, Brokat S, Mahmoodzadeh S, Brozova E, Davidson MM, Ruiz Noppinger P, Grohé C, Regitz-Zagrosek V: Sex-specific pathways in early cardiac response to pressure overload in mice. Mol Med. 2008, 86 (9): 1013-24. 10.1007/s00109-008-0385-4.
Fermin DR, Barac A, Lee S, Polster SP, Hannenhalli S, Bergemann TL, Grindle S, Dyke DB, Pagani F, Miller LW, Tan S, Dos Remedios C, Cappola TP, Margulies KB, Hall JL: Sex and age dimorphism of myocardial gene expression in nonischemic human heart failure. Circ Cardiovasc Genet. 2008, 1 (2): 117-25. 10.1161/CIRCGENETICS.108.802652.
Wells HF, Buzby JC: Dietary assessment of major trends in US food consumption (1970-2005). USDA Economic Information Bulletin. 2008, 33:
Kavanagh K, Jones KL, Sawyer J, Kelley K, Carr JJ, Wagner JD, Rudel LL: Trans fat diet induces abdominal obesity and changes in insulin sensitivity in monkeys. Obesity (Silver Spring). 2007, 15 (7): 1675-84. 10.1038/oby.2007.200.
Mozaffarian D, Pischon T, Hankinson SE, Rifai N, Joshipura K, Willett WC, Rimm EB: Dietary intake of trans fatty acids and systemic inflammation in women. Am J Clin Nutr. 2004, 79: 606-612.
Woodside JV, McKinley MC, Young IS: Saturated and trans fatty acids and coronary heart disease. Curr Atheroscler Rep. 2008, 10 (6): 460-466. 10.1007/s11883-008-0072-5.
Risérus U, Willett WC, Hu FB: Dietary fats and prevention of type 2 diabetes. Prog Lipid Res. 2009, 48 (1): 44-51. 10.1016/j.plipres.2008.10.002.
Mozaffarian D, Katan MB, Ascherio A, Stampfer MJ, Willett WC: Trans fatty acids and cardiovascular disease. N Engl J Med. 2006, 354: 1601-1613. 10.1056/NEJMra054035.
Ascherio A, Willett WC: Health effects of trans fatty acids. Am J Clin Nutr. 1997, 66 (4 Suppl): 1006S-1010S.
Mozaffarian D, Clarke R: Quantitative effects on cardiovascular risk factors and coronary heart disease risk of replacing partially hydrogenated vegetable oils with other fats and oils. Eur J Clin Nutr. 2009, 63 (Suppl 2): S22-33.
Barrett-Connor E: Sex differences in coronary heart disease Why are women so superior? The 1995 Ancel Keys Lecture. Circulation. 1997, 95 (1): 252-64.
Damanhori A, Al Khaja K, Sequeira RP: Gender-based treatment outcomes in diabetic hypertension. J Postgrad Med. 2008, 54 (4): 252-8. 10.4103/0022-3859.41433.
Silander K, Alanne M, Kristiansson K, Saarela O, Ripatti S, Auro K, Karvanen J, Kulathinal S, Niemelä M, Ellonen P, Vartiainen E, Jousilahti P, Saarela J, Kuulasmaa K, Evans A, Perola M, Salomaa V, Peltonen L: Gender differences in genetic risk profiles for cardiovascular disease. PLoS One. 2008, 3 (10): e3615-10.1371/journal.pone.0003615.
Maletínská L, Maixnerová J, Pirník Z, Kiss A, Zelezná B: Comparison of the obesity phenotypes related to monosodium glutamate effect on arcuate nucleus and/or the high fat diet feeding in C57Bl/6 and NMRI mice. Physiol Res. 2008, 57 (5): 727-34.
Hamaoka K, Sawada T: Hypoplastic heart induced by neonatal hypothalamic lesions in mice. J Mol Cell Cardiol. 1987, 19 (8): 741-9. 10.1016/S0022-2828(87)80385-1.
Pesini P, Rois JL, Menéndez L, Vidal S: The neonatal treatment of rats with monosodium glutamate induces morphological changes in the subfornical organ. Anat Histol Embryol. 2004, 33 (5): 273-7. 10.1111/j.1439-0264.2004.00547.x.
Larsen PJ, Mikkelsen JD, Jessop D, Lightman SL, Chowdrey HS: Neonatal monosodium glutamate treatment alters both the activity and the sensitivity of the rat hypothalamo-pituitary-adrenocortical axis. J Endocrinol. 1994, 141 (3): 497-503. 10.1677/joe.0.1410497.
Yu T, Zhao Y, Shi W, Ma R, Yu L: Effects of maternal oral administration of monosodium glutamate at a late stage of pregnancy on developing mouse fetal brain. Brain Res. 1997, 747: 195-206. 10.1016/S0006-8993(96)01181-X.
Beyreuther K, Biesalski HK, Fernstrom JD, Grimm P, Hammes WP, Heinemann U, Kempski O, Stehle P, Steinhart H, Walker R: Consensus meeting: monosodium glutamate - an update. Eur J Clin Nutr. 2007, 61: 304-13. 10.1038/sj.ejcn.1602526.
Rol de Lama MA, Pérez-Romero A, Hermanussen M, Ariznavarreta C, Tresguerres JA: Sexual dimorphism in growth as measured by microknemometry: different responses to GH deficiency and exogenous GH administration. Neuroendocrinology. 1998, 68 (3): 210-9. 10.1159/000054368.
Nakanishi Y, Tsuneyama K, Fujimoto M, Salunga TL, Nomoto K, An JL, Takano Y, Iizuka S, Nagata M, Suzuki W, Shimada T, Aburada M, Nakano M, Selmi C, Gershwin ME: Monosodium glutamate (MSG): a villain and promoter of liver inflammation and dysplasia. J Autoimmun. 2008, 30 (1-2): 42-50. 10.1016/j.jaut.2007.11.016.
Han HC, Austin KJ, Nathanielsz PW, Ford SP, Nijland MJ, Hansen TR: Maternal nutrient restriction alters gene expression in the ovine fetal heart. J Physiol. 2004, 558 (Pt 1): 111-21.
Collison KS, Maqbool Z, Saleh SM, Inglis A, Makhoul NJ, Bakheet R, Al-Johi M, Al-Rabiah R, Zaidi MZ, Al-Mohanna FA: Effect of dietary monosodium glutamate on trans fat-induced nonalcoholic fatty liver disease. J Lipid Res. 2009, 50 (8): 1521-37. 10.1194/jlr.M800418-JLR200.
Seidelmann SB, De Luca C, Leibel RL, Breslow JL, Tall AR, Welch CL: Quantitative trait locus mapping of genetic modifiers of metabolic syndrome and atherosclerosis in low-density lipoprotein receptor-deficient mice: identification of a locus for metabolic syndrome and increased atherosclerosis on chromosome 4. Arterioscler Thromb Vasc Biol. 2005, 25 (1): 204-10.
Ellmers LJ, Scott NJ, Piuhola J, Maeda N, Smithies O, Frampton CM, Richards AM, Cameron VA: Npr1-regulated gene pathways contributing to cardiac hypertrophy and fibrosis. J Mol Endocrinol. 2007, 38 (1-2): 245-57.
Suzuki YJ: Stress-induced activation of GATA-4 in cardiac muscle cells. Free Radic Biol Med. 2003, 34 (12): 1589-98. 10.1016/S0891-5849(03)00208-9.
Kim Y, Phan D, van RE, Wang DZ, McAnally J, Qi X, Richardson JA, Hill JA, Bassel-Duby R, Olson EN: The MEF2D transcription factor mediates stress-dependent cardiac remodeling in mice. J Clin Invest. 2008, 118: 124-132. 10.1172/JCI33255.
Park HJ, Begley U, Kong D, Yu H, Yin L, Hillgartner FB, Osborne TF, Galper JB: Role of sterol regulatory element binding proteins in the regulation of Galpha(i2) expression in cultured atrial cells. Circ Res. 2002, 91 (1): 32-7. 10.1161/01.RES.0000026502.79063.66.
Dahl CP, Gullestad L, Fevang B, Holm AM, Landrø L, Vinge LE, Fiane AE, Sandberg WJ, Otterdal K, Frøland SS, Damås JK, Halvorsen B, Aukrust P, Øie E, Yndestad A: Increased expression of LIGHT/TNFSF14 and its receptors in experimental and clinical heart failure. Eur J Heart Fail. 2008, 10 (4): 352-9. 10.1016/j.ejheart.2008.02.010.
Hirota T, Okano T, Kokame K, Shirotani-Ikejima H, Miyata T, Fukada Y: Glucose down-regulates Per1 and Per2 mRNA levels and induces circadian gene expression in cultured Rat-1 fibroblasts. J Biol Chem. 2002, 277 (46): 44244-51. 10.1074/jbc.M206233200.
Awad AB, Chattopadhyay JP: Alteration of rat heart sarcolemma lipid composition by dietary elaidic acid. J Nutr. 1983, 113 (4): 913-20.
Lawson LD, Kummerow FA: Beta-Oxidation of the coenzyme A esters of elaidic, oleic, and stearic acids and their full-cycle intermediates by rat heart mitochondria. Biochim Biophys Acta. 1979, 573 (2): 245-54.
Yu W, Liang X, Ensenauer RE, Vockley J, Sweetman L, Schulz H: Leaky beta-oxidation of a trans-fatty acid: incomplete beta-oxidation of elaidic acid is due to the accumulation of 5-trans-tetradecenoyl-CoA and its hydrolysis and conversion to 5-trans-tetradecenoylcarnitine in the matrix of rat mitochondria. J Biol Chem. 2004, 279 (50): 52160-7. 10.1074/jbc.M409640200.
Hwang LL, Wang CH, Li TL, Chang SD, Lin LC, Chen CP, Chen CT, Liang KC, Ho IK, Yang WS, Chiou LC: Sex differences in high-fat diet-induced obesity, metabolic alterations and learning, and synaptic plasticity deficits in mice. Obesity (Silver Spring). 2010, 18 (3): 463-9. 10.1038/oby.2009.273.
Clough RW, Aravich PF, Sladek CD: Monosodium glutamate neurotoxicity: a sex-specific impairment of blood pressure but not vasopressin in developing rats. Brain Res Bull. 1986, 17 (1): 51-8. 10.1016/0361-9230(86)90160-7.
Miyazaki M, Jacobson MJ, Man WC, Cohen P, Asilmaz E, Friedman JM, Ntambi JM: Identification and characterization of murine SCD4, a novel heart-specific stearoyl-CoA desaturase isoform regulated by leptin and dietary factors. Journal of Biological Chemistry. 2003, 278: 33904-33911. 10.1074/jbc.M304724200.
Paroo Z, Haist JV, Karmazyn M, Noble EG: Exercise improves postischemic cardiac function in males but not females: consequences of a novel sex-specific heat shock protein 70 response. Circ Res. 2002, 90 (8): 911-7. 10.1161/01.RES.0000016963.43856.B1.
Haynes SG, Feinleib M: Women, work and coronary heart disease: prospective findings from the Framingham heart study. Am J Public Health. 1980, 70: 133-141. 10.2105/AJPH.70.2.133.
Mendelsohn ME, Karas RH: The protective effects of estrogen on the cardiovascular system. N Engl J Med. 1999, 340: 1801-1811. 10.1056/NEJM199906103402306.
Shih JH, Michalowska AM, Dobbin K, Ye Y, Qiu TH, Green JE: Effects of pooling mRNA in microarray class comparisons. Bioinformatics. 2004, 20 (18): 3318-25. 10.1093/bioinformatics/bth391.
Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC: Homeostasis model assessment: insulin resistance and ß-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985, 28: 412-419. 10.1007/BF00280883.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf J, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 2: 249-264.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc Series B (Methodological). 1995, 57 (1): 289-300.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4 (5): P3-10.1186/gb-2003-4-5-p3.
We thank Rhea Mondreal and Rana Al-Rabiah for technical assistance. Thanks to the Comparative Medicine Department, KFSH&RC for animal husbandry. Our gratitude goes to Mr Hakim Al-Enazi for his unparalleled help in coordinating research resources. This project was funded in part by KFSH&RC project#2060 007. This research was supported by grant #08-MED490-20 from the National Comprehensive Plan for Science and Technology (NCPST), Kingdom of Saudi Arabia.
KSC devised the study, analyzed the data, wrote the manuscript and contributed to the figures and tables. MZZ analysed the data and contributed to the figures and tables; MS gave advice over statistical analysis; ZM performed the RNA extraction and Affymetrix hybridization; SMS performed the QRT-PCR; AI and NJM bred the experimental subjects, directed animal husbandry and performed all biochemical analysis, and FAA oversaw the project, contributed to the study design and critiqued the manuscript. All authors have read and approve this manuscript.
Electronic supplementary material
Additional File 1:Figure S1. Histology of cardiac tissue. A: transverse section and B: longitudinal sections of hematoxin and eosin-stained cardiac tissue from mice in control, MSG, TFA and TFA+MSG diet groups. Scale is x790. (PDF 180 KB)
Additional file 2:Table S2a. Differentially expressed cardiac genes in males compared to females (p ≤0.01) with a fold change of ≥ 1.5. Table S2b. Differentially expressed cardiac genes in females compared to males (p ≤0.01) with a fold change of ≥ 1.5. (PDF 2 MB)
Additional file 3:Table S3. Gene ontologies & pathways enriched for differentially expressed genes comparing all males to all females. (PDF 118 KB)
Additional file 4:Table S4. Differentially expressed genes in either males or females for the comparison TFA vs Control with respect to diet and sex (P < 0.01) and fold change ≥ ± 1.5. (PDF 363 KB)
Additional file 5:Table S5. Differentially expressed genes in either males or females for the comparison MSG vs Control with respect to diet and sex (P < 0.01) and fold change ≥ ± 1.5. (PDF 266 KB)
Additional file 6:Table S6. Gene ontologies enriched for differentially expressed genes comparing MSG to Control diet in males and females. (PDF 160 KB)
Additional file 7:Table S7. Differentially expressed genes in either males or females for the comparison TFA+MSG vs TFA with respect to diet and sex (P < 0.01) and fold Change ≥ ± 1.5. (PDF 280 KB)
Additional file 8:Table S8. Gene ontologies enriched for differentially expressed genes comparing TFA+MSG to TFA diet in males and females. (PDF 149 KB)
Additional file 9:Table S9. Intensities of genes/ESTs differentially regulated only in males with a fold change of ≥ ± 1.5 for either of the comparisons Control vs MSG, Control vs TFA, or TFA vs TFA+MSG. (PDF 188 KB)
Additional file 10:Table S10. Intensities of genes/ESTs differentially regulated only in females with a fold change of ≥ ± 1.5 for either of the comparisons Control vs MSG, Control vs TFA, or TFA vs TFA+MSG. (PDF 150 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Collison, K.S., Zaidi, M.Z., Maqbool, Z. et al. Sex-dimorphism in Cardiac Nutrigenomics: effect of Trans fat and/or Monosodium Glutamate consumption. BMC Genomics 12, 555 (2011). https://doi.org/10.1186/1471-2164-12-555
- Sexual Dimorphism
- Diet Group
- Trans Fatty Acid
- Monosodium Glutamate
- Cardiac Gene