Single-cell RNA sequencing analysis of mouse follicular somatic cells†

Long non-coding ribonucleic acid SNHG18 induced human granulosa cell apoptosis via disruption of glycolysis in ovarian aging

BACKGROUND

In-depth understanding of dynamic expression profiles of human granulosa cells (GCs) during follicular development will contribute to the diagnostic and targeted interventions for female infertility. However, genome-scale analysis of long non-coding ribonucleic acid (lncRNA) in GCs across diverse developmental stages is challenging. Meanwhile, further research is needed to determine how aberrant lncRNA expression participates in ovarian diseases.

METHODS

Granulosa cell-related lncRNAs data spanning five follicular development stages were retrieved and filtered from the NCBI dataset (GSE107746). Stage-specific lncRNA expression patterns and mRNA-lncRNA co-expression networks were identified with bioinformatic approaches. Subsequently, the expression pattern of SNHG18 was detected in GCs during ovarian aging. And SNHG18 siRNA or overexpression plasmids were transfected to SVOG cells in examining the regulatory roles of SNHG18 in GC proliferation and apoptosis. Moreover, whether PKCɛ/SNHG18 signaling take part in GC glycolysis via ENO1 were verified in SVOG cells.

RESULTS

We demonstrated that GC-related lncRNAs were specifically expressed across different developmental stages, and coordinated crucial biological functions like mitotic cell cycle and metabolic processes in the folliculogenesis. Thereafter, we noticed a strong correlation of PRKCE and SNHG18 expression in our analysis. With downregulated SNHG18 of GCs identified in the context of ovarian aging, SNHG18 knockdown could further induce cell apoptosis, retard cell proliferation and exacerbate DNA damage in SVOG cell. Moreover, downregulated PKCɛ/SNHG18 pathway interrupted the SVOG cell glycolysis by lowering the ENO1 expression.

CONCLUSIONS

Altogether, our results revealed that folliculogenesis-related lncRNA SNHG18 participated in the pathogenesis of ovarian aging, which may provide novel biomarkers for ovarian function and new insights for the infertility treatment.

Single-cell and spatiotemporal profile of ovulation in the mouse ovary

Ovulation is a spatiotemporally coordinated process that involves several tightly controlled events, including oocyte meiotic maturation, cumulus expansion, follicle wall rupture and repair, and ovarian stroma remodeling. To date, no studies have detailed the precise window of ovulation at single-cell resolution. Here, we performed parallel single-cell RNA-seq and spatial transcriptomics on paired mouse ovaries across an ovulation time course to map the spatiotemporal profile of ovarian cell types. We show that major ovarian cell types exhibit time-dependent transcriptional states enriched for distinct functions and have specific localization profiles within the ovary. We also identified gene markers for ovulation-dependent cell states and validated these using orthogonal methods. Finally, we performed cell-cell interaction analyses to identify ligand-receptor pairs that may drive ovulation, revealing previously unappreciated interactions. Taken together, our data provides a rich and comprehensive resource of murine ovulation that can be mined for discovery by the scientific community.

Transcriptomic signature of luteinized cumulus cells of oocytes developing to live birth after women received intracytoplasmic sperm injection

OBJECTIVE

To compare the transcriptome of human cumulus cells (CCs) from oocytes with different outcomes (pregnancy yes/no, live birth [LB] yes/no), to identify noninvasive biomarkers for oocyte selection as well as new therapeutic targets to increase LB rates from assisted reproductive technologies (ART).

DESIGN

Retrospective observational study.

SETTINGS

This study was conducted at a University Hospital in Switzerland.

PATIENTS

Subfertile couples undergoing controlled ovarian superstimulation and intracytoplasmic sperm injection with subsequent unbiopsied embryo transfer below the female age of 43 years.

INTERVENTION(S)

RNA sequencing of CCs from oocytes results in a pregnancy, no pregnancy, LB, or no LB.

MAIN OUTCOME MEASURES

Differential gene expression (DEG) between CCs of oocytes results in "no pregnancy" vs. "pregnancy" and "pregnancy only" vs. "live birth."

RESULTS

Although RNA sequencing did not reveal DEGs when comparing the transcriptomic profiles of the groups "no pregnancy" with "pregnancy," we identified 139 DEGs by comparing "pregnancy only" with "live birth," of which 28 belonged to clusters relevant to successful ART outcomes (i.e., CTGF, SERPINE2, PCK1, HHIP, HS3ST, and BIRC5). A functional enrichment analysis revealed that the transcriptome of CCs associated with LB depicts pathways of extracellular matrix, inflammatory cascades leading to ovulation, cell patterning, proliferation, and differentiation, and silencing pathways leading to apoptosis.

CONCLUSION

We identified a CCs transcriptomic profile associated with LB after embryo transfer that, after further validation, could serve to predict successful ART outcomes. The definition of relevant pathways of CCs related to oocyte competency contributes to a broader understanding of the cumulus oocyte complex and helps identify further therapeutic targets for improving ART success.

Unraveling the complexity of polycystic ovary syndrome with animal models

Polycystic ovary syndrome (PCOS) is a highly familial and heritable endocrine disorder. Over half of the daughters born to women with PCOS may eventually develop their own PCOS-related symptoms. Progress in the treatment of PCOS is currently hindered by the complexity of its clinical manifestations and incomplete knowledge of its etiopathogenesis. Various animal models, including experimentally induced, naturally occurring, and spontaneously arising ones, have been established to emulate a wide range of phenotypical and pathological traits of human PCOS. These studies have led to a paradigm shift in understanding the genetic, developmental, and evolutionary origins of this disorder. Furthermore, emerging evidence suggests that animal models are useful in evaluating state-of-the-art drugs and treatments for PCOS. This review aims to provide a comprehensive summary of recent studies of PCOS in animal models, highlighting the power of these disease models in understanding the biology of PCOS and aiding high-throughput approaches.

A spatiotemporal molecular atlas of the ovulating mouse ovary

Ovulation is essential for reproductive success, yet the underlying cellular and molecular mechanisms are far from clear. Here, we applied high-resolution spatiotemporal transcriptomics to map out cell type- and ovulation stage-specific molecular programs as function of time during follicle maturation and ovulation in mice. Our analysis revealed dynamic molecular transitions within granulosa cell types that occur in tight coordination with mesenchymal cell proliferation. We identified molecular markers for the emerging cumulus cell fate during the preantral-to-antral transition. We describe transcriptional programs that respond rapidly to ovulation stimulation and those associated with follicle rupture, highlighting the prominent roles of apoptotic and metabolic pathways during the final stages of follicle maturation. We further report stage-specific oocyte-cumulus cell interactions and diverging molecular differentiation in follicles approaching ovulation. Collectively, this study provides insights into the cellular and molecular processes that regulate mouse ovarian follicle maturation and ovulation with important implications for advancing therapeutic strategies in reproductive medicine.

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.

Single-cell sequencing reveals the reproductive variations between primiparous and multiparous Hu ewes

BACKGROUND

In the modern sheep production systems, the reproductive performance of ewes determines the economic profitability of farming. Revealing the genetic mechanisms underlying differences in the litter size is important for the selection and breeding of highly prolific ewes. Hu sheep, a high-quality Chinese sheep breed, is known for its high fecundity and is often used as a model to study prolificacy traits. In the current study, animals were divided into two groups according to their delivery rates in three consecutive lambing seasons (namely, the high and low reproductive groups with ≥ 3 lambs and one lamb per season, n = 3, respectively). The ewes were slaughtered within 12 h of estrus, and unilateral ovarian tissues were collected and analyzed by 10× Genomics single-cell RNA sequencing.

RESULTS

A total of 5 types of somatic cells were identified and corresponding expression profiles were mapped in the ovaries of each group. Noticeably, the differences in the ovary somatic cell expression profiles between the high and low reproductive groups were mainly clustered in the granulosa cells. Furthermore, four granulosa cell subtypes were identified. GeneSwitches analysis revealed that the abundance of JPH1 expression and the reduction of LOC101112291 expression could lead to different evolutionary directions of the granulosa cells. Additionally, the expression levels of FTH1 and FTL in mural granulosa cells of the highly reproductive group were significantly higher. These genes inhibit necroptosis and ferroptosis of mural granulosa cells, which helps prevent follicular atresia.

CONCLUSIONS

This study provides insights into the molecular mechanisms underlying the high fecundity of Hu sheep. The differences in gene expression profiles, particularly in the granulosa cells, suggest that these cells play a critical role in female prolificacy. The findings also highlight the importance of genes such as JPH1, LOC101112291, FTH1, and FTL in regulating granulosa cell function and follicular development.

High-quality single-cell transcriptomics from ovarian histological sections during folliculogenesis

High-quality, straightforward single-cell RNA sequencing (RNA-seq) with spatial resolution remains challenging. Here, we developed DRaqL (direct RNA recovery and quenching for laser capture microdissection), an experimental approach for efficient cell lysis of tissue sections, directly applicable to cDNA amplification. Single-cell RNA-seq combined with DRaqL allowed transcriptomic profiling from alcohol-fixed sections with efficiency comparable with that of profiling from freshly dissociated cells, together with effective exon-exon junction profiling. The combination of DRaqL with protease treatment enabled robust and efficient single-cell transcriptome analysis from formalin-fixed tissue sections. Applying this method to mouse ovarian sections, we were able to predict the transcriptome of oocytes by their size and identified an anomaly in the size-transcriptome relationship relevant to growth retardation of oocytes, in addition to detecting oocyte-specific splice isoforms. Furthermore, we identified differentially expressed genes in granulosa cells in association with their proximity to the oocytes, suggesting distinct epigenetic regulations and cell-cycle activities governing the germ-soma relationship. Thus, DRaqL is a versatile, efficient approach for high-quality single-cell RNA-seq from tissue sections, thereby revealing histological heterogeneity in folliculogenic transcriptome.

FSHR-mTOR-HIF1 signaling alleviates mouse follicles from AMPK-induced atresia

The majority of activated ovarian follicles undergo atresia during reproductive life in mammals, and only a small number of follicles are ovulated. Though hormone treatment has been widely used to promote folliculogenesis, the molecular mechanism behind follicle selection and atresia remains under debate due to inconsistency among investigation models. Using a high-throughput molecular pathology strategy, we depicted a transcriptional atlas of mouse follicular granulosa cells (GCs) under physiological condition and obtained molecular signatures in healthy and atresia GCs during development. Functional results revealed hypoxia-inducible factor 1 (HIF1) as a major effector downstream of follicle-stimulating hormone (FSH), and HIF1 activation is essential for follicle growth. Energy shortage leads to prevalent AMP-activated protein kinase (AMPK) activation and drives follicular atresia. FSHR-mTOR-HIF1 signaling helps follicles escape from the atresia fate, while energy stress persists. Our work provides a comprehensive understanding of the molecular network behind follicle selection and atresia under physiological condition.

Maternal RNA binding protein with multiple splicing 2 (RBPMS2) is involved in mouse blastocyst formation through the bone morphogenetic protein pathway

RESEARCH QUESTION

Is early embryo development in mice influenced by RNA binding protein with multiple splicing 2 (RBPMS2), a maternal factor that accumulates and is stored in the cytoplasm of mature oocytes?

DESIGN

The expression patterns of RBPMS2 in mouse were analysed using quantitative real-time PCR (qRT PCR) and immunofluorescence staining. The effect of knockdown of RBPMS2 on embryo development was evaluated through a microinjection of specific morpholino or small interfering RNA. RNA sequencing was performed for mechanistic analysis. The interaction between RBPMS2 and the bone morphogenetic protein (BMP) pathway was studied using BMP inhibitor and activator. The effect on the localization of E-cadherin was determined by immunofluorescence staining.

RESULTS

Maternal protein RBPMS2 is highly expressed in mouse oocytes, and knockdown of RBPMS2 inhibits embryo development from the morula to the blastocyst stage. Mechanistically, RNA sequencing showed that the differentially expressed genes were enriched in the transforming growth factor-β (TGF-β) signalling pathway. BMPs are members of the TGF-β superfamily of growth factors. It was found that the addition of BMP inhibitor to the culture medium led to a morula-stage arrest, similar to that seen in RBPMS2 knockdown embryos. This morula-stage arrest defect caused by RBPMS2 knockdown was partially rescued by BMP activator. Furthermore, the localization of E-cadherin to the membrane was impaired in response to a knockdown of RBPMS2 or inhibition of the BMP pathway.

CONCLUSION

This study suggests that RBPMS2 activates the BMP pathway and thus influences the localization of E-cadherin, which is important for early mouse embryo development during blastocyst formation.

Deciphering Cellular Heterogeneity and Communication Patterns in Porcine Antral Follicles by Single-Cell RNA Sequencing

The antral follicle stage is a critical period in mammalian oocyte maturation, marked by complex interactions between oocyte development and neighboring granulosa cells. Understanding the heterogeneity and communication patterns of granulosa cells within antral follicles is crucial for deciphering their roles in follicle development and oocyte maturation. Here, we employed single-cell RNA-sequencing to explore the molecular and cellular characteristics of porcine antral follicles. Our analysis revealed distinct subpopulations within mural and cumulus granulosa cells, indicating diverse cellular states and functions within the follicles. Functional enrichment analysis unveiled the involvement of specific subpopulations in steroid biosynthesis, cumulus expansion, and cellular communication. Moreover, comparing mature and less mature follicles highlighted differences in cell distribution and functions, indicating developmental-specific variations. Our findings shed light on the intricate cellular heterogeneity and communication network within porcine antral follicles, providing valuable insights into the regulation of follicle development and oocyte maturation in pigs. These results hold promise for improving pig reproductive efficiency and advancing human reproductive medicine.

A single-cell atlas of the aging murine ovary

Ovarian aging leads to diminished fertility, dysregulated endocrine signaling, and increased chronic disease burden. These effects begin to emerge long before follicular exhaustion. Around 35 years old, women experience a sharp decline in fertility, corresponding to declines in oocyte quality. Despite a growing body of work, the field lacks a comprehensive cellular map of the transcriptomic changes in the aging ovary to identify early drivers of ovarian decline. To fill this gap, we performed single-cell RNA sequencing on ovarian tissue from young (3-month-old) and reproductively aged (9-month-old) mice. Our analysis revealed a doubling of immune cells in the aged ovary, with lymphocyte proportions increasing the most, which was confirmed by flow cytometry. We also found an age-related downregulation of collagenase pathways in stromal fibroblasts, which corresponds to rises in ovarian fibrosis. Follicular cells displayed stress response, immunogenic, and fibrotic signaling pathway inductions with aging. This report raises provides critical insights into mechanisms responsible for ovarian aging phenotypes.

A Spatiotemporal Molecular Atlas of the Ovulating Mouse Ovary

Ovulation is essential for reproductive success, yet the underlying cellular and molecular mechanisms are far from clear. Here, we applied high-resolution spatiotemporal transcriptomics to map out cell-type- and ovulation-stage-specific molecular programs as function of time during follicle maturation and ovulation in mice. Our analysis revealed dynamic molecular transitions within granulosa cell types that occur in tight coordination with mesenchymal cell proliferation. We identified new molecular markers for the emerging cumulus cell fate during the preantral-to-antral transition. We describe transcriptional programs that respond rapidly to ovulation stimulation and those associated with follicle rupture, highlighting the prominent roles of apoptotic and metabolic pathways during the final stages of follicle maturation. We further report stage-specific oocyte-cumulus cell interactions and diverging molecular differentiation in follicles approaching ovulation. Collectively, this study provides insights into the cellular and molecular processes that regulate mouse ovarian follicle maturation and ovulation with important implications for advancing therapeutic strategies in reproductive medicine.

Advances in Oocyte Maturation In Vivo and In Vitro in Mammals

The quality and maturation of an oocyte not only play decisive roles in fertilization and embryo success, but also have long-term impacts on the later growth and development of the fetus. Female fertility declines with age, reflecting a decline in oocyte quantity. However, the meiosis of oocytes involves a complex and orderly regulatory process whose mechanisms have not yet been fully elucidated. This review therefore mainly focuses on the regulation mechanism of oocyte maturation, including folliculogenesis, oogenesis, and the interactions between granulosa cells and oocytes, plus in vitro technology and nuclear/cytoplasm maturation in oocytes. Additionally, we have reviewed advances made in the single-cell mRNA sequencing technology related to oocyte maturation in order to improve our understanding of the mechanism of oocyte maturation and to provide a theoretical basis for subsequent research into oocyte maturation.

RNA sequencing and expression analysis reveal a role for Lhx9 in the haploinsufficient adult mouse ovary

Understanding the molecular pathways that underpin ovarian development and function is vital for improving the research approaches to investigating fertility. Despite a significant improvement in our knowledge of molecular activity in the ovary, many questions remain unanswered in the quest to understand factors influencing fertility and ovarian pathologies such as cancer. Here, we present an investigation into the expression and function of the developmental transcription factor LIM Homeobox 9 (LHX9) in the adult mouse ovary. We have characterized Lhx9 expression in several cell types of the mature ovary across follicle stages. To evaluate possible LHX9 function in the adult ovary, we investigated ovarian anatomy and transcription in an Lhx9^+/- knockout mouse model displaying subfertility. Despite a lack of gross anatomical differences between genotypes, RNA-sequencing found that 90 differentially expressed genes between Lhx9^{+/ -} and Lhx9^+/+ mice. Gene ontology analyses revealed a reduced expression of genes with major roles in ovarian steroidogenesis and an increased expression of genes associated with ovarian cancer. Analysis of the ovarian epithelium revealed Lhx9^{+/ -} mice have a disorganized epithelial phenotype, corresponding to a significant increase in epithelial marker gene expression. These results provide an analysis of Lhx9 in the adult mouse ovary, suggesting a role in fertility and ovarian epithelial cancer.

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.

Oocyte quality following in vitro follicle development†

In vitro follicle development (IVFD) is an adequate model to obtain basic knowledge of folliculogenesis and provides a tool for ovarian toxicity screening. IVFD yielding competent oocytes may also offer an option for fertility and species preservation. To promote follicle growth and oocyte maturation in vitro, various culture systems are utilized for IVFD in rodents, domestic animals, wild animals, nonhuman primates, and humans. Follicle culture conditions have been improved by optimizing gonadotropin levels, regulatory factors, nutrient supplements, oxygen concentration, and culture matrices. This review summarizes quality assessment of oocytes generated from in vitro-developed antral follicles from the preantral stage, including oocyte epigenetic and genetic profile, cytoplasmic and nuclear maturation, preimplantation embryonic development following in vitro fertilization, as well as pregnancy and live offspring after embryo transfer. The limitations of oocyte quality evaluation following IVFD and the gaps in our knowledge of IVFD to support proper oocyte development are also discussed. The information may advance our understanding of the requirements for IVFD, with a goal of producing competent oocytes with genetic integrity to sustain embryonic development resulting in healthy offspring.