Morphometric and gene expression analyses of stromal expansion during development of the bovine fetal ovary

Analysis of genetic diversity and selection signals in Chaling cattle of southern China using whole-genome scan

China has diversified resources of indigenous cattle, which are classified into Northern, Central, and Southern groups according to their geographical distribution. Chaling cattle belong to Southern group. This breed is famous for the production of good quality meat with elite meat grades. To analyze the genetic diversity of Chaling cattle, 20 samples were sequenced using whole-genome resequencing technology, along with 138 published whole-genome sequencing data of Indian indicine cattle, Chinese indicine cattle, East Asian taurine cattle, Eurasian taurine cattle, and European taurine cattle as control. It was found that Chaling cattle originated from Chinese indicine cattle. The genetic diversity of Chaling cattle is higher than that of Indian indicine cattle, East Asian taurine cattle, Eurasian taurine cattle, and European taurine cattle, but lower than that of Chinese indicine cattle and Xiangxi cattle. Annotating the selection signals obtained by composite likelihood ratio, θπ, F_ST , π-ratio, and XP-EHH methods, several genes associated with immunity, heat tolerance, reproduction, growth, and meat quality showed strong selection signals. In general, this study provides a theoretical basis for analyzing the genetic mechanism of Chaling cattle with excellent adaptability, rough feeding tolerance, good immune performance, and good meat quality. This work lays a foundation for genetic breeding of Chaling cattle in future.

Expression of PCOS candidate genes in bovine fetal and adult ovarian somatic cells

Polycystic ovary syndrome (PCOS) is an endocrine metabolic disorder that appears to have a genetic predisposition and a fetal origin. The fetal ovary has two major somatic cell types shown previously to be of different cellular origins, different morphologies and to differentially express 15 genes. We isolated the somatic gonadal ridge epithelial-like (GREL) cells (n = 7) and ovarian fetal fibroblasts (n = 6) by clonal expansion. Using qRT-PCR, we compared the gene expression levels of PCOS candidate genes with previous data on the expression levels in whole fetal ovaries across gestation. We also compared these levels with those in bovine adult ovarian cells including fibroblasts (n = 4), granulosa cells (n = 5) and surface epithelial cells (n = 5). Adult cell types exhibited clear differences in the expression of most genes. In fetal ovarian cells, DENND1A and ERBB3 had significantly higher expression in GREL cells. HMGA2 and TGFB1I1 tended to have higher expression in fetal fibroblasts than GREL cells. Another 19 genes did not exhibit differences between GREL cells and fetal fibroblasts and FBN3, FSHB, LHCGR, FSHR and ZBTB16 were very lowly expressed in GREL cells and fibroblasts. The culture of fetal fibroblasts in EGF-containing medium resulted in lower expression of NEIL2, but higher expression of MAPRE1 compared to culture in the absence of EGF. Thus, the two fetal ovarian somatic cell types mostly lacked differential expression of PCOS candidate genes.

A single-cell atlas of the cycling murine ovary

The estrous cycle is regulated by rhythmic endocrine interactions of the nervous and reproductive systems, which coordinate the hormonal and ovulatory functions of the ovary. Folliculogenesis and follicle progression require the orchestrated response of a variety of cell types to allow the maturation of the follicle and its sequela, ovulation, corpus luteum formation, and ovulatory wound repair. Little is known about the cell state dynamics of the ovary during the estrous cycle and the paracrine factors that help coordinate this process. Herein, we used single-cell RNA sequencing to evaluate the transcriptome of >34,000 cells of the adult mouse ovary and describe the transcriptional changes that occur across the normal estrous cycle and other reproductive states to build a comprehensive dynamic atlas of murine ovarian cell types and states.

Key signalling pathways underlying the aetiology of polycystic ovary syndrome

Polycystic ovary syndrome (PCOS) is a common endocrine condition characterised by a range of reproductive, endocrine, metabolic and psychological abnormalities. Reports estimate that around 10% of women of reproductive age are affected by PCOS, representing a significant prevalence worldwide, which poses a high economic health burden. As the origin of PCOS remains largely unknown, there is neither a cure nor mechanism-based treatments leaving patient management suboptimal and focused solely on symptomatic treatment. However, if the underlying mechanisms underpinning the development of PCOS were uncovered then this would pave the way for the development of new interventions for PCOS. Recently, there have been significant advances in our understanding of the underlying pathways likely involved in PCOS pathogenesis. Key insights include the potential involvement of androgens, insulin, anti-Müllerian hormone and transforming growth factor beta in the development of PCOS. This review will summarise the significant scientific discoveries on these factors that have enhanced our knowledge of the mechanisms involved in the development of PCOS and discuss the impact these insights may have in shaping the future development of effective strategies for women with PCOS.

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

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.

Ovarian Biomechanics: From Health to Disease

Biomechanics is a physical phenomenon which mainly related with deformation and movement of life forms. As a mechanical signal, it participates in the growth and development of many tissues and organs, including ovary. Mechanical signals not only participate in multiple processes in the ovary but also play a critical role in ovarian growth and normal physiological functions. Additionally, the involvement of mechanical signals has been found in ovarian cancer and other ovarian diseases, prompting us to focus on the roles of mechanical signals in the process of ovarian health to disease. This review mainly discusses the effects and signal transduction of biomechanics (including elastic force, shear force, compressive stress and tensile stress) in ovarian development as a regulatory signal, as well as in the pathological process of normal ovarian diseases and cancer. This review also aims to provide new research ideas for the further research and treatment of ovarian-related diseases.

OUP accepted manuscript

STUDY QUESTION

Could changes in transforming growth factor β (TGFβ) signalling during foetal ovary development alter the expression of polycystic ovary syndrome (PCOS) candidate genes leading to a predisposition to PCOS?

SUMMARY ANSWER

TGFβ signalling molecules are dynamically expressed during foetal ovary development and TGFβ1 inhibits expression of the androgen receptor (AR) and 7 (INSR, C8H9orf3, RAD50, ERBB3, NEIL2, IRF1 and ZBTB16) of the 25 PCOS candidate genes in foetal ovarian fibroblasts in vitro, whilst increasing expression of the AR cofactor TGFβ-induced transcript 1 (TGFB1I1 or Hic5).

WHAT IS KNOWN ALREADY

The ovarian stroma arises from the mesonephros during foetal ovary development. Changes in the morphology of the ovarian stroma are cardinal features of PCOS. The ovary is more fibrous and has more tunica and cortical and subcortical stroma. It is not known why this is and when this arises. PCOS has a foetal origin and perhaps ovarian stroma development is altered during foetal life to determine the formation of a polycystic ovary later in life. PCOS also has a genetic origin with 19 loci containing 25 PCOS candidate genes. In many adult tissues, TGFβ is known to stimulate fibroblast replication and collagen deposition in stroma, though it has the opposite effect in the non-scaring foetal tissues. Our previous studies showed that TGFβ signalling molecules [TGFβs and their receptors, latent TGFβ binding proteins (LTBPs) and fibrillins, which are extracellular matrix proteins that bind LTBPs] are expressed in foetal ovaries. Also, we previously showed that TGFβ1 inhibited expression of AR and 3 PCOS candidate genes (INSR, C8H9orf3 and RAD50) and stimulated expression of TGFB1I1 in cultured foetal ovarian fibroblasts.

STUDY DESIGN, SIZE, DURATION

We used Bos taurus for this study as we can ethically collect foetal ovaries from across the full 9-month gestational period. Foetal ovaries (62–276 days, n = 19) from across gestation were collected from pregnant B. taurus cows for RNA-sequencing (RNA-seq) analyses. Foetal ovaries from B. taurus cows were collected (160–198 days, n = 6) for culture of ovarian fibroblasts.

PARTICIPANTS/MATERIALS, SETTING, METHODS

RNA-seq transcriptome profiling was performed on foetal ovaries and the data on genes involved in TGFβ signalling were extracted. Cells were dispersed from foetal ovaries and fibroblasts cultured and treated with TGFβ1. The effects of TGFβ regulation on the remaining eight PCOS candidate genes not previously studied (ERBB3, MAPRE1, FDFT1, NEIL2, ARL14EP, PLGRKT, IRF1 and ZBTB16) were examined.

MAIN RESULTS AND THE ROLE OF CHANCE

Many TGFβ signalling molecules are expressed in the foetal ovary, and for most, their expression levels increased accross gestation (LTBP1/2/3/4, FBN1, TGFB2/3, TGFBR2/3 and TGFB1I1), while a few decreased (FBN3, TGFBR3L, TGFBI and TGFB1) and others remained relatively constant (TGFBRAP1, TGFBR1 and FBN2). TGFβ1 significantly decreased expression of PCOS candidate genes ERBB3, NEIL2, IRF1 and ZBTB16 in cultured foetal ovarian fibroblasts.

LARGE SCALE DATA

The FASTQ files, normalized data and experimental information have been deposited in the Gene Expression Omnibus (GEO) accessible by accession number GSE178450.

LIMITATIONS, REASONS FOR CAUTION

Regulation of PCOS candidate genes by TGFβ was carried out in vitro and further studies in vivo are required. This study was carried out in bovine where foetal ovaries from across all of the 9-month gestational period were available, unlike in the human where it is not ethically possible to obtain ovaries from the second half of gestation.

WIDER IMPLICATIONS OF THE FINDINGS

From our current and previous results we speculate that inhibition of TGFβ signalling in the foetal ovary is likely to (i) increase androgen sensitivity by enhancing expression of AR, (ii) increase stromal activity by stimulating expression of COL1A1 and COL3A1 and (iii) increase the expression of 7 of the 25 PCOS candidate genes. Thus inhibition of TGFβ signalling could be part of the aetiology of PCOS or at least the aetiology of polycystic ovaries.

STUDY FUNDING/COMPETING INTEREST(S)

Funding was received from Adelaide University China Fee Scholarship (M.L.), Australian Research Training Program (R.A.) and the Faculty of Health and Medical Science Divisional Scholarship (R.A.), Adelaide Graduate Research Scholarships (R.A. and N.A.B.), Australia Awards Scholarship (M.D.H.), Robinson Research Institute Career Development Fellowship (K.H.) and Building On Ideas Grant (K.H.), National Health and Medical Research Council of Australia Centre for Research Excellence in the Evaluation, Management and Health Care Needs of Polycystic Ovary Syndrome (N.A.B., M.D.H. and R.J.R.; GTN1078444) and the Centre for Research Excellence on Women’s Health in Reproductive life (R.A., R.J.R. and K.H.; GTN1171592) and the UK Medical Research Council (R.A.A.; grant no. G1100357). The funders did not play any role in the study design, data collection and analysis, decision to publish or preparation of the manuscript. The 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.

Role of Fibulins in Embryonic Stage Development and Their Involvement in Various Diseases

The extracellular matrix (ECM) plays an important role in the evolution of early metazoans, as it provides structural and biochemical support to the surrounding cells through the cell–cell and cell–matrix interactions. In multi-cellular organisms, ECM plays a pivotal role in the differentiation of tissues and in the development of organs. Fibulins are ECM glycoproteins, found in a variety of tissues associated with basement membranes, elastic fibers, proteoglycan aggregates, and fibronectin microfibrils. The expression profile of fibulins reveals their role in various developmental processes such as elastogenesis, development of organs during the embryonic stage, tissue remodeling, maintenance of the structural integrity of basement membrane, and elastic fibers, as well as other cellular processes. Apart from this, fibulins are also involved in the progression of human diseases such as cancer, cardiac diseases, congenital disorders, and chronic fibrotic disorders. Different isoforms of fibulins show a dual role of tumor-suppressive and tumor-promoting activities, depending on the cell type and cellular microenvironment in the body. Knockout animal models have provided deep insight into their role in development and diseases. The present review covers details of the structural and expression patterns, along with the role of fibulins in embryonic development and disease progression, with more emphasis on their involvement in the modulation of cancer diseases.

Analysis of expression of candidate genes for polycystic ovary syndrome in adult and fetal human and fetal bovine ovaries†

Polycystic ovary syndrome (PCOS) appears to have a genetic predisposition and a fetal origin. We compared the expression levels of 25 PCOS candidate genes from adult control and PCOS human ovaries (n = 16) using microarrays. Only one gene was potentially statistically different. Using qRT-PCR, expression of PCOS candidate genes was examined in bovine fetal ovaries from early stages when they first developed stroma through to completion of development (n = 27; 60-270 days of gestation). The levels of ERBB3 mRNA negatively correlated with gestational age but positively with HMGA2, FBN3, TOX3, GATA4, and DENND1A.X1,2,3,4, previously identified as correlated with each other and expressed early. PLGRKT and ZBTB16, and less so IRF1, were also correlated with AMH, FSHR, AR, INSR, and TGFB1I1, previously identified as correlated with each other and expressed late. ARL14EP, FDFT1, NEIL2, and MAPRE1 were expressed across gestation and not correlated with gestational age as shown previously for THADA, ERBB4, RAD50, C8H9orf3, YAP1, RAB5B, SUOX, and KRR1. LHCGR, because of its unusual bimodal expression pattern, had some unusual correlations with other genes. In human ovaries (n = 15; <150 days of gestation), ERBB3.V1 and ERBB3.VS were expressed and correlated negatively with gestational age and positively with FBN3, HMGA2, DENND1A.V1,3,4, DENND1A.V1-7, GATA4, and FSHR, previously identified as correlated with each other and expressed early. Thus, the general lack of differential expression of candidate genes in adult ovaries contrasting with dynamic patterns of gene expression in fetal ovaries is consistent with a vulnerability to disturbance in the fetal ovary that may underpin development of PCOS.

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.

Formation of the Bovine Ovarian Surface Epithelium during Fetal Development

When first formed, the ovary only has an established epithelium at its base or hilum. Later, an epithelium is established around the rest of the ovary. To examine this further, we conducted scanning electron microscopy of the surface of bovine fetal ovaries and immunohistochemistry of ovarian cross-sections. From the earliest time point, the cells on the surface of the base or hilum of the ovary were cuboidal. On the remainder of the ovary, the surface was more irregular. By mid-development, the surface was covered completely with either a stratified or simple epithelium of cuboidal cells. Clefts were observed in the surface and appeared to form due to the expansion of stroma surrounding each open ovigerous cord, elevating the areas surrounding each cord, while leaving the opening of the cord to form the base of each cleft. The continued expansion of the surrounding stroma below the surface appeared not only to close the ovigerous cords from the surface but to compress the clefts into the shape of a groove. Later, most of the ovarian surface was covered with a simple cuboidal epithelium. The changes to the ovarian surface during fetal development coincide with the remodeling of the stroma and cords below.

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.

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

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.