Transcriptome profiling of granulosa cells from bovine ovarian follicles during atresia

Background: The major function of the ovary is to produce oocytes for fertilisation. Oocytes mature in follicles surrounded by nurturing granulosa cells and all are enclosed by a basal lamina. During growth, granulosa cells replicate and a large fluid-filled cavity (the antrum) develops in the centre. Only follicles that have enlarged to over 10 mm can ovulate in cows. In mammals, the number of primordial follicles far exceeds the numbers that ever ovulate and atresia or regression of follicles is a mechanism to regulate the number of oocytes ovulated and to contribute to the timing of ovulation. To better understand the molecular basis of follicular atresia, we undertook transcriptome profiling of granulosa cells from healthy (n = 10) and atretic (n = 5) bovine follicles at early antral stages (< 5 mm). Results: Principal Component Analysis (PCA) and hierarchical classification of the signal intensity plots for the arrays showed primary clustering into two groups, healthy and atretic. These analyses and size-frequency plots of coefficients of variation of signal intensities revealed that the healthy follicles were more heterogeneous. Examining the differentially-expressed genes the most significantly affected functions in atretic follicles were cell death, organ development, tissue development and embryonic development. The overall processes influenced by transcription factor gene TP53 were predicted to be activated, whereas those of MYC were inhibited on the basis of known interactions with the genes in our dataset. The top ranked canonical pathway contained signalling molecules common to various inflammatory/fibrotic pathways such as the transforming growth factor-beta and tumour necrosis factor-alpha pathways. The two most significant networks also reflect this pattern of tissue remodelling/fibrosis gene expression. These networks also contain molecules which are present in the canonical pathways of hepatic fibrosis/hepatic stellate cell activation and transforming growth factor-beta signalling and were up regulated. Conclusions: Small healthy antral follicles, which have a number of growth outcomes, exhibit greater variability in gene expression, particularly in genes associated with cell division and other growth-related functions. Atresia, on the other hand, not only involves cell death but clearly is an active process similar to wound healing.Nicholas Hatzirodos, Katja Hummitzsch, Helen F Irving-Rodgers, Margaret L Harland, Stephanie E Morris and Raymond J Rodger

The global effect of follicle-stimulating hormone and tumour necrosis factor α on gene expression in cultured bovine ovarian granulosa cells

Background Oocytes mature in ovarian follicles surrounded by granulosa cells. During follicle growth, granulosa cells replicate and secrete hormones, particularly steroids close to ovulation. However, most follicles cease growing and undergo atresia or regression instead of ovulating. To investigate the effects of stimulatory (follicle-stimulating hormone; FSH) and inhibitory (tumour necrosis factor alpha; TNFα) factors on the granulosa cell transcriptome, bovine ovaries were obtained from a local abattoir and pools of granulosa cells were cultured in vitro for six days under defined serum-free conditions with treatments present on days 3–6. Initially dose–response experiments (n = 4) were performed to determine the optimal concentrations of FSH (0.33 ng/ml) and TNFα (10 ng/ml) to be used for the microarray experiments. For array experiments cells were cultured under control conditions, with FSH, with TNFα, or with FSH plus TNFα (n = 4 per group) and RNA was harvested for microarray analyses. Results Statistical analysis showed primary clustering of the arrays into two groups, control/FSH and TNFα/TNFα plus FSH. The effect of TNFα on gene expression dominated that of FSH, with substantially more genes differentially regulated, and the pathways and genes regulated by TNFα being similar to those of FSH plus TNFα treatment. TNFα treatment reduced the endocrine activity of granulosa cells with reductions in expression of FST, INHA, INBA and AMH. The top-ranked canonical pathways and GO biological terms for the TNFα treatments included antigen presentation, inflammatory response and other pathways indicative of innate immune function and fibrosis. The two most significant networks also reflect this, containing molecules which are present in the canonical pathways of hepatic fibrosis/hepatic stellate cell activation and transforming growth factor β signalling, and these were up regulated. Upstream regulator analyses also predicted TNF, interferons γ and β1 and interleukin 1β. Conclusions In vitro, the transcriptome of granulosa cells responded minimally to FSH compared with the response to TNFα. The response to TNFα indicated an active process akin to tissue remodelling as would occur upon atresia. Additionally there was reduction in endocrine function and induction of an inflammatory response to TNFα that displays features similar to immune cells

Transcriptome profiling of granulosa cells of bovine ovarian follicles during growth from small to large antral sizes

Background: At later stages of folliculogenesis, the mammalian ovarian follicle contains layers of epithelial granulosa cells surrounding an antral cavity. During follicle development granulosa cells replicate, secrete hormones and support the growth of the oocyte. In cattle, the follicle needs to grow > 10 mm in diameter to allow an oocyte to ovulate, following which the granulosa cells cease dividing and differentiate into the specialised cells of the corpus luteum. To better understand the molecular basis of follicular growth and granulosa cell maturation, we undertook transcriptome profiling of granulosa cells from small ( 10 mm, n = 4) healthy bovine follicles using Affymetrix microarrays (24,128 probe sets).Results: Principal component analysis for the first two components and hierarchical clustering showed clustering into two groups, small and large, with the former being more heterogeneous. Size-frequency distributions of the coefficient of variation of the signal intensities of each probe set also revealed that small follicles were more heterogeneous than the large. IPA and GO enrichment analyses revealed that processes of axonal guidance, immune signalling and cell rearrangement were most affected in large follicles. The most important networks were associated with: (A) Notch, SLIT/ROBO and PI3K signalling, and (B) ITGB5 and extracellular matrix signalling through extracellular signal related kinases (ERKs). Upstream regulator genes which were predicted to be active in large follicles included STAT and XBP1. By comparison, developmental processes such as those stimulated by KIT, IHH and MEST were most active in small follicles. MGEA5 was identified as an upstream regulator in small follicles. It encodes an enzyme that modifies the activity of many target proteins, including those involved in energy sensing, by removal of N-acetylglucosamine from serine and threonine residues.Conclusions: Our data suggest that as follicles enlarge more genes and/or pathways are activated than are inactivated, and gene expression becomes more uniform. These findings could be interpreted that either the cells in large follicles are more uniform in their gene expression, or that follicles are more uniform or a combination of both and that additional factors, such as LH, are additionally controlling the granulosa cells. © 2014 Hatzirodos et al.; licensee BioMed Central Ltd.Nicholas Hatzirodos, Helen F Irving-Rodgers, Katja Hummitzsch, Margaret L Harland, Stephanie E Morris and Raymond J Rodger

Transcriptome profiling of the theca interna from bovine ovarian follicles during atresia

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.Nicholas Hatzirodos, Helen F. Irving-Rodgers, Katja Hummitzsch, Raymond J. Rodger

Transcriptome comparisons identify new cell markers for theca interna and granulosa cells from small and large antral ovarian follicles

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.Nicholas Hatzirodos, Katja Hummitzsch, Helen F. Irving-Rodgers, Raymond J. Rodger

Exercise and insulin resistance in PCOS: muscle insulin signalling and fibrosis

OBJECTIVE:Mechanisms of insulin resistance in polycystic ovary syndrome (PCOS) remain ill-defined, contributing to sub-optimal therapies. Recognising skeletal muscle plays a key role in glucose homeostasis we investigated early insulin signalling, its association with aberrant transforming growth factor β (TGFβ) regulated tissue fibrosis. We also explored the impact of aerobic exercise on these molecular pathways. METHODS:A secondary analysis from a cross-sectional study was undertaken in women with (n=30) or without (n=29) PCOS across lean and overweight BMIs. A subset of participants with (n=8) or without (n=8) PCOS who were overweight completed 12-weeks of aerobic exercise training. Muscle was sampled before and 30 min into a euglycaemic-hyperinsulinaemic clamp pre- and post-training. RESULTS:We found reduced signalling in PCOS of mechanistic target of rapamycin (mTOR). Exercise training augmented but did not completely rescue this signalling defect in women with PCOS. Genes in the TGFβ signalling network were upregulated in skeletal muscle in the overweight women with PCOS but were unresponsive to exercise training except for genes encoding LOX, collagen 1 and 3. CONCLUSIONS:We provide new insights into defects in early insulin signalling, tissue fibrosis, and hyperandrogenism in PCOS-specific insulin resistance in lean and overweight women. PCOS-specific insulin-signalling defects were isolated to mTOR, while gene expression implicated TGFβ ligand regulating a fibrosis in the PCOS-obesity synergy in insulin resistance and altered responses to exercise. Interestingly, there was little evidence for hyperandrogenism as a mechanism for insulin resistance

Could perturbed fetal development of the ovary contribute to the development of polycystic ovary syndrome in later life?

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

Regulation of fibrillins and modulators of TGFβ in fetal bovine and human ovaries

Fibrillins 1–3 are stromal extracellular matrix proteins that play important roles in regulating TGFβ activity, which stimulates fibroblasts to proliferate and synthesize collagen. In the developing ovary, the action of stroma is initially necessary for the formation of ovigerous cords and subsequently for the formation of follicles and the surface epithelium of the ovary. FBN3 is highly expressed only in early ovarian development and then it declines. In contrast, FBN1 and 2 are upregulated in later ovarian development. We examined the expression of FBN1–3 in bovine and human fetal ovaries. We used cell dispersion and monolayer culture, cell passaging and tissue culture. Cells were treated with growth factors, hormones or inhibitors to assess the regulation of expression of FBN1–3. When bovine fetal ovarian tissue was cultured, FBN3 expression declined significantly. Treatment with TGFβ-1 increased FBN1 and FBN2 expression in bovine fibroblasts, but did not affect FBN3 expression. Additionally, in cultures of human fetal ovarian fibroblasts (9–17 weeks gestational age), the expression of FBN1 and FBN2 increased with passage, whereas FBN3 dramatically decreased. Treatment with activin A and a TGFβ family signaling inhibitor, SB431542, differentially regulated the expression of a range of modulators of TGFβ signaling and of other growth factors in cultured human fetal ovarian fibroblasts suggesting that TGFβ signaling is differentially involved in the regulation of ovarian fibroblasts. Additionally, since the changes in FBN1–3 expression that occur in vitro are those that occur with increasing gestational age in vivo, we suggest that the fetal ovarian fibroblasts mature in vitro.Nicole A Bastian, Rosemary A Bayne, Katja Hummitzsch, Nicholas Hatzirodos, Wendy M Bonner, Monica D Hartanti, Helen F Irving-Rodgers, Richard A Anderson and Raymond J Rodger

Morphometric analyses and gene expression related to germ cells, gonadal ridge epithelial-like cells and granulosa cells during development of the bovine fetal ovary

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.Katja Hummitzsch, Nicholas Hatzirodos, Helen F. Irving-Rodgers, Monica D. Hartanti, Viv E. A. Perry, Richard A. Anderson, Raymond J. Rodger

Transcript abundance of stromal and thecal cell related genes during bovine ovarian development

<div><p>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 <i>Bos taurus</i> cattle. RNA was isolated from ovaries for quantitative RT-PCR. Expression of the majority of genes in TGFβ signalling, stromal transcription factors (<i>NR2F2</i>, <i>AR)</i>, and some stromal matrix genes (<i>COL1A1</i>, <i>COL3A1</i> and <i>FBN1</i>, but not <i>FBN3</i>) showed a positive linear increase with gestational age. Expression of genes associated with follicles (<i>INSL3</i>, <i>CYP17A1</i>, <i>CYP11A1</i> and <i>HSD3B1</i>), was low until mid-gestation and then increased with gestational age. <i>LHCGR</i> showed an unusual bimodal pattern; high levels in the first and last trimesters. <i>RARRES1</i> and <i>IGFBP3</i> 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 (<i>GATA4</i>, <i>FBN3</i>, <i>LHCGR</i>, <i>CYP19A1</i>, <i>ESR2</i>, <i>OCT4</i>, <i>DSG2</i>, <i>TGFB1</i>, <i>CCND2</i>, <i>LGR5</i>, <i>NR5A1</i>) were characterised by high expression only in the first trimester. Cluster 2 genes (<i>FSHR</i>, <i>INSL3</i>, <i>HSD3B1</i>, <i>CYP11A1</i>, <i>CYP17A1</i>, <i>AMH</i>, <i>IGFBP3</i>, <i>INHBA</i>) were highly expressed in the third trimester and largely associated with follicle function. Cluster 3 (<i>COL1A1</i>, <i>COL3A1</i>, <i>FBN1</i>, <i>TGFB2 TGFB3</i>, <i>TGFBR2</i>, <i>TGFBR3</i>, <i>LTBP2</i>, <i>LTBP3</i>, <i>LTBP4</i>, <i>TGFB1I1</i>, <i>ALDH1A1</i>, <i>AR</i>, <i>ESR1</i>, <i>NR2F2</i>) 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 (<i>TGFB1I1</i>; hic5). <i>GATA4</i>, <i>FBN3</i> and <i>LHCGR</i>, 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.</p></div