Skip to main content
PLOS ONE logoLink to PLOS ONE
. 2022 Jul 8;17(7):e0268467. doi: 10.1371/journal.pone.0268467

Isolation, culture, and characterisation of bovine ovarian fetal fibroblasts and gonadal ridge epithelial-like cells and comparison to their adult counterparts

Menghe Liu 1, Katja Hummitzsch 1, Nicole A Bastian 1, Monica D Hartanti 1,2, Qianhui Wan 1, Helen F Irving-Rodgers 1,3, Richard A Anderson 4, Raymond J Rodgers 1,*
Editor: Rajakumar Anbazhagan5
PMCID: PMC9269465  PMID: 35802560

Abstract

During ovarian development, gonadal ridge epithelial-like (GREL) cells arise from the epithelial cells of the ventral surface of the mesonephros. They ultimately develop into follicular granulosa cells or into ovarian surface epithelial cells. Stromal fibroblasts arise from the mesonephros and penetrate the ovary. We developed methods for isolating and culturing fetal ovarian GREL cells and ovarian fibroblasts by expansion of colonies without passage. In culture, these two cell types were morphologically different. We examined the expression profile of 34 genes by qRT-PCR, of which 24 genes had previously been studied in whole fetal ovaries. Expression of nine of the 10 newly-examined genes in fetal ovaries correlated with gestational age (MUC1, PKP2, CCNE1 and CCNE2 negatively; STAR, COL4A1, GJA1, LAMB2 and HSD17B1 positively). Comparison between GREL cells and fetal fibroblasts revealed higher expression of KRT19, PKP2, OCLN, MUC1, ESR1 and LGR5 and lower expression of GJA1, FOXL2, NR2F2, FBN1, COL1A1, NR5A1, CCND2, CCNE1 and ALDH1A1. Expression of CCND2, CCNE1, CCNE2, ESR2 and TGFBR1 was higher in the fetal fibroblasts than in adult fibroblasts; FBN1 was lower. Expression of OCLN, MUC1, LAMB2, NR5A1, ESR1, ESR2, and TGFBR3 was lower in GREL cells than ovarian surface epithelial cells. Expression of KRT19, DSG2, PKP2, OCLN, MUC1, FBN1, COL1A1, COL3A1, STAR and TGFBR2 was higher and GJA1, CTNNB1, LAMB2, NR5A1, CYP11A1, HSD3B1, CYP19A1, HSD17B1, FOXL2, ESR1, ESR2, TGFBR3 and CCND2 was lower in GREL cells compared to granulosa cells. TGFβ1 altered the expression of COL1A1, COL3A1 and FBN1 in fetal fibroblasts and epidermal growth factor altered the expression of FBN1 and COL1A1. In summary, the two major somatic cell types of the developing ovary have distinct gene expression profiles. They, especially GREL cells, also differ from the cells they ultimately differentiate in to. The regulation of cell fate determination, particularly of the bi-potential GREL cells, remains to be elucidated.

Introduction

Understanding the development of the mammalian ovary is important for identifying the origins of different ovarian cell types and ovarian diseases. Despite much research in this area, the knowledge of ovarian development is still somewhat disjointed, particularly across species. It is commonly accepted that the ovary arises from the mesonephros [1, 2], but the origins of somatic ovarian cells is not clear. For instance, granulosa cells that are epithelial cells that surround each oocyte in follicles have been proposed to arise from the rete ovarii [35], the ovarian surface epithelium [6, 7] or a recently identified cell type, the gonadal ridge epithelial-like (GREL) cell [8]. Based upon morphological observations GREL cells are not only a precursor of granulosa cells but also of ovarian surface epithelial cells [8, 9]. Penetration of stromal fibroblasts from the mesonephros into the primordial ovary has consistently been observed in a number of species [6, 8, 10, 11].

In the ovary, granulosa cells contribute supportive functions to the oocyte via gap junctions [1214] and they communicate with the thecal cells [15, 16] that surround them. They play an important role in follicle selection and follicle development [17, 18]. Furthermore, enzymes which are essential to convert androgens into oestrogens [19, 20] are highly expressed in granulosa cells late in follicular development. On the other hand, the ovarian surface epithelium is a single layer of epithelial cells covering the surface of the ovary. It is involved in the recovery of the ovarian surface following the rupture of follicles at ovulation by cell proliferation. Lgr5-positive stem cells have been identified in the ovarian surface epithelium where they function as stem/progenitor cells for the recovery of the ovarian surface after ovulation [21, 22]. The differences between the two cell types which GREL cells differentiate into could be the result of regulation by oocyte-secreted factors such as growth differentiation factor 9 (GDF9) and bone morphogenetic protein 15 (BMP15) which play essential roles in regulating the proliferation and differentiation of granulosa cells [23]. The ovarian surface epithelium, on the other hand, is separated from oocytes by stroma and basal lamina and is likely not exposed to oocyte-secreted factors.

Stroma in the ovary provides structural support to this organ, but it also plays active roles. During the development of the ovary, stroma from the mesonephros penetrates into the primordial fetal ovary, and branches as it does so. This branching thus partitioning the GREL cells and oogonia thereby creates the ovigerous cords. The stroma further partitions the ovigerous cords into smaller groups of cells which become primordial follicles containing an oocyte, derived from oogonia, and granulosa cells, derived from GREL cells. The stroma penetrates to just below the surface of the ovary and spreads laterally and the GREL cells on the surface differentiate into ovarian surface epithelium [8, 24]. The stroma continues to play significant roles in adult ovaries. In vitro culture of follicles demonstrates that the strength and rigidity of an artificial matrix can influence the growth of follicles [25], which suggests that in vivo the strength and rigidity of stroma surrounding follicles can influence follicle development. The fibroblast is the major cell type in stroma, and hormones and growth factors can regulate the proliferation and extracellular matrix (ECM) production of fibroblasts [26].

To conduct the first examination of GREL cells and to compare them to fetal fibroblasts, we examined the expression of a range of genes related to epithelial cells, stromal cells, steroidogenesis, signalling receptors, cell cycle and germ/stem cells. We also examined the expression of these genes in the adult cells that the GREL cells and fetal fibroblast differentiate in to. We also examined potential regulators of expression of stromal-related genes in fetal fibroblast.

Materials and methods

Collection of bovine fetal and adult ovaries

Fetuses from pregnant Bos taurus cows and ovaries from non-pregnant Bos taurus cows were collected at local abattoirs (Thomas Foods International, Murray Bridge, SA, Australia and Strath Meats, Strathalbyn, SA, Australia) and transported on ice in Hank’s balanced-salt solution (HBSS) with calcium and magnesium (HBSS+/+) to the laboratory. The gestational ages of the fetuses were estimated from the crown-rump length (CRL) [27]. Fetuses with a CRL < 8 cm underwent a sex determination as previously reported [8]. Samples of whole fetal ovaries throughout gestation (n = 27) were from our previous study [28] and collected from The Midfield Group, Warrnambool, VIC, Australia. The bovine ovaries used in this study were scavenged from animals that were being processed for human consumption and were not owned by the authors or their institutions. As such the University of Adelaide’s Animal Ethics Committee only requires notification of this.

Isolation of fetal ovarian GREL cells and fibroblasts

Ovaries were dissected from the fetus and rinsed once with 70% ethanol and twice with HBSS+/+. Then the fetal ovaries were fully minced and digested in 10 ml HBSS+/+ containing 1 mg/ml collagenase type I (GIBCO, Carlsbad, CA, USA) and 50 μg/ml DNase I (Sigma-Aldrich, New South Wales, Australia) at 37°C and shaking at 100 rpm. Minced ovaries from fetuses (CRL ≤ 20 cm) were digested for 30 min and ovaries from fetuses (CRL > 20 cm) were digested for 45 min. The cell suspensions were then passed through cell strainers and centrifuged at 1200 rpm for 5 min and the supernatants were discarded. Subsequently, the cell pellets were washed with 5 ml HBSS+/+ and digested with 1 ml red blood cell lysis solution (Sigma-Aldrich, New South Wales, Australia; Cat # R7757) by gently pipetting up and down for 45 sec. Thereafter, cell pellets were resuspended in 1 ml DMEM/F12 medium (Sigma-Aldrich, New South Wales, Australia; Cat # D8900) containing 5% fetal calf serum (Sigma-Aldrich, New South Wales, Australia), 100 IU/ml penicillin, 0.01% streptomycin sulfate (GIBCO, Carlsbad, CA, USA; Cat # 15140122), 0.1% fungizone (GIBCO, Carlsbad, CA, USA. Cat # 15290–018) and 3 ng/ml human epidermal growth factor (EGF; Sigma-Aldrich, New South Wales, Australia; Cat # E9644). Total cell number was determined using trypan blue method and a Neubauer haemocytometer. Then 250 cells were seeded per wellin 24-well plates which were pre-coated with collagen type I (Sigma-Aldrich, New South Wales, Australia) at a density of 10 μg/cm2. The cells were then incubated at 5% CO2 and 38.5°C. The culture medium was changed every two days and the growth of purified GREL or fibroblast clusters was checked every day by light microscopy (Olympus, Tokyo, Japan). Clusters with characteristic morphology were observed from day 2 onwards (S1 Fig). The majority of wells contained cultures of mixed cell types with fibroblast clusters and only one or two large purified GREL cell clusters in each 24-well plate. The largest purified and healthy clusters of GREL cells or fetal fibroblasts were selected on day 3 and day 4 (S1 Fig). In case the clusters were surrounded by some cells of another cell type (S1 Fig), these cells were discarded by scraping with a sterile blunt-ended glass Pasteur pipette and subsequent culture medium change. When the cell clusters were confluent (S2 and S3 Figs), GREL cells and fibroblasts were collected. To collect pure GREL cell or fibroblast cultures for subsequent RNA extraction, the clusters were scraped with a blunt-ended Pasteur pipette and the culture media with the cells transferred into sterile Eppendorf tubes. After centrifugation at 1500 rpm for 5 min, the cell pellets were stored at -80°C until required.

Isolation of adult ovarian fibroblasts

After collection, any connective tissue surrounding the adult ovaries was discarded and the ovaries washed with 70% ethanol and HBSS+/+. Individual ovaries were cut into quarters and tissue areas devoid of visible follicles was cut to a depth of about 0.5 mm from the surface and minced with sterile razor blades and digested with HBSS+/+ containing 1 mg/ml collagenase type I and 50 μg/ml DNase I at 37°C shaking at 100 rpm for 1 to 2 h. Subsequently, cell suspensions were passed through 70 μm cell strainers (Sartorius, Dandenong South, VIC, Australia) and centrifuged at 1500 rpm for 5 min. The cell pellet was resuspended in 5 ml 0.05% trypsin/EDTA (GIBCO, Carlsbad, CA, USA/Life Technologies Australia Pty Ltd, Tullamarine, Vic, Australia) in HBSS without calcium or magnesium (HBSS-/-) for digestion with pipetting, and the digestion was stopped with DMEM/F12 containing 5% fetal calf serum, 100 IU/ml penicillin, 0.01% streptomycin sulfate, 0.1% fungizone and 3 ng/ml EGF after 5 min. Total cell number was determined as described above for fetal cells and the adult fibroblasts seeded into 24-well plates and cultured in the same medium and way as the fetal ovarian cells. Purified clusters were collected and cell pellets were stored at -80°C for subsequent use.

Collection of adult ovarian surface epithelial cells

Individual adult ovaries were placed into a 60 mm culture dish and the connective tissues around the ovary was removed with forceps and scissors. Subsequently, ovaries were moved to another dish with 12 ml DMEM/F12 containing 5% fetal calf serum, 100 IU/ml penicillin, 0.01% streptomycin sulfate and 0.1% fungizone and the surface of each ovary was scraped with the blunt side of a scalpel blade. Ovarian surface epithelium from five to seven ovaries were pooled and transferred into a 15 ml centrifuge tube and centrifuged at 1500 rpm for 5 min. Supernatant was discarded and cell pellets were stored at -80°C until required. In total 5 pools were collected.

Collection of adult granulosa cells

Granulosa cells were isolated from the largest healthy antral follicle (9–12 mm) from each adult ovarian pair. The diameters of the follicles were measured with callipers. The follicles were cut out with scissors and transferred to a dish containing 1 ml HBSS-/-. The follicles were cut open with scissors and granulosa cells scraped from the follicle wall with a blunt-ended Pasteur pipette. The cell suspension from each individual follicle was transferred into a 15 ml tube containing HBSS-/- and centrifuged at 1500 rpm for 5 min. Then the cell pellets were resuspended in 1 ml HBSS-/- and the ratio of alive/dead cells determined using the trypan blue method and a Neubauer haemocytometer. Cell suspensions containing more than 70% live cells were centrifuged at 1500 rpm for 5 min, and cell pellets stored at -80°C for subsequent use.

Treatment of fetal fibroblasts

Fetal fibroblasts were from a previous study and had been cultured and treated as described previously [26]. Note that these cells were not cultured with EGF during their isolation culture. Briefly, bovine fetal fibroblasts were seeded in 24-well plates at a density of 30,000 cells/well in DMEM/F12 containing 5% fetal calf serum, 1% penicillin and streptomycin sulphate and 0.1% fungizone. When cells were 60–70% confluent, cells were treated with different agents for 18 h. After 18 h, cells were harvested and stored at -80°C. Expression of NR2F2, FBN1, COL1A1 and COL3A1 was subsequently measured. Cells that were treated with 24 different agents as listed in S1 Table were from weeks 13, 14, 17, 19 and 33 of gestation (n = 5) and those treated with TGFβ were from gestational weeks 19–26 (n = 6).

RNA isolation, cDNA synthesis and quantitative real-time PCR

RNA from granulosa cells, ovarian surface epithelial cells and adult fibroblasts was extracted with TRIzolTM (Thermo Fisher Scientific, Waltham, MA, USA) according to manufacturer’s instructions. RNA from GREL cells and fetal fibroblasts was extracted with RNAqueousTM- Micro Total RNA Isolation kit. Ten μg of RNA of each sample was treated with DNase I (Promega Australia, Alexandria, NSW, Australia/Life Technologies Australia Pty Ltd, Tullamarine, Vic, Australia) for 20 min at 37°C. Two hundred ng of DNAse-treated RNA was used for cDNA synthesis as previously described [8, 29]. Based on available RNA sequences in NCBI, PCR primers were designed to span introns using Primer3 plus and Net primer software (PREMIER Biosoft Palo Alto, CA, USA) and primers are listed in S2 Table. Primer combinations were tested as previously described [29]. No addition of reverse transcriptase was used to confirm no amplification of any potential genomic DNA contamination. The amplification conditions were 95°C for 15 sec, then 60°C for 60 sec for 40 cycles using Rotor Gene 6000 cycler (Q series, Qiagen GmbH, Hilden, Germany). Ct values were determined using Rotor Gene 6000 software at a threshold of 0.05 normalised fluorescent unit. Gene expression values were determined using 2-ΔCt method and ribosomal protein L32 (RPL32) and peptidylprolyl isomerase A (PPIA) were used as housekeeping genes. The fold change of gene expression in fetal fibroblasts treated with different agents were presented as 2-ΔΔCt data with the same housekeeping genes in 2-ΔCt method. We analysed the expression levels of 34 genes relating to cell junctions, extracellular matrix, enzymes and transcription factors in order to characterise the cultured GREL cells and fetal fibroblasts.

Statistical analysis

All statistical analyses were carried out using Microsoft Office Excel 365 (Microsoft Redmond, WA, USA) and IBM SPSS Statistics for Windows, version 25 (IBM Corp., Armonk, N.Y., USA). Scatter plots and columns were made using GraphPad Prism version 8.00 (GraphPad Software Inc., La Jolla, CA, USA). Pearson correlation test was used to analyse the correlation between gene expression levels and gestational age. The differences between GREL cells and fetal fibroblasts, the differences between fetal fibroblasts and adult fibroblasts, and the differences between fetal fibroblasts cultured with or without EGF were analysed with Student’s t-test. To determine significant differences among adult fibroblasts, granulosa cells and ovarian surface epithelial cells, and differences between GREL cells with adult ovarian surface epithelial cells or granulosa cells, one-way ANOVA and Tukey’s post-hoc tests were conducted. To compare the difference between the control and the treatments in fetal fibroblasts, one-way ANOVA with Dunnett’s post-hoc tests were conducted.

Results

Morphology of GREL cells and fetal fibroblasts in culture

Early in culture (day 2 and 3), GREL cells in clusters appeared to be of irregularly polygonal shape, with clear visible nuclei, and in some areas tightly packed (S1A and S1C Fig), as previously observed [8]. In the following days, cells appeared to be more tightly connected (S1E and S1F Fig) and remained so until the cells were ready for harvesting (Fig 1A and S2 Fig). In contrast, fetal fibroblasts showed the typical spindle shaped appearance as also seen in adult fibroblasts (Fig 1B, S1B and S1D Fig, S3 Fig). The purified GREL and fetal fibroblast clusters demonstrated that the current method of isolating GREL and fetal fibroblast was feasible. In order to explore more characters of GREL cells and fetal fibroblasts and to make comparisons with their adult counterparts, we conducted analysis of expression of genes including those involved in cellular junctions, extracellular matrix, hormone production, transcription factors, oestrogen and TGFβ signalling, cell cycle and stem and germ cells.

Fig 1. Morphology of fetal somatic cells.

Fig 1

Representative micrograph image of (A) GREL cells and (B) fetal fibroblasts. The cell types were isolated from fetal ovaries with gestational ages of 56 days and 73 days, respectively. Bars: 100 μm.

Presentation of gene expression data

In Figs 2 to 7 we present expression data of the 34 genes examined and then group them based upon the functionally of genes for each Figure. For each gene the expression across gestation is presented, and for 24 of the 34 genes this had been published previously and is acknowledged and reprinted in the current Figs for ease of comparison. The data from the other 10 genes are new data presented in this manuscript. We then present the expression in fetal fibroblasts and GREL cells and the statistical comparison between them. We then present the expression in adult fibroblasts, granulosa cells and ovarian surface epithelial cells and the statistical comparison between them. The data for each gene from gestation and from isolated cells in Figs 2 to 7 is plotted on the same scale for direct comparisons. The mean levels of expression in isolated cells and statistical comparisons of fetal and adult cells in presented in Table 1. The individual data of expression in fetal fibroblasts and GREL cells across gestation is presented in a supplementary figure. The data with correlations between the 10 new genes across gestation and age of gestation is presented in a supplementary table and cannot be combined with previously published data on the other 24 genes as the RNA used for the 10 new genes is from a different collection of fetal ovaries.

Fig 2. mRNA expression levels of genes associated with cellular junctions.

Fig 2

Scatter plots in each subfigure show gene expression levels in fetal ovaries at different gestational ages (n = 27), and scatter plots for KRT19 (A), DSG2 (B) and CTNNB1 (G) are from our previous publications [30, 31]. Gene expression levels in cultured GREL cells (n = 7; gestational ages are 51, 56, 73, 110, 110, 127 and 177 days) and fetal fibroblasts (n = 6; gestational ages are 51, 73, 101, 110, 110 and 177 days) are presented as mean ± SEM (normalised to PPIA and RPL32) in the middle graph of each subfigure. The differences between GREL cells and fetal fibroblasts were analysed with Student’s t-test. Data of gene expression in adult fibroblasts (n = 4), granulosa cells (GC; n = 5) and ovarian surface epithelium (OSE; n = 5) are presented as mean ± SEM (normalised to PPIA and RPL32) in the graph on the right side of each subfigure. Significant differences among the three adult cell types were determined by one-way ANOVA with Tukey’s post-hoc test. * P < 0.05, ** P < 0.01, *** P < 0.001.

Fig 7. mRNA expression levels of stem and germ cell genes.

Fig 7

Scatter plots in each subfigure show gene expression levels in fetal ovaries at different gestational ages (n = 27) and are from our previous publications [30, 31]. Gene expression levels in cultured GREL cells (n = 7; gestational ages are 51, 56, 73, 110, 110, 127 and 177 days) and fetal fibroblasts (n = 6; gestational ages are 51, 73, 101, 110, 110 and 177 days) are presented as mean ± SEM (normalised to PPIA and RPL32) in the middle graph of each subfigure. The differences between GREL cells and fetal fibroblasts were analysed with Student’s t-tests. Data of gene expression in adult fibroblasts (n = 4), granulosa cells (GC; n = 5) and ovarian surface epithelium (OSE; n = 5) are presented as mean ± SEM (normalised to PPIA and RPL32) in the graph on the right side of each subfigure. Significant differences among the three adult cell types were determined by one-way ANOVA with Tukey’s post-hoc test. * P < 0.05, ** P < 0.01, *** P < 0.001.

Table 1. Relative gene expression levels (2-ΔCt) are presented as mean ± SEM for each cell type.

Comparisons are between fetal fibroblasts (n = 6) with adult fibroblasts (n = 4), and then fetal GREL cells (n = 7) with adult ovarian surface epithelial cells (n = 5) or granulosa cells (n = 5).

Gene Fetal fibroblasts Adult fibroblasts Fetal GREL cells Adult ovarian surface epithelial cells Adult granulosa cells
KRT19 0.542±0.203 0.286±0.120 5.044±1.313 3.090±0.515 0.003±0.001**
DSG2 0.059±0.024 0.027±0.011 0.132±0.048 0.041±0.005 0.000±0.000*
PKP2 0.022±0.005 0.01±0.002 0.080±0.024 0.106±0.013 0.012±0.003***
OCLN 0.001±0.000 0.006±0.005 0.017±0.005 0.272±0.057*** 0.003±0.000**
MUC1 0.004±0.001 0.005±0.003 0.009±0.002 0.180±0.027*** 0.002±0.001***
GJA1 2.640±0.332 3.732±0.798 1.461±0.319 1.507±0.287 17.049±5.306**
CTNNB1 0.281±0.047 0.177±0.023 0.290±0.064 0.504±0.071 0.922±0.302*
LAMB2 0.193±0.026 0.130±0.037 0.165±0.037 0.856±0.191** 0.650±0.107*
FBN1 0.913±0.103 1.425±0.161* 0.374±0.104 0.211±0.029 0.012±0.003**
COL1A1 16.111±2.289 14.009±2.382 4.653±1.984 3.907±0.640 0.043±0.021***
COL3A1 6.116±1.553 4.484±0.453 3.613±1.339 2.186±0.381 0.009±0.003*
COL4A1 1.797±0.845 0.329±0.081 0.752±0.418 0.301±0.129 0.441±0.127
NR5A1 0.061±0.009 0.028±0.006* 0.016±0.007 0.120±0.021** 0.416±0.105***
STAR 0.038±0.011 0.007±0.002 0.018±0.005 0.009±0.002 0.002±0.001***
CYP11A1 0.005±0.002 0.039±0.017 0.001±0.001 0.127±0.077*** 3.758±1.128***
HSD3B1 0.000±0.000 0.000±0.000 0.000±0.000 0.005±0.002*** 0.316±0.116***
CYP19A1 0.001±0.000 0.000±0.000 0.000±0.000 0.056±0.051* 10.307±2.106***
HSD17B1 0.000±0.000 0.000±0.000 0.000±0.000 0.003±0.001*** 0.328±0.157***
INHBA 0.002±0.000 0.001±0.000 0.001±0.000 0.008±0.005 0.005±0.002
NR2F2 0.559±0.128 0.455±0.096 0.135±0.037 0.172±0.024 0.050±0.016
FOXL2 0.051±0.008 0.024±0.006* 0.012±0.006 0.028±0.006 1.015±0.280***
ESR1 0.001±0.001 0.002±0.001 0.014±0.007 0.570±0.094** 0.116±0.010*
ESR2 0.007±0.001 0.003±0.000** 0.009±0.006 0.022±0.003* 0.144±0.039***
TGFBR1 0.148±0.033 0.06±0.008* 0.153±0.036 0.136±0.013 0.088±0.015
TGFBR2 0.105±0.026 0.05±0.011 0.129±0.037 0.533±0.076* 0.007±0.006***
TGFBR3 0.195±0.058 0.065±0.038 0.125±0.027 0.483±0.096** 0.325±0.083*
CCND2 0.096±0.011 0.018±0.003** 0.038±0.011 0.111±0.032 0.480±0.156**
CCNE1 0.038±0.005 0.009±0.002** 0.018±0.004 0.025±0.011 0.024±0.006
CCNE2 0.039±0.005 0.005±0.001** 0.024±0.005 0.022±0.007 0.060±0.0170
ALDH1A1 0.095±0.026 0.102±0.030 0.019±0.003 0.038±0.011 0.008±0.002
LGR5 0.000±0.000 0.000±0.000 0.001±0.000 0.045±0.007* 0.000±0.000
OCT4 0.012±0.001 0.007±0.001* 0.011±0.003 0.039±0.006** 0.014±0.003
DAZL 0.003±0.001 0.001±0.000 0.001±0.000 0.003±0.001 0.001±0.000
VASA 0.001±0.000 0.000±0.000* 0.001±0.000 0.001±0.001 0.000±0.000

Student’s t-tests were used to analyse the differences between fetal fibroblasts and adult fibroblasts. One-way ANOVA with Tukey’s post-hoc test was used to analyse the differences between GREL cells with adult ovarian surface epithelial cells or granulosa cells

* P < 0.05

** P < 0.01

*** P < 0.001 report significant differences between the corresponding fetal cells with the adult cells.

Cellular junctions

In the fetal ovary, KRT19 (Fig 2A), OCLN (Fig 2D and S3 Table), CTNNB1 and (Fig 2G) were expressed constantly throughout the gestation. Expression of DSG2 (Fig 2B), PKP2 (Fig 2C and S3 Table) and MUC1 (Fig 2E and S3 Table) were negatively correlated with gestational age, however, expression of GJA1 (Fig 2F and S3 Table) was positively correlated with gestational age. In GREL cells and fetal fibroblasts, expression of KRT19 (Fig 2A), PKP2 (Fig 2C), OCLN (Fig 2D) and MUC1 (Fig 2E) was significantly higher in GREL cells than fetal fibroblasts. DSG2 (Fig 2B) did not exhibit significant differences between GREL cells and fetal fibroblasts, even though most GREL cells exhibited higher DSG2 expression (S4A2 Fig) than fetal fibroblasts from the same fetus. GJA1 (Fig 2F) had significantly higher expression in fetal fibroblasts. Expression of CTNNB1 (Fig 2G) did not exhibit significant differences (P > 0.05) between GREL cells and fetal fibroblasts. In adult cells, expression of KRT19 (Fig 2A), PKP2 (Fig 2C), OCLN (Fig 2D) and MUC1 (Fig 2E) was greater in ovarian surface epithelial cells, whereas GJA1 (Fig 2F) was higher in granulosa cells. DSG2 (Fig 2B) was more highly expressed in ovarian surface epithelial cells and adult fibroblasts than in granulosa cells. Additionally expression of CTNNB1 (Fig 2G) was not significantly different between the three cell types.

Extracellular matrix

In the fetal ovary, expression of LAMB2 (Fig 3B and S3 Table), FBN1 (Fig 3C), COL1A1 (Fig 3D), COL3A1 (Fig 3E) and COL4A1 (Fig 3F and S3 Table) were positively correlated with the gestational age. FBN1 (Fig 3C) and COL1A1 (Fig 3D) levels were significantly higher in fetal fibroblasts than in GREL cells. LAMB2 (Fig 3B), COL3A1 (Fig 3E) and COL4A1 (Fig 3F) expression did not exhibit significant differences between GREL cells and fetal fibroblasts. In adult cells, adult fibroblasts had significant higher expression of FBN1 (Fig 3C) and COL1A1 (Fig 3D) and significantly lower expression of LAMB2 (Fig 3B) than ovarian surface epithelial cells and granulosa cells. Ovarian epithelial cells and adult fibroblast had higher expression of COL3A1 (Fig 3E) than granulosa cells. COL4A1 (Fig 3F) did not have significant differences among the three cell types.

Fig 3. mRNA expression levels of extracellular matrix genes.

Fig 3

Scatter plots in each subfigure show gene expression levels in fetal ovaries at different gestational ages (n = 27), and scatter plots for FBN1 (B), COL1A1 (C) and COL3A1 (D) are from our previous publications [30, 31]. Gene expression levels in cultured GREL cells (n = 7; gestational ages are 51, 56, 73, 110, 110, 127 and 177 days) and fetal fibroblasts (n = 6; gestational ages are 51, 73, 101, 110, 110 and 177 days) are presented as mean ± SEM (normalised to PPIA and RPL32) in the middle graph of each subfigure. The differences between GREL cells and fetal fibroblasts were analysed with Student’s t-tests. Data of gene expression in adult fibroblasts (n = 4), granulosa cells (GC; n = 5) and ovarian surface epithelium (OSE; n = 5) are presented as mean ± SEM (normalised to PPIA and RPL32) in the graph on the right side of each subfigure. Significant differences among the three adult cell types were determined by one-way ANOVA with Tukey’s post-hoc test. * P < 0.05, ** P < 0.01, *** P < 0.001.

Genes associated with hormone production

In the fetal ovary, expression levels of NR5A1 (Fig 4A), CYP11A1 (Fig 4C) and HSD3B1 (Fig 3D) were not significantly correlated with gestational age. STAR (Fig 4B and S3 Table), HSD17B1 (Fig 4F and S3 Table) and INHBA (Fig 4G) were positively correlated with gestational age but CYP19A1 (Fig 4E) was negatively correlated with gestational age. Fetal fibroblasts had higher expression of NR5A1 (Fig 4A) than GREL cells. Likewise, CYP11A1 (Fig 4C) and INHBA (Fig 4G) had also a significantly higher expression in fetal fibroblasts, even though it was very low. STAR did not show significance differences between GREL cells and fetal fibroblasts (Fig 4B), even if most GREL cells had higher STAR expression than fetal fibroblasts of the same age (S4C2 Fig). Expression of HSD3B1 (Fig 4D), CYP19A1 (Fig 4E) and HSD17B1 (Fig 4F) exhibited no significant differences between GREL cells and fetal fibroblasts, and they were expressed at very low levels. In adult cells, CYP11A1 (Fig 4C), HSD3B1 (Fig 4D) and CYP19A1 (Fig 4E) were more highly expressed in granulosa cells. NR5A1 (Fig 3A) was more highly expressed in both granulosa cells and ovarian surface epithelium than adult fibroblasts. Ovarian surface epithelial cells and adult fibroblast had higher expression of STAR (Fig 4B) than granulosa cells. HSD17B1 (Fig 4F) and INHBA (Fig 4G) did not exhibit significant differences among the three cell types.

Fig 4. mRNA expression levels of genes involved in hormone production.

Fig 4

Scatter plots in each subfigure show gene expression levels in fetal ovaries at different gestational ages (n = 27), and scatter plots for NR5A1 (A), CYP11A1(C), HSD3B1 (D), CYP19A1 (E) and INHBA (G) are from our previous publications [30, 31]. Gene expression levels in cultured GREL cells (n = 7; gestational ages are 51, 56, 73, 110, 110, 127 and 177 days) and fetal fibroblasts (n = 6; gestational ages are 51, 73, 101, 110, 110 and 177 days) are presented as mean ± SEM (normalised to PPIA and RPL32) in the middle graph of each subfigure. The differences between GREL cells and fetal fibroblasts were analysed with Student’s t-tests. Data of gene expression in adult fibroblasts (n = 4), granulosa cells (GC; n = 5) and ovarian surface epithelium (OSE; n = 5) are presented as mean ± SEM (normalised to PPIA and RPL32) in the graph on the right side of each subfigure. Significant differences among the three adult cell types were determined by one-way ANOVA with Tukey’s post-hoc test. * P < 0.05, *** P < 0.001.

Transcription factors, oestrogen receptors and TGFβ signalling pathways

In the fetal ovary, NR2F2 (Fig 5A), ESR1 (Fig 5C), TGFBR2 (Fig 5F) and TGFBR3 (Fig 5G) were positively correlated with gestational age but ESR2 (Fig 5D) was negatively correlated with gestational age. TGFBR1 (Fig 5E) and FOXL2 (Fig 5B) were expressed constantly during the gestation. In GREL cells and fetal fibroblasts, ESR1 (Fig 5C) was expressed significantly higher in GREL cells. ESR2 (Fig 5D) was hardly detectable in GREL cells or fetal fibroblasts. FOXL2 (Fig 5B) and NR2F2 (Fig 5A) was expressed significantly higher in fetal fibroblasts. GREL cells and fetal fibroblasts did not exhibit significant differences in expression of TGFBR1 (Fig 5E), TGFBR2 (Fig 5F) and TGFBR3 (Fig 5G). In adult cells, ESR1 (Fig 5C), TGFBR1 (Fig 5E) and TGFBR2 (Fig 5F) were more highly expressed in ovarian surface epithelial cells and ESR2 (Fig 5D) and FOXL2 (Fig 5B) were more highly expressed in granulosa cells. TGFBR3 (Fig 5G) was expressed more in both surface epithelial cells and granulosa cells than in adult fibroblasts. Adult fibroblasts had significant higher expression of NR2F2 (Fig 5A) than ovarian surface epithelial cells and granulosa cells.

Fig 5. mRNA expression levels of transcription factors, oestrogen receptors and TGFβ signalling pathways.

Fig 5

Scatter plots left in each subfigure of gene expression levels in fetal ovaries at different gestational ages (n = 27) are from our previous publications [30, 31]. Gene expression levels in cultured GREL cells (n = 7; gestational ages are 51, 56, 73, 110, 110, 127 and 177 days) and fetal fibroblasts (n = 6; gestational ages are 51, 73, 101, 110, 110 and 177 days) are presented as mean ± SEM (normalised to PPIA and RPL32) in the middle graph of each subfigure. The differences between GREL cells and fetal fibroblasts were analysed with Student’s t-tests. Data of gene expression in adult fibroblasts (n = 4), granulosa cells (GC; n = 5) and ovarian surface epithelium (OSE; n = 5) are presented as mean ± SEM (normalised to PPIA and RPL32) in the graph on the right side of each subfigure. Significant differences among the three adult cell types were determined by one-way ANOVA with Tukey’s post-hoc test. * P < 0.05, ** P < 0.01, *** P < 0.001.

Cell cycle genes

In the fetal ovary, expression of CCND2 (Fig 6A), CCNE1 (Fig 6B and S3 Table) and CCNE2 (Fig 6C and S3 Table) was negatively correlated with gestational age. CCNE1 and CCNE2 levels reached to a peak around 100 days of gestation before declining. In GREL cells and fetal fibroblasts, expression of CCND2 (Fig 6A) and CCNE1 (Fig 6B) was significantly higher in fetal fibroblasts, but CCNE2 (Fig 6C) did not exhibit significant differences between the two cell types. In adult cells, expression of CCND2 (Fig 6A) was higher in granulosa cells. CCNE2 (Fig 6C) was more highly expressed in both granulosa cells and ovarian surface epithelium than in adult fibroblasts, however, there was no statistically significant difference between adult ovarian surface epithelial cells and adult fibroblasts. Expression of CCNE1 (Fig 6B) was not significantly difference among the three adult cell types.

Fig 6. mRNA expression levels of cell cycle genes.

Fig 6

Scatter plots in each subfigure show gene expression levels in fetal ovaries at different gestational ages (n = 27), and the scatter plot for CCND2 (A) is from our previous publication [31]. Gene expression levels in cultured GREL cells (n = 7; gestational ages are 51, 56, 73, 110, 110, 127 and 177 days) and fetal fibroblasts (n = 6; gestational ages are 51, 73, 101, 110, 110 and 177 days) are presented as mean ± SEM (normalised to PPIA and RPL32) in the middle graph of each subfigure. The differences between GREL cells and fetal fibroblasts were analysed with Student’s t-tests. Data of gene expression in adult fibroblasts (n = 4), granulosa cells (GC; n = 5) and ovarian surface epithelium (OSE; n = 5) are presented as mean ± SEM (normalised to PPIA and RPL32) in the graph on the right side of each subfigure. Significant differences among the three adult cell types were determined by one-way ANOVA with Tukey’s post-hoc test. * P < 0.05, ** P < 0.01, *** P < 0.001.

Stem or germ cell genes

In the fetal ovary, expression of ALDH1A1 (Fig 7A) was positively correlated with gestational age but OCT4 (Fig 7C) was negatively correlated with gestational age. LGR5 (Fig 7B), DAZL (Fig 7D) and VASA (Fig 7E) were constantly expressed throughout the gestation. In GREL cells and fetal fibroblasts, expression of ALDH1A1 (Fig 7A) was significantly higher in fetal fibroblasts but LGR5 (Fig 7B) was significantly higher in GREL cells. OCT4 (Fig 7C), DAZL (Fig 7D) and VASA (Fig 7E) were hardly detected in GREL cells and fetal fibroblasts as expected. In adult cells, LGR5 (Fig 7B) and OCT4 (Fig 7C) had higher expression in the ovarian surface epithelium. Similarly, ALDH1A1 (Fig 7A) had higher expression in both ovarian surface epithelium and fibroblasts than granulosa cells. DAZL (Fig 7D) and VASA (Fig 7E) were hardly detectable in the three cell types.

Fetal versus adult cells

CCND2, CCNE1, CCNE2, ESR2 and TGFBR1 were more highly expressed in the fetal fibroblasts than in adult fibroblasts (Table 1) whereas FBN1 was lower in the fetal fibroblasts. Expression of OCLN, MUC1, LAMB2, NR5A1, ESR1, ESR2 and TGFBR3 was lower in GREL cells than ovarian surface epithelial cells (Table 1). Expression of NR5A1, CYP11A1, HSD3B1, CYP19A1 and HSD17B1 was statistically lower in GREL cells than ovarian surface epithelial cells but all were much lower than granulosa cells (Table 1). Expression of KRT19, DSG2, PKP2, OCLN, MUC1, FBN1, COL1A1, COL3A1, STAR and TGFBR2 was higher in GREL cells compared to granulosa cells whereas expression of GJA1, CTNNB1, LAMB2, NR5A1, CYP11A1, HSD3B1, CYP19A1, HSD17B1, FOXL2, ESR1, ESR2, TGFBR3 and CCND2 was lower (Table 1).

Regulation of genes in fetal fibroblasts

Fetal fibroblasts previously isolated and cultured without EGF [26] and treated with 24 of hormones, stimulators and growth factors were analysed for the expression of NR2F2, FBN1, COL1A1 and COL3A1 which are commonly considered as genes expressed by ovarian fibroblasts. Results showed that many of the treatments had no significant effect on the expression of those genes (S5 Fig). However, results of treatment of doses of TGFβ1 (5 or 20 ng/ml) on fetal fibroblasts (n = 6 from the second trimester of gestation) showed significantly increased FBN1 expression in fetal fibroblast and significantly decreased the expression of COL1A1 and COL3A1 (Fig 8).

Fig 8. Normalised gene expression genes in fetal fibroblasts.

Fig 8

Fibroblasts from bovine fetal ovaries (19–26 weeks, n = 6) were treated with 5 or 20 ng/ml TGFβ1 for 18 h. Fold changes of gene expression to the control groups are presented as mean ± SEM (normalised to PPIA and RPL32). Significant differences between groups were determined by one-way ANOVA with Tukey’s post-hoc test. ** P < 0.01, *** P < 0.001.

We compared the fetal fibroblast cultured in EGF-containing medium (n = 6 from 51, 56, 73, 110, 127 and 177 days of gestation) isolated in the current experiments with fetal fibroblasts cultured in EGF-free medium (n = 11 from the first, second and third trimester of gestation) from a previous study [26]. The result showed that EGF significantly increased FBN1 expression, but decreased the gene expression of COL1A1 (S6 Fig).

Discussion

In this study we developed methods for isolating and culturing fetal ovarian GREL cells and fetal ovarian fibroblasts by removal of contaminating cells during culture and expansion from single colonies without passage. We had to adopt this approach because we had not identified a specific antigen on the surface of GREL cells to allow us to utilise more direct cell isolation methods such as the use of magnetic beads or fluorescence-activated cell sorting. We investigated the expression of a variety of genes to identify properties of GREL cells and fetal fibroblasts and examined differences with their adult counterparts. It should be acknowledged that culturing of GREL cells could have contributed to differences observed with the granulosa cells or ovarian surface epithelial cells which were not cultured. Furthermore, we additionally carried out a series of treatments to identify possible regulators of extracellular matrix in fetal fibroblasts. This is the first study of GREL cells since their identification.

Of the 10 newly-examined genes and the 24 previously-examined genes in whole fetal ovaries across gestation [30, 31] that we examined here in in-vitro cultured GREL cells and fetal fibroblasts, many were expressed early in gestation (MUC1, PKP2, CCNE1, CCNE2, CCND2, CYP19A1, DSG2, ESR2 and OCT4) and were negatively correlated with gestational age. Others were expressed late in gestation (ALDH1A1, COL1A1, COL3A1, ESR1, FBN1, INHBA, NR2F2, TGFBR2, TGFBR3, STAR, COL4A1, GJA1, LAMB2 and HSD17B1) and were positively correlated with gestational age. The expression of NR2F2, FBN1, COL1A1 and COL3A1 that were all positively correlated with the gestational age, could be as a result of the increasing volume of stroma in the fetal ovary [32]. Many of the 10 newly-examined genes were either negatively or positively correlated with each other. These patterns of gene expression and correlations between genes have been observed previously [2831] and reinforce the major developmental processes that the ovary undergoes during its development [8].

Many earlier studies of fetal ovaries focused on germ cells and often refer to the other cells as somatic cells which to some degree creates the impression that they are one cell type. However, our earlier studies identified that stromal fibroblasts and GREL cells have different cellular origins and exist in the ovary in different cellular compartments separated by a basal lamina [8]. In the current study the morphology of the cells, when cultured, differed substantially. The fibroblasts had a classical elongated shape and the GREL cells were of irregular polygonal shape and were tightly packed together in colonies. Their gene expression profiles were also different. Fifteen out of the 34 genes examined were differentially expressed with GREL cells having higher expression of KRT19, PKP2, OCLN, MUC1, ESR1 and LGR5 and lower expression of GJA1, FOXL2, NR2F2, FBN1, COL1A1, NR5A1, CCND2, CCNE1 and ALDH1A1. That FOXL2 was higher in the fibroblasts is at odds with our previous immunohistochemistry data [8]. Our previous microarray data demonstrated that GREL cells in culture had higher signal intensity of KRT19, DSG2 and PKP2 than cultured fibroblasts isolated from the adult ovary [8] and this is consistent our current on KRT19 and PKP2. There was a trend for higher DSG2 expression in GREL cells, but not significantly so. The elevated levels of NR2F2, FBN1 and COL1A1 in stromal fibroblasts is in agreement with previous immunolocalisation studies [8, 32]. Thus fetal fibroblasts and GREL cells are clearly different somatic cell types.

During development of the ovary the GREL cells, being the only somatic cells on the surface, mature into surface epithelial cells once the stroma has penetrated to below the surface and spread laterally below it [8, 9]. It is when the stroma underlies the cells on the surface that GREL cells form into a mature epithelium overlaying an epithelial basal lamina that separates the epithelium from the underlying stroma. The position of basal lamina and stroma on the basal side of the surface epithelial cells is important for the apical-basal polarity of epithelia and much has been written about the effects of stroma on adjacent epithelial cells both in organ development and in diseases processes such as cancer [33]. From the current study, the cells in the adult ovary most closely related in phenotype to the GREL cells were the surface epithelial cells. However, the adult surface epithelial cells had higher expression of OCLN, MUC1, LAMB2, NR5A1, ESR1, ESR2 and TGFBR3. DSG2, KRT19 and PKP2 were also highly expressed in GREL cells and ovarian surface epithelial cells. It has been found that cytokeratins, desmogleins and plakophilins bind to E-cadherin to facilitate formation of desmosomes [3436], suggesting that desmosomes exist in both, GREL cells and the ovarian surface epithelium. The higher expression of OCLN suggests that the epithelium on the surface of the adult ovary has matured with tight junctions between the cells and the higher expression of MUC1 suggests that they are actively and directionally secreting mucins. Mucins are expressed in epithelia contacting relatively harsh microenvironments and form a mucin layer on the apical surface of those epithelia [37, 38]. As such the adult surface epithelium would have apical-basal polarity cued by its interface with underlying stroma. The GREL cells in the developing ovary would not have developed that polarity until such time as the stroma has completed its penetration and expanded laterally below the surface. This is completed by the end of gestation [32].

GREL cells also have the bi-potential ability and develop into granulosa cells of the follicles [8]. They are the only somatic cells in the ovigerous cords. The other somatic cells of the ovary are outside the cords in stroma and are separated from the inside of the cords by a basal lamina. In the current study there were many genes that differed between the GREL cells and granulosa cells suggesting that the GREL cells would need to change considerably to form granulosa cells. It would not just be a process of gaining functions as GREL cells would need to reduce expression of KRT19, DSG2 and PKP2 and thus differentiating into granulosa cells may be more complex than the differentiation into ovarian surface epithelium. However, this might possibly not be the case as the comparison we made was with granulosa cells from mature antral follicles. By this stage the granulosa cells are replicating and forming a multilayered membrana granulosa and are hormonally active. In fact five of the genes differentially expressed between GREL cells and mature granulosa cells have been previously identified to be differentially expressed (COL4A1, INHBA, STAR, CYP11PA1 and CYP19A1), even from small antral to large antral follicle sizes [39]. In hindsight, a better comparison would have been with granulosa cells derived from primordial follicles, which we have never been successful at isolating.

The fetal ovarian fibroblasts were similar to adult fibroblasts in their gene expression profiles. However, in order to isolate them free from GREL cells we cultured them and compared them with adult fibroblasts which were also cultured. We have previously observed that ovarian fetal fibroblasts from fetuses in their first trimester of gestation express and secrete fibrillin 3 [8, 40], but they lose that ability rapidly in culture and instead express fibrillin 1 [26]. However, the cultured fetal fibroblasts did retain some fetal characteristics as TGFβ1 decreased expression of COL1A1 and COL3A1, as happens with fetal fibroblasts and which is the opposite of what adult fibroblasts do [41, 42]. We also examined the effects of 24 hormones, stimulators and growth factors and analysed the expression of genes expressed in fibroblasts (NR2F2, FBN1, COL1A1 and COL3A1). Only TGFβ1 altered the expression of three of these genes. We were also able to compare the expression patterns of genes in fibroblasts isolated in a previous experiment which were cultured with no added EGF [26] with those isolated in the current experiments isolated with EGF. The result showed that EGF significantly increased FBN1 expression, but decreased the gene expression of COL1A1. Thus fetal fibroblasts appear to have limited number of regulators.

In this study NR2F2, FBN1, COL1A1 and COL1A3 were used as the marker genes for the stromal fibroblasts in the ovary [30, 43, 44]. NR2F2 encodes COUP-TF2 and it is abundantly expressed in the human fetal ovary, playing a role as a sex-determining factor in the female gonad [45]. FBN1, COL1A1, COL1A3 which encode fibrillin 1, alpha 1 chain of collagen type I and alpha 1 chain of collagen type Ⅲ, respectively, are three components in ECM. Fibrillin 1 is a glycoprotein essential for the formation of microfibrils, which appear to provide a scaffold for the tissues [46]. Moreover, fibrillin networks exhibit regulatory functions on cell signalling and behaviour [47]. Collagens interact with fibres in ECM and regulate the cell morphology, growth and migration [48, 49], and collagen type I and III are secreted by fibroblasts [48, 50]. The high expression of NR2F2, FBN1 and COL1A1 in fetal fibroblasts was consistent with the higher expression in adult fibroblasts than ovarian surface epithelial cells and granulosa cells. LAMB2 encodes the beta 2 chain of laminin, which is a major component in basal laminas and essential to cell adhesion [51, 52]. Expression of LAMB2 and COL4A1 increased slightly across gestation suggesting that more basal lamina is present in late gestational fetal ovaries, which was consistent with the formation and growth of follicles and maturation of the ovarian surface epithelium.

In genes involved in steroidogenesis, NR5A1 encodes steroidogenic factor 1, and it could regulate genes encoding enzymes for steroidogenesis by binding to the promoter regions of genes encoding steroidogenesis enzymes, such as CYP11A1, CYP19A1 and HSD3B1 [53, 54]. STAR, also regulated by steroidogenic factor 1 [54], encodes steroidogenic acute regulatory protein, which is a transport protein for the transfer of cholesterol for the production of progesterone [55]. Our current data showed that NR5A1 was predominantly expressed in adult epithelial cells and granulosa cells, and this was consistent with the higher NR5A1 expression in GREL cells than fetal fibroblasts, suggesting that GREL cells exhibited a similar pattern with its adult counterparts in NR5A1 expression. CYP11A1, HSD3B1, CYP19A1 and HSD17B1 encode different enzymes in steroidogenesis, and they were predominantly expressed in granulosa cells. However, they were only minimally expressed in GREL cells.

Previous studies have demonstrated that oestrogens promote proliferation of ovarian surface epithelium [56] and granulosa cells [57]. ESR1 and ESR2 encode oestrogen receptor α and β, respectively. Our current data demonstrated that ESR1 was mostly expressed in the ovarian surface epithelial cells and ESR2 was predominantly expressed in granulosa cells, consistent with previous studies [58, 59]. The very low levels of expression of ESR1 and ESR2 in fetal and adult fibroblasts suggested that oestrogens do not directly regulate fibroblasts in ovarian development or function. Moreover, the weak expression of ESR1 and ESR2 in GREL cells suggested that GREL cells in culture are not regulated by oestrogens.

TGFβ signalling regulates processes such as follicle growth, oocyte meiotic maturation and luteinisation in the ovary [60, 61]. TGFBR1 and TGFBR2 encode TGFβ receptors 1 and 2, and have different affinities to TGFβ1, TGFβ2 and TGFβ3 [6264]. Unlike these receptors, the type 3 receptor, encoded by TGFBR3, does not have any kinase activity, and it functions as an accessory receptor stabilising the complexes between type 1 and 2 receptors and modulating the activity of TGFβ ligands [65, 66]. Compared with the fetal ovary, GREL cells and fetal fibroblasts in culture had higher expression levels of TGFBR1, however, the two fetal ovarian cells did not have higher expression of TGFBR2 or TGFBR3 than the fetal ovary. Moreover, the higher expression levels of TGFBR1, TGFBR2 and TGFBR3 in ovarian surface epithelial cells suggests that they could be more sensitive to ligands of the TGFβ family than granulosa cells. Similarly, this scenario also happens in other genes, such as KRT19, FBN1 and STAR, which implies the difference of granulosa cell from ovarian surface epithelial cell and GREL cell. Oocyte-derived factors, GDF9 and BMP15, regulate proliferation and differentiation of granulosa cells [23]. Therefore, oocyte derived factors would be more necessary in the differentiation of GREL cell into granulosa cell.

CCND2, CCNE1 and CCNE2 encode cyclin D2, cyclin E1 and cyclin E2, respectively, and activate cyclin-dependent kinases, which regulate cell cycle G1/S transition [67, 68] and increase cell proliferation [6971]. Our previous study demonstrated that the proliferating index was decreasing in both, the non-stromal component and the stroma of the fetal ovary, and the expression of CCND2 was decreasing across gestation [31, 32]. Consistent with the previous studies, CCNE1 and CCNE2 in the fetal ovary were also negatively correlated with the gestational age.

Aldehyde dehydrogenase 1 family member A1 is encoded by ALDH1A1 and functions as an enzyme in the biosynthesis of retinoic acid which regulates cell functions, such as self-renewal and differentiation. Furthermore, ALDH1A1 regulates both normal stem cells and tumor-initiating cells, and it is often used as a marker for stem cells [72, 73]. Leucine-rich repeat-containing G-protein coupled receptor 5 (LGR5) is encoded by LGR5 and it plays an important role in cell self-renewal [74, 75]. LGR5-positive stem cells have been identified in the mouse ovarian surface epithelium [21]. OCT4 is required for the survival of primordial germ cells [76], and its expression decreases at the onset of meiotic prophase I [77, 78]. Depending on the species, VASA and DAZL are essential for the survival or migration of primordial germ cell [79, 80]. The observed low expression of ALDH1A1, LGR5, OCT4, DAZL and VASA in GREL cells and fetal fibroblasts indicated that the two cell types probably do not have properties of stem/germ cells.

Conclusions

The two major somatic cell types of the developing ovary, GREL cells and fetal fibroblasts, have different cellular origins and exhibited different morphologies in cell culture. They also have distinct gene expression profiles. They also differ from the cells they ultimately differentiate in to. GREL cells had more epithelial properties which were closer to ovarian surface epithelial cells and less so to granulosa cells. Considering the regulation of the oocyte on granulosa cells, perhaps oocyte secreted factors might be necessary to establish the granulosa cell lineage. Fetal fibroblasts differed minimally from adult fibroblasts apart from higher expression of cell cycle genes. TGFβ1 and EGF may remodel ECM by regulating gene expression in fetal fibroblasts. In summary, whilst we can conclude that there are at least two somatic cells in the ovary with very different properties that also distinguish them from their adult counterparts, the regulation of cell fate determination, particularly of the bi-potentiality of GREL cells, remains to be elucidated.

Supporting information

S1 Fig. Behaviour of bovine fetal ovarian cells cultured on collagen type I coated plates.

After 2 days in culture, sporadic GREL cell clusters (A) and fetal fibroblasts (B) could be observed. GREL clusters grew slowly and formed larger clusters (C) in the subsequent culture days. GREL cells appeared polygonal, and in some areas tightly packed, as previously observed [8]. (D) Fetal fibroblasts grew faster than GREL cells and showed the typical spindle shaped appearance as seen in adult fibroblasts. (E) Alongside wells with pure GREL cells and pure fibroblasts, wells with both cell types could be observed. In these GREL cells formed clearly distinguishable clusters surrounded by fibroblasts. A sterile blunt-ended glass pipette was used to remove the fibroblasts and continue the culture to obtain GREL cells (F). Gestational ages were: (A) 51 and (B-F) 110 days. Bars: A-B, D-F = 100 μm; C = 200 μm.

(PDF)

S2 Fig. Representative micrograph images of each isolate of GREL cells used for qRT-PCR.

GREL cell clusters were grown on collagen type I coated plates. All 7 GREL cell samples collected for qRT-PCR were from passage 0. Gestational ages were (A) 51, (B) 56, (C) 73, (D, E) 110, (F) 127 and (G) 177 days. Bars: (A-G) = 100 μm.

(PDF)

S3 Fig. Representative micrograph images of each isolate of ovarian fetal fibroblast cells used for qRT-PCR.

Fetal fibroblasts were grown on collagen type I coated plates and collected on the same day as the corresponding GREL cells for gene expression analysis. All 6 fibroblast samples were from passage 0. Gestational ages were (A) 51, (B) 56, (C) 73, (D, E) 110, (F) 127 and (G) 177 days. Bars: (A-F) = 100 μm.

(PDF)

S4 Fig. Normalised gene expression of genes in GREL cells and fetal fibroblasts.

Gene expression levels in cultured GREL (n = 7) and fetal fibroblasts (n = 6) were normalised to PPIA and RPL32, and they are presented as scatter plots.

(PDF)

S5 Fig. Normalised gene expression of stromal genes in fetal fibroblasts in the presence of different treatments.

Fibroblasts from bovine fetal ovary (n = 5) were cultured with different treatments for 18 h. Fold changes of gene expression to the control groups are presented as mean ± SEM (normalised to PPIA and RPL32). Significantly different results among groups were determined by one-way ANOVA with Tukey’s post-hoc test.

(PDF)

S6 Fig. mRNA expression levels of stromal genes in fetal fibroblasts cultured in presence or absence of EGF.

Columns in white colour and black colour show gene expression in fetal fibroblasts cultured in the absence of EGF (n = 11) in a previous study [26] and in the presence of EGF (n = 7) in the current study, respectively. Data of gene expression are presented as mean ± SEM (normalised to PPIA and RPL32). Significantly different results between groups were determined by Student’s t-test. ** P < 0.01.

(PDF)

S1 Table. Treatments of cultured bovine fetal ovarian fibroblasts.

(PDF)

S2 Table. List of genes and primers used for qRT-PCR.

(PDF)

S3 Table. Pearson correlation coefficients (R) of expression levels of 10 newly-examined genes in bovine fetal ovaries and gestational age (n = 27).

(PDF)

Acknowledgments

We thank Thomas Foods International (Murray Bridge, SA), Strath Meats (Strathalbyn, SA) and The Midfield Group (Warrnambool, VIC) for providing bovine tissues and Mrs. Wendy M. Bonner for the collection of bovine fetal and adult ovaries.

Data Availability

All the data that pertain to the results are presented in the text or supplementary figures and tables.

Funding Statement

Funding was received from Adelaide University China Fee Scholarships (ML and QW; https://www.csc.edu.cn/chuguo), Adelaide Graduate Research Scholarship (NAB; http://www.adelaide.edu.au/), Australia Awards Scholarship (MDH; https://www.dfat.gov.au/people-to-people/australia-awards/australia-awards-scholarships), Robinson Research Institute Career Development Fellowship (KH; https://www.adelaide.edu.au/robinson-research-institute/) and Building On Ideas Grant (KH; https://www.adelaide.edu.au/robinson-research-institute/), National Health and Medical Research Council of Australia Centre for Research Excellence in the Evaluation, Management and Health Care Needs of Polycystic Ovary Syndrome (RJR; GTN1078444; https://www.nhmrc.gov.au/) and The Centre for Research Excellence on Women’s Health in Reproductive life (RJR and KH; GTN1171592; https://www.nhmrc.gov.au/) and the UK Medical Research Council (RAA; grant No. G1100357; https://mrc.ukri.org/). The funders did not play any role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

  • 1.Maheshwari A, Fowler PA. Primordial follicular assembly in humans—revisited. Zygote. 2008;16(4):285–96. doi: 10.1017/S0967199408004802 . [DOI] [PubMed] [Google Scholar]
  • 2.Wilhelm D, Palmer S, Koopman P. Sex determination and gonadal development in mammals. Physiol Rev. 2007;87(1):1–28. doi: 10.1152/physrev.00009.2006 . [DOI] [PubMed] [Google Scholar]
  • 3.Byskov AG, Lintern-Moore S. Follicle formation in the immature mouse ovary: the role of the rete ovarii. J Anat. 1973;116(2):207–17. . [PMC free article] [PubMed] [Google Scholar]
  • 4.Loffler S, Horn LC, Weber W, Spanel-Borowski K. The transient disappearance of cytokeratin in human fetal and adult ovaries. Anat Embryol (Berl). 2000;201(3):207–15. Epub 2000/02/09. doi: 10.1007/s004290050019 . [DOI] [PubMed] [Google Scholar]
  • 5.Byskov AG. Does the rete ovarii act as a trigger for the onset of meiosis? Nature. 1974;252(5482):396–7. Epub 1974/11/29. doi: 10.1038/252396a0 . [DOI] [PubMed] [Google Scholar]
  • 6.Sawyer HR, Smith P, Heath DA, Juengel JL, Wakefield SJ, McNatty KP. Formation of ovarian follicles during fetal development in sheep. Biol Reprod. 2002;66(4):1134–50. Epub 2002/03/22. doi: 10.1095/biolreprod66.4.1134 . [DOI] [PubMed] [Google Scholar]
  • 7.Juengel JL, Sawyer HR, Smith PR, Quirke LD, Heath DA, Lun S, et al. Origins of follicular cells and ontogeny of steroidogenesis in ovine fetal ovaries. Mol Cell Endocrinol. 2002;191(1):1–10. doi: 10.1016/s0303-7207(02)00045-x . [DOI] [PubMed] [Google Scholar]
  • 8.Hummitzsch K, Irving-Rodgers HF, Hatzirodos N, Bonner W, Sabatier L, Reinhardt DP, et al. A new model of development of the mammalian ovary and follicles. PLoS One. 2013;8(2):e55578. doi: 10.1371/journal.pone.0055578 ; PubMed Central PMCID: PMC3567121. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 9.Hartanti MD, Hummitzsch K, Bonner WM, Bastian NA, Irving-Rodgers HF, Rodgers RJ. Formation of the bovine ovarian surface epithelium during fetal development. J Histochem Cytochem. 2020;68(2):113–26. Epub 2019/12/20. doi: 10.1369/0022155419896797 ; PubMed Central PMCID: PMC7003494. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 10.Liu CF, Liu C, Yao HH. Building pathways for ovary organogenesis in the mouse embryo. Curr Top Dev Biol. 2010;90:263–90. doi: 10.1016/S0070-2153(10)90007-0 ; PubMed Central PMCID: PMC3400115. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 11.Heeren A, van Iperen L, Klootwijk DB, de Melo Bernardo A, Roost MS, Gomes Fernandes MM, et al. Development of the follicular basement membrane during human gametogenesis and early folliculogenesis. BMC Dev Biol. 2015;15(1):4. doi: 10.1186/s12861-015-0054-0 . [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 12.Kidder GM, Mhawi AA. Gap junctions and ovarian folliculogenesis. Reproduction. 2002;123(5):613–20. doi: 10.1530/rep.0.1230613 . [DOI] [PubMed] [Google Scholar]
  • 13.Carabatsos MJ, Sellitto C, Goodenough DA, Albertini DF. Oocyte-granulosa cell heterologous gap junctions are required for the coordination of nuclear and cytoplasmic meiotic competence. Dev Biol. 2000;226(2):167–79. doi: 10.1006/dbio.2000.9863 . [DOI] [PubMed] [Google Scholar]
  • 14.Shuhaibar LC, Egbert JR, Norris RP, Lampe PD, Nikolaev VO, Thunemann M, et al. Intercellular signaling via cyclic GMP diffusion through gap junctions restarts meiosis in mouse ovarian follicles. Proc Natl Acad Sci U S A. 2015;112(17):5527–32. doi: 10.1073/pnas.1423598112 ; PubMed Central PMCID: PMC4418852. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 15.Young JM, McNeilly AS. Theca: the forgotten cell of the ovarian follicle. Reproduction. 2010;140(4):489–504. doi: 10.1530/REP-10-0094 . [DOI] [PubMed] [Google Scholar]
  • 16.Liu C, Peng J, Matzuk MM, Yao HH. Lineage specification of ovarian theca cells requires multicellular interactions via oocyte and granulosa cells. Nat Commun. 2015;6:6934. Epub 2015/04/29. doi: 10.1038/ncomms7934 ; PubMed Central PMCID: PMC4413950. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 17.Zeleznik AJ. The physiology of follicle selection. Reprod Biol Endocrinol. 2004;2:31. doi: 10.1186/1477-7827-2-31 ; PubMed Central PMCID: PMC442133. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 18.Holesh JE, Bass AN, Lord M. Physiology, Ovulation. StatPearls. Treasure Island (FL) 2021. [PubMed] [Google Scholar]
  • 19.Mehta RV, Malcom PJ, Chang RJ. The effect of androgen blockade on granulosa cell estradiol production after follicle-stimulating hormone stimulation in women with polycystic ovary syndrome. J Clin Endocrinol Metab. 2006;91(9):3503–6. doi: 10.1210/jc.2006-0752 . [DOI] [PubMed] [Google Scholar]
  • 20.Kakuta H, Iguchi T, Sato T. The Involvement of Granulosa Cells in the Regulation by Gonadotropins of Cyp17a1 in Theca Cells. In Vivo. 2018;32(6):1387–401. doi: 10.21873/invivo.11391 ; PubMed Central PMCID: PMC6365763. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 21.Ng A, Tan S, Singh G, Rizk P, Swathi Y, Tan TZ, et al. Lgr5 marks stem/progenitor cells in ovary and tubal epithelia. Nature cell biology. 2014;16(8):745–57. Epub 2014/07/07. doi: 10.1038/ncb3000 . [DOI] [PubMed] [Google Scholar]
  • 22.Flesken-Nikitin A, Hwang CI, Cheng CY, Michurina TV, Enikolopov G, Nikitin AY. Ovarian surface epithelium at the junction area contains a cancer-prone stem cell niche. Nature. 2013;495(7440):241–5. Epub 2013/03/08. doi: 10.1038/nature11979 . [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 23.McNatty KP, Moore LG, Hudson NL, Quirke LD, Lawrence SB, Reader K, et al. The oocyte and its role in regulating ovulation rate: a new paradigm in reproductive biology. Reproduction. 2004;128(4):379–86. Epub 2004/09/30. doi: 10.1530/rep.1.00280 . [DOI] [PubMed] [Google Scholar]
  • 24.Hartanti MD, Hummitzsch K, Irving-Rodgers HF, Bonner WM, Copping KJ, Anderson RA, et al. Morphometric and gene expression analyses of stromal expansion during development of the bovine fetal ovary. Reprod Fertil Dev. 2019;31(3):482–95. doi: 10.1071/RD18218 . [DOI] [PubMed] [Google Scholar]
  • 25.Hornick JE, Duncan FE, Shea LD, Woodruff TK. Isolated primate primordial follicles require a rigid physical environment to survive and grow in vitro. Hum Reprod. 2012;27(6):1801–10. Epub 2012/03/30. doi: 10.1093/humrep/der468 ; PubMed Central PMCID: PMC3357191. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 26.Bastian NA, Bayne RA, Hummitzsch K, Hatzirodos N, Bonner WM, Hartanti MD, et al. Regulation of fibrillins and modulators of TGFbeta in fetal bovine and human ovaries. Reproduction. 2016;152(2):127–37. Epub 2016/05/26. doi: 10.1530/REP-16-0172 . [DOI] [PubMed] [Google Scholar]
  • 27.Russe I. Oogenesis in cattle and sheep. Bibl Anat. 1983;24:77–92. Epub 1983/01/01. . [PubMed] [Google Scholar]
  • 28.Hartanti MD, Rosario R, Hummitzsch K, Bastian NA, Hatzirodos N, Bonner WM, et al. Could perturbed fetal development of the ovary contribute to the development of polycystic ovary syndrome in later life? PLOS ONE. 2020;15(2):e0229351. doi: 10.1371/journal.pone.0229351 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 29.Liu M, Hummitzsch K, Hartanti MD, Rosario R, Bastian NA, Hatzirodos N, et al. Analysis of expression of candidate genes for polycystic ovary syndrome in adult and fetal human and fetal bovine ovariesdagger. Biol Reprod. 2020;103(4):840–53. doi: 10.1093/biolre/ioaa119 . [DOI] [PubMed] [Google Scholar]
  • 30.Hatzirodos N, Hummitzsch K, Irving-Rodgers HF, Breen J, Perry VEA, Anderson RA, et al. Transcript abundance of stromal and thecal cell related genes during bovine ovarian development. PLoS One. 2019;14(3):e0213575. doi: 10.1371/journal.pone.0213575 eCollection 2019. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 31.Hummitzsch K, Hatzirodos N, Irving-Rodgers HF, Hartanti MD, Perry VE, Anderson RA, et al. Morphometic analyses and gene expression related to germ cells, GREL cells and granulosa cells during development of the bovine fetal ovary. Plos One. 2019;14:e0214130. doi: 10.1371/journal.pone.0214130 [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 32.Hartanti MD, Hummitzsch K, Irving-Rodgers HF, Bonner WM, KJ C, Anderson RA, et al. Morphometric and gene expression analyses of stromal expansion during bovine fetal ovarian development. Reprod Fert Devel. 2019;31(3):482–95. doi: 10.1071/RD18218 [DOI] [PubMed] [Google Scholar]
  • 33.Cunha GR. Mesenchymal-epithelial interactions: past, present, and future. Differentiation. 2008;76(6):578–86. Epub 2008/06/19. doi: 10.1111/j.1432-0436.2008.00290.x . [DOI] [PubMed] [Google Scholar]
  • 34.Moch M, Schwarz N, Windoffer R, Leube RE. The keratin-desmosome scaffold: pivotal role of desmosomes for keratin network morphogenesis. Cell Mol Life Sci. 2020;77(3):543–58. doi: 10.1007/s00018-019-03198-y ; PubMed Central PMCID: PMC7010626. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 35.Shafraz O, Rubsam M, Stahley SN, Caldara AL, Kowalczyk AP, Niessen CM, et al. E-cadherin binds to desmoglein to facilitate desmosome assembly. eLife. 2018;7. doi: 10.7554/eLife.37629 ; PubMed Central PMCID: PMC6066328. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 36.Gosavi P, Kundu ST, Khapare N, Sehgal L, Karkhanis MS, Dalal SN. E-cadherin and plakoglobin recruit plakophilin3 to the cell border to initiate desmosome assembly. Cell Mol Life Sci. 2011;68(8):1439–54. doi: 10.1007/s00018-010-0531-3 . [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 37.Chauhan SC, Kumar D, Jaggi M. Mucins in ovarian cancer diagnosis and therapy. Journal of ovarian research. 2009;2:21. doi: 10.1186/1757-2215-2-21 ; PubMed Central PMCID: PMC2804676. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 38.Gipson IK, Spurr-Michaud S, Tisdale A, Menon BB. Comparison of the transmembrane mucins MUC1 and MUC16 in epithelial barrier function. PLoS One. 2014;9(6):e100393. doi: 10.1371/journal.pone.0100393 ; PubMed Central PMCID: PMC4072602. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 39.Hatzirodos N, Irving-Rodgers HF, Hummitzsch K, Harland ML, Morris SE, Rodgers RJ. Transcriptome profiling of granulosa cells of bovine ovarian follicles during growth from small to large antral sizes. BMC Genomics. 2014;15(1):24. doi: 10.1186/1471-2164-15-24 . [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 40.Prodoehl MJ, Irving-Rodgers HF, Bonner WM, Sullivan TM, Micke GC, Gibson MA, et al. Fibrillins and latent TGFbeta binding proteins in bovine ovaries of offspring following high or low protein diets during pregnancy of dams. Mol Cell Endocrinol. 2009;307(1–2):133–41. doi: 10.1016/j.mce.2009.03.002 . [DOI] [PubMed] [Google Scholar]
  • 41.Rolfe KJ, Irvine LM, Grobbelaar AO, Linge C. Differential gene expression in response to transforming growth factor-beta1 by fetal and postnatal dermal fibroblasts. Wound Repair Regen. 2007;15(6):897–906. doi: 10.1111/j.1524-475X.2007.00314.x . [DOI] [PubMed] [Google Scholar]
  • 42.Rolfe KJ, Grobbelaar AO. A review of fetal scarless healing. ISRN Dermatol. 2012;2012:698034. doi: 10.5402/2012/698034 ; PubMed Central PMCID: PMC3362931. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 43.Fan X, Bialecka M, Moustakas I, Lam E, Torrens-Juaneda V, Borggreven NV, et al. Single-cell reconstruction of follicular remodeling in the human adult ovary. Nat Commun. 2019;10(1):3164. doi: 10.1038/s41467-019-11036-9 ; PubMed Central PMCID: PMC6639403. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 44.Wagner M, Yoshihara M, Douagi I, Damdimopoulos A, Panula S, Petropoulos S, et al. Single-cell analysis of human ovarian cortex identifies distinct cell populations but no oogonial stem cells. Nat Commun. 2020;11(1):1147. doi: 10.1038/s41467-020-14936-3 ; PubMed Central PMCID: PMC7052271. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 45.Bashamboo A, Eozenou C, Jorgensen A, Bignon-Topalovic J, Siffroi JP, Hyon C, et al. Loss of Function of the Nuclear Receptor NR2F2, Encoding COUP-TF2, Causes Testis Development and Cardiac Defects in 46,XX Children. Am J Hum Genet. 2018;102(3):487–93. doi: 10.1016/j.ajhg.2018.01.021 ; PubMed Central PMCID: PMC5985285. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 46.Thomson J, Singh M, Eckersley A, Cain SA, Sherratt MJ, Baldock C. Fibrillin microfibrils and elastic fibre proteins: Functional interactions and extracellular regulation of growth factors. Seminars in cell & developmental biology. 2019;89:109–17. doi: 10.1016/j.semcdb.2018.07.016 ; PubMed Central PMCID: PMC6461133. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 47.Schrenk S, Cenzi C, Bertalot T, Conconi MT, Di Liddo R. Structural and functional failure of fibrillin1 in human diseases (Review). Int J Mol Med. 2018;41(3):1213–23. doi: 10.3892/ijmm.2017.3343 . [DOI] [PubMed] [Google Scholar]
  • 48.Frantz C, Stewart KM, Weaver VM. The extracellular matrix at a glance. J Cell Sci. 2010;123(Pt 24):4195–200. doi: 10.1242/jcs.023820 ; PubMed Central PMCID: PMC2995612. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 49.Fassett J, Tobolt D, Hansen LK. Type I collagen structure regulates cell morphology and EGF signaling in primary rat hepatocytes through cAMP-dependent protein kinase A. Mol Biol Cell. 2006;17(1):345–56. doi: 10.1091/mbc.e05-09-0871 ; PubMed Central PMCID: PMC1345672. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 50.Mariggio MA, Cassano A, Vinella A, Vincenti A, Fumarulo R, Lo Muzio L, et al. Enhancement of fibroblast proliferation, collagen biosynthesis and production of growth factors as a result of combining sodium hyaluronate and aminoacids. Int J Immunopathol Pharmacol. 2009;22(2):485–92. doi: 10.1177/039463200902200225 . [DOI] [PubMed] [Google Scholar]
  • 51.Aumailley M. The laminin family. Cell Adh Migr. 2013;7(1):48–55. doi: 10.4161/cam.22826 ; PubMed Central PMCID: PMC3544786. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 52.Hohenester E. Structural biology of laminins. Essays Biochem. 2019;63(3):285–95. doi: 10.1042/EBC20180075 ; PubMed Central PMCID: PMC6744579. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 53.Parker KL, Schimmer BP. Steroidogenic factor 1: a key determinant of endocrine development and function. Endocr Rev. 1997;18(3):361–77. doi: 10.1210/edrv.18.3.0301 . [DOI] [PubMed] [Google Scholar]
  • 54.Meinsohn MC, Smith OE, Bertolin K, Murphy BD. The Orphan Nuclear Receptors Steroidogenic Factor-1 and Liver Receptor Homolog-1: Structure, Regulation, and Essential Roles in Mammalian Reproduction. Physiol Rev. 2019;99(2):1249–79. doi: 10.1152/physrev.00019.2018 . [DOI] [PubMed] [Google Scholar]
  • 55.Payne AH, Hales DB. Overview of steroidogenic enzymes in the pathway from cholesterol to active steroid hormones. Endocr Rev. 2004;25(6):947–70. doi: 10.1210/er.2003-0030 . [DOI] [PubMed] [Google Scholar]
  • 56.Bai W, Oliveros-Saunders B, Wang Q, Acevedo-Duncan ME, Nicosia SV. Estrogen stimulation of ovarian surface epithelial cell proliferation. In Vitro Cell Dev Biol Anim. 2000;36(10):657–66. doi: . [DOI] [PubMed] [Google Scholar]
  • 57.Robker RL, Richards JS. Hormone-induced proliferation and differentiation of granulosa cells: a coordinated balance of the cell cycle regulators cyclin D2 and p27Kip1. Mol Endocrinol. 1998;12(7):924–40. Epub 1998/07/11. doi: 10.1210/mend.12.7.0138 . [DOI] [PubMed] [Google Scholar]
  • 58.Wright JW, Jurevic L, Stouffer RL. Dynamics of the primate ovarian surface epithelium during the ovulatory menstrual cycle. Hum Reprod. 2011;26(6):1408–21. doi: 10.1093/humrep/der057 ; PubMed Central PMCID: PMC3096557. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 59.Chu S, Nishi Y, Yanase T, Nawata H, Fuller PJ. Transrepression of estrogen receptor beta signaling by nuclear factor-kappab in ovarian granulosa cells. Mol Endocrinol. 2004;18(8):1919–28. doi: 10.1210/me.2004-0021 . [DOI] [PubMed] [Google Scholar]
  • 60.Yu C, Zhou JJ, Fan HY. Studying the Functions of TGF-beta Signaling in the Ovary. Methods in molecular biology (Clifton, NJ). 2016;1344:301–11. doi: 10.1007/978-1-4939-2966-5_19 . [DOI] [PubMed] [Google Scholar]
  • 61.Yang J, Zhang Y, Xu X, Li J, Yuan F, Bo S, et al. Transforming growth factor-beta is involved in maintaining oocyte meiotic arrest by promoting natriuretic peptide type C expression in mouse granulosa cells. Cell Death Dis. 2019;10(8):558. doi: 10.1038/s41419-019-1797-5 ; PubMed Central PMCID: PMC6646305. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 62.Qian SW, Burmester JK, Tsang ML, Weatherbee JA, Hinck AP, Ohlsen DJ, et al. Binding affinity of transforming growth factor-beta for its type II receptor is determined by the C-terminal region of the molecule. J Biol Chem. 1996;271(48):30656–62. doi: 10.1074/jbc.271.48.30656 . [DOI] [PubMed] [Google Scholar]
  • 63.De Crescenzo G, Pham PL, Durocher Y, O’Connor-McCourt MD. Transforming growth factor-beta (TGF-beta) binding to the extracellular domain of the type II TGF-beta receptor: receptor capture on a biosensor surface using a new coiled-coil capture system demonstrates that avidity contributes significantly to high affinity binding. J Mol Biol. 2003;328(5):1173–83. doi: 10.1016/s0022-2836(03)00360-7 . [DOI] [PubMed] [Google Scholar]
  • 64.Groppe J, Hinck CS, Samavarchi-Tehrani P, Zubieta C, Schuermann JP, Taylor AB, et al. Cooperative assembly of TGF-beta superfamily signaling complexes is mediated by two disparate mechanisms and distinct modes of receptor binding. Mol Cell. 2008;29(2):157–68. doi: 10.1016/j.molcel.2007.11.039 . [DOI] [PubMed] [Google Scholar]
  • 65.Vander Ark A, Cao J, Li X. TGF-beta receptors: In and beyond TGF-beta signaling. Cell Signal. 2018;52:112–20. doi: 10.1016/j.cellsig.2018.09.002 . [DOI] [PubMed] [Google Scholar]
  • 66.Kaivo-oja N, Jeffery LA, Ritvos O, Mottershead DG. Smad signalling in the ovary. Reprod Biol Endocrinol. 2006;4:21. doi: 10.1186/1477-7827-4-21 ; PubMed Central PMCID: PMC1459162. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 67.Quelle DE, Ashmun RA, Shurtleff SA, Kato JY, Bar-Sagi D, Roussel MF, et al. Overexpression of mouse D-type cyclins accelerates G1 phase in rodent fibroblasts. Genes & development. 1993;7(8):1559–71. doi: 10.1101/gad.7.8.1559 . [DOI] [PubMed] [Google Scholar]
  • 68.Sherr CJ, Roberts JM. CDK inhibitors: positive and negative regulators of G1-phase progression. Genes & development. 1999;13(12):1501–12. doi: 10.1101/gad.13.12.1501 . [DOI] [PubMed] [Google Scholar]
  • 69.Duronio RJ, Xiong Y. Signaling pathways that control cell proliferation. Cold Spring Harb Perspect Biol. 2013;5(3):a008904. doi: 10.1101/cshperspect.a008904 ; PubMed Central PMCID: PMC3578363. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 70.Keyomarsi K, Herliczek TW. The role of cyclin E in cell proliferation, development and cancer. Prog Cell Cycle Res. 1997;3:171–91. doi: 10.1007/978-1-4615-5371-7_14 . [DOI] [PubMed] [Google Scholar]
  • 71.Diril MK, Ratnacaram CK, Padmakumar VC, Du T, Wasser M, Coppola V, et al. Cyclin-dependent kinase 1 (Cdk1) is essential for cell division and suppression of DNA re-replication but not for liver regeneration. Proc Natl Acad Sci U S A. 2012;109(10):3826–31. doi: 10.1073/pnas.1115201109 ; PubMed Central PMCID: PMC3309725. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 72.Tomita H, Tanaka K, Tanaka T, Hara A. Aldehyde dehydrogenase 1A1 in stem cells and cancer. Onco_target. 2016;7(10):11018–32. doi: 10.18632/onco_target.6920 ; PubMed Central PMCID: PMC4905455. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 73.Vassalli G. Aldehyde Dehydrogenases: Not Just Markers, but Functional Regulators of Stem Cells. Stem Cells Int. 2019;2019:3904645. doi: 10.1155/2019/3904645 ; PubMed Central PMCID: PMC6348814. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 74.Birchmeier W. Stem cells: Orphan receptors find a home. Nature. 2011;476(7360):287–8. doi: 10.1038/476287a . [DOI] [PubMed] [Google Scholar]
  • 75.Van der Flier LG, Sabates-Bellver J, Oving I, Haegebarth A, De Palo M, Anti M, et al. The Intestinal Wnt/TCF Signature. Gastroenterology. 2007;132(2):628–32. doi: 10.1053/j.gastro.2006.08.039 . [DOI] [PubMed] [Google Scholar]
  • 76.Kehler J, Tolkunova E, Koschorz B, Pesce M, Gentile L, Boiani M, et al. Oct4 is required for primordial germ cell survival. EMBO Rep. 2004;5(11):1078–83. doi: 10.1038/sj.embor.7400279 ; PubMed Central PMCID: PMC1299174. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 77.Pesce M, Wang X, Wolgemuth DJ, Scholer H. Differential expression of the Oct-4 transcription factor during mouse germ cell differentiation. Mech Dev. 1998;71(1–2):89–98. doi: 10.1016/s0925-4773(98)00002-1 . [DOI] [PubMed] [Google Scholar]
  • 78.Yeom YI, Fuhrmann G, Ovitt CE, Brehm A, Ohbo K, Gross M, et al. Germline regulatory element of Oct-4 specific for the totipotent cycle of embryonal cells. Development. 1996;122(3):881–94. doi: 10.1242/dev.122.3.881 . [DOI] [PubMed] [Google Scholar]
  • 79.Li M, Hong N, Xu H, Yi M, Li C, Gui J, et al. Medaka vasa is required for migration but not survival of primordial germ cells. Mech Dev. 2009;126(5–6):366–81. doi: 10.1016/j.mod.2009.02.004 . [DOI] [PubMed] [Google Scholar]
  • 80.Zagore LL, Sweet TJ, Hannigan MM, Weyn-Vanhentenryck SM, Jobava R, Hatzoglou M, et al. DAZL Regulates Germ Cell Survival through a Network of PolyA-Proximal mRNA Interactions. Cell Rep. 2018;25(5):1225–40 e6. doi: 10.1016/j.celrep.2018.10.012 ; PubMed Central PMCID: PMC6878787. [DOI] [PMC free article] [PubMed] [Google Scholar]

Decision Letter 0

Rajakumar Anbazhagan

8 Feb 2022

PONE-D-21-29973Characterisation of bovine ovarian fetal fibroblasts and gonadal ridge epithelial-like cells and comparison to their adult counterpartsPLOS ONE

Dear Dr. Rodgers,

Thank you for submitting your manuscript to PLOS ONE. After careful consideration, we feel that it has merit but does not fully meet PLOS ONE’s publication criteria as it currently stands. Therefore, we invite you to submit a revised version of the manuscript that addresses the points raised during the review process.

Please submit your revised manuscript by Mar 25 2022 11:59PM. If you will need more time than this to complete your revisions, please reply to this message or contact the journal office at plosone@plos.org. When you're ready to submit your revision, log on to https://www.editorialmanager.com/pone/ and select the 'Submissions Needing Revision' folder to locate your manuscript file.

Please include the following items when submitting your revised manuscript:

  • A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). You should upload this letter as a separate file labeled 'Response to Reviewers'.

  • A marked-up copy of your manuscript that highlights changes made to the original version. You should upload this as a separate file labeled 'Revised Manuscript with Track Changes'.

  • An unmarked version of your revised paper without tracked changes. You should upload this as a separate file labeled 'Manuscript'.

If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter. Guidelines for resubmitting your figure files are available below the reviewer comments at the end of this letter.

If applicable, we recommend that you deposit your laboratory protocols in protocols.io to enhance the reproducibility of your results. Protocols.io assigns your protocol its own identifier (DOI) so that it can be cited independently in the future. For instructions see: https://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols. Additionally, PLOS ONE offers an option for publishing peer-reviewed Lab Protocol articles, which describe protocols hosted on protocols.io. Read more information on sharing protocols at https://plos.org/protocols?utm_medium=editorial-email&utm_source=authorletters&utm_campaign=protocols.

We look forward to receiving your revised manuscript.

Kind regards,

Academic Editor

PLOS ONE

Journal Requirements:

When submitting your revision, we need you to address these additional requirements.

1. Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at 

https://journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and 

https://journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf

2. Thank you for stating the following in the Acknowledgments/ Funding Section of your manuscript: 

Funding was received from Adelaide University China Fee Scholarships (ML and QW; https://www.csc.edu.cn/chuguo), Adelaide Graduate Research Scholarship (NAB; http:// www.adelaide.edu.au/), Australia Awards Scholarship (MDH; https://www.dfat.gov.au/people-to-people/australia-awards/australia-awards-scholarships), Robinson Research Institute Career Development Fellowship (KH; https://www.adelaide.edu.au/robinson-research-institute/) and Building On Ideas Grant (KH; https://www.adelaide.edu.au/robinson-research-institute/), National Health and Medical Research Council of Australia Centre for Research Excellence in the Evaluation, Management and Health Care Needs of Polycystic Ovary Syndrome (RJR; GTN1078444; https://www. nhmrc.gov.au/) and The Centre for Research Excellence on Women’s Health in Reproductive life (RJR and KH; GTN1171592; https://www. nhmrc.gov.au/) and the UK Medical Research Council (RAA; grant No. G1100357; https://mrc.ukri.org/). The funders did not play any role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Please note that funding information should not appear in the Acknowledgments section or other areas of your manuscript. We will only publish funding information present in the Funding Statement section of the online submission form. 

Please remove any funding-related text from the manuscript and let us know how you would like to update your Funding Statement. Currently, your Funding Statement reads as follows: 

Funding was received from Adelaide University China Fee Scholarships (ML and QW; https://www.csc.edu.cn/chuguo), Adelaide Graduate Research Scholarship (NAB; http:// www.adelaide.edu.au/), Australia Awards Scholarship (MDH; https://www.dfat.gov.au/people-to-people/australia-awards/australia-awards-scholarships), Robinson Research Institute Career Development Fellowship (KH; https://www.adelaide.edu.au/robinson-research-institute/) and Building On Ideas Grant (KH; https://www.adelaide.edu.au/robinson-research-institute/), National Health and Medical Research Council of Australia Centre for Research Excellence in the Evaluation, Management and Health Care Needs of Polycystic Ovary Syndrome (RJR; GTN1078444; https://www. nhmrc.gov.au/) and The Centre for Research Excellence on Women’s Health in Reproductive life (RJR and KH; GTN1171592; https://www. nhmrc.gov.au/) and the UK Medical Research Council (RAA; grant No. G1100357; https://mrc.ukri.org/). The funders did not play any role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Please include your amended statements within your cover letter; we will change the online submission form on your behalf. 

3. Thank you for stating the following in the Competing Interests section: 

RAA reports consultancy work for Ferring, Merck, IBSA, Roche Diagnostics, NeRRe Therapeutics and Sojournix Inc. The other authors of this manuscript have nothing to declare and no conflict of interest that could be perceived as prejudicing the impartiality of the research reported.

Please confirm that this does not alter your adherence to all PLOS ONE policies on sharing data and materials, by including the following statement: "This does not alter our adherence to  PLOS ONE policies on sharing data and materials.” (as detailed online in our guide for authors http://journals.plos.org/plosone/s/competing-interests).  If there are restrictions on sharing of data and/or materials, please state these. Please note that we cannot proceed with consideration of your article until this information has been declared. 

Please include your updated Competing Interests statement in your cover letter; we will change the online submission form on your behalf.

[Note: HTML markup is below. Please do not edit.]

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

Reviewer #4: Partly

**********

2. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

Reviewer #4: Yes

**********

3. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

Reviewer #4: No

**********

4. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

Reviewer #4: Yes

**********

5. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: In this manuscript “Characterisation of bovine ovarian fetal fibroblasts and gonadal ridge epithelial-like cells and comparison to their adult counterparts' ' reports importance of cell types in ovary development. This study can be helpful to have more knowledge on ovary development and ovary diseases. Article has details of methods and results, which are presented well in the manuscript. This article has made careful and further analysis of ovary development from the previous Plos One articles from these authors. This topic seems to help readers connect with the author's previous study in the features cell types in ovary development. I would recommend publication of this study.

Minor comment:

In line 496 authors should correct “week” to weak.

Authors can add a flowchart how to collect and maintain fetal fibroblasts and GREL1 cells from ovary tissues.

Reviewer #2: Reviewer’s comment

In the study titled ‘Characterization of bovine ovarian fetal fibroblasts and gonadal ridge epithelial-like cells and comparison to their adult counterparts’ the authors Liu et al isolated fibroblasts and gonadal ridge epithelial-like cells (GREL) from fetal ovary of pregnant cow. They have also isolated fibroblasts, surface epithelial cells, and granulosa cells from adult ovaries from non-pregnant cows. The culture of cells was done without passage, and based on the morphology, the contaminant cell if any, were removed by scraping. The cells were harvested and frozen for gene expression studies. To identify possible regulators of extracellular matrix, the fetal ovary fibroblasts were treated with hormones, stimulators, and growth factors. The expression profile of various genes was studied.

The two somatic cells, GREL cells and fetal fibroblasts, were distinguishable in terms of cellular origins, morphology, and gene expression profile in cell culture. These cells are also distinguishable from their adult counterparts. TGFβ1 and EGF may remodel ECM by regulating gene expression in fetal fibroblasts.

The article is a preliminary work in understanding the development of the ovary. The authors have put a good effort in preparing the manuscript. However, for the readers, it may a bit complicated. Hence need to be simplified. The following suggestions might increase the value of the manuscript.

1. The authors may consider changing the title to have a clear differentiation from their previous publication (Hummitzsch et al. 2019. Morphometric analyses and gene expression related to germ cells, GREL cells, and granulosa cells during development of the bovine fetal ovary). The suggested new title is ‘Isolation, culture, and characterization of bovine fetal ovarian fibroblasts and gonadal ridge epithelial-like cells compared to their adult counterparts.

2. The authors may consider simplifying the abstract. As in the discussion, group names of the genes may be mentioned.

3. Line 53: In culture, the cell types were morphologically different. - Rephrase the sentence. ‘In culture, these two cell types were morphologically different.

4. Published studies from the literature on isolating the ovarian cells are missing in the introduction. Studies in other species also can be included.

5. The methodology can be explained with reasons such that readers will understand the purpose.

6. The basis of selecting the genes can be explained in methods.

7. More details on samples collection can be included in methods (how many days of pregnancy, age of cow etc), although some of them have been explained in figure /table legends.

8. Line 126-137 After digesting with collagenase, how do you remove the undigested cell debris? A low-speed spin first or filtration?

9. PCR is a sensitive technique. What are the precautions taken to prevent fetal ovary fibroblasts and GREL cells from getting contaminated with other cells?

10. Any specific advantage of washing the tissues with 70% alcohol over washing it with media with antibiotics and antifungals?

11. Fig S4: X-axes title to be included.

12. Fig S5: Is there any statistically significant expression? If yes, the significance can be indicated by asterisks* as in Fig S6.

13. Figure S1 can be moved to the main body, or a representative figure of various cells may be included.

Reviewer #3: Characterisation of bovine ovarian fetal fibroblasts and gonadal ridge epithelial like cells and comparison to their adult counterparts

Menghe Liu et al

This paper seeks to isolate and characterise fibroblasts and GREL cells, the 2 major components of the fetal ovarian stroma. In the context of ovarian development, the developmental origins of health and disease (notably PCOS) the work is particularly relevant and important and, perhaps understated somewhat in the introduction.

The work is of a good standard, the results are of interest and build on our knowledge in this area. The manuscript is in general well crafted and put together. There are occasional lapses in grammar but these do not impact on the message.

Introduction

Lines 84-88. Give the discussion later around GREL cell differentiation into surface epithelial cells and granulosa cells I believe some mention or reference to the oocytes impacts of granulosa cell function/development would be appropriate (e.g. McNatty et al, Reproduction vol 128, 2004).

Minor issues

Lines94-95. Grammar could be addressed. Remove into or replace with thereby. ?

Line 105. A little understated I would have thought

Materials and methods. What I have struggled to find is the number of fetuses that were collected for the isolation of fibroblasts and GREL cells and importantly the age range of these fetuses. This should be clearly laid out in the M&M somewhere or in the figure legends.

Line 141 Grammar could be addressed here.

Lines 157-171. Isolation of adult ovarian fibroblasts. It is unclear to me and would be interesting to know whether the tissue to be minced and subsequently cultured contained any follicles

Line 254-255. While some of this data has been previously published it is my opinion that inclusion of this data is warranted

Presentation of gene expression data. Can any inference be drawn from the differences in expression levels between fetal ovaries and the isolated fetal cells, could this be an “artifact” of culture?. It should be acknowledged that for the adult cells fibroblasts have been cultured whereas both granulosa cells and OSE have not. Therefore differences between fibroblasts and granulosa or OSE may be attributed the culture of the cells. The same could be argued for differences between GREL cells and granulosa or OSE. If the authors do not accept this then their reason should be outlined in the discussion.

Line 268 and 272. it is OCLN not OCLD?

Lines 301-302. At the low levels of expression seen is the data reliable enough to support a significant difference, has the PCR been validated at these levels of expression.

Lines 407-408.Would the authors care to speculate of the difference in FOXL2 expression between this work and the previous IHC work. Could this be a function of culturing cells?.

Lines 421-422 and 438-440. This would not be surprising given that the granulosa cells have been under the influence of oocytes for some time.

Line 478. This would be consistent with the formation and growth of follicles

Line 497. The caveat with this statement ( and a number of others in the discussion) is that we are talking about expression in cultured cells. This should really be mentioned in the discussion.

Lines 503-505. This scenario is true for other genes as well (eg FBN1, STAR and KRT19). What is the implication of this pattern of expression.

Lines 533-534. I expected to see more around this statement in the discussion.

Figures and Legends. These are acceptable, well formatted and convey the intended message well. The exception is Supplementary Fig 1C which has poor uneven lighting and should be replaced

Discussion

Line 279. GJA1 was higher in granulosa cells not less.

Line 330-331. It would appear that for CCNE1 and CCNE2 the pattern of expression over gestation does not necessarily follow a linear pattern. The possibility should be acknowledged in the results or addressed with a more complex analysis.

Line 335. The difference in CCNE2 expression between adult OSE and adult fibroblasts is not indicated in the actual figure.

Reviewer #4: The manuscript "Characterisation of bovine ovarian fetal fibroblasts and gonadal ridge epithelial-like

cells and comparison to their adult counterparts" is basically a study where differential genes that are expressed in bovine fetal and adult tissue are compared after specific cell types have been cultured to isolate them.

General comments- A major problem with the manuscript is that it describes lists of genes different but does not explain why the differences in genes between each isolated cell population is important or what the different genes might confer to their functionality. Also if these genes are going to be utilized as markers of specific genes at specific ages (which would be a good enough function) then the processes they are involved in should be better described and there also needs to be better functional data to support these claims. Also the differences in the genes may also demonstrate why these cells are morphologically different and then how that helps their function.

As is-- the manuscript leaves the reviewer questioning what each gene is: As there are no genes that are defined and then what the function of the gene is--and how the group of genes that are differentially expressed--what does that mean? The authors appear to be demonstrating how different and similar fetal tissues are to adult but they need to develop a graphical abstract or diagram to demonstrate how these different fetal tissues may ultimately differentiate into adult--or what genes demonstrate different cell lineages (what genes are markers or are enriched for specific genes and that makes them have a specific lineage or type of cell). However, no identification of processes or function is made and the reader is left to see what genes are different but not really have a good idea of why these differences may be important.

Also no functional data other that treatment of TGFbeta and EGF is conducted. The outcome of these treatments are not fully incorporated into what the relevance of the differential gene expression observed means. Aspects of the paper are interesting but it does remain a characterization study throughout the paper and lacks some functional aspects that are necessary to give the paper more significance to its findings.

**********

6. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #2: No

Reviewer #3: No

Reviewer #4: No

[NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files.]

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com/. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email PLOS at figures@plos.org. Please note that Supporting Information files do not need this step.

PLoS One. 2022 Jul 8;17(7):e0268467. doi: 10.1371/journal.pone.0268467.r002

Author response to Decision Letter 0


22 Mar 2022

Note the line numbers refer to the track changes version of the manuscript.

Reviewer #1:

COMMENT: In line 496 authors should correct “week” to weak.

REPLY: Done.

COMMENT: Authors can add a flowchart how to collect and maintain fetal fibroblasts and GREL1 cells from ovary tissues.

REPLY: We prefer not to do this and keep the methods presented simply.

Reviewer #2:

COMMENT: 1. The authors may consider changing the title to have a clear differentiation from their previous publication (Hummitzsch et al. 2019. Morphometric analyses and gene expression related to germ cells, GREL cells, and granulosa cells during development of the bovine fetal ovary). The suggested new title is ‘Isolation, culture, and characterization of bovine fetal ovarian fibroblasts and gonadal ridge epithelial-like cells compared to their adult counterparts.

REPLY: Done.

COMMENT: 2. The authors may consider simplifying the abstract. As in the discussion, group names of the genes may be mentioned.

REPLY: We think listing gene names in the abstract would be helpful to the reader. No changes were made.

COMMENT: 3. Line 53: In culture, the cell types were morphologically different. - Rephrase the sentence. ‘In culture, these two cell types were morphologically different.

REPLY: Done (Line 52).

COMMENT: 4. Published studies from the literature on isolating the ovarian cells are missing in the introduction. Studies in other species also can be included.

REPLY: Whilst this is an important point it is something that could be done later in a methods review. To do this would potentially detract from the lines of reasoning that we have in our introduction.

COMMENT: 5. The methodology can be explained with reasons such that readers will understand the purpose.

REPLY: Done. We now say ‘We had to adopt this approach because we had not identified a specific antigen on surface of GREL cells to allow us to utilise more direct cell isolation methods such as the use of magnetic beads or fluorescence-activated cell sorting.’ (Line 399-402)

COMMENT: 6. The basis of selecting the genes can be explained in methods.

REPLY: Done. We now say “We analysed the expression levels of 34 genes relating to cell junctions, extracellular matrix, enzymes and transcription factors in order to characterise the cultured GREL cells and fetal fibroblasts” (Line 231-233).

COMMENT: 7. More details on samples collection can be included in methods (how many days of pregnancy, age of cow etc), although some of them have been explained in figure /table legends.

REPLY: We have now added the days of pregnancy into the relevant figure legends. (Pages 26-28)

COMMENT: 8. Line 126-137 After digesting with collagenase, how do you remove the undigested cell debris? A low-speed spin first or filtration?

REPLY: We now say “The cell suspensions were then passed through cell strainers and centrifuged at 1200 rpm for 5 min and the supernatants were discarded” (Line 142-143).

COMMENT: 9. PCR is a sensitive technique. What are the precautions taken to prevent fetal ovary fibroblasts and GREL cells from getting contaminated with other cells?

REPLY: When cDNA was synthesised, RNA was isolated from only one cell type at one time. When cDNA samples were added into the PCR systems, cDNA samples from different cell types was not added at one time.

COMMENT: 10. Any specific advantage of washing the tissues with 70% alcohol over washing it with media with antibiotics and antifungals?

REPLY: Since our samples were collected from the local abattoir, washing the tissue with 70% alcohol has a greater ability to reduce contamination in the subsequent cell culture.

COMMENT: 11. Fig S4: X-axes title to be included.

REPLY: The X-axes title has been added in Fig S4.

COMMENT: 12. Fig S5: Is there any statistically significant expression? If yes, the significance can be indicated by asterisks* as in Fig S6.

REPLY: There is no statistical significant differences, and this has been described in the results (Line 384).

COMMENT: 13. Figure S1 can be moved to the main body, or a representative figure of various cells may be included.

REPLY: We now include a representative figure (Fig 1)

Reviewer #3:

COMMENT: 1. Lines 84-88. Give the discussion later around GREL cell differentiation into surface epithelial cells and granulosa cells I believe some mention or reference to the oocytes impacts of granulosa cell function/development would be appropriate (e.g. McNatty et al, Reproduction vol 128, 2004).

REPLY: Now we say “The differences between the two cell types which GREL cells differentiate into could be the result of regulation by oocyte-secreted factors such as growth differentiation factor 9 (GDF9) and bone morphogenetic protein 15 (BMP15) which play essential roles in regulating the proliferation and differentiation of granulosa cells [23]. The ovarian surface epithelium, on the other hand, is separated from oocytes by stroma and basal lamina and is likely not exposed to oocyte-secreted factors” in Line 93-100.

COMMENT: 2. Lines 94-95. Grammar could be addressed. Remove into or replace with thereby. ?

REPLY: Done (Line 103-104).

COMMENT: 3. Line 105. A little understated I would have thought

REPLY: We now say ‘To conduct the first examination of GREL cells and to compare them to fetal fibroblasts,…..’ (Line 114-115)

COMMENT: 4. Materials and methods. What I have struggled to find is the number of fetuses that were collected for the isolation of fibroblasts and GREL cells and importantly the age range of these fetuses. This should be clearly laid out in the M&M somewhere or in the figure legends.

REPLY: The ages of fetuses that were collected for fibroblasts and GREL cells have been added in the figure legends. Pages 26, 27 and 28.

COMMENT: 5. Line 141 Grammar could be addressed here.

REPLY: Done. We now say “Then 250 cells were seeded per well in 24-well plates which were …” (Line 151-152).

COMMENT: 6. Lines 157-171. Isolation of adult ovarian fibroblasts. It is unclear to me and would be interesting to know whether the tissue to be minced and subsequently cultured contained any follicles.

REPLY: We have added some additional information to indicate that we avoided areas with follicles. We now say ‘Individual ovaries were cut into quarters and tissue areas devoid of visible follicles was cut to a depth of about 0.5 mm from the surface and minced with sterile razor blades…’

COMMENT: 7. Line 254-255. While some of this data has been previously published it is my opinion that inclusion of this data is warranted

REPLY: The data are shown in the figs. Its just the statement on line 268-269 that gave the wrong message and we have amended it.

COMMENT: 8. Presentation of gene expression data. Can any inference be drawn from the differences in expression levels between fetal ovaries and the isolated fetal cells, could this be an “artifact” of culture? It should be acknowledged that for the adult cells fibroblasts have been cultured whereas both granulosa cells and OSE have not. Therefore differences between fibroblasts and granulosa or OSE may be attributed the culture of the cells. The same could be argued for differences between GREL cells and granulosa or OSE. If the authors do not accept this then their reason should be outlined in the discussion.

REPLY: Yes we accept this comment. We now say in the opening paragraph of the Discussion ‘In this study we developed methods for isolating and culturing fetal ovarian GREL cells and fetal ovarian fibroblasts by removal of contaminating cells during culture and expansion from single colonies without passage. We had to adopt this approach because we had not identified a specific antigen on the surface of GREL cells to allow us to utilise more direct cell isolation methods such as the use of magnetic beads or fluorescence-activated cell sorting. We investigated the expression of a variety of genes to identify properties of GREL cells and fetal fibroblasts and examined differences with their adult counterparts. It should be acknowledged that culturing of GREL cells could have contributed to differences observed with the granulosa cells or ovarian surface epithelial cells which were not cultured’.

COMMENT: 8. Line 268 and 272. it is OCLN not OCLD?

REPLY: This error has been corrected throughout the manuscript.

COMMENT: 9. Lines 301-302. At the low levels of expression seen is the data reliable enough to support a significant difference, has the PCR been validated at these levels of expression.

REPLY: The Ct values of CYP11A1 and INHBA of the fetal fibroblast samples were all below 30. However, in the PCR examination of the blank controls of the two genes, Ct values were above 35 or no Ct values were obtained. Therefore, we believed that the two genes were expressed in fetal fibroblasts, even if they were in very low levels.

COMMENT: 10. Lines 407-408.Would the authors care to speculate of the difference in FOXL2 expression between this work and the previous IHC work. Could this be a function of culturing cells?.

REPLY: We are not sure why this is the case at this stage.

COMMENT: 11. Lines 421-422 and 438-440. This would not be surprising given that the granulosa cells have been under the influence of oocytes for some time.

REPLY: We agree with this opinion. The addition of sentences to the introduction address this issue directly (Line 93-100).

COMMENT: 12. Line 478. This would be consistent with the formation and growth of follicles

REPLY: We agree with this opinion, and this statement has been added in Line 501-502.

COMMENT: 13. Line 497. The caveat with this statement ( and a number of others in the discussion) is that we are talking about expression in cultured cells. This should really be mentioned in the discussion.

REPLY: This has been mentioned in Line 409, 431, 432, 521 and 527.

COMMENT: 14. Lines 503-505. This scenario is true for other genes as well (eg FBN1, STAR and KRT19). What is the implication of this pattern of expression.

REPLY: We think that the big differences between granulosa cell and other cell types (GREL cell, ovarian surface epithelial cell and fibroblast) imply the oocyte-derived regulation affects the differentiation from GREL cell into granulosa cell. Another statement has been added in the discussion in Line 531-535.

COMMENT: 15. Lines 533-534. I expected to see more around this statement in the discussion.

REPLY: We have added more statement in Line 531-535.

COMMENT: 16. Figures and Legends. These are acceptable, well formatted and convey the intended message well. The exception is Supplementary Fig 1C which has poor uneven lighting and should be replaced.

REPLY: Done.

COMMENT: 17. Line 279. GJA1 was higher in granulosa cells not less.

REPLY: This error has been corrected in Line 294.

COMMENT: 18. Line 330-331. It would appear that for CCNE1 and CCNE2 the pattern of expression over gestation does not necessarily follow a linear pattern. The possibility should be acknowledged in the results or addressed with a more complex analysis.

REPLY: We chose linear models so that we only identified strong relationships. There were many identified. No changes made.

COMMENT: 19. Line 335. The difference in CCNE2 expression between adult OSE and adult fibroblasts is not indicated in the actual figure.

REPLY: There is no statistically significant difference between adult ovarian surface epithelial cells and adult fibroblasts. We now say ‘….however, there was no statistically significant difference between adult ovarian surface epithelial cells and adult fibroblasts.’ (Line 353-354).

Reviewer #4:

COMMENT A major problem with the manuscript is that it describes lists of genes but does not explain why the differences in genes between each isolated cell population is important or what the different genes might confer to their functionality. Also if these genes are going to be utilized as markers of specific genes at specific ages (which would be a good enough function) then the processes they are involved in should be better described and there also needs to be better functional data to support these claims. Also the differences in the genes may also demonstrate why these cells are morphologically different and then how that helps their function.

As is the manuscript leaves the reviewer questioning what each gene is: As there are no genes that are defined and then what the function of the gene is--and how the group of genes that are differentially expressed--what does that mean? The authors appear to be demonstrating how different and similar fetal tissues are to adult but they need to develop a graphical abstract or diagram to demonstrate how these different fetal tissues may ultimately differentiate into adult--or what genes demonstrate different cell lineages (what genes are markers or are enriched for specific genes and that makes them have a specific lineage or type of cell). However, no identification of processes or function is made and the reader is left to see what genes are different but not really have a good idea of why these differences may be important.

REPLY: The selection of the genes emanates from our previous publication (Hummitzsch et al, 2013) where we first describe GREL cells, their morphology and some microarray data. We also chose genes known to be expressed in fibroblasts. In this study, we developed a method to isolate and culture GREL cells, and the result showed that these GREL cells were mostly consistent with our previous study. The fibroblasts expressed genes as expected. This this validates our cell isolation and then we compared the cells with themselves and with adult counterparts. Thus providing new knowledge.

COMMENT: Also no functional data other that treatment of TGFbeta and EGF is conducted. The outcome of these treatments are not fully incorporated into what the relevance of the differential gene expression observed means. Aspects of the paper are interesting but it does remain a characterization study throughout the paper and lacks some functional aspects that are necessary to give the paper more significance to its findings.

REPLY: We do not disagree that the study is mostly a characterization study. However, this is the first time to show that GREL cells are a real entity. They are very different somatic cell to fibroblasts and they are not granulosa cells nor surface epithelial cells. It is the first time they have been isolated and studied.

Attachment

Submitted filename: reply to reviewers v5.docx

Decision Letter 1

Rajakumar Anbazhagan

12 Apr 2022

PONE-D-21-29973R1Isolation, culture, and characterisation of bovine ovarian fetal fibroblasts and gonadal ridge epithelial-like cells and comparison to their adult counterpartsPLOS ONE

Dear Dr. Rodgers,

Thank you for submitting your manuscript to PLOS ONE. After careful consideration, we feel that it has merit but does not fully meet PLOS ONE’s publication criteria as it currently stands. Therefore, we invite you to submit a revised version (Minor Revision) of the manuscript that addresses the points raised during the review process.

Please submit your revised manuscript by May 27 2022 11:59PM. If you will need more time than this to complete your revisions, please reply to this message or contact the journal office at plosone@plos.org. When you're ready to submit your revision, log on to https://www.editorialmanager.com/pone/ and select the 'Submissions Needing Revision' folder to locate your manuscript file.

Please include the following items when submitting your revised manuscript:

  • A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). You should upload this letter as a separate file labeled 'Response to Reviewers'.

  • A marked-up copy of your manuscript that highlights changes made to the original version. You should upload this as a separate file labeled 'Revised Manuscript with Track Changes'.

  • An unmarked version of your revised paper without tracked changes. You should upload this as a separate file labeled 'Manuscript'.

If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter. Guidelines for resubmitting your figure files are available below the reviewer comments at the end of this letter.

If applicable, we recommend that you deposit your laboratory protocols in protocols.io to enhance the reproducibility of your results. Protocols.io assigns your protocol its own identifier (DOI) so that it can be cited independently in the future. For instructions see: https://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols. Additionally, PLOS ONE offers an option for publishing peer-reviewed Lab Protocol articles, which describe protocols hosted on protocols.io. Read more information on sharing protocols at https://plos.org/protocols?utm_medium=editorial-email&utm_source=authorletters&utm_campaign=protocols.

We look forward to receiving your revised manuscript.

Kind regards,

Academic Editor

PLOS ONE

Journal Requirements:

Please review your reference list to ensure that it is complete and correct. If you have cited papers that have been retracted, please include the rationale for doing so in the manuscript text, or remove these references and replace them with relevant current references. Any changes to the reference list should be mentioned in the rebuttal letter that accompanies your revised manuscript. If you need to cite a retracted article, indicate the article’s retracted status in the References list and also include a citation and full reference for the retraction notice.

[Note: HTML markup is below. Please do not edit.]

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. If the authors have adequately addressed your comments raised in a previous round of review and you feel that this manuscript is now acceptable for publication, you may indicate that here to bypass the “Comments to the Author” section, enter your conflict of interest statement in the “Confidential to Editor” section, and submit your "Accept" recommendation.

Reviewer #2: All comments have been addressed

Reviewer #3: All comments have been addressed

Reviewer #4: All comments have been addressed

Reviewer #5: (No Response)

**********

2. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #2: Yes

Reviewer #3: Yes

Reviewer #4: Yes

Reviewer #5: Partly

**********

3. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #2: Yes

Reviewer #3: Yes

Reviewer #4: Yes

Reviewer #5: Yes

**********

4. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #2: Yes

Reviewer #3: Yes

Reviewer #4: Yes

Reviewer #5: Yes

**********

5. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #2: Yes

Reviewer #3: Yes

Reviewer #4: Yes

Reviewer #5: No

**********

6. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #2: The author’s response and changes made in the manuscript are satisfactory. The manuscript is improved now.

Reviewer #3: I am comfortable that the questions raised have been addressed to a satisfactory level. While I agree with reviewer 4 that more functional studies would have been useful i believe there is enough new and interesting material to warrant publication as a largely characterization study.

Reviewer #4: All concerns by this reviewer were met by authors. While a graphical abstract would help tell the story this reviewer is not going to make it a requirement for publication.

Reviewer #5: Isolation, culture, and Characterization of bovine ovarian fetal fibroblasts and gonadal ridge epithelial-like cells and comparison to their counterparts

Main Theme of the manuscript:

In this article, the authors have isolated gonadal ridge epithelial-like cells (GREL) and examined the expression of various genes in GREL cells along with other counterpart cells such as bovine ovarian fetal fibroblasts, adult fibroblasts, granulosa cells (GC), and ovarian surface cells (OSE) by qRT-PCR. The gene expression analysis shows that GREL cells have a differential pattern of gene expression compared with other cells. I appreciate the author's attempts at the basic characterization of GREL cells for the first time. However, this article is lacking some details as mentioned below.

Major comments:

This manuscript is mainly focused on just technical details of isolation, culture, and gene expression analysis. The authors listed out the differential expression of a list of genes but didn’t explain what it could be the effect of these differential gene expressions make them unique in terms of functionality of the GREL cells. To justify the title of the manuscript, authors must fully characterize the cells by performing an analysis of cell lineage markers, cell doubling time, cell cycle analysis, etc., and compare them with other counterpart cells. This would demonstrate why these cells are morphologically and functionally different from other cells. Since it is just isolation, culture, and gene expression analysis, I feel it could be a short communication or a technical report instead of a full-length research article.

Minor comments:

1. There are a few minor typos and grammatical errors throughout. The author must check the grammatical errors thoroughly.

2. It would give a good look at the gene expression data if it was shown in a heat-map instead of bar charts. Just one glance would give a complete picture of up and down-regulated genes across the various cell types. For example, the heat-map shows down-regulated genes in green and up-regulated genes in red in the heat-map.

3. What is the reason for adding 3 ng/ml EGF while isolating GREL and fibroblast cells?

4. Why is the incubation temperature at 38.5 c? is there any specific reason to keep them at that temperature?

5. How pure these cells are? how did the authors make sure that they didn’t pick up other cells while scrapping off the colonies? since PCR is a very sensitive technique, authors may get the expression of the genes from the cross-contaminated cells as well. Did the authors include fibroblast lineage marker genes as a control to check the purity of the GREL cells in the PCR panel?

6. Did the authors check the genomic DNA contamination by including non-RT controls? If yes it must be included in the methods section.

7. It is interesting to see the effect of TGF-b, and EGF treatment on gene expression but the authors didn’t explain the functional outcome. The manuscript lacks explanations for functional aspects although they are significant findings.

8. X and Y axis titles must be present for every graph.

9. In the figure, S4 scatter plot labels are missing, what are the blue and orange dots represent? What does n mean? the number of colonies or number of samples? please mention it in the figure legend.

**********

7. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #2: No

Reviewer #3: No

Reviewer #4: No

Reviewer #5: Yes: Vijaya kumar Pidugu

[NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files.]

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com/. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email PLOS at figures@plos.org. Please note that Supporting Information files do not need this step.

PLoS One. 2022 Jul 8;17(7):e0268467. doi: 10.1371/journal.pone.0268467.r004

Author response to Decision Letter 1


20 Apr 2022

COMMENT: This manuscript is mainly focused on just technical details of isolation, culture, and gene expression analysis. The authors listed out the differential expression of a list of genes but didn’t explain what it could be the effect of these differential gene expressions make them unique in terms of functionality of the GREL cells. To justify the title of the manuscript, authors must fully characterize the cells by performing an analysis of cell lineage markers, cell doubling time, cell cycle analysis, etc., and compare them with other counterpart cells. This would demonstrate why these cells are morphologically and functionally different from other cells. Since it is just isolation, culture, and gene expression analysis, I feel it could be a short communication or a technical report instead of a full-length research article.

REPLY: It should be appreciated that this is the first time GREL cells have been isolated and cultured and as such this a substantial contribution to field of fetal ovary development. We examined expression of 34 genes and compared those with their counterparts. This is a necessary first step in understanding these cells. Hence at this time there are no lineage markers discovered much less are they available to be monitored.

If you feel it can be a short communication then we think the editor should decide this.

COMMENT 1. There are a few minor typos and grammatical errors throughout. The author must check the grammatical errors thoroughly.

REPLY: Thank you for the suggestion. We have proofed it again.

COMMENT 2. It would give a good look at the gene expression data if it was shown in a heat-map instead of bar charts. Just one glance would give a complete picture of up and down-regulated genes across the various cell types. For example, the heat-map shows down-regulated genes in green and up-regulated genes in red in the heat-map.

REPLY: Yes the heat maps could be useful, but we went to a lot of work to ensure we also had data on expression during the continuum of gestation and those data do not lend themselves readily to a heat map, hence we did not include one.

COMMENT 3. What is the reason for adding 3 ng/ml EGF while isolating GREL and fibroblast cells?

REPLY: We found that GREL cells, being an epithelial cell, needed EGF to maintain the morphological phenotype.

COMMENT 4. Why is the incubation temperature at 38.5 c? is there any specific reason to keep them at that temperature?

REPLY: The bovine body temperature is 38.5 C.

COMMENT 5. How pure these cells are? how did the authors make sure that they didn’t pick up other cells while scrapping off the colonies? since PCR is a very sensitive technique, authors may get the expression of the genes from the cross-contaminated cells as well. Did the authors include fibroblast lineage marker genes as a control to check the purity of the GREL cells in the PCR panel?

REPLY: Yes its possible that any contamination could lead to a PCR signal. However, we were using quantitative PCR and not using it for detection. So the comment about the sensitivity of PCR is not so relevant in this situation. Also Figs S2 and S3 show photos of each preparation of cells that were actually used in these experiments and clearly they are highly enriched for their respective cell types.

In any event we said ‘In this study NR2F2, FBN1, COL1A1 and COL1A3 were used as the marker genes for the stromal fibroblasts in the ovary’ and these were higher in fibroblasts than GREL cells. We are confident of our results.

COMMENT 6. Did the authors check the genomic DNA contamination by including non-RT controls? If yes it must be included in the methods section.

REPLY: Yes. We assumed this was a routine control that everyone did. In any event as we reported we used primers spanning introns. We now include this information in the methods (line 207).

COMMENT 7. It is interesting to see the effect of TGF-b, and EGF treatment on gene expression but the authors didn’t explain the functional outcome. The manuscript lacks explanations for functional aspects although they are significant findings.

REPLY: We said ‘TGFβ1 and EGF may remodel ECM by regulating gene expression in fetal fibroblasts.’ Given that we do not know about the in vivo expression of these molecules we feel in premature to speculate too much at this stage about the meaning of the results for the in vivo situation.

COMMENT 8. X and Y axis titles must be present for every graph.

REPLY: We prefer to have the X-axis labelled for the Figs as a whole not as individual parts of the Fig.

9. In the figure, S4 scatter plot labels are missing, what are the blue and orange dots represent? What does n mean? the number of colonies or number of samples? please mention it in the figure legend.

REPLY: These are now clarified in the legend.

Attachment

Submitted filename: Reply to reviewer.docx

Decision Letter 2

Rajakumar Anbazhagan

1 May 2022

Isolation, culture, and characterisation of bovine ovarian fetal fibroblasts and gonadal ridge epithelial-like cells and comparison to their adult counterparts

PONE-D-21-29973R2

Dear Dr. Rodgers,

We’re pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it meets all outstanding technical requirements.

Within one week, you’ll receive an e-mail detailing the required amendments. When these have been addressed, you’ll receive a formal acceptance letter and your manuscript will be scheduled for publication.

An invoice for payment will follow shortly after the formal acceptance. To ensure an efficient process, please log into Editorial Manager at http://www.editorialmanager.com/pone/, click the 'Update My Information' link at the top of the page, and double check that your user information is up-to-date. If you have any billing related questions, please contact our Author Billing department directly at authorbilling@plos.org.

If your institution or institutions have a press office, please notify them about your upcoming paper to help maximize its impact. If they’ll be preparing press materials, please inform our press team as soon as possible -- no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org.

Kind regards,

Academic Editor

PLOS ONE

Additional Editor Comments (optional):

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. If the authors have adequately addressed your comments raised in a previous round of review and you feel that this manuscript is now acceptable for publication, you may indicate that here to bypass the “Comments to the Author” section, enter your conflict of interest statement in the “Confidential to Editor” section, and submit your "Accept" recommendation.

Reviewer #2: All comments have been addressed

Reviewer #4: All comments have been addressed

Reviewer #5: All comments have been addressed

**********

2. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #2: Yes

Reviewer #4: Yes

Reviewer #5: Yes

**********

3. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #2: Yes

Reviewer #4: Yes

Reviewer #5: Yes

**********

4. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #2: Yes

Reviewer #4: Yes

Reviewer #5: Yes

**********

5. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #2: Yes

Reviewer #4: Yes

Reviewer #5: Yes

**********

6. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #2: Based on the reviewers comments the authors have made changes to the manuscript. Now the manuscript looks good.

The manuscript can be considered for publication.

Reviewer #4: (No Response)

Reviewer #5: (No Response)

**********

7. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #2: No

Reviewer #4: No

Reviewer #5: No

Acceptance letter

Rajakumar Anbazhagan

29 Jun 2022

PONE-D-21-29973R2

Isolation, culture, and characterisation of bovine ovarian fetal fibroblasts and gonadal ridge epithelial-like cells and comparison to their adult counterparts

Dear Dr. Rodgers:

I'm pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now with our production department.

If your institution or institutions have a press office, please let them know about your upcoming paper now to help maximize its impact. If they'll be preparing press materials, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information please contact onepress@plos.org.

If we can help with anything else, please email us at plosone@plos.org.

Thank you for submitting your work to PLOS ONE and supporting open access.

Kind regards,

PLOS ONE Editorial Office Staff

on behalf of

Dr. Rajakumar Anbazhagan

Academic Editor

PLOS ONE

Associated Data

    This section collects any data citations, data availability statements, or supplementary materials included in this article.

    Supplementary Materials

    S1 Fig. Behaviour of bovine fetal ovarian cells cultured on collagen type I coated plates.

    After 2 days in culture, sporadic GREL cell clusters (A) and fetal fibroblasts (B) could be observed. GREL clusters grew slowly and formed larger clusters (C) in the subsequent culture days. GREL cells appeared polygonal, and in some areas tightly packed, as previously observed [8]. (D) Fetal fibroblasts grew faster than GREL cells and showed the typical spindle shaped appearance as seen in adult fibroblasts. (E) Alongside wells with pure GREL cells and pure fibroblasts, wells with both cell types could be observed. In these GREL cells formed clearly distinguishable clusters surrounded by fibroblasts. A sterile blunt-ended glass pipette was used to remove the fibroblasts and continue the culture to obtain GREL cells (F). Gestational ages were: (A) 51 and (B-F) 110 days. Bars: A-B, D-F = 100 μm; C = 200 μm.

    (PDF)

    S2 Fig. Representative micrograph images of each isolate of GREL cells used for qRT-PCR.

    GREL cell clusters were grown on collagen type I coated plates. All 7 GREL cell samples collected for qRT-PCR were from passage 0. Gestational ages were (A) 51, (B) 56, (C) 73, (D, E) 110, (F) 127 and (G) 177 days. Bars: (A-G) = 100 μm.

    (PDF)

    S3 Fig. Representative micrograph images of each isolate of ovarian fetal fibroblast cells used for qRT-PCR.

    Fetal fibroblasts were grown on collagen type I coated plates and collected on the same day as the corresponding GREL cells for gene expression analysis. All 6 fibroblast samples were from passage 0. Gestational ages were (A) 51, (B) 56, (C) 73, (D, E) 110, (F) 127 and (G) 177 days. Bars: (A-F) = 100 μm.

    (PDF)

    S4 Fig. Normalised gene expression of genes in GREL cells and fetal fibroblasts.

    Gene expression levels in cultured GREL (n = 7) and fetal fibroblasts (n = 6) were normalised to PPIA and RPL32, and they are presented as scatter plots.

    (PDF)

    S5 Fig. Normalised gene expression of stromal genes in fetal fibroblasts in the presence of different treatments.

    Fibroblasts from bovine fetal ovary (n = 5) were cultured with different treatments for 18 h. Fold changes of gene expression to the control groups are presented as mean ± SEM (normalised to PPIA and RPL32). Significantly different results among groups were determined by one-way ANOVA with Tukey’s post-hoc test.

    (PDF)

    S6 Fig. mRNA expression levels of stromal genes in fetal fibroblasts cultured in presence or absence of EGF.

    Columns in white colour and black colour show gene expression in fetal fibroblasts cultured in the absence of EGF (n = 11) in a previous study [26] and in the presence of EGF (n = 7) in the current study, respectively. Data of gene expression are presented as mean ± SEM (normalised to PPIA and RPL32). Significantly different results between groups were determined by Student’s t-test. ** P < 0.01.

    (PDF)

    S1 Table. Treatments of cultured bovine fetal ovarian fibroblasts.

    (PDF)

    S2 Table. List of genes and primers used for qRT-PCR.

    (PDF)

    S3 Table. Pearson correlation coefficients (R) of expression levels of 10 newly-examined genes in bovine fetal ovaries and gestational age (n = 27).

    (PDF)

    Attachment

    Submitted filename: reply to reviewers v5.docx

    Attachment

    Submitted filename: Reply to reviewer.docx

    Data Availability Statement

    All the data that pertain to the results are presented in the text or supplementary figures and tables.


    Articles from PLoS ONE are provided here courtesy of PLOS

    RESOURCES

      NODES
    admin 4
    Idea 5
    idea 5
    INTERN 3
    Note 16
    Project 4
    twitter 2