Disruption of postnatal folliculogenesis and development of ovarian tumor in a mouse model with aberrant transforming growth factor beta signaling

Background Transforming growth factor beta (TGFB) superfamily signaling is implicated in the development of sex cord-stromal tumors, a category of poorly defined gonadal tumors. The aim of this study was to determine potential effects of dysregulated TGFB signaling in the ovary using Cre recombinase driven by growth differentiation factor 9 (Gdf9) promoter known to be expressed in oocytes. Methods A mouse model containing constitutively active TGFBR1 (TGFBR1CA) using Gdf9-iCre (termed TGFBR1-CAG9Cre) was generated. Hematoxylin and eosin (H & E) staining, follicle counting, and immunohistochemistry and immunofluorescence analyses using antibodies directed to Ki67, forkhead box L2 (FOXL2), forkhead box O1 (FOXO1), inhibin alpha (INHA), and SRY (sex determining region Y)-box 9 were performed to determine the characteristics of the TGFBR1-CAG9Cre ovary. Terminal deoxynucleotidyl transferase (TdT) labeling of 3’-OH ends of DNA fragments, real-time PCR, and western blotting were used to examine apoptosis, select gene expression, and TGFBR1 activation. RNAscope in situ hybridization was used to localize the expression of GLI-Kruppel family member GLI1 (Gli1) in ovarian tumor tissues. Results TGFBR1-CAG9Cre females were sterile. Sustained activation of TGFBR1 led to altered granulosa cell proliferation evidenced by high expression of Ki67. At an early age, these mice demonstrated follicular defects and development of ovarian granulosa cell tumors, which were immunoreactive for granulosa cell markers including FOXL2, FOXO1, and INHA. Further histochemical and molecular analyses provided evidence of overactivation of TGFBR1 in the granulosa cell compartment during ovarian pathogenesis in TGFBR1-CAG9Cre mice, along with upregulation of Gli1 and Gli2 and downregulation of Tgfbr3 in ovarian tumor tissues. Conclusions These results reinforce the role of constitutively active TGFBR1 in promoting ovarian tumorigenesis in mice. The mouse model created in this study may be further exploited to define the cellular and molecular mechanisms of TGFB/activin downstream signaling in granulosa cell tumor development. Future studies are needed to test whether activation of TGFB/activin signaling contributes to the development of human granulosa cell tumors. Electronic supplementary material The online version of this article (10.1186/s12958-017-0312-z) contains supplementary material, which is available to authorized users.


Background
Transforming growth factor beta (TGFB) superfamily members regulate a variety of biological processes in a contextually dependent manner via the interaction with membrane associated TGFBR1/TGFBR2 serine/ threonine kinase receptor complexes and downstream SMAD proteins [1]. A growing body of research documents critical roles of TGFB signaling in reproductive development, function, and dysfunction [2,3]. TGFB signaling is known for its tumor suppressive function, particularly in epithelial cancer cells [4,5]. However, TGFB signaling also promotes cancer progression through modulating tumor cell invasion and metastasis [6]. Recently, we have shown that TGFB pathway activation plays a crucial role in pancreatic tumor initiation, through its capacity to induce Acinar-to-Ductal metaplasia, providing a favourable environment for KRAS G12D -dependent carcinogenesis [7].
Disruption of genes encoding TGFB superfamily signaling mediators [e.g., inhibin alpha (Inha), bone morphogenetic protein (BMP)-responsive SMADs (i.e., Smad1/5), and BMP type 1 receptors (i.e., Bmpr1a and Bmpr1b)] results in the development of sex cordstromal tumors [8][9][10]. Studies have revealed tumorigenic function of activin-responsive SMAD3 and tumor suppressive roles of inhibins and BMPs in the gonad [8][9][10][11]. A cornerstone study has demonstrated that nearly all adult granulosa cell tumors (GCTs), a major subtype of sex cord-stromal tumors, harbor a somatic mutation of forkhead box L2 (FOXL2; 402C→G; C134W) [12]. However, the specific role of FOXL2 mutation and its associated tumorigenic signals in GCTs remain to be clarified. Therefore, it is imperative to define the mechanism underlying the initiation and progression of ovarian GCTs, with an ultimate goal of developing new therapeutic strategies for this type of poorly characterized tumor.
As the functional unit of the ovary, a follicle consists of granulosa and theca somatic cells and a germ cell, the oocyte. Through secretion of paracrine growth factors, the oocyte regulates multiple functions of granulosa cells including, but not limited to, apoptosis, proliferation and differentiation, steroidogenesis, and metabolism [13,14]. It is noteworthy that dysregulation of phosphoinositide-3-kinase (PI3K)-AKT signaling in the mouse oocyte impairs normal ovarian function and causes ovarian GCT formation [15]. In an early report, we have shown that constitutive activation of TGFB signaling in mouse granulosa cells promotes ovarian tumorigenesis [16]. To further explore the tumorigenic function of dysregulated TGFB signaling in distinct cellular compartments of the ovary, we created a mouse model that harbors constitutively active TGFBR1 using growth differentiation factor 9 (Gdf9)-iCre known to be expressed in the oocyte.

Histology and follicle counting
Hematoxylin and eosin (H & E) staining was conducted using standard procedures. To perform follicle counting, ovaries were first serially sectioned (5 μm) and stained with periodic acid Schiff 's (PAS) and hematoxylin. Follicles were counted from every 5th section. Follicle classification was based on previously established criteria [23,24]. Primordial and primary follicles were counted, regardless of the status of the oocyte nucleus. For secondary follicles, only those with visible oocyte nucleus in the section were counted as previously described [23]. Since the oocyte diameter for primordial and primary follicles averages~14 μm, the same follicle could be counted every two sections. Hence, the estimated follicle numbers per ovary were calculated by multiplying the cumulative number of each ovary by 5 and then dividing by 2 as described [23].

In situ apoptosis analysis
A commercially available In Situ Apoptosis Detection Kit (Abcam) was used to detect apoptosis in ovarian samples from control and TGFBR1-CA G9Cre mice at postnatal day 3 (PD3), PD5, PD7, PD12, and PD21. Briefly, paraffin sections were deparaffinized in xylene and rehydrated in a graded alcohol series before the treatment with Proteinase K. Then 3% H 2 O 2 was added to inactivate endogenous peroxidases. Apoptotic cells were labeled with terminal deoxynucleotidyl transferase (TdT) that catalyzes the addition of biotin-labeled deoxynucleotides, followed by incubation with streptavidinhorseradish peroxidase (HRP) conjugate. Positive controls where tissue sections were treated with DNase I and negative controls where TdT was substituted with water were included. The signal was detected using 3,3′diaminobenzidine (DAB) substrate and sections counterstained with Methyl Green.

X-gal staining
As described [25], ovarian samples were first fixed in 2% paraformaldehyde-0.2% glutaraldehyde in a 0.1 M phosphate buffer solution (pH 7.4), and then stained with X-gal staining buffer containing 1 mg/ml X-gal, 5 mM potassium ferricyanide, and 5 mM potassium ferrocyanide. After the staining procedure, samples were post-fixed with 10% neutral buffered formalin for histological analysis. Fast red (Vector Laboratories) was used to counterstain the nucleus.

RNAscope analysis
RNAscope 2.5 HD Reagent Kit and Gli1 probe (Catalog no. 311001) were purchased from Advanced Cell Diagnostics. Positive control peptidylprolyl isomerase B (Ppib) and negative control dihydrodipicolinate reductase (DapB) probes were included. RNAscope in situ hybridization analysis was performed using formalin-fixed paraffin-embedded sections based on the manufacturer instructions. Hybridization signals were developed and visualized using DAB. Sections were counterstained using hematoxylin and mounted with Permount Mounting Medium.

Statistical analyses
A two-tailed t-test (unpaired) was used to determine the difference of means between two groups. Data are mean ± standard error of the mean (s.e.m). Statistical significance was reported at * P < 0.05, * * P < 0.01, or * * * P < 0.001.

Results
Development of ovarian tumors in female mice harboring TGFBR1 CA conditional allele and Gdf9-iCre The TGFBR1 CA allele was engineered to contain three missense mutations to constitutively activate TGFBR1 and prevent the inhibitory effect of FK506 binding protein 1A (FKBP12) as previously described [29,30]. To determine the potential phenotype of TGFBR1-CA G9Cre mice, we first examined their fertility. During a 3-month fertility test, we found that TGFBR1-CA G9Cre females (n = 3) were sterile, in contrast to wild type mice (8.25 ± 0.26 pups/litter and 1.22 ± 0.07 litter/month; n = 6). To define the cause of infertility, we examined the morphology of the ovary. Remarkably, we found the development of gross ovarian tumors in TGFBR1-CA G9Cre mice at the age of 2-3 months. A summary of tumor development in mice at the age of 8-10 weeks is listed in Table 1. Histological analysis using H & E staining suggested the development of hemorrhagic GCTs in TGFBR1-CA G9Cre mice compared with controls ( Fig. 1a-d). The neoplastic cells appeared mitotically active ( Fig. 1e and f ). This observation was confirmed by high expression of Ki67, a cell proliferation marker, compared with controls ( Fig.  1g and h). Figure 1i shows ovarian tumors in TGFBR1-CA G9Cre mice at the age of 7 months.

Perturbed follicular development at an early stage
To study the potential effect of TGFBR1 activation on follicular development, we performed follicle counting using ovaries from control and TGFBR1-CA G9Cre mice at PD5 and PD7. The results did not reveal a significant difference in the number of primordial follicles between control and TGFBR1-CA G9Cre mice at PD5, although a reduction of primary follicles was detected in the TGFBR1-CA G9Cre mice (Fig. 2a). At PD7, the number of primordial follicles, primary follicles, and secondary follicles was reduced in TGFBR1-CA G9Cre ovaries compared with controls (Fig. 2b). Meanwhile, large follicles or histologically abnormal follicle-like structures were found in the ovaries of TGFBR1-CA G9Cre mice ( Fig. 2c  and d). Further double immunofluorescence staining of ACTA2 (a theca marker) and DDX4 (an oocyte marker) using PD12 ovaries independently confirmed these findings (Fig. 3a-f; dotted yellow lines). Ovarian histology of control and TGFBR1-CA G9Cre mice at the age of PD7 and PD21 is shown in Additional file 1: Figure S1. These results suggest that constitutive activation of TGFBR1 using Gdf9-iCre disrupts ovarian follicle development, and the reduced follicle numbers in TGFBR1-CA G9Cre mice may be associated with destruction/loss of follicles resulting from tumor initiation/development.

Molecular analysis of ovarian tumor type
To define the molecular characteristics of ovarian tumors in TGFBR1-CA G9Cre mice, we performed immunostaining to examine the expression of a granulosa cell lineage marker FOXL2 [31] and three other granulosa cell-expressed proteins, FOXO1, INHA, and AMH. DDX4, a germ cell marker, was also included. The localization of FOXL2 (Fig. 4a), INHA (Fig. 4c), FOXO1 (Fig. 4e), AMH (Fig. 4g), and DDX4 ( Fig. 4i) was detected in the granulosa cell or oocyte compartment of control ovaries, while ovarian tumor tissues from TGFBR1-CA G9Cre mice were immunoreactive with FOXL2 ( Fig. 4b), INHA (Fig. 4d), and FOXO1 (Fig. 4f ), supporting the development of GCTs in these mice. However, immunoreactive signals of AMH were low to undetectable in these tumors (Fig. 4h), consistent with the reported lack of AMH expression in GCTs resulting from conditional deletion of FOXO1/3 and PTEN [32]. The tumor nodules did not express DDX4 (Fig. 4j). Representative negative controls using rabbit and goat IgGs are respectively shown in Fig. 4k and l. Because GCTs can express Sertoli cell genes [32], we examined whether SOX9 was detectable in TGFBR1-CA G9Cre tumors. Our results showed heterogeneous expression of SOX9 in tumors from TGFBR1-CA G9Cre mice (Fig. 5d-i; arrows) versus controls (Fig. 5a-c). Taken together, these findings support the histological observation that sustained activation of TGFBR1 using Gdf9-iCre leads to the formation of sex cord-stromal tumors reminiscent of GCTs.

TGFBR1-CA G9Cre mice did not show increased apoptosis in the oocyte
It has been suggested that loss of oocytes during follicular development may alter the differentiation and cell fate of ovarian granulosa cells [33]. To determine whether TGFBR1 activation promotes oocyte apoptosis, we examined fragmentation/damage on DNA suggestive of apoptosis using an In situ Apoptosis Detection Kit and ovarian samples collected at PD3-PD21. To validate the experimental procedure, we showed that ovarian samples of positive controls treated with DNase I to generate fragmented DNA contained strong immunoreactive signals (Fig. 6a), while negative controls showed only background staining (Fig. 6b). The analysis did not show alterations of oocyte apoptosis in TGFBR1-CA G9Cre ovaries at PD3-PD12 ( Fig. 6c and d and data for PD5, PD7, and PD12 not shown). At PD21, however, signals were observed in ovarian somatic cells within some follicles in TGFBR1-CA G9Cre mice compared with controls ( Fig. 6e and f ). Thus, these results indicate that oocyte apoptosis is not a main contributing factor to the overall ovarian pathology observed in the TGFBR1-CA G9Cre mouse model.

Evidence of TGFBR1 activation in ovarian granulosa cells of TGFBR1-CA G9Cre mice
In line with the expression of Gdf9-iCre in primordial follicles from PD3 [17], TGFBR1 CA transcripts were highly expressed in the ovaries of TGFBR1-CA G9Cre mice (Fig. 7a), coinciding with increased expression of Smad7, a TGFB target gene (Fig. 7b). Moreover, the mRNA expression of Inha was increased and Zp3 reduced at PD7 (Fig. 7c and d), corroborating altered granulosa cell and oocyte properties. Somewhat unexpected, our further analysis showed that ovarian tumor tissues which barely contained oocytes (Fig. 4j) readily expressed TGFBR1 CA , along with increased levels of phospho-SMAD2/3 and reduced expression of HSD3B in comparison to controls (Fig. 7e). Activation of TGFB signaling was also examined using ovarian samples from 1-month-old mice (Additional file 2: Figure S2). These results suggest overactivation of TGFB signaling in GCT cells and alteration of cell differentiation.
The aforementioned findings raised the question on whether overactivation of TGFBR1 in granulosa cells conferred the development of ovarian tumors in TGFBR1-CA G9Cre mice. To verify the expression of Gdf9-iCre in female germ cells, we performed X-gal staining using ovaries from Rosa26/Gdf9-iCre mice and demonstrated predominant Cre activity in oocytes from the primordial stage (Fig. 7f-h). However, minor sporadic signals were observed in some follicles ( Fig. 7g;  arrows), and the ovarian surface also appeared to be stained (Fig. 7h). As a quality control, ovaries from Rosa26/Zp3-Cre mice were included, where X-gal staining was exclusively observed in growing oocytes (Additional file 3: Figure S3 A and B). Negative controls using ovaries from Rosa26 reporter mice alone did not show X-gal staining (Fig. 7i). Moreover, higher levels of  (a and b). Panel f is a higher magnification image for the boxed region in panel (e). g and h Immunohistochemical analysis of Ki67 using 8-week-old TGFBR1-CA G9Cre (g) and control (h) ovaries. Experiment was performed using ABC method, and signals were developed using NovaRED Peroxidase Substrate Kit. Sections were counterstained with hematoxylin. Scale bar is representatively depicted in (a) and equals 12.5 μm (f), 25 μm (e, g, and h), 50 μm (b), and 250 μm (a, c, and d). H & E staining and immunohistochemistry were conducted using 3-5 independent samples per group. i Gross ovarian tumor morphology of a 7-month-old mouse. Yellow arrows denote the ovary and ovarian tumors in the control and TGFBR1-CA G9Cre mice, respectively Gli1 and Gli2 transcription factors and lower levels of Tgfbr3 were found in TGFBR1-CA G9Cre ovaries versus controls by real-time PCR (Fig. 7j). RNAscope in situ hybridization analysis further revealed the strong localization of Gli1 mRNA to GCT tissues (Fig. 7l), in contrast to the restricted expression pattern of Gli1 mRNA to theca layers in control ovaries (Fig. 7k).
Representative positive and negative controls for the RNAscope experiment are respectively depicted in Fig. 7m and n.

Discussion
GCTs are the major type of sex cord-stromal tumors, and both genomic and genetic factors are involved in Fig. 2 Altered follicular development in TGFBR1-CA G9Cre mice. a and b Follicle counts of control and TGFBR1-CA G9Cre ovaries at PD5 (a) and PD7 (b). Data are mean ± s.e.m. n = 3. * P < 0.05. Ns, not significant. c and d Immunohistochemical localization of INHA in PD7 control and TGFBR1-CA G9Cre ovaries. Arrows indicate abnormal follicle structures. Experiment was performed using ABC method, and signals were developed using NovaRED Peroxidase Substrate Kit. Sections were counterstained with hematoxylin. Four independent samples per group were used for immunohistochemical analyses. Scale bar is representatively depicted in (c) and equals 50 μm (c and d) Fig. 3 Immunofluorescence analysis of follicular defects in TGFBR1-CA G9Cre mice. a-f Double immunofluorescence of DDX4 (red) and ACTA2 (green) using PD12 control (a-c) and TGFBR1-CA G9Cre (d-f) mice. Three independent samples per group were analyzed using immunohistochemistry and/or immunofluorescence. Note the presence of large follicles or abnormal follicle-like structures (dotted yellow lines) in the ovaries of TGFBR1-CA G9Cre mice in comparison with age-matched controls. Scale bar is representatively depicted in (a), and equals 100 μm (a-f) the pathogenesis of these tumors [34,35]. Human GCTs can be divided into adult and juvenile types; the former represents more than 90% of all GCTs and mainly occurs in perimenopausal or postmenopausal women around 50-54 years of age [36]. The Juvenile GCTs are relatively rare and mainly occur in prepubertal girls [37]. Although the 5-year survival rate in stage I patients is generally high, a poor prognosis occurs in patients with advanced stages of disease [38][39][40]. Of note, human adult GCTs have a high risk of recurrence [41], leading to resistance to chemotherapy and death. Because of the rarity of this category of tumors, mouse models are beneficial to study the etiopathology of these tumors.
Several reports using genetically modified mouse models highlight the importance of TGFB superfamily signaling components, such as INHA, SMAD3, SMAD1/5, and BMPR1A/BMPR1B, in sex cordstromal tumor development [8-11, 42, 43]. A causal link between the mRNA expression levels of TGFBR1 and the development of GCTs in humans has not been established [43]. Growing evidence supports that activation of TGFB/activin signaling may represent a driving force of GCT development in mice [44]. It has also been found that TGFB/activin signaling is active in human GCTs [32,45]. Notably, nearly all adult GCTs bear a somatic missense mutation of FOXL2 C→G (C134W) [12], which may alter activin/TGFB and BMP signaling activity and the proliferation and differentiation status of granulosa cells [46][47][48]. In an effort to define the role of aberrant activation of TGFB signaling in the pathogenesis of ovarian tumors, we utilized a mouse model harboring constitutively active TGFBR1. Our results corroborate the role of TGFB signaling in ovarian tumorigenesis evidenced by the finding that constitutive activation of TGFB signaling in ovarian somatic cells leads to the formation of ovarian malignancies that phenocopy GCTs in several perspectives [16]. Of note, sustained activation of TGFBR1 led to the phosphorylation of SMAD2/3, which are downstream signaling elements shared by TGFBs and activins. Therefore, our results do not discount the role of activin signaling in GCT development. Indeed, our early studies have shown that activin signaling is important in promoting ovarian sex cord-stromal tumor development [11,42]. Collectively, our approach to manipulate TGFBR1 activity and SMAD2/3 activation represents a valuable tool to study ovarian GCT development.
Oocyte-somatic cell communication is critical for follicular development. Through the secretion of growth factors, some of which are TGFB superfamily proteins, oocytes regulate folliculogenesis in a bi-directional , and DDX4 (i and j) are shown. k and l Negative controls using isotype-matched rabbit and goat IgGs. Experiment was performed using ABC method, and signals were developed using NovaRED Peroxidase Substrate Kit. Sections were counterstained with hematoxylin. Five independent samples per group were utilized in this analysis. Scale bar is representatively depicted in (a) and equals 50 μm (a-l) paradigm [13,14]. The involvement and potential role of oocyte-granulosa cell regulatory loop in sex cordstromal tumor development are poorly understood. A recent study has provided functional evidence that constitutive activation of PI3K-AKT signaling in the oocyte promotes GCT formation [15], suggesting that dysregulation of major growth regulatory pathways in the oocyte impacts granulosa cell growth and differentiation. Genetic evidence suggests that SMAD4dependent canonical TGFB signaling in oocytes is largely dispensable for female fertility [49]. However, it has not been determined whether unopposed TGFB signaling in the oocyte is detrimental to ovarian development and function. In the current study, ovarian tumors developed in mice harboring TGFBR1 CA and Gdf9-iCre. By performing histological and molecular analyses, we demonstrated that the ovarian neoplasms were sex cord-stromal tumors which expressed granulosa cell markers including FOXL2, INHA, and FOXO1. An elegant review by Pitman and colleagues suggests that premature loss of ovarian germ cells may promote the formation of ovarian neoplasms [33]. However, the tumor phenotypes appear to be associated with the type of mutations and the onset of oocyte loss [33]. In this study, the specific cause of follicle/oocyte reduction and its potential involvement in ovarian tumor development remain unknown. The observation that overactivation of TGFBR1 using Gdf9-iCre did not increase oocyte apoptosis suggests that the reduction of follicle/oocyte numbers in TGFBR1-CA G9Cre mice may be associated with the disruptive effect (e.g., impairment/destruction) of tumor formation and development.
The development of GCTs in TGFBR1-CA G9Cre mice begs the question of how overactivation of TGFBR1 links to ovarian tumor formation. As a first step, expression of Gdf9-iCre was verified in our model. Our X-gal staining showed predominant Cre activity in the oocyte. Somewhat unexpectedly, we also observed minor sporadic signals in some follicles. Of note, in TGFBR1-CA G9Cre mice, Gdf9-iCre was transmitted maternally due to the extremely low efficiency to generate TGFBR1 CA flox/+ ; Gdf9-iCre mice using male breeders (i.e., Gdf9-iCre males) to transmit Gdf9-iCre. In our breeding strategy, Gdf9-iCre females were crossed with TGFBR1 CA flox/+ male mice to produce TGFBR1 CA flox/+ ; Gdf9-iCre females. The reason for the inability of obtaining TGFBR1 CA flox/+ ; Gdf9-iCre mice using males to transmit Gdf9-iCre was not clear. Since constitutively active TGFBR1 in ovarian somatic cells is a strong driver of ovarian GCTs in mice [16], it was possible that even low to negligible Cre activity in a subset of granulosa cells might be permissive for TGFBR1 CA activation, leading to GCT development. Further supporting the contribution of overactivation of TGFBR1 in granulosa cells to ovarian tumorigenesis in the TGFBR1-CA G9Cre model, mice harboring TGFBR1 CA and Zp3-Cre, which is expressed in growing oocytes, but not non-growing oocytes of primordial follicles [19], did not develop ovarian tumors, regardless of the parental origin of Cre transmission (Y. Gao and Q. Li, unpublished observation). We previously found that ovarian GCTs induced by overactivation of TGFBR1 in granulosa cells express higher levels of Gli1 and Gli2 transcription factors and lower levels of Tgfbr3 compared with normal ovaries [16]. Interestingly, a similar transcript expression pattern of these genes was found in the ovaries of TGFBR1-CA G9Cre mice. These findings collectively suggest that ovarian tumor formation in TGFBR1-CA G9Cre mice may be attributable, at least partially, to overactivation of TGFBR1 in the granulosa cell compartment. Since the TGFBR1 CA is tagged with HA, we have tested the utility of a number of commercially available antibodies directed to HA in immunohistochemical/immunofluorescence applications, with the aim of identifying malignant granulosa cells harboring TGFBR1 overactivation. Although some antibodies performed reasonably in western blotting analysis, none of them generated reproducible and convincing results in immunohistochemical/immunofluorescence assays. Future work using Rosa26/ TGFBR1 CA /Gdf9-iCre mice may help elucidate the cellular origin of TGFBR1 overactivation and its contribution to ovarian GCT formation. It is worthwhile mentioning that a potential contribution of overactivation of TGFBR1 in the primordial oocyte or in both primordial oocytes and granulosa cells to GCT development could not be excluded. Further clarification of this question relays on the technical capability of specifically manipulating TGFB signaling in the oocyte of primordial follicles.
In summary, this study has created a mouse model of GCTs with defined disease onset that can be further exploited to study the role of TGFB signaling in ovarian tumor development. Of note, our studies were performed using mice and extrapolation of the findings to human GCTs needs further investigation. Identification of molecular markers for GCTs would benefit early diagnosis and treatment. To date, there are no reliable biomarkers for GCTs, although serum levels of estradiol, inhibin, and AMH have been extensively investigated and appear to correlate with GCTs in some cases [50]. Ongoing studies are to identify molecular signatures of GCTs in our model during tumor initiation and development.

Conclusions
Results of the current study reinforce the role of constitutively active TGFBR1 in promoting ovarian tumorigenesis in mice. The mouse model created in this study may be further exploited to define the cellular and molecular mechanisms of TGFB/activin downstream signaling in GCT development. Future studies are needed to test whether activation of TGFB/activin signaling contributes to the development of human GCTs.

Additional files
Additional file 1: Figure S1. Histological and morphological analysis of ovaries from immature TGFBR1-CA G9Cre and control mice. (A-H) Periodic acid Schiff's staining of ovarian samples from TGFBR1-CA G9Cre and control mice at PD7 and PD21. Panels (B, D, F, and H) represent higher magnification images for the corresponding panels (A, C, E, and G). Scale bar is representatively Fig. 7 Evidence of TGFBR1 activation in ovarian granulosa cells of TGFBR1-CA G9Cre mice. a-d Real-time PCR analysis of expression of TGFBR1 CA , Smad7, Inha, and Zp3 in ovaries from control and TGFBR1-CA G9Cre mice at PD3 and PD7. Real-time PCR was performed using ΔΔCT method. Data are mean ± s.e.m. n = 4-5. * P < 0.05, ** P < 0.01, and *** P < 0.001. Ns, not significant. e Western blotting analysis of TGFBR1 CA , phospho-SMAD2/3, and HSD3B using ovaries from 2-month-old control and TGFBR1-CA G9Cre mice. TGFBR1 CA was detected using an anti-HA antibody. n = 3. Each lane represents an independent sample. f-i X-gal staining using ovaries from Rosa26/Gdf9-iCre mice (f-h) and Rosa26 control mice (i). Panels g and h are higher magnification images of two different fields of panel (f). At least 3 independent samples per group were used. j Real-time PCR analysis of the expression of Gli1, Gli2, and Tgfbr3 using ovaries from 8-week-old control and TGFBR1-CA G9Cre mice. Data are mean ± s.e.m. n = 4-5. * P < 0.05 and *** P < 0.001. k-n RNAscope in situ hybridization analysis of Gli1 mRNA distribution using 8-week-old control (k) and TGFBR1-CA G9Cre ovaries (l). n = 4. Positive and negative controls using TGFBR1 CA flox/+ ovaries were shown in (m) and (n), respectively. Sections were counterstained with hematoxylin. Scale bar is representatively shown in (f) and equals 25 μm (g, h, and k-n) and 100 μm (f and i) shown in (A) and equals 25 μm (B, D, F, and H) and 100 μm (A, C, E, and G). (I) Reproductive tract of TGFBR1-CA G9Cre and control mice at PD21. OV, ovary; Ut, uterus. Scale bar = 10 mm (TIFF 3522 kb) Additional file 2: Figure S2. Western blotting analysis of phospho-SMAD2/ 3 and TGFBR1 CA using ovaries from 1-month-old control and TGFBR1-CA G9Cre mice. Note that a pronounced increase in phospho-SMAD2 but not phospho-SMAD3 was observed at this stage, suggesting that SMAD3 activation may have a later onset or is masked by high levels of phospho-SMAD3 within control ovaries at this stage. TGFBR1 CA was detected using an anti-HA antibody. ACTB was included as internal control. n = 3-4. Each lane represents an independent sample (TIFF 366 kb) Additional file 3: Figure S3.