- Research article
- Open Access
Gene expression profiles in rat mesenteric lymph nodes upon supplementation with Conjugated Linoleic Acid during gestation and suckling
BMC Genomics volume 12, Article number: 182 (2011)
Diet plays a role on the development of the immune system, and polyunsaturated fatty acids can modulate the expression of a variety of genes. Human milk contains conjugated linoleic acid (CLA), a fatty acid that seems to contribute to immune development. Indeed, recent studies carried out in our group in suckling animals have shown that the immune function is enhanced after feeding them with an 80:20 isomer mix composed of c9,t11 and t10,c12 CLA. However, little work has been done on the effects of CLA on gene expression, and even less regarding immune system development in early life.
The expression profile of mesenteric lymph nodes from animals supplemented with CLA during gestation and suckling through dam's milk (Group A) or by oral gavage (Group B), supplemented just during suckling (Group C) and control animals (Group D) was determined with the aid of the specific GeneChip® Rat Genome 230 2.0 (Affymettrix). Bioinformatics analyses were performed using the GeneSpring GX software package v10.0.2 and lead to the identification of 89 genes differentially expressed in all three dietary approaches. Generation of a biological association network evidenced several genes, such as connective tissue growth factor (Ctgf), tissue inhibitor of metalloproteinase 1 (Timp1), galanin (Gal), synaptotagmin 1 (Syt1), growth factor receptor bound protein 2 (Grb2), actin gamma 2 (Actg2) and smooth muscle alpha actin (Acta2), as highly interconnected nodes of the resulting network. Gene underexpression was confirmed by Real-Time RT-PCR.
Ctgf, Timp1, Gal and Syt1, among others, are genes modulated by CLA supplementation that may have a role on mucosal immune responses in early life.
Food components play a role in influencing, either directly or indirectly (through hormonal regulation), the expression of genes encoding for proteins involved in energy metabolism, cell differentiation and growth and immune responses. More specifically, diet exerts diverse effects on the development of the immune system, even at the level of gene regulation . It is known that polyunsaturated fatty acids (PUFAs) can modulate the expression of a variety of genes encoding for cytokines, adhesion molecules, and inflammatory proteins [2, 3]. This fact seems to be very important during early life, since docosahexanoic and arachidonic acids were reported to participate in the development of the neonate immune system, although their proportion among total fatty acids in human breast milk is very low . Human milk contains PUFAs such as conjugated linoleic acid (CLA), among others, that seem to contribute to immune development [5–8].
CLA is a class of positional and geometric conjugated dienoic isomers of linoleic acid, among which, cis9,trans11 (c9,t11) predominate, accounting for 83% to 100% of total CLA present in milk [5–7], whereas trans10,cis12 (t10,c12) exist in lower proportion [9, 10]. However, it is the intake of food of ruminant origin which determines the total concentration of CLA in dam's milk . Many other beneficial physiological effects have also been attributed to CLA, including reduced body fat, and inhibition of carcinogenesis, atherosclerosis, and diabetes [11–13].
Existing data regarding the effects of CLA on the immune system show great variability, mainly due to differences in the animal species used, the length of the supplementation period, and the differences in the isomer mixtures used in the experimental approach. In this direction, recent studies in suckling animals showed that the immune function is enhanced after feeding with an 80:20 isomer mix of c9,t11 and t10,c12 CLA [14, 15]. Specifically, sera IgG concentration and IgM in vitro production by splenocytes are increased after CLA supplementation during suckling. CLA downmodulatory effects on lymphoproliferation were only observed after an extra week of diet . The immune effects of CLA were also described in adult rats receiving this mixture since pregnancy . However, little work has been done on the effects of CLA on gene expression, and even less regarding the development of the immune system in early life.
Nutritional genomics is a result of the genetic revolution experienced over the past 10 years. Nutrigenomics deals with the interactions between dietary components and the genome and the resulting changes in proteins and other metabolites. On the other hand, nutrigenetics aims to understand the gene-based differences in response to dietary components and to develop nutraceuticals that are the most compatible with the health status of individuals based on their genetic makeup . The number of successful examples of transcriptome, proteome, and metabolome profiling as tools for evaluating the cellular responses to nutrients and identifying their molecular targets, has grown significantly. The use of high-density microarrays is a useful approach to estimate correlations among genes, which in turn can become the basis of transcriptional networks. The availability of microarrays for a number of model systems allows the quantification of relative transcript abundances in a comprehensive fashion. Despite the relatively simplistic nature of correlation measurements, they reflect an integrative view of gene-gene interactions in any given system, pointing out general structure characteristics of transcriptional interaction networks .
Additionally, nutrigenomic approaches have been undertaken to get further insight on the molecular understanding of mechanisms triggered by nutritional interventions. Diets enriched in different long-chain polyunsaturated fatty acids (LC-PUFAs) have been tested in rat nutritional intervention models. One report revealed steaoryl-CoA desaturase as an enzyme target for an arachidonate-enriched diet . In another study, Berger et al.  looked at transcriptional effects of these LC-PUFA-enriched diets on murine hepatic and hippocampal gene expression. Additionally, the beneficial effect of LC-PUFAs was assessed by a nutrigenomics experiment designed to understand the mechanisms by which these lipids induce and control gene signalling involved in carcinogenesis .
With all these concepts in mind and based on previous studies, we hypothesized that CLA intake during developmental phases would exert some influencing effect, among others, on genes involved in the regulation of the immune system. The aim of our study was the evaluation, by using whole genome microarrays, of the effects of dietary supplementation with an 80:20 isomer mix of c9,t11 and t10,c12 CLA, on mesenteric lymph nodes (MLN) gene expression, during gestation and/or suckling. The list of common genes differentially expressed in the three dietary interventions was used to construct a Biological Association Network (BAN). This approach allowed us to obtain a global view of gene expression in MLN, formed by a collection of nodes with different degrees of interrelationship, that could be used to explain the molecular mechanisms triggered by CLA.
Effect of CLA supplementation on rat MLN gene expression profiles
Wistar rats were subjected to a dietary supplementation during gestation and/or suckling with an 80:20 isomer mix of c9,t11 and t10,c12 CLA (Figure 1). CLA transfer was confirmed at day 21, when the proportion of CLA present in pup's plasma was around two times higher in groups A and B than that in group C, while values on reference group were very low . At day 21, MLN were obtained and total RNA was prepared. The expression profile of MLN from control animals and supplemented with CLA was determined with the aid of the specific GeneChip® Rat Genome 230 2.0, which includes more than 28,000 rat genes. Data from these microarrays were uploaded in the database repository of Gene Expression Omnibus (GEO, ) and can be accessed through series accession number GSE23004. Upon normalization and statistical filtering using GeneSpring GX software package v10.0.2, lists of differentially expressed genes by 2-fold were built as described in Methods, and presented as additional files 1, 2 and 3. CLA supplementation during gestation and suckling (groups A and B) modulated the expression of 1332 genes whereas the dietary supplementation only during suckling modulated the expression of 517 genes. Moreover, CLA supplementation during gestation and suckling through dam's milk (group A) decreased the expression of 1001 genes and up-regulated 103 genes in MLN. The dietary supplementation during gestation and suckling by oral gavage (group B) decreased the expression of 257 genes and up-regulated 371 genes, whereas supplementation only during suckling (group C) decreased the expression of 271 genes and up-regulated 351 genes in MLN.
Identification of common genes modulated by CLA supplementation
Venn diagrams were used to compare the lists of 2-fold differentially expressed genes obtained for the three experimental groups of animals (Figure 2). Lists of upregulated and dowregulated genes were compared separately in order to find the same expression patterns between the dietary interventions compared (e.g. genes upregulated in both). Experimental groups A and B shared 84 upregulated and 31 downregulated genes, groups A and C displayed 80 upregulated and 21 downregulated genes in common and 276 upregulated and 117 downregulated genes were found in common for both groups B and C. A total of 89 genes were found to be differentially expressed in all three conditions due to CLA supplementation (76 up plus 13 down). These 89 genes are presented as Table 1. We also analyzed global changes in gene expression independently of the type of regulation. We observed that 34 genes displayed downregulation in group A (CLA transfer through milk) and upregulation in groups B/C (CLA transfer by oral gavage) (Data not shown). This behavior suggests a role for the way of CLA administration in modulating gene expression.
Detection of gene nodes upon BAN generation
A BAN was generated as described in Methods with the list of differentially expressed genes in common among the three groups of animals. This type of graphical representation evidenced several genes, such as connective tissue growth factor (Ctgf), tissue inhibitor of metalloproteinase 1 (Timp1), galanin (Gal), synaptotagmin 1 (Syt1), growth factor receptor bound protein 2 (Grb2), actin gamma 2 (Actg2) and smooth muscle alpha actin (Acta2), as highly interconnected nodes of the resulting network (Figure 3).
Validation of node-genes selected from the BAN
We proceeded to validate the differential expression of the node-genes to verify the changes in their mRNA levels. Real-time PCR offers a non hybridization-based detection and was chosen as complementary to microarrays. The genes selected from the BAN representation were analyzed by RT-Real time PCR under the same experimental conditions as for the microarray analyses. Results are shown in Figure 4, where the levels of mRNA determined both in the microarrays and by RT-Real Time PCR for each gene in each experimental group are represented. In MLN, the mRNA levels for Gal and Actg2 were decreased by around 99%, and those of Timp1 and Syt1 decreased by 90-95% upon CLA administration in all groups of animals tested (p < 0.001). Treatment with CLA also decreased mRNA expression of Ctgf and Acta2 by around 80% (p < 0.01). Grb2 showed no changes in gene expression when determined by RT-Real Time PCR (p > 0.05, data not shown), at variance with the data obtained from the microarrays. These results confirmed the RNA data obtained in the screening performed using the microarrays for all the selected genes with the exception of Grb2.
In this work, we analyzed the gene expression profile of MLN cells from suckling rats receiving an 80:20 c9,t11:t10,c12 CLA mix during early life. This particular mixture enriched with the c9,t11 isomer was used due to its resemblance with that in human breast milk . Previous studies have demonstrated that this mix modulates immune development in suckling rats receiving it during gestation/suckling . The main effect of this type of diet was found regarding immunoglobulin production at systemic and intestinal level [14, 15]. However, the precise mechanism of action of this PUFA was not established. Thus, we used microarrays to get further insight in the effects of CLA supplementation on MLN gene expression. A large number of genes were modulated when CLA was administered during gestation and suckling (groups A and B, 5 wk). However, the use of Venn Diagrams evidenced that many genes were differentially expressed in both groups B and C, and less common genes were found between group A and the others. This behavior could be due to the different nutritional intervention that was performed. In the group A, CLA was added to the diet of the dams during suckling period and therefore CLA was incorporated in breast milk and transferred to the suckling pups. However, in groups B and C, CLA was directly administered by oral gavage to the suckling pups. This fact influences the proportion between the two CLA isomers that get to pups, as the ratio of c9,t11 CLA/t10,c12 CLA in pups's plasma was 83:14 for group A, but 93:7 for group B and 94:6 for group C as previously described .
The list of differentially expressed genes in the three conditions was used to generate a BAN. Several node genes, namely Ctgf, Timp1, Gal, Syt1, Actg2 and Acta2, were selected from this BAN, and changes in their mRNA expression were confirmed by RT-Real time PCR. Some of these node-genes may explain the behavior of CLA as immune modulator during early age.
Connective tissue growth factor (Ctgf) is a secreted extracellular matrix-associated protein that modulates many cellular functions, including proliferation, migration, adhesion, and extracellular matrix production, and it is involved in several biological and pathological processes . Dietary supplementation with CLA reduced mRNA levels of Ctgf in MLN in our three experimental approaches, indicating that downregulation of Ctgf is independent of the length and via of CLA supplementation. Ctgf is mainly regulated at the transcriptional level by bioactive lipids, as well as TGFβ and downstream mitogen-activated protein (MAP) kinase signaling .
We have previously reported that young (28-day-old) and adult animals fed with the same CLA mixture showed a reduction in the proliferation activity of spleen and MLN cells [16, 17]. Given that Ctgf promotes DNA synthesis in chondrocytes, osteoblasts and fibroblasts , the observed decrease in proliferation could be mediated through Ctgf downregulation.
Tissue inhibitors of metalloproteinases (Timp) are multifunctional proteins including four members (Timp1-4). Timp1 is expressed and secreted by a variety of cell types and is present in most tissues [27, 28]. Timp1 is able to promote cell growth, and its expression levels inversely correlate with the susceptibility to induction of apoptosis (reviewed in ).
Until now, two studies have reported different effects of CLA on Timp1 expression. On one hand, 1% CLA diet in a mice model of metastatic mammary tumor increased mRNA and protein levels of Timp1 , whereas an in vitro study using pure CLA isomers (c9,t11 CLA and t10,c12 CLA) on a macrophage cell line did not modify Timp1 gene expression . Thus, the effect of CLA on Timp1 expression may depend on the tissue microenvironment and the experimental approach used. In our study, CLA administration reduced mRNA levels for Timp1 in all three groups of animals. As Timp1 may stimulate cell proliferation by using several distinct signalling pathways, such as MAPK, the inhibitory effect of CLA on this gene may reduce lymphoproliferation. Thus, both Timp1 and Ctgf would play a role in a same signaling pathway that could explain the downmodulation of the lymphoproliferative response observed in 28-day-old animals fed with CLA from birth .
The neuropeptide Galanin (Gal) was also downregulated at the mRNA level by early CLA administration, independently of length, duration and via of supplementation. This peptide exerts a wide range of effects, not only in the central nervous system  but also in the enteric nervous system . Although Gal is mainly produced by neurons, it has been described that lymphocytes and macrophages can also synthesize it at a lower level . A clear effect of neuropeptides such as somatostatin (SOM) or vasoactive intestinal peptide (VIP), including gut-associated immune system (GALT), has been already reported on the immune system . When these neuropeptides are found in the intestine, they seem to affect the proliferation, differentiation and function of immune cells . Thus, the effect of CLA on Gal expression could explain the immunoenhancing effects of CLA on developing GALT .
There is no data available on the effects of Gal on lymphoproliferation. SOM and VIP are mainly inhibitors of lymphocyte proliferation, apparently by reducing interleukin-2 production , whereas calcitonin-gene-related-peptide (CGRP), substance P (SP) or met-enkephalin stimulate lymphoproliferation [35, 36]. Gal-mediated effects through its receptor Galr1 led to antiproliferative effects, whereas those through Galr2 and Galr3 were proliferative. All three receptors are expressed in immune cells . When the same CLA mix was administered to young rats supplemented during suckling  or until adult age , a downmodulatory activity on proliferation in both spleen and MLN was observed, suggesting that this effect may be mediated in part by modulation of Gal expression. Thus, CLA is able to reduce Gal promoting activity on lymphoproliferation, indicating that Gal would be acting mainly through its Galr2 and Galr3 receptors in that environment.
During early life CLA is able to enhance the production of immunoglobulins  and more specifically IgA in the intestine . Enhancement of intestinal IgA concentration by CLA could be due to the decrease of Gal expression in MLN, as it is the case for SOM in murine spleen, Peyer patches and MLN . Moreover, as neuropeptides influence cytokine production, downregulation of Gal by CLA may be the responsible effector of the changes in cytokine concentrations found in previous studies [14, 15]. However, the effect on Gal found here in early life, and their parallel immune effects, may not be present at older ages, since the maturation state of lymphocytes seem to affect the response to neuropeptides and cytokines due to their differential neuropeptidergic receptor expression . This fact is in agreement with the effect described for CLA on developing immune system and immunoglobulin production [15, 17], which is not observed later in life.
Synaptotagmins (Syts) constitute a family of proteins that act as Ca2+-sensors for regulated exocytosis and endocytosis [39–41]. Synaptotagmin proteins were first detected in the nervous system, but their expression has also been reported in other tissues. Thus, synaptotagmins have been proposed to be widespread regulatory proteins that confer Ca2+-sensitivity in a wide range of fusion events . Indeed, Syts have been described to play a role in the maintenance of plasma membrane integrity , in cell resealing , in chemotactic cell migration , in insulin release from the pancreatic β-cells [44, 45], in FGF-1 secretion in response to heat shock  and in degranulation of mast cells .
CLA supplementation caused downregulation of Syt1 by about 2.3 fold in the three dietary approaches. Syt1 is mainly localized in synaptic vesicles, where it functions as a Ca2+-sensor for exocytosis . Indeed, Syt1 has been proposed to play specific roles in synaptic vesicles trafficking, in facilitating docking , and in vesicle fusion . Although it has been traditionally proposed to be a neuron-specific isoform, Syt1 has been demonstrated to have important roles in membrane trafficking in non-neural and non-hormonal cells [51, 52].
Syt1 was previously associated with CLA supplementation, as its expression was significantly decreased in mice receiving a diet enriched in the t10,c12 CLA isomer compared with c9,t11 CLA . Syt1 downregulation was associated with the presence of hepatic fat, suggesting that changes in expression of Syt1 may indicate that other genes than the traditional are involved in lipid handling in hepatic cells, and probably in MLN. Syt1 has also been shown to play a role in neurotransmitter release [53, 54]. As described above, we found the neuropeptide Gal underexpressed upon CLA supplementation. Since it has been described that Gal is secreted by enteric nerves , we propose that downregulation of Syt1 upon CLA supplementation would lead to a decrease in Gal secretion, thus enhancing the negative effects that downregulation of Gal would have on lymphoproliferation.
Actins are highly conserved proteins that are involved in various types of cell motility, and maintenance of the cytoskeleton. In vertebrates, three main groups of actin isoforms, alpha, beta and gamma have been identified. The alpha actins are found in muscle tissues and are a major constituent of the contractile apparatus. The beta and gamma actins co-exist in most cell types as components of the cytoskeleton, and as mediators of internal cell motility. There are three α-actins (skeletal, cardiac, and smooth muscle), one β-actin (β-nonmuscle), and two γ-actins (γ-smooth muscle and γ-nonmuscle) . The cytoskeleton constitutes a network that serves to maintain cell shape and to regulate dynamic cellular functions. It provides essential scaffolding for the localization and activation or inhibition of diverse cytoplasmic signaling molecules, as well as anchors for motor proteins that are necessary for intracellular transport and cell division [57, 58]. In our microarray experiments we found that expression of smooth muscle γ-actin (Actg2), found in enteric tissues, and α2-actin (Acta2) were decreased upon dietary supplementation with CLA.
Alpha-smooth muscle actin expression is regulated by TGF-β and TNF-α. Changes in alpha-smooth muscle actin, collagen, and fibronectin expression result in decreased contraction and stiffness of collagen matrices . The effect of TGF-β has been related to a conserved TGF-β control element (TCE) within the 5'-region of alpha-smooth muscle actin, also known as smooth muscle cell differentiation marker. TCE could mediate both transcriptional activation and repression in cultured smooth muscle cells through interaction with members of the zinc finger Kruppel-like transcription factor family . On the other hand, it has been reported that TGF-β1 is able to significantly increase Acta2 expression in both normal and cancer associated fibroblasts cultures , which are believed to promote tumor growth and progression. Overall, the downregulation of actins upon CLA supplementation would affect both the contractile apparatus and the cytoskeleton with a final effect on cell motility. However its relationship with the development of the immune system remains unclear. Deficiencies in proteins that regulate the cortical actin cytoskeleton have increasingly been associated with immunodeficiency and autoimmune or autoinflammatory disease, indicating a critical role for these regulators in immune response and tolerance (reviewed in ).
Besides the modulatory action of CLA found on the above node genes, other genes with importance on the immune system were differentially modulated and would require further studies. This is the case of for the receptors for interleukin-7 (Il7r) and chemokines (i.e. Cxcr7) which are involved in mucosal immune responses and chemotaxis, respectively [62, 63].
In summary, by using rat whole genome microarrays, we determined the changes in gene expression induced by an 80:20 c9,t11:t10,c12 CLA mix in MLN from suckling rats. The generation of a biological association network allowed the identification of specific node genes that might be involved in immune responses. We conclude that Ctgf, Timp1, Gal and Syt1, among others, are gene target candidates modulated by CLA which may explain the effect of this PUFA on mucosal immune responses in early life.
Wistar rats at 7 days of gestation were obtained from Harlan (Barcelona, Spain) and housed in individual cages under conditions of controlled temperature and humidity in a 12 h:12 h light:dark cycle, with access to food and water ad libitum. The rats were monitored daily and allowed to deliver at term. The day of birth was reported as day 1. Pups, unified in litters of ten per lactating dam, had free access to the nipples and rat diet. Animals were daily identified and weighed, and handling was done in the same time range to avoid the influence of biological rhythms. On day 21 (weaning day), pups were anaesthetized with ketamine/xylazine (ketamine 90 mg/kg plus xylacine 10 mg/kg of rat weight) and sacrificed by humanitarian methods. MLN were obtained and were flash-frozen in liquid N2 and stored immediately at -80°C until processing. Studies were performed in accordance with the institutional guidelines for the care and use of laboratory animals established by the Ethical Committee for Animal Experimentation of the University of Barcelona and approved by the Catalonian Government (CEEA 303/05, UB/DMA 3242).
Animals were randomly distributed in 4 dietary groups, according to total period of CLA supplementation and administration route used in the pups . Total period of CLA supplementation (TPS) is shown in parentheses in the experimental design (Figure 1). Pups from dams fed with a 1% CLA diet during the two last weeks of gestation and throughout the suckling period constituted Group A. These pups received CLA through the dam's milk during suckling (TPS 5 wk). Pups from dams fed only during gestation (two last weeks) with a 1% CLA diet and receiving a daily CLA supplement by oral gavage throughout the suckling period (TPS 5 wk) constituted group B. Pups from dams fed with a standard diet during the two last weeks of gestation and suckling and receiving CLA by daily oral gavage throughout the suckling period (TPS 3 wk) constituted Group C. Group D, pups from dams fed standard diet throughout the study, constituted the reference diet group (TPS 0 wk).
In this later group, CLA levels present in milk samples were very low for the c9,t11-CLA isomer (0.02 ± 0.00% in total milk fatty acids) and undetectable for the t10,c12-CLA isomer .
Dietary CLA supplementation
The standard diet corresponded to the American Institute of Nutrition (AIN)-93G formulation , containing 7% soybean oil. A 1% CLA diet was obtained from modified standard flour AIN-513 (Harlan) containing 10 g CLA/kg as previously described . Thus, the supplemented diet contained 6% soybean oil plus 1% CLA oil. The CLA isomer mixture used was 80% c9,t11 and 20% t10,c12 from the total CLA isomers in oil. The 1% CLA diet in suckling animals corresponded to a daily administration of 1.5 mg CLA oil provided/g rat from day 1 to 21. Low-capacity syringes (Hamilton Bonaduz AG) adapted to oral 25- or 23-gauge gavage tubes, 27 mm in length (ASICO) were used for oral administration before and after day 5, respectively. To allow gastric emptying, litters were separated from dams 1 h before oral supplementation . CLA arrival to pups was confirmed by its quantification in 21-day-old plasma; moreover, immunomodulatory action of CLA in these animals by means of evaluation of sera Ig concentration, and in vitro Ig production and proliferation by isolated splenocytes was also previously observed .
Three animals were randomly selected from the 10 that constituted each experimental group. Total RNA from MLN of each animal was prepared using the RNAeasy Lipid Tissue Mini kit (Qiagen) following the recommendations of the manufacturer. Briefly, tissue samples were thawed and homogenized in QIAzol Lysis Reagent. After addition of chloroform, the homogenate was separated into aqueous and organic phases by centrifugation. The aqueous phase was applied to the RNeasy spin column and RNA was eluted in RNase-free water. An additional step of phenol-chloroform extraction and ethanol precipitation was performed to ensure the purity of the RNA samples from MLN. Ribosomal RNA band integrity was assessed on an Agilent BioAnalyzer 2100 using an RNA Nano LabChip (Agilent Technologies).
Gene expression was determined by hybridization to the GeneChip® Rat Genome 230 2.0 (Affymettrix), that allows the simultaneous analysis of the expression level of over 30,000 transcripts and variants from over 28,000 well-substantiated rat genes. Labeling, hybridization and detection were carried out following the manufacturer's specifications. Triplicate samples were hybridized for each experimental condition.
Microarray data analysis
Gene expression analysis was carried out with GeneSpring GX v10.0.2 software (Agilent Technologies), using the latest gene annotations available. This software package allows multifilter comparisons using data from different experiments to perform the normalization, generation of restriction lists and the functional classification of the differentially expressed genes. All the samples were normalized against the median of the control samples (Group D). The expression of each gene is reported as the ratio of the value obtained after each condition relative to control condition after normalization of the data. Lists of differentially expressed genes were generated using data from the three independent experiments for each condition. A first filter was applied to select the genes that displayed a p-value of less than 0.05. The output of this analysis was then filtered by fold expression, generating lists of differentially expressed genes by 2-fold for each of the experimental groups. These lists were split in two others of upregulated or downregulated genes. Comparisons of the lists of upregulated genes among them were performed by Venn Diagrams in GeneSpring GX. Lists of downregulated genes were also compared among them using the same approach. This procedure allowed us to find differentially expressed genes that followed the same pattern (e.g. upregulated or downregulated) in common among the experimental conditions.
Biological Association Networks generation
The list of common genes differentially expressed by 2-fold with a p-value < 0.05 in the tree groups of animals was used to construct a biological association network (BAN) using the Pathway Architect software integrated within GeneSpring GX. Briefly, the Pathway Architect software generates interaction networks starting with the genes in a given list (entities) taking into account the information present in a database of known molecular interactions. The lists correspond to the collection of differentially expressed genes under specific conditions. The database of molecular interactions is composed by more than 1.6 million interactions divided into different classes (binding, regulation, promoter binding, transport, metabolism, protein metabolism and expression). The interactions are extracted from literature using a Natural Language Processing tool run on Medline Abstracts (NLP references), plus those obtained from external curated databases like BIND  and MINT . Interactions in the interaction database are scored into 5 different categories: maximum, high, medium, low and minimal. Curated interactions (BIND and MINT sources) get the Maximum quality score as do any interactions which have at least 3 NLP references. Pathway Architect gathers all that information to construct novel views as to how the entities in a list could be interacting with each other, even including entities not present in the original list (neighbors resulting from the expanded interaction). Customized analyses were performed to select relevance interaction networks with an associated high confidence index since such networks are likely to mirror biological significance. One-step expansion (expand network) of the original set of entities with maximum score interaction were then analyzed by setting an advanced filter that included the categories of binding, expression, metabolism, promoter binding, protein modification and regulation. This procedure gives a final view formed by a collection of nodes with different degrees of interrelationship. A number of gene products from the original list that were not significantly connected with the other members or neighbors were removed from the final view.
RT-Real Time PCR
cDNA was synthesized in a total volume of 20 μl by mixing 1 μg of total RNA, 125 ng of random hexamers (BioTools), in the presence of 75 mM KCl, 3 mM MgCl2, 10 mM dithiothreitol, 20 units of RNAsin (Promega), 0.5 mM dNTPs (BioTools), 200 units of M-MLV reverse transcriptase (Invitrogen) and 50 mM Tris-HCl buffer, pH 8.3. The reaction mixture was incubated at 37°C for 60 min. The cDNA product was used for subsequent amplification by Real Time-PCR. The expression levels of outlier genes differentially expressed in the microarrays was determined in an ABI Prism 7000 Sequence Detection System (Applied Biosystems) using 3 μl of the cDNA mixture and the Assays-on-demand Rn00583681_m1 for Gal, Rn00587558_m1 for Timp1, Rn00573960_g1 for Ctgf, Rn00578230_m1 for Grb2, Rn00436862_m1 for Syt1, Rn00563662_m1 for Actg2, Rn01759928_g1 for Acta2 and Rn01432775_m1 for Aprt (all from Applied Biosystems). Aprt mRNA was used as an endogenous control. The reaction was performed following the manufacturers recommendations. Fold-changes in gene expression were calculated using the standard ΔΔCt method.
Sanderson IR, Naik S: Dietary regulation of intestinal gene expression. Annu Rev Nutr. 2000, 20: 311-338. 10.1146/annurev.nutr.20.1.311.
Wallace FA, Miles EA, Evans C, Stock TE, Yaqoob P, Calder PC: Dietary fatty acids influence the production of Th1- but not Th2-type cytokines. J Leukoc Biol. 2001, 69: 449-457.
Calder PC, Field CJ: Fatty acids, inflamation and immunity. Nutrition and Immune Function. Edited by: Calder PC, Field CJ, Gill HS. 2002, Southhampton (UK): CAB International, 57-92. full_text.
Field CJ, Clandinin MT, Van Aerde JE: Polyunsaturated fatty acids and T-cell function: implications for the neonate. Lipids. 2001, 36: 1025-1032. 10.1007/s11745-001-0813-6.
Mcguire MK, Park Y, Behre RA, Harrison LY, Shultz TD, Mcguire MA: Conjugated linoleic acid concentrations of human milk and infant formula. Nutr Res. 1997, 17: 1277-1283. 10.1016/S0271-5317(97)00111-5.
Jensen RG, Lammi-Keefe C: The anticarcinogenic conjugated fatty acid c9, t11-c18:2, or rumenic acid, in human milk: amounts and effects. Adv Exp Med Biol. 2001, 501: 153-156.
Luna P, Juarez M, de la Fuente MA: Fatty acid and conjugated linoleic acid isomer profiles in human milk fat. Eur J Lipid Sci Technol. 2007, 109: 1160-1166. 10.1002/ejlt.200700139.
Chin SF, Storkson JM, Albright KJ, Cook ME, Pariza MW: Conjugated linoleic acid is a growth factor for rats as shown by enhanced weight gain and improved feed efficiency. J Nutr. 1994, 124: 2344-2349.
Pariza MW, Park Y, Cook ME: Mechanisms of action of conjugated linoleic acid: evidence and speculation. Proc Soc Exp Biol Med. 2000, 223: 8-13. 10.1046/j.1525-1373.2000.22302.x.
Hayashi AA, de Medeiros SR, Carvalho MH, Lanna DP: Conjugated linoleic acid (CLA) effects on pups growth, milk composition and lipogenic enzymes in lactating rats. J Dairy Res. 2007, 74: 160-166. 10.1017/S002202990600224X.
Visonneau S, Cesano A, Tepper SA, Scimeca JA, Santoli D, Kritchevsky D: Conjugated linoleic acid suppresses the growth of human breast adenocarcinoma cells in SCID mice. Anticancer Res. 1997, 17: 969-973.
Ip MM, Masso-Welch PA, Shoemaker SF, Shea-Eaton WK, Ip C: Conjugated linoleic acid inhibits proliferation and induces apoptosis of normal rat mammary epithelial cells in primary culture. Exp Cell Res. 1999, 250: 22-34. 10.1006/excr.1999.4499.
Kelley NS, Hubbard NE, Erickson KL: Conjugated linoleic acid isomers and cancer. J Nutr. 2007, 137: 2599-2607.
Perez-Cano FJ, Ramirez-Santana C, Molero-Luis M, Castell M, Rivero M, Castellote C, Franch A: Mucosal IgA increase in rats by continuous CLA feeding during suckling and early infancy. J Lipid Res. 2009, 50: 467-476. 10.1194/jlr.M800356-JLR200.
Ramirez-Santana C, Perez-Cano FJ, Castellote C, Castell M, Rivero M, Rodriguez-Palmero M, Franch A: Higher immunoglobulin production in conjugated linoleic acid-supplemented rats during gestation and suckling. Br J Nutr. 2009, 102: 858-868. 10.1017/S0007114509301592.
Ramírez-Santana C, Castellote C, Castell M, Moltó-Puigmartí C, Rivero M, Pérez-Cano FJ, Franch A: Enhancement of antibody synthesis in rats by feeding cis-9,trans-11 conjugated linoleic acid during early life. J Nutr Biochem. 2010 , 22: 495-501. 10.1016/j.jnutbio.2010.04.005 .
Ramirez-Santana C, Castellote C, Castell M, Rivero M, Rodriguez-Palmero M, Franch A, Perez-Cano FJ: Long-term feeding of the cis-9,trans-11 isomer of conjugated linoleic acid reinforces the specific immune response in rats. J Nutr. 2009, 139: 76-81.
Kaput J: Nutrigenomics research for personalized nutrition and medicine. Curr Opin Biotechnol. 2008, 19: 110-120. 10.1016/j.copbio.2008.02.005.
Hubner N, Wallace CA, Zimdahl H, Petretto E, Schulz H, Maciver F, Mueller M, Hummel O, Monti J, Zidek V, et al: Integrated transcriptional profiling and linkage analysis for identification of genes underlying disease. Nat Genet. 2005, 37: 243-253. 10.1038/ng1522.
Mutch DM, Grigorov M, Berger A, Fay LB, Roberts MA, Watkins SM, Williamson G, German JB: An integrative metabolism approach identifies stearoyl-CoA desaturase as a target for an arachidonate-enriched diet. FASEB J. 2005, 19: 599-601. 10.1096/fj.05-3911rev.
Berger A, Mutch DM, German JB, Roberts MA: Dietary effects of arachidonate-rich fungal oil and fish oil on murine hepatic and hippocampal gene expression. Lipids Health Dis. 2002, 1: 2-10.1186/1476-511X-1-2.
Anderle P, Farmer P, Berger A, Roberts MA: Nutrigenomic approach to understanding the mechanisms by which dietary long-chain fatty acids induce gene signals and control mechanisms involved in carcinogenesis. Nutrition. 2004, 20: 103-108. 10.1016/j.nut.2003.09.018.
Gene Expression Omnibus. [http://www.ncbi.nlm.nih.gov/geo]
Rachfal AW, Brigstock DR: Structural and functional properties of CCN proteins. Vitam Horm. 2005, 70: 69-103. 10.1016/S0083-6729(05)70003-0.
Cicha I, Goppelt-Struebe M: Connective tissue growth factor: context-dependent functions and mechanisms of regulation. Biofactors. 2009, 35: 200-208. 10.1002/biof.30.
Chen CC, Lau LF: Functions and mechanisms of action of CCN matricellular proteins. Int J Biochem Cell Biol. 2009, 41: 771-783. 10.1016/j.biocel.2008.07.025.
Welgus HG, Stricklin GP: Human skin fibroblast collagenase inhibitor. Comparative studies in human connective tissues, serum, and amniotic fluid. J Biol Chem. 1983, 258: 12259-12264.
Nuttall RK, Sampieri CL, Pennington CJ, Gill SE, Schultz GA, Edwards DR: Expression analysis of the entire MMP and TIMP gene families during mouse tissue development. FEBS Lett. 2004, 563: 129-134. 10.1016/S0014-5793(04)00281-9.
Stetler-Stevenson WG: Tissue inhibitors of metalloproteinases in cell signaling: metalloproteinase-independent biological activities. Sci Signal. 2008, 1: re6-10.1126/scisignal.127re6.
Hubbard NE, Lim D, Erickson KL: Conjugated linoleic acid alters matrix metalloproteinases of metastatic mouse mammary tumor cells. J Nutr. 2007, 137: 1423-1429.
Ringseis R, Schulz N, Saal D, Eder K: Troglitazone but not conjugated linoleic acid reduces gene expression and activity of matrix-metalloproteinases-2 and -9 in PMA-differentiated THP-1 macrophages. J Nutr Biochem. 2008, 19: 594-603. 10.1016/j.jnutbio.2007.08.003.
McCown TJ: Adeno-associated virus vector-mediated expression and constitutive secretion of galanin suppresses limbic seizure activity. Neurotherapeutics. 2009, 6: 307-311. 10.1016/j.nurt.2009.01.004.
Genton L, Kudsk KA: Interactions between the enteric nervous system and the immune system: role of neuropeptides and nutrition. Am J Surg. 2003, 186: 253-258. 10.1016/S0002-9610(03)00210-1.
Pascual DW, Stanisz AM, Bienenstock J, Bost KL: Neural intervention in mucosal immunity. Mucosal Immunology. Edited by: Ogra PL, Mestecky J, Lamm ME, Strober W, Bienenstock J, McGhee JR. 1999, San Diego: Academic Press, 631-642. 2
Gabrilovac J, Marotti T: Gender-related differences in murine T- and B-lymphocyte proliferative ability in response to in vivo [Met(5)]enkephalin administration. Eur J Pharmacol. 2000, 392: 101-108. 10.1016/S0014-2999(00)00118-7.
Del Rio M, Hernanz A, de la Fuente M: Bombesin, gastrin-releasing peptide, and neuromedin C modulate murine lymphocyte proliferation through adherent accessory cells and activate protein kinase C. Peptides. 1994, 15: 15-22. 10.1016/0196-9781(94)90164-3.
Sugimoto T, Seki N, Shimizu S, Kikkawa N, Tsukada J, Shimada H, Sasaki K, Hanazawa T, Okamoto Y, Hata A: The galanin signaling cascade is a candidate pathway regulating oncogenesis in human squamous cell carcinoma. Genes Chromosomes Cancer. 2009, 48: 132-142. 10.1002/gcc.20626.
Stanisz AM, Befus D, Bienenstock J: Differential effects of vasoactive intestinal peptide, substance P, and somatostatin on immunoglobulin synthesis and proliferations by lymphocytes from Peyer's patches, mesenteric lymph nodes, and spleen. J Immunol. 1986, 136: 152-156.
Tucker WC, Chapman ER: Role of synaptotagmin in Ca2+-triggered exocytosis. Biochem J. 2002, 366: 1-13.
Chakrabarti S, Kobayashi KS, Flavell RA, Marks CB, Miyake K, Liston DR, Fowler KT, Gorelick FS, Andrews NW: Impaired membrane resealing and autoimmune myositis in synaptotagmin VII-deficient mice. J Cell Biol. 2003, 162: 543-549. 10.1083/jcb.200305131.
Brose N, Petrenko AG, Sudhof TC, Jahn R: Synaptotagmin: a calcium sensor on the synaptic vesicle surface. Science. 1992, 256: 1021-1025. 10.1126/science.1589771.
McNeil PL, Steinhardt RA: Loss, restoration, and maintenance of plasma membrane integrity. J Cell Biol. 1997, 137: 1-4. 10.1083/jcb.137.1.1.
Colvin RA, Means TK, Diefenbach TJ, Moita LF, Friday RP, Sever S, Campanella GS, Abrazinski T, Manice LA, Moita C: Synaptotagmin-mediated vesicle fusion regulates cell migration. Nat Immunol. 11: 495-502. 10.1038/ni.1878.
Lang J, Fukuda M, Zhang H, Mikoshiba K, Wollheim CB: The first C2 domain of synaptotagmin is required for exocytosis of insulin from pancreatic beta-cells: action of synaptotagmin at low micromolar calcium. EMBO J. 1997, 16: 5837-5846. 10.1093/emboj/16.19.5837.
Gao Z, Reavey-Cantwell J, Young RA, Jegier P, Wolf BA: Synaptotagmin III/VII isoforms mediate Ca2+-induced insulin secretion in pancreatic islet beta -cells. J Biol Chem. 2000, 275: 36079-36085. 10.1074/jbc.M004284200.
Mouta Carreira C, LaVallee TM, Tarantini F, Jackson A, Lathrop JT, Hampton B, Burgess WH, Maciag T: S100A13 is involved in the regulation of fibroblast growth factor-1 and p40 synaptotagmin-1 release in vitro. J Biol Chem. 1998, 273: 22224-22231. 10.1074/jbc.273.35.22224.
Textor B, Licht AH, Tuckermann JP, Jessberger R, Razin E, Angel P, Schorpp-Kistner M, Hartenstein B: JunB is required for IgE-mediated degranulation and cytokine release of mast cells. J Immunol. 2007, 179: 6873-6880.
Sudhof TC: Synaptotagmins: why so many?. J Biol Chem. 2002, 277: 7629-7632. 10.1074/jbc.R100052200.
Reist NE, Buchanan J, Li J, DiAntonio A, Buxton EM, Schwarz TL: Morphologically docked synaptic vesicles are reduced in synaptotagmin mutants of Drosophila. J Neurosci. 1998, 18: 7662-7673.
Geppert M, Goda Y, Hammer RE, Li C, Rosahl TW, Stevens CF, Sudhof TC: Synaptotagmin I: a major Ca2+ sensor for transmitter release at a central synapse. Cell. 1994, 79: 717-727. 10.1016/0092-8674(94)90556-8.
Rastaldi MP, Armelloni S, Berra S, Calvaresi N, Corbelli A, Giardino LA, Li M, Wang GQ, Fornasieri A, Villa A, Heikkila E, Soliymani R, Boucherot A, Cohen CD, Kretzler M, Nitsche A, Ripamonti M, Malgaroli A, Pesaresi M, Forloni GL, Schlöndorff D, Holthofer H, D'Amico G: Glomerular podocytes contain neuron-like functional synaptic vesicles. FASEB J. 2006, 20: 976-978. 10.1096/fj.05-4962fje.
Guillen N, Navarro MA, Arnal C, Noone E, Arbones-Mainar JM, Acin S, Surra JC, Muniesa P, Roche HM, Osada J: Microarray analysis of hepatic gene expression identifies new genes involved in steatotic liver. Physiol Genomics. 2009, 37: 187-198. 10.1152/physiolgenomics.90339.2008.
Koh TW, Bellen HJ: Synaptotagmin I, a Ca2+ sensor for neurotransmitter release. Trends Neurosci. 2003, 26: 413-422. 10.1016/S0166-2236(03)00195-4.
Littleton JT, Bai J, Vyas B, Desai R, Baltus AE, Garment MB, Carlson SD, Ganetzky B, Chapman ER: synaptotagmin mutants reveal essential functions for the C2B domain in Ca2+-triggered fusion and recycling of synaptic vesicles in vivo. J Neurosci. 2001, 21: 1421-1433.
Hoyle CH, Burnstock G: Galanin-like immunoreactivity in enteric neurons of the human colon. J Anat. 1989, 166: 23-33.
Wickramarachchi DC, Theofilopoulos AN, Kono DH: Immune pathology associated with altered actin cytoskeleton regulation. Autoimmunity. 43: 64-75. 10.3109/08916930903374634.
Janmey PA: The cytoskeleton and cell signaling: component localization and mechanical coupling. Physiol Rev. 1998, 78: 763-781.
Chesarone MA, Goode BL: Actin nucleation and elongation factors: mechanisms and interplay. Curr Opin Cell Biol. 2009, 21: 28-37. 10.1016/j.ceb.2008.12.001.
Goldberg MT, Han YP, Yan C, Shaw MC, Garner WL: TNF-alpha suppresses alpha-smooth muscle actin expression in human dermal fibroblasts: an implication for abnormal wound healing. J Invest Dermatol. 2007, 127: 2645-2655. 10.1038/sj.jid.5700890.
Liu Y, Sinha S, Owens G: A transforming growth factor-beta control element required for SM alpha-actin expression in vivo also partially mediates GKLF-dependent transcriptional repression. J Biol Chem. 2003, 278: 48004-48011. 10.1074/jbc.M301902200.
Casey TM, Eneman J, Crocker A, White J, Tessitore J, Stanley M, Harlow S, Bunn JY, Weaver D, Muss H, Plaut K: Cancer associated fibroblasts stimulated by transforming growth factor beta1 (TGF-beta 1) increase invasion rate of tumor cells: a population study. Breast Cancer Res Treat. 2008, 110: 39-49. 10.1007/s10549-007-9684-7.
Fagarasan S, Kawamoto S, Kanagawa O, Suzuki K: Adaptive immune regulation in the gut: T cell-dependent and T cell-independent IgA synthesis. Annu Rev Immunol. 2010, 28: 243-273. 10.1146/annurev-immunol-030409-101314.
Neusser MA, Kraus AK, Regele H, Cohen CD, Fehr T, Kerjaschki D, Wuthrich RP, Penfold ME, Schall T, Segerer S: The chemokine receptor CXCR7 is expressed on lymphatic endothelial cells during renal allograft rejection. Kidney Int. 2010, 77: 801-808. 10.1038/ki.2010.6.
Reeves PG, Nielsen FH, Fahey GC: AIN-93 purified diets for laboratory rodents: final report of the American Institute of Nutrition ad hoc writing committee on the reformulation of the AIN-76A rodent diet. J Nutr. 1993, 123: 1939-1951.
Perez-Cano FJ, Marin-Gallen S, Castell M, Rodriguez-Palmero M, Rivero M, Franch A, Castellote C: Bovine whey protein concentrate supplementation modulates maturation of immune system in suckling rats. Br J Nutr. 2007, 98 (Suppl 1): S80-84.
Bader GD, Donaldson I, Wolting C, Ouellette BF, Pawson T, Hogue CW: BIND--The Biomolecular Interaction Network Database. Nucleic Acids Res. 2001, 29: 242-245. 10.1093/nar/29.1.242.
Zanzoni A, Montecchi-Palazzi L, Quondam M, Ausiello G, Helmer-Citterich M, Cesareni G: MINT: a Molecular INTeraction database. FEBS Lett. 2002, 513: 135-140. 10.1016/S0014-5793(01)03293-8.
This work was supported in part by SAF08-43 and Generalitat de Catalunya (SGR05-00833 and SGR09-118). The oil used in the study was a gift from Loders Croklaan, Lipid Nutrition, Wormerveer, The Netherlands.
The authors declare that they have no competing interests.
CC, FJPC, MR, AF and VN designed the study and supervised the experimental work. ES, CRS, FJPC, CC, AF and VN performed the experimental work. ES, CJC, and FJPC analyzed the data. FJPC, AF, ES and VN wrote the manuscript, with input from all authors. All authors read and approved the final manuscript
Electronic supplementary material
Additional file 1:List of genes differentially expressed by 2-fold in experimental Group A (CLA supplementation during gestation and during suckling through dam's milk). Excel file containing the list of 2-fold differentially expressed genes in group A with respect to the control (Group D) generated using GeneSpring software. It includes the gene symbol of all genes and the associated description. The absolute fold change values relative to the control group and the type of regulation (Up or Down) are provided. The differentially expressed transcripts corresponding to open reading frames, transcribed sequences, cDNA clones or hypothetical genes were removed. (XLS 158 KB)
List of genes differentially expressed by 2-fold in experimental Group B (CLA supplementation during gestation and during suckling through oral gavage)
Additional file 2:. Excel file containing the list of 2-fold differentially expressed genes in group B with respect to the control generated as described in Additional file 1. (XLS 98 KB)
List of genes differentially expressed by 2-fold in experimental Group C (CLA supplementation during suckling through oral gavage)
Additional file 3:. Excel file containing the list of 2-fold differentially expressed genes in group C with respect to the control generated as described in Additional file 1. (XLS 88 KB)
About this article
Cite this article
Selga, E., Pérez-Cano, F.J., Franch, À. et al. Gene expression profiles in rat mesenteric lymph nodes upon supplementation with Conjugated Linoleic Acid during gestation and suckling. BMC Genomics 12, 182 (2011). https://doi.org/10.1186/1471-2164-12-182
- Conjugate Linoleic Acid
- Connective Tissue Growth Factor
- Mesenteric Lymph Node
- Conjugate Linoleic Acid Isomer
- Suckling Period