Search

During maturation oocytes at the germinal vesicle (GV) stage progress to metaphase II (MII). However, during in vitro maturation a proportion often fail to progress. To understand these processes, we employed RNA sequencing to examine the transcriptome profile of these three groups of oocytes from the pig. We compared our findings with similar public oocyte data from humans. The transcriptomes in oocytes that failed to progress was similar to those that did. We found in both species, the most upregulated genes in MII oocytes were associated with chromosome segregation and cell cycle processes, while the most down regulated genes were relevant to ribosomal and mitochondrial pathways. Moreover, those genes involved in chromosome segregation during GV to MII transition were conserved in pig and human. We also compared MII and GV oocyte transcriptomes at the isoform transcript level in both species. Several thousands of genes (including DTNBP1, MAPK1, RAB35, GOLGA7, ATP1A1 and ATP2B1) identified as not different in expression at a gene transcript level were found to have differences in isoform transcript levels. Many of these genes were involved in ATPase-dependent or GTPase-dependent intracellular transport in pig and human, respectively. In conclusion, our study suggests the failure to progress to MII in vitro may not be regulated at the level of the genome and that many genes are differentially regulated at the isoform level, particular those involved ATPase- or GTPase-dependent intracellular transport., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Tang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The ovarian KGN granulosa-like tumour cell line is commonly used as a model for human granulosa cells, especially since it produces steroid hormones. To explore this further, we identified genes that were differentially expressed by KGN cells compared to primary human granulosa cells using three public RNA sequence datasets. Of significance, we identified that the expression of the antioxidant gene TXNRD1 (thioredoxin reductase 1) was extremely high in KGN cells. This is ominous since cytochrome P450 enzymes leak electrons and produce reactive oxygen species during the biosynthesis of steroid hormones. Gene Ontology (GO) analysis identified steroid biosynthetic and cholesterol metabolic processes were more active in primary granulosa cells, whilst in KGN cells, DNA processing, chromosome segregation and kinetochore pathways were more prominent. Expression of cytochrome P450 cholesterol side-chain cleavage (CYP11A1) and cytochrome P450 aromatase (CYP19A1), which are important for the biosynthesis of the steroid hormones progesterone and oestrogen, plus their electron transport chain members (FDXR, FDX1, POR) were measured in cultured KGN cells. KGN cells were treated with 1 mM dibutyryl cAMP (dbcAMP) or 10 μM forskolin, with or without siRNA knockdown of TXNRD1. We also examined expression of antioxidant genes, H2O2 production by Amplex Red assay and DNA damage by γH2Ax staining. Significant increases in CYP11A1 and CYP19A1 were observed by either dbcAMP or forskolin treatments. However, no significant changes in H2O2 levels or DNA damage were found. Knockdown of expression of TXNRD1 by siRNA blocked the stimulation of expression of CYP11A1 and CYP19A1 by dbcAMP. Thus, with TXNRD1 playing such a pivotal role in steroidogenesis in the KGN cells and it being so highly overexpressed, we conclude that KGN cells might not be the most appropriate model of primary granulosa cells for studying the interplay between ovarian steroidogenesis, reactive oxygen species and antioxidants., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Tang et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

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., Competing Interests: 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. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

Polycystic ovary syndrome (PCOS) affects around 10% of young women, with adverse consequences on fertility and cardiometabolic outcomes. PCOS appears to result from a genetic predisposition interacting with developmental events during fetal or perinatal life. We hypothesised that PCOS candidate genes might be expressed in the fetal ovary when the stroma develops; mechanistically linking the genetics, fetal origins and adult ovarian phenotype of PCOS. In bovine fetal ovaries (n = 37) of 18 PCOS candidate genes only SUMO1P1 was not expressed. Three patterns of expression were observed: early gestation (FBN3, GATA4, HMGA2, TOX3, DENND1A, LHCGR and FSHB), late gestation (INSR, FSHR, and LHCGR) and throughout gestation (THADA, ERBB4, RAD50, C8H9orf3, YAP1, RAB5B, SUOX and KRR1). A splice variant of FSHB exon 3 was also detected early in the bovine ovaries, but exon 2 was not detected. Three other genes, likely to be related to the PCOS aetiology (AMH, AR and TGFB1I1), were also expressed late in gestation. Significantly within each of the three gene groups, the mRNA levels of many genes were highly correlated with each other, despite, in some instances, being expressed in different cell types. TGFβ is a well-known stimulator of stromal cell replication and collagen synthesis and TGFβ treatment of cultured fetal ovarian stromal cells inhibited the expression of INSR, AR, C8H9orf3 and RAD50 and stimulated the expression of TGFB1I1. In human ovaries (n = 15, < 150 days gestation) many of the same genes as in bovine (FBN3, GATA4, HMGA2, FSHR, DENND1A and LHCGR but not TOX3 or FSHB) were expressed and correlated with each other. With so many relationships between PCOS candidate genes during development of the fetal ovary, including TGFβ and androgen signalling, we suggest that future studies should determine if perturbations of these genes in the fetal ovary can lead to PCOS in later life., Competing Interests: The authors have declared that no competing interests exist.

Cells on the surface of the mesonephros give rise to replicating Gonadal Ridge Epithelial-Like (GREL) cells, the first somatic cells of the gonadal ridge. Later germ cells associate with the GREL cells in the ovigerous cords, and the GREL cells subsequently give rise to the granulosa cells in follicles. To examine these events further, 27 bovine fetal ovaries of different gestational ages were collected and prepared for immunohistochemical localisation of collagen type I and Ki67 to identify regions of the ovary and cell proliferation, respectively. The non-stromal cortical areas (collagen-negative) containing GREL cells and germ cells and later in development, the follicles with oocytes and granulosa cells, were analysed morphometrically. Another set of ovaries (n = 17) were collected and the expression of genes associated with germ cell lineages and GREL/granulosa cells were quantitated by RT-PCR. The total volume of non-stromal areas in the cortex increased significantly and progressively with ovarian development, plateauing at the time the surface epithelium developed. However, the proportion of non-stromal areas in the cortex declined significantly and progressively throughout gestation, largely due to a cessation in growth of the non-stroma cells and the continued growth of stroma. The proliferation index in the non-stromal area was very high initially and then declined substantially at the time follicles formed. Thereafter, it remained low. The numerical density of the non-stromal cells was relatively constant throughout ovarian development. The expression levels of a number of genes across gestation either increased (AMH, FSHR, ESR1, INHBA), declined (CYP19A1, ESR2, ALDH1A1, DSG2, OCT4, LGR5) or showed no particular pattern (CCND2, CTNNB1, DAZL, FOXL2, GATA4, IGFBP3, KRT19, NR5A1, RARRES1, VASA, WNT2B). Many of the genes whose expression changed across gestation, were positively or negatively correlated with each other. The relationships between these genes may reflect their roles in the important events such as the transition of ovigerous cords to follicles, oogonia to oocytes or GREL cells to granulosa cells., Competing Interests: The authors have declared that no competing interests exist.

Movement and expansion of mesonephric-derived stroma appears to be very important in the development of the ovary. Here, we examined the expression of 24 genes associated with stroma in fetal ovaries during gestation (n = 17; days 58-274) from Bos taurus cattle. RNA was isolated from ovaries for quantitative RT-PCR. Expression of the majority of genes in TGFβ signalling, stromal transcription factors (NR2F2, AR), and some stromal matrix genes (COL1A1, COL3A1 and FBN1, but not FBN3) showed a positive linear increase with gestational age. Expression of genes associated with follicles (INSL3, CYP17A1, CYP11A1 and HSD3B1), was low until mid-gestation and then increased with gestational age. LHCGR showed an unusual bimodal pattern; high levels in the first and last trimesters. RARRES1 and IGFBP3 also increased with gestational age. To relate changes in gene expression in stromal cells with that in non stromal cells during development of the ovary we combined the data on the stromal genes with another 20 genes from non stromal cells published previously and then performed hierarchical clustering analysis. Three major clusters were identified. Cluster 1 genes (GATA4, FBN3, LHCGR, CYP19A1, ESR2, OCT4, DSG2, TGFB1, CCND2, LGR5, NR5A1) were characterised by high expression only in the first trimester. Cluster 2 genes (FSHR, INSL3, HSD3B1, CYP11A1, CYP17A1, AMH, IGFBP3, INHBA) were highly expressed in the third trimester and largely associated with follicle function. Cluster 3 (COL1A1, COL3A1, FBN1, TGFB2 TGFB3, TGFBR2, TGFBR3, LTBP2, LTBP3, LTBP4, TGFB1I1, ALDH1A1, AR, ESR1, NR2F2) had much low expression in the first trimester rising in the second trimester and remaining at that level during the third trimester. Cluster 3 contained members of two pathways, androgen and TGFβ signalling, including a common member of both pathways namely the androgen receptor cofactor TGFβ1 induced transcript 1 protein (TGFB1I1; hic5). GATA4, FBN3 and LHCGR, were highly correlated with each other and were expressed highly in the first trimester during stromal expansion before follicle formation, suggesting that this could be a critical phase in the development of the ovarian stroma., Competing Interests: The authors have declared that no competing interests exist.

In vitro culture of ovarian granulosa cells and theca cells has been very important for our understanding of their function and regulation. One of the most eagerly sought attributes of cell culture is the use of chemically-defined conditions. However, even under such in vitro conditions cell behaviour could differ from the in vivo situation because of differences in oxygen tension, nutrients, adhesion matrix and other factors. To examine this further we compared the transcriptomes of both granulosa cells and cells from the theca interna that were cultured in what are arguably the best in vitro conditions for maintaining the 'follicular' phenotypes of both tissue types, as displayed by their respective freshly-isolated counterparts. The array data analysed are from recently published data and use the same sizes of bovine follicles (small antral 3-6 mm) and the same Affymetrix arrays. We conducted analysis using Partek, Ingenuity Pathway Analysis and GOEAST. Principal Component Analysis (PCA) and hierarchical clustering clearly separated the in vivo from the in vitro groups for both cells types and transcriptomes were more homogeneous upon culture. In both cell cultures behaviours associated with cell adhesion, migration and interaction with matrix or substrate were more abundant. However, the pathways involved generally differed between the two cell types. With the thecal cultures a gene expression signature of an immune response was more abundant, probably by leukocytes amongst the cells cultured from the theca interna. These results indicate differences between in vivo and in vitro that should be considered when interpreting in vitro data.

In studies using isolated ovarian granulosa and thecal cells it is important to assess the degree of cross contamination. Marker genes commonly used for granulosa cells include FSHR, CYP19A1 and AMH while CYP17A1 and INSL3 are used for thecal cells. To increase the number of marker genes available we compared expression microarray data from isolated theca interna with that from granulosa cells of bovine small (n = 10 for both theca and granulosa cells; 3-5 mm) and large (n = 4 for both theca and granulosa cells, > 9 mm) antral follicles. Validation was conducted by qRT-PCR analyses. Known markers such as CYP19A1, FSHR and NR5A2 and another 11 genes (LOC404103, MGARP, GLDC, CHST8, CSN2, GPX3, SLC35G1, CA8, CLGN, FAM78A, SLC16A3) were common to the lists of the 50 most up regulated genes in granulosa cells from both follicle sizes. The expression in theca interna was more consistent than in granulosa cells between the two follicle sizes. Many genes up regulated in theca interna were common to both sizes of follicles (MGP, DCN, ASPN, ALDH1A1, COL1A2, FN1, COL3A1, OGN, APOD, COL5A2, IGF2, NID1, LHFP, ACTA2, DUSP12, ACTG2, SPARCL1, FILIP1L, EGFLAM, ADAMDEC1, HPGD, COL12A1, FBLN5, RAMP2, COL15A1, PLK2, COL6A3, LOXL1, RARRES1, FLI1, LAMA2). Many of these were stromal extracellular matrix genes. MGARP, GLDC, CHST8, GPX3 were identified as new potential markers for granulosa cells, while FBLN5, OGN, RAMP2 were significantly elevated in the theca interna.

The theca interna is a specialized stromal layer that envelops each growing ovarian follicle. It contains capillaries, fibroblasts, immune cells and the steroidogenic cells that synthesize androgens for conversion to estradiol by the neighboring granulosa cells. During reproductive life only a small number of follicles will grow to a sufficient size to ovulate, whereas the majority of follicles will undergo regression/atresia and phagocytosis by macrophages. To identify genes which are differentially regulated in the theca interna during follicular atresia, we undertook transcriptome profiling of the theca interna from healthy (n = 10) and antral atretic (n = 5) bovine follicles at early antral stages (<5 mm). Principal Component Analyses and hierarchical classification of the signal intensity plots for the arrays showed primary clustering into two groups, healthy and atretic. A total of 543 probe sets were differentially expressed between the atretic and healthy theca interna. Further analyses of these genes by Ingenuity Pathway Analysis and Gene Ontology Enrichment Analysis Toolkit software found most of the genes being expressed were related to cytokines, hormones and receptors as well as the cell cycle and DNA replication. Cell cycle genes which encode components of the replicating chromosome complex and mitotic spindle were down-regulated in atretic theca interna, whereas stress response and inflammation-related genes such as TP53, IKBKB and TGFB1 were up-regulated. In addition to cell cycle regulators, upstream regulators that were predicted to be inhibited included Retinoblastoma 1, E2 transcription factor 1, and hepatocyte growth factor. Our study suggests that during antral atresia of small follicles in the theca interna, arrest of cell cycle and DNA replication occurs rather than up- regulation of apoptosis-associated genes as occurs in granulosa cells.

The theca interna layer of the ovarian follicle forms during the antral stage of follicle development and lies adjacent to and directly outside the follicular basal lamina. It supplies androgens and communicates with the granulosa cells and the oocyte by extracellular signaling. To better understand developmental changes in the theca interna, we undertook transcriptome profiling of the theca interna from small (3-5 mm, n = 10) and large (9-12 mm, n = 5) healthy antral bovine follicles, representing a calculated >7-fold increase in the amount of thecal tissue. Principal Component Analysis and hierarchical classification of the signal intensity plots for the arrays showed no clustering of the theca interna samples into groups depending on follicle size or subcategories of small follicles. From the over 23,000 probe sets analysed, only 76 were differentially expressed between large and small healthy follicles. Some of the differentially expressed genes were associated with processes such as myoblast differentiation, protein ubiquitination, nitric oxide and transforming growth factor β signaling. The most significant pathway affected from our analyses was found to be Wnt signaling, which was suppressed in large follicles via down-regulation of WNT2B and up-regulation of the inhibitor FRZB. These changes in the transcriptional profile could have been due to changes in cellular function or alternatively since the theca interna is composed of a number of different cell types it could have been due to any systematic change in the volume density of any particular cell type. However, our study suggests that the transcriptional profile of the theca interna is relatively stable during antral follicle development unlike that of granulosa cells observed previously. Thus both the cellular composition and cellular behavior of the theca interna and its contribution to follicular development appear to be relatively constant throughout the follicle growth phase examined.

Ovarian follicular granulosa cells surround and nurture oocytes, and produce sex steroid hormones. It is believed that during development the ovarian surface epithelial cells penetrate into the ovary and develop into granulosa cells when associating with oogonia to form follicles. Using bovine fetal ovaries (n = 80) we identified a novel cell type, termed GREL for Gonadal Ridge Epithelial-Like. Using 26 markers for GREL and other cells and extracellular matrix we conducted immunohistochemistry and electron microscopy and chronologically tracked all somatic cell types during development. Before 70 days of gestation the gonadal ridge/ovarian primordium is formed by proliferation of GREL cells at the surface epithelium of the mesonephros. Primordial germ cells (PGCs) migrate into the ovarian primordium. After 70 days, stroma from the underlying mesonephros begins to penetrate the primordium, partitioning the developing ovary into irregularly-shaped ovigerous cords composed of GREL cells and PGCs/oogonia. Importantly we identified that the cords are always separated from the stroma by a basal lamina. Around 130 days of gestation the stroma expands laterally below the outermost layers of GREL cells forming a sub-epithelial basal lamina and establishing an epithelial-stromal interface. It is at this stage that a mature surface epithelium develops from the GREL cells on the surface of the ovary primordium. Expansion of the stroma continues to partition the ovigerous cords into smaller groups of cells eventually forming follicles containing an oogonium/oocyte surrounded by GREL cells, which become granulosa cells, all enclosed by a basal lamina. Thus in contrast to the prevailing theory, the ovarian surface epithelial cells do not penetrate into the ovary to form the granulosa cells of follicles, instead ovarian surface epithelial cells and granulosa cells have a common precursor, the GREL cell.