|
|
||||||||
1 Endocrine Research Laboratory, Aurora St. Luke's Medical Center, Milwaukee, Wisconsin
2 Department of Medicine, Medical College of Wisconsin, Milwaukee, Wisconsin
3 Department of Biology, Boston University, Boston, Massachusetts
| ABSTRACT |
|---|
|
|
|---|
steroidogenesis; neonatal; mRNA; adrenal cortex; corticosterone; ACTH
| INTRODUCTION |
|---|
|
|
|---|
In the neonatal rat, there is a unique developmental aspect to the maturation of the adrenocortical response to stressors such as hypoxia (30, 34). We have previously shown that hypoxia from birth to postnatal day (PD) seven significantly increased plasma corticosterone and adrenal sensitivity to ACTH in suckling rat pups (25, 26). This increase was largely independent of increased plasma ACTH concentrations and was partially attenuated by the administration of guanethidine, a postganglionic sympathetic nerve terminal blocker that prevents the release of norepinephrine from nerve terminals and affects adrenocortical responsiveness (27, 33). Furthermore, the increase in corticosterone occurred during the so-called stress hyporesponsive period of the neonatal rat (30, 34, 36). This period of relative suppression of the HPA axis is thought to play a protective role during a critical period of brain development.
The molecular and cellular aspects of hypoxia-induced increases in plasma corticosterone have also been explored. We have shown that hypoxia elicited small but significant increases in steroidogenic acute regulatory (StAR) protein expression (25). Other studies from our laboratory measured unique changes in the lipid and fatty acid profile of the hypoxic adrenal (6). Interestingly, hypoxia significantly increased the total concentration of cholesterol ester (CE) in the adrenal, with specific increases in CE-associated 20:4n6 (arachidonic acid) concentrations. Previous studies using semiquantitative PCR revealed no changes in the expression of genes encoding proteins involved in the steroidogenic pathway (25, 26). It was therefore pertinent to address changes in gene expression using more advanced screening approaches.
The present study examined the effects of hypoxia from birth to PD7 on adrenal gene expression using microarray analysis. This approach was used to search for hypoxia-induced changes in the expression of novel genes not necessarily considered a component of the steroid synthetic pathway. Microarray results considered highly significant were verified with real-time PCR. In addition, expression of genes involved with the steroidogenic pathway was also analyzed by real-time PCR. Food intake of nursing dams was measured daily to examine the interaction between the anorectic effects of hypoxia and HPA axis responsiveness in pups. Nursing normoxic dams were pair-fed an amount of food that matched the intake of hypoxic dams. Plasma ACTH and corticosterone, as well as changes in adrenal gene expression, were measured in the three treatment groups (Normoxia, Hypoxia, and Normoxia Pair Fed) as functional endpoints of HPA axis activity. We hypothesized that hypoxia from birth would increase the expression of genes involved in steroidogenesis or the steps that lead up to steroidogenesis. We also hypothesized that maternal food restriction per se may contribute to changes in plasma ACTH or corticosterone.
| METHODS |
|---|
|
|
|---|
Total RNA for microarray analysis.
Total RNA from pooled adrenal glands (n = 4 pooled samples per Normoxia or Hypoxia) was isolated with TRIzol reagent (Invitrogen Life Technologies, Carlsbad, CA) and column purified (Qiagen, Valencia, CA) as described previously (16). RNA quality was assessed spectrophotometrically on the basis of the A260/A280 ratio. All RNA samples were checked for integrity of 18S and 28S RNA by gel electrophoresis.
Microarray analysis.
Changes in adrenal gene expression due to hypoxia were measured by microarray analysis, as described previously (16). In brief, double-stranded DNA was synthesized from 10 µg of total RNA using a Superscript cDNA Synthesis Kit (Invitrogen). Biotin-labeled cRNA was generated by transcription with T7 polymerase and purified on RNeasy affinity columns (Qiagen). Fragmented, biotinylated cRNA (10 µg), along with hybrid controls (Affymetrix, Santa Clara, CA), were hybridized to the Affymetrix Rat Genome U34A GeneChip array containing probes for 15,923 transcripts. Arrays were washed and stained and then scanned at 488 nm in a G2500A GeneArray Scanner (Agilent, Palo Alto, CA).
Microarray data quantification, normalization, and analysis.
Evaluation of data from microarray analysis was performed as described in detail previously (16). In brief, scanned images were quantified by GeneChip Operating Software 1.1 (Affymetrix). Signal intensities were used to determine overall expression level and a detection confidence score. Fold-change in expression was calculated from the average signal intensity of each group. A modified version of the false-discovery rate (FDR) method was used to assess significance of results (16). To identify genes significantly affected by hypoxia, genes with an FDR
0.25 and a minimum fold-change of ±1.3 were selected. Functional categories used to sort genes according to biological function were derived from the Expression Analysis Systematic Explorer software. More specific descriptions of gene function were obtained from the National Center for Biotechnology Information Entrez Gene Database (http://www.ncbi.nlm.nih.gov). The data comply with the MIAME standard, and the accession numbers are GSM143946 through GSM143953.
Real-time PCR analysis.
Total RNA for real-time PCR was isolated from a separate set of pooled adrenal tissue [one litter (24 adrenals) per pooled sample; n = 4 pooled samples per treatment group] using the RNeasy Mini Protocol (Qiagen). The concentration of RNA was quantified using the absorbance value at 260 nm and the quality of the sample preparation was assessed using the A260/A280 ratio. All RNA preparations were diluted to a final concentration of 10 ng/µl. Real-time PCR was performed using the Taqman One-Step RT-PCR Protocol [Applied Biosystems (ABI), Foster City, CA]. Premade primers and probes were purchased from ABI (see Table 1), and expression of 18S RNA was used as an endogenous control. The final reaction volume of 25 µl consisted of 1x AmpliTaq Gold DNA Polymerase mix, 1x RT enzyme mix containing MultiScribe Reverse Transcriptase and RNase Inhibitor, 1x primer/probe mix, and 50 ng of total RNA. Amplification and detection were performed with the ABI Prism 7300 Sequence Detection System with the following thermal cycler conditions: 48°C for 30 min (RT), 95°C for 10 min, and 40 cycles at 95°C for 0.25 min and 60°C for 1 min. Each measurement was carried out in triplicate or quadruplicate. Differences in gene expression, expressed as fold-change, were calculated using the 2
Ct method using 18S as the internal control.
|
Statistics.
Data for food intake measurements were analyzed by two-way repeated measures analysis of variance (ANOVA) with P < 0.05 considered significant. Plasma hormone data were analyzed with one-way ANOVA with P < 0.05 considered significant. The
Ct values obtained from real-time PCR, those used for calculation of fold-change differences (expressed here as percent of control), were also used for statistical analyses. PCR data were analyzed with one-way ANOVA with P < 0.05 considered significant. All post hoc analyses were performed by Student-Newman-Keuls method for multiple comparisons (SigmaStat 2.03). Explanation of statistical analyses used for data generated by microarray may be found above.
| RESULTS |
|---|
|
|
|---|
|
|
|
Table 3 lists real-time PCR results, expressed as a percent of control (Normoxia). Calculated
Ct values are also included in Table 3. Real-time PCR was initially used to verify selected microarray results. We selected 12 genes on the basis of their possible contribution to the functional endpoint of the study (i.e., increased plasma corticosterone). Three of the 12 PCR results did not agree with those indicated by microarray (Gmpr, Rgs2, and Cpt2). This discordance was confirmed by a second PCR assay for each of the three genes. The other nine results from real-time PCR showed trends that agreed with the microarray; however, only two of them reached statistical significance (Cyp1b1 and Nsf).
|
|
| DISCUSSION |
|---|
|
|
|---|
There was a small but significant increase in plasma ACTH concentrations in hypoxic pups in the current study. We have previously shown either no or only a small increase in ACTH using our experimental model (2527). We are confident that the large increase in plasma corticosterone is due to increased adrenal sensitivity to ACTH, at least partly mediated by the sympathetic nervous system (27). In addition, changes in adrenal gene expression described below were probably not mediated by small increases in plasma ACTH, since similar changes were not observed with pharmacological doses of ACTH in neonatal rat pups (16). Therefore, the effects described next appear to be primarily mediated by hypoxia.
Microarray results.
In an effort to further the understanding of hypoxia-induced increases in plasma corticosterone, a broad approach (i.e., microarray) was used to analyze changes in adrenal gene expression. A similar approach has been used to study the effects of chronic ACTH treatment in rat pups (16). None of the genes that encode proteins of the corticosterone synthetic pathway were significantly affected by hypoxia according to this method. Chronic ACTH treatment in rat pups was also without effect on these genes (16). Two genes whose products are involved in further metabolism of corticosterone, Srd5a1 (steroid 5
-reductase) and LOC191574 (3
HSD), were significantly increased and decreased by hypoxia, respectively. An increase in plasma concentrations of 5
-dihydrocorticosterone (5
THB, the product of steroid 5
-reductase) could influence physiological responses to hypoxia, since it has been shown that 5
THB binds to and activates the glucocorticoid receptor (20). A decrease in the expression of Hsd11b2 (11ßHSD2) was also detected and may indicate decreased conversion of corticosterone to its inactive metabolite in the adrenal cortex.
Hypoxia affected the expression of a number of genes encoding proteins that take part in intermediary metabolism and oxidative phosphorylation (microarray data). Expression of Aldoc (aldolase C) and Eno1 (enolase 1
), enzymes involved in the glycolytic pathway, was significantly increased. In addition, the expression levels of Got1 (glutamate oxaloacetate transaminase 1), Cpt2 (carnitine palmitoytransferase 2), Fmo1 (flavin-containing monooxygenase 1), and Cox4b (cytochrome oxidase c, subunit 4b), all involved in mitochondrial respiration, were increased by hypoxia. These findings suggest that hypoxia may increase substrate available for mitochondrial respiration, thus fostering the concomitant increase in steroid production. Recent studies provided evidence that Leydig cell steroidogenesis requires actively respiring mitochondria that maintain the mitochondrial membrane potential (
m) (1).
Interestingly, the expression of Ucp1 (uncoupling protein 1) was significantly decreased in response to hypoxia. Ucp1 dissipates the proton gradient required for mitochondrial ATP synthesis and is involved in thermogenesis but is normally associated with mitochondria of adipose, heart, and muscle tissue (11, 21). Our detection of Ucp1 mRNA in the neonatal rat adrenal may be a novel finding. A decrease in Ucp1 expression, along with the changes noted above, suggests that neonatal adrenal mitochondria may operate more efficiently when exposed to a moderately low oxygen environment. This phenomenon could augment steroid production as proposed above. An argument against this mitochondrial theory would be the nearly eightfold decrease in Cox8h (cytochrome c oxidase, subunit VIII-H); however, these findings must be interpreted with caution. It is difficult to assign cause and effect when dealing with a multitude of changes in gene expression as in the present study.
Since whole adrenal glands were pooled for microarray and real-time PCR analyses, a brief discussion of genes linked to adrenal neural and adrenomedullary function is pertinent. It is also pertinent since we have previously shown that hypoxia-induced increases in plasma corticosterone were attenuated with guanethidine administration, which blocks postganglionic sympathetic innervation of the adrenal cortex (27, 33). Subsequent studies failed to measure any hypoxia-induced changes in plasma catecholamine concentrations, although there was a tendency toward increased plasma epinephrine (5). Hypoxia increased the expression of Dbh (dopamine ß-hydroxylase), Nsf (N-ethylmaleimide sensitive fusion protein), and Slc18a2 (solute carrier family 18, member 2) in the present study. Dbh encodes the enzyme responsible for the conversion of dopamine to norepinephrine (NE), and an increase in its expression may indicate hypoxia-induced adrenal NE synthesis. The proteins encoded by Slc18a2 and Nsf function in vesicular monoamine transport, a prerequisite to the exocytic pathway involved in monoamine release (2, 7, 12, 29). Taken together with our previous studies with guanethidine (27), these findings may indicate increased adrenocortical sympathetic input and/or increased chromaffin cell catecholamine production during hypoxia. These changes, in concert with increases in steroidogenic gene expression, may contribute to the stimulation of corticosterone production. An interesting possibility is that increased corticosterone during hypoxia may alter adrenomedullary gene expression and function via paracrine effects (37). Additional studies will be required to evaluate this possible mechanism.
Real-time PCR vs. microarray: concordant results.
A dozen genes that were significantly affected by hypoxia, according to microarray results, were chosen for confirmation with real-time PCR. Genes were selected based on involvement in intra- and intercellular pathways that may have an influence on adrenal steroid production (Table 3). We confirmed that hypoxia significantly increased the expression of Cyp1b1 and Nsf (N-ethylmaleimide sensitive factor). Cyp1b1 is probably not involved directly in steroidogenesis (3); however, it may be involved in O2 metabolism as it does alter the metabolism of aromatic hydrocarbons. The stimulation of Cyp1b1 expression during hypoxia overrides the normal developmental suppression of the expression of this gene, as shown previously (3). There was a tendency for increased expression of Xdh, Srd5a1, and Eno1 according to real-time PCR, but the differences were not considered significant. There was also a tendency for decreased expression of Cdh1, Pla2g4a, Aox1, and Cyp11b3. While agreeing with the trends measured by microarray, these PCR results did not reach statistical significance.
Real-time PCR vs. microarray: discordant results.
Comparison of microarray and real-time PCR results revealed discordance for three of the 12 genes selected (Gmpr, Rgs2, and Cpt2). In the case of Rgs2, microarray results indicated a decrease in expression (1.5-fold), while real-time PCR revealed no change. The expression level of Cpt2 was increased by hypoxia according to the microarray, but it was decreased when assayed with real-time PCR. In the third case, microarray detected a significant increase in Gmpr expression, while real-time PCR results indicated a decrease, albeit not significant. We verified real-time PCR by running a second assay for each gene. The discordance of these results may be explained by the decreased sensitivity of the microarray per se or by the relatively low expression levels of these genes.
Expression of steroidogenic genes and functional correlates.
Using real-time PCR, a more sensitive method of detection, we measured hypoxia-induced changes in the expression of genes involved in adrenocortical function. Hypoxia from birth to PD7, compared with Normoxia (food ad libitum), tended to increase the expression of Star, Cyp11a, and Ldlr. The only significant increases in expression were with Cyp21a1 and Cyp1b1. It is possible that the measured increase in Cyp21a1 mRNA expression was due to the detection of genomic DNA (see Table 1). These results corroborate the finding of increased plasma corticosterone concentrations and increased sensitivity to ACTH (25, 27) and indicate that increased steroid production due to hypoxia is at least partly mediated by changes in adrenal gene expression. We have previously shown that zona fasciculata cells from 7-day-old pups exposed to hypoxia from birth exhibit increased corticosterone production (compared with cells from normoxic pups) when studied in vitro (24). These increases were probably driven by increased early pathway activity (i.e., pregnenolone production) (26). This fits with the current findings of increased expression of Ldlr, Star, and Cyp11a mRNA, as these genes encode proteins involved in cholesterol uptake, transport into mitochondria, and metabolism to pregnenolone, respectively.
Expression of hypoxia-responsive genes.
Hypoxia had no effect on the expression of Hif1a (hypoxia-inducible factor 1
) or Vegf (vascular endothelial growth factor), although the expression of Vegf tended to be increased by hypoxia. The fact that Hif1a mRNA expression was not altered by hypoxia does not indicate that HIF-1
protein expression also remained unaffected. It is the protein itself that is stabilized in response to acute decreases in oxygen concentration (31), and the transcriptional effects of this stabilization may normalize at some point during the hypoxic treatment used in the present study. The Vegf gene is a target of the HIF-1
protein (31), and expression of Vegf mRNA would be expected to increase in hypoxic tissues to maximize perfusion. Hypoxia has been shown to increase adrenal blood flow in dogs (4), and the present results with Vegf expression, though not statistically significant, may support this phenomenon.
Hypoxia and maternal food restriction.
Our laboratory has consistently shown that hypoxia from birth to PD7 attenuates weight gain in rat pups (5, 23, 27). The nature of our experimental model necessitates placing the dam in the hypoxic chamber with her pups. Therefore, examination of the effects of hypoxia on the dam is necessary. It was not surprising to find that hypoxia-induced anorexia was most profound at the onset of the hypoxic exposure (PD1PD2). Hypoxic dam food intake remained significantly lower than normoxic controls throughout the study, even though hypoxic dams increased their food intake daily from PD36. Pups of normoxic dams pair-fed to match hypoxic dam intake did not exhibit any changes in plasma ACTH or corticosterone. Body weights of these pups were not different from those reared by dams fed ad libitum. This confirmed that the effects of hypoxia on the neonatal HPA axis were not due to reduced maternal food intake.
When the anorexic effect of hypoxia on the dam was accounted for (Hypoxia vs. Normoxia Pair Fed) with the real-time PCR results, the magnitude of the increases in gene expression became greater and more genes (Star, Cyp21a1, Ldlr, and Cyp1b1) were significantly different. This phenomenon suggests that hypoxia-induced decreases in maternal food intake did not play a major role in the modulation of adrenal function in pups. Pups of pair fed dams did not exhibit elevated plasma corticosterone concentrations, the functional endpoint of adrenocortical function. Maternal food restriction in lactating dams has been shown to decrease the plasma concentration of corticosterone-binding globulin (CBG) in pups and was thought to play a reinforcing role in the suppression of HPA axis activity in the pups (18). We have not found any changes in CBG concentrations in pups in our experimental model of neonatal hypoxia (25).
In summary, the present study used microarray and real-time PCR analyses to examine changes in adrenal gene expression in the hypoxic neonatal rat. We also measured the effect of hypoxia on maternal food intake, accounting for the anorectic effects of hypoxia. Hypoxia from birth to 7 days of age significantly increased plasma corticosterone and elicited a small but significant increase in plasma ACTH. Matching normoxic dam food intake to that of hypoxic dams had no effect on plasma ACTH or corticosterone concentrations in 7-day-old pups. Real-time PCR analysis measured a significant increase in the expression of Cyp21a1 in the hypoxic adrenal. The effects of hypoxia became more pronounced when accounting for hypoxia-induced anorexia (e.g., significant increases in Star, Cyp21a1, and Ldlr expression). Microarray results revealed significant changes in the expression of genes involved in both adrenocortical and adrenomedullary function. We postulate that hypoxia induces shifts in adrenal cellular function such as mitochondrial respiration and vesicular trafficking that, together, act in concert with augmented cholesterol uptake to increase sensitivity to ACTH and corticosterone output.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Article published online before print. See web site for date of publication (http://physiolgenomics.physiology.org).
| REFERENCES |
|---|
|
|
|---|
-reduced glucocorticoids, novel endogenous activators of the glucocorticoid receptor. J Biol Chem 279: 2290822912, 2004.This article has been cited by other articles:
![]() |
H. Raff, R. Molthen, W. Pan, and A. J. Kastin Commentaries on viewpoint: effect of altitude on leptin, does it go up or down? J Appl Physiol, November 1, 2008; 105(5): 1686 - 1687. [Full Text] [PDF] |
||||
![]() |
E. D. Bruder, J. K. Taylor, K. J. Kamer, and H. Raff Development of the ACTH and corticosterone response to acute hypoxia in the neonatal rat Am J Physiol Regulatory Integrative Comp Physiol, October 1, 2008; 295(4): R1195 - R1203. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |