| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Division of Neuroscience (X.Q.X., K.L.G., M.S.S.), Oregon National Primate Research Center; Biostatistics & Bioinformatics Shared Resource (S.Y.L., S.M.); and Department of Physiology and Pharmacology (M.S.S.), Oregon Health & Science University, Beaverton, Oregon 97006
Address all correspondence and requests for reprints to: Dr. M. Susan Smith, Division of Neuroscience, Oregon National Primate Research Center, Oregon Health & Science University, 505 Northwest 185th Avenue, Beaverton, Oregon 97006. E-mail: smithsu{at}ohsu.edu.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
| Materials and Methods |
|---|
|
|
|---|
Microdissection of ARH/ventromedial nucleus of the hypothalamus (VMH), RNA isolation and purification
Tissues for microarray or real-time PCR were collected by rapid decapitation. Brains were quickly removed, frozen on dry ice, and stored at 80 C until the microdissection of hypothalamic nuclei using a micropunch technique. Frozen coronal brain sections (300 µm) were cut in a cryostat at 9 C, then thaw-mounted onto glass slides and refrozen. Using a 1-mm-diameter circular punch (Biomedical Research Instruments, Inc., Malden, MA), tissues were obtained from the sections with coordinates relative to bregma between 2.0 and 4.0, according to the rat atlas. The micropunches from these areas include all of the ARH and the ventral part of the VMH. For microarray, four ARH/VMH pools (each pool contained the punched tissues from three animals) were prepared for the Cont and Lac+8 groups. An additional five ARH/VMH pools (each pool contained the punched tissues from two animals) were prepared from Cont, Lac+8, and Lac+0 groups for real-time PCR. The pools were homogenized in TRIzol reagent (Invitrogen, Carlsbad, CA), and total RNA was isolated according to the manufacturers specifications, as previously described (12). RNA was further purified using an affinity resin column (RNeasy Mini Kit; QIAGEN, Valencia, CA). The quality and concentration of the RNA were determined by measuring the absorbance at 260 and 280 nm, and RNA integrity was confirmed by an Agilent 2100 Bioanalyzer (Agilent Technologies, Inc., Palo Alto, CA).
Microarray procedures
Microarray assays are performed in the Affymetrix Microarray Core (Affymetrix, Inc., Santa Clara, CA), a unit of the OHSU Gene Microarray Shared Resource, using procedures adapted from the Affymetrix Gene-Chip Expression Analysis Technical Manual, revision 3 (http://www.affymetrix.com/support/technical/manual/expression_manual.affx). Briefly, mRNA is amplified and labeled from 4 µg total RNA, and mRNA is converted to double-stranded cDNA using Superscript Reverse Transcriptase (Invitrogen) and an oligo-dT primer linked to a T7 RNA polymerase binding site sequence (Integrated DNA Technologies, Inc., Coralville, IA). The amplified, labeled cRNA (the target) is produced in an in vitro transcription reaction using T7 RNA polymerase and biotin-CTP (GeneChip IVT Labeling Kit, Affymetrix). After removal of free nucleotides, target yield is measured by UV260 absorbance. Using microfluidics separation technology, 200 ng of each labeled target is analyzed alongside a control cRNA target on an RNA NanoChip processed on the 2100 Bioanalyzer (Agilent). Target quality is based on cRNA yield and size distribution. Samples that fail yield and size distribution thresholds are discarded or relabeled. Borderline or atypical targets are applied to test arrays for supplemental quality assessment before hybridization to a GeneChip genome array. Labeled target is fragmented at 95 C in the presence of high magnesium concentration. The fragmented material is combined with biotinylated hybridization control oligomer and biotinylated control cRNAs for the bacterial genes BioB, BioC, BioD, and CreX (Affymetrix) in hybridization buffer. Ten micrograms of target are hybridized with the RAE 230 A GeneChip (http://www.affymetrix.com/products/arrays/specific/rat230.affx) overnight, followed by washing, staining with streptavidin-phycoerythrin (Molecular Probes, Inc., Eugene, OR), signal amplification with biotinylated antistreptavidin antibody (Vector Laboratories, Inc., Burlingame, CA), and a final staining step on the Fluidics Station 450 (Affymetrix). The distribution of fluorescent material on the processed array is determined using the GeneArray laser scanner (Affymetrix). Image inspection is performed manually immediately after each scan.
Data measurement
Scans were performed with an Affymetrix GeneChip scanner, and the expression value for each gene was calculated using the Affymetrix Microarray Analysis Suite 5.0 (MAS5.0), computing the expression intensities in signal units defined by the software. After the initial analysis, the absolute analyses were rerun using global scaling to an average target intensity of 300. The scaling allows for the direct comparison of hybridization values from the different targets analyzed. For each analysis, scaled or unscaled, the parameters
1 and
2 were set to 0.1 and 0.15, respectively. These parameters set the point at which a probe set is called present, marginal, or undetectable. This call is based on the detection P value of the probe set as determined by the software.
Affymetrix comparison algorithms incorporated in the MAS5.0 software were used to derive metrics that identify differences between two-probe arrays for every probe set. We performed the following procedures to identify genes affected in each condition. First, a pair-wise comparison analysis was performed among the four hybridizations prepared from Lac+8 (used as experimental arrays) and the four Cont (used as baseline arrays). Among the 16 possible comparisons between four Cont pools and four Lac+8 pools, eight comparisons were randomly selected, which provides a very high level of representation and avoids bias. For the high-stringency analysis, all eight of the eight randomly selected comparisons (100%) had to show differential expression during lactation using a signed-rank test by the MAS software (P < 0.045). As an additional stringency filter, genes with expression ratios lower than ±1.3 were eliminated. However, due to the inherent variability of microarray analysis, it was likely that some important genes would not meet these criteria. Therefore, we performed a moderate-stringency analysis in which differential expression had to be observed in six of eight comparisons (75%) between Cont pools and Lac+8 pools. Because this decreased stringency increased the likelihood of identifying false positives, genes with expression ratios lower than ±1.4 were eliminated.
We used NetAffx (http://www.affymetrix.com), an interactive database, to retrieve the probe information, functional annotation, map location, and other related information for each gene. Genes were manually grouped into functional clusters based on their GeneOntology annotations (primarily GO: biological process) and PubMed literature.
Real-time PCR
Methods used for real-time PCR are similar to those previously described (12, 13). RNA samples were prepared for real-time PCR by random-primed reverse transcription reaction using random hexamer primers (Promega Corp., Madison, WI) and 1 µg RNA. The sequences of primers and probes used are summarized in Table 1
. The reverse transcription reaction was then diluted 1:50 for PCR analysis. The amplification was performed as follows: 2 min at 50 C, 10 min at 95 C, then 40 cycles each at 95 C for 15 sec and 60 C for 60 sec in the ABI/Prism 7900 Sequences Detector System (Applied Biosystems, Foster City, CA). After PCR was completed, baseline and threshold values were set to optimize the amplification plot, and the data were exported to an Excel spreadsheet. Standard curves were drawn on the basis of the log of the input RNA vs. the critical threshold cycle, which is the cycle in which the fluorescence of the sample was greater than the threshold of baseline fluorescence. These standard curves allowed for the critical threshold values to be converted to relative RNA concentrations for each sample. 18S RNA amplifications were conducted with the Pre-Developed TaqMan Assay Reagent (Applied Biosystems). The primers and probe for G protein-coupled receptor 88 (GPR88) were provided by Custom TaqMan Gene Expression from Applied Biosystems (Rn00581712_m1), and other primers and probes were designed using the Primer Express software from Applied Biosystems (Table 1
).
|
| Results |
|---|
|
|
|---|
= 0.001), given that eight of the 12 genes were identified by both methods. Expectedly, the majority of affected genes encoded proteins involved in mechanisms of neurotransmission, signal transduction, growth modulation, transporters, and structure remodeling.
|
New genes identified as being differentially expressed during lactation
Neurotransmission-related genes.
In ARH/VMH punches, oxytocin transcript increased about 103% during lactation (Table 2
). Oxytocin has been classically identified in magnocellular neurons of the PVH and supraoptic nucleus (SON) (18). Although it is highly unlikely that these micropunches included any portion of the PVH or the main body of the SON, it is possible that the accessory SON was included, thus explaining the presence of oxytocin. The increase in oxytocin fits with all we know about changes in this system during lactation (18). Galanin, an important peptide neuromodulator involved in stimulating appetite and neuroendocrine regulation (19), was reduced by 36% (Table 2
). Finally, endothelin-converting enzyme-like 1 (ECEL1), an enzyme of the M13 family of endopeptidases and previously known as XCE (20), was decreased by 50% during lactation (Table 2
).
Growth factor-related genes.
The most surprising finding of this study was the large increase in IGF binding protein 3 (IGFBP3) (Table 2
) during lactation, which is an important modulator of IGF action (21). VGF, a secreted polypeptide that is synthesized by neurons and is abundant in the hypothalamus (22), was reduced by 70% during lactation (Table 2
).
Signaling-related genes.
CR16, an SH3 domain binding protein, which is encoded by a gene previously cloned as a glucocorticoid-regulated mRNA from a rat cDNA library (23), was increased by almost 80% during lactation; and similarly, serum/glucocorticoid-induced kinase (SGK) (24) was increased by 60% (Table 2
).
Other genes.
The glial-specific structural protein, glial fibrillary acidic protein (GFAP), was down-regulated by 40% during lactation (Table 2
). The metabolic transporter, Slc16a1 (monocarboxylate transporter 1, MCT1), which is mainly involved in the transport of lactate across cell membranes, showed a moderate reduction (35%, Table 2
).
Additional genes identified as being differentially expressed during lactation: moderate-stringency analysis
Because the high-stringency analysis likely overlooked genes that were truly differentially expressed, the stringency was reduced by requiring significant changes in only six of eight random comparisons between Cont and Lac+8 pools, instead of eight of eight comparisons as described above. This analysis resulted in identifying an additional 10 genes (Table 3
) and three transcribed sequences (expressed sequence tags). Two of these genes, dynorphin and cocaine- and amphetamine-regulated transcript (CART), are known as satiety agents (25, 26). Dynorphin, one of endogenous opioids expressed in the ARH, was down-regulated by 60%, consistent with previous reports (25), whereas CART was reduced by 45% during lactation. In addition, remarkable reductions in endothelin receptor type B (EDNRB), GHRH, and fatty acid binding protein 7 and increases in sulfotransferase family 1A, connective tissue growth factor (CTGF), and GPR88 were observed during lactation (Table 3
).
|
|
| Discussion |
|---|
|
|
|---|
Methodology and validation of microarray analysis
The micropunches of the ARH/VMH used in this study included the whole rostro-caudal extent of the ARH and portions of the VMH, containing a very heterogenous population of neurons. Therefore, regional specific changes in gene expression may not have been detected in these pools. In addition, these studies were conducted in OVX females to minimize the effects of differences in ovarian sex steroid hormone levels on ARH gene expression. For these studies, we also pooled samples from three animals to obtain sufficient mRNA for microarray analysis. Thus, the observed changes are more representative of the general population. In addition, it should be noted that, in the analysis of the microarray data, we used a very conservative approach to detect the significant changes. The filter to exclude genes that were significantly variable among replicates in the same experimental group likely eliminates some genes that are truly differentially expressed, at the same time providing a robust stringency that limits the selection of false positives. For example, both neurokinin B and signal transducer and activator of transcription 3 tended to be down-regulated during lactation. However, they were not selected by either high-stringency or moderate-stringency analysis due to the variability specifically in the lactation group. When analyzed by real-time PCR, neurokinin B was reduced by 50% (n = 5; P < 0.05), whereas signal transducer and activator of transcription 33 remained unaltered during lactation (data not shown). Another limitation of the current study is the selection of the gene chip, which lacked several interesting genes involved in the regulation of energy balance, such as POMC, AGRP, SOCS3 (suppressor of cytokine signaling-3), and MC3 (melanocortin-3) receptor. An important validation of our analysis was the identification of changes in expression of genes that previously had been shown to be altered during lactation. The increase in ARH NPY during lactation has been firmly established (7, 9, 14) and is thought to be an important component in the hyperphagia and suppression of gonadotropin secretion. The increase in ARH enkephalin gene expression during lactation has been proposed to be involved in the hyperprolactinemia and the suppression of gonadotropin secretion (15). TH, the rate-limiting enzyme for dopamine synthesis and a marker of tuberoinfundibular dopaminergic neurons in the ARH, was down-regulated over 100% in lactation. The decrease in dopamine activity has been shown to be associated with suckling-induced prolactin secretion (16, 17). Because dynorphin, along with other opioids in the ARH, induces satiety, its reduction during lactation is consistent with the hyperphagia and negative energy balance (25).
Identification of new genes in the ARH that are differentially expressed and may contribute to the hyperphagia during lactation
During lactation, 2- to 4-fold increases in food and water intake are induced to meet the dramatic increase in the nutrient requirements of the mother for milk production (2). This is likely accomplished through the changing pattern of orexigenic and anorexigenic signals (4, 5, 6, 25, 29). In the current study, we confirm the changes in some of these compounds using microarray, as well as identify several new genes that may be involved in the hyperphagia.
The reduction in CART in the ARH during lactation is consistent with changes observed in response to food deprivation and its role as a potent anorexic agent (26, 30). CART expression in the ARH is limited to POMC neurons; POMC expression is also down-regulated during lactation (4). Decreased CART may promote the chronic increase in food intake needed to counter the vast efflux of nutrients from the mammary gland for milk production during lactation.
The reduction in galanin observed in the ARH during lactation does not fit with its role as a stimulator of food intake. Administration of galanin into the third ventricle induces a moderate increase in food intake (19). However, it is difficult to discern the effects of endogenous galanin expression on feeding behavior, because galanin-overexpressing and knockout mice both have normal body weights and feeding patterns (31). It is likely that the reduced galanin expression in ARH/VMH during lactation is more closely related to alterations in other neuroendocrine functions instead of food intake (32). Galanin is coexpressed in GHRH neurons in the ARH (33), and its decrease during lactation correlates with the decrease in GHRH expression.
The glucocorticoid, corticosterone, is a factor mediating the regulation of food intake and body adiposity. Interestingly, in rodents, most forms of genetic obesity are associated with increased glucocorticoid levels, whereas a lack of glucocorticoids is linked to hypophagia and reduced body weight (34). It has been reported that lactating animals have elevated levels of corticosterone and ACTH (35). In addition to the glucocorticoids and mineralocorticoids themselves, our study shows that some of their responsive elements, SGK (36) and CR16 (23), were also significantly increased in the ARH/VMH, which may contribute to the hyperphagia.
The role of VGF in the regulation of food intake may be complex. In the normal mouse, VGF is coexpressed in POMC neurons, but expression is increased in NPY neurons in response to fasting (22). VGF knockout mice are lean, hypermetabolic, and infertile (37), suggesting a critical role in the regulation of energy balance and an overall involvement in stimulating food intake. This role does not fit with VGF being inhibited in the ARH during lactation. However, the VGF gene is expressed extensively throughout the CNS, so it is difficult to determine where it acts to modulate energy homeostasis.
Monocarboxylates, such as lactate and pyruvate, play a central role in cellular metabolism. Essential to these roles is their rapid transport across the plasma membrane, which is regulated by recently identified MCTs. MCT1 is strongly expressed by neurons and astrocytes, which may have important implications for delivery of fuel to the brain (38). In the present study, we found a moderate inhibition of MCT1 expression in the ARH/VMH during lactation. This may be a very important adaptation to support the hyperphagia associated with lactation. Lactate is a potent inhibitor of food intake (39, 40), and our studies and others have shown that lactate production by skeletal muscle is greatly increased during lactation (12). To counteract the elevated levels of circulating lactate, a decrease in MCT1 expression in the ARH would prevent a lactate-induced hypophagia.
Somatotropic signaling is another critical component of the anorexia system in the hypothalamic regulation of energy homeostasis. Levels of IGFBP3 are important indicators of states of altered GH secretion and IGF action (21, 41). The high level of expression of IGFBP3 in ARH/VMH during lactation, by binding IGF, would result in a sustained low level of IGF. Given that intracerebroventricular injection of IGF results in a decrease in food intake and body weight (41), low levels of IGF would act to sustain the appetite signal. IGF and GH levels have been shown to be low during lactation in some species (42). Our data showed that GHRH gene expression in the ARH was significantly decreased during lactation (Table 3
). The reduced GHRH is most likely due to the increase in ARH NPY, because increases in NPY have been shown to inhibit GHRH (43).
Taken together, the sustained hyperphagia during lactation could be achieved through increases in the orexigenic peptides, NPY and AGRP, as well as glucocorticoid signals, and the decreases in the anorexigenic signals, CART, dynorphin, and POMC. At the systemic level, the suppression in serum leptin would enhance all of the alterations in the various appetite-regulating neuropeptide systems that are observed in the ARH (44). It is also likely that the changes in ARH/VMH gene expression involved in the sustained hyperphagia are linked to the suppression of GnRH activity and cyclic reproductive function that is characteristic of lactation.
New differentially regulated genes in ARH that may be involved in several adaptations occurring during lactation
In addition to increased appetite and cessation of reproductive cyclicity, there are other important adaptations during lactation that are under the influence of complex neuroendocrine control mechanisms (2, 18, 45, 46, 47). The decrease in the endothelin system (ECEL1 and EDNRB) during lactation is likely associated with a decrease in endothelin production, although the role of endothelin as a neurotransmitter is presently not clearly understood (20). Another interesting differentially regulated gene identified by these studies is GPR88, an orphan G protein-coupled receptor, although the ligand for this receptor is yet to be identified. Furthermore, removal of the pups for 48 h fully reversed the changes in GPR88 occurring during lactation, indicating that this change correlates with the hyperphagia and suppressed reproductive function that occurs during lactation. Interestingly, several structural proteins, GFAP, and CTGF, as well as CR16, were differentially expressed during lactation, indicating that extracellular matrix and synaptic remodeling may be occurring in ARH/VMH (47).
In summary, by using Affymetrix microarray analysis and a stringent method of data analysis, we acquired a reliable profile of the lactation-induced transcriptional changes in the ARH/VMH areas. These differentially regulated gene transcripts during lactation are involved in a variety of functions, such as metabolic homeostasis, milk production, suppression of cyclic reproductive function, maternal behavior, and structure remodeling. Some of the changes include previously unheralded transcription factors. These data identify additional candidate genes that may be involved in ARH-mediated changes in neuroendocrine function during lactation. Additional studies are required to establish the specific cell types within the ARH/VMH that express these genes and to determine their functional significance.
| Acknowledgments |
|---|
| Footnotes |
|---|
First Published Online July 7, 2005
1 X.Q.X. and K.L.G. contributed equally as senior authors. ![]()
Abbreviations: AGRP, Agouti-related protein; ARH, arcuate nucleus of hypothalamus; CART, cocaine- and amphetamine-regulated transcript; Cont, control group; CTGF, connective tissue growth factor; ECEL1, endothelin-converting enzyme-like 1; EDNRB, endothelin receptor type B; GFAP, glial fibrillary acidic protein; GPR88, G protein-coupled receptor 88; IGFBP3, IGF binding protein 3; Lac+8, suckled until d 11; Lac+0, removed at 1500 h on d 9; MCT, monocarboxylate transporter; NPY, neuropeptide Y; OVX, ovariectomy/ovariectomized; POMC, proopiomelanocortin; PVH, paraventricular nucleus; SGK, serum/glucocorticoid-induced kinase; SON, supraoptic nucleus; TH, tyrosine hydroxylase; VMH, ventromedial nucleus of the hypothalamus.
Received May 9, 2005.
Accepted for publication June 30, 2005.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. Arumugam, D. Fleenor, and M. Freemark Lactogenic and Somatogenic Hormones Regulate the Expression of Neuropeptide Y and Cocaine- and Amphetamine-Regulated Transcript in Rat Insulinoma (INS-1) Cells: Interactions with Glucose and Glucocorticoids Endocrinology, January 1, 2007; 148(1): 258 - 267. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Endocrinology | Endocrine Reviews | J. Clin. End. & Metab. |
| Molecular Endocrinology | Recent Prog. Horm. Res. | All Endocrine Journals |