Gene expression profiles in rat mesenteric lymph nodes upon supplementation with Conjugated Linoleic Acid during gestation and suckling
Department of Biochemistry and Molecular Biology, School of Pharmacy, University of Barcelona, Barcelona, Spain
Department of Physiology, School of Pharmacy, University of Barcelona, Barcelona, Spain
CIBER in Epidemiology and Public Health (CIBERESP), Barcelona, Spain
Ordesa Group, Research Department, Scientific Park of Barcelona, Barcelona, Spain
Abstract
Background
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.
Results
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.
Conclusions
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.
Background
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
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
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
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
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
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.
Results
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
Figure 1
Experimental design
Experimental design. Total period of CLA supplementation (TPS) is shown in parentheses for each group of animals. Pups from dams fed with 1% CLA diet during the last two weeks of gestation and throughout the suckling period constituted Group A. Pups from dams fed only during the last two weeks of gestation with a 1% CLA diet and receiving CLA by daily oral gavage throughout the suckling period represented Group B. Pups from dams fed with a standard diet during the last two weeks of gestation and suckling and receiving CLA by daily oral gavage throughout the suckling period corresponded to Group C. Pups from dams fed with a standard diet throughout the study, constituted the reference diet group our Group D.
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.
Click here for file
Additional file 2
List of genes differentially expressed by 2-fold in experimental Group B (CLA supplementation during gestation and during suckling through oral gavage). 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
Click here for file
Additional file 3
List of genes differentially expressed by 2-fold in experimental Group C (CLA supplementation during suckling through oral gavage). 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
Click here for file
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
Figure 2
Venn Diagram of differentially expressed genes
Venn Diagram of differentially expressed genes. Venn diagrams in GeneSpring GX were used to compare the lists of genes differentially expressed by 2-fold (p < 0.05) in each of the three experimental groups of animals. The diagrams show the number of genes that belong to each of the individual lists, the genes in common between each pair of lists and the genes in common among all three lists (in the center of the representation) for each type of regulation (A: Upregulation and B: Downregulation).
Table 1
Common genes among the three experimental groups supplemented with CLA.
Gene Symbol
Description
Gene Ontology Biological Process
Ratio A
Ratio B
Ratio C
Actg2
actin, gamma 2, smooth muscle, enteric
70.79 D
97.48 D
105.85 D
Cnn1
calponin 1, basic, smooth muscle
actin cytoskeleton organization
17.18 D
20.13 D
19.90 D
Myh11
myosin, heavy chain 11, smooth muscle
smooth muscle contraction
16.49 D
16.57 D
19.68 D
Pcp4
Purkinje cell protein 4
12.37 D
12.81 D
12.94 D
Tpm2
tropomyosin 2
regulation of ATPase activity
7.89 D
7.26 D
6.84 D
Krt19
keratin 19
response to estrogen stimulus
7.78 D
5.59 D
9.12 D
Tpm1
Tropomyosin 1, alpha
muscle contraction
7.38 D
8.09 D
8.27 D
Cfb
complement factor B
complement activation
7.11 D
4.99 D
7.82 D
Acta1
actin, alpha 1, skeletal muscle
cell growth
7.10 D
7.54 D
7.18 D
Gal
galanin prepropeptide
inflammatory response
6.38 D
7.03 D
7.72 D
Rnase4
ribonuclease, RNase A family 4
proton transport
5.38 D
2.53 D
4.06 D
Upk1b
uroplakin 1B
epithelial cell differentiation
5.36 D
4.17 D
5.51 D
Rnase4
ribonuclease, RNase A family 4
proton transport
5.36 D
2.38 D
3.58 D
Mylk
myosin light chain kinase
5.04 D
4.11 D
3.23 D
Tpm1
tropomyosin 1, alpha
muscle contraction
4.97 D
5.34 D
5.80 D
Tnfaip6
tumor necrosis factor alpha induced 6
cell adhesion
4.94 D
3.35 D
4.39 D
Klhl23
kelch-like 23
4.62 D
2.99 D
2.93 D
Tnfrsf11b
tumor necrosis factor receptor 11b
signal transduction
4.37 D
2.02 D
2.46 D
Fhl1
four and a half LIM domains 1
cell differentiation
4.22 D
2.56 D
2.73 D
Synm
synemin, intermediate filament protein
3.96 D
3.69 D
4.74 D
Tmem100
transmembrane protein 100
3.76 D
3.06 D
3.71 D
Smoc2
SPARC related modular calcium binding 2
cell-substrate adhesion
3.68 D
5.14 D
7.01 D
Pln
phospholamban
calcium ion transport
3.68 D
3.66 D
3.89 D
Tm4sf1
transmembrane 4 L six family member 1
3.61 D
2.09 D
2.35 D
Acin1
apoptotic chromatin condensation inducer 1
chromosome condensation
3.59 U
3.09 U
2.69 U
Grem2
gremlin 2, cysteine knot superfamily, homolog
BMP signaling pathway
3.52 D
3.38 D
4.83 D
Pik3c2b
phosphoinositide-3-kinase, class 2, beta
cell communication
3.45 U
3.42 U
2.95 U
Fbxl22
F-box and leucine-rich repeat protein 22
3.36 D
3.72 D
3.67 D
Nupr1
nuclear protein 1
acute inflammatory response
3.31 D
3.00 D
3.01 D
Nov
nephroblastoma overexpressed gene
regulation of cell growth
3.27 D
3.10 D
3.48 D
Scg2
secretogranin II
MAPKKK cascade
3.22 D
3.16 D
3.15 D
Rbm5
RNA binding motif protein 5
nuclear mRNA splicing
3.14 U
3.99 U
4.04 U
Tpm1
tropomyosin 1, alpha
muscle contraction
3.11 D
3.29 D
3.27 D
Tpm1
tropomyosin 1, alpha
muscle contraction
3.11 D
4.63 D
4.57 D
Bzrap1
benzodiazapine receptor associated protein 1
3.08 U
3.61 U
3.25 U
Zeb2
Zinc finger E-box binding homeobox 2
regulation of transcription
3.07 U
5.80 U
2.06 U
Adamts1
ADAM metallopeptidase with thrombospondin 1
proteolysis
3.01 D
2.26 D
2.34 D
Myl9
myosin, light chain 9, regulatory
3.01 D
5.94 D
5.54 D
Gap43
growth associated protein 43
regulation of cell growth
2.98 D
3.34 D
3.22 D
Aldh1a1
aldehyde dehydrogenase 1 family, A1
metabolic process
2.94 D
2.77 D
2.91 D
Krt8
keratin 8
apoptosis
2.93 D
3.97 D
3.54 D
Gstm5
glutathione S-transferase, mu 5
metabolic process
2.86 D
3.14 D
2.81 D
Sh3bgr
SH3 domain binding glutamic acid-rich
2.85 D
3.32 D
3.40 D
Prph
peripherin
cytoskeleton organization
2.81 D
3.59 D
4.12 D
Ctgf
connective tissue growth factor
cartilage condensation
2.80 D
3.59 D
5.12 D
Sparcl1
SPARC-like 1
2.74 D
2.00 D
2.65 D
Pgm5
phosphoglucomutase 5
carbohydrate metabolism
2.74 D
3.86 D
4.72 D
Klhl23
kelch-like 23
2.74 D
2.26 D
2.16 D
Fhl1
four and a half LIM domains 1
cell differentiation
2.73 D
2.74 D
3.13 D
Acta2
smooth muscle alpha-actin
muscle contraction
2.70 D
4.63 D
2.86 D
Msln
mesothelin
cell adhesion
2.69 D
2.58 D
3.63 D
Sstr3
somatostatin receptor 3
signal transduction
2.67 U
3.11 U
2.90 U
Ppp1r14a
protein phosphatase 1, regulatory subunit 14A
phosphorylation
2.65 D
4.99 D
4.83 D
Timp1
TIMP metallopeptidase inhibitor 1
cell activation
2.64 D
4.04 D
4.54 D
Bcl2l2//Pabpn1
Bcl2-like 2//poly(A) binding protein nuclear 1
apoptosis
2.64 U
2.40 U
2.40 U
Ash1l
Ash1 (absent, small, or homeotic)-like
2.58 U
3.43 U
3.85 U
Mylk
myosin light chain kinase
2.57 D
3.10 D
3.94 D
Adh1
alcohol dehydrogenase 1
retinoid metabolism
2.52 D
2.52 D
2.56 D
Ptgis
prostaglandin I2 synthase
prostaglandin biosynthesis
2.48 D
2.92 D
3.33 D
Stmn2
stathmin-like 2
intracellular signaling
2.45 D
2.40 D
2.33 D
Gpd1
glycerol-3-phosphate dehydrogenase 1
gluconeogenesis
2.43 D
2.40 D
2.79 D
Pdlim3
PDZ and LIM domain 3
actin filament organization
2.42 D
2.33 D
2.50 D
Cgnl1
cingulin-like 1
2.41 D
2.01 D
2.15 D
Pla2g2a
phospholipase A2, group IIA
phospholipid metabolism
2.41 D
2.93 D
3.05 D
Ikzf2
IKAROS family zinc finger 2
2.39 U
5.10 U
4.82 U
Rgs4
regulator of G-protein signaling 4
signal transduction
2.39 D
2.27 D
2.55 D
Argbp2
Arg/Abl-interacting protein ArgBP2
intracellular signaling
2.36 D
2.30 D
2.45 D
Fgf13
fibroblast growth factor 13
MAPKKK cascade
2.35 D
2.80 D
2.86 D
Igf2
insulin-like growth factor 2
cell proliferation
2.30 D
4.82 D
6.01 D
Syt1
synaptotagmin I
transport
2.29 D
2.46 D
2.30 D
Hspb1
heat shock protein 1
response to heat
2.24 D
3.69 D
5.15 D
Gja5
gap junction protein, alpha 5
cell communication
2.24 D
3.04 D
2.93 D
Schip1
schwannomin interacting protein 1
2.23 D
2.92 D
2.80 D
Efemp1
EGF-containing fibulin-like extracellular matrix 1
2.23 D
2.41 D
3.39 D
Cxcr7
chemokine (C-X-C motif) receptor 7
signal transduction
2.22 D
2.24 D
3.13 D
Phemx
pan hematopoietic expression
2.21 U
2.56 U
2.61 U
Gpc3
glypican 3
regulation of growth
2.20 D
3.02 D
4.39 D
Leng8
leukocyte receptor cluster member 8
2.20 U
2.34 U
2.10 U
Crispld2
cysteine-rich secretory protein LCCL domain containing 2
lung development
2.20 D
3.82 D
3.77 D
Tagln
transgelin
cytoskeleton organization
2.20 D
3.49 D
3.54 D
Grb2
growth factor receptor bound protein 2
MAPKKK cascade
2.19 U
2.42 U
2.25 U
Parva
parvin, alpha
cell adhesion
2.16 D
2.04 D
2.30 D
Hist1h4b
histone cluster 1, H4b
nucleosome assembly
2.15 U
2.44 U
2.08 U
Des
desmin
2.12 D
3.69 D
4.05 D
Cryab
crystallin, alpha B
glucose metabolic process
2.09 D
3.37 D
3.18 D
Enpp3
ectonucleotidepyrophosphatase phosphodiesterase 3
phosphate metabolism
2.05 D
2.21 D
2.52 D
Bag2
Bcl2-associated athanogene 2
apoptosis
2.04 D
3.04 D
3.12 D
Wfdc1
WAP four-disulfide core domain 1
regulation of cell growth
2.01 D
2.66 D
2.80 D
Cd9
CD9 molecule
cell adhesion
2.01 D
2.24 D
2.12 D
The lists of genes differentially expressed by 2-fold with a p-value < 0.05 obtained for each experimental group were compared using Venn Diagrams in GeneSpring GX software v 10.0.2. The table shows the list of genes differentially expressed in all three conditions, and includes the Gene Symbol for all genes, their associated description and one of the Gene Onthology categories to which the genes belong according to GeneSpring GX. The ratio columns correspond to the absolute fold change in expression for the genes in each experimental group (A, B or C) relative to the control group (D) and the type of regulation (U, upregulation; D, downregulation)
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
Figure 3
BAN of differentially expressed genes
BAN of differentially expressed genes. The list of common genes among the three experimental groups of animals was the starting point to construct a BAN using the Pathway Architect software within GeneSpring GX. Expanded networks were constructed by setting an advanced filter that included the categories of binding, expression, metabolism, promoter binding, protein modification and regulation (see inset legend). Only proteins are represented. Genes present in the original list of common genes are encircled in blue, whereas the other genes were added by the software from the interactions database. The BAN presented shows some highly interconnected node genes that were object of further studies (pointed with arrows).
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
Figure 4
mRNA levels of selected node genes
mRNA levels of selected node genes. The mRNA levels of the indicated genes were determined for the four groups of animals by RT-Real Time PCR as described in Methods. Empty bars represent the mRNA levels for each specific gene as determined in the microarray experiments, filled bars represent the mRNA levels for each specific gene as determined by RT-Real Time PCR. ** p < 0.01 and *** p < 0.001.
Discussion
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
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
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
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
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
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
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
During early life CLA is able to enhance the production of immunoglobulins
Synaptotagmins (Syts) constitute a family of proteins that act as Ca2+-sensors for regulated exocytosis and endocytosis
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
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
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)
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
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
Conclusions
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.
Methods
Animals
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
Experimental Design
Animals were randomly distributed in 4 dietary groups, according to total period of CLA supplementation and administration route used in the pups
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
RNA extraction
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).
Microarrays
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
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.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
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
Acknowledgements
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.
Dietary regulation of intestinal gene expression
Dietary fatty acids influence the production of Th1- but not Th2-type cytokines
Fatty acids, inflamation and immunity
Polyunsaturated fatty acids and T-cell function: implications for the neonate
Conjugated linoleic acid concentrations of human milk and infant formula
The anticarcinogenic conjugated fatty acid c9, t11-c18:2, or rumenic acid, in human milk: amounts and effects
Fatty acid and conjugated linoleic acid isomer profiles in human milk fat
Conjugated linoleic acid is a growth factor for rats as shown by enhanced weight gain and improved feed efficiency
Mechanisms of action of conjugated linoleic acid: evidence and speculation
Conjugated linoleic acid (CLA) effects on pups growth, milk composition and lipogenic enzymes in lactating rats
Conjugated linoleic acid suppresses the growth of human breast adenocarcinoma cells in SCID mice
Conjugated linoleic acid inhibits proliferation and induces apoptosis of normal rat mammary epithelial cells in primary culture
Conjugated linoleic acid isomers and cancer
Mucosal IgA increase in rats by continuous CLA feeding during suckling and early infancy
Higher immunoglobulin production in conjugated linoleic acid-supplemented rats during gestation and suckling
Enhancement of antibody synthesis in rats by feeding cis-9,trans-11 conjugated linoleic acid during early life
Long-term feeding of the cis-9,trans-11 isomer of conjugated linoleic acid reinforces the specific immune response in rats
Nutrigenomics research for personalized nutrition and medicine
Integrated transcriptional profiling and linkage analysis for identification of genes underlying disease
An integrative metabolism approach identifies stearoyl-CoA desaturase as a target for an arachidonate-enriched diet
Dietary effects of arachidonate-rich fungal oil and fish oil on murine hepatic and hippocampal gene expression
Nutrigenomic approach to understanding the mechanisms by which dietary long-chain fatty acids induce gene signals and control mechanisms involved in carcinogenesis
Gene Expression Omnibus
Structural and functional properties of CCN proteins
Connective tissue growth factor: context-dependent functions and mechanisms of regulation
Functions and mechanisms of action of CCN matricellular proteins
Human skin fibroblast collagenase inhibitor. Comparative studies in human connective tissues, serum, and amniotic fluid
Expression analysis of the entire MMP and TIMP gene families during mouse tissue development
Tissue inhibitors of metalloproteinases in cell signaling: metalloproteinase-independent biological activities
Conjugated linoleic acid alters matrix metalloproteinases of metastatic mouse mammary tumor cells
Troglitazone but not conjugated linoleic acid reduces gene expression and activity of matrix-metalloproteinases-2 and -9 in PMA-differentiated THP-1 macrophages
Adeno-associated virus vector-mediated expression and constitutive secretion of galanin suppresses limbic seizure activity
Interactions between the enteric nervous system and the immune system: role of neuropeptides and nutrition
Neural intervention in mucosal immunity
Gender-related differences in murine T- and B-lymphocyte proliferative ability in response to in vivo [Met(5)]enkephalin administration
Bombesin, gastrin-releasing peptide, and neuromedin C modulate murine lymphocyte proliferation through adherent accessory cells and activate protein kinase C
The galanin signaling cascade is a candidate pathway regulating oncogenesis in human squamous cell carcinoma
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
Role of synaptotagmin in Ca2+-triggered exocytosis
Impaired membrane resealing and autoimmune myositis in synaptotagmin VII-deficient mice
Synaptotagmin: a calcium sensor on the synaptic vesicle surface
Loss, restoration, and maintenance of plasma membrane integrity
Synaptotagmin-mediated vesicle fusion regulates cell migration
The first C2 domain of synaptotagmin is required for exocytosis of insulin from pancreatic beta-cells: action of synaptotagmin at low micromolar calcium
Synaptotagmin III/VII isoforms mediate Ca2+-induced insulin secretion in pancreatic islet beta -cells
S100A13 is involved in the regulation of fibroblast growth factor-1 and p40 synaptotagmin-1 release in vitro
JunB is required for IgE-mediated degranulation and cytokine release of mast cells
Synaptotagmins: why so many?
Morphologically docked synaptic vesicles are reduced in synaptotagmin mutants of Drosophila
Synaptotagmin I: a major Ca2+ sensor for transmitter release at a central synapse
Glomerular podocytes contain neuron-like functional synaptic vesicles
Microarray analysis of hepatic gene expression identifies new genes involved in steatotic liver
Synaptotagmin I, a Ca2+ sensor for neurotransmitter release
synaptotagmin mutants reveal essential functions for the C2B domain in Ca2+-triggered fusion and recycling of synaptic vesicles in vivo
Galanin-like immunoreactivity in enteric neurons of the human colon
Immune pathology associated with altered actin cytoskeleton regulation
The cytoskeleton and cell signaling: component localization and mechanical coupling
Actin nucleation and elongation factors: mechanisms and interplay
TNF-alpha suppresses alpha-smooth muscle actin expression in human dermal fibroblasts: an implication for abnormal wound healing
A transforming growth factor-beta control element required for SM alpha-actin expression in vivo also partially mediates GKLF-dependent transcriptional repression
Cancer associated fibroblasts stimulated by transforming growth factor beta1 (TGF-beta 1) increase invasion rate of tumor cells: a population study
Adaptive immune regulation in the gut: T cell-dependent and T cell-independent IgA synthesis
The chemokine receptor CXCR7 is expressed on lymphatic endothelial cells during renal allograft rejection
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
Bovine whey protein concentrate supplementation modulates maturation of immune system in suckling rats
BIND--The Biomolecular Interaction Network Database
MINT: a Molecular INTeraction database