| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Department of Medicine and Committee on Cancer Biology, University of Chicago, Chicago, Illinois 60637
Address all correspondence and requests for reprints to: Suzanne D. Conzen, M.D., Department of Medicine and Ben May Department for Cancer Research, University of Chicago, 5841 South Maryland Avenue, MC 2115, Chicago, Illinois 60637. E-mail: sdconzen{at}uchicago.edu.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
As an approach to dissecting GCs influence on gene transcription networks in different cell types, our laboratory is examining gene expression data over time after GC treatment. In this way we hope to untangle direct transcriptional events after glucocorticoid receptor (GR) activation, as well as to uncover novel mechanisms of indirect GR-mediated regulation of gene expression. Studying gene expression at several time points has made it clear that there are distinctive groups of genes with varying temporal patterns of gene expression after GR activation. By clustering genes of similar patterns (8), we can predict possible regulatory pathways downstream of GR activation. For example, many genes whose expression levels change at a similar time point may be direct targets of a common GR-regulated transcription factor. By identifying those transcription factors and looking for common target genes, we are beginning to construct regulatory networks downstream of GR activation. These networks may help to uncover GRs indirect mechanisms for regulating gene expression, including the mechanisms underlying GC-dependent down-regulation of proinflammatory genes.
In a previous paper (9), we investigated the role of the transcription factor forkhead box O3a (FOXO3a) in GR-mediated cell survival signaling. GCs induce serum and GC regulated kinase (SGK) (10), a serine threonine kinase that is subsequently phosphorylated and activated in a phosphatidylinositol 3-kinase dependent fashion. After transcriptional induction by GR, activated SGK-1 phosphorylates FOXO3a, resulting in FOXO3a nuclear exclusion and inactivation (9, 11). Overexpression of FOXO3a also limits the ability of dexamethasone (Dex) to protect mammary epithelial cells from apoptosis. We have previously identified genes down-regulated at least 30% by GCs that contain putative FOXO3a binding sites in the 1000-bp upstream of the transcriptional start site (TSS) to begin to identify novel indirect targets of GR repression. More recently, the chromatin immunoprecipitation microarray (ChIP-chip) assay has emerged as a genome-wide tool for studying transcription factor-chromatin interaction. In the current study, we additionally use FOXO3a ChIP-chip results to refine our identification of FOXO3a target genes. Interestingly, one of the most significantly down-regulated genes also found using the FOXO3a ChIP-chip assay was urokinase plasminogen activator (uPA or PLAU), a very well-established down-regulated target of GC action that is not a known FOXO3a target. These results suggest a novel mechanism whereby GCs indirectly down-regulate uPA.
uPA plays an important role in several diseases, most notably cancer. uPA is a serine protease that generates plasmin, causing matrix metalloproteinase activation and degradation of the extracellular matrix. In breast cancer, uPA and the plasminogen system are thought to degrade the extracellular matrix, mediating tumor cell invasion (12). Furthermore, elevated uPA activity correlates with poor breast cancer prognosis (13). GCs suppress uPA mRNA expression and activity in breast cancer cells (14, 15). Down-regulation of uPA via inhibition of activator protein-1 (AP-1) and nuclear factor
-B (NF-
B) function (16, 17) has been long postulated as the major mechanism. In this study we identify FOXO3a inactivation as an additional important mechanism for uPA repression in breast epithelial cells.
| Materials and Methods |
|---|
|
|
|---|
Both MCF10A-Myc mammary epithelial and MDA-MB-231 breast cancer cell lines were chosen because they express abundant endogenous GR despite the absence of estrogen, progesterone, and Her-2/neu receptors.
Retroviruses were made by transient transfection of retroviral vectors into amphotropic Phoenix cells (a gift of Dr. Gary Nolan, Stanford University, Palo Alto, CA) using Polyfect Transfection Reagent (QIAGEN, Inc., Valencia, CA) per manufacturers instructions. Retroviruses were used to create stable MCF10A-Myc or MDA-MB-231 cells expressing pLPCX alone, pLPCX- wild-type (WT)-FOXO3a, pLPCX-triple mutant (TM) FOXO3a, scrambled sequence (SS) short hairpin RNA (shRNA), or SGK shRNA; colonies were selected with puromycin as described previously (18). Fugene HD (Roche Diagnostics Corp., Roche Applied Science, Indianapolis, IN) was used following the manufacturers protocol with a 5:2 Fugene to DNA ratio to create C-28/I2 cells transiently expressing pLPCX alone, or pLPCX containing WT-FOXO3a, or TM-FOXO3a.
Plasmids
WT-FOXO3a-hemagglutinin (HA), TM-FOXO3a-HA were obtained from Dr. Michael E. Greenberg (Harvard Medical School, Boston, MA) and inserted into pLPCX retroviral vectors as previously described (9). SGK shRNA and SS shRNA constructs were created in the retroviral expression vector RNAi-Ready pSIREN-RetroQ (BD Biosciences, Palo Alto, CA) as described previously (9).
Drug treatment
Dex (Sigma-Aldrich) was dissolved in 100% ethanol to make stock solutions. Before treatment, cell culture dishes were washed twice with 1x PBS, placed in media without serum, and starved 24 h for MDA-MB-231 cells or C-28/I2 cells or 72 h for MCF10A-Myc and SK-BR-3 cells before drug treatment. Cells were treated for the time specified in the experiment with Dex (10–6 M) or an equal volume of ethanol.
Time course microarray data
The procedure for cRNA preparation and hybridization to high-density oligonucleotide arrays was described previously (19). The experiment was previously performed (9) on MCF10A-Myc cells starved for 72 h, followed by treatment with either vehicle (ethanol), Dex (10–6 M), or concomitant Dex/RU486 (10–7 M) for 0.5, 2, 4, or 24 h. Preparation of the microarrays (Affymetrix human genechip HG-U133A; Affymetrix, Inc., Santa Clara, CA) was performed at the University of Chicago Microarray Core Facility. The genechip used contains 22,215 probe sets representing approximately 16,000 human genes. The experiment was performed three independent times. Fold-change values (Dex vs. vehicle) were determined after analysis using Robust Multiarray Average algorithm (40).
ChIP-chip analysis
Ethanol and Dex treated immunoprecipitated chromatin and cellular input DNA were produced as in the conventional chromatin immunoprecipitation (ChIP) assay. These were then amplified, fragmented, and labeled according to the Affymetrix ChIP Protocol. These samples were then hybridized to the Affymetrix GeneChip Human Promoter 1.0R Array (P/N 900775) and scanned at the University of Chicago Functional Genomics Core Facility.
Probe signals were analyzed using the software Model-based Analysis of Tiling-array (MAT) (20) to detect enriched regions of transcription factor binding based on National Center for Biotechnology Informations build 36 of the human genome. Data were normalized using the corresponding input sample. A window size of 400 bp was set based on observed DNA shearing. A threshold of P < 10–3 was used to determine which regions were occupied by FOXO3a, and the resulting browser extensible display (".bed") file was analyzed by the cis-regulatory element annotation system (21) to generate a list of genes with putative FOXO3a chromatin interaction in their core promoter (–5000/+500). Only regions that appeared to be occupied by FOXO3a in the ethanol vs. the Dex samples were identified in the final list.
Quantitative real-time RT-PCR
Total RNA was extracted, using the RNeasy Mini kit (QIAGEN) from cells after treatment. cDNA was synthesized from RNA with TaqMan RT reagents (Applied Biosystems, Foster City, CA) per manufacturers instructions in a 7300 Real-Time PCR System (Applied Biosystems). Quantitative real-time RT-PCRs were performed as previously described (9). The following primers were used: uPA (PLAU, NM_002658), 5'-ACGCAAGGGGAGATGAAG-3' (forward) and 5'-TCAGCAAGGCAATGTCGT-3' (reverse); and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (NM_002046), 5'-GAGTCAACGGATTTGGTCGT-3' (forward) and 5'-TTGATTTTGGAGGGATCTCG-3' (reverse). GAPDH was amplified as an internal control. The samples were loaded in triplicate for each primer pair, and the results of each sample were normalized to GAPDH. Fold change was calculated as a ratio of treatment over control (vehicle). Agarose DNA gels were run, or Dissociated Curve analysis was performed to verify the absence of nonspecific products.
ChIP assay
MCF10A-Myc cells were starved for 72 h and then treated with either vehicle (ethanol) or Dex (10–6 M). With slight modifications, manufacturers directions in the Chromatin Immunoprecipitation Assay Kit (Upstate Biotechnology, Waltham, MA) were followed. Cells were fixed with formaldehyde for 20 min, and lysates were sonicated on ice four times for 30 sec each. Lysates were precleared for 1 h, after which 1% of the lysates were taken for input. Lysates were incubated with anti-FOXO3a antibody (Santa Cruz Biotechnology, Inc., Santa Cruz, CA) or antirabbit IgG (Santa Cruz Biotechnology). DNA was purified with a PCR Purification Kit (QIAGEN). PCR was performed on the purified DNA and analyzed on agarose gels. Primers testing a known region of TNF-related apoptosis inducing ligand (TRAIL) were used as a positive control: forward, 5'-CCTGGGCGATAAAGTGAGAT-3'; reverse, 5'-GGCCCAGCTGTATGTTGTCT-3'.
To test the uPA binding sites, three primer pairs were designed to test the different regions of the uPA promoter. Because the resolution of a ChIP assay is limited by the size of DNA shearing, it is unlikely that the ChIP assay could distinguish further among these four binding sites. The uPA primers used were: uPA primer pair 1, forward 5'-TCTGTTTATGTTCTGGTCACGTTT-3' and reverse 5'-CACTGTCTTCTGGAGAGACTTCTG-3'; uPA primer pair 2, forward 5'-TGGAATGGGCTATACAGAGATG-3' and reverse 5'-TGGAGGACAAAATAAGAGGAAGA-3'; and uPA primer pair 3, forward 5'-ATCAAGCCCTGTCAAAAACG-3' and reverse 5'-ACCAGTGGATGGAACAAAGC-3'.
Western blot analysis
MCF10A-Myc cells, MDA-MB-231 cells, or C-28/I2 cells were starved of growth factors for times previously described, and then treated with Dex (10–6 M) or vehicle (ethanol) for 24 h. Cell lysates were harvested with 2x Laemmli lysis buffer or as directed for the ChIP assay. Equal amounts of protein were resolved by SDS-PAGE, transferred to a nitrocellulose membrane, blocked by incubation in 5% skim milk in 0.1% Tween 20 in Tris-buffered saline, and probed with a 1:1000 dilution of anti-SGK-1 (22), anti-FOXO3a (Santa Cruz Biotechnology), anti-P-FOXO3a (Santa Cruz Biotechnology), anti-uPA (American Diagnostica Inc., Stamford, CT), anti-Na+/K+ ATPase (Affinity BioReagents, Inc., Golden, CO), anti-
-tubulin (Calbiochem, San Diego, CA), or anti-β-Actin (Sigma-Aldrich). After washing with 0.1% Tween 20 in Tris-buffered saline and incubation with peroxidase-conjugated goat antirabbit or antimouse secondary antibody (1:5000 dilution), the membranes were treated with enhanced chemiluminescent staining (Amersham Biosciences Inc., Piscataway, NJ) per manufacturers instructions before film development.
Luciferase constructs
A 1270-bp fragment of the uPA promoter (–1252/+18) was generated with Kpn-1 and Bgl-II restriction sites using standard nested PCR and inserted into the luciferase vector pGL3-basic (Promega Corp., Madison, WI). The primers used were: left outer primer (5'-CTGCAAAGGAAGGAGAAGTCAGGGCAAG-3'); right outer primer (5'-CACAGAAGTCTCTCCAGAAGACAGTGGGTCA-3'); left inner primer (5'-GGGGTACCCAGAGATGCGTGTGTGTGTG-3'); and right inner primer (5'-GAAGATCTCTGCGGTCTCCGCACTGT-3').
The fragment of the human uPA gene promoter was then mutated to introduce an A to G mutation (mut) in the core (AAACA) FOXO3a consensus binding site (AGACA) by site-directed mutagenesis (QuikChange site-directed mutagenesis kit; Stratagene, La Jolla, CA) per manufacturers instructions using the primers 5'-GTTAAGAGTGTAAAAGGACAGACAGTGAAAAATAAATCTTCC-3' (forward) and 5'-GGAAGATTTATTTTTCACTGTCTGTCCTTTTACACTCTTAAC-3' (reverse). Underlining indicates which nucleotide was mutated. The resulting vectors (WT- and mut-uPA-pGL3) were sequenced for validation.
Luciferase assay
Thirty-five millimeter plates of SK-BR-3 cells were transfected with 0.9 µg WT-uPA-pGL3, mut-uPA-pGL3, or pGL3 vector alone together with 0.5 µg of either pLPCX, WT-FOXO3a, or TM-FOXO3a, and 0.1 µg of a pCMV-β-galactosidase-encoding vector (gift of Dr. Geoffrey Greene, University of Chicago, Chicago, IL). Forty-eight hours after transfection using PolyFect (QIAGEN) per manufacturers instructions, each transfection was split into multiple wells to provide triplicates for luciferase measurement. The protocol for the Luciferase Assay Kit (Promega) was followed to collect and prepare the luciferase lysates. Luminescence was measured using a Wallace luminescence counter (Beckman Coulter, Inc., Fullerton, CA), and β-galactosidase activity was measured separately at 420 nm. The relative luciferase activity in each condition was normalized to β-galactosidase activity.
Statistical analysis
For all gene expression and luciferase assays, statistical significance between two conditions was determined by an unpaired two-tailed t test. Results were considered statistically significant if the P value was less than 0.05.
| Results |
|---|
|
|
|---|
|
We then compared (Fig. 1B
) our list of significantly down-regulated genes after GC treatment in the same MCF10A-Myc cell line (9) with the list of genes identified by ChIP-chip. Although the ChIP-chip technique provides evidence for complex formation between the immunoprecipitated transcription factor and DNA, the functional significance of this interaction is not always clear (23). By examining only those genes down-regulated after GR activation (potentially with reduced transcription secondary to SGK-1-mediated inactivation of FOXO3a), we attempted to reduce false positives and to increase physiologically relevant FOXO3a target genes. These data revealed 11 genes (Fig. 1C
) as likely regulated by the GR through an indirect mechanism involving FOXO3a inactivation. These genes met several stringent criteria to be considered as putative FOXO3a targets, suggesting that they deserved further study. The known FOXO3a target TRAIL appeared on the list, whereas IGF binding protein-3, previously identified as a FOXO3a target (9), was absent. Therefore, although this list may be useful as a starting point for future investigation, it is unlikely to encompass all functional FOXO3a targets due to the stringency of the ChIP-chip assay.
Of the genes identified, uPA was one of the most highly down-regulated genes, and also contains one of the most significant P values using MAT analysis of the ChIp-chip data. Because uPA plays a well-established role in both cancer and arthritis, we chose to study it as a potential FOXO3a target.
Dex down-regulates uPA mRNA expression at 24 h
To confirm the expression microarray data previously generated (9) for uPA (Fig. 2A
), semiquantitative real-time RT-PCR was used to examine uPA mRNA after treatment with Dex (10–6 M) or vehicle (ethanol) for 24 h in MCF10A-Myc cells. We also examined the MDA-MB-231 breast cancer cell line as, well as the C-28/I2 chondrocyte cell line because of the role uPA plays in both cancer and arthritis.
|
Conventional ChIP confirms uPA promoter occupancy by FOXO3a in breast cells
Figure 3A
presents the 2000-bp promoter region upstream of the uPA TSS, with the corresponding ChIP-chip data shown above the uPA promoter. In the ChIP-chip data, the vertical lines represent a score calculated for a window around each microarray probe by the MAT software, and the solid horizontal bar with an asterisk is the region determined by MAT to be a significant FOXO3a binding site (P = 2.74 x 10–4) in the vehicle but not the Dex-treated immunoprecipitates. Three primer pairs were then chosen to examine by conventional ChIP assay the four predicted FOXO3a binding sites (A–D) found in this region (Fig. 3A
).
|
As a positive control, we used PCR primers for a region of the TRAIL promoter previously shown to be bound by FOXO3a in the absence of Dex (26). For this primer pair, we confirmed FOXO3a promoter occupancy in vehicle but not Dex-treated immunoprecipitates (Fig. 3C
). In addition, uPA primer 1 did not amplify the far upstream uPA region, suggesting that this region is not a FOXO3a target. uPA primers 2 and 3 did exhibit vehicle-specific FOXO3a promoter occupancy. Together, the ChIP and ChIP-chip experiments suggest a region (–1262/–770 nt) of the uPA promoter that is occupied by FOXO3a in MCF10A-Myc cells.
Constitutively active FOXO3a increases uPA mRNA and protein steady-state levels and inhibits Dex-mediated uPA repression
Having identified putative FOXO3a binding sites in the uPA promoter, we then asked whether or not FOXO3a plays a functional role in uPA expression. MCF10A-Myc and MDA-MB-231 cells were stably transfected with either vector alone (pLPCX), pLPCX-WT-FOXO3a, or a pLPCX-constitutively active FOXO3a (TM-FOXO3a), in which three FOXO3a phosphorylation sites (Thr-32, Ser-253, and Ser-315) were mutated to alanines, thereby preventing phosphorylation and nuclear exclusion (27). C-28/I2 cells were transiently transfected with the same vectors.
In all three transfected cell lines, significant levels of ectopic FOXO3a proteins were observed (Fig. 4A
). As expected, there was significant FOXO3a phosphorylation in the cells expressing WT-FOXO3a, but not in the cells expressing TM-FOXO3a. We found that ectopic expression of either WT- or TM-FOXO3a significantly increased the steady-state levels of uPA mRNA in MCF10A-Myc and MDA-MB-231, but somewhat surprisingly, not in C-28/I2 cells (Fig. 4A
). This suggests that FOXO3a increases uPA expression in a cell type-dependent fashion.
|
We next examined whether or not WT-FOXO3a or TM-FOXO3a expression could reverse Dex-mediated down-regulation of uPA mRNA in these cell lines. We hypothesized that unlike the breast cell lines, constitutively active TM-FOXO3a would not affect GC-mediated uPA down-regulation in the C-28/I2 cells because TM-FOXO3a did not induce uPA expression in this cell type. Indeed, we found that the constitutively active TM-FOXO3a reversed the repression by GCs only in MCF10A-Myc and MDA-MB-231 cells (Fig. 4C
). Furthermore, GCs likely inactivated ectopic WT-FOXO3a activity through induction of SGK-1, whereas TM-FOXO3a remained active. In summary, these results suggest that FOXO3a induces uPA expression in breast epithelial and cancer cells, and that FOXO3a inactivation contributes to GC-mediated down-regulation of uPA mRNA and protein steady-state levels.
After GR activation, SGK-1 expression contributes to uPA down-regulation
SGK-1 and Akt-1 can coordinately mediate FOXO3a phosphorylation and inhibit FOXO3a activity by causing its nuclear exclusion (11). After Dex treatment we have previously observed phosphorylation and nuclear exclusion of FOXO3a that appears to be dependent on SGK-1 induction alone (9). To test the hypothesis that GC-mediated down-regulation of uPA mRNA requires SGK-1, we examined whether or not inhibiting SGK-1 expression affects uPA down-regulation. MCF10A-Myc and MDA-MB-231 cells were stably transfected with either SS shRNA or SGK shRNA. Western analysis demonstrated that SGK-1 was efficiently decreased in both cell lines (Fig. 5A
).
|
2- to 3-fold) uPA mRNA (Fig. 5A
We also examined uPA protein expression in the MDA-MB-231 cells by Western blot. In agreement with our mRNA results, intracellular pro-uPA protein was significantly decreased by Dex in control shRNA cells, but not in SGK shRNA cells (Fig. 5B
). The near full reversal of the Dex-mediated down-regulation of uPA protein, in comparison to the partial reversal of uPA mRNA levels (Fig. 5A
), suggests that SGK-1 may also have posttranscriptional effects on uPA expression.
FOXO3a can activate the uPA promoter-luciferase construct through a conserved FOXO3a binding sequence
All of the four identified putative FOXO3a sites in the 2000-bp region upstream of the uPA TSS were in the region identified by ChIP-chip (Fig. 3A
) and verified by conventional ChIP (Fig. 3C
). We further studied the central site (Fig. 3A
, site "C") to determine whether or not FOXO3a binds to and activates the uPA promoter. Luciferase constructs containing a 1270-bp (–1252/+18, Fig. 3A
) region of the uPA promoter were generated either containing the WT or a single nucleotide mutation in the core sequence (Fig. 6A
). The WT or mutant construct was then cotransfected with a pCMV-β-galactosidase reporter into SK-BR-3 breast cancer cells along with either a pLPCX control, WT-FOXO3a, or TM-FOXO3a vector.
|
Overall, these data support the hypothesis that SGK-1-mediated phosphorylation and inactivation of FOXO3a are important mechanisms contributing to decreased uPA expression after GR activation in breast epithelial cells. Interestingly, GCs appear to use a common mechanism, i.e. transcription factor inhibition (AP-1, NF-
B, and FOXO3a), to down-regulate genes in a cell type-specific manner (Fig. 7
).
|
| Discussion |
|---|
|
|
|---|
uPA expression also appears to correlate with the clinical severity of arthritis (29). GCs are used to alleviate arthritic symptoms and also down-regulate uPA in synovial cells (30), chondrocytes (25), and mononuclear cells (31), potentially accounting for their effectiveness in treatment. Because of this we examined the mechanism of uPA regulation by GCs in a human chondrocyte cell line. Interestingly, our data suggest that unlike breast cancer cells, uPA does not appear to be a FOXO3a target in chondrocytes. This is consistent with evidence suggesting that FOXO activity at target promoters is cell type dependent (32).
In cancer therapy, SGK-1 (33) and uPA (34) are considered to be potential therapeutic targets. Our data suggest that inhibiting SGK-1 might actually result in higher uPA levels, which could promote cancer cell invasion under certain conditions. Increased uPA production and concomitant cellular invasion might also be predicted upon inhibiting Akt-1, another phosphatidylinositol 3-kinase pathway member that phosphorylates and inactivates FOXO3a (35). In fact, consistent with the hypothesis that Akt-1 activation could result in down-regulated uPA expression via FOXO3a inactivation, Akt-1 activity inhibits tumor cell invasion in an in vitro assay (36). Akt-1s inhibitory effects on tumor cell invasion have been suggested to result from its ability to cause the ubiquitin modification and degradation of nuclear factor of activated T cells (37). In turn, nuclear factor of activated T-cell inactivation results in the down-regulation of the pro-invasion gene cyclooxygenase-2 (38). Our data suggest that the observed inhibition of invasion after Akt-1 expression could also be a result of decreased uPA expression. Furthermore, in an in vivo mouse model, constitutively active Akt-1 expression appears to have opposing effects on mammary gland tumor growth and metastatic progression: Akt-1 accelerates tumorigenesis, while suppressing tumor invasiveness (39).
The difference in mechanism of uPA regulation by GCs in breast epithelial cells vs. chondrocytes suggests that confirming bona fide FOXO3a target genes will require cell type-specific verification. Furthermore, the discovery that FOXO3a contributes to uPA regulation implies that in addition to its well-described role in promoting apoptosis in cancer cells, FOXO3a may also play a role in tumor invasion.
| Acknowledgments |
|---|
| Footnotes |
|---|
Disclosure Statement: The authors have nothing to disclose.
First Published Online January 31, 2008
Abbreviations: AP-1, Activator protein-1; ChIP, chromatin immunoprecipitation; ChIP-chip, chromatin immunoprecipitation microarray; Dex, dexamethasone; FOXO3a, forkhead box O3a; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; GC, glucocorticoid; GR, glucocorticoid receptor; MAT, Model-based Analysis of Tiling-array; mut, mutation; NF-
B, nuclear factor
-B; SGK, serum and glucocorticoid regulated kinase; shRNA, short hairpin RNA; SS, scrambled sequence; TM, triple mutant; TRAIL, TNF-related apoptosis inducing ligand; TSS, transcriptional start site; uPA, urokinase plasminogen activator; WT, wild type.
Received August 8, 2007.
Accepted for publication January 22, 2008.
| References |
|---|
|
|
|---|
B-mediated protease production by oral cancer. Biochim Biophys Acta 1586:11–22[Medline]
. Inflamm Res 46:60–64[CrossRef][Medline]
EF1 to activate growth suppressive genes in B lymphocytes. J Immunol 176:2711–2721
) can accelerate ErbB-2-mediated mammary tumorigenesis but suppresses tumor invasion. Cancer Res 64:3171–3178
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 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 |