Human granulosa cells function as innate immune cells executing an inflammatory reaction during ovulation: a microarray analysis

Progesterone and 17-hydroxy-progesterone concentrations in follicular fluid and serum reflect their production in granulosa and theca cells

RESEARCH QUESTION

How is the production of progesterone (P4) and 17-hydroxy-P4 (17-OH-P4) regulated between theca cells and granulosa cells during the follicular phase, during ovulation and after transformation into a corpus luteum?

DESIGN

Three cohorts were examined: (i) 31 women undergoing natural and stimulated cycles, with serum hormone measurements taken every 3 days; (ii) 50 women undergoing ovarian stimulation, with hormone concentrations in serum and follicular fluid assessed at five time points during final follicle maturation; and (iii) 12 women undergoing fertility preservation, with hormone concentrations evaluated via the follicular fluid of small antral follicles.

RESULTS

In the early follicular phase, theca cells primarily synthesized 17-OH-P4 while granulosa cells produced limited P4, maintaining the P4:17-OH-P4 ratio <1. As follicles reached follicle selection at a diameter of approximately 10 mm, P4 synthesis in granulosa cells was up-regulated, but P4 was mainly accumulated in follicular fluid. During final maturation, enhanced activity of the enzyme HSD3B2 in granulosa cells enhanced P4 production, with the P4:17-OH-P4 ratio increasing to >1. The concentration of 17-OH-P4 in the luteal phase was similar to that in the follicular phase, but P4 production increased in the luteal phase, yielding a P4:17-OH-P4 ratio significantly >1.

CONCLUSIONS

The P4:17-OH-P4 ratio reflects the activity of granulosa cells and theca cells during the follicular phase and following luteinization in the corpus luteum. Managing the function of granulosa cells is key for reducing the concentration of P4 during ovarian stimulation, but the concerted action of FSH and LH on granulosa cells during the second half of the follicular phase makes this complex.

Leptin signalling regulates transcriptional differences in granulosa cells from genetically obese mice but not the activation of NLRP3 inflammasome

Obesity is associated with increased ovarian inflammation and the establishment of leptin resistance. We presently investigated the role of impaired leptin signalling on transcriptional regulation in granulosa cells (GCs) collected from genetically obese mice. Furthermore, we characterised the association between ovarian leptin signalling, the activation of the NOD-like receptor protein 3 (NLRP3) inflammasome and macrophage infiltration in obese mice. After phenotype characterisation, ovaries were collected from distinct group of animals for protein and mRNA expression analysis: (i) mice subjected to a diet-induced obesity (DIO) protocol, where one group was fed a high-fat diet (HFD) and another a standard chow diet (CD) for durations of 4 or 16 weeks; (ii) mice genetically deficient in the long isoform of the leptin receptor (ObRb; db/db); (iii) mice genetically deficient in leptin (ob/ob); and (iv) mice rendered pharmacologically hyperleptinemic (LEPT). Next, GCs from antral follicles isolated from db/db and ob/ob mice were subjected to transcriptome analysis. Transcriptional analysis revealed opposing profiles in genes associated with steroidogenesis and prostaglandin action between the genetic models, despite the similarities in body weight. Furthermore, we observed no changes in the mRNA and protein levels of NLRP3 inflammasome components in the ovaries of db/db mice or in markers of M1 and M2 macrophage infiltration. This contrasted with the downregulation of NLRP3 inflammasome components and M1 markers in ob/ob and 16-wk HFD-fed mice. We concluded that leptin signalling regulates NLRP3 inflammasome activation and the expression of M1 markers in the ovaries of obese mice in an ObRb-dependent and ObRb-independent manner. Furthermore, we found no changes in the expression of leptin signalling and NLRP3 inflammasome genes in GCs from db/db and ob/ob mice, which was associated with no effects on macrophage infiltration genes, despite the dysregulation of genes associated with steroidogenesis in homozygous obese db/db. Our results suggest that: (i) the crosstalk between leptin signalling, NLRP3 inflammasome and macrophage infiltration takes place in ovarian components other than the GC compartment; and (ii) transcriptional changes in GCs from homozygous obese ob/ob mice suggest structural rearrangement and organisation, whereas in db/db mice the impairment in steroidogenesis and secretory activity.

The intrafollicular concentrations of biologically active cortisol in women rise abruptly shortly before ovulation and follicular rupture

STUDY QUESTION

What is the temporal activity and the concentration in follicular fluid (FF) of the anti-inflammatory steroid cortisol during the ovulatory process in humans?

SUMMARY ANSWER

Intrafollicular concentrations of cortisol become massively upregulated close to ovulation concomitant with an exceptionally high biological activity securing a timely and efficient termination of inflammatory processes.

WHAT IS KNOWN ALREADY

Ovulation has been described as a local, controlled inflammatory process resulting in the degeneration of the follicle wall which facilitate oocyte extrusion. Ovulation also affects the glucocorticoid metabolism of granulosa cells (GCs) and although de novo synthesis of cortisol only occurs in the adrenal cortex, the mid-cycle surge has been shown to induce a change from high expression of HSD11B2, inactivating cortisol to cortisone, to high expression of HSD11B1 which reversibly catalyses cortisol production from cortisone. Furthermore, high concentrations of progesterone and 17OH-progesterone within follicles may cause dislodging of cortisol from cortisol binding protein (CBP) thereby activating the biological activity of cortisol.

STUDY DESIGN, SIZE, DURATION

This prospective cohort study included 50 women undergoing fertility treatment according to a standard antagonist protocol at a university hospital-affiliated fertility clinic in Denmark.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Women donated FF and GCs from one follicle for research purpose aspirated at one of four time points during the process of final maturation of follicles: T = 0 h, T = 12 h, T = 17 h, T = 32 h. A second sample was collected at oocyte pick up at T = 36 h. The concentration of cortisol and cortisone together with a range of sex steroids was measured by LC-MS/MS in FF collected at the five time points mentioned above. Whole genome microarray data, validated by q-PCR analysis, was used to evaluate gene expression of CYP11B1, CYP21A2, HSD11B1, HSD11B2, and NR3C1 in GCs at the same time points.

MAIN RESULTS AND THE ROLE OF CHANCE

The concentration of cortisol was significantly increased from a few nM at 0 h to around 100-140 nM (P ≤ 0.0001) at 32-36 h, whilst cortisone was almost constant from 0 to 17 h at a concentration of between 90 and 100 nM being significantly reduced to 25-40 nM (P ≤ 0.0001) at 32-36 h. This was paralleled by a 690-fold upregulation of HSD11B1 from 0 to 12 h increasing to a more than 20.000-fold change at 36 h. HSD11B2 was quickly downregulated 15- to 20-fold after ovulation induction. Concentrations of progesterone and 17OH-progesterone increased during the ovulatory process to high levels which in essence displaces cortisol from its binding protein CBP due to similar binding affinities. Furthermore, a significant decrease in 11-deoxycortisol expression was seen, but CYP11B1 expression was below detection limit in GCs.

LIMITATIONS, REASONS FOR CAUTION

The study included women undergoing ovarian stimulation and results may differ from the natural cycle. More observations at each specific time point may have strengthened the conclusions. Furthermore, we have not been able to measure the actual active biological concentration of cortisol.

WIDER IMPLICATIONS OF THE FINDINGS

For the first time, this study collectively evaluated the temporal pattern of cortisol and cortisone concentrations during human ovulation, rendering a physiological framework for understanding potential dysregulations in the inflammatory reaction of ovulation.

STUDY FUNDING/COMPETING INTEREST(S)

This research was supported by the University Hospital of Copenhagen, Rigshospitalet, and Novo Nordisk Foundation grant number NNF21OC00700556. Interreg V ÔKS through ReproUnion (www.reprounion.eu); Region Zealand Research Foundation. The funders had no role in study design, collection of data, analyses, writing of the article, or the decision to submit it for publication. The authors have no conflicts of interest to declare.

TRIAL REGISTRATION NUMBER

N/A.

Germline cell de novo mutations and potential effects of inflammation on germline cell genome stability

Uncontrolled pathogenic genome mutations in germline cells might impair adult fertility, lead to birth defects or even affect the adaptability of a species. Understanding the sources of DNA damage, as well as the features of damage response in germline cells are the overarching tasks to reduce the mutations in germline cells. With the accumulation of human genome data and genetic reports, genome variants formed in germline cells are being extensively explored. However, the sources of DNA damage, the damage repair mechanisms, and the effects of DNA damage or mutations on the development of germline cells are still unclear. Besides exogenous triggers of DNA damage such as irradiation and genotoxic chemicals, endogenous exposure to inflammation may also contribute to the genome instability of germline cells. In this review, we summarized the features of de novo mutations and the specific DNA damage responses in germline cells and explored the possible roles of inflammation on the genome stability of germline cells.

Toxicological Mechanisms and Potencies of Organophosphate Esters in KGN Human Ovarian Granulosa Cells as Revealed by High-throughput Transcriptomics

Despite the growing number of studies reporting potential risks associated with exposure to organophosphate esters (OPEs), their molecular mechanisms of action remain poorly defined. We used the high-throughput TempO-Seq™ platform to investigate the effects of frequently detected OPEs on the expression of ∼3000 environmentally responsive genes in KGN human ovarian granulosa cells. Cells were exposed for 48 h to one of five OPEs (0.1 to 50 μM): tris(methylphenyl) phosphate (TMPP), isopropylated triphenyl phosphate (IPPP), tert-butylphenyl diphenyl phosphate (BPDP), triphenyl phosphate (TPHP), or tris(2-butoxyethyl) phosphate (TBOEP). The sequencing data indicate that four OPEs induced transcriptional changes, whereas TBOEP had no effect within the concentration range tested. Multiple pathway databases were used to predict alterations in biological processes based on differentially expressed genes. At lower concentrations, inhibition of the cholesterol biosynthetic pathway was the predominant effect of OPEs; this was likely a consequence of intracellular cholesterol accumulation. At higher concentrations, BPDP and TPHP had distinct effects, primarily affecting pathways involved in cell cycle progression and other stress responses. Benchmark concentration (BMC) modelling revealed that BPDP had the lowest transcriptomic point of departure. However, in vitro to in vivo extrapolation modeling indicated that TMPP was bioactive at lower concentrations than the other OPEs. We conclude that these new approach methodologies provide information on the mechanism(s) underlying the effects of data-poor compounds and assist in the derivation of protective points of departure for use in chemical read-across and decision-making.

Comprehensive transcriptomic analysis of long non-coding RNAs in bovine ovarian follicles and early embryos

Long non-coding RNAs (lncRNAs) have been the subject of numerous studies over the past decade. First thought to come from aberrant transcriptional events, lncRNAs are now considered a crucial component of the genome with roles in multiple cellular functions. However, the functional annotation and characterization of bovine lncRNAs during early development remain limited. In this comprehensive analysis, we review lncRNAs expression in bovine ovarian follicles and early embryos, based on a unique database comprising 468 microarray hybridizations from a single platform designed to target 7,724 lncRNA transcripts, of which 5,272 are intergenic (lincRNA), 958 are intronic, and 1,524 are antisense (lncNAT). Compared to translated mRNA, lncRNAs have been shown to be more tissue-specific and expressed in low copy numbers. This analysis revealed that protein-coding genes and lncRNAs are both expressed more in oocytes. Differences between the oocyte and the 2-cell embryo are also more apparent in terms of lncRNAs than mRNAs. Co-expression network analysis using WGCNA generated 25 modules with differing proportions of lncRNAs. The modules exhibiting a higher proportion of lncRNAs were found to be associated with fewer annotated mRNAs and housekeeping functions. Functional annotation of co-expressed mRNAs allowed attribution of lncRNAs to a wide array of key cellular events such as meiosis, translation initiation, immune response, and mitochondrial related functions. We thus provide evidence that lncRNAs play diverse physiological roles that are tissue-specific and associated with key cellular functions alongside mRNAs in bovine ovarian follicles and early embryos. This contributes to add lncRNAs as active molecules in the complex regulatory networks driving folliculogenesis, oogenesis and early embryogenesis all of which are necessary for reproductive success.

A single-cell gene expression atlas of human follicular aspirates: Identification of leukocyte subpopulations and their paracrine factors

Leukocytes are in situ regulators critical for ovarian function. However, little is known about leukocyte subpopulations and their interaction with follicular cells in ovulatory follicles, especially in humans. Single-cell RNA sequencing (scRNA-seq) was performed using follicular aspirates obtained from four IVF patients and identified 13 cell groups: one granulosa cell group, one thecal cell group, 10 subsets of leukocytes, and one group of RBC/platelet. RNA velocity analyses on five granulosa cell populations predicted developmental dynamics denoting two projections of differentiation states. The cell type-specific transcriptomic profiling analyses revealed the presence of a diverse array of leukocyte-derived factors that can directly impact granulosa cell function by activating their receptors (e.g., cytokines and secretory ligands) and are involved in tissue remodeling (e.g., MMPs, ADAMs, ADAMTSs, and TIMPs) and angiogenesis (e.g., VEGFs, PGF, FGF, IGF, and THBS1) in ovulatory follicles. Consistent with the findings from the scRNA-seq data, the leukocyte-specific expression of CD68, IL1B, and MMP9 was verified in follicle tissues collected before and at defined hours after hCG administration from regularly cycling women. Collectively, this study demonstrates that this data can be used as an invaluable resource for identifying important leukocyte-derived factors that promote follicular cell function, thereby facilitating ovulation and luteinization in women.

The effect of C–C motif chemokine ligand 2 supplementation on in vitro maturation of porcine cumulus-oocyte complexes and subsequent developmental competence after parthenogenetic activation

Porcine embryos are used for a variety of applications. However, the maturation rate in vitro remains low, and novel in vitro maturation (IVM) techniques that facilitate the collection of mature oocytes are necessary. C-C motif chemokine ligand 2 (CCL2) is a key periovulatory chemokine present in cumulus-oocyte complexes (COCs). We aimed to examine the effects of CCL2 supplementation during IVM on oocyte maturation and embryonic development. The CCL2 concentration was significantly higher in porcine follicular fluid (pFF) derived from follicles &gt;8 mm in size than in pFF derived from smaller follicles. There was a significant increase in CCL2 mRNA levels in all follicular cells after IVM compared with that before IVM. We analyzed the localization of CCL2 and its receptor, the CCL2 receptor, in follicular cells. During IVM, different concentrations of CCL2 were added to COCs cultured in a maturation medium. After IVM, the group treated with 100 ng/mL CCL2 showed significantly higher metaphase II rates than the control group. All CCL2-treatment groups showed a significant increase in intracellular glutathione levels and a significant decrease in reactive oxygen species levels, compared to the control. In CCs treated with 100 ng/mL CCL2, the mRNA levels of BAX, CASP3, and NPR2 were significantly decreased. Furthermore, the mRNA levels of SOD1, SOD2, and CD44 were significantly increased. In oocytes treated with 10 ng/mL CCL2, mRNA levels of BAX and CASP3 were significantly decreased, whereas, NRF2 and NPM2 were significantly increased. ERK1 exhibited significantly increased mRNA expression in both CCs and oocytes treated with 10 ng/mL CCL2. The protein expression ratio of phosphorylated ERK1/2 to total ERK1/2 was significantly increased in CCs treated with 10 ng/mL CCL2. After parthenogenetic activation, cleavage rates were significantly improved in the 100 ng/mL CCL2 treatment group, and blastocyst formation rates were significantly enhanced in the 10 ng/mL CCL2 treatment group. Overall, our results suggest that IVM medium along with CCL2 improves porcine oocyte maturation and the development of parthenogenetically-activated embryos.

AREG upregulates secreted protein acidic and rich in cysteine expression in human granulosa cells

The secreted protein acidic and rich in cysteine (SPARC) is a secreted glycoprotein and the expression of ovarian SPARC peaks during ovulation and luteinization. Besides, SPARC expression was induced by human chorionic gonadotropin (hCG) in rat granulosa cells. Amphiregulin (AREG) is the most abundant epidermal growth factor receptor (EGFR) ligand expressed in human granulosa cells and follicular fluid. AREG mediates the physiological functions of luteinizing hormone (LH)/hCG in the ovary. However, to date, the biological function of SPARC in the human ovary remains undetermined, and whether AREG regulates SPARC expression in human granulosa cells is unknown. In this study, we show that AREG upregulated SPARC expression via EGFR in a human granulosa-like tumor cell line, KGN. Treatment of AREG activated ERK1/2, JNK, p38 MAPK, and PI3K/AKT signaling pathways and all of them were required for the AREG-induced SPARC expression. Using RNA-sequencing, we identified that steroidogenic acute regulatory protein (StAR) was a downstream target gene of SPARC. In addition, we demonstrated that SPARC mRNA levels were positively correlated with the levels of StAR mRNA in the primary culture of human granulosa cells. Moreover, SPARC protein levels were positively correlated with progesterone levels in follicular fluid of in vitro fertilization patients. This study provides the regulatory role of AREG on the expression of SPARC and reveals the novel function of SPARC in progesterone production in granulosa cells.

A single-cell transcriptomic atlas characterizes cell types and their molecular features in yak ovarian cortex

The ovary as one of the most dynamic organs produces steroids to orchestrate female secondary sexual characteristics, harbors ovarian reserve for oocytes, releases mature oocytes for fertilization, and maintains pregnancy. Yak (Bos grunniens) is the only bovid animal that can adapt to the harsh climatic conditions on the Qinghai-Tibetan Plateau (altitudes of over 3000 m above sea level). However, the cellular atlas is composed of oocytes and other somatic cells, and their individual molecular characteristics remain to be elucidated in the yak ovary. Here, single-cell RNA sequencing (scRNA-seq) was performed to delineate the molecular signature of various cell types in the yak ovarian cortex. A cellular atlas of yak ovarian cortex was constructed successfully on the basis of the differentially expressed genes (DEGs) from the distinct cell types and their functional enrichment analysis, comprising endothelial cells, nature kill cells, stromal cells, smooth muscle cells, oocytes, macrophages, epithelial cells, and granulosa cells. Meanwhile, the signature genes were determined based on their expression specificity in each cell type. A cell-to-cell communication network was built in light of the differentially overexpressed ligand and receptor genes from each cell type. Further, the oocytes were subdivided into four subtypes based on their individual DEGs and the functional enrichment of the DEGs. FST and TOP2A were identified as maker genes for oocytes by immunostaining in the yak ovarian cortex. The cellular atlas reveals the biological characteristics of the ovarian cortex at the cellular molecular level and provides insights into female reproductive biology via cellular communications in the yak.

Oncostatin M and its receptor in women with polycystic ovary syndrome and association with assisted reproductive technology outcomes

The role of adipokines in ovarian-related disorders such as polycystic ovary syndrome (PCOS) has been reported. However, the involvement of Oncostatin M (OSM), a recently identified adipokine, in ovarian function is unknown. Therefore, we investigated the association of the OSM signaling pathway with ovarian functions and PCOS pathogenesis. This case-control study enrolled 30 PCOS and 30 healthy women who underwent the intracytoplasmic sperm injection procedure. OSM and OSM receptor (OSMR) levels were evaluated in the follicular fluid (FF). Moreover, the expression of insulin receptor substrates (IRS1 and IRS2), OSM, OSMR, suppressor of cytokine signaling 3 (SOCS3), and androgen receptor (AR) genes were analyzed in the isolated cumulus cells (CCs). For the in-vitro experiment, the effect of recombinant OSM on the expression of related genes in isolated CCs was analyzed. Follicular concentrations of OSM and OSMR were significantly lower in PCOS (123.91±48.58 pg/ml and 0.93±0.35 ng/ml, respectively) compared to control women (283.53 ± 96.62 pg/ml and 1.45 ± 0.18 ng/ml, respectively; p < 0.001) and were positively correlated with the oocyte maturation (r = 0.611 and r = 0.611, respectively) and fertilization (r = 0.592 and r = 0.627, respectively) rates in the PCOS group. Furthermore, the SOCS3 expression was upregulated about eight times in PCOS patients compared to the controls (p < 0.05). The treatment of cells with recombinant OSM significantly increased SOCS3, OSMR, IRS-1, and -2 expression and decreased AR expression. The decreased levels of OSM and its receptor in PCOS patients, possibly mediated by SOCS3, could negatively affect oocyte maturation and fertilization rates.

Differential proteomic analysis demonstrates follicle fluid participate immune reaction and protein translation in yak

Background

Ovarian follicle fluid (FF) as a microenvironment surrounding oocyte plays critical roles in physio-biochemical processes of follicle development and oocyte maturation. It is hypothesized that proteins in yak FF participate in the physio-biochemical pathways. The primary aims of this study were to find differentially expressed proteins (DEPs) between mature and immature FF, and to elucidating functions of the mature and immature FF in yak.

Results

The mature and immature FF samples were obtained from three healthy yaks that were nonpregnant, aged from four to five years, and free from any anatomical reproductive disorders. The FF samples were subjected to mass spectrometry with the isobaric tags for relative and absolute quantification (iTRAQ). The FF samples went through correlation analysis, principle component analysis, and expression pattern analysis based on quantification of the identified proteins. Four hundred sixty-three DEPs between mature and immature FF were identified. The DEPs between the mature and immature FF samples underwent gene ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG), and protein-protein interaction (PPI) analysis. The DEPs highly expressed in the mature FF mainly took parts in the complement and coagulation cascades, defense response, acute-phase response, response to other organism pathways to avoid invasion of exogenous microorganisms. The complement activation pathway contains eight DEPs, namely C2, C5, C6, C7, C9, C4BPA, CFH, and MBL2. The three DEPs, CATHL4, CHGA, and PGLYRP1, take parts in defense response pathway to prevent invasion of exogenetic microorganism. The coagulation cascades pathway involves many coagulation factors, such as F7, F13A1, FGA, FGB, FGG, KLKB1, KNG1, MASP1, SERPINA1, and SERPIND1. While the DEPs highly expressed in the immature FF participated in protein translation, peptide biosynthetic process, DNA conformation change, and DNA geometric change pathways to facilitate follicle development. The translation pathway contains many ribosomal proteins, such as RPL3, RPL5, RPS3, RPS6, and other translation factors, such as EIF3J, EIF4G2, ETF1, MOV10, and NARS. The DNA conformation change and DNA geometric change involve nine DEPs, DDX1, G3BP1, HMGB1, HMGB2, HMGB3, MCM3, MCM5, MCM6, and RUVBL2. Furthermore, the expressed levels of the main DEPs, C2 and SERPIND1, were confirmed by western blot.

Conclusions

The differential proteomics revealed the up-regulated DEPs in mature FF take parts in immunoreaction to prevent invasion of microorganisms and the up-regulated DEPs in immature FF participate in protein synthesis, which may improve our knowledge of the follicular microenvironment and its biological roles for reproductive processes in yak. The DEPs, C2 and SERPIND1, can be considered as protein markers for mature yak follicle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-021-03097-0.

High-definition ultrasonography in the evaluation of the reproductive tract of bitches during the follicular phase of the estrous cycle

The aim was to sonographically evaluate the reproductive tract of bitches during the follicular phase of the estrous cycle using High Density (HD) ultrasonic techniques. Females (n = 8) were evaluated at five different times throughout the follicular phase, as determined by vaginal cytology and blood progesterone concentrations. Ultrasonic exams were performed using the ACUSON S2000/SIEMENS device utilizing a multifrequency HD transducer (5.5-18 MHz). Videos of the ovaries were obtained and recordings were evaluated using a DICOM viewer software for counting and measuring the ovarian structures, which were assigned to groups based on diameter in mm: G1: ≤ 1; G2: from 1.01 to 3.5; G3 from 3.51 to 5.5; G4: from 5.51 to 10. There was a greater uterine thickness with the progression of the follicular phase (P < 0.05). Six distinct regions were identified in the uterine wall. The ovarian dimensions increased (P < 0.05) as stage of the follicular phase advanced. There was fluid detected around the ovaries after ovulation. There was a characteristic fat tissue hyperechogenicity around the ovaries at all timepoints. There was a difference in the number of ovarian structures of each dimension group at each time there were assessments (P < 0.05). There was a difference in diameter of the largest ovarian structure and in average value of wall thickness at all timepoints when there were evaluations (P < 0.05). The HD ultrasonography technique provides for excellent image resolution, allowing for a more precise characterization of the bitch's reproductive structures and changes occurring during the follicular phase of the estrous cycle.

WNT5A Enhances LH-Mediated Expression of HAS2 in Granulosa Cells

In adult ovary, WNT5A is involved in follicular responses to gonadotropins and necessary for ovarian follicle development. However, the mechanism by which gonadotropins regulate WNT5A and the role of WNT5A in modulating follicular responses to gonadotropins are unclear. In mice, we discovered that the expression of Wnt5a was increased in granulosa cells of mouse ovaries during ovulation, and regulated by gonadotropin-activated intra-ovarian cytokine interleukin 6 (IL6). Using human granulosa-like KGN cells, we confirmed that forskolin plus phorbol myristate acetate (PMA) which mimic the luteinizing hormone (LH) action induced the expression of WNT5A and cumulus expansion gene HAS2. However, this effect was suppressed by a NF-κB pathway inhibitor. Inhibition of NF-κB pathway also blocked the activation of WNT5A signaling components ROR2 and JNK. Moreover, exogenous WNT5A enhanced the expression of HAS2 in KGN cells through JNK and AKT signaling pathways. Knockdown of WNT5A expression by siRNA disrupted LH-mediated expression of HAS2. Our findings indicate that WNT5A could be a fine tuner for LH-induced HAS2 expression in ovarian granulosa cells.

Expression Profile of New Gene Markers and Signaling Pathways Involved in Immunological Processes in Human Cumulus-Oophorus Cells

The function of the immune system extends from defense against external pathogens to the recognition and elimination of mutated or dying cells, aiding elimination of malignant potential and/or maintaining homeostasis. The many cell types of the immune system secrete a broad range of factors to enable cellular signaling that is vital to physiological processes. Additionally, in the ovary, follicular selection and maturation, as well as ovulation, are directly regulated by the nearby immune cells. Additionally, ovulation and rupture of the follicle have been observed to resemble a local inflammatory response. Cells of the cumulus-oocyte complex (COC) show evolving gene expression profiles throughout the oocytes' lifespan, including genes associated with immunological processes. Analysis of these genes allows the identification of useful molecular markers, as well as highlighting gene functions and interactions in these cells. Cumulus cells were obtained from hormonally stimulated patients undergoing an in vitro fertilization procedure and studied under long-term culture conditions. The microarray technique made it possible to compare the level of CCs' gene expression on the 1st, 7th, 15th and 30th day of cultivation. Additionally, RNA microarray analysis was performed to map gene expression in these cells, associated with immunological processes and associated cytokine signaling. Subsequently, the use of DAVID software allowed us to identify the "defense response to other organism", "defense response", "defense response to virus", "cytokine secretion", "cytokine production" and "cytokine-mediated signaling pathway" GO BP terms, as well as allowing further analysis of the most differentially expressed genes associated with these processes. Of the 122 genes involved, 121 were upregulated and only one was downregulated. The seven most upregulated genes related to the abovementioned terms were ANXA3, IFIT1, HLA-DPA1, MX1, KRT8, HLA-DRA and KRT18. Therefore, genes involved in immunological defense processes are upregulated in CC cultures and could serve as useful molecular markers of growth and development in the COC, as well as the proliferation of granulosa and cumulus cells.

Follicular hormone dynamics during the midcycle surge of gonadotropins in women undergoing fertility treatment

Changes in concentrations of intra-follicular hormones during ovulation are important for final oocyte maturation and endometrial priming to ensure reproductive success. As no human studies have investigated these changes in detail, our objective was to describe the dynamics of major follicular fluid (FF) hormones and transcription of steroidogenic enzymes and steroid receptors in human granulosa cells (GCs) during ovulation. We conducted a prospective cohort study at a public fertility clinic in 2016-2018. Fifty women undergoing ovarian stimulation for fertility treatment were included. From each woman, FF and GCs were collected by transvaginal ultrasound-guided follicle puncture of one follicle at two specific time points during ovulation, and the study covered a total of five time points: before ovulation induction (OI), 12, 17, 32 and 36 h after OI. Follicular fluid concentrations of oestradiol, progesterone, androstenedione, testosterone, 17-hydroxyprogesterone, anti-Mullerian hormone, inhibin A and inhibin B were measured using ELISA assays, and a statistical mixed model was used to analyse differences in hormone levels between time points. Gene expression of 33 steroidogenic enzymes and six hormone receptors in GCs across ovulation were assessed by microarray analysis, and selected genes were validated by quantitative reverse transcription PCR. We found that concentrations of oestradiol, testosterone, progesterone, AMH, inhibin A and inhibin B (P < 0.001) and gene expression of 12 steroidogenic enzymes and five receptors (false discovery rate < 0.0001) changed significantly during ovulation. Furthermore, we found parallel changes in plasma hormones. The substantial changes in follicular hormone production during ovulation highlight their importance for reproductive success.

Two waves of transcriptomic changes in periovulatory human granulosa cells

STUDY QUESTION

How does the human granulosa cell (GC) transcriptome change during ovulation?

SUMMARY ANSWER

Two transcriptional peaks were observed at 12 h and at 36 h after induction of ovulation, both dominated by genes and pathways known from the inflammatory system.

WHAT IS KNOWN ALREADY

The crosstalk between GCs and the oocyte, which is essential for ovulation and oocyte maturation, can be assessed through transcriptomic profiling of GCs. Detailed transcriptional changes during ovulation have not previously been assessed in humans.

STUDY DESIGN, SIZE, DURATION

This prospective cohort study comprised 50 women undergoing fertility treatment in a standard antagonist protocol at a university hospital-affiliated fertility clinic in 2016-2018.

PARTICIPANTS/MATERIALS, SETTING, METHODS

From each woman, one sample of GCs was collected by transvaginal ultrasound-guided follicle aspiration either before or 12 h, 17 h or 32 h after ovulation induction (OI). A second sample was collected at oocyte retrieval, 36 h after OI. Total RNA was isolated from GCs and analyzed by microarray. Gene expression differences between the five time points were assessed by ANOVA with a random factor accounting for the pairing of samples, and seven clusters of protein-coding genes representing distinct expression profiles were identified. These were used as input for subsequent bioinformatic analyses to identify enriched pathways and suggest upstream regulators. Subsets of genes were assessed to explore specific ovulatory functions.

MAIN RESULTS AND THE ROLE OF CHANCE

We identified 13 345 differentially expressed transcripts across the five time points (false discovery rate, <0.01) of which 58% were protein-coding genes. Two clusters of mainly downregulated genes represented cell cycle pathways and DNA repair. Upregulated genes showed one peak at 12 h that resembled the initiation of an inflammatory response, and one peak at 36 h that resembled the effector functions of inflammation such as vasodilation, angiogenesis, coagulation, chemotaxis and tissue remodelling. Genes involved in cell-matrix interactions as a part of cytoskeletal rearrangement and cell motility were also upregulated at 36 h. Predicted activated upstream regulators of ovulation included FSH, LH, transforming growth factor B1, tumour necrosis factor, nuclear factor kappa-light-chain-enhancer of activated B cells, coagulation factor 2, fibroblast growth factor 2, interleukin 1 and cortisol, among others. The results confirmed early regulation of several previously described factors in a cascade inducing meiotic resumption and suggested new factors involved in cumulus expansion and follicle rupture through co-regulation with previously described factors.

LARGE SCALE DATA

The microarray data were deposited to the Gene Expression Omnibus (www.ncbi.nlm.nih.gov/gds/, accession number: GSE133868).

LIMITATIONS, REASONS FOR CAUTION

The study included women undergoing ovarian stimulation and the findings may therefore differ from a natural cycle. However, the results confirm significant regulation of many well-established ovulatory genes from a series of previous studies such as amphiregulin, epiregulin, tumour necrosis factor alfa induced protein 6, tissue inhibitor of metallopeptidases 1 and plasminogen activator inhibitor 1, which support the relevance of the results.

WIDER IMPLICATIONS OF THE FINDINGS

The study increases our understanding of human ovarian function during ovulation, and the publicly available dataset is a valuable resource for future investigations. Suggested upstream regulators and highly differentially expressed genes may be potential pharmaceutical targets in fertility treatment and gynaecology.

STUDY FUNDING/COMPETING INTEREST(S)

The study was funded by EU Interreg ÔKS V through ReproUnion (www.reprounion.eu) and by a grant from the Region Zealand Research Foundation. None of the authors have any conflicts of interest to declare.

Serum amyloid A in polycystic ovary syndrome

Polycystic ovary syndrome (PCOS), i.e., anovulation, hyperandrogenemia and polycystic ovary, is an endocrine-metabolic disease affecting reproductive aged women. Women with PCOS are likely to develop obesity, dyslipidemia, type 2 diabetes mellitus (T2DM) and cardiovascular diseases at a younger age. Despite high frequency and severe disease burden, the pathophysiological mechanisms of PCOS remain poorly defined and correspondingly have no therapeutic options. Emerging evidence has demonstrated that PCOS is accompanied with low-grade chronic inflammation and biomarkers thereof. Interestingly, serum amyloid A (SAA) has recently been identified as a potential marker of infection and inflammation and a number of studies have reported an association with PCOS. In this review, we explore the relationship between SAA and hyperandrogenemia, inflammation, obesity and insulin resistance, and provide convincing evidence for SAA as a potential inflammatory biomarker in PCOS.

Effects of flame retardants on ovarian function

Flame retardants have been added to a variety of consumer products and are now found ubiquitously throughout the environment. Epidemiological, in vivo, and in vitro studies have shown that polybrominated diphenyl ether (PBDE) flame retardants may have a negative impact on human health; this has resulted in their phase-out and replacement by alternative flame retardants, such as hexabromocyclododecane (HBCDD), tetrabromobisphenol A (TBBPA), and organophosphate esters (OPEs). Evidence suggests that some of these chemicals induce ovarian dysfunction and thus may be detrimental to female fertility; however, the effects of many of these flame retardants on the ovary remain unclear. In this review, we present an overview of the effects of brominated and organophosphate ester flame retardants on ovarian function and discuss the possible mechanisms which may mediate these effects.

Transcriptome analysis of buffalo granulosa cells in three dimensional culture systems

Hanging drop (HD) three-dimensional (3D) culture model for buffalo granulosa cells (GC) was reported to mimic the preovulatory stage of ovarian follicles in our previous study. To further verify its reliability, the present study attempted a comparative transcriptome profile of buffalo GC freshly isolated from ovarian follicles (<8 mm diameter) (FC) and their cultures in normal culture dish (ND or 2D), polyHEMA coated dish (PH) and HD culture systems (3D). Out of 223 significantly (-log2 fold change: >3; p < .0005; false discovery rate [FDR]: <0.1) differentially expressed genes (SDEGs) among different culture systems, 137 were found unannotated, and 94, 29, and 66 were exclusively expressed in FC, PH, and HD, respectively. However, on eliminating the fixed points of p values and FDR from the entire raw data, only 11 genes related to long noncoding RNA, 12 genes related to luteinization, and 3 genes related to follicular maturation were exclusively expressed in FC, PH, and HD culture systems, respectively. The quantitative real time-PCR validation and the next generation sequencing data had more than 90% correlation. Bioinformatics analyses of the exclusively expressed SDEG revealed that the freshly aspirated GCs were a true representative of GCs from small follicles (<8 mm diameter), the GC spheroids under PH maintained mitochondrial function, and those cultured in HD system for 6 days simulated the inflammatory milieu required for ovulation. Therefore, the comparative transcriptome profile also reinforced that HD culture system is better in vitro culture method than the other methods analyzed in this study for buffalo GC.

Polybrominated Diphenyl Ethers in Human Follicular Fluid Dysregulate Mural and Cumulus Granulosa Cell Gene Expression

Polybrominated diphenyl ethers (PBDEs), a major class of flame retardants incorporated into numerous consumer products, leach out into dust resulting in widespread exposure. There is evidence from in vitro and in vivo animal studies that PBDEs affect ovarian granulosa cell function and follicular development, yet human studies of their association with female infertility are inconclusive. Here, we tested the hypothesis that exposure to the PBDEs in follicular fluid is associated with dysregulation of gene expression in the mural and cumulus granulosa cells collected from women undergoing in vitro fertilization by intracytoplasmic sperm injection. The median concentration of the ∑ 10PBDEs detected in the follicular fluid samples (n = 37) was 15.04 pg/g wet weight. RNA microarray analyses revealed that many genes were differentially expressed in mural and cumulus granulosa cells. Highest vs lowest quartile exposure to the Σ 10PBDEs or to 2 predominant PBDE congeners, BDE-47 or BDE-153, was associated with significant effects on gene expression in both cell types. Mural granulosa cells were generally more sensitive to PBDE exposure compared to cumulus cells. Overall, gene expression changes associated with BDE-47 exposure were similar to those for ∑ 10PBDEs but distinct from those associated with BDE-153 exposure. Interestingly, exposure to BDE-47 and ∑ 10PBDEs activated the expression of genes in pathways that are important in innate immunity and inflammation. To the best of our knowledge, this is the first demonstration that exposure to these environmental chemicals is associated with the dysregulation of pathways that play an essential role in ovulation.

Effect of luteinizing hormone concentration on transcriptome and subcellular organelle phenotype of ovarian granulosa cells

RESEARCH QUESTION

Granulosa cells (GCs) surrounding oocytes are crucial for follicular growth, oocyte development, ovulation, and luteinization under the dynamic co-stimulation of follicle stimulating hormone (FSH) and luteinizing hormone (LH). This study aimed to investigate the effect of LH levels on GCs in preovulatory follicles under gonadotropin releasing hormone antagonist-based ovarian stimulation. In vitro experiments were also conducted to study the direct effect of LH on GCs.

METHODS

Twelve infertile women were divided into low (L), medium (M), and high (H) LH groups according to their serum LH levels during ovarian stimulation. RNA-sequencing (RNA-seq) was conducted to examine the transcriptome profiles of GCs obtained from the above patients during the oocyte retrieval. The activity of mitochondrial dehydrogenase was measured under the stimulation of recombinant LH (rLH) concentration gradient combined with recombinant FSH. The ultrastructures of subcellular organelles were observed.

RESULTS

Bioinformatic analyses showed that compared with the M group, molecule and pathway changes in the L group and in the H group were similar. In cultured GCs, both insufficient and excessive rLH impaired the activity of mitochondrial dehydrogenase. With the medium rLH concentration, numerous cell connections and abundant mitochondria and liposomes were observed. Compared with the medium concentration, GCs showed smaller and rounder mitochondria, more autophagosomes, and massive organelles damages with excessive rLH, and swollen, circular, or forked mitochondria were observed with inadequate rLH.

CONCLUSIONS

RNA-seq provided a novel spectrum of transcriptome characteristics of GCs potentially affected by serum LH levels during ovarian stimulation. In vitro, rLH could directly affect GCs at the subcellular level.

Patients who failed to conceive following an in vitro fertilization cycle can be clustered into different failure causes using gene expression hierarchical analysis†

The patient's response to an IVF stimulation protocol is highly variable and thus difficult to predict. When a cycle fails, there are often no apparent or obvious reasons to explain the failure. Having clues on what went wrong during stimulation could serve as a basis to improve and personalize the next protocol. This exploratory study aimed to investigate if it is possible to distinguish different failure causes or different follicular responses in a population of nonpregnant IVF patients. Using qRT-PCR, we analyzed a panel of genes indicative of different failure causes in patients who did not achieve pregnancy following an IVF cycle. For each patient, a pool of follicular cells from all aspirated follicles was used as a sample which gives a global picture of the patient's ovary and not a specific picture of each follicle. We performed hierarchical clustering analysis to split the patients according to the gene expression pattern. Hierarchical analysis showed that the population of nonpregnant IVF patients could be divided into three clusters. Gene expression was significantly different, and each cluster displayed a particular gene expression pattern. Follicular cells from patients in clusters 1, 2 and 3 displayed respectively a pattern of gene expression related to large incompetent follicles with a higher apoptosis (over matured), to follicles not ready to ovulate (under mature) and to an excess of inflammation with no visible symptoms. This study reinforces the idea that women often have different response to the same protocol and would benefit from more personalized treatments.

Expression and Contribution of NLRP3 Inflammasome During the Follicular Development Induced by PMSG

Follicular development and following ovulation induced by luteinizing hormone (LH) surge are critical for ovarian functions, but the molecular mechanism regulating ovarian ovulation attracts more attention and remains mainly unknown. Recent researches on the nucleotide leukin rich polypeptide 3 (NLRP3) inflammasome shred light on it. Given pregnant mare serum gonadotropin (PMSG) can not only trigger the follicular development, but also induce the following ovulation, the present study therefore examined that expression and localization of NLRP3 inflammasome through immunohistochemistry and Western blotting during the follicular development induced by PMSG. The results showed expressions of NLRP3 and the adaptor protein apoptosis-associated speck-like protein (ASC) significantly increased in the outside of intrafollicular fluid, further analysis found that caspase-1 was activated and IL-1β production was also upregulated after 52 h-treatment of PMSG. Furthermore, a significant increase of ovulation-related genes, hypoxia inducible factor (HIF)-1α and endothelin (ET)-1, was found after 52 h-treatment of PMSG. To our knowledge, it is the first time to clearly indicated the activation of NLRP3 inflammasome may contribute to the ovulation of PMSG-treated ovaries, which will help to further clarify the ovulatory mechanism in mammals.

Contribution of the immune system to follicle differentiation, ovulation and early corpus luteum formation

Much of what we know about the involvement of the immune system in periovulatory follicle differentiation, ovulation and subsequent formation of the corpus luteum in cattle is drawn from the findings of studies in several mammalian livestock species. By integrating published histological data from cattle, sheep and pigs and referring back to the more comprehensive knowledge bank that exists for mouse and humans we can sketch out the key cells of the immune system and the cytokines and growth factors that they produce that are involved in follicle differentiation and luteinization, ovulation and early follicle development. These contributions are reviewed and the key findings, discussed.

Progressive changes in human follicular fluid composition over the course of ovulation: quantitative proteomic analyses

Follicular fluid (FF) acts as a vehicle for paracrine signalling between somatic cells of the follicle and the oocyte. To investigate changes in the protein composition of FF during ovulation, we conducted a prospective cohort study including 25 women undergoing fertility treatment. Follicular fluid was aspirated either before or 12, 17, 32 or 36 h after induction of ovulation (five patients per time point). Liquid chromatography-mass spectrometry was used to identify and quantify FF proteins. In total, 400 proteins were identified and the levels of 40 proteins changed significantly across ovulation, evaluated by analysis of covariance (adjusted p < 0.05) and on-off expression patterns. The majority peaked after 12-17 h, e.g., AREG (p < 0.0001), TNFAIP6 (p < 0.0001), and LDHB (p = 0.0316), while some increased to peak after 36 h e.g., ACPP (p < 0.0001), TIMP1 (p < 0.0001) and SERPINE1 (p = 0.0002). Collectively, this study highlights proteins and pathways of importance for ovulation and oocyte competence in humans.

Single-cell reconstruction of follicular remodeling in the human adult ovary

The ovary is perhaps the most dynamic organ in the human body, only rivaled by the uterus. The molecular mechanisms that regulate follicular growth and regression, ensuring ovarian tissue homeostasis, remain elusive. We have performed single-cell RNA-sequencing using human adult ovaries to provide a map of the molecular signature of growing and regressing follicular populations. We have identified different types of granulosa and theca cells and detected local production of components of the complement system by (atretic) theca cells and stromal cells. We also have detected a mixture of adaptive and innate immune cells, as well as several types of endothelial and smooth muscle cells to aid the remodeling process. Our results highlight the relevance of mapping whole adult organs at the single-cell level and reflect ongoing efforts to map the human body. The association between complement system and follicular remodeling may provide key insights in reproductive biology and (in)fertility.

Perturbations in Lineage Specification of Granulosa and Theca Cells May Alter Corpus Luteum Formation and Function

Anovulation is a major cause of infertility, and it is the major leading reproductive disorder in mammalian females. Without ovulation, an oocyte is not released from the ovarian follicle to be fertilized and a corpus luteum is not formed. The corpus luteum formed from the luteinized somatic follicular cells following ovulation, vasculature cells, and immune cells is critical for progesterone production and maintenance of pregnancy. Follicular theca cells differentiate into small luteal cells (SLCs) that produce progesterone in response to luteinizing hormone (LH), and granulosa cells luteinize to become large luteal cells (LLCs) that have a high rate of basal production of progesterone. The formation and function of the corpus luteum rely on the appropriate proliferation and differentiation of both granulosa and theca cells. If any aspect of granulosa or theca cell luteinization is perturbed, then the resulting luteal cell populations (SLC, LLC, vascular, and immune cells) may be reduced and compromise progesterone production. Thus, many factors that affect the differentiation/lineage of the somatic cells and their gene expression profiles can alter the ability of a corpus luteum to produce the progesterone critical for pregnancy. Our laboratory has identified genes that are enriched in somatic follicular cells and luteal cells through gene expression microarray. This work was the first to compare the gene expression profiles of the four somatic cell types involved in the follicle-to-luteal transition and to support previous immunofluorescence data indicating theca cells differentiate into SLCs while granulosa cells become LLCs. Using these data and incorporating knowledge about the ways in which luteinization can go awry, we can extrapolate the impact that alterations in the theca and granulosa cell gene expression profiles and lineages could have on the formation and function of the corpus luteum. While interactions with other cell types such as vascular and immune cells are critical for appropriate corpus luteum function, we are restricting this review to focus on granulosa, theca, and luteal cells and how perturbations such as androgen excess and inflammation may affect their function and fertility.