Dietary effects of arachidonate-rich fungal oil and fish oil on murine hepatic and hippocampal gene expression

Background The functions, actions, and regulation of tissue metabolism affected by the consumption of long chain polyunsaturated fatty acids (LC-PUFA) from fish oil and other sources remain poorly understood; particularly how LC-PUFAs affect transcription of genes involved in regulating metabolism. In the present work, mice were fed diets containing fish oil rich in eicosapentaenoic acid and docosahexaenoic acid, fungal oil rich in arachidonic acid, or the combination of both. Liver and hippocampus tissue were then analyzed through a combined gene expression- and lipid- profiling strategy in order to annotate the molecular functions and targets of dietary LC-PUFA. Results Using microarray technology, 329 and 356 dietary regulated transcripts were identified in the liver and hippocampus, respectively. All genes selected as differentially expressed were grouped by expression patterns through a combined k-means/hierarchical clustering approach, and annotated using gene ontology classifications. In the liver, groups of genes were linked to the transcription factors PPARα, HNFα, and SREBP-1; transcription factors known to control lipid metabolism. The pattern of differentially regulated genes, further supported with quantitative lipid profiling, suggested that the experimental diets increased hepatic β-oxidation and gluconeogenesis while decreasing fatty acid synthesis. Lastly, novel hippocampal gene changes were identified. Conclusions Examining the broad transcriptional effects of LC-PUFAs confirmed previously identified PUFA-mediated gene expression changes and identified novel gene targets. Gene expression profiling displayed a complex and diverse gene pattern underlying the biological response to dietary LC-PUFAs. The results of the studied dietary changes highlighted broad-spectrum effects on the major eukaryotic lipid metabolism transcription factors. Further focused studies, stemming from such transcriptomic data, will need to dissect the transcription factor signaling pathways to fully explain how fish oils and arachidonic acid achieve their specific effects on health.


Background
The presence and abundance of long chain polyunsaturated fatty acids (LC-PUFA) in humans diets have been the basis of scientific investigations for decades. Studies have yet to resolve all the mechanisms of action of LC-PUFA, in particular the remarkable phenotypic differences between populations that consume different quantities of fish and fish oil in their diets [1][2][3][4][5]. Dietary LC-PUFA affect processes including: growth, neurological development, lean and fat mass accretion, reproduction, innate and acquired immunity, infectious pathologies of viruses, bacteria and parasites; and the incidence and severity of virtually all chronic and degenerative diseases including cancer, atherosclerosis, stroke, arthritis, diabetes, osteoporosis, and neurodegenerative, inflammatory, and skin diseases [6][7][8][9][10][11]. Due to lack of mechanistic knowledge, predictions for optimal intake of nutritionally important LC-PUFA remain vague and dietary recommendations simply relate to prevention of overt deficiency. There is not even scientific consensus on whether LC-PUFA present in human milk should be included in infant formulas [12].
The most often cited action of n3 LC-PUFA is that they induce alterations in eicosanoids by antagonizing the primary substrate of eicosanoid metabolism: arachidonic acid (AA) [13]. Certainly, n3 LC-PUFA replace AA in cellular membranes to some extent and alter eicosanoid products, but is this the primary mechanism underlying the benefits of consuming fish? The n3 LC-PUFA alter postprandial lipoprotein production and clearance, and lower circulating triglyceride-rich lipoproteins [14]. N3 LC-PUFA induce alterations in membrane properties including ion transport, trafficking and vesicular transport, and alter the transcription of genes encoding proteins in lipid biosynthesis, lipid desaturation, and lipid oxidation [15][16][17][18][19].
Fish oil contains, among others, docosahexaenoic acid (DHA; 22:6n3) and eicosapentaenoic acid (EPA; 20:5n3) as important components. The concentration of DHA in the retina and neurological membranes suggests a specific role for DHA in the brain. However, neither the mechanism of uptake nor the precise role of DHA in the brain is currently well understood [20]. Rigorous elimination of DHA and its precursor 18:3n3 from animal diets leads to limited, but significant, phenotypic effects [21]. The basic biochemical pathway of DHA biosynthesis from its precursor 18:3n3 is still being resolved [22]. The EPA component of fish oil is linked to anti-inflammatory, antithrombotic, and generally, eicosanoid-mediated effects by antagonizing AA metabolism.
Arachidonic acid (AA; 20:4n6), either formed from n6 precursors in mammals or obtained directly from the diet, is believed to exert most of its biological activity via con-version to various eicosanoids. N6 LC-PUFA, and especially AA, are also now recognized to elicit significant physiological effects when consumed in the diet. The lack of significant commercial sources of AA has limited the study of this fatty acid (FA) as a dietary ingredient until relatively recently. Due to the commercialization of fungal biomass production to produce very high levels of AA as a storage oil, it is now possible to study this FA.
Until recently, LC-PUFA effects on gene transcription were thought primarily to be mediated by a single subfamily of orphan nuclear receptors -peroxisome proliferator activated receptors (PPARs); however, it is now becoming evident that FAs can affect many different genes either via direct interactions or indirectly through additional transcription factors including hepatic nuclear-4a (HNF-4a), nuclear factor kb (NF-kb), retinoid X receptor a (RXRa), sterol regulatory element binding protein-1c (SREBP-1c), and liver X receptors (LXRa and LXRb) [17]. Indeed, the enthusiasm for uncovering the biological pathways underlying the beneficial actions of LC-PUFA has revealed a story that is increasingly more complex than originally supposed [23], thereby making microarray technology an ideal platform to further decipher the many roles of these nutritional lipids.
In the present work, mice were fed control diets adequate in 18:2n6 and 18:3n3 but lacking LC-PUFA; or a diet enriched in fungal oil (FUNG) enriched in AA, the n6 elongation and desaturation product of 18:2n6; or a diet containing fish oil (FISH) enriched in 22:6n3 and 20:5n3, the major n3 LC-PUFA elongation and desaturation products of 18:3n3; or a diet containing both fungal oil and fish oil (FUNG+FISH). Thereafter, gene expression profiling was performed on two organs in parallel with quantitative metabolic profiling of a broad spectrum of liver FAs. Tissues of interest were liver, the major lipid metabolizing tissue; and brain, the major neurological tissue and principal site for accumulation and functionality of DHA (and other fish oil components). The hippocampus brain region was specifically chosen because of its importance in memory and learning [24]. Liver and brain are metabolically connected organs as both AA and DHA cross the blood brain barrier via several possible carriers. Diets were fed for 57 days to study chronic and sustained alterations in control of gene expression rather than acute effects. modulation (p = 0.17). Feed intake data showed no significant differences amongst groups (classic ANOVA considering the cages and diets as cofactors; p > 0.05). There were also no obvious differences in brain and liver weights (p > 0.05). By these criteria, differences in gene expression and metabolic profiles described herein were judged to be linked to the FAs present in FUNG and FISH oils on healthy, normal weight animals.
The incorporation of 20:4n6 and 22:4n6 following feeding of the AA-rich FUNG diet, and the incorporation of 20:5n3, 22:5n3, and 22:6n3 following feeding of the EPA/ DHA-rich FISH diet, into whole brain, hippocampal, and hepatic phospholipid (PL) pools (Table 1) was further evidence the diets were consumed, resulting in significant changes to PL acyl composition. Fish oil feeding is well established to increase n3 LC-PUFA in hepatic PL as well as whole brain PL [18]; whereas deficiency of n3 LC-PUFA can decrease 22:6n-3 in hippocampal PL [25]. DHA present in the brain can be formed from precursors in the liver or astrocytes, and intact DHA passes the blood brain barrier via selective transport mechanisms [26]. These data confirm that diets did result in actual increases in the brain of the appropriate FAs and set the stage for interpreting transcriptional responses observed in each organ. Importantly, our lipid profiling data revealed that FUNG and FISH FAs were not identically incorporated into individual PL in the two organs studied (partly shown in Table  1), and this may affect subsequent transcription factor sig-naling. For example, there was consistently more AA accumulating in various hepatic PL than in the hippocampus following AA feeding; and more DHA accumulated in hepatic cardiolipin (CL) than in hippocampal CL, whereas the reverse was true for phosphatidylserine/phosphatidylinositol (PS/PI). Thus, there are several explanations for why FUNG, FISH, and FUNG+FISH feeding resulted in different transcriptional profiles in liver and hippocampus, as discussed subsequently.

Gene Expression analysis
Differentially expressed genes induced by the dietary treatments were selected using a limit fold change (LFC) model, in addition to a more conventional confidence interval test [27]. Genes that satisfied the 5% LFC model and which lay above the 99.9% confidence interval of pooled replicates were selected as differentially expressed. It should also be noted that for the liver data, inter-individual variation was examined. All genes selected as differentially expressed using pooled data in the liver also laid outside the 99.0% confidence interval when compared to individual mouse variability, thereby lending additional confidence to the selection method used with both the liver and hippocampus data. According to these criteria, which involves examining all possible pair-wise comparisons between the four diets, 329 and 356 differentially expressed genes were selected in the liver and hippocampus, respectively. In the present document, genes that are discussed are followed by parentheses indicating fold changes according to the following convention: FUNG versus . Those fold changes that are significant by the 5% LFC model are further identified with an asterisk (*). It is important to note that the nature of the LFC model allows differentially expressed genes to be selected based on differences between the treatments, and therefore explains why some genes discussed in the text have no fold changes marked with an asterisk (as significance was observed between treatments and not versus the control diet). All liver genes discussed in the present manuscript are found in Figure 1 (green dots) and lie within the range of validation data (black squares). The high degree of concordance between RT-PCR and microarray data (86%-reported in detail else-where [27]) establishes a certain level of confidence in the results discussed throughout this manuscript.
To understand the trends in transcriptional response within and between the two organs studied, genes were annotated using GO classifications; a dynamic, structured, yet precisely defined system which characterizes genes based on defined categories of molecular function, biological process, and cellular component [28]. Since the GO classifications contain several thousand descriptors, many differentially expressed genes had a unique annotation. Record counts for most GO identifiers were not high, although a number of significant trends were observable. Venn diagrams were constructed for the top 10 GO descriptors for each category, according to record count, to

Figure 1
Graphical representation of the trend and spread of the hepatic microarray data. Blue dots indicate the highest fold change for each probe set on the Affymetrix Mu11K chip, where the highest fold change was identified by examining all fold changes in every possible pairwise comparison between diets. Those probe sets that were identified as absent (by Affymetrix software) across all dietary treatments were treated as unchanged and therefore not considered in the limit fold change (LFC) selection model. The trend within this dataset indicates that as absolute expression increases (along the x-axis), fold changes decrease, i.e. the chance of seeing large fold changes with highly expressed genes is minimal. The 5% LFC line of best fit (red line) takes this observed trend into account and identifies the upper 5% of those probe sets with the greatest fold changes within predefined bins (model parameters fully explained in [27]). The 329 probe sets that lie above the 5% LFC line are those genes defined as differentially regulated. Green circles represent the hepatic genes discussed in the current manuscript and black squares represent those genes validated by real time PCR. When examining the spread of the genes discussed and validated, it is apparent that genes were selected across the entire range of absolute values and not concentrated in one area. Furthermore, as the validation data (which had 86% concordance with microarray data; see [27]) lies across the entire range of absolute expression levels, a high degree of confidence is established for the genes discussed within the present manuscript.
visualize broad transcriptional responses induced by dietary LC-PUFA ( Figure 2). The interpretation of transcriptomic data, especially that stemming from a nutritional study such as this, benefits greatly from such a classification system. The diverse whole-body effects of LC-PUFA render individual gene interpretation difficult; however, considering groups of genes reveals preliminary information concerning the biological mechanisms responding to the addition of dietary LC-PUFA.
The GO analysis and specific gene expression patterns were interpreted to indicate that LC-PUFA-enriched oils affected liver metabolism through the regulation of genes directly involved in FA metabolism and lipid binding/ transport (Figure 2), as previously reported in rats [18]. When examining GO classifications for the hippocampus, such descriptors as cell signaling and cell cycle control are seen in response to dietary LC-PUFA. This result was expected, as the two organs have very different biochemical roles and import dietary FA via distinctly different means. The tissue functions were reflected in the transcriptomic profile. Although there were relatively few overlapping functions between the two organs, the most prominent overlapping function indicated that LC-PUFA was affecting genes related to transcription and DNA/RNA-binding.

Hepatic gene expression analysis and transcription factor signatures
It is difficult to interpret transcriptional profiling data alone, since dimerization, post-translational modifications, relative quantities of co-activators and co-repressors, and the formation of multiple transcriptional complexes binding to common promoter regions ultimately determines the binding activation potential of transcription factors. Nevertheless, clustering analysis of differentially regulated genes revealed numerous genes known to be regulated by specific transcription factors, thereby enabling several transcription factor signatures to be established. These signatures then provide clues concerning the various pathways through which FA may effect gene transcription. Additionally, FUNG and FISH diets regulated some transcription factors not (Table 2) previ-

Figure 2
Venn diagrams depicting the top 10 GO functional categories in the hippocampus and liver. Examining the top 10 GO classifications for cellular component (CC), biological process (BP), and molecular function (MF) indicated that functional commonalties (overlapping circles) in response to LC-PUFAs exist between the two organs, implying that molecular responses to dietary LC-PUFA occurs via common biological pathways rather than through identical genes. Additional organ-specific functions are also indicated (individual circles).

genes in response to dietary LC-PUFAs in (A) hippocampus and (B) liver. Table includes gene categories discussed in the manuscript. Genes are grouped by functional category, sorted alphabetically, and annotated with NCBI LocusLink. Average difference intensity (ADI) represents gene absolute expression, and was used to calculate fold changes. Differentially expressed genes were clustered with a combined k-means/hierarchical approach. Asterisk (*) indicates ADIs significantly different from control diet, using 5% LFC gene selection model equations (LFC hippocampus = 1.36 + (90/min ADI); LFC liver = 1.52+(100/min ADI)). (Continued)
ously associated with these dietary oils, thereby indicating novel LC-PUFA mediated molecular mechanisms of action. These include: Spot14 (Thrsp), transcription factor 7 T-cell specific (Tcf7), transcription factor 2 (Tcf2), nuclear transcription factor-Y b(NF-Y b; Nfyb transcript), D site albumin promoter binding protein (Dbp), general transcription factor II (Gtf2I), activating transcription factor 5 (Atf5), Msx-interacting-zinc finger (Miz1), and pre B-cell leukemia transcription factor 1 (Pbx1). In the present work, the focus will be on transcription factors previously identified as mediators of LC-PUFA actions, such as PPARs, SREBP, and HNFs.

Hepatic peroxisome proliferator-activated receptors (PPARs)
PPARs are nuclear hormone receptors that are activated by micromolar concentrations of various FAs and FA analogues, such as fibrates and thiazolidinediones [29]. This subfamily can be divided into three isotypes, designated PPARa, PPARd and PPARg, each with tissue-specific expression. PPARs form heterodimers with retinoid X receptor (RXR), and interact with the peroxisome proliferator response element (PPREs) in gene promoters. PPREs are direct repeats (DR) of a hexanucleotide sequence AGGT-CA separated by one nucleotide and are therefore referred to as a DR-1 response element [30]. PPARa and PPARg play critical roles in the catabolism and storage of FA, whereas the function of PPARd is less certain. PPARa is the predominant PPAR subtype expressed in liver, as confirmed presently. The other two subtypes were not likely expressed in the mouse liver, as negative values were obtained with the Affymetrix software. PPARa is reported to regulate genes involved in FA transport, synthesis and oxidation, glucose and lipid metabolism, ketogenesis and D5, D6, and D9-desaturation [17].
In spite of dramatic differences in the expression of a variety of PPAR-responsive genes, PPAR gene expression itself
A recent study examined effects of PPARa-specific induction in murine livers using cDNA microarrays [30]. Despite the different strain of mice examined (C57BL/6J) and the different genomic technology utilized for the examination of differentially expressed genes (cDNA microarrays), several genes were found differentially regulated in both experiments, indicating that the PPAR agonist and LC-PUFA affect gene expression through similar pathways. Genes such as Cd36, Cyp2b9, Hmgcs2, and Ech1 were induced more than 4-fold when Wy-14,643 (a PPA-Ra-specific agonist) was administered to mice for 2 weeks. These transcripts were similarly regulated in the current study, providing further indication that many of the differentially expressed genes were likely PPARa-activated.
Preliminary evidence suggests preferred ligands for PPAR activation are 18 C and 20 C FAs, as FAs shorter than 14 C or longer than 20 C will not fit correctly into the PPAR binding pocket and stabilize the AF-2 helix [17]. The net effect would be that PPAR co-activators are not recruited to relax nucleosomal DNA and allow transcription machinery to interact with various promoter elements. Therefore DHA (22:6:n3) may require retroconversion via boxidation to an 18 C or 20 C LC-PUFA prior to activating PPARs [17,33]. Without the use of purified FA it is difficult to distinguish whether the DHA in our fish oil preparation was retroconverted since the fish oil also provided the retroconversion products, 20:5n3 and 22:5n3 (Table 3). Likewise, AA was not provided as a pure FA component making comparisons to the control to detect retroconversion difficult. However, earlier work clearly shows that both DHA and AA can be retroconverted to 2 carbon shortened intermediates.

Hepatic sterol regulatory element-binding proteins (SREB-Ps)
SREBPs are critically important basic helix-loop helix (bHLH) transcription factors activated in response to variations in intracellular cholesterol levels that induce transcription of genes with sterol regulatory element (SRE) and E-box containing promoters. Within the SREBP promoter element, the liver X receptor (LXR)-responsive elements (LXRE) may have an important role [34]. LC-PUFA may affect SREBP-1c mediated transcription via several mechanisms. One mechanism is to antagonize activation of the LXR receptor leading to fewer LXR/RXR heterodimer complexes to bind the LXRE in the SREBP-1c promoter and less SREBP-1c-mediated transcription of SRE-containing genes [34]. The PPARa signaling cascade intertwines with SREBP pathways since LC-PUFA also activate PPARa, which can induce LXRa expression [17].
In the present experiment, saturated FA replaced AA and DHA to equilibrate total fat content. Yoshikawa et al. [34] found AA to be more potent than DHA (saturated and monounsaturated being ineffective) at inhibiting activation of SREBP-1c promoter and LXRE enhancer constructs, induced by co-expression of LXRa or -b. An increase in transcription of SREBP-1 may be associated with increased levels of the active SREBP-1 protein.
Several genes that may contain SRE promoters such as SREBP-1 and Acly, and Elovl2 and Cd36, were found in the same cluster. However, similar to PPRE-containing genes, not all SRE-containing genes co-clustered, indicating that not all SRE-containing genes were regulated identically by the three LC-PUFA diets. Overall, the number of regulated genes with SREs and PPREs indicates the breadth of molecular responses to dietary LC-PUFA, many of which appear to be co-orchestrated not only by these two transcription factors, but by additional elements as well.
Hepatic SREBP-1 was down regulated with FISH and FUNG+FISH, but importantly, not affected with FUNG (-1.1, -2.1*, -1.8*). Thus, chronic AA feeding alone does not appear, in our experiments, to interact with the SREBP-1 signaling cascade. It is not clear how AA is handled differently in the liver to prevent its interaction with SREBP. The probe set in the Affymetrix Murine 11 k chip could not differentiate between the -1a and -1c splicing variants of SREBP-1. There was however high sequence identity (93%) between the probe set sequence and rat ADD1 mR-NA, the murine homologue of human SREBP-1c. The SREBP-1c/ADD1 isoform regulates transcription of genes involved in FA biosynthesis and desaturation [37].
SREBP-1c, in addition to interacting with SREs, binds Ebox containing promoters [38]. E-box containing promoters were not examined so the current data analysis could have underestimated SREBP-1-linked effects. Some E-box or bHLH binding factors activate the same genes as SREBP-1. For example, both upstream stimulating factor (USF; an E-box binding factor) and SREBP bind the FA synthase (Fasn) E-box [38]. Hepatic Usf2 transcription was significantly down regulated with FISH, but not affected by FUNG or the combination diet using microarrays, but this result could not fully be confirmed using RT-PCR (1.2, -5.6*, 1.1 with microarray; 2.5, -1.1, 2.0 with RT-PCR).
The NF-Y transcription factor adds further complexity to transcriptional regulation since it can interact and modulate the expression of SREBP-1c target genes containing Ybox motifs, such as those found in Fasn [37]. Mutating the NF-Y binding site in Fasn relieved LC-PUFA-induced Fasn down regulation [39]. Nfyb was differentially expressed after dietary LC-PUFA addition (1.8, -1.6, 1.6). Although no significant changes were observed between dietary treatments and control diet, a significant difference was observed between the FUNG and FISH diets. Thus, LC-PUFA could have induced a joint SREBP-1c/NFY effect on common promoter elements.

Hepatic nuclear factors (HNFs)
Hepatic nuclear factors have been recognized to be transcription factors mediating the transcription of a variety of genes putatively involved in energy metabolism. HNFs may compete with PPARs for binding to DR-1 elements. Thus, LC-PUFA can activate PPARa/RXRa, leading to displacement of HNF4a from PPRE; or be converted to CoA derivatives, which in turn can suppress HNF4a transcriptional activity. Alternatively, these transcription factors may respond differently to chronic versus acute administration of LC-PUFA.
Hnf4a was not differentially regulated in the current experiment. The related transcription factor HNF3g (winged helix protein, fork head, Foxa3) was however down regulated -2.0* fold with FUNG+FISH, but not with the individual LC-PUFA. HNF3g protein may be involved in etiology of obesity, hyperlipidemia, and diabetes, by regulating glucagon transcription, insulin resistance, and pancreatic g cell function [40,41]. To our knowledge, the regulation of hepatic HNF3g expression by LC-PUFA has not been previously reported.
HNF3g-signaling may regulate transcription of Pck1 (up regulated ca. 2* fold with each diet, presently), and G6pc (differentially regulated in the current experiment) [42]. HNF3g is also reported to increase transcription of several genes that were not differentially regulated in the present experiment including, transferrin, tyrosine amino transferase [43], and IPF-1/PDX-1 [40].

Hepatic fatty acid synthesis
In the cytoplasm, FA synthesis is controlled via Acly, acetyl CoA carboxylase (ACC), and Fasn. The Acly reaction enhances FA synthesis by providing more acetyl CoA and NADPH substrates and enhancing ACC activity. The malonyl CoA generated in the ACC reaction then inhibits FA oxidation by reversibly inhibiting Cpt1.
It has previously been shown that the Acly gene in human cultured hepatic cells has a regulatory region from -61 to -49 in the 5' flanking region similar to that of Fasn, which is responsive to insulin and LC-PUFA [44]. In the present experiment, there was coordinated down regulation of both Acly (between -1.8* and -8.9* fold regulation for four different probe sets) and Fasn (-1.9*, -1.3, -5.4*; FUNG and FISH were not statistically significant from one another, but each was significantly different with the combination diet). It is intriguing that, in general, both Acly and Fasn were more strongly down regulated with the combination diet than the individual LC-PUFA. This synergy represents just one of the many unexpected findings revealed with a transcriptomic approach. Although ACC was not regulated at the transcriptional level, cytoplasmic citrate (generated from the Acly reaction) is an allosteric activator of ACC, so less cytoplasmic citrate would lead to less activation of ACC [45]. Transcript levels of hepatic acetyl CoA synthetase 1 AMP forming (acetate:CoA ligase AMP forming; Acas1) were consistently down regulated by LC-PUFA (-2.5*, -1.7*, -4.4*). This enzyme typically has specificity for 2-3 C FA chain lengths. In the liver, intestine, adipose, and mammary gland, the cytoplasmic form synthesizes acetyl CoA for lipogenesis [46]. It is regulated by SREBP-1a, -1c-and -2, as well as SP1 and SP3 binding to GC-boxes; the promoter also contains an E-box [47]. It is down regulated with fasting, and induced with re-feeding; such dietary regulation does not exist in SREBP-1 knockout mice [48]. More generally, acyl CoA synthase is reported to be PPARa-activated [49]. The observed down regulation is consistent with less acetyl CoA available for the ACC reaction, as well as a plethora of other reactions. Overall, it is plausible that both FUNG and FISH feeding decreased FA synthesis in the current experiment.

Hepatic FA oxidaton
As described above, Acly expression (and possibly ACC activity) was significantly down regulated by LC-PUFAs, likely leading to less FA synthesis. Malonyl CoA inhibits FA oxidation by reversibly inhibiting Cpt1. Cpt1 transfers acyl residues from CoA to carnitine in the intermembrane space of the inner mitochondrial membrane, and has a regulatory role in controlling hepatic b-oxidation [50]. Cpt2 transfers carnitinoyl residues to CoA in the mitochondrial matrix for b-oxidation. Carnitine acetyltransferase (Crat) may have the same role as CPT2, but with specificity for 2-10 C units.  [51]. The present data generally suggest LC-PUFA increased FA oxidation by down regulating Acly; and up regulating Cpt-1 and -2, Ech1 and Acadm. However, FUNG, FISH, and FUNG+FISH did not affect regulation of these transcripts identically.

Hepatic carbohydrate metabolism
Phosphoenolpyruvate carboxykinase (Pck1) is the ratelimiting enzyme of gluconeogenesis, and most of its activity is regulated transcriptionally. Feeding of all three LC-PUFA diets led to an increase in hepatic Pck1 transcription, indicating gluconeogenesis was favored over glycolysis. Specifically, four probe sets with high sequence similarity to rat Pck1 (accession numbers AA106463, AA110781, AA080172, AA063800) demonstrated differential gene expression for the mouse Pck1 homologue. The signaling pathways affecting Pck1 transcription are discussed in the section "cAMP signaling/links to carbohydrate metabolism".
G6pc converts glucose-6 phosphate to free glucose during gluconeogenesis. G6pc was unaffected with FUNG, significantly up regulated 1.6* fold with FISH (consistent with up regulated gluconeogenesis), but down regulated -1.9* fold with the combination diet (1.1, 1.6*, -1.9*). Rat G6pc has HNF4a consensus binding sites [52]. In HeLa cells and HepG2 hepatoma cells, PUFA down regulated G6pc, possibly via inhibition of HNF 4a promoter activity [52]. In our current mouse experiment, the transcriptional effects of LC-PUFA on G6pc may thus have been HNF4aor PPARa mediated.

Hepatic cAMP signaling/links to carbohydrate metabolism
cAMP is an important signaling molecule that can be found in specific cellular domains [53]. Levels are controlled through synthesis via adenylate cyclase, and degradation via phosphodiesterases (PDEs). Calmodulin PDEs (CaM-PDEs) are sensitive to calmodulin and calcium up regulation [54]. cAMP can then affect cell signaling via interactions with the regulatory subunits of cAMP activated kinases, including protein kinase A, leading to their subsequent activation, and via binding of cAMP response element binding proteins (CREBs) to cAMP response elements (CREs) or cAMP response units (CRU) consisting of CREs.
In the liver, cAMP regulates Pck1 transcription via the binding of CREB, as well as CCAAT/enhancer binding protein (C/EBP), to the CRU. The Pck1 promoter can thus exist in altered states of cAMP responsiveness, depending on the transcription factors bound to CRU elements [55]. In hepatocytes, SREBP-1c can block transcriptional induction of cAMP on PEPCK-C metabolic action [56].
In the present study, liver Cam-Pde 1c was not regulated by FUNG and FISH feeding alone, but was significantly up regulated 6.0* fold with FUNG+FISH together. A calcium/ calmodulin-dependent protein kinase 2 (CKII; Camk2b) was not affected with FUNG feeding or the combination diet, but was significantly down regulated -7.6* fold with FISH alone. CKII is a calmodulin target in the brain, involved in synaptic plasticity. In the liver, CKII may be associated with SAG protein (sensitive to apoptosis gene) and phosphorylate it. SAG in turn has a role in preventing oxidation by inhibiting cytochrome c release and caspase activation. Further, CKII may activate CPT-1 via cytoskeletal phosphorylation of cytokeratin. Thus, an activation of CKII may lead to increased b-oxidation and antioxidant protection. Feeding of all three LC-PUFA diets led to an approximately 2* fold increase in hepatic Pck1 for each probe set identified.
Overall, LC-PUFA may up regulate Pck1, the rate-limiting enzyme in gluconeogenesis, via cAMP-, SREBP-1-and/or PPAR-signaling [57]. This has the overall effect of limiting fat accumulation and shunting metabolic flux to gluconeogenesis. In a molecular sense, FUNG and FISH feeding mimics (or accentuates) the fasted metabolic state.

Hepatic Cytochrome p450 metabolism of LC-PUFA
There are at least 17 cytochrome p450 gene families (Cyps), with four hepatic types (Cyps 1-4) that metabolize foreign compounds and lipophilic substrates, including hormones, FAs, drugs and xenobiotics. Cyp 1 transcription is activated when cytosolic hydrocarbons activate the Ah receptor which then translocates to the nucleus and heterodimerizes with Arnt, leading to binding to an upstream xenobiotic response element on Cyp 1. Cyps 2-4 are induced through nuclear hormone receptors. Phenobarbitol and phenobarbitol-like substances activate CAR-RXR heterodimers (CAR, constitutive androstane receptor), that bind DR-4 response elements on the Cyp 2 promoter. Cyps 3 and 4 are activated by PXR-RXR (PXR, pregnane X receptor) and PPARa-RXR, respectively, also binding DR-4 response elements in both cases.
Nine hepatic cytochrome p450 family members were differentially regulated by LC-PUFA (Cyp1a2, Cyp2a4, Cyp2b9, Cyp3a11, Cyp3a16, Cyp4a10, Cyp4a14, Cypf13, and Por). Two broad clusters (hierarchical clustering alone) and 3 specific clusters (hierarchical-and k-means clustering) of co-expression exist within this subset, indicating transcriptional co-regulation (Table 2). Furthermore, additional specific clusters were close to each other, indicating only a fine split in expression behavior across diets.

Hepatic desaturation of LC-PUFA
As previously noted, SCD-1 is an important enzyme whose transcription is mediated by both PPAR-and SREBP-mediated signaling [59]. SCD-1 introduces cis double bonds at position C9 in FA chains. It has preferred substrate specificity for C18:0, but will also desaturate other saturated chain lengths. FUNG+FISH showed a significant down regulation of transcripts (-1.1, -1.0, -1.8* with microarray; 2.5*, -1.4, -2.7* with RT-PCR). In this instance, the microarray indicated no change with FUNG; however, RT-PCR indicated an up regulation. This may be expected as RT-PCR is more sensitive than microarrays. It is important to note however, that changes deemed significant by the microarray were also confirmed to be significant by RT-PCR, as is the case for FUNG+FISH induced changes in Scd1. The transcriptional changes in Scd1 expression were in some cases consistent with a marker for apparent decreased enzymatic activity, namely the ratios of FA in hepatic PL pools with a D9 double bond/same FA without a D9 double bond (Table 4). For example, in PS/ PI and phosphatidylethanolamine (PE), the ratio of 16:1n7/16:0 and 18:1n9/18:0 was lowest for the FUNG+FISH diet.

Hepatic elongation of LC-PUFA
There are suggestions in the literature that FA elongation may also be controlled by transcription factor signaling. The elongase Elovl2 may contain an SRE, and therefore be under the transcriptional control of SREBP [35]). Due to observations in the metabolite data, which suggests a repression of elongase activity, the dietary regulation of Elovl2 and Elovl3 was focused on. Changes in genes with known elongase activity were examined. The yeast gene Elo1 is associated with elongating FAs with 14 to 16 C; yet, this activity is often masked by cytosolic FA synthase activity [60]. Evidence shows Elovl2 elongates FA to 24C; and Elovl3 elongates FA to 26C, with preference for converting FA from 24C to 26C, in yeast [61].
The probe set for Elovl2 was down regulated, particularly with FISH and the combined diet (-1.3, -1.7*, -1.8*), consistent with some of the lipid data described below; whereas the probe set for Elovl3 was not affected by FUNG, up regulated with FISH, and down regulated, although not significantly, with the combination diet (-1.2, 1.9*, -1. 6). An NCBI BLAST search identified the murine homologue of yeast ELO2 (Elovl2) to be Ssc2 (Sequence similarity to Cig30), predominating in liver and testes [62]; complementation studies in yeast mutants indicated the functional equivalent of ELO3 (Elovl3) is Ssc1 [35]. Substrates of mouse Ssc2 were AA and EPA (DHA was not tested) [35]. Thus, FUNG and FISH feeding may transcriptionally regulate their subsequent conversion to elongated products.
Our metabolic data showed evidence that elongase activities may have been altered by LC-PUFA feeding. There was reduced conversion of 18 to 16C FA following LC-PUFA feeding. Specifically, the ratio of C18:0/16:0 in phosphatidylcholine (PC), PS/PI, PE and ether PE was decreased by FISH (

Hepatic Lipid Transport of FAs
Transcription factors such as PPARs and SREBPs are also known to have an effect on various lipid transport genes. Hepatic Abca1, Abca7 and Abce1 were differentially regulated; however, their roles in the liver are not precisely known. Several members of the ATP-binding cassette (ABC) superfamily, such as Abca1, are induced via PPARg and LXRa activation [63]. Abca1 has been shown to play a role in reverse cholesterol pathways and may be linked to HDL deficiency in metabolic syndromes [64]. In macrophages, 18:2n6 and AA increased ABCA1 protein degradation without affecting transcriptional levels [65]. In the current study, Abca1 was significantly decreased by FUNG and FISH.
FA translocase (Cd36) was significantly up regulated with FISH alone (-1.1, 2.4*, -1.5). Cd36 co-clustered with E-FABP, Apoe, Apoc1, and Abce1 transcript, all involved in FA and sterol transport. CD36 protein mediates long chain free FA uptake in adipocytes and other tissues, including liver. In murine pre-adipocytes, it is transcriptionally up regulated by long chain free FAs via a PPARamediated process [66]. Approaches using antisense technology, null mice, transgenic mice, and microarrays revealed important additional roles of CD36 protein as a receptor for thrombospondin, collagen type I, and oxidized LDL on macrophages [67]. The up regulation of Cd36 with FISH, and antagonism in presence of FUNG, merits further study.
FA binding proteins (FABPs) are small cytosolic proteins with virtually identical backbone structures that facilitate the solubility and intracellular transport of FAs. E-FABP is found in liver and other organs [68], but it has been most studied in psoriatic skin lesions [69] where it has putative roles in FA cytoplasmic transport, skin barrier function and keratinocyte differentiation [70]. Although links between PPARs and E-FABP have not currently been identified, this gene clustered with other known PPAR-affected genes suggesting a possibly important regulatory link. Purified E-FABP has five fold higher affinity for 18:0, than for 18:1n9 and 20:4n6 (22:6n3 not studied) [69]. Herein, hepatic E-FABP expression (Fabp5 transcript) was down regulated by all three diets (-5.5*, -8.4*, -13.6* with microarray; -2.9*, -8.5*, -16.4* with RT-PCR;). The significance of this large decrease in hepatic expression by both FUNG and FISH, and particularly the combination of the two, is a potentially important finding whose physiological consequence must be further investigated. Whether LC-PUFA induces the same decrease in Fabp5 in psoriatic skin is not known to our knowledge. N3 LC-PUFA, particularly 20:5n3, are known to be of clinical benefit in treating psoriasis [71].

Hepatic cholesterol metabolism, ketogenesis, and apolipoproteins
Mitochondrial Hmgcs2, which is involved in ketogenesis, was up regulated, particularly with the combination diet (1.2, 1.3, 1.6*). Mitochondrial Hmgcs2 is known to be induced by PPARa activation in experiments with hepatocytes [72,73]. The cytosolic form involved in cholesterol biosynthesis was not affected. Fpps was consistently down regulated (-2.4*, -1.6*, -1.9*) across diets, which could lead to reduced cholesterol synthesis; cholesterol synthesis and plasma cholesterol levels were not determined. Oils rich in DHA, such as fish oil, have cholesterol and triacylglycerol lowering properties [74].
Apoe and Apoc1, and Apoa4 and Apoc2, were similarly co-expressed/co-clustered across diets. When clustered with genes in the more general "transport" category, categories of co-expressed genes with underlying roles in lipid and sterol transport were evident. For example, Apoe and Apoc1 co-clustered with Abce1 and Cd36; and Apoa4 and Apoc2 co-clustered with Abca1 and Fabp5.
FUNG and FISH alone moderately down regulated Apoa4, whereas the combination diet resulted in a large down regulation (-2.6*, -3.2*, -17.2* with microarray; -1.3, -3.0, -4.9* with RT-PCR). A similar trend was seen for Apoa4 precursor transcript (-3.2*, -4.2*, -9.7*). ApoAIV protein is an anti-atherogenic component of high-density lipoprotein (HDL) particles and is also found in very lowdensity lipoproteins and chylomicrons [75]. It is involved in cholesterol absorption; low-density lipoprotein (LDL) clearance as a cofactor for hepatic lipase; increasing dolichol acylation; and increasing PE, but not PC, hydrolysis. FUNG and FISH may decrease Apoa4 expression through PPARa [76,77], and in the current experiment, LC-PUFA may have affected PPARa transactivation and signaling. ApoAIV activity could also be affected by diet induced structural changes to membrane PL in lipoprotein particles, as Apoa4 intercalates with PC in model membranes [78]. The biological consequences of reduced Apoa4 transcription, particularly with the combination diet, merits further investigation. Rodents fed low cholesterol diets, as in the present study, carry most of their cholesterol in the HDL fraction, and under these conditions, fish oil feeding reduces cholesterol from the HDL fraction, so the above mentioned changes in apolipoproteins, although fascinating, must be interpreted within the context of the chosen models.

FISH and FUNG-mediated effects in the Hippocampus
In the present study, 356 differentially expressed genes were identified in the hippocampus using a 5% LFC model. The trend in the data was similar to that seen with liver data, as shown in Figure 1 (see Materials and Methods for details). These data revealed that the consumption of dietary oils led to differences in the expression of a broad range of hippocampal genes that have been previously indicated to be involved in cell signaling and signal transduction, including cation transport, calcium binding, protein phosphorylation, protein kinases, and ATP-and GTP-binding as revealed by examining GO identifiers across all of the diet-regulated genes ( Figure 2). In striking contrast to the liver, a number of hippocampal LC-PUFAinduced differentially expressed genes had annotations suggesting immune related functions, as well as general defense responses specific to the hippocampus, including: "MHC-interacting proteins, lymphocyte antigens, immunoglobulins, and defense response." Additionally, there were numerous ontological assignments for genes which suggested dietary LC-PUFA feeding affected cell cycle control including the following categories: "cytoskeleton, cell cycle control, and cell cycle regulator".
In the hippocampus, dietary LC-PUFA activated a different set of transcription factors than observed in the liver, including: nuclear factor of k light chain gene enhancer in B cells 1 p105 (Nfkb1); nescient helix loop helix 1 (Nhlh1); E2F transcription factor 3 (E2f3); basic transcription factor 3 (Btf3); inhibitor of DNA binding 1 (Idb1); D site albumin promoter binding protein (Dbp); paired box gene 9 (Pax9); and the ELK1 member of the ETS oncogene family. Although dietary LC-PUFA activated numerous transcriptional factors in both liver and hippocampus, it is a noteworthy difference that the transcriptional signatures for PPAR, SREBP, and HNF could not be detected in the hippocampal dataset.
The interpretation of all 356 hippocampal transcripts is beyond the scope of the current manuscript, but all raw data is deposited at the NCBI gene expression Omnibus [79]. The present data provides evidence that LC-PUFA were incorporated into brain and hippocampal tissues (as previously discussed) and thereby influenced novel signaling pathways via receptor interactions. Some of these signaling pathways are highlighted below.

Hippocampal signaling pathways affected by dietary LC-PUFA
Transthyretin signaling Thyroid hormones including transthyretin are generally involved in neuronal proliferation and differentiation, and are required for normal cytoskeletal assembly. Transthyretin is specifically involved in transporting thyroxine (T4) and retinol-binding protein in cerebrospinal fluid and brain serum. AA has been found to bind transthyretin and inhibit thyroxin associations [80]. Transthyretin may also sequester b amyloid, thereby having neuroprotective properties. Dietary FISH had the effect of decreasing Ttr; whereas, FUNG and particularly the combination diet, increased Ttr (1.6, -1.9*, 2.2*). Similar to our findings in mice, feeding rats fish oil significantly decreased Ttr expression in whole brain 2.9 fold [18]. Interestingly, feeding of Ginkgo biloba extract (publicized to affect mental function) in mice increased hippocampal Ttr [81]. Levels of transthyretin, and associations of transthyretin with thyroxine and retinol-binding protein should be assessed in future LC-PUFA feeding experiments.
Alpha-2,8-linked disialic acid residues are found in glycoproteins [88]. During learning, there is a transient increase in neuronal polysialyation in the dentate gyrus of the hippocampus. This has been associated with selective retention and/or elimination of synapses that are transiently overproduced during memory consolidation [89]. In rat hippocampus, during development, PKCd can negatively regulate polysialyltransferase activity and neural cell adhesion molecule (NCAM) polysialylation state [89]. NCAM is implicated in neural differentiation and cellular plasticity [90]. Our observed decrease in Prkcd expression may have affected hippocampal polysialyltransferase activity. It is not known if Prkcd regulates expression of Siat8e.

Proopiomelanocortin (POMC)
POMC is a hypothalamic neuropeptide, and a target for leptin. POMC is cleaved post-translationally producing bioactive peptides (adrenocorticotrophin, b-endorphin and a-, band g-melanocyte stimulating hormones) that interact with neurons expressing melanocortin receptors, leading to appetite suppression [91]. Feeding high fat diets is known to increase Pomc expression [92]. Herein, FUNG feeding resulted in a 6.0* fold increase in hippocampal Pomc, whereas FISH and the combination diet had smaller, non-statistically significant effects (6.0*, 1.6, 1.4). The physiological effect of potentially increasing POMC protein with FUNG feeding is not known. As a next step, levels of POMC protein and downstream POMC peptides should be quantified

Conclusions and key findings
The molecular effects of dietary fish oil fed in reasonable quantities over prolonged periods has been found in humans and animals to be related to alterations in health outcome. In the present study, mice were fed relevant quantities of fish oil and AA rich fungal oils and the responses were studied in murine liver and hippocampus using a combined gene expression-and lipid-profiling strategy. The simple addition of these oils to the habitual diet of mice influenced hundreds of gene products [79]. These same conclusions regarding energy expenditure were reached in experiments assessing enzymatic activities, and transcriptional changes in rat liver using Northern blots, where fish oil feeding: increased palmitoyl CoA mitochondrial and peroxisomal oxidation rate; increased Cpt1, Cpt2, D 3 , D 2 -enoyl CoA isomerase, and Ech1 expression and activity; and decreased expression and activity of FA synthase, malic enzyme, glucose 6-phosphate dehydrogenase, and pyruvate kinase, compared to palm and/ or safflower oil feeding [93]. Likewise, in rat white retroperitoneal adipose tissue, fish oil feeding also decreased Fasn expression [94]. Intriguingly, the hormone leptin, which also acts through SREBP-1 and PPARa, affected many transcripts similarly to LC-PUFA (e.g., leptin decreased Gpam and Apoa4 transcripts; and up regulated Ech1, mitochondrial Hmgcs2, and Cyp4a14 transcripts [95].
Chronic feeding of fish and AA-enriched fungal oils broadly affected genes regulated by PPARa, HNFa, and SREBP-1 in the liver; transcription factors implicated in lipid metabolism as well as other processes. Only with future studies in which PPAR-specific agonists and antagonists are included, along with receptor antagonists, will the precise role of PPAR-signaling in response to LC-PUFA feeding be elucidated. The overlapping actions/effects of these various transcription factors present a challenge that will require tools, such as genomics and metabolomics, to provide a stepping stone towards the complete understanding of LC-PUFA metabolism.
There were striking differences in the genes and transcriptional factors activated by LC-PUFA in the liver and hippocampus, despite these organs being metabolically linked, as AA and DHA cross the blood brain barrier. There are at least three reasons for these inter-organ differences. Firstly, the dietary FAs were not identically incorporated into both tissues, as indicated by tissue lipid profiling. Secondly, different transcription factors exist and were activated in the two organs, which would affect downstream signaling. Thirdly, there is the possibility, that even when the same transcription factors were activated in both organs, factors including dimerization partners, multiple binding of various transcription factors, post-translational modifications, and signaling kinetics could affect the transcriptional profile.
Currently, nutritional FA recommendations are made largely by examining levels of key LC-PUFA in plasma; and less often, red blood cells, platelets, and biopsied tissue. Transcriptional gene profiles have not been considered in making such recommendations. The present work convincingly shows the potential importance of examining transcriptional profiles in addition to FA profiles; and the importance of examining more than one tissue, in making such LC-PUFA nutritional and clinical recommendations.

Experimental diets
All diets contained 90% fat-free AIN93G rodent diet in powder form (Lot 9350-5, Dyets, Inc., Bethlehem, PA), 0.4% milk fat, 1.2% palm olein, 1.9% Trisun sunflower oil, 1.5% soybean oil and 2.1-5.1% medium chain triacylglycerol oil (all of the above ingredients except the AIN93G were from Nestlé affiliated companies) (Tables 5  and 6). Part of the medium chain triacylglycerol oil in control diet was replaced with: 1.1% fungal oil (providing 0.5 dietary wt% AA and 1.0en%) in the FUNG diet; 1.9% fish oil (providing 0.5 dietary wt% DHA and 1.0en%) in the FISH diet; and 1.1% fungal oil and 1.9% fish oil in the combined FUNG+FISH diet. The levels of AA and DHA chosen were similar to that previously shown in rats to affect neurotransmitter levels and behavior [96]; and were not excessively high, being 2-3 fold higher than that recommended for human babies who have a slower D6 desaturase activity than rodents. Food was maintained at -80°C in daily aliquots under liquid N 2 , thawed each afternoon before administration to mice, and uneaten food was discarded daily.

Mice, feeding, and dissection conditions
Prior to actual experimentation, test dissections were made on identical mice to perfect the dissection techniques at Porsolt & Partners Pharmacology (Paris, France).
Male Rj:NMRI mice from Elevage Janvier, Le Genest-Saint-Isle France, weighing 10-11 g at delivery and 33-51 grams on experimental d 42, were housed 10 per cage in wire cages with bedding and normal light cycle. Mice received ad libitum quantities of distilled water and purified powdered diets (7.5 g/mouse) in ceramic cups (10/group) from d 1-57 or 58, and were thereafter sacrificed in the morning by cervical dislocation over these last 2 experi- mental days. Food was not removed the night before, and since nocturnal mice are known to eat during the normal period dark cycle employed, the mice could have been fasted 0-6 h at sacrifice. Mice were thus not completely fasted at sacrifice, a complete fast being 12-16 h.
After feeding diets to 10 mice per group, 5 mice were randomly selected for inclusion in hepatic and hippocampal gene expression analysis. From these 5 mice, livers were dissected according to standard protocols (Porsolt & Partners Pharmacology, Boulogne-Billancourt, France), which consisted of carefully dissecting out the liver, rinsing with 70% EtOH to remove adhering blood, trimming of adhering connective tissue with scissors, cutting into 100-150 mg subsections, freezing in liquid N 2 , and storing at -80°C until gene expression and lipid analysis. From these same 5 mice, the hippocampus was dissected from the mice according to standard protocols (Porsolt & Partners Pharmacology, Boulogne-Billancourt, France), rinsed as described for livers, cut into 17-50 mg subsections, and frozen as described for livers. Hepatic and hippocampal lipid analysis was performed according to standard methods at Lipomics, Inc. (U.S.A) on individual mice (n = 5).

Nucleic acid preparation
Liver tissue was extracted from 5 individual mice per group and extracted separately using Qiagen RNeasy mini-kits (Basel, Switzerland) according to manufacturer instructions, except that RNeasy columns were impregnated with DNase I (Roche, Basel, Switzerland) to remove any genomic DNA contamination.
For the group fed a control diet, two independent hepatic pools (n = 2) of 10 mg each were created, consisting of equal amounts (2 mg) of total RNA from each of the 5 aforementioned mice. This was done to assess measurement variability. Secondly, a separate piece of liver was excised from three of the individual control mice (n = 3) to additionally and independently extract a second set of RNA; these three mice were chosen at random from the group of five. This second round of extraction was used to independently assess the variability of normal control mice housed under the prevailing conditions of the study.
For each of the other dietary groups, a single pool (n = 1) of 10 mg total hepatic RNA was created for comparison against the control. Hippocampal RNA was prepared as described above for hepatic RNA, except only 1 pooled control RNA sample from 5 mice was prepared; hence All RNA samples were quantified with the RiboGreen RNA Quantification Kit (Molecular Probes, Eugene Oregon), then analyzed via agarose gel electrophoresis for intact 18S and 28S rRNA. All study samples were judged to contain sufficient amounts of high-quality RNA for hybridization to GeneChips. As another quality control measure, prior to hybridization with Affymetrix Gene-Chips (Affymetrix, Inc., Santa Clara, CA), we confirmed that pooled samples gave strong signals using Affymetrix test chips.

Gene expression analysis using the murine 11k GeneChip cRNA preparation
For the control diet group, each of three individual mouse liver RNA samples (n = 3) were hybridized to Affymetrix GeneChips. For both the control and LC-PUFA diets, 15 mg pooled total RNA was the starting material for hybridization to the Mu11kA & B GeneChips. Pooled liver controls were evaluated in duplicate on the Affymetrix probe array cartridge to estimate experimental variation in gene chip hybridization, as previously described in detail [27].
In all cases, total RNA was converted to biotinylated cR-NA, hybridized in the Affymetrix probe array cartridge, stained, and then quantified. First and second strand cDNA synthesis was performed using the SuperScript Choice System (Invitrogen AG, Basel, Switzerland), according to manufacturer instructions, but using an oligo-dT primer containing a T7 RNA polymerase binding site. Labeled cRNA was prepared with the MEGAscript In Vitro Transcription kit (Ambion, Inc., Austin, TX). Biotinylated CTP and UTP (Enzo Biochem Inc., NY) were used together with unlabeled dNTPs in the reaction, and unincorporated nucleotides were removed with RNeasy columns.
Array hybridization and scanning cRNA (10 mg) was fragmented at 94°C for 35 min in buffer containing 40 mM Tris-acetate pH 8.1, 100mM KOAc, 30 mM MgOAc. Prior to hybridization, fragmented cRNA in a 6X SSPE-T hybridization buffer (6X: 1M NaCl, 10mM Tris pH 7.6, 0.005% Triton), was heated to 95°C for 5 min, cooled to 40°C and loaded onto an Affymetrix probe array cartridge. The probe array was incubated for 16 h at 40°C at constant rotation (60 rpm), then exposed to 10 washes in 6X SSPE-T at 25°C followed by 4 washes in 0.5X SSPE-T at 50°C. Biotinylated cRNA was stained with 10 g/ mL streptavidin-phycoerythrin conjugate (Molecular Probes) in 6X SSPE-T for 30 min at 25°C followed by 10 washes in 6xSSPE-T at 25°C. Probe arrays were scanned at 560 nm using a confocal laser scanning microscope (made for Affymetrix by Hewlett-Packard). Readings from the quantitative scanning were analyzed with Affymetrix Gene Expression Analysis Software.

Selection of differentially regulated genes and data analysis
A mathematical method was developed and applied to the raw GeneChip data for the selection of differentially regulated genes. This method moves beyond setting a single fold change cut-off by considering fold change in the context of absolute expression [27]. In brief, this method includes: (A) the determination of the upper 5% of fold changes within narrow bins of absolute expression levels to determine a limit fold change (LFC) function; and (B) subsequent ranking of genes by a combined fold change/ absolute difference calculation. Herein, absolute expression is referred to as average difference intensity (ADI) in all mathematical equations. The function describing the upper 5% of hippocampal and liver data was described with the following equations: LFC hippocampus = 1.36 + (90/min ADI); and LFC liver = 1.52+(100/min ADI), respectively.
The selection of differentially expressed genes was corroborated in three ways for liver samples: (1) assessing the measurement variability and (2) assessing individual mouse variability in order to establish confidence intervals with which to calculate the significance of differences with pooled diet response groups and (3) confirmation by RT-PCR for a sub-selection of genes. Replicate analysis established that all genes selected by this method also laid above a more conventional 99.9% confidence interval, when compared to the variation of pooled replicates (n = 2). Subsequent analysis of the variation between individual mice (n = 3) indicated that this variation is comparable to the pooled variation and that all genes selected as differentially expressed lied above the 99.0% confidence interval for individual mouse variability. The gene expression means of the three individual mice were calculated and compared to pooled data in a scatter plot. A correlation co-efficient of 0.93 was obtained indicating that the pools are accurate representations of the individual mice used to create them. Lastly, validation studies showed up to 85.7% concordance of microarray data with real time polymerase chain reaction (RT-PCR) across the range of absolute expression for specific genes selected as differentially regulated [27]. All RT-PCR data, including information about the low false positive rate, for all of the liver samples mentioned in the present work, can be found in the aforementioned reference.
The 5% LFC method indicated that 513 genes in the hippocampus and 489 genes in the liver were differentially regulated by LC-PUFA from the original 13,179 genes represented on the Mu11K Affymetrix GeneChip. Absence/ presence calls were not considered a priori in selecting significantly regulated genes, but were a post-selection criterion. After removing those genes that were identified as absent across all dietary conditions, 356 hippocampal and 329 hepatic differentially expressed probe sets remained in the list for further analysis and interpretation. Of these, 221 and 214 genes represented annotated sequences in the hippocampus and liver, respectively.
Fold changes are always listed in the respective order: FUNG, FISH, and FUNG+FISH without repeating these labels; and negative signs indicate down regulation. To determine if the ADI for a treatment group was significantly different from the control diet, the 5% LFC best-fit equations were utilized: LFC hippocampus = 1.36 + (90/min ADI); LFC liver = 1.52+(100/min ADI). Groups significantly different from the control, as determined by a pair-wise comparison, are labeled with an asterisk in the text and Table 2. The mathematical model utilized also selected some genes due to significant differences between experimental treatments, even when no significant differences versus the control were observed. For reasons of simplicity, significant differences between experimental treatments are not indicated in Table 2, but are discussed in the text.

Data analysis
Selected genes were grouped by expression patterns through a combined k-means and hierarchical clustering approach (Spotfire, Inc., Sommerville, MA) by setting the maximum number of clusters in either method to 5 (a compromise between not having too many or too few clusters) and then graphing on a 2-D gene plot (Figures not shown; hierarchical cluster numbers (h1-h5), and kmeans cluster number (k1-5) included in Table 2). Genes were then annotated using gene ontology (GO) classifications; a dynamic, structured, yet precisely defined system which characterizes genes based on defined categories of molecular function, biological process, and cellular component [28]. Finally, genes were assigned names through NCBI Locuslink (Table 2).

Gene Confirmation
A subset of the genes selected by the 5% LFC model were examined by RT-PCR using SYBR ® green dye (Molecular Probes) which interacts with double-stranded nucleotide sequences. Data represented means of triplicates for each experimental treatment using pooled RNA samples (n = 5). Amplification was performed using an ABI 5700 machine (Applied Biosystems, Foster City, CA, USA) with a hot start at 95°C for 10 minutes, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min for denaturation, annealing and elongation. Data was normalized to b-actin and GAPDH, which were not statistically, differentially expressed. A students t-test (a = 0.05) was used to confirm significantly differentially expressed genes following RT-PCR.

Lipid analysis
Following sacrifice of the animals, the livers (250 mg piece), hippocampus (from half of the mice), and whole brains (from the other half of mice) from each animal were rapidly frozen in liquid N 2 and stored at -80°C. Each whole brain weighed between 395-454 mg. A slurry of each whole brain was made by thawing the tissue samples at room temperature for the minimum time necessary, adding 12.6 mL of CHCL3:MeOH (2:1, v/v), briefly homogenizing on ice, and re-freezing the samples under N 2 at -80°C (concentration was ca. 30 mg brain/mL solvent). Livers, hippocampus, and whole brain (4 individual samples/tissue) were quantitatively analyzed for FA distribution (mg FA/g tissue) in each PL class, by Lipomics Technologies, Inc. (West Sacramento, CA). Brain samples were also separately extracted, and N-acylethanolamines (NAEs) and monoacylglycerols (MAGs) purified and characterized as previously described in detail [84].