Gene Expression Changes in the Septum: Possible Implications for MicroRNAs in Sculpting the Maternal Brain

Abstract

The transition from the non-maternal to the maternal state is characterized by a variety of CNS alterations that support the care of offspring. The septum (including lateral and medial portions) is a brain region previously linked to various emotional and motivational processes, including maternal care. In this study, we used microarrays (PLIER algorithm) to examine gene expression changes in the septum of postpartum mice and employed gene set enrichment analysis (GSEA) to identify possible regulators of altered gene expression. Genes of interest identified as differentially regulated with microarray analysis were validated with quantitative real-time PCR. We found that fatty acid binding protein 7 (Fabp7) and galanin (Gal) were downregulated, whereas insulin-like growth factor binding protein 3 (Igfbp3) was upregulated in postpartum mice compared to virgin females. These genes were previously found to be differentially regulated in other brain regions during lactation. We also identified altered expression of novel genes not previously linked to maternal behavior, but that could play a role in postpartum processes, including glutamate-ammonia ligase (Glul) and somatostatin receptor 1 (Sstr1) (both upregulated in postpartum). Genes implicated in metabolism, cell differentiation, or proliferation also exhibited altered expression. Unexpectedly, enrichment analysis revealed a high number of microRNAs, transcription factors, or conserved binding sites (177 with corrected P-value <0.05) that were significantly linked to maternal upregulated genes, while none were linked to downregulated genes. MicroRNAs have been linked to placenta and mammary gland development, but this is the first indication they may also play a key role in sculpting the maternal brain. Together, this study provides new insights into genes (along with possible mechanisms for their regulation) that are involved in septum-mediated adaptations during the postpartum period.

About the Authors:

Changjiu Zhao

* E-mail: [email protected]

Affiliation: Department of Zoology, University of Wisconsin-Madison, Madison, Wisconsin, United States of America

Michael C. Saul

Affiliation: Department of Zoology, University of Wisconsin-Madison, Madison, Wisconsin, United States of America

Terri Driessen

Affiliation: Department of Zoology, University of Wisconsin-Madison, Madison, Wisconsin, United States of America

Stephen C. Gammie

Affiliations Department of Zoology, University of Wisconsin-Madison, Madison, Wisconsin, United States of America, Neuroscience Training Program, University of Wisconsin-Madison, Madison, Wisconsin, United States of America

Introduction

A variety of behavioral, metabolic, structural, and neuroendocrine alterations occur during the transition from a virgin to a lactating state [1]–[4]. These physiological adaptations during lactation are believed to be critically important for the survival and development of the offspring. For example, in rodents, within hours of parturition, the mother retrieves the scattered pups, gathers them together in the nest, and adopts a nursing posture over the pups to permit suckling [5]. Adaptive changes in several neuroendocrine systems during lactation are also evident. For example, the release of both adrenocorticotropic hormone and corticosterone in response to various stressors is markedly suppressed [6], [7]. A striking change in metabolic aspect during lactation is the negative energy balance resulting from the increased energy demand and decreased adaptive thermogenesis [8]. This energy-saving mechanism facilitates the availability of energy for milk production.

The mammalian septum is a heterogeneous forebrain structure that is divided into two main subdivisions, the medial septum (MS) and the lateral septum (LS), and has been linked to various aspects of maternal care, including offspring protection and pup retrieval [9]–[12]. Neuroanatomically, MS is densely interconnected with the hippocampus while LS has extensive reciprocal connections with numerous brain regions known to regulate affect and motivation, such as hypothalamus, amygdala, bed nucleus of the stria terminalis, periaqueductal gray, ventral tegmental area, and raphe nuclei [13]–[15]. Neurochemically, a high number of neurons in LS contain GABA and in subregions of LS, neurons also express other neurotransmitters, such as neurotensin, enkephalin, galanin, and somatostatin [13], [14]. Neurons within MS are predominantly GABAergic and/or cholinergic [16]. LS has long been considered to play a critical role in regulating multiple affective, behavioral and cognitive processes, such as fear, anxiety, depression, aggression, maternal behavior, and social recognition [11], [15], [17]–[19]. MS functions mainly as a region that modulates processes related to attention and memory [20], [21].

Alterations in neuronal activity and gene expression in several brain regions during the postpartum period have been reported. For example, during lactation, neuronal activity of neuropeptide Y (NPY) and agouti-related protein are increased in the arcuate nucleus of hypothalamus (Arc), a site known to be a core feeding center, whereas proopiomelanocortin is reduced [22], [23]. These changes in neuropeptide activity have been proposed to be important in integrating Arc-mediated food intake and energy balance. Expression changes of some genes in forebrain structures during lactation have been demonstrated to be involved in the control of maternal/parental and/or reproductive behaviors [24]–[26]. However, possible alterations in gene expression in the septum that may contribute to the septum-mediated adaptations during lactation remain largely unexplored. In this study, we employed microarray analysis to identify genes in the septum with altered expression during the postpartum period.

As part of this study, we were also interested in how identified gene expression differences may have occurred. Producing the maternal brain is a complex process and is associated with various contributing factors including hormonal changes and sensory input from different events, such as the mating, parturition, nursing, and additional interactions with offspring [27]–[29]. There are likely large numbers of transcription factors that contribute to the gene expression changes but which have not to date been evaluated. Newer data mining tools, such as gene set enrichment analysis (GSEA), provide a unique approach to gain insight into how altered gene expression occurs [30]. The databases include for each gene information on all known binding sites in that region, including those for a wide range of transcription factors and microRNAs (miRNAs), that may be involved in gene expression. MiRNAs are of interest because they have been implicated in peripheral alterations, such as placenta and mammary gland development [31], [32] in the maternal female, but to date no roles for miRNAs have been examined for the maternal brain. In this study, we used GSEA to gain insights into how large scale changes may have been occurring during pregnancy and early lactation with a subset of transcriptional regulators.

Results

Genes Identified as being Differentially Expressed in the Septum of Lactating Versus Virgin Mice

Using PLIER analysis with a P<0.0025 cutoff value, we identified 116 genes that displayed significant changes in expression (Table 1). Two genes of interest, Gal (P = 0.018, fold change is 0.83) and Sstrl (P = 0.004, fold change is 1.22) were not included in Table 1 because their P values were more than 0.0025. The full list of all 35,557 targets, their relative expression and P value ranking using PLIER algorithm is presented in Table S1. The genes identified as being differentially expressed in lactating versus virgin mice were distributed across a number of categories (Table 1). The function/category of each gene was sorted out individually based on the biological process using PubMed and GenBank databases.

[Figure omitted. See PDF.]

Table 1. List of genes showing highest significant differences in gene expression between lactating and virgin mice when analyzed with the probe logarithmic intensity error statistics (PLIER) (P<0.0025).

https://doi.org/10.1371/journal.pone.0038602.t001

Verification of Microarray Results by Quantitative Real-time PCR Analysis

Consistent with the findings from microarray analysis, mRNA levels of three representative genes, Glul (P = 0.023), Igfbp3 (P = 0.004) and Sstr1 (P = 0.048) were identified as being significantly up-regulated in maternal mice as compared to virgin mice, whereas Fabp7 (P<0.001) and Gal (P = 0.023) were significantly down-regulated (Fig. 1). While Npy1r mRNA levels were found to be significantly increased in lactating versus virgin mice in the microarray analysis, we did not confirm a difference between the two groups using highly sensitive real-time PCR technique (Fig. 1).

[Figure omitted. See PDF.]

Figure 1. Quantitative real-time PCR analysis of Fabp7, Gal, Glul, Igfbp3, Npy1r and Sstr1 expression in the septum.

Relative expression distribution of mRNA (y-axis) represented as a ratio of maternal versus virgin mice, was calculated using Ywhaz and CycA as reference genes, and shown by box-and-whisker plot as medians (dashed lines), interquartile range (boxes) and ranges (whiskers). Ratios over one indicate genes with higher expression in maternal relative to virgin mice, and ratios of less than one indicate genes with lower expression in maternal as opposed to virgin mice. Confirming array results, postpartum mice exhibited increased Glul, Igfbp3, Sstr1and decreased Fabp7 and Gal mRNA levels relative to virgin control mice. Npy1 was not confirmed. *P<0.05 maternal mice versus virgin control.

https://doi.org/10.1371/journal.pone.0038602.g001

Gene Set Enrichment Analysis Results

When evaluating genes upregulated in postpartum mice, 177 unique regions involved in transcriptional regulation reached significance (corrected P-value <0.05). The majority of these were miRNA binding sites (67%; 119 out of 177). For example, the site with the fifth highest p-value was a binding site for miRNA, Mir-183. For some binding sites, multiple miRNAs were listed. For example, for the 10th most significant region is the conserved binding site of AGCACTT, which can be bound by multiple miRNAs, including Mir-93 and Mir-302B. A partial (Table 2) and full list (Table S2) of the GSEA results are provided. The binding sites for 29 different transcription factors were identified as significantly upregulating gene expression in maternal mice (Table 3, Table S2). The database includes target data from multiple species and in some cases a given transcription factor (e.g., YY1) reached significance with multiple species. Here we only counted each transcription factor once and only list those associated with mammals (e.g., mouse, human). Twenty seven additional conserved binding sites in close proximity to genes were also identified as being linked with upregulated expression in maternal mice (Table S2). For these regions, no known transcription factor or miRNA has been identified that target these regions.

[Figure omitted. See PDF.]

Table 2. List of all miRNAs found to be significantly linked with upregulated gene expression in postpartum mice using GSEA.

https://doi.org/10.1371/journal.pone.0038602.t002

[Figure omitted. See PDF.]

Table 3. List of transcription factors (and corrected p-values) that are significantly linked to upregulation of subsets of genes in maternal mice using GSEA.

https://doi.org/10.1371/journal.pone.0038602.t003

Discussion

This study used Affymetrix microarray analysis and real-time PCR techniques to examine CNS gene expression changes in the septum during lactation. We identified 116 genes (P<0.0025 as analyzed by PLIER) and 2844 genes (using a lower cutoff of P<0.05) which were distributed across a number of categories in biological process, as being differentially expressed in lactating versus virgin mice. Some of the genes identified are consistent with previous studies testing gene expression changes during lactation. Our study also found a list of novel genes that may play a role in the transition from a non-parental to a parental state. Further, using data mining approaches and GSEA we identified a number of miRNAs and transcription factors that may play a role in producing the maternal brain.

Identification of Expression Changes During Lactation Consistent with Previous Studies

In the septum, mRNA of Igfbp3 increased about 30% during lactation (Table 1), and was confirmed by qPCR (Fig. 1). This result is consistent with previous findings that Igfbp3 was increased in the arcuate nucleus of hypothalamus (Arc)/ventromedial nucleus of hypothalamus (VMH) [33], and in large hypothalamic CNS regions in previous maternal array studies [34]. Igfbp3 is a member of the superfamily of insulin-like growth factor (IGF) binding proteins and is expressed in the brain. In vivo studies show that Igfbp3 inhibits IGF action by modulating the bioavailability and distribution of IGFs [35]. One possibility is that increased Igfbp3 in the septum during lactation would result in reduced IGF action, but this would need to be determined. Given that IGF negatively regulates food intake and body weight [36], the increased expression of Igfbp3 would be expected to increase food intake through an IGF-mediated mechanism, consistent with the low levels of IGF in association with the metabolic adaptations that occur during lactation [37]. Based on its role in neurogenesis [38], Igfbp3 has been linked to the involvement in certain disease processes, such as ischemia, hypoxia and autism [39], [40]. There are multiple pathways by which Igfbp3 could possibly interact with maternal circuitry.

Fabp7 transcript was reduced about 30% during lactation in this study and was also found to be significantly decreased in maternal rodents in Arc [33] and hypothalamic CNS regions [34] in previous array studies. The mammalian Fabp7 is expressed in radial glial cells, astrocytes, and neuronal cell precursors, functioning as a modulator of cell differentiation, proliferation, and neurogenesis [41]. Interestingly, Fabp7 was preferentially expressed in numerous astrocytes in the amygdala and septum, regions known to be critically involved in the regulation of fear and anxiety [42]. In addition, previous studies show that this gene modulates the processes of sleep, memory, and anxiety [42], [43]. Several lines of evidence associate Fabp7 with neuropsychiatric diseases such as schizophrenia [44]. Fabp7-deficient mice displayed enhanced fear memory and anxiety in adulthood [42]. Together with the functional significance of the septum in fear, anxiety and memory [15], [20], [21], these findings suggest the possibility that Fabp7 plays an important role in emotional changes that occur during the postpartum period.

We found a decrease of galanin in the gene array, which was confirmed by qPCR. Significant reduction of galanin during lactation has been observed in previous array studies [33], [34]. Galanin is a 29 amino-acid (30 in humans) neuropeptide and is richly expressed in the CNS. Although galanin was not listed in the top identified genes in terms of significance (P = 0.018), it was selected as a gene of interest since this peptide has been known to regulate a variety of behavioral processes, including food intake, sexual behavior, learning, memory, reward, cognition, sleep, seizure, as well as emotion/mood-related behaviors, such as stress, fear, anxiety and depression [45]–[47]. Given that LS has many roles in common with galanin in the processes of fear, anxiety and depression, galanin may be an interesting candidate gene for regulation of LS-mediated affective processes.

Identification of Novel Genes During Lactation

The upregulation of Glul mRNA in the septum of maternal mice was noteworthy for a few reasons. The astrocytic Glul or enzyme glutamine synthetase (GS), is a key enzyme involved in the glutamate-glutamine cycle. Glutamate is metabolized mainly to glutamine by the catalyzation of Glul. Further, glutamine may also be utilized for de novo synthesis of GABA, a major inhibitory neurotransmitter. Thus, Glul modulates levels of glutamate, glutamine, and GABA in the brain [48], and influences GABAergic neurotransmission [49]. As GABA signaling regulates maternal defense [11], it may be that Glul plays a role in the emergence of offspring protection during the postpartum period. Research on the role of excitatory neurotransmitter glutamate in maternal behavior has so far received little attention, but future work evaluating glutamate and maternal care appears to be merited.

We identified increased Sstr1 transcript in the septum during lactation. The neuropeptide, somatostatin, has been linked to hormone release, cell proliferation, and many pathophysiological processes of brain disorders [50]. Somatostatin exerts its actions by interacting with specific G-protein coupled receptors. Sstr1 acts as an inhibitory autoreceptor located on somatostatin neurons and is involved in the modulation of the ultradian release of growth hormone from pituitary [51]. As somatostatin has been shown to regulate GABA receptor functions [52], and GABA signaling is involved in maternal defense, it will be interesting to know whether and how Sstr1 regulates this aspect of maternal care.

It should be noted that many genes implicated in metabolism, cell differentiation, or proliferation exhibited altered expression pattern during lactation, suggesting that septum may be an important neural site for lactation-induced metabolic and structural alterations. As can be seen in Table 1, 18 out of 116 genes were linked to cell division, death, or proliferation, while about 16% were metabolism-related genes. However, the mechanism by which these genes regulate metabolic and structural adaptations during lactation remains to be evaluated. This study provides potential candidate genes for future research to answer these questions.

Genes Previously Shown to be Differentially Expressed During Lactation Not Identified in the Present Study

Expression of neuropeptide Y, preproenkephalin, and oxytocin has been reported to be increased during lactation [22], [24], [33], [34], whereas proopiomelanocortin was down-regulated [53]. Although these changes were not observed in this study, one possible explanation for negative results is a region-specific effect as these neuropeptides were examined in other brain regions such as hypothalamus and Arc, while the septum was tested in this study. Alternatively, as discussed below, using the whole septum as a target, rather than subregions of septum, may contribute to the lack of difference. In addition to the above two possible factors, time lapse between the parturition and the gene analysis may be another contributing factor. Previous studies have shown that gene expression in specific brain regions undergoes a time-course change from parturition through lactation [24], [54], [55]. For example, preproenkephaline mRNA in the anterior arcuate nucleus of lactating versus virgin females significantly increased on postpartum day 10, but did not differ on postpartum day 3 [24], suggesting that gene expression is dynamic and alters during the postpartum period. Future study will be directed towards understanding the time-course effect on gene expression during the postpartum period.

Gene Set Enrichment Analysis

One surprising result from this study was the finding of a large number of miRNAs and transcription factors (N = 177) that were significantly linked to the upregulation of genes in maternal mice, whereas no significant links were found with downregulated maternal expression. Why so many regulators for upregulation, but not one for downregulation were found is not clear, but the finding may reflect the underlying processes of developing the maternal brain. It was also surprising to find a high number of miRNAs that were linked to upregulated gene expression. miRNAs are mostly linked with downregulation of gene expression as they can bind to mRNA transcripts, leading to degradation or suppressed translation [56]. Therefore, one explanation could be that in the process of developing the maternal brain, a number of miRNAs are themselves downregulated and this, in turn, leads to elevated expression of select genes. MiRNAs appear to play important roles in mammalian females in both developments of the placenta and mammary glands [31], [32]. Thus, an important role for miRNAs in reproduction has already been documented. To our knowledge, this study is the first to suggest that miRNAs may also play an important role in producing the maternal brain, but direct studies are needed to address this possibility.

The identification of 29 transcription factors, including Elk-1, YY1, and HSF1, that were linked to upregulated genes provides possible new insight into how the maternal brain may be produced. Although key roles for estrogen and progesterone receptors have been found in producing the maternal brain [27], [28], steroid receptors were not identified here, with the exception of estrogen related receptor 1 (Err1). One explanation is that steroid receptors alter expression of other transcription factors and the ones identified here may be key agents for some of the final changes in maternal gene expression in septum. Some of the transcription factors identified by GSEA were also found to have altered expression from the PLIER array analysis (P<0.05), such as Hsf1. It is not necessarily surprising that transcription factors or miRNAs linked by GSEA to altered expression are themselves not found to be differentially expressed because their action could have occurred prior to the day of tissue collection (e.g., during pregnancy or early postpartum).

It was also noteworthy that some conserved binding sites were found that were linked to maternally upregulated genes, although there are currently no known factors that bind to these sites. Future studies that determine whether or how these sites are acted upon will be critical for developing an understanding of the relevance of these binding sites.

Methodological Considerations and Limitations

A number of factors contribute to the formation of the maternal brain, including the experience of mating, pregnancy, parturition, lactation, and the sensory input from the pups. All these factors orchestrate many of changes in gene expression that occur during the postpartum period. This study seeks to explore the results of this constellation of experiences by examining the differences in gene expression between virgin mice, who have never been exposed to any of these experiences, and lactating mice, who have been exposed to all of them. Thus, the observed changes in gene expression in postpartum female mice do not simply reflect maternal care or lactation, but rather the integration of multiple factors. In this study, we did not discriminate between virgins at different stages of the estrus cycle although estrous stage can influence gene expression [57]. Virgin mice were either in estrus, diestrus I and diestrus II (periods of relatively consistent estradiol and progesterone levels), but not proestrus. In this way, only differences large enough to rise above any noise created by varying estrus states appear in the results. While this noise may obscure differences in gene expression, the present approach serves to highlight general differences in virgin and postpartum mice. However, subtler difference between the virgin and lactating mice may be missed. It should be noted that virgin female mice were fed a breeder diet to provide identical food to both groups. This feeding paradigm may have had a somewhat different effect on gene expression changes relative to regular rodent chow as breeder diet differentially modulates expression of genes implicated in food intake and energy balance due to its high fat content [58], [59].

One limitation of using whole septum is that we are unable to determine whether changes occur throughout septum or are specific to subregions. Future studies examining subregions can address this issue. Also, if changes in gene expression occur in opposite directions within these subregions, then the difference may be missed as false negatives. A related general issue is that microarray analysis can produce false negatives. The chips target a subset of regions along a given mRNA transcript and it is possible with different targeting that significant differences may have been found. Although different algorithms produce different false negative rates [60], [61], the approach used here, PLIER, is the current recommended algorithm as it introduces a higher sensitivity to changes in abundance for targets near background [60]. One example of a possible false negative is that we recently glutamic acid decarboxylase (GAD67) mRNA in the septum to be upregulated in lactating versus virgin mice (unpublished observation), while no difference was observed in this study. Microarray technology has become a useful tool for analyzing gene expression profiling and in identifying new genes or molecular pathways involved in multiple processes [62]–[64]. Five out of six genes identified as being differentially expressed (except Npy1r) in this study were confirmed by real-time PCR, validating our microarray analysis as reliable means to evaluate gene expression.

Conclusions

Using high density oligonucleotide microarrays, we observed a distinct gene expression profile in the septum of maternal mice. Given that the identified genes may play a role in diverse biological processes, this study provides new insights into the genes linked to the septum and brain changes that promote the emergence of maternal care during the postpartum period. Our findings suggest a potential key role for miRNAs in creating the maternal brain and could open an important new area of research. Future studies could involve directly evaluating and manipulating miRNAs to determine the relative importance of miRNAs in producing a maternal brain.

[Figure omitted. See PDF.]

Figure 2. Schematic representation of the brain regions (boxed areas) dissected for gene array analysis.

Distance from Bregma in the rostrocaudal planes is indicated. Reprinted and modified from The Allen Mouse Brain Atlas (reference atlas version 1, 2008). Abbreviations: LSc, caudal part of lateral septal nucleus; LSr, rostral part of lateral septal nucleus; MS, medial septal nucleus; VL, lateral ventricle; aco, anterior commissure, olfactory limb.

https://doi.org/10.1371/journal.pone.0038602.g002

[Figure omitted. See PDF.]

Table 4. Primers of interest and reference genes used for quantitative real-time PCR assay.

https://doi.org/10.1371/journal.pone.0038602.t004

Materials and Methods

Animals

Untested nulliparous female mice from a line of mice previously selected for high maternal defense (original stock was outbred hsd:ICR mice) (Mus domesticus) (Harlan, Madison, WI) were used. Selected mice were used because they exhibit the reliable emergence of a number of maternal characteristics, including nursing and offspring protection, and thus provide a solid platform for comparing maternal and non-maternal brains. As for any work that uses one strain, there is a limitation in knowing whether results apply to other strains. Only with comparisons across strains and species can reliable alterations with the postpartum period be identified, but here we provide a start for examining septal changes in gene expression with maternal behavior in mice. All animals were age matched (∼70 days old at time of dissection). For mating, females were housed with a breeder male (hsd:ICR strain) for 2 weeks. At the same time, virgin females were co-housed to provide similar social stimuli. When breeder males were removed, all females (pregnant and virgin) were housed individually and provided precut nesting material until dissections. Thus, virgin and postpartum females experienced similar levels of co-housing and single housing. The timing of cohousing and isolation was performed to minimize the effects of isolation-induced stress. Polypropylene cages were changed once weekly, but when pups were born (postpartum Day 0), cages were not changed for the postpartum female or the age matched virgin control for the remainder of the experiment. All animals were housed in the same room and cages of virgin and lactating females were alternated with one another on the same shelves. A 14∶10 light/dark cycle with lights on at 06∶00 h CST was used. Female mice were given ad lib access to breeder chow (Harlan) and tap water. All procedures were carried out in strict accordance with the guidelines of the National Institutes of Health Guide for the Care and Use of Laboratory Animals, and studies were approved by the University of Wisconsin Animal Care and Use Committee.

Tissue Collection and RNA Extraction

On postpartum Days 6–8, brains were removed from lactating females between 10∶00 and 12∶00 h. Brains from age-matched, virgin females were also removed on the same day and dissections were alternated between the two groups so that an equal number of dissections from each group were obtained. The animals were lightly anaesthetized with isoflurane, killed by cervical dislocation and then decapitated. Following decapitation, virgin females were examined for stage of estrous cycle using a vaginal lavage [65], [66]. Female mice in vaginal proestrus were not used because this stage is associated with mating in mice. For the microarray study, the stages of the six virgin mice were estrous, diestrous I, and diestrous II (N = 2 each). Although estrous stage can influence gene expression [57], the use of mice from different estrous states (excluding proestrus) provides a generic view of the virgin female and should provide insights into reliable differences between maternal and virgin mice. The whole brain was removed, snap frozen in isopentane on dry ice, and then stored at −80°C until sliced. Brain sections were sliced in a cryostat (Leica, CM1850, Bannockburn, IL, USA) with the first two sections at 300 microns and the last section at 200 microns and then mounted on glass slides. Target tissue was captured using a micropunch technique [67]. Microdissection of frozen brain sections was made with Brain Punch Set (Stoelting, Wood Dale, IL, USA) under a dissecting microscope. The whole septum including lateral and medial parts (Fig. 2) was collected bilaterally from Bregma 1.10 to 0.14 and pooled, so that each mouse provided one sample of the septum. Microdissections from ten animals in each group were flash frozen on dry ice and stored at −80°C until being processed for either gene array analysis or real-time PCR. Total RNA was extracted and prepared in matched pairs (virgin versus lactating) using an Aurum Total RNA Fatty and Fibrous Tissue Kit (Bio-Rad, Hercules, CA, USA) according to the manufacturer’s instructions. Following isolation, the integrity of RNA was assessed using Agilent RNA 6000 Nano Chips with Agilent Bioanalyzer 2100 (Agilent Technologies, Palo Alto, CA). The purity of RNA was tested, and the yield of RNA was determined using NanoDrop 1000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA). Purified total RNA was stored at −80°C until processed.

High-Density Oligonucleotide Array Hybridization

Six out of ten samples in each group (n = 6 per group) were randomly chosen for use in the microarray experiment. Microarray analysis was performed with the GeneChip Mouse Gene 1.0 ST array (Affymetrix, Santa Clara, CA, USA) using targets derived from total RNA isolated from septum as described above. cDNA for array hybridization was synthesized from 200 ng of total RNA using an Ambion GeneChip WT Expression Kit (Ambion, Austin, TX, USA) according to the manufacturer’s specifications. Briefly, total RNA was used to synthesize double-stranded cDNA, which was then used as a template for single-stranded cRNA synthesis. This cRNA was in turn used as a template for a second round of single-stranded cDNA synthesis, and the resultant DNA-RNA hybrids were then degraded using RNase H. Amplified cDNA was then fragmented and biotinylated using an Affymetrix WT Terminal Labeling Kit (Affymetrix, Santa Clara, CA, USA) according to the manufacturer’s specifications. Fragmented and labeled cDNA samples were hybridized with the arrays at 45°C for 16 hours. The hybridized arrays were then washed and stained according to manufacturer specifications, and arrays were scanned at 570 nm on an Affymetrix GC3000 G7 Scanner. Data were extracted and processed from scans using Affymetrix Command Console v. 3.1.1.1229. cDNA synthesis, fragmentation, labeling, array hybridization, staining, and scanning were performed by the Gene Expression Center at the University of Wisconsin–Madison.

Probeset Level Summarization and Microarray Statistical Analysis

Probeset level summarization and data normalization were performed with the PLIER algorithm with GC bin background correction using Affymetrix Power Tools v. 1.12.0. The microarray data discussed in this publication, both raw and summarized, have been deposited in NCBI’s Gene Expression Omnibus [68], and are accessible through GEO Series accession number GSE30836 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE30836). GEO reporting fulfills the requirements of the MIAME. Inferential statistics for differential expression between maternal and virgin samples were calculated using the array-specific empirical Bayesian implementation of ANOVA in the BioConductor package limma v.3.6.9 [69]. Both nominal and false discovery rate corrected p-values were calculated, and fold-change differences for each gene were calculated in Excel by dividing the limma-calculated average maternal expression by the limma-calculated average virgin expression.

Verification of Microarray Results with Quantitative Real-time PCR (qPCR)

To confirm gene expression results obtained from microarray analysis, six genes: fatty acid binding protein 7 (Fabp7), galanin (Gal), glutamate-ammonia ligase (Glul)/glutamine synthetase (GS), insulin-like growth factor binding protein 3 (Igfbp3), neuropeptide Y receptor 1 (Npy1r), and somatostatin receptor 1 (Sstr1) were evaluated using qPCR. These genes were chosen because of their possible involvement in multiple lactation-induced processes and the function of the genes is described in the Discussion. A SuperScript III First-Strand Synthesis System for RT-PCR (Invitrogen, Carlsbad, CA, USA) was used to reverse transcribe 100 ng of RNA to cDNA in an Eppendorf MasterCycler Personal PCR machine (Eppendorf, Hamburg, Germany) using poly-T 20mer primers. The cDNA was then amplified using a SsoFast EvaGreen Supermix kit (Bio-Rad, Hercules, CA, USA) in a StepOnePlus real-time PCR system (Applied Biosystems, Foster City, CA, USA). The amplification mixtures (20 µL) contained 1x SsoFast EvaGreen Supermix, 160 ng template cDNA, and 500 nM forward and reverse primers. Each sample was run in triplicate and standard amplification procedures were used. The cycling profile is as follows: an initial melting step at 95°C for 30 Sec followed by 40 cycles of a 95°C melting step for 5 sec, a 58°C annealing step except Fabp7 (at 60°C) for 20 sec, and a 72°C elongation step for 20 sec. Primers for genes of interest and reference genes (Table 4) were designed and screened for specificity using NCBI Primer-BLAST. Ywhaz and CycA were used as reference genes, as they are found to be among the most stable genes in rodent brain [70], [71] and found not to be different between lactating and virgin females in our arrays results. Following amplification, a standard curve was generated to assess the empirical PCR reaction efficiency, and a dissociation curve analysis was performed to insure specificity of PCR products. Ct values were calculated using the StepOnePlus software. The expression ratio of mRNA of genes in lactating relative to virgin (normalized to the reference genes Ywhaz and CycA) was calculated using a relative expression software tool REST 2009, which corrects for empirical PCR efficiency, allows for the use of multiple reference genes, and utilizes a randomization test of significance [72]. The N’s used for qPCR were as follows: virgin, N = 10; maternal, N = 10. Four females were added to each group from the array mice and the virgin states were estrous (N = 2), diestrous I (N = 4), and diestrous II (N = 4).

Gene Set Enrichment Analysis (GSEA)

GSEA is a technique that calculates up or down regulation of precompiled lists of genes and was performed using GSEA software v. 2.07 [30]. The specific gene sets described here came from the MSigDB C3 computationally-curated gene set of miRNA and transcription factor binding sites. For transcription factor binding sites, data from multiple species were used, including mice, humans, and rats. For miRNA, only human data sets were available, but binding sites for given miRNAs in mice and humans are highly conserved [73]. The analysis used z-transformed values for PLIER-normalized array intensities, t-tests to calculate individual gene inferential statistics, gene set permutation, and a weighted test statistic. The GSEA output provides both nominal and corrected p-values for each target.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Full list of all targets, their relative expression and P value ranking using PLIER algorithm.

https://doi.org/10.1371/journal.pone.0038602.s001

(XLS)

Table S2.

List of target sites linked significantly with upregulation of subsets of genes in postpartum mice by GSEA.

https://doi.org/10.1371/journal.pone.0038602.s002

(DOC)

Acknowledgments

The authors wish to thank Sharon Stevenson for administrative help, Wayne Davis for microarray technical assistance and Kate Skogen and Jeff Alexander for animal care.

Author Contributions

Conceived and designed the experiments: CZ SG. Performed the experiments: CZ MS TD. Analyzed the data: CZ MS SG. Contributed reagents/materials/analysis tools: SG. Wrote the paper: CZ SG.

Citation: Zhao C, Saul MC, Driessen T, Gammie SC (2012) Gene Expression Changes in the Septum: Possible Implications for MicroRNAs in Sculpting the Maternal Brain. PLoS ONE7(6): e38602. https://doi.org/10.1371/journal.pone.0038602

References

1. Montagnese CM, Poulain DA, Vincent JD, Theodosis DT (1987) Structural plasticity in the rat supraoptic nucleus during gestation, post-partum lactation and suckling-induced pseudogestation and lactation. J Endocrinol 115: 97–105.CM MontagneseDA PoulainJD VincentDT Theodosis1987Structural plasticity in the rat supraoptic nucleus during gestation, post-partum lactation and suckling-induced pseudogestation and lactation.J Endocrinol11597105

2. Russell JA, Douglas AJ, Ingram CD (2001) Brain preparations for maternity–adaptive changes in behavioral and neuroendocrine systems during pregnancy and lactation. An overview. Prog Brain Res 133: 1–38.JA RussellAJ DouglasCD Ingram2001Brain preparations for maternity–adaptive changes in behavioral and neuroendocrine systems during pregnancy and lactation. An overview.Prog Brain Res133138

3. Neumann ID (2001) Alterations in behavioral and neuroendocrine stress coping strategies in pregnant, parturient and lactating rats. Prog Brain Res 133: 143–152.ID Neumann2001Alterations in behavioral and neuroendocrine stress coping strategies in pregnant, parturient and lactating rats.Prog Brain Res133143152

4. Levy F, Gheusi G, Keller M (2011) Plasticity of the parental brain: a case for neurogenesis. J Neuroendocrinol 23: 984–993.F. LevyG. GheusiM. Keller2011Plasticity of the parental brain: a case for neurogenesis.J Neuroendocrinol23984993

5. Rosenblatt JS, Lehrman DS (1963) Maternal behavior in the laboratory rat. In: Rheingold HL, editor. Maternal behavior in mammals. New York: John Wiley and Sons. pp. 8–57.JS RosenblattDS Lehrman1963Maternal behavior in the laboratory rat.HL RheingoldMaternal behavior in mammalsNew YorkJohn Wiley and Sons857

6. Lightman SL (1992) Alterations in hypothalamic-pituitary responsiveness during lactation. Ann N Y Acad Sci 652: 340–346.SL Lightman1992Alterations in hypothalamic-pituitary responsiveness during lactation.Ann N Y Acad Sci652340346

7. Windle RJ, Wood S, Shanks N, Perks P, Conde GL, et al. (1997) Endocrine and behavioural responses to noise stress: comparison of virgin and lactating female rats during non-disrupted maternal activity. J Neuroendocrinol 9: 407–414.RJ WindleS. WoodN. ShanksP. PerksGL Conde1997Endocrine and behavioural responses to noise stress: comparison of virgin and lactating female rats during non-disrupted maternal activity.J Neuroendocrinol9407414

8. Smith MS, Grove KL (2002) Integration of the regulation of reproductive function and energy balance: lactation as a model. Front Neuroendocrinol 23: 225–256.MS SmithKL Grove2002Integration of the regulation of reproductive function and energy balance: lactation as a model.Front Neuroendocrinol23225256

9. Zhao C, Li M (2010) c-Fos identification of neuroanatomical sites associated with haloperidol and clozapine disruption of maternal behavior in the rat. Neuroscience 166: 1043–1055.C. ZhaoM. Li2010c-Fos identification of neuroanatomical sites associated with haloperidol and clozapine disruption of maternal behavior in the rat.Neuroscience16610431055

10. Fleischer S, Slotnick BM (1978) Disruption of maternal behavior in rats with lesions of the septal area. Physiol Behav 21: 189–200.S. FleischerBM Slotnick1978Disruption of maternal behavior in rats with lesions of the septal area.Physiol Behav21189200

11. Lee G, Gammie SC (2009) GABA(A) receptor signaling in the lateral septum regulates maternal aggression in mice. Behav Neurosci 123: 1169–1177.G. LeeSC Gammie2009GABA(A) receptor signaling in the lateral septum regulates maternal aggression in mice.Behav Neurosci12311691177

12. D’Anna KL, Gammie SC (2009) Activation of corticotropin-releasing factor receptor 2 in lateral septum negatively regulates maternal defense. Behav Neurosci 123: 356–368.KL D’AnnaSC Gammie2009Activation of corticotropin-releasing factor receptor 2 in lateral septum negatively regulates maternal defense.Behav Neurosci123356368

13. Risold PY, Swanson LW (1997a) Connections of the rat lateral septal complex. Brain Res Brain Res Rev 24: 115–195.PY RisoldLW Swanson1997aConnections of the rat lateral septal complex.Brain Res Brain Res Rev24115195

14. Risold PY, Swanson LW (1997b) Chemoarchitecture of the rat lateral septal nucleus. Brain Res Brain Res Rev 24: 91–113.PY RisoldLW Swanson1997bChemoarchitecture of the rat lateral septal nucleus.Brain Res Brain Res Rev2491113

15. Sheehan TP, Chambers RA, Russell DS (2004) Regulation of affect by the lateral septum: implications for neuropsychiatry. Brain Res Brain Res Rev 46: 71–117.TP SheehanRA ChambersDS Russell2004Regulation of affect by the lateral septum: implications for neuropsychiatry.Brain Res Brain Res Rev4671117

16. Jakab RL, Leranth C (1995) Septum. In: Paxinos G, editor. The Rat Nervous System. San Diego: Academic Press. pp. 405–442.RL JakabC. Leranth1995Septum.G. PaxinosThe Rat Nervous SystemSan DiegoAcademic Press405442

17. Everts HG, Koolhaas JM (1999) Differential modulation of lateral septal vasopressin receptor blockade in spatial learning, social recognition, and anxiety-related behaviors in rats. Behav Brain Res 99: 7–16.HG EvertsJM Koolhaas1999Differential modulation of lateral septal vasopressin receptor blockade in spatial learning, social recognition, and anxiety-related behaviors in rats.Behav Brain Res99716

18. Lukas M, Bredewold R, Landgraf R, Neumann ID, Veenema AH (2011) Early life stress impairs social recognition due to a blunted response of vasopressin release within the septum of adult male rats. Psychoneuroendocrinology 36: 843–853.M. LukasR. BredewoldR. LandgrafID NeumannAH Veenema2011Early life stress impairs social recognition due to a blunted response of vasopressin release within the septum of adult male rats.Psychoneuroendocrinology36843853

19. Veenema AH, Beiderbeck DI, Lukas M, Neumann ID (2010) Distinct correlations of vasopressin release within the lateral septum and the bed nucleus of the stria terminalis with the display of intermale aggression. Horm Behav 58: 273–281.AH VeenemaDI BeiderbeckM. LukasID Neumann2010Distinct correlations of vasopressin release within the lateral septum and the bed nucleus of the stria terminalis with the display of intermale aggression.Horm Behav58273281

20. Klein RL, Hirko AC, Meyers CA, Grimes JR, Muzyczka N, et al. (2000) NGF gene transfer to intrinsic basal forebrain neurons increases cholinergic cell size and protects from age-related, spatial memory deficits in middle-aged rats. Brain Res 875: 144–151.RL KleinAC HirkoCA MeyersJR GrimesN. Muzyczka2000NGF gene transfer to intrinsic basal forebrain neurons increases cholinergic cell size and protects from age-related, spatial memory deficits in middle-aged rats.Brain Res875144151

21. Rudick CN, Gibbs RB, Woolley CS (2003) A role for the basal forebrain cholinergic system in estrogen-induced disinhibition of hippocampal pyramidal cells. J Neurosci 23: 4479–4490.CN RudickRB GibbsCS Woolley2003A role for the basal forebrain cholinergic system in estrogen-induced disinhibition of hippocampal pyramidal cells.J Neurosci2344794490

22. Smith MS (1993) Lactation alters neuropeptide-Y and proopiomelanocortin gene expression in the arcuate nucleus of the rat. Endocrinology 133: 1258–1265.MS Smith1993Lactation alters neuropeptide-Y and proopiomelanocortin gene expression in the arcuate nucleus of the rat.Endocrinology13312581265

23. Chen P, Li C, Haskell-Luevano C, Cone RD, Smith MS (1999) Altered expression of agouti-related protein and its colocalization with neuropeptide Y in the arcuate nucleus of the hypothalamus during lactation. Endocrinology 140: 2645–2650.P. ChenC. LiC. Haskell-LuevanoRD ConeMS Smith1999Altered expression of agouti-related protein and its colocalization with neuropeptide Y in the arcuate nucleus of the hypothalamus during lactation.Endocrinology14026452650

24. Ottinger MA, Rosewell KL, Weiland NG, Margaretten KT, Wise PM (1995) Effect of lactation on hypothalamic preproenkephalin gene expression. J Neuroendocrinol 7: 341–346.MA OttingerKL RosewellNG WeilandKT MargarettenPM Wise1995Effect of lactation on hypothalamic preproenkephalin gene expression.J Neuroendocrinol7341346

25. Broad KD, Kendrick KM, Sirinathsinghji DJ, Keverne EB (1993) Changes in oxytocin immunoreactivity and mRNA expression in the sheep brain during pregnancy, parturition and lactation and in response to oestrogen and progesterone. J Neuroendocrinol 5: 435–444.KD BroadKM KendrickDJ SirinathsinghjiEB Keverne1993Changes in oxytocin immunoreactivity and mRNA expression in the sheep brain during pregnancy, parturition and lactation and in response to oestrogen and progesterone.J Neuroendocrinol5435444

26. Kuroda KO, Tachikawa K, Yoshida S, Tsuneoka Y, Numan M (2011) Neuromolecular basis of parental behavior in laboratory mice and rats: with special emphasis on technical issues of using mouse genetics. Prog Neuropsychopharmacol Biol Psychiatry 35: 1205–1231.KO KurodaK. TachikawaS. YoshidaY. TsuneokaM. Numan2011Neuromolecular basis of parental behavior in laboratory mice and rats: with special emphasis on technical issues of using mouse genetics.Prog Neuropsychopharmacol Biol Psychiatry3512051231

27. Giordano AL, Siegel HI, Rosenblatt JS (1989) Nuclear estrogen receptor binding in the preoptic area and hypothalamus of pregnancy-terminated rats: correlation with the onset of maternal behavior. Neuroendocrinology 50: 248–258.AL GiordanoHI SiegelJS Rosenblatt1989Nuclear estrogen receptor binding in the preoptic area and hypothalamus of pregnancy-terminated rats: correlation with the onset of maternal behavior.Neuroendocrinology50248258

28. Numan M, Roach JK, del Cerro MC, Guillamon A, Segovia S, et al. (1999) Expression of intracellular progesterone receptors in rat brain during different reproductive states, and involvement in maternal behavior. Brain Res 830: 358–371.M. NumanJK RoachMC del CerroA. GuillamonS. Segovia1999Expression of intracellular progesterone receptors in rat brain during different reproductive states, and involvement in maternal behavior.Brain Res830358371

29. Kinsley CH, Amory-Meyer E (2011) Why the maternal brain? J Neuroendocrinol 23: 974–983.CH KinsleyE. Amory-Meyer2011Why the maternal brain?J Neuroendocrinol23974983

30. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, et al. (2005) Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A 102: 15545–15550.A. SubramanianP. TamayoVK MoothaS. MukherjeeBL Ebert2005Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.Proc Natl Acad Sci U S A1021554515550

31. Prieto DM, Markert UR (2011) MicroRNAs in pregnancy. J Reprod Immunol 88: 106–111.DM PrietoUR Markert2011MicroRNAs in pregnancy.J Reprod Immunol88106111

32. Silveri L, Tilly G, Vilotte JL, Le Provost F (2006) MicroRNA involvement in mammary gland development and breast cancer. Reprod Nutr Dev 46: 549–556.L. SilveriG. TillyJL VilotteF. Le Provost2006MicroRNA involvement in mammary gland development and breast cancer.Reprod Nutr Dev46549556

33. Xiao XQ, Grove KL, Lau SY, McWeeney S, Smith MS (2005) Deoxyribonucleic acid microarray analysis of gene expression pattern in the arcuate nucleus/ventromedial nucleus of hypothalamus during lactation. Endocrinology 146: 4391–4398.XQ XiaoKL GroveSY LauS. McWeeneyMS Smith2005Deoxyribonucleic acid microarray analysis of gene expression pattern in the arcuate nucleus/ventromedial nucleus of hypothalamus during lactation.Endocrinology14643914398

34. Gammie SC, Hasen NS, Awad TA, Auger AP, Jessen HM, et al. (2005) Gene array profiling of large hypothalamic CNS regions in lactating and randomly cycling virgin mice. Brain Res Mol Brain Res 139: 201–211.SC GammieNS HasenTA AwadAP AugerHM Jessen2005Gene array profiling of large hypothalamic CNS regions in lactating and randomly cycling virgin mice.Brain Res Mol Brain Res139201211

35. Clemmons DR, Underwood LE (1991) Nutritional regulation of IGF-I and IGF binding proteins. Annu Rev Nutr 11: 393–412.DR ClemmonsLE Underwood1991Nutritional regulation of IGF-I and IGF binding proteins.Annu Rev Nutr11393412

36. Lauterio TJ, Marson L, Daughaday WH, Baile CA (1987) Evidence for the role of insulin-like growth factor II (IGF-II) in the control of food intake. Physiol Behav 40: 755–758.TJ LauterioL. MarsonWH DaughadayCA Baile1987Evidence for the role of insulin-like growth factor II (IGF-II) in the control of food intake.Physiol Behav40755758

37. Butler ST, Marr AL, Pelton SH, Radcliff RP, Lucy MC, et al. (2003) Insulin restores GH responsiveness during lactation-induced negative energy balance in dairy cattle: effects on expression of IGF-I and GH receptor 1A. J Endocrinol 176: 205–217.ST ButlerAL MarrSH PeltonRP RadcliffMC Lucy2003Insulin restores GH responsiveness during lactation-induced negative energy balance in dairy cattle: effects on expression of IGF-I and GH receptor 1A.J Endocrinol176205217

38. Kalluri HS, Dempsey RJ (2011) IGFBP-3 inhibits the proliferation of neural progenitor cells. Neurochem Res 36: 406–411.HS KalluriRJ Dempsey2011IGFBP-3 inhibits the proliferation of neural progenitor cells.Neurochem Res36406411

39. Gluckman P, Klempt N, Guan J, Mallard C, Sirimanne E, et al. (1992) A role for IGF-1 in the rescue of CNS neurons following hypoxic-ischemic injury. Biochem Biophys Res Commun 182: 593–599.P. GluckmanN. KlemptJ. GuanC. MallardE. Sirimanne1992A role for IGF-1 in the rescue of CNS neurons following hypoxic-ischemic injury.Biochem Biophys Res Commun182593599

40. Mills JL, Hediger ML, Molloy CA, Chrousos GP, Manning-Courtney P, et al. (2007) Elevated levels of growth-related hormones in autism and autism spectrum disorder. Clin Endocrinol (Oxf) 67: 230–237.JL MillsML HedigerCA MolloyGP ChrousosP. Manning-Courtney2007Elevated levels of growth-related hormones in autism and autism spectrum disorder.Clin Endocrinol (Oxf)67230237

41. Arai Y, Funatsu N, Numayama-Tsuruta K, Nomura T, Nakamura S, et al. (2005) Role of Fabp7, a downstream gene of Pax6, in the maintenance of neuroepithelial cells during early embryonic development of the rat cortex. J Neurosci 25: 9752–9761.Y. AraiN. FunatsuK. Numayama-TsurutaT. NomuraS. Nakamura2005Role of Fabp7, a downstream gene of Pax6, in the maintenance of neuroepithelial cells during early embryonic development of the rat cortex.J Neurosci2597529761

42. Owada Y, Abdelwahab SA, Kitanaka N, Sakagami H, Takano H, et al. (2006) Altered emotional behavioral responses in mice lacking brain-type fatty acid-binding protein gene. Eur J Neurosci 24: 175–187.Y. OwadaSA AbdelwahabN. KitanakaH. SakagamiH. Takano2006Altered emotional behavioral responses in mice lacking brain-type fatty acid-binding protein gene.Eur J Neurosci24175187

43. Gerstner JR, Vanderheyden WM, Shaw PJ, Landry CF, Yin JC (2011) Fatty-acid binding proteins modulate sleep and enhance long-term memory consolidation in Drosophila. PLoS One 6: e15890.JR GerstnerWM VanderheydenPJ ShawCF LandryJC Yin2011Fatty-acid binding proteins modulate sleep and enhance long-term memory consolidation in Drosophila.PLoS One6e15890

44. Gross L (2007) A candidate gene for a biological marker of schizophrenia in mice. PLoS Biol 5: e320.L. Gross2007A candidate gene for a biological marker of schizophrenia in mice.PLoS Biol5e320

45. Karlsson RM, Holmes A (2006) Galanin as a modulator of anxiety and depression and a therapeutic target for affective disease. Amino Acids 31: 231–239.RM KarlssonA. Holmes2006Galanin as a modulator of anxiety and depression and a therapeutic target for affective disease.Amino Acids31231239

46. Hokfelt T, Broberger C, Diez M, Xu ZQ, Shi T, et al. (1999) Galanin and NPY, two peptides with multiple putative roles in the nervous system. Horm Metab Res 31: 330–334.T. HokfeltC. BrobergerM. DiezZQ XuT. Shi1999Galanin and NPY, two peptides with multiple putative roles in the nervous system.Horm Metab Res31330334

47. Wrenn CC, Crawley JN (2001) Pharmacological evidence supporting a role for galanin in cognition and affect. Prog Neuropsychopharmacol Biol Psychiatry 25: 283–299.CC WrennJN Crawley2001Pharmacological evidence supporting a role for galanin in cognition and affect.Prog Neuropsychopharmacol Biol Psychiatry25283299

48. Fonnum F, Paulsen RE (1990) Comparison of transmitter amino acid levels in rat globus pallidus and neostriatum during hypoglycemia or after treatment with methionine sulfoximine or gamma-vinyl gamma-aminobutyric acid. J Neurochem 54: 1253–1257.F. FonnumRE Paulsen1990Comparison of transmitter amino acid levels in rat globus pallidus and neostriatum during hypoglycemia or after treatment with methionine sulfoximine or gamma-vinyl gamma-aminobutyric acid.J Neurochem5412531257

49. Cremer CM, Bidmon HJ, Gorg B, Palomero-Gallagher N, Escobar JL, et al. (2010) Inhibition of glutamate/glutamine cycle in vivo results in decreased benzodiazepine binding and differentially regulated GABAergic subunit expression in the rat brain. Epilepsia 51: 1446–1455.CM CremerHJ BidmonB. GorgN. Palomero-GallagherJL Escobar2010Inhibition of glutamate/glutamine cycle in vivo results in decreased benzodiazepine binding and differentially regulated GABAergic subunit expression in the rat brain.Epilepsia5114461455

50. Olias G, Viollet C, Kusserow H, Epelbaum J, Meyerhof W (2004) Regulation and function of somatostatin receptors. J Neurochem 89: 1057–1091.G. OliasC. ViolletH. KusserowJ. EpelbaumW. Meyerhof2004Regulation and function of somatostatin receptors.J Neurochem8910571091

51. Stroh T, van Schouwenburg MR, Beaudet A, Tannenbaum GS (2009) Subcellular dynamics of somatostatin receptor subtype 1 in the rat arcuate nucleus: receptor localization and synaptic connectivity vary in parallel with the ultradian rhythm of growth hormone secretion. J Neurosci 29: 8198–8205.T. StrohMR van SchouwenburgA. BeaudetGS Tannenbaum2009Subcellular dynamics of somatostatin receptor subtype 1 in the rat arcuate nucleus: receptor localization and synaptic connectivity vary in parallel with the ultradian rhythm of growth hormone secretion.J Neurosci2981988205

52. Cammalleri M, Cervia D, Langenegger D, Liu Y, Dal Monte M, et al. (2004) Somatostatin receptors differentially affect spontaneous epileptiform activity in mouse hippocampal slices. Eur J Neurosci 20: 2711–2721.M. CammalleriD. CerviaD. LangeneggerY. LiuM. Dal Monte2004Somatostatin receptors differentially affect spontaneous epileptiform activity in mouse hippocampal slices.Eur J Neurosci2027112721

53. Mann PE, Rubin BS, Bridges RS (1997) Differential proopiomelanocortin gene expression in the medial basal hypothalamus of rats during pregnancy and lactation. Brain Res Mol Brain Res 46: 9–16.PE MannBS RubinRS Bridges1997Differential proopiomelanocortin gene expression in the medial basal hypothalamus of rats during pregnancy and lactation.Brain Res Mol Brain Res46916

54. Mann PE, Bridges RS (2002) Prolactin receptor gene expression in the forebrain of pregnant and lactating rats. Brain Res Mol Brain Res 105: 136–145.PE MannRS Bridges2002Prolactin receptor gene expression in the forebrain of pregnant and lactating rats.Brain Res Mol Brain Res105136145

55. da Costa AP, Wood S, Ingram CD, Lightman SL (1996) Region-specific reduction in stress-induced c-fos mRNA expression during pregnancy and lactation. Brain Res 742: 177–184.AP da CostaS. WoodCD IngramSL Lightman1996Region-specific reduction in stress-induced c-fos mRNA expression during pregnancy and lactation.Brain Res742177184

56. Shenouda SK, Alahari SK (2009) MicroRNA function in cancer: oncogene or a tumor suppressor? Cancer Metastasis Rev 28: 369–378.SK ShenoudaSK Alahari2009MicroRNA function in cancer: oncogene or a tumor suppressor?Cancer Metastasis Rev28369378

57. Romano GJ, Mobbs CV, Howells RD, Pfaff DW (1989) Estrogen regulation of proenkephalin gene expression in the ventromedial hypothalamus of the rat: temporal qualities and synergism with progesterone. Brain Res Mol Brain Res 5: 51–58.GJ RomanoCV MobbsRD HowellsDW Pfaff1989Estrogen regulation of proenkephalin gene expression in the ventromedial hypothalamus of the rat: temporal qualities and synergism with progesterone.Brain Res Mol Brain Res55158

58. Mercer JG, Lawrence CB, Atkinson T (1996) Regulation of galanin gene expression in the hypothalamic paraventricular nucleus of the obese Zucker rat by manipulation of dietary macronutrients. Brain Res Mol Brain Res 43: 202–208.JG MercerCB LawrenceT. Atkinson1996Regulation of galanin gene expression in the hypothalamic paraventricular nucleus of the obese Zucker rat by manipulation of dietary macronutrients.Brain Res Mol Brain Res43202208

59. Schaffhauser AO, Madiehe AM, Braymer HD, Bray GA, York DA (2002) Effects of a high-fat diet and strain on hypothalamic gene expression in rats. Obes Res 10: 1188–1196.AO SchaffhauserAM MadieheHD BraymerGA BrayDA York2002Effects of a high-fat diet and strain on hypothalamic gene expression in rats.Obes Res1011881196

60. Seo J, Hoffman EP (2006) Probe set algorithms: is there a rational best bet? BMC Bioinformatics 7: 395.J. SeoEP Hoffman2006Probe set algorithms: is there a rational best bet?BMC Bioinformatics7395

61. Rosati B, Grau F, Kuehler A, Rodriguez S, McKinnon D (2004) Comparison of different probe-level analysis techniques for oligonucleotide microarrays. Biotechniques 36: 316–322.B. RosatiF. GrauA. KuehlerS. RodriguezD. McKinnon2004Comparison of different probe-level analysis techniques for oligonucleotide microarrays.Biotechniques36316322

62. Uriguen L, Arteta D, Diez-Alarcia R, Ferrer-Alcon M, Diaz A, et al. (2008) Gene expression patterns in brain cortex of three different animal models of depression. Genes Brain Behav 7: 649–658.L. UriguenD. ArtetaR. Diez-AlarciaM. Ferrer-AlconA. Diaz2008Gene expression patterns in brain cortex of three different animal models of depression.Genes Brain Behav7649658

63. Renn SC, Aubin-Horth N, Hofmann HA (2004) Biologically meaningful expression profiling across species using heterologous hybridization to a cDNA microarray. BMC Genomics 5: 42.SC RennN. Aubin-HorthHA Hofmann2004Biologically meaningful expression profiling across species using heterologous hybridization to a cDNA microarray.BMC Genomics542

64. Schumer M, Krishnakant K, Renn SC (2011) Comparative gene expression profiles for highly similar aggressive phenotypes in male and female cichlid fishes (Julidochromis). J Exp Biol 214: 3269–3278.M. SchumerK. KrishnakantSC Renn2011Comparative gene expression profiles for highly similar aggressive phenotypes in male and female cichlid fishes (Julidochromis).J Exp Biol21432693278

65. Drazen DL, Klein SL, Burnett AL, Wallach EE, Crone JK, et al. (1999) Reproductive function in female mice lacking the gene for endothelial nitric oxide synthase. Nitric Oxide 3: 366–374.DL DrazenSL KleinAL BurnettEE WallachJK Crone1999Reproductive function in female mice lacking the gene for endothelial nitric oxide synthase.Nitric Oxide3366374

66. Marcondes FK, Bianchi FJ, Tanno AP (2002) Determination of the estrous cycle phases of rats: some helpful considerations. Brazilian Journal of Biology 62: 609–614.FK MarcondesFJ BianchiAP Tanno2002Determination of the estrous cycle phases of rats: some helpful considerations.Brazilian Journal of Biology62609614

67. Makino S, Gold PW, Schulkin J (1994) Effects of corticosterone on CRH mRNA and content in the bed nucleus of the stria terminalis; comparison with the effects in the central nucleus of the amygdala and the paraventricular nucleus of the hypothalamus. Brain Res 657: 141–149.S. MakinoPW GoldJ. Schulkin1994Effects of corticosterone on CRH mRNA and content in the bed nucleus of the stria terminalis; comparison with the effects in the central nucleus of the amygdala and the paraventricular nucleus of the hypothalamus.Brain Res657141149

68. Edgar R, Domrachev M, Lash AE (2002) Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res 30: 207–210.R. EdgarM. DomrachevAE Lash2002Gene Expression Omnibus: NCBI gene expression and hybridization array data repository.Nucleic Acids Res30207210

69. Smyth GK (2004) Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol 3: Article3.GK Smyth2004Linear models and empirical bayes methods for assessing differential expression in microarray experiments.Stat Appl Genet Mol Biol3Article3

70. Nelissen K, Smeets K, Mulder M, Hendriks JJ, Ameloot M (2010) Selection of reference genes for gene expression studies in rat oligodendrocytes using quantitative real time PCR. J Neurosci Methods 187: 78–83.K. NelissenK. SmeetsM. MulderJJ HendriksM. Ameloot2010Selection of reference genes for gene expression studies in rat oligodendrocytes using quantitative real time PCR.J Neurosci Methods1877883

71. Bonefeld BE, Elfving B, Wegener G (2008) Reference genes for normalization: a study of rat brain tissue. Synapse 62: 302–309.BE BonefeldB. ElfvingG. Wegener2008Reference genes for normalization: a study of rat brain tissue.Synapse62302309

72. Pfaffl MW, Horgan GW, Dempfle L (2002) Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucleic Acids Res 30: e36.MW PfafflGW HorganL. Dempfle2002Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR.Nucleic Acids Res30e36

73. Moss EG, Tang L (2003) Conservation of the heterochronic regulator Lin-28, its developmental expression and microRNA complementary sites. Dev Biol 258: 432–442.EG MossL. Tang2003Conservation of the heterochronic regulator Lin-28, its developmental expression and microRNA complementary sites.Dev Biol258432442

74. Caldwell KK, Sheema S, Paz RD, Samudio-Ruiz SL, Laughlin MH, et al. (2008) Fetal alcohol spectrum disorder-associated depression: evidence for reductions in the levels of brain-derived neurotrophic factor in a mouse model. Pharmacol Biochem Behav 90: 614–624.KK CaldwellS. SheemaRD PazSL Samudio-RuizMH Laughlin2008Fetal alcohol spectrum disorder-associated depression: evidence for reductions in the levels of brain-derived neurotrophic factor in a mouse model.Pharmacol Biochem Behav90614624