Single-cell transcriptome landscape of ovarian cells during primordial follicle assembly in mice

Inter- and trans-generational impacts of environmental exposures on the germline resolved at the single-cell level

Reproduction is a remarkably intricate process involving the interaction of multiple cell types and organ systems unfolding over long periods of time and that culminates with the production of gametes. The initiation of germ cell development takes place during embryogenesis but only completes decades later in humans. The complexity inherent to reproduction and its study has long hampered our ability to decipher how environmental agents disrupt this process. Single-cell approaches provide an opportunity for a deeper understanding of the action of toxicants on germline function and analyze how the response to their exposure is differentially distributed across tissues and cell types. In addition to single-cell RNA sequencing, other single-cell or nucleus level approaches such as ATAC-sequencing and multi-omics have expanded the strategies that can be implemented in reproductive toxicological studies to include epigenomic and the nuclear transcriptomic data. Here we will discuss the current state of single-cell technologies and how they can best be utilized to advance reproductive toxicological studies. We will then discuss case studies in two model organisms (Caenorhabditis elegans and mouse) studying different environmental exposures (alcohol and e-cigarettes respectively) to highlight the value of single-cell and single-nucleus approaches for reproductive biology and reproductive toxicology.

Endometrial stem cells alleviate cisplatin-induced ferroptosis of granulosa cells by regulating Nrf2 expression

BACKGROUND

Premature ovarian failure (POF) caused by cisplatin is a severe and intractable sequela for young women with cancer who received chemotherapy. Cisplatin causes the dysfunction of granulosa cells and mainly leads to but is not limited to its apoptosis and autophagy. Ferroptosis has been also reported to participate, while little is known about it. Our previous experiment has demonstrated that endometrial stem cells (EnSCs) can repair cisplatin-injured granulosa cells. However, it is still unclear whether EnSCs can play a repair role by acting on ferroptosis.

METHODS

Western blotting and quantitative reverse-transcription polymerase chain reaction (qRT-PCR) were applied to detect the expression levels of ferroptosis-related genes. CCK-8 and 5-Ethynyl-2'-deoxyuridine (EdU) assays were used to evaluate cell viability. Transmission electron microscopy (TEM) was performed to detect ferroptosis in morphology. And the extent of ferroptosis was assessed by ROS, GPx, GSSG and MDA indicators. In vivo, ovarian morphology was presented by HE staining and the protein expression in ovarian tissue was detected by immunohistochemistry.

RESULTS

Our results showed that ferroptosis could occur in cisplatin-injured granulosa cells. Ferroptosis inhibitor ferrostatin-1 (Fer-1) and EnSCs partly restored cell viability and mitigated the damage of cisplatin to granulosa cells by inhibiting ferroptosis. Moreover, the repair potential of EnSCs can be markedly blocked by ML385.

CONCLUSION

Our study demonstrated that cisplatin could induce ferroptosis in granulosa cells, while EnSCs could inhibit ferroptosis and thus exert repair effects on the cisplatin-induced injury model both in vivo and in vitro. Meanwhile, Nrf2 was validated to participate in this regulatory process and played an essential role.

Single-cell profiling reveals transcriptome dynamics during bovine oocyte growth

BACKGROUND

Mammalian follicle development is characterized by extensive changes in morphology, endocrine responsiveness, and function, providing the optimum environment for oocyte growth, development, and resumption of meiosis. In cattle, the first signs of transcription activation in the oocyte are observed in the secondary follicle, later than during mouse and human oogenesis. While many studies have generated extensive datasets characterizing gene expression in bovine oocytes, they are mostly limited to the analysis of fully grown and matured oocytes. The aim of the present study was to apply single-cell RNA sequencing to interrogate the transcriptome of the growing bovine oocyte from the secondary follicle stage through to the mid-antral follicle stage.

RESULTS

Single-cell RNA-seq libraries were generated from oocytes of known diameters (< 60 to > 120 μm), and datasets were binned into non-overlapping size groups for downstream analysis. Combining the results of weighted gene co-expression network and Trendy analyses, and differently expressed genes (DEGs) between size groups, we identified a decrease in oxidative phosphorylation and an increase in maternal -genes and transcription regulators across the bovine oocyte growth phase. In addition, around 5,000 genes did not change in expression, revealing a cohort of stable genes. An interesting switch in gene expression profile was noted in oocytes greater than 100 μm in diameter, when the expression of genes related to cytoplasmic activities was replaced by genes related to nuclear activities (e.g., chromosome segregation). The highest number of DEGs were detected in the comparison of oocytes 100-109 versus 110-119 μm in diameter, revealing a profound change in the molecular profile of oocytes at the end of their growth phase.

CONCLUSIONS

The current study provides a unique dataset of the key genes and pathways characteristic of each stage of oocyte development, contributing an important resource for a greater understanding of bovine oogenesis.

The SLC19A1-AS/miR-1343/WNT11 axis is a novel positive regulatory ceRNA network governing goat granulosa cell proliferation

Long noncoding RNAs (lncRNAs), as competitive endogenous RNAs (ceRNAs), can directly or indirectly affect the proliferation and apoptosis of granulosa cells by regulating microRNA (miRNA) pathways. A ceRNA network of the SLC19A1-AS-miR-1343-WNT11 axis was constructed via comprehensive transcriptome sequencing of ovaries from goats with various fertility levels to further elucidate the function and regulatory mechanism of SLC19A1-AS in modulating miR-1343 and WNT11 during granulosa cell proliferation and apoptosis. Subsequent validation experiments were conducted in vitro using granulosa cells. In these experiments, we performed RNA immunoprecipitation (RIP) and identified SLC19A1-AS as a ceRNA in goat granulosa cells that promoted proliferation. Through bioinformatics prediction, luciferase reporter gene assays, and RNA pulldown assays, we confirmed that SLC19A1-AS acts as a sponge for miR-1343, preventing its binding to WNT11 mRNA and thereby increasing the expression of WNT11. This interaction also influenced the proliferation and apoptosis of granulosa cells. Our study systematically validated the biological function of the lncRNA-miRNA-mRNA ceRNA network in goat ovaries and revealed the potential regulatory mechanism by which SLC19A1-AS functions as a ceRNA in granulosa cells. These findings are expected to provide an important experimental foundation for further elucidating the physiological regulatory network of the ovary and contributing to reproductive health in goats.

A matter of new life and cell death: programmed cell death in the mammalian ovary

BACKGROUND

The mammalian ovary is a unique organ that displays a distinctive feature of cyclic changes throughout the entire reproductive period. The estrous/menstrual cycles are associated with drastic functional and morphological rearrangements of ovarian tissue, including follicular development and degeneration, and the formation and subsequent atrophy of the corpus luteum. The flawless execution of these reiterative processes is impossible without the involvement of programmed cell death (PCD).

MAIN TEXT

PCD is crucial for efficient and careful clearance of excessive, depleted, or obsolete ovarian structures for ovarian cycling. Moreover, PCD facilitates selection of high-quality oocytes and formation of the ovarian reserve during embryonic and juvenile development. Disruption of PCD regulation can heavily impact the ovarian functions and is associated with various pathologies, from a moderate decrease in fertility to severe hormonal disturbance, complete loss of reproductive function, and tumorigenesis. This comprehensive review aims to provide updated information on the role of PCD in various processes occurring in normal and pathologic ovaries. Three major events of PCD in the ovary-progenitor germ cell depletion, follicular atresia, and corpus luteum degradation-are described, alongside the detailed information on molecular regulation of these processes, highlighting the contribution of apoptosis, autophagy, necroptosis, and ferroptosis. Ultimately, the current knowledge of PCD aberrations associated with pathologies, such as polycystic ovarian syndrome, premature ovarian insufficiency, and tumors of ovarian origin, is outlined.

CONCLUSION

PCD is an essential element in ovarian development, functions and pathologies. A thorough understanding of molecular mechanisms regulating PCD events is required for future advances in the diagnosis and management of various disorders of the ovary and the female reproductive system in general.

Melatonin alleviates the toxic effect of di(2-ethylhexyl) phthalate on oocyte quality resulting from CEBPB suppression during primordial follicle formation

Presently, the exposure of plasticizers to humans and animals occurs daily, which pose a potential threat to reproductive health. In the present study, a pregnant mouse model exposed to di(2-ethylhexyl) phthalate (DEHP, one of the most common plasticizers) and melatonin was established, and the single-cell transcriptome technology was applied to investigate the effects of melatonin in ovarian cells against DEHP. Results showed that DEHP markedly altered the gene expression pattern of ovarian cells, and severely weakened the histone methylation modification of oocytes. The administration of melatonin recovered the expression of LHX8 and SOHLH1 proteins that essential for primordial follicle formation, and increased the expression of CEBPB, as well as key genes of histone methylation modification (such as Smyd3 and Kdm5a). In addition, the ovarian damage caused by DEHP was also relieved after the overexpression of CEBPB, which suggested melatonin could improve primordial follicle formation progress via enhancing CEBPB expression in mice. Besides, the apoptosis of ovarian cells induced by DEHP also was diminished by melatonin. The study provides evidence of melatonin preventing the damage mediated by plasticizers on the reproductive system in females and CEBPB may serve as a downstream target factor of melatonin in the process.

Mapping Cell Atlases at the Single-Cell Level

Recent advancements in single-cell technologies have led to rapid developments in the construction of cell atlases. These atlases have the potential to provide detailed information about every cell type in different organisms, enabling the characterization of cellular diversity at the single-cell level. Global efforts in developing comprehensive cell atlases have profound implications for both basic research and clinical applications. This review provides a broad overview of the cellular diversity and dynamics across various biological systems. In addition, the incorporation of machine learning techniques into cell atlas analyses opens up exciting prospects for the field of integrative biology.

Data Analysis Pipeline for scRNA-seq Experiments to Study Early Oogenesis

Germ cells as the means for the transmission of genetic information between generations have been a hot topic of research for decades. The analysis of the transcriptomes, that is of the RNA transcripts produced by the genotype at a given time, of germ cells and the surrounding somatic cells, is essential to unravel the cellular and molecular processes regulating gametogenesis. However, the asynchronized differentiation of germ cells and high cellular heterogeneity in the developing ovary or testis represent two unsurmountable challenges for delineating the transcription regulation mechanism of germ cells using traditional bulk RNA sequencing. By performing single-cell RNA sequencing (scRNA-seq), it is now possible to dissect the transcriptome of germ cell development at single-cell resolution, and apply powerful bioinformatics methods to translate raw sequencing data into meaningful information. Here, using the 10× Genomic platform and the most widely cited bioinformatics tools, we describe how to analyze early female germ cell development using scRNA-seq data generated from mouse E11.5 to E14.5 ovaries. This pipeline will provide a guide for exploring the processes of early germ cell development at single-cell resolution.

NADase CD38 is a key determinant of ovarian aging

The ovary ages earlier than most other tissues, yet the underlying mechanisms remain elusive. Here a comprehensive analysis of transcriptomic landscapes in different organs in young and middle-aged mice revealed that the ovaries showed earlier expression of age-associated genes, identifying increased NADase CD38 expression and decreased NAD+ levels in the ovary of middle-aged mice. Bulk and single-cell RNA sequencing revealed that CD38 deletion mitigated ovarian aging, preserving fertility and follicle reserve in aged mice by countering age-related gene expression changes and intercellular communication alterations. Mechanistically, the earlier onset of inflammation induced higher expression levels of CD38 and decreased NAD+ levels in the ovary, thereby accelerating ovarian aging. Consistently, pharmacological inhibition of CD38 enhanced fertility in middle-aged mice. Our findings revealed the mechanisms underlying the earlier aging of the ovary relative to other organs, providing a potential therapeutic target for ameliorating age-related female infertility.

A single-cell atlas of the aging mouse ovary

Ovarian aging leads to diminished fertility, dysregulated endocrine signaling and increased chronic disease burden. These effects begin to emerge long before follicular exhaustion. Female humans experience a sharp decline in fertility around 35 years of age, which corresponds 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 mouse 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 provides critical insights into mechanisms responsible for ovarian aging phenotypes. The data can be explored interactively via a Shiny-based web application.

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.

Single-cell RNA-seq identified novel genes involved in primordial follicle formation

Introduction

The number of primordial follicles (PFs) in mammals determines the ovarian reserve, and impairment of primordial follicle formation (PFF) will cause premature ovarian insufficiency (POI).

Methods

By analyzing public single-cell RNA sequencing performed during PFF on mice and human ovaries, we identified novel functional genes and novel ligand-receptor interaction during PFF. Based on immunofluorescence and in vitro ovarian culture, we confirmed mechanisms of genes and ligand-receptor interaction in PFF. We also applied whole exome sequencing (WES) in 93 cases with POI and whole genome sequencing (WGS) in 465 controls. Variants in POI patients were further investigated by in silico analysis and functional verification.

Results

We revealed ANXA7 (annexin A7) and GTF2F1 (general transcription factor IIF subunit 1) in germ cells to be novel potentially genes in promoting PFF. Ligand Mdk (midkine) in germ cells and its receptor Sdc1 (syndecan 1) in granulosa cells are novel interaction crucial for PFF. Based on immunofluorescence, we confirmed significant up-regulation of ANXA7 in PFs compared with germline cysts, and uniform expression of GTF2F1, MDK and SDC1 during PFF, in 25 weeks human fetal ovary. In vitro investigation indicated that Anxa7 and Gtf2f1 are vital for mice PFF by regulating Jak/Stat3 and Jnk signaling pathways, respectively. Ligand-receptor (Mdk-Sdc1) are crucial for PFF by regulating Pi3k-akt signaling pathway. Two heterozygous variants in GTF2F1, and one heterozygous variants in SDC1 were identified in cases, but no variant were identified in controls. The protein level of GTF2F1 or SDC1 in POI cases are significantly lower than that of controls, indicating the pathogenic effects of the two genes on ovarian function were dosage dependent.

Discussion

Our study identified novel genes and novel ligand-receptor interaction during PFF, and further expanding the genetic architecture of POI.

A spatiotemporal gene expression and cell atlases of the developing rat ovary

Normal ovarian development is necessary for the production of healthy oocytes. However, the characteristics of oocytes development at different stages and the regulatory relationship between oocytes and somatic cells remain to be fully explained. Here, we combined scRNA-seq and spatial transcriptomic sequencing to profile the transcriptomic atlas of developing ovarian of the rat. We identified four components from developing granulosa cells including cumulus, primitive, mural, and luteal cells, and constructed their differential transcriptional regulatory networks. Several novel growth signals from oocytes to cumulus cells were identified, such as JAG1-NOTCH2 and FGF9-FGFR2. Moreover, we observed three cumulus sequential phases during follicle development determined by the key transcriptional factors in each cumulus phase (Bckaf1, Gata6, Cebpb, etc.), as well as the potential pinpointed roles of macrophages in luteal regression. Altogether, the single-cell spatial transcriptomic profile of the ovary provides not only a new research dimension for temporal and spatial analysis of ovary development, but also valuable data resources and a research basis for in-depth excavation of the mechanisms of mammalian ovary development.

Maternal vitamin B1 is a determinant for the fate of primordial follicle formation in offspring

The mediation of maternal-embryonic cross-talk via nutrition and metabolism impacts greatly on offspring health. However, the underlying key interfaces remain elusive. Here, we determined that maternal high-fat diet during pregnancy in mice impaired preservation of the ovarian primordial follicle pool in female offspring, which was concomitant with mitochondrial dysfunction of germ cells. Furthermore, this occurred through a reduction in maternal gut microbiota-related vitamin B1 while the defects were restored via vitamin B1 supplementation. Intriguingly, vitamin B1 promoted acetyl-CoA metabolism in offspring ovaries, contributing to histone acetylation and chromatin accessibility at the promoters of cell cycle-related genes, enhancement of mitochondrial function, and improvement of granulosa cell proliferation. In humans, vitamin B1 is downregulated in the serum of women with gestational diabetes mellitus. In this work, these findings uncover the role of the non-gamete transmission of maternal high-fat diet in influencing offspring oogenic fate. Vitamin B1 could be a promising therapeutic approach for protecting offspring health.

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.

Cyanidin-3-O-glucoside Mitigates the Ovarian Defect Induced by Zearalenone via p53-GADD45a Signaling during Primordial Follicle Assembly

Zearalenone (ZEN) is well known as a kind of endocrine disruptor whose exposure is capable of causing reproductive toxicity in animals. Cyanidin-3-O-glucoside (C3G) is a derivative of cyanidin and owns multiple biofunctions, and prior efforts have suggested that C3G has therapeutic actions for reproductive diseases. In this article, a ZEN exposure model during primordial follicle assembly was constructed using the in vitro culture platform of neonatal mouse ovaries. We investigated the protective effect of C3G on ZEN-induced ovarian toxicity during primordial follicle assembly in mice, as well as its potential mechanism. Interestingly, we observed that C3G could effectively protect the ovary from ZEN damage, mainly by restoring primordial follicle assembly, which upregulated the expression of LHX8 and SOHLH1 proteins and relieved ZEN-induced DNA damage. Next, to explore the mechanism by which C3G rescued ZEN-induced injury, we performed RNA sequencing (RNA-seq). The bioinformatic analysis illustrated that the rescue pathway of C3G was associated with p53-Gadd45a signaling and cell cycle. Then, western blotting and flow cytometry results revealed that C3G restored the expression levels of cyclin-dependent kinase 6 (CDK6) and cyclin D2 (CCND2) and regulated the ovarian cell cycle to normal. In conclusion, our findings manifested that C3G could alleviate ZEN-induced primordial follicle assembly impairment by restoring the cell cycle involved in p53-GADD45a signaling.

Cross-species analysis of transcriptome emphasizes a critical role of TNF-α in mediating MAP2K7/AKT2 signaling in zearalenone-induced apoptosis

Zearalenone (ZEN) is a widespread and transgenerational toxicant that can cause serious reproductive health risks, which poses a potential threat to global agricultural production and human health; its estrogenic activity can lead to reproductive toxicity through the induction of granulosa cell apoptosis. Herein, comparative transcriptome analysis, single-cell transcriptome analysis, and weighted gene co-expression network analysis (WGCNA) combined with gene knockout in vivo and RNA interference in vitro were used to comprehensively describe the damage caused by ZEN exposure on ovarian granulosa cells. Comparative transcriptome analysis and WGCNA suggested that the tumor necrosis factor (TNF)-α-mediated mitogen-activated protein kinase 7 (MAP2K7)/ AKT serine/threonine kinase 2 (AKT2) axis was disordered after ZEN exposure in porcine granulosa cells (pGCs) and mouse granulosa cells (mGCs). In vivo gene knockout and in vitro RNA interference verified that TNF-α-mediated MAP2K7/AKT2 was the guiding signal in ZEN-induced apoptosis in pGCs and mGCs. Moreover, single-cell transcriptome analysis showed that ZEN exposure could induce changes in the TNF signaling pathway in offspring. Overall, we concluded that the TNF-α-mediated MAP2K7/AKT2 axis was the main signaling pathway of ZEN-induced apoptosis in pGCs and mGCs. This work provides new insights into the mechanism of ZEN toxicity and provides new potential therapeutic targets for the loss of livestock and human reproductive health caused by ZEN.

Epigenetic programming in the ovarian reserve

The ovarian reserve defines female reproductive lifespan, which in humans spans decades. The ovarian reserve consists of oocytes residing in primordial follicles arrested in meiotic prophase I and is maintained independent of DNA replication and cell proliferation, thereby lacking stem cell-based maintenance. Largely unknown is how cellular states of the ovarian reserve are established and maintained for decades. Our recent study revealed that a distinct chromatin state is established during ovarian reserve formation in mice, uncovering a novel window of epigenetic programming in female germline development. We showed that an epigenetic regulator, Polycomb Repressive Complex 1 (PRC1), establishes a repressive chromatin state in perinatal mouse oocytes that is essential for prophase I-arrested oocytes to form the ovarian reserve. Here we discuss the biological roles and mechanisms underlying epigenetic programming in ovarian reserve formation, highlighting current knowledge gaps and emerging research areas in female reproductive biology.

The transcriptome reveals the molecular regulatory network of primordial follicle depletion in obese mice

OBJECTIVE

To explore the dynamic transcriptional regulatory network of primordial follicle fate in obese mice to elucidate the potential mechanism of primordial follicle depletion.

DESIGN

Experimental study and transcriptomic analysis.

ANIMALS

Healthy (n=15) and obese (n=15) female mice.

INTERVENTIONS

Six-week-old CD-1 mice were divided into healthy and high-fat diet groups and fed continuously for 12 weeks. The diet of healthy mice contained 10% fat. The diet of high-fat mice contained 60% fat.

MAIN OUTCOME MEASURES

Primordial to primary follicle transition rate, gene expression changes, enriched Kyoto Encyclopedia of Genes and Genomes pathways, and ferroptosis.

RESULTS

Primordial follicle depletion was increased in the ovaries of obese mice. We found that deposited fat around primordial and primary follicles of obese mice was higher than that for healthy mice. The proliferation of granulosa cells around primary follicles was increased in obese mice. In addition, we uncovered specific gene signatures associated with the primordial to primary follicle transition (PPT) in obese mice using laser capture microdissection RNA sequencing analysis. Gene set enrichment analysis indicated that ferroptosis, cell oxidation, vascular endothelial growth factor, and mammalian target of rapamycin signaling were increased significantly in the primordial follicles of obese mice. Notably, the ferritin, acyl CoA synthetase long-chain family member 4, and solid carrier family 7 member 11 associated proteins of the ferroptosis signaling pathway were significantly increased in the PPT phase of obese mice.

CONCLUSION

Our work suggests that ferroptosis is a key pathway activated within immature ovarian follicles in the context of obesity and that the process may be involved in the physiological regulation of the PPT as well.

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.

LSM14B controls oocyte mRNA storage and stability to ensure female fertility

Controlled mRNA storage and stability is essential for oocyte meiosis and early embryonic development. However, how to regulate mRNA storage and stability in mammalian oogenesis remains elusive. Here we showed that LSM14B, a component of membraneless compartments including P-body-like granules and mitochondria-associated ribonucleoprotein domain (MARDO) in germ cell, is indispensable for female fertility. To reveal loss of LSM14B disrupted primordial follicle assembly and caused mRNA reduction in non-growing oocytes, which was concomitant with the impaired assembly of P-body-like granules. 10× Genomics single-cell RNA-sequencing and immunostaining were performed. Meanwhile, we conducted RNA-seq analysis of GV-stage oocytes and found that Lsm14b deficiency not only impaired the maternal mRNA accumulation but also disrupted the translation in fully grown oocytes, which was closely associated with dissolution of MARDO components. Moreover, Lsm14b-deficient oocytes reassembled a pronucleus containing decondensed chromatin after extrusion of the first polar body, through compromising the activation of maturation promoting factor, while the defects were restored via WEE1/2 inhibitor. Together, our findings reveal that Lsm14b plays a pivotal role in mammalian oogenesis by specifically controlling of oocyte mRNA storage and stability.

Purinergic receptor P2X7 activates NOX2/JNK signaling to participate in granulosa cell inflammation and apoptosis in polycystic ovary syndrome

Increasing evidence shows that polycystic ovary syndrome (PCOS) is often accompanied by an inflammatory response, hence, appropriately managing granulosa cell inflammation is critical to regaining ovarian function in PCOS. In this study, the differential levels of purinergic receptor P2X7 between the control and PCOS samples in the dataset GSE34526 were assessed, then PCOS mouse models were established. Following evaluating the fluctuations in hormone levels, inflammatory cytokines, and P2X7, mice received treatment with the P2X7 antagonist A740003. Its effects on hormones, inflammation, apoptosis, and NOX2 signaling in mice were examined. Afterward, primary mouse granulosa cells were isolated, and the mediating role of NOX2 signaling in the P2X7 regulatory pathway was confirmed by transfection of NOX2 overexpression plasmids. The results demonstrated that P2X7 was significantly elevated in the PCOS samples in the dataset. Compared with the control group, PCOS mice had significant differences in the follicle-stimulating hormone, luteinizing hormone, testosterone, anti-Müllerian hormone, inflammatory factors, and P2X7. Treatment with A740003 partially restored these parameter levels, including NOX2 signaling. Based on in vitro experiments on primary mouse granulosa cells, the above findings were re-verified, and the overexpression of NOX2 could reverse the regulatory function of P2X7. The present study highlights that P2X7 level increases in PCOS, and inhibition of P2X7 can reduce disease symptoms. It is involved in inflammation and apoptosis in granulosa cells through NOX2/JNK signaling.

Current progress on in vitro differentiation of ovarian follicles from pluripotent stem cells

Mammalian female reproduction requires a functional ovary. Competence of the ovary is determined by the quality of its basic unit-ovarian follicles. A normal follicle consists of an oocyte enclosed within ovarian follicular cells. In humans and mice, the ovarian follicles are formed at the foetal and the early neonatal stage respectively, and their renewal at the adult stage is controversial. Extensive research emerges recently to produce ovarian follicles in-vitro from different species. Previous reports demonstrated the differentiation of mouse and human pluripotent stem cells into germline cells, termed primordial germ cell-like cells (PGCLCs). The germ cell-specific gene expressions and epigenetic features including global DNA demethylation and histone modifications of the pluripotent stem cells-derived PGCLCs were extensively characterized. The PGCLCs hold potential for forming ovarian follicles or organoids upon cocultured with ovarian somatic cells. Intriguingly, the oocytes isolated from the organoids could be fertilized in-vitro. Based on the knowledge of in-vivo derived pre-granulosa cells, the generation of these cells from pluripotent stem cells termed foetal ovarian somatic cell-like cells was also reported recently. Despite successful in-vitro folliculogenesis from pluripotent stem cells, the efficiency remains low, mainly due to the lack of information on the interaction between PGCLCs and pre-granulosa cells. The establishment of in-vitro pluripotent stem cell-based models paves the way for understanding the critical signalling pathways and molecules during folliculogenesis. This article aims to review the developmental events during in-vivo follicular development and discuss the current progress of generation of PGCLCs, pre-granulosa and theca cells in-vitro.

Single-cell transcriptomics allows novel insights into the endocrine-disrupting chemicals induced mammalian reproductive disorder

Increasing environmental pollution has led to the spread of many endocrine-disrupting chemicals (EDCs) around the world, which are toxic substances in the form of compounds that pose a great threat to the reproductive health of mammals and become a potential cause of many reproductive function-related diseases. In the past decade, the rapid development of single-cell RNA sequencing (scRNA-seq) technology has greatly promoted the study of the toxic mechanisms of EDCs in the mammalian reproductive system, including DEHP, ZEN, BPA, and BDE47. These studies aim to resolve the interference of EDCs in critical stages of reproductive development, including prepubertal and pubertal in males, meiosis I and early follicle formation in females. This paper introduces the sequencing process and analysis methods of current mainstream scRNA-seq technology, systematically reviews the outstanding contributions and specific research ideas of this technology in the study of reproductive system toxicity, lists representative cases of using this technology to explore reproductive damage caused by EDCs, and summarizes in detail the connection between environmental pollution and reproductive development disorders. It provides an important theoretical basis and direction for further exploring the mechanism of damage to the physiological functions of toxic substances on the reproductive system and the prevention and treatment of reproductive diseases.

Review: Environment of the ovulatory follicle: modifications and use of biotechnologies to enhance oocyte competence and increase fertility in cattle

The oocyte is the basis of life, supporting development from a fertilized cell to an independent multicellular organism. The oocyte's competence to drive the first cell cycles postfertilization are critical to embryonic survival and subsequent successful pregnancy. Coupled with the complex processes of follicle assembly, activation, differentiation, growth, and terminal maturation, oocyte developmental competence is gradually acquired during oocyte growth and meiotic maturation. Most reproduction management technologies and interventions are centered around these highly coordinated processes, targeting the ovarian follicle and the oocyte within. Thus, our objective was to highlight key aspects of oocyte and follicle development in cattle, and to discuss recent advances in oocyte and follicle-centered reproductive biotechnologies.

Single cell epigenomic and transcriptomic analysis uncovers potential transcription factors regulating mitotic/meiotic switch

In order to reveal the complex mechanism governing the mitotic/meiotic switch in female germ cells at epigenomic and genomic levels, we examined the chromatin accessibility (scATAC-seq) and the transcriptional dynamics (scRNA-seq) in germ cells of mouse embryonic ovary between E11.5 to 13.5 at single-cell resolution. Adopting a strict transcription factors (TFs) screening framework that makes it easier to understand the single-cell chromatin signature and a TF interaction algorithm that integrates the transcript levels, chromatin accessibility, and motif scores, we identified 14 TFs potentially regulating the mitotic/meiotic switch, including TCFL5, E2F1, E2F2, E2F6, E2F8, BATF3, SP1, FOS, FOXN3, VEZF1, GBX2, CEBPG, JUND, and TFDP1. Focusing on TCFL5, we constructed Tcfl5^+/- mice which showed significantly reduced fertility and found that decreasing TCFL5 expression in cultured E12.5 ovaries by RNAi impaired meiotic progression from leptotene to zygotene. Bioinformatics analysis of published results of the embryonic germ cell transcriptome and the finding that in these cells central meiotic genes (Stra8, Tcfl5, Sycp3, and E2f2) possess open chromatin status already at the mitotic stage together with other features of TCFL5 (potential capability to interact with core TFs and activate meiotic genes, its progressive activation after preleptotene, binding sites in the promoter region of E2f2 and Sycp3), indicated extensive amplification of transcriptional programs associated to mitotic/meiotic switch with an important contribution of TCFL5. We conclude that the identified TFs, are involved in various stages of the mitotic/meiotic switch in female germ cells, TCFL5 primarily in meiotic progression. Further investigation on these factors might give a significant contribution to unravel the molecular mechanisms of this fundamental process of oogenesis and provide clues about pathologies in women such as primary ovarian insufficiency (POI) due at least in part to meiotic defects.

Reproductive genomics of the mouse: implications for human fertility and infertility

Genetic analyses of mammalian gametogenesis and fertility have the potential to inform about two important and interrelated clinical areas: infertility and contraception. Here, we address the genetics and genomics underlying gamete formation, productivity and function in the context of reproductive success in mammalian systems, primarily mouse and human. Although much is known about the specific genes and proteins required for meiotic processes and sperm function, we know relatively little about other gametic determinants of overall fertility, such as regulation of gamete numbers, duration of gamete production, and gamete selection and function in fertilization. As fertility is not a binary trait, attention is now appropriately focused on the oligogenic, quantitative aspects of reproduction. Multiparent mouse populations, created by complex crossing strategies, exhibit genetic diversity similar to human populations and will be valuable resources for genetic discovery, helping to overcome current limitations to our knowledge of mammalian reproductive genetics. Finally, we discuss how what we know about the genomics of reproduction can ultimately be brought to the clinic, informing our concepts of human fertility and infertility, and improving assisted reproductive technologies.

M1AP interacts with the mammalian ZZS complex and promotes male meiotic recombination

Following meiotic recombination, each pair of homologous chromosomes acquires at least one crossover, which ensures accurate chromosome segregation and allows reciprocal exchange of genetic information. Recombination failure often leads to meiotic arrest, impairing fertility, but the molecular basis of recombination remains elusive. Here, we report a homozygous M1AP splicing mutation (c.1074 + 2T > C) in patients with severe oligozoospermia owing to meiotic metaphase I arrest. The mutation abolishes M1AP foci on the chromosome axes, resulting in decreased recombination intermediates and crossovers in male mouse models. M1AP interacts with the mammalian ZZS (an acronym for yeast proteins Zip2-Zip4-Spo16) complex components, SHOC1, TEX11, and SPO16. M1AP localizes to chromosomal axes in a SPO16-dependent manner and colocalizes with TEX11. Ablation of M1AP does not alter SHOC1 localization but reduces the recruitment of TEX11 to recombination intermediates. M1AP shows cytoplasmic localization in fetal oocytes and is dispensable for fertility and crossover formation in female mice. Our study provides the first evidence that M1AP acts as a copartner of the ZZS complex to promote crossover formation and meiotic progression in males.

A mouse model mimicking gender-affirming treatment with pubertal suppression followed by testosterone in transmasculine youth

STUDY QUESTION

Can mice serve as a translational model to examine the reproductive consequences of pubertal suppression with GnRH agonist (GnRHa) followed by testosterone (T) administration, a typical therapy in peripubertal transmasculine youth?

SUMMARY ANSWER

An implanted depot with 3.6 mg of GnRHa followed by T enanthate at 0.45 mg weekly can be used in peripubertal female mice for investigating the impact of gender-affirming hormone therapy in transmasculine youth.

WHAT IS KNOWN ALREADY

There is limited knowledge available in transgender medicine to provide evidence-based fertility care, with the current guidelines being based on the assumption of fertility loss. We recently successfully developed a mouse model to investigate the reproductive consequences of T therapy given to transgender men. On the other hand, to our knowledge, there is no mouse model to assess the reproductive outcomes in peripubertal transmasculine youth.

STUDY DESIGN, SIZE, DURATION

A total of 80 C57BL/6N female mice were used in this study, with n = 7 mice in each experimental group.

PARTICIPANTS/MATERIALS, SETTING, METHODS

We first assessed the effectiveness of GnRHa in arresting pubertal development in the female mice. In this experiment, 26-day-old female mice were subcutaneously implanted with a GnRHa (3.6 mg) depot. Controls underwent a sham surgery. Animals were euthanized at 3, 9, 21 and 28 days after the day of surgery. In the second experiment, we induced a transmasculine youth mouse model. C57BL/6N female mice were subcutaneously implanted with a 3.6 mg GnRHa depot on postnatal day 26 for 21 days and this was followed by weekly injections of 0.45 mg T enanthate for 6 weeks. The control for the GnRH treatment was sham surgery and the control for T treatment was sesame oil vehicle injections. Animals were sacrificed 0.5 weeks after the last injection. The data collected included the day of the vaginal opening and first estrus, daily vaginal cytology, weekly and terminal reproductive hormones levels, body/organ weights, ovarian follicular distribution and corpora lutea (CL) counts.

MAIN RESULTS AND THE ROLE OF CHANCE

GnRHa implanted animals remained in persistent diestrus and had reduced levels of FSH (P = 0.0013), LH (P = 0.0082) and estradiol (P = 0.0155), decreased uterine (P < 0.0001) and ovarian weights (P = 0.0002), and a lack of CL at 21 days after GnRHa implantation. T-only and GnRHa+T-treated animals were acyclic throughout the treatment period, had sustained elevated levels of T, suppressed LH levels (P < 0.0001), and an absence of CL compared to controls (P < 0.0001). Paired ovarian weights were reduced in the T-only and GnRHa+T groups compared with the control and GnRHa-only groups.

LARGE SCALE DATA

N/A.

LIMITATIONS, REASONS FOR CAUTION

Although it is an appropriate tool to provide relevant findings, precaution is needed to extrapolate mouse model results to mirror human reproductive physiology.

WIDER IMPLICATIONS OF THE FINDINGS

To our knowledge, this study describes the first mouse model mimicking gender-affirming hormone therapy in peripubertal transmasculine youth. This model provides a tool for researchers studying the effects of GnRHa-T therapy on other aspects of reproduction, other organ systems and transgenerational effects. The model is supported by GnRHa suppressing puberty and maintaining acyclicity during T treatment, lower LH levels and absence of CL. The results also suggest GnRHa+T therapy in peripubertal female mice does not affect ovarian reserve, since the number of primordial follicles was not affected by treatment.

STUDY FUNDING/COMPETING INTEREST(S)

This work was supported by the Michigan Institute for Clinical and Health Research grants KL2 TR 002241 and UL1 TR 002240 (C.D.C.); National Institutes of Health grants F30-HD100163 and T32-HD079342 (H.M.K.); University of Michigan Office of Research funding U058227 (A.S.); American Society for Reproductive Medicine/Society for Reproductive Endocrinology and Infertility grant (M.B.M.); and National Institutes of Health R01-HD098233 (M.B.M.). The University of Virginia Center for Research in Reproduction Ligand Assay and Analysis Core Facility was supported by the Eunice Kennedy Shriver NICHD/NIH grants P50-HD028934 and R24-HD102061. The authors declare that they have no competing interests.

Single-Cell Transcriptomics Analysis Reveals a Cell Atlas and Cell Communication in Yak Ovary

Yaks (Bos grunniens) are the only bovine species that adapt well to the harsh high-altitude environment in the Qinghai-Tibetan plateau. However, the reproductive adaptation to the climate of the high elevation remains to be elucidated. Cell composition and molecular characteristics are the foundation of normal ovary function which determines reproductive performance. So, delineating ovarian characteristics at a cellular molecular level is conducive to elucidating the mechanism underlying the reproductive adaption of yaks. Here, the single-cell RNA-sequencing (scRNA-seq) was employed to depict an atlas containing different cell types with specific molecular signatures in the yak ovary. The cell types were identified on the basis of their specifically expressed genes and biological functions. As a result, a cellular atlas of yak ovary was established successfully containing theca cells, stromal cells, endothelial cells, smooth muscle cells, natural killer cells, macrophages, and proliferating cells. A cell-to-cell communication network between the distinct cell types was constructed. The theca cells were clustered into five subtypes based on their biological functions. Further, CYP11A1 was confirmed as a marker gene for the theca cells by immunofluorescence staining. Our work reveals an ovarian atlas at the cellular molecular level and contributes to providing insights into reproductive adaption in yaks.

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.

Dynamic Autophagy Map in Mouse Female Germ Cells Throughout the Fetal to Postnatal Life

Autophagy plays vital roles in mouse female germ cells, but the potential mechanism is largely unknown. In this study, by interrogating single-cell RNA-seq dataset, we investigated the dynamic expression of autophagy-related genes in seven types of germ cells (mitosis, pre-leptotene, leptotene, zygotene, pachytene, diplotene, and dictyate) and discovered stage-specific autophagy-related genes. Using immunofluorescence (IF) and transmission electron microscopy (TEM), autophagy activity and autophagosome numbers were revealed from mitosis to follicular assembly (E12.5 (embryonic day 12.5) to P5 (postnatal day 5)). Furthermore, single-sample gene set enrichment analysis (ssGSEA) was performed to validate the autophagy kinetics from E12.5 to P5. Our study proved that the mitosis, diplotene, and dictyate female germ cells had relatively higher autophagy activity among the seven subtypes. In summary, our work provided an autophagy map, suggesting that autophagy was complicated in mouse female germ cell development from the fetal to postnatal life, which paved a new insight for deciphering the autophagy regulatory networks for cell-fate transition and female infertility issues like primary ovarian insufficiency (POI).

Single-cell profiling of lncRNAs in human germ cells and molecular analysis reveals transcriptional regulation of LNC1845 on LHX8

Non-coding RNAs exert diverse functions in many cell types. In addition to transcription factors from coding genes, non-coding RNAs may also play essential roles in shaping and directing the fate of germ cells. The presence of many long non-coding RNAs (lncRNAs) which are specifically expressed in the germ cells during human gonadal development were reported and one divergent lncRNA, LNC1845, was functionally characterized. Comprehensive bioinformatic analysis of these lncRNAs indicates that divergent lncRNAs occupied the majority of female and male germ cells. Integrating lncRNA expression into the bioinformatic analysis also enhances the cell-type classification of female germ cells. Functional dissection using in vitro differentiation of human pluripotent stem cells to germ cells revealed the regulatory role of LNC1845 on a transcription factor essential for ovarian follicle development, LHX8, by modulating the levels of histone modifications, H3K4me3 and H3K27Ac. Hence, bioinformatical analysis and experimental verification provide a comprehensive analysis of lncRNAs in developing germ cells and elucidate how an lncRNA function as a cis regulator during human germ cell development.

H3K4me3 as a target of di(2-ethylhexyl) phthalate (DEHP) impairing primordial follicle assembly

Di (2-ethylhexyl) phthalate (DEHP) is a widely used plastics additive that growing evidence indicates as endocrine disruptor able to negatively affect various reproductive processes both in female and male animals, including humans. However, the precise molecular mechanism of such actions is not completely understood. In the present study, scRNA-seq was performed on the ovaries of offspring from mothers exposed to DEHP from 16.5 days post coitum to 3 days post-partum, when the primordial follicle (PF) stockpile is established. While the histological observations of the offspring ovaries from DEHP exposed mothers confirmed previous data about a distinct reduction of oocytes enclosed in PFs. Focusing on oocytes, scRNA-seq analyses showed that the genes that mostly changed by DEHP were enriched GO terms related to histone H3-K4 methylation. Moreover, we observed H3K4me3 level, an epigenetics modification of H3 that is crucial for chromatin transcription, decreased by 40.28% (P < 0.01) in DEHP-treated group compared with control. When the newborn ovaries were cultured in vitro, the DEHP effects were abolished by tamoxifen (an estrogen receptor antagonist) or overexpression of Smyd3 (one specific methyltransferase of H3K4me3), in particular, the percentage of oocyte enclosed in PF was increased by 15.39% in DEHP plus Smyd3 overexpression group than of DEHP group (P < 0.01), which was accompanied by the upregulation of H3K4me3. Collectively, the present results discover Smyd3-H3K4me3 as a novel target of the deleterious ER-mediated effect of DEHP on PF formation during early folliculogenesis in the mouse and highlight epigenetics changes as prominent targets of endocrine disruptors like DEHP.

Transcriptional profiling of the developing rat ovary following intrauterine exposure to the endocrine disruptors diethylstilbestrol and ketoconazole

Exposure to endocrine-disrupting chemicals (EDCs) during development may cause reproductive disorders in women. Although female reproductive endpoints are assessed in rodent toxicity studies, a concern is that typical endpoints are not sensitive enough to detect chemicals of concern to human health. If so, measured endpoints must be improved or new biomarkers of effects included. Herein, we have characterized the dynamic transcriptional landscape of developing rat ovaries exposed to two well-known EDCs, diethylstilbestrol (DES) and ketoconazole (KTZ), by 3' RNA sequencing. Rats were orally exposed from day 7 of gestation until birth, and from postnatal day 1 until days 6, 14 or 22. Three exposure doses for each chemical were used: 3, 6 and 12 µg/kg bw/day of DES; 3, 6, 12 mg/kg bw/day of KTZ. The transcriptome changed dynamically during perinatal development in control ovaries, with 1137 differentially expressed genes (DEGs) partitioned into 3 broad expression patterns. A cross-species deconvolution strategy based on a mouse ovary developmental cell atlas was used to map any changes to ovarian cellularity across the perinatal period to allow for characterization of actual changes to gene transcript levels. A total of 184 DEGs were observed across dose groups and developmental stages in DES-exposed ovaries, and 111 DEGs in KTZ-exposed ovaries across dose groups and developmental stages. Based on our analyses, we have identified new candidate biomarkers for female reproductive toxicity induced by EDC, including Kcne2, Calb2 and Insl3.

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.

Multi-omics analysis reveals that iron deficiency impairs spermatogenesis by gut-hormone synthesis axis

Considering that research has mainly focussed on how excessive iron supplementation leads to reproductive cytotoxicity, there is a lack of in-depth research on reproductive system disorders caused by iron deficiency. To gain a better understanding of the effects of iron deficiency on the reproductive system, especially spermatogenesis, we first constructed a mouse model of iron deficiency. We employed multi-omic analysis, including transcriptomics, metabolomics, and microbiomics, to comprehensively dissect the impact of iron deficiency on spermatogenesis. Moreover, we verified our findings in detail using western blot, immunofluorescence, immunohistochemistry, qRT-PCR and other techniques. Microbiomic analysis revealed altered gut microbiota in iron-deficient mice, and functional predictive analysis showed that gut microbiota can regulate spermatogenesis. The transcriptomic data indicated that iron deficiency directly alters expression of meiosis-related genes. Transcriptome data also revealed that iron deficiency indirectly regulates spermatogenesis by affecting hormone synthesis, findings confirmed by metabolomic data, western blot and immunofluorescence. Interestingly, competing endogenous RNA networks also play a vital role in regulating spermatogenesis after iron deficiency. Taken together, the data elucidate that iron deficiency impairs spermatogenesis and increases the risk of male infertility by affecting hormone synthesis and promoting gut microbiota imbalance.

Single-cell transcriptomic profiling to evaluate the effects of Di(2-ethylhexyl)phthalate exposure on early meiosis of female mouse germ cells

Di(2-ethylhexyl)phthalate (DEHP) has proven characteristics of an endocrine-disrupting compound (EDC), which can threaten the reproductive health of humans and other animals. In mammals, a series of chromosomal events occur during the meiotic stage of oocytes. External toxins may enter the body and cause infertility and other related diseases. Therefore, it is crucial to explore the influence of DEHP exposure on the molecular mechanism of germ cell meiosis. We used single-cell RNA sequencing (scRNA-seq) to analyse the ovaries of foetal mice at embryonic day 12.5 (E12.5) and E14.5 after maternal DEHP exposure. DEHP exposure further activated the pathways related to DNA repair in germ cells, increased the expression of genes related to DNA damage and changed the developmental trajectory of germ cells. DEHP exposure may affect the proliferation of pregranulosa (PG) cells. Moreover, DEHP exposure altered the signal transduction between PG cells and germ cells. We showed that DEHP affects meiosis by causing DNA damage in oocytes and disrupting the signal transduction between PG cells and germ cells. These results provide a strong theoretical basis for the prevention and treatment of DEHP-mediated female reproductive health problems.

Single-nucleus RNA Sequencing reveals the mechanism of cigarette smoke exposure on diminished ovarian reserve in mice

The systematic toxicological mechanism of cigarette smoke (CS) on ovarian reserve has not been extensively investigated. Female 8-week-old C57BL/6 mice at peak fertility were exposed to CS or indoor air only for 30 days (100 mice per group) and the effects of CS on ovarian reserve were assessed using Single-Nucleus RNA Sequencing (snRNA-seq). In addition, further biochemical experiments, including immunohistochemical staining, ELISA, immunofluorescence staining, transmission electron microscopy, cell counting kit-8 assay, flow cytometry analysis, senescence-associated β-galactosidase staining, and western blotting, were accomplished to confirm the snRNA-seq results. We identified nine main cell types in adult ovaries and the cell-type-specific differentially expressed genes (DEGs) induced by CS exposure. Western blot results verified that down-regulation of antioxidant genes (Gpx1 and Wnt10b) and the steroid biosynthesis gene (Fdx1) occurred in both ovarian tissue and human granulosa cell-like tumor cell line (KGN cells) after CS exposure. Five percent cigarette smoke extract (CSE) effectively stimulated the production of reactive oxygen species (ROS), DNA damage, cellular senescence and markedly inhibited KGN cell proliferation by inducing G1-phase cell cycle arrest. Moreover, down-regulation of Gja1, Lama1 and the Ferroptosis indicator (Gpx4) in granulosa cells plays a significant role in ultrastructural changes in the ovary induced by CS exposure. These observations suggest that CS exposure impaired ovarian follicle reserve might be caused by REDOX imbalance in granulosa cells. The current study systematically determined the damage caused by CS in mouse ovaries and provides a theoretical basis for early clinical prediction, diagnosis and intervention of CS exposure-associated primary ovarian insufficiency (POI), and is of great significance in improving female reproductive health.

BNC1 deficiency-triggered ferroptosis through the NF2-YAP pathway induces primary ovarian insufficiency

Primary ovarian insufficiency (POI) is a clinical syndrome of ovarian dysfunction characterized by premature exhaustion of primordial follicles. POI causes infertility, severe daily life disturbances and long-term health risks. However, the underlying mechanism remains largely unknown. We previously identified a Basonuclin 1 (BNC1) mutation from a large Chinese POI pedigree and found that mice with targeted Bnc1 mutation exhibit symptoms of POI. In this study, we found that BNC1 plays key roles in ovarian reserve and maintaining lipid metabolism and redox homeostasis in oocytes during follicle development. Deficiency of BNC1 results in premature follicular activation and excessive follicular atresia. Mechanistically, BNC1 deficiency triggers oocyte ferroptosis via the NF2-YAP pathway. We demonstrated that pharmacologic inhibition of YAP signaling or ferroptosis significantly rescues Bnc1 mutation-induced POI. These findings uncover a pathologic mechanism of POI based on BNC1 deficiency and suggest YAP and ferroptosis inhibitors as potential therapeutic targets for POI.

Primary ovarian insufficiency (POI) is a clinical syndrome of ovarian dysfunction that results in infertility. Here they show that BCN1 mutation results in premature ovarian follicle activation and atresia through dysregulation of ferroptosis.

Effects and potential mechanism of Ca2+/calmodulin‑dependent protein kinase II pathway inhibitor KN93 on the development of ovarian follicle

Adequate regulation of the speed of follicular development has been reported to prolong the reproductive life of the ovary. The aim of the present study was to assess the potential effects and mechanism of the Ca²⁺/calmodulin-dependent protein kinase II (CaMKII) pathway on the development of ovarian follicle. In the present study, the expression of CaMKII was measured in the ovary of mice at different developmental stages by immunofluorescence, confirming that CaMKII has a role in follicular development. Subsequently, the 17.5 days post-coitus (dpc) embryonic ovaries were collected and cultured with KN93 for 4 days in vitro. It was revealed that KN93 inhibited the development of follicles, where it reduced the expression levels of oocyte and granulosa cell markers DEAD-box helicase 4 (DDX4) and forkhead box L2 (FOXL2). These results suggested that KN93 could delay follicular development. Proteomics technology was then used to find that 262 proteins of KN93 treated 17.5 dpc embryonic ovaries were significantly altered after in vitro culture. Bioinformatics analysis was used to analyze these altered proteins. In total, four important Kyoto Encyclopedia of Genes and Genome pathways, namely steroid biosynthesis, p53 signaling pathway and retinol metabolism and metabolic pathways, were particularly enriched. Further analysis revealed that the upregulated proteins NADP-dependent steroid dehydrogenase-like (Nsdhl), lanosterol synthase (Lss), farnesyl-diphosphate farnesyltransferase 1 (Fdft1), cytochrome P450 family 51 family A member 1 (Cyp51a1), hydroxymethylglutaryl-CoA synthase 1 (Hmgcs1), fatty acid synthase (Fasn) and dimethylallyltranstransferase (Fdps) were directly interacting with each other in the four enriched pathways. In summary, the potential mechanism of KN93 in slowing down follicular development most likely lies in its inhibitory effects on CaMKII, which upregulated the expression of Nsdhl, Lss, Fdft1, Cyp51a1, Hmgcs1, Fasn and Fdps. This downregulated the expression of oocyte and granulosa cell markers DDX4 and FOXL2 in the follicles, thereby delaying follicular development. Overall, these results provide novel insight into the potential mechanism by which KN93 and CaMKII can delay follicular development.

A New Understanding, Guided by Single-Cell Sequencing, of the Establishment and Maintenance of the Ovarian Reserve in Mammals

BACKGROUND

Oocytes are a finite and non-renewable resource that are maintained in primordial follicle structures. The ovarian reserve is the totality of primordial follicles, present from birth, within the ovary and its establishment, size, and maintenance dictates the duration of the female reproductive lifespan. Understanding the cellular and molecular dynamics relevant to the establishment and maintenance of the reserve provides the first steps necessary for modulating both individual human and animal reproductive health as well as population dynamics.

SUMMARY

This review details the key stages of establishment and maintenance of the ovarian reserve, encompassing germ cell nest formation, germ cell nest breakdown, and primordial follicle formation and activation. Furthermore, we spotlight several formative single-cell sequencing studies that have significantly advanced our knowledge of novel molecular regulators of the ovarian reserve, which may improve our ability to modulate female reproductive lifespans.

KEY MESSAGES

The application of single-cell sequencing to studies of ovarian development in mammals, especially when leveraging genetic and environmental models, offers significant insights into fertility and its regulation. Moreover, comparative studies looking at key stages in the development of the ovarian reserve across species has the potential to impact not just human fertility, but also conservation biology, invasive species management, and agriculture.

Bone morphogenetic protein 2- and estradiol-17β-induced changes in ovarian transcriptome during primordial follicle formation†

Estradiol-17β has been shown to promote primordial follicle formation and to involve bone morphogenetic protein 2 (BMP2) as a downstream effector to promote primordial follicle in hamsters. However, the molecular mechanism whereby these factors regulate ovarian somatic cells to pre-granulosa cells transition leading to primordial follicle formation remains unclear. The objective of this study was to determine whether BMP2 and/or estradiol-17β would regulate the expression of specific ovarian transcriptome during pre-granulosa cells transition and primordial follicle formation in the mouse ovary. BMP2 mRNA level increased during the period of primordial follicle formation with the concurrent presence of BMP2 protein in ovarian somatic cells. Estradiol-17β but not BMP2 exposure led to increased expression of ovarian BMP2 messenger RNA (mRNA), and the effect of estradiol-17β could not be suppressed by 4-[6-[4-(1-Piperazinyl)phenyl]pyrazolo[1,5-a]pyrimidin-3-yl]quinoline dihydrochloride (LDN) 193189. BMP2 or estradiol-17β stimulated primordial follicle formation without inducing apoptosis. Ribonucleic acid-sequence analysis (RNA-seq) of ovaries exposed to exogenous BMP2 or estradiol-17β revealed differential expression of several thousand genes. Most of the differentially expressed genes, which were common between BMP2 or estradiol-17β treatment demonstrated concordant changes, suggesting that estradiol-17β and BMP2 affected the same set of genes during primordial follicle formation. Further, we have identified that estradiol-17β, in cooperation with BMP2, could affect the expression of three major transcription factors, GATA binding protein 2, GATA binding protein 4 and Early growth response 2, and one serine protease, hepsin, in pre-granulosa cells during primordial follicle formation. Taken together, results of this study suggest that estradiol-17β and BMP2 may regulate ovarian gene expression that promote somatic cells to pre-granulosa cells transition and primordial follicle formation in the mouse ovary.

Maternal ochratoxin A exposure impairs meiosis progression and primordial follicle formation of F1 offspring

Ochratoxin A (OTA), a mycotoxin produced by Aspergillus and Penicillium fungi, widely contaminates feed, food and their raw materials. OTA has been proved to have hepatotoxicity and nephrotoxicity. Its reproductive toxicity needs to be further explored. We found that OTA inhibited the progression of meiosis, keeping more germ cells at leptotene and zygotene. Furthermore, OTA impaired primordial follicle formation, keeping more germ cells in cysts. Increased γH2AX suggested that DNA damage occurred both at the stages of meiosis and primordial follicle formation. The expression of RAD51 increased with the concentration of OTA at the stage of meiosis, while decreased later, suggesting the activated DNA repair induced by DNA damage then inhibited by persistent and excessive stress of DNA damage, which further induced apoptosis. DEGs caused by OTA were also mainly enriched in DNA damage and repair through RNA-seq analysis. Higher level of reactive oxygen species (ROS) and increased degree of oxidative damage marker 8-OHdG were both found in the ovaries exposed to OTA. We concluded that maternal OTA exposure affected meiosis progression and primordial follicle formation via oxidative damage and DNA repair. Clarification of the mechanism of OTA will contribute to the development of more effective detoxification strategies.

Advances in application of single-cell RNA sequencing in cardiovascular research

Single-cell RNA sequencing (scRNA-seq) provides high-resolution information on transcriptomic changes at the single-cell level, which is of great significance for distinguishing cell subtypes, identifying stem cell differentiation processes, and identifying targets for disease treatment. In recent years, emerging single-cell RNA sequencing technologies have been used to make breakthroughs regarding decoding developmental trajectories, phenotypic transitions, and cellular interactions in the cardiovascular system, providing new insights into cardiovascular disease. This paper reviews the technical processes of single-cell RNA sequencing and the latest progress based on single-cell RNA sequencing in the field of cardiovascular system research, compares single-cell RNA sequencing with other single-cell technologies, and summarizes the extended applications and advantages and disadvantages of single-cell RNA sequencing. Finally, the prospects for applying single-cell RNA sequencing in the field of cardiovascular research are discussed.

Dysfunctional Ovarian Stem Cells Due to Neonatal Endocrine Disruption Result in PCOS and Ovarian Insufficiency in Adult Mice

Polycystic ovarian syndrome (PCOS) is a common global cause of anovulatory infertility but underlying etiology leading to PCOS still remains elusive. Fetal and perinatal endocrine disruption reportedly affects germ cell nests (GCN) breakdown, meiosis, and primordial follicle (PF) assembly with unassembled oocytes in neonatal ovaries. We recently reported that very small embryonic-like stem cells (VSELs) and ovarian stem cells (OSCs) express ERα, ERβ and FSHR, undergo distinct cyclic changes and neo-oogenesis encompassing GCN formation, meiosis, and primordial follicle (PF) assembly on regular basis in adult mice ovaries and these GCN are arrested in pre-meiotic or early meiotic stage in aged ovaries. Present study was undertaken to evaluate whether neonatal exposure to endocrine disruption (estradiol E2 or diethylstilbestrol DES) affects ovarian stem cells and their differentiation (neo-oogenesis) and PF assembly in adult 100 days old ovaries. Neonatal exposure to E2 resulted in typical features of PCOS including hyperandrogenism, infertility, increased stromal compartment, absent corpus lutea, and cystic follicles whereas DES treated ovaries showed rapid recruitment of follicles in young ovaries and multi-ovular/cystic follicles. Ovary surface epithelial cells smears showed large numbers of growth-arrested GCN in zygotene/pachytene with increased expression of Mlh-1 and Scp-1 suggesting defects at synapsis and recombination stages during prophase-1 of meiosis. Being immortal and expression of ERα and ERβ makes VSELs directly vulnerable to carry developmental endocrine insults to adult life. Dysfunction of VSELs/OSCs possibly results in oocyte defects observed in our study in PCOS/POI besides the widely reported defects in granulosa cells.

The Requirement of Ubiquitin C-Terminal Hydrolase L1 (UCHL1) in Mouse Ovarian Development and Fertility †

Ubiquitin C-Terminal Hydrolase L1 (UCHL1) is a de-ubiquitinating enzyme enriched in neuronal and gonadal tissues known to regulate the cellular stores of mono-ubiquitin and protein turnover. While its function in maintaining proper motor neuron function is well-established, investigation into its role in the health and function of reproductive processes is only just beginning to be studied. Single-cell-sequencing analysis of all ovarian cells from the murine perinatal period revealed that Uchl1 is very highly expressed in the developing oocyte population, an observation which was corroborated by high levels of oocyte-enriched UCHL1 protein expression in oocytes of all stages throughout the mouse reproductive lifespan. To better understand the role UCHL1 may be playing in oocytes, we utilized a UCHL1-deficient mouse line, finding reduced number of litters, reduced litter sizes, altered folliculogenesis, morphologically abnormal oocytes, disrupted estrous cyclicity and apparent endocrine dysfunction in these animals compared to their wild-type and heterozygous littermates. These data reveal a novel role of UCHL1 in female fertility as well as overall ovarian function, and suggest a potentially essential role for the ubiquitin proteasome pathway in mediating reproductive health. Summary sentence: Ubiquitin C-Terminal Hydrolase L1 (UCHL1) is required for proper ovarian folliculogenesis, estrous cyclicity, and fertility in the female mouse.

Impaired primordial follicle assembly in offspring ovaries from zearalenone-exposed mothers involves reduced mitochondrial activity and altered epigenetics in oocytes

Previous works have shown that zearalenone (ZEA), as an estrogenic pollutant, has adverse effects on mammalian folliculogenesis. In the present study, we found that prolonged exposure of female mice to ZEA around the end of pregnancy caused severe impairment of primordial follicle formation in the ovaries of newborn mice and altered the expression of many genes in oocytes as revealed by single-cell RNA sequencing (scRNA-seq). These changes were associated with morphological and molecular alterations of mitochondria, increased autophagic markers in oocytes, and epigenetic changes in the ovaries of newborn mice from ZEA-exposed mothers. The latter increased expression of HDAC2 deacetylases was leading to decreased levels of H3K9ac and H4K12ac. Most of these modifications were relieved when the expression of  Hdac2 in newborn ovaries was reduced by RNA interference during in vitro culture in the presence of ZEA. Such changes were also alleviated in offspring ovaries from mothers treated with both ZEA and the coenzyme Q10 (CoQ10), which is known to be able to restore mitochondrial activities. We concluded that impaired mitochondrial activities in oocytes caused by ZEA are at the origin of metabolic alterations that modify the expression of genes controlling autophagy and primordial follicle assembly through changes in epigenetic histones.