Integrated Analysis of Transcriptome and Histone Modifications in Granulosa Cells During Ovulation in Female Mice

Granulosa cell expression of Fos is critical for regulating ovulatory gene expressions in the mouse ovary

A previous study showed that female Fos null mice fail to ovulate even when given gonadotropins, suggesting that ovarian expression of Fos is critical for successful ovulation. However, the expression of FOS and function of FOS have not been determined in the mouse ovary. FOS, a member of the Fos family (Fos, Fosb, Fosl1, and Fosl2), functions as a transcription factor by forming a heterodimer complex with a member of Jun family (Jun, Junb, and Jund). This study demonstrated rapid increases in Fos, along with other Fos and Jun family members, after hCG administration in the ovary of immature PMSG-primed mice and after the LH surge in naturally cycling animals. ChIP-seq analysis identified 1965 FOS-binding genes in granulosa cells collected at 3 h post-hCG, including Pgr, Ptgs2, Tnfiap6, and Edn2, genes known to be involved in the ovulatory process. When super-ovulation was induced, the number of oocytes released was significantly reduced in Esr2cre/+-driven granulosa cell-specific Fos knockout (gcFosKO) mice. This reduction was accompanied by lower expressions of Pgr, Ptgs2, Ptgs1, and Edn2 in preovulatory follicles of gcFosKO mice compared to those in control littermates. In addition, gcFosKO mice showed a trend toward a decreased average litter size. Together, the present study indicates that the preovulatory induction of Fos expression is crucial for increasing the expression of key ovulatory genes, yet the role of FOS may be partially substituted by other Fos and Jun family members induced in the preovulatory follicle in the gcFosKO mouse ovary.

New insights into the ovulatory process in the human ovary

BACKGROUND

Successful ovulation is essential for natural conception and fertility. Defects in the ovulatory process are associated with various conditions of infertility or subfertility in women. However, our understanding of the intra-ovarian biochemical mechanisms underlying this process in women has lagged compared to our understanding of animal models. This has been largely due to the limited availability of human ovarian samples that can be used to examine changes across the ovulatory period and delineate the underlying cellular/molecular mechanisms in women. Despite this challenge, steady progress has been made to improve our knowledge of the ovulatory process in women by: (i) collecting granulosa cells across the IVF interval, (ii) creating a novel approach to collecting follicular cells and tissues across the periovulatory period from normally cycling women, and (iii) developing unique in vitro models to examine the LH surge or hCG administration-induced ovulatory changes in gene expression, the regulatory mechanisms underlying the ovulatory changes, and the specific functions of the ovulatory factors.

OBJECTIVE AND RATIONALE

The objective of this review is to summarize findings generated using in vivo and in vitro models of human ovulation, with the goal of providing new insights into the mechanisms underlying the ovulatory process in women.

SEARCH METHODS

This review is based on the authors' own studies and a search of the relevant literature on human ovulation to date using PubMed search terms such as 'human ovulation EGF-signaling', 'human ovulation steroidogenesis', 'human ovulation transcription factor', 'human ovulation prostaglandin', 'human ovulation proteinase', 'human ovulation angiogenesis' 'human ovulation chemokine', 'human ovulatory disorder', 'human granulosa cell culture'. Our approach includes comparing the data from the authors' studies with the existing microarray or RNA-seq datasets generated using ovarian cells obtained throughout the ovulatory period from humans, monkeys, and mice.

OUTCOMES

Current findings from studies using in vivo and in vitro models demonstrate that the LH surge or hCG administration increases the expression of ovulatory mediators, including EGF-like factors, steroids, transcription factors, prostaglandins, proteolytic systems, and other autocrine and paracrine factors, similar to those observed in other animal models such as rodents, ruminants, and monkeys. However, the specific ovulatory factors induced, their expression pattern, and their regulatory mechanisms vary among different species. These species-specific differences stress the necessity of utilizing human samples to delineate the mechanisms underlying the ovulatory process in women.

WIDER IMPLICATIONS

The data from human ovulation in vivo and in vitro models have begun to fill the gaps in our understanding of the ovulatory process in women. Further efforts are needed to discover novel ovulatory factors. One approach to address these gaps is to improve existing in vitro models to more closely mimic in vivo ovulatory conditions in humans. This is critically important as the knowledge obtained from these human studies can be translated directly to aid in the diagnosis of ovulation-associated pathological conditions, for the development of more effective treatment to help women with anovulatory infertility or, conversely, to better manage ovulation for contraceptive purposes.

REGISTRATION NUMBER

N/A.

LH-stimulated periodic lincRNA HEOE regulates follicular dynamics and influences estrous cycle and fertility via miR-16-ZMAT3 and PGF2α in pigs

Disruption of the estrous cycle affects fertility and reproductive health. Follicular dynamics are key to the regularity of the estrous cycle. We identified a novel lincRNA, HEOE, showing significant upregulation in the ovaries during the estrus phase across various pig breeds. Functional analysis revealed that HEOE is responsive to luteinizing hormone (LH) stimulation, modulating transcriptional suppression and alternative splicing in ovarian granulosa cells (GCs). This leads to increased GC apoptosis and inhibition of proliferation. Mechanistically, HEOE inhibits miR-16 maturation in the nucleus, and sequesters miR-16 in the cytoplasm, thereby collectively reducing miR-16's inhibition on ZMAT3, enhancing the expression of ZMAT3, a key factor in the p53 pathway and alternative splicing, thereby regulating follicular development. This effect was validated in both mice and pig follicles. Persistent overexpression or suppression of HEOE throughout the estrous cycle impairs cycle regularity and reduces litter size. These outcomes are associated with HEOE reduced follicular PGF2α levels and modulation of the cAMP signaling pathway. Our data, combined with public databases, indicate that the high expression of HEOE during the estrus phase is crucial for maintaining the estrous cycle. HEOE is a potential therapeutic target for regulating fertility and ensuring estrous cycle regularity in pigs.

Differential gene expression in decidualized human endometrial stromal cells induced by different stimuli

Decidualization can be induced by culturing human endometrial stromal cells (ESCs) with several decidualization stimuli, such as cAMP, medroxyprogesterone acetate (MPA) or Estradiol (E2). However, it has been unclear how decidualized cells induced by different stimuli are different. We compared transcriptomes and cellular functions of decidualized ESCs induced by different stimuli (MPA, E2 + MPA, cAMP, and cAMP + MPA). We also investigated which decidualization stimulus induces a closer in vivo decidualization. Differentially expressed genes (DEGs) and altered cellular functions by each decidualization stimuli were identified by RNA-sequence and gene-ontology analysis. DEGs was about two times higher for stimuli that use cAMP (cAMP and cAMP + MPA) than for stimuli that did not use cAMP (MPA and E2 + MPA). cAMP-using stimuli altered the cellular functions including angiogenesis, inflammation, immune system, and embryo implantation whereas MPA-using stimuli (MPA, E2 + MPA, and cAMP + MPA) altered the cellular functions associated with insulin signaling. A public single-cell RNA-sequence data of the human endometrium was utilized to analyze in vivo decidualization. The altered cellular functions by in vivo decidualization were close to those observed by cAMP + MPA-induced decidualization. In conclusion, decidualized cells induced by different stimuli have different transcriptome and cellular functions. cAMP + MPA may induce a decidualization most closely to in vivo decidualization.

DNA methylation as a window into female reproductive aging

People with ovaries experience reproductive aging as their reproductive function and system declines. This has significant implications for both fertility and long-term health, with people experiencing an increased risk of cardiometabolic disorders after menopause. Reproductive aging can be assessed through markers of ovarian reserve, response to fertility treatment or molecular biomarkers, including DNA methylation. Changes in DNA methylation with age associate with poorer reproductive outcomes, and epigenome-wide studies can provide insight into genes and pathways involved. DNA methylation-based epigenetic clocks can quantify biological age in reproductive tissues and systemically. This review provides an overview of hallmarks and theories of aging in the context of the reproductive system, and then focuses on studies of DNA methylation in reproductive tissues.

Analysis of cell-cell interaction between mural granulosa cells and cumulus granulosa cells during ovulation using single-cell RNA sequencing data of mouse ovary

Purpose

We investigated the interactions between mural granulosa cells (MGCs) and cumulus granulosa cells (CGCs) during ovulation after the LH surge.

Methods

We performed clustering, pseudotime, and interactome analyses utilizing reported single-cell RNA sequencing data of mouse ovary at 6 h after eCG-hCG injection.

Results

Clustering analysis classified granulosa cells into two distinct populations, MGCs and CGCs. Pseudotime analysis divided granulosa cells into before and after the LH surge, and further divided them into two branches, the ovulatory MGCs and the ovulatory CGCs. Interactome analysis was performed to identify the interactions between MGCs and CGCs. Twenty-six interactions were acting from CGCs toward MGCs, involving ovulation and steroidogenesis. Thirty-six interactions were acting from MGCs toward CGCs, involving hyaluronan synthesis. There were 25 bidirectional interactions, involving the EGFR pathway. In addition, we found three novel interactions: Ephrins-Ephs pathway and Wnt-Lrp6 pathway from CGCs to MGCs, associated with steroidogenesis and lipid transport, respectively, and TGF-β-TGFBR1 pathway from MGCs to CGCs, associated with hyaluronan synthesis.

Conclusions

MGCs and CGCs interact with each other in the preovulatory follicle after the LH surge, and their interactions have roles in corpus luteum formation, oocyte maturation, and follicle rupture.

Cellular atlases of ovarian microenvironment alterations by diet and genetically-induced obesity

Obesity, which can arise from genetic or environmental factors, has been shown to cause serious damages to the reproductive system. The ovary, as one of the primary regulators of female fertility, is a complex organ comprised of heterogeneous cell types that work together to maintain a normal ovarian microenvironment (OME). Despite its importance, the effect of obesity on the entire ovary remains poorly documented. In this study, we performed ovary single-cell and nanoscale spatial RNA sequencing to investigate how the OME changed under different kinds of obesity, including high-fat diet (HFD) induced obesity and Leptin ablation induced obesity (OB). Our results demonstrate that OB, but not HFD, dramatically altered the proportion of ovarian granulosa cells, theca-interstitial cells, luteal cells, and endothelial cells. Furthermore, based on the spatial dynamics of follicular development, we defined four subpopulations of granulosa cell and found that obesity drastically disrupted the differentiation of mural granulosa cells from small to large antral follicles. Functionally, HFD enhanced follicle-stimulating hormone (FSH) sensitivity and hormone conversion, while OB caused decreased sensitivity, inadequate steroid hormone conversion, and impaired follicular development. These differences can be explained by the differential expression pattern of the transcription factor Foxo1. Overall, our study provides a powerful and high-resolution resource for profiling obesity-induced OME and offers insights into the diverse effects of obesity on female reproductive disorders.

LH/hCG Regulation of Circular RNA in Mural Granulosa Cells during the Periovulatory Period in Mice

Ovarian follicles undergo a series of dynamic changes following the ovulatory surge of luteinizing hormone including cumulus expansion, oocyte maturation, ovulation, and luteinization. Post-transcriptional gene regulatory events are critical for mediating LH follicular responses, and among all RNA isoforms, circular RNA (circRNA) is one of the most abundant forms present in cells, yet they remain the least studied. Functionally, circRNA can act as miRNA sponges, protein sponges/decoys, and regulators of transcription and translation. In the context of ovarian follicular development, the identity and roles of circRNA are relatively unknown. In the present study, high throughput RNA sequencing of granulosa cells immediately prior to and 4-h after the LH/hCG surge identified 42,381 circRNA originating from 7712 genes. A total of 54 circRNA were identified as differentially expressed between 0-h and 4-h time points (Fold Change ± 1.5, FDR ≤ 0.1), among them 42 circRNA were upregulated and 12 circRNA were downregulated. All differentially expressed circRNA between the 0-h and 4-h groups were subjected to circinteractome analysis and identified networks of circRNA-protein and circRNA-miRNA were further subjected to "micro-RNA target filter analysis" in Ingenuity Pathway Analyses, which resulted in the identification of miRNA targeted mRNAs. A comparison of these circRNA target mRNAs with LH-induced mRNAs identified Runx2, Egfr, Areg, Sult1el, Cyp19a1, Cyp11a1, and Hsd17b1 as targets of circKif2, circVcan, circMast4, and circMIIt10. These newly identified LH/hCG-induced circRNA, their target miRNA and protein networks provide new insights into the complex interactions associated with periovulatory follicular development.

Intrafollicular Retinoic Acid Signaling Is Important for Luteinizing Hormone-Induced Oocyte Meiotic Resumption

It has been clear that retinoic acid (RA), the most active vitamin A (VA) derivative, plays a central role in governing oocyte meiosis initiation. However, it has not been functionally determined if RA participates in luteinizing hormone (LH)-induced resumption from long-lasting oocyte meiotic arrest, which is essential for haploid oocyte formation. In the present study, using well-established in vivo and in vitro models, we identified that intrafollicular RA signaling is important for normal oocyte meiotic resumption. A mechanistic study indicated that mural granulosa cells (MGCs) are the indispensable follicular compartment for RA-prompted meiotic resumption. Moreover, retinoic acid receptor (RAR) is essential for mediating RA signaling to regulate meiotic resumption. Furthermore, we found zinc finger protein 36 (ZFP36) is the transcriptional target of RAR. Both RA signaling and epidermal growth factor (EGF) signaling were activated in MGCs in response to LH surge, and two intrafollicular signalings cooperate to induce rapid Zfp36 upregulation and Nppc mRNA decrease, which is critical to LH-induced meiotic resumption. These findings extend our understanding of the role of RA in oocyte meiosis: RA not only governs meiotic initiation but also regulates LH-induced meiotic resumption. We also emphasize the importance of LH-induced metabolic changes in MGCs in this process.

Application of Single-Cell RNA Sequencing in Ovarian Development

The ovary is a female reproductive organ that plays a key role in fertility and the maintenance of endocrine homeostasis, which is of great importance to women's health. It is characterized by a high heterogeneity, with different cellular subpopulations primarily containing oocytes, granulosa cells, stromal cells, endothelial cells, vascular smooth muscle cells, and diverse immune cell types. Each has unique and important functions. From the fetal period to old age, the ovary experiences continuous structural and functional changes, with the gene expression of each cell type undergoing dramatic changes. In addition, ovarian development strongly relies on the communication between germ and somatic cells. Compared to traditional bulk RNA sequencing techniques, the single-cell RNA sequencing (scRNA-seq) approach has substantial advantages in analyzing individual cells within an ever-changing and complicated tissue, classifying them into cell types, characterizing single cells, delineating the cellular developmental trajectory, and studying cell-to-cell interactions. In this review, we present single-cell transcriptome mapping of the ovary, summarize the characteristics of the important constituent cells of the ovary and the critical cellular developmental processes, and describe key signaling pathways for cell-to-cell communication in the ovary, as revealed by scRNA-seq. This review will undoubtedly improve our understanding of the characteristics of ovarian cells and development, thus enabling the identification of novel therapeutic targets for ovarian-related diseases.

Effects of Melatonin on the Transcriptome of Human Granulosa Cells, Fertilization and Blastocyst Formation

Melatonin is a promising reagent that can improve assisted reproductive technology (ART) outcomes in infertility patients. However, melatonin is not effective for all infertile patients, and it remains unclear for which patients melatonin would be effective. This study examined the effects of melatonin on ART outcomes and examined its mechanisms. Melatonin increased the fertilization rate in patients whose fertilization rates in the previous cycle were less than 50%, but not in patients whose fertilization rates were more than 50% in the previous cycle. Melatonin increased the blastocyst formation rate in patients whose embryo development rates in the previous cycle were less than 50%, but not in patients whose embryo development rates were more than 50% in the previous cycle. To clarify its mechanisms, transcriptome changes by melatonin treatment in granulosa cells (GCs) of the patients were examined by RNA-sequence. Melatonin treatment altered the transcriptomes of GCs of patients with poor ART outcomes so that they were similar to the transcriptomes of patients with good ART outcomes. The altered genes were associated with the inhibition of cell death and T-cell activity, and the activation of steroidogenesis and angiogenesis. Melatonin treatment was effective for patients with poor fertilization rates and poor embryo development rates in the previous ART cycle. Melatonin alters the GCs transcriptome and, thus, their functions, and this could improve the oocyte quality, leading to good ART outcomes.

The granulosa cell response to luteinizing hormone is partly mediated by YAP1-dependent induction of amphiregulin

Background

The LH surge is a pivotal event that triggers multiple key ovarian processes including oocyte maturation, cumulus expansion, follicular wall rupture and luteinization of mural granulosa and theca cells. Recently, LH-dependent activation of the Hippo signaling pathway has been shown to be required for the differentiation of granulosa cells into luteal cells. Still, the precise interactions between Hippo and LH signaling in murine granulosa cells remain to be elucidated.

Methods

To detect the expression of effectors of the Hippo pathway, western blot, immunohistochemical and RT-qPCR analyses were performed on granulosa cells treated with LH in vitro or isolated from immature mice treated with eCG and hCG. Cultured granulosa cells were pretreated with pharmacologic inhibitors to identify the signaling pathways involved in Hippo regulation by LH. To study the roles of Yap1 and Taz in the regulation of the LH signaling cascade, RT-qPCR and microarray analyses were done on granulosa cells from Yap1^f/f;Taz^f/f mice treated with an adenovirus to drive cre expression. RT-qPCR was performed to evaluate YAP1 binding to the Areg promoter following chromatin immunoprecipitation of granulosa cells collected from mice prior to or 60 min following hCG treatment.

Results

Granulosa cells showed a transient increase in LATS1, YAP1 and TAZ phosphorylation levels in response to the ovulatory signal. This Hippo activation by LH was mediated by protein kinase A. Furthermore, Yap1 and Taz are required for the induction of several LH target genes such as Areg, Pgr and Ptgs2, and for the activation of the ERK1/2 pathway. Consistent with these results, there was a substantial overlap between genes that are upregulated by LH and those that are downregulated following loss of Yap1/Taz, highlighting a major role for Hippo in mediating LH actions in the ovulation process. Finally, we showed that there is a marked recruitment of YAP1 to the Areg promoter of granulosa cells in response to hCG stimulation.

Conclusions

Overall, these results indicate that Hippo collaborates with the cAMP/PKA and ERK1/2 pathways to participate in the precise regulation of the LH cascade, and that Areg, as a direct transcriptional target of YAP1, is involved in mediating its actions in the ovary.