| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Departments of Pathology (K.H.B., S.O., M.M.M.), Molecular and Human Genetics (K.H.B., M.M.M.), and Molecular and Cellular Biology (M.M.M.), Baylor College of Medicine, Houston, Texas 77030; and Department of Pharmacology, Case Western Reserve University (G.E.O., J.H.N.), Cleveland, Ohio 44106
Address all correspondence and requests for reprints to: Dr. Martin M. Matzuk, Department of Pathology, One Baylor Plaza, Baylor College of Medicine, Houston, Texas 77030. E-mail: mmatzuk{at}bcm.tmc.edu.
| Abstract |
|---|
|
|
|---|
knockout mouse models. Oligonucleotide array hybridization experiments using Affymetrix GeneChip technology and NIA 15K murine cDNA microarray studies identified hundreds of transcripts differentially expressed compared with wild type, over 30 of which were selected for further characterization by Northern blot analyses. Additionally, we performed in situ hybridization studies to localize 10 mRNAs, melanocyte-specific gene 1, amino acid transporter SN2, overexpressed and amplified in teratocarcinoma (Bcat1), Forkhead box protein FOXO1, 24p3, vascular cell adhesion molecule, epiregulin, Bcl2-like10, PC3B, and retinoblastoma binding protein 7. These 10 genes have expression patterns and postulated functions suggesting that they mediate important processes in the physiology and pathology of ovarian and adrenal tissue. | Introduction |
|---|
|
|
|---|
Completion of follicular development depends on the pituitary gonadotropins, FSH and LH (3). FSH and LH are glycoprotein
:ß heterodimers that share a common
-subunit in noncovalent association with a unique ß-subunit. Knockout mice homozygous for a null allele at the FSHß-encoding locus (Fshb-/-) exhibit female infertility due to a block in follicle development before the climax of granulosa cell proliferation and antrum formation that typify preovulatory follicles (4). Although these mice develop high levels of circulating LH, granulosa cells within these arrested follicles fail to become responsive to LH, do not express markers associated with luteinization, and are defective in the production of steroidogenic enzymes and peptide hormones (5). FSH receptor (Fshr) knockout females display a similar block in follicle development and also exhibit high levels of circulating LH (6, 7). It is believed that LH and LH-induced androgen synthesis in the Fshr model contribute to the development of ovarian sex cord stromal tumors, which are seen in 92% of these mice by 12 months of age. In these tumors there is a loss of granulosa cell proliferation control programs accompanied by an up-regulation of Sertoli cell markers, Mullerian inhibiting substance, and the GATA4 transcription factor (8).
Gonadal peptide hormones, inhibins and activins, normally elaborated by granulosa cells in females and Sertoli cells in males, comprise an endocrine feedback loop to the pituitary. Inhibins are
:ß heterodimers (
:ßA;
:ßB), and activins are ß:ß dimers (ßA:ßA; ßB:ßB; ßA:ßB), named for their effects, which inhibit and activate, respectively, pituitary production of FSH. Knockout mice with a deletion of the inhibin
-subunit locus, and thus deficient in inhibins, develop mixed granulosa/Sertoli cell tumors in the ovaries and testes (9). The progression of tumor development in these mice is associated with a cachexia-like wasting syndrome and high serum concentrations of the FSH and activins (9, 10). Together FSH and LH are required for the development of tumors in these mice (11), whereas FSH and sex steroids are important regulators of the process (Refs. 12 and 13 and our unpublished data). Castrated inhibin
knockout mice develop steroidogenic tumors of the adrenal cortex with similar histological features (10).
Granulosa cells normally up-regulate LH receptor expression in large preovulatory follicles and terminally differentiate or luteinize in response to the surge of circulating LH at ovulation (14, 15). Together with luteinizing follicular thecal cells, these cells form the corpus luteum and express factors critical to the maintenance of early pregnancy. Although this physiological LH exposure induces cell cycle withdrawal in the granulosa population and increased expression of cell cycle inhibitors, female LH-overexpressing transgenic mice develop tumors of the granulosa cell lineage in some genetic backgrounds (16, 17, 18). In these mice, a chimeric bovine LHß-human chorionic gonadotropin (hCG) transgene is driven in a tissue-specific manner to achieve a 10-fold increase in circulating LH immunoreactivity. The tumor phenotype in this model can be averted by administrating ovulatory doses of hCG, diverting granulosa cells to instead form luteomas that have a distinct histological appearance and transcriptome profile compared with the granulosa cell tumors (19).
All of the conditions modeled in the lines of transgenic mice considered above have analogies in human medicine (3). Women who are homozygous or compound heterozygous for inactivating FSHß ligand (frameshift/truncation and missense) or FSH receptor (missense) mutations exhibit normal preantral follicle development, but antral stage follicles capable of ovulation do not form (20). Granulosa cell tumors comprise nearly 10% of ovarian cancer cases in women, and they exhibit propensity for invasiveness, metastasis, and recurrence. Gain of function mutations of the LH receptor have been described in families with autosomal dominant precocious puberty (21), and a nonconservative amino acid change in inhibin
(Ala257Thr) has been associated with several cases of premature ovarian failure with elevated serum gonadotropins (22). However, to date neither activating mutations in gonadotropin receptors nor loss of inhibin function have been associated with human sex cord stromal tumor development.
To better appreciate the molecular mechanisms underlying ovarian and adrenal responses to gonadotropins and inhibins and the biological bases for the phenotypes of these transgenic models, we have undertaken a series of mRNA expression profiling experiments using Affymetrix 11K oligonucleotide array GeneChips and NIA 15K cDNA microarray slides. These approaches have revealed several novel candidate genes with potentially pivotal roles in ovarian and adrenal function and pathophysiology.
| Materials and Methods |
|---|
|
|
|---|
knockout (Inha-/-) mice and genotyping protocols have been described previously (4, 9, 23).
For RNA extraction, ovaries and adrenal glands were collected from 1316 wild-type C57BL/6J/129S6/SvEv randomly cycling females, 23 months old, and pooled unless otherwise indicated. Ovarian tumors were collected from individual inhibin
knockout females and were not pooled. These mice were killed at approximately 8 wk of age when demonstrating signs of a cachexia-like wasting syndrome associated with tumor development (9). To study adrenal tumor gene expression, castrated inhibin
null females were killed upon weight loss at 2235 wk of age (10). Tissues were immediately homogenized in RNA-STAT (Leedo Medical Laboratories, Houston, TX) and processed in accordance with the manufacturers protocol or were stored in RNAlater (Ambion, Austin, TX) before RNA extraction. For in situ hybridization experiments, ovaries, testes, or ovarian tumors were collected from mice, 23 months old, and immediately submerged in 4% paraformaldehyde fixative before processing and wax embedding.
Affymetrix gene chip mRNA expression analyses
Total RNA was used as a template for cDNA synthesis and biotinylated antisense cRNA probe preparation as described by the manufacturers of the SuperScript System kit (Invitrogen, Carlsbad, CA) and the ENZO BioArray HighYield RNA labeling kit (Enzo Diagnostics, Farmingdale, NY). Unincorporated nucleotides were removed from the riboprobe preparation using the RNeasy Mini kit (Qiagen, Valencia, CA). The integrity of the riboprobe was checked by gel electrophoresis.
The Affymetrix murine 11K oligonucleotide array set (Mu11KsubA, Mu11KsubB) was hybridized, washed, and scanned using Affymetrix equipment and protocols (Affymetrix, Santa Clara, CA). Data analysis was performed using Affymetrix software with the decision matrix default parameters recommended. For the purposes of this manuscript, only 3-fold or greater comparative changes in mRNA expression were considered significant and included in Tables 24![]()
![]()
below. More information detailing our protocols and data analyses and a listing of the genes surveyed in this analysis are available at http://www.affymetrix.com.
|
|
|
knockout female; these RNA samples were identical to three of those used for the Affymetrix analysis.
Wild-type (Cy3) and mutant (Cy5) probes were hybridized to the NIA 15K cDNA clone set (http://lgsun.grc.nia.nih.gov/cDNA/15k.html) at the Population Center Gene Array Facility (http://depts.washington.edu/popctrma/index.shtml). After washing, slides were scanned using a GenePix scanner (Axon Instruments, Union City, CA), and returned values were analyzed with GeneSpring software (Silicon Genetics, Redwood City, CA). For the purposes of this manuscript, only changes in mRNA expression more than 2-fold (ratio of medians, >2 or <0.5) with less than 20% difference between the ratio of the medians, median of ratios, and the normalized ratio of medians were considered significant and included in Tables 5
and 6
below.
|
|
knockout mice were reverse transcribed to cDNAs using the SuperScript RT kit (Invitrogen, Carlsbad, CA). The coding sequences of select genes were used to find primers 1922 bp in length with predicted Tms of 60 C using the Cprimer program (http://iubio. bio.indiana.edu/soft/molbio/mac/). Sequences were then amplified by PCR using an annealing temperature equivalent to or 5 C lower than the lower predicted primer Tm, ligated into pGEM(R)-T or T-easy vector (Promega, Madison, WI), and checked by bidirectional sequencing before being used as templates to synthesize Northern blot and in situ hybridization probes. The sequences amplified are listed in Table 1
|
For in situ hybridization, [35S]UTP antisense and sense probes were transcribed from the template plasmids (Table 1
) using T7 and SP6 polymerases (Promega, Madison, WI). Paraffin-embedded samples were cut into 5-µm-thick sections, which were dewaxed, pretreated, hybridized, and washed as previously detailed (27). Probe hybridization was detected by autoradiography using the NTB-2 emulsion (Eastman Kodak, Rochester, NY). Slides were developed as described by the manufacturer and were counterstained with hematoxylin.
Serum testosterone analysis
Ten randomly cycling, wild-type and Fshb-/- females, 34 months old, were anesthetized by isoflurane inhalation (Abbott Laboratories, North Chicago, IL), and blood was recovered by closed cardiac puncture. Serum was separated in Microtainer tubes (BD Biosciences, Franklin Lakes, NJ) and stored at -80 C before analysis. Testosterone measurements were made at the University of Virginia Ligand Core Facility (Specialized Cooperative Centers Program in Reproduction Research) and are reported as the average ± SE.
| Results and Discussion |
|---|
|
|
|---|
Profiling ovarian mRNA expression changes in FSHß knockout mice treated with pregnant mare serum gonadotropin (PMSG)
To identify potential gonadotropin-responsive gene products, we collected ovaries from 3-wk-old Fshb null females 8 h after administering 7.5 IU PMSG, an FSH receptor agonist that functions well regardless of the presence of endogenous FSH (4). Riboprobe prepared from total RNA was then hybridized to Affymetrix GeneChip oligonucleotide arrays, and the hybridization signals were compared with results obtained using untreated 3-wk-old Fshb-/- mice. The results of this comparison are shown in Table 3
. Several of the identified changes in gene expression are consistent with our previous results (Table 2
), including the values obtained for 24p3, inhibin/activin ßB, prostaglandin D synthetase, P450 aromatase, scavenger receptor class B type I, and steroidogenic acute regulatory protein. Thus, in select cases, aberrant gene expression in the Fshb knockout ovary can be readily reversed with PMSG administration.
Although PMSG is chiefly used to mimic FSH function in superovulatory regimes in rodents, it also has been demonstrated to activate the LH receptor (28, 29). Thus, several genes identified in this screen may be up-regulated as part of a periovulatory sequence in the ovary, as has been reported for epiregulin and steroidogenic acute regulatory protein, which are also induced after administration of the LH analog, hCG, to wild-type mice (30). Consistent with this, we found that up-regulation of epiregulin (Ereg), lymphocyte antigen 6 complex (Ly6d), and mouse 3T3 double minute 1 were recapitulated in wild-type mice only when PMSG was followed by treatment with hCG (see below).
Affymetrix profiling of mRNA expression in tumors of inhibin
knockout and LHß-overexpressing mice
Because it is feasible to recover large amounts of total RNA from the ovarian tumors that develop in the inhibin
knockout and LHß-overexpressing mouse models, and because it seemed conceivable that variation in gene expression patterns would probably exist between tumor samples and even within different regions of an individual sample, we performed Affymetrix analyses on individual tumors. A comparison analysis of one of the tumor samples recovered from an inhibin
knockout demonstrated hundreds of mRNA species with differential expression (>3-fold) compared with wild-type ovary (240 transcripts identified). The serial addition of a second inhibin
knockout tumor sample (262 transcripts identified), allowed us to compile a list of 118 transcripts with common directions of regulation in both samples. Table 4
lists the common expression changes of 74 known genes shared between the two inhibin
knockout tumor samples (excludes EST and duplicate sequences). We annotated this table by adding expression changes observed in the tumor transcriptome of an LHß-overexpressing mouse.
After gonadectomy, inhibin
knockout mice develop sex steroidogenic tumors of the adrenal cortex. Although these adrenal tumors do not result from metastases, they bear histological resemblance to granulosa cell tumors and perhaps share a similar means of pathogenesis (10). We used the Affymetrix system to assess gene expression in a unilateral adrenal tumor recovered from a 22-wk-old Inha-/- female after ovariectomy at 5 wk of age. A column in Table 4
indicates those gene products identified as present (P) vs. those identified as absent (A) in the absolute analysis of this sample. We speculate that gene products expressed in both gonadal and adrenal tumor samples have the potential to participate in transformation.
cDNA microarray analyses
To confirm the results of our Affymetrix studies and to identify additional changes in ovarian gene expression in FSHß and inhibin
knockouts, samples of ovarian RNA from 2- to 3-month-old, wild-type and Fshb-/- mice and ovarian tumor RNA from an inhibin
knockout female were studied by cDNA microarray analyses. Total RNA was used to transcribe fluorescently labeled cDNA, which was then hybridized to the NIA 15K murine microarray glass slides. Hybridization signals were scanned, and Cy5(mutant)/Cy3(wild-type) intensity ratios were compared. As expected, the majority of gene products and all loading controls demonstrated essentially equivalent hybridization signals and a Cy5/Cy3 ratio of approximately 1.0. For the purposes of this report, comparisons yielding fold changes with absolute values less than 2 as well as data obtained for the expression of ESTs have been omitted. The results of this analysis are listed in Tables 5
and 6
.
Comparing quantitative technologies: array analyses vs. Northern blot phosphorimaging
The array analyses described above provided us with quantitative data reflecting changes in the expression of many mRNA species in the ovaries of FSHß and inhibin
knockout mice vs. wild-type mice or in young FSHß knockout mice treated with PMSG. To better appreciate which of these findings were reproducible by standard molecular techniques and to what extent such analyses would provide comparable fold change values, we selected 32 genes for Northern blot expression analyses with phosphorimaging quantitation (26). The fold change values of the Affymetrix analyses with streptavidin conjugate amplification and/or of the cDNA microarray analysis as well as the Northern blot phosphorimaging fold changes are depicted in Table 7
, and autoradiographs are shown in Figs. 14![]()
![]()
![]()
. Adult wild-type ovary, Fshb-/- ovary, and wild-type adrenal gland RNA samples were extracted from pooled tissues. Ovarian and adrenal tumors from inhibin
knockout mice were not pooled. To further investigate gene products identified as up-regulated by gene chip analysis after PMSG treatment, we performed Northern blot analyses on RNA of pooled ovaries of 3- to 4-wk-old mice both 46 h after PMSG treatment (7.5 IU) and after 5 h of subsequent hCG treatment (5 IU). Although the majority of high throughput gene expression profiling findings were verifiable, some results were not readily reproduced in related samples or were not reproduced in the tissue RNA sample first profiled. This underscores the need to interpret array data with caution even when similar profiling results are obtained by related experimental approaches or replicate sampling.
|
|
|
|
|
|
The amino acid transporter SN2 was identified in these analyses as being up-regulated in ovaries and ovarian tumors of Fshb-/- and Inha-/- mice, respectively (Tables 5
and 6
). These findings were confirmed by Northern blot analyses, which also demonstrated up-regulation of SN2 mRNA expression in adrenal tumors of Inha-/- mice (Figs. 1
, 3
, and 4
). In situ hybridization showed that ovarian SN2 mRNA is localized exclusively to granulosa cells of growing follicles and is not detectable in corpora lutea (Fig. 5
, C and D). Thus, the increased level of ovarian SN2 mRNA in Fshb null mice reflects in part the relative enrichment of nonluteinized granulosa cells rather than direct gene regulation. SN2 is a neutral amino acid transporter that is known to be highly expressed in many tissues in rats and may be important in tumor development, because it makes glutamine available for catabolism or ammonia export (34, 35). Based on this, it seems that ovarian SN2 is transcribed in an FSH-independent manner and encodes a protein with potential roles in regulating or supporting the proliferation of granulosa cells.
Amino acid metabolism is regulated within cells by specific processing enzymes, including the branched chain aminotransferases (BCATs), which facilitate the first catabolic steps in the processing of leucine, isoleucine, and valine. One BCAT-encoding gene, Bcat1 or ECA39, is tightly linked to the Kras2 oncogene in both the mouse and human genomes. It was originally identified as being overexpressed and amplified in teratocarcinomas (36) and was further described as being a direct transcriptional target of the c-myc oncogene (37, 38). Bcat1 mRNA was identified in our microarray studies as being up-regulated in ovarian tumors, although by Northern blot analysis we found that this was true in adrenal tumors and was not readily reproducible in ovarian tumors (Table 5
and Figs. 3
and 4
). Despite these findings, which suggest that BCAT1 promotes cell proliferation, our in situ hybridization shows that Bcat1 mRNA is highly expressed in terminally differentiated luteal cells (Fig. 5
, E and F). There is other evidence that BCAT1 may function dichotomously in modulating cell cycle progression, as loss of the yeast ortholog leads to faster cell growth and a shortened G1 phase (39). Although their specific functions remain unclear, our identification of SN2 and BCAT1 calls attention to the pathways regulating granulosa cell amino acid import and metabolism and how they may participate as determinants of granulosa cell division and differentiation.
The transcription factor FOXO1 (FKHR) is classified as a Forkhead family member based on its three
-helix, two loop DNA-binding domain. There is accumulating evidence that these Forkhead family members are important regulators of cell death and cell cycle progression. Our observations of Foxo1 mRNA expression in oocytes of early follicle stages and in the granulosa cells of secondary and preovulatory follicles, but not corpora lutea (Fig. 5
, G and H), corroborate findings recently reported by Richards et al. (40). Based on granulosa cell culture studies, this group reported that FOXO1 function seems to be antagonized at the transcriptional and posttranslational levels by FSH, which down-regulates the expression of Foxo1 mRNA and causes the phosphorylation and nuclear to cytoplasmic sequestration of the FOXO1 protein. Consistent with this premise, we observed high levels of Foxo1 mRNA expression in granulosa cells of FSHß knockout ovaries, although Foxo1 was not identified in our gene chip analyses for the subtle increase in expression seen in the absence of FSH (Fig. 1
), which we ascribe to be largely due to enrichment of this cell population rather than direct gene derepression. Rather, Foxo1 mRNA was called to our attention because of its high level of expression in granulosa cell tumors in inhibin
knockout females despite high levels of circulating FSH, and this was confirmed by Northern blot analysis and in situ hybridization (Fig. 4
and data not shown).
Interstitial cell effects of the gonadotropins
Just as these studies identified genes expressed in granulosa cells, as expected, a surprising number of interstitial cell gene products were brought to our attention. Rather than representing direct targets of FSH or PMSG, we hypothesize that transcriptional control of these genes reflects interstitial cell responses to elevated levels of LH in Fshb-/- mice or occurs in response to steroid production in the granulosa cell compartment. Two mRNAs encoding 24p3 and vascular cell adhesion molecule 1 (VCAM1), identified in the Affymetrix array analysis as dramatically up-regulated in the ovaries of FSHß-deficient mice, were studied by in situ hybridization analyses. Both transcripts were highly expressed in the ovarian interstitial cell population in Fshb-/- mice, and no signal was detected in granulosa cells or oocytes of ovarian follicles (Fig. 5
, I and J). No expression was detectable in ovaries from age-matched, wild-type mice.
The 24p3 gene encodes a hydrophobic molecule transporter protein of the lipocalin family. Originally called an oncogene up-regulated in simian virus 40-transformed cells, 24p3 mRNA has been shown to be induced in response to glucocorticoid and retinoic acid (41) and is induced by steroid hormone signaling in the mouse uterus (42). 24p3 is also expressed in the epithelial cells of the mouse epididymis, and the protein associates with the acrosomes of caudal spermatozoa (43). We find by oligonucleotide array and Northern blot analyses that 24p3 is highly expressed in ovaries of FSHß knockout mice, although it is virtually undetectable in the wild-type ovary (Table 2
and Fig. 1
). In situ hybridization reveals that the mRNA species accumulates within the interstitial cell population between the arrested follicles (Fig. 5I
); no signal was detected in wild-type ovary. Notably, 24p3 protein exhibits closest homology to prostaglandin D synthase, which was also identified in the Affymetrix analysis as being up-regulated in the absence of FSH (Table 2
); both are down-regulated with PMSG administration (Table 3
).
VCAM1 (CD106) is a transmembrane protein belonging to the immunoglobulin superfamily, which is known to mediate cell-cell interactions by binding to
4 integrins and has been implicated in inflammatory processes (44). Other binding partners for VCAM1 have not been described, and the significance of its expression in the Fshb-/- ovary is not clear, since integrin
4 is termed absent by Affymetrix analysis in both wild-type ovary (45) and Fshb-/- ovary. Targeted disruption of Vcam1 in vivo results in lethality between d 1012 of mouse embryonic development due to defects in the placenta and heart, and this has precluded studies of its functions in other tissues (46, 47). Defects in VCAM1 expression are also seen in the embryonic lethal phenotypes of knockouts lacking SMAD1 and the Forkhead transcription factor FOXF1, although whether these changes reflect direct Vcam1 gene regulation is not clear (48, 49). Like 24p3, in situ hybridization reveals VCAM1 mRNA within the interstitial cell population of Fshb-/- ovaries (Fig. 5J
).
Two additional interstitial cell expressed genes, Ereg and Ly6d, were identified in these studies as being up-regulated in Fshb-/- mice given PMSG. Epiregulin is a member of the epidermal growth factor family and may have effects that either inhibit or promote cell cycle progression depending on its context. Ereg has been identified previously as a granulosa cell expressed gene that is responsive to gonadotropins, and it has been speculated that epiregulin may be involved in enhancing cyclooxygenase 2 mRNA expression during the periovulatory period (30). Our in situ hybridization confirms that epiregulin transcript is expressed in granulosa cells of preovulatory follicles after hCG administration and also demonstrates high levels of expression in the thecal/interstitial cell population (Fig. 5
, K and L). Based on the intensity of the in situ hybridization signal, Ereg does not seem to be differentially regulated between mural and cumulus granulosa cell populations. Ereg mRNA expression is not detectable in corpora lutea, indicating that the transcript is cleared upon ovulation and the differentiation of luteal cells.
LY6D (or ThB) is a cell surface glycoprotein that is involved in cell-cell interactions (50). Although its normal functions in vivo are not well understood, differential expression of Ly6d in tumor cell populations seems to be a determinant of how cancers interact with endothelial cells and blood vessel precursors to ultimately regulate the potential for tumor invasiveness, metastasis, and angiogenesis (51). By Northern blot and in situ hybridization analysis, Ly6d mRNA is weakly detected in the interstitial cells of wild-type ovaries after treatment with hCG (Fig. 2
and data not shown). Ly6d is highly expressed in the outermost layer of the adrenal cortex of wild-type mice, and this expression is lost during adrenal tumorigenesis in the inhibin
knockout mouse model (Fig. 4
and data not shown).
Androgen biosynthesis in the FSHß knockout mouse
Affymetrix oligonucleotide hybridization analyses of Fshb-/- ovaries indicated aberrant expression levels of a number of mRNAs encoding proteins associated with steroidogenic pathways, some of which we had previously studied (5). Additionally, this study identified two androgen-responsive genes encoding seminal vesicle F protein and seminal vesicle secretory protein as being up-regulated in the ovaries of Fshb-/- mice. Given the high levels of circulating LH in these mice (4) and their relative lack of aromatase mRNA for converting ovarian androgens to estrogens (5), we hypothesized that the lack of FSH would result in elevations in circulating androgens. To test this, we measured circulating testosterone in age-matched, adult wild-type and FSHß knockout mice and found that knockouts have a 4-fold increase in circulating testosterone over wild-type mice (0.94 ± 0.24 vs. 0.21 ± 0.07 ng/ml). Corroborating these findings, elevated testosterone levels have been recently reported for FSH receptor knockout mice (52). Thus, these studies indicate a potential for elevated androgens to contribute to the etiology of FSH deficiency syndromes in women.
Cell cycle control in oocytes: potential roles for Bcl210, PC3B, and Rbb7
In situ hybridization analyses demonstrated oocyte expression of mRNAs encoding the cell cycle regulators BCL2-like-10 (also known as BCLB, BOO, and DIVA), PC3B (also known as B cell translocation gene 4), and retinoblastoma binding protein 7 (RBBP7; also known as mRbAp46). The first two transcripts were identified as being down-regulated in tumor samples from inhibin
knockout mice; the last was identified as up-regulated in FSHß knockout ovaries. Their relative expression in each case reflects changes in tissue composition between these knockout models; there are no oocytes in advanced tumors, and there is a relative increase in oocytes in Fshb-/- ovaries (see below). In oocytes, these gene products may function to inhibit apoptotic programs and maintain meiotic arrest.
Bcl2l10 encodes an apoptotic protease-activating factor 1-binding protein that prevents apoptotic protease-activating factor 1 from activating downstream caspase enzymes (53, 54). Bcl2l10 mRNA is highly expressed in the mammalian ovary and epididymis based on our analyses and previous descriptions of the gene (53, 54). The observation of ovarian expression led to the speculation that BCL2-like-10 is a regulator of apoptotic programs that mediate follicle atresia. Our in situ hybridization analysis reveals that Bcl2l10 mRNA is expressed preferentially in oocytes, and that this expression persists in oocytes of primary to preovulatory stage follicles (Fig. 5M
). This expression was also clearly visible in ovaries of Gdf9-/- and Fshb-/- mice (data not shown). Notably, there is no expression detectable in oocytes of primordial follicles; thus, BCL2-like-10 is more likely to influence oocyte apoptosis after follicle recruitment than the survival of primordial germ cells. Although a recent knockout study indicates that BCL2-like-10 is not alone necessary for ovarian development or function (55), it is possible that double mutant studies will indicate important, although redundant, in vivo roles for this protein. For example, BCL2-Like-10 may cooperate with BCLX or BCL2, both of which are important apoptotic inhibitors in mouse oocytes based on loss of function and gain of function mouse models (56, 57, 58).
Our in situ hybridization studies confirm the oocyte expression of PC3B, which was previously reported based on EST database and RT-PCR analyses (59). PC3B is a member of the PC3/B cell translocation gene family identified by EST database sequence analyses (59). Related gene products can inhibit cyclin D1 expression in a retinoblastoma-dependent manner and therefore limit the propensity of a cell to proceed from G1 to S phase (60). Similarly, PC3B is able to induce G1 arrest when expressed in NIH-3T3 cells (59). In addition to being expressed in oocytes, PC3B mRNA has been detected in testis and differentiating olfactory epithelium (59). We localized PC3B mRNA to growing oocytes from the primary stage throughout the remainder of follicle development (Fig. 5N
and data not shown). There is no detectable expression in other ovarian cell types.
Finally, a third gene, Rbbp7, was identified in our studies as being preferentially expressed in mouse oocytes. Like Bcl2l10 and PC3B mRNAs, Rbbp7 transcript is detectable in the cytoplasm of oocytes at the primary follicle stage and in oocytes of more advanced ovarian follicles. It is also highly expressed in the arrested primary follicles of Gdf9 null ovaries (Fig. 5
, O and P). RBBP7 was originally identified from HeLa cell lysates based on its binding affinity for RB, and ribonuclease protection assays indicated that it is expressed widely in mouse tissues, with highest expression in testis and ovary/uterus (61). RBBP7 (RbAp46) is closely related to RBBP4 (RbAp48), and both homologs are highly conserved between mouse and humans. Although this is the first report of Rbbp7 expression in mammalian oocytes, the Xenopus RbAp46/48 ortholog is known to complex with a 57-kDa histone deacetylase expressed only during oogenesis and early cleavage stage embryos (62). The functions of RBBP7 have not yet been identified in vivo, but may prove important in mediating chromatin configuration and accessibility to transcription factors in oocytes and early embryos.
Implications for future expression profiling endeavors in reproductive tissues
In a previous study, we used Affymetrix GeneChip analysis to define the direct downstream targets of GDF9 in granulosa cells (63). In these experiments we isolated granulosa cells from PMSG-treated females, placed them in culture, and maintained them for 5 h in the presence or absence of recombinant GDF9. These studies examined an isolated population of cells at one time point, the only difference between the two samples being GDF9 stimulation. We were able to confirm known targets of GDF9 (e.g. cyclooxygenase 2 and hyaluronan synthase 2) (64) and also identify novel downstream targets of GDF9 (e.g. pentraxin 3 and TNF-induced protein 6). In the present study, tissues composed of multiple cell lineages were used in their entirety, and several genes were identified as abnormally expressed in the absence of either FSH or inhibins. Although some of these genes may represent direct targets of FSH or participants in the tumor process, in situ hybridization analyses demonstrated that many transcripts were instead increased or decreased in representation because of structural alterations in the ovaries. Specifically, FSHß knockout ovaries have a block in folliculogenesis before antral follicle development, resulting in a relative increase in the number of secondary follicles and oocytes compared with wild-type adult ovaries, in which much of the ovary is filled with functional or regressing corpora lutea. The result is that genes expressed in oocytes, such as Zp1 and Rbbp7, are termed up-regulated in the absence of FSH. Similarly, markers of luteinization, including calcyclin and LH receptor, are termed down-regulated, as luteinization is precluded without FSH (5). Likewise, comparing ovarian tumors of inhibin
knockout mice and wild-type ovaries yields a host of down-regulated gene products specific to oocytes (e.g. Gdf9, Zp2, Zp3, Bcl2l10, and PC3B) and luteinized cells (e.g. calcyclin and LH receptor), which are not found in ovarian tumors. Thus, whether similar studies are being planned to investigate in vivo gene expression in the ovary or in other tissues, it is critical to perform in situ hybridizations to determine the principle cell population expressing a gene of interest. In addition, short-term in vitro experiments with single cell types and single ligands can be expected to simplify the results obtained and identify target genes more readily.
The data in this study have been compiled using two different expression profiling methodologies, Affymetrix oligonucleotide arrays and cDNA microarrays, and we find the results of these analyses to be complementary in that each has revealed different transcripts to be aberrantly expressed in analyses of the same RNA sample. Reasons for this include differences in the sequences represented on each array and may also extend to different sources of technical variability, as different strategies are employed by each to synthesize riboprobe, different conditions are used for the hybridization, and different quantifications and mathematical analyses to identify regulated transcripts are used. For each type of analysis, we did not perform replicate analyses on the same RNA sample, and data indicate that for genes that are highly or even moderately expressed provide good reproducibility in such analyses (65). We followed up specific changes by Northern blot analysis rather than combining multiple array studies for samples that were relatively homogeneous, such as ovarian samples pooled from a group of wild-type mice or FSHß knockout mice. We relied on the fact that these samples were pooled from a large group of genetically identical and age-matched mice to minimize changes in gene expression in response to stress effects (known to vary more widely in comparative studies of nonovarian tissues) (65) or other causes of variability from individual to individual. This concern presents another rationale for single cell analyses, as different samples of cultured cells demonstrate variability in gene chip studies comparable to those attributable to technical variability from chip to chip (65). We adopted a different strategy in analyzing tumor samples from inhibin
knockout mice, profiling two ovarian tumors, and subsequently annotating our results with gene changes shared by analysis of an adrenal tumor from an inhibin
knockout mouse and a study of an ovarian tumor from an LHß-overexpressing mouse. Using two tumor samples enabled us to narrow a list of approximately 250 gene products defined as regulated in either sample to a list half this size of gene expression changes shared by both samples. Consistent with our hypothesis that there are commonalities in the pathophysiology of granulosa cell tumors of different mouse models and between these ovarian tumors and the adrenal tumors that occur in the absence of inhibins, we find several examples of shared gene expression changes. Together, these studies suggest several avenues for collecting and verifying expression profiling data.
In conclusion, this work has identified a myriad of gene products, many of which have not been previously characterized in the mammalian ovary or adrenal gland or implicated in ovarian tumorigenesis. We focused more in-depth in situ hybridization analyses of gene expression on transcripts encoding putative cell cycle regulators and have uncovered mRNAs preferentially expressed in every ovarian cell type. Only a few similar studies have been undertaken and reported to date, and none with a more extensive series of follow-up expression studies (19, 66, 67). Together, these analyses will prompt future investigations into how specific gene products contribute to normal and aberrant ovarian functioning. Some of the differential gene expressions uncovered may prove important in the endocrinopathies modeled by our transgenic mice, whereas others, perhaps equally relevant to our understanding of ovarian biology, reflect the results of alterations in tissue composition between these models.
| Acknowledgments |
|---|
| Footnotes |
|---|
Abbreviations: Bcat, Branched chain aminotransferase 1; Ereg, epiregulin; EST, expressed sequence tag; Gapd, glyceraldehyde 3-phosphate dehydrogenase; hCG, human chorionic gonadotropin; Ly6d, lymphocyte antigen 6 complex; Msg1, melanocyte-specific gene 1; PMSG, pregnant mare serum gonadotropin; RBBP7, retinoblastoma binding protein 7; VCAM1, vascular cell adhesion molecule 1.
Received April 16, 2003.
Accepted for publication June 17, 2003.
| References |
|---|
|
|
|---|
-Inhibin is a tumor-suppressor gene with gonadal specificity in mice. Nature 360:313319[CrossRef][Medline]