Transcriptome Landscape of Human Folliculogenesis Reveals Oocyte and Granulosa Cell Interactions

Downregulation of glucose-energy metabolism via AMPK signaling pathway in granulosa cells of diminished ovarian reserve patients

Glucose metabolism plays a crucial role in the function of granulosa cells (GCs) and the development of follicles. In cases of diminished ovarian reserve (DOR), alterations in these processes can impact female fertility. This study aims to investigate changes in glucose-energy metabolism in GCs of young DOR patients aged 20 to 35 years and their correlation with the onset and progression of DOR. 72 DOR cases and 75 women with normal ovarian reserve (NOR) as controls were included based on the POSEIDON and Bologna criteria. Samples of GCs and follicular fluid (FF) were collected for a comprehensive analysis involving transcriptomics, metabolomics, RT-qPCR, JC-1 staining, and flow cytometry. The study identified differentially expressed genes and metabolites in GCs of DOR and NOR groups, revealing 7 common pathways related to glucose-energy metabolism, along with 11 downregulated genes and 14 metabolites. Key substances in the glucose-energy metabolism pathway, such as succinate, lactate, NADP, ATP, and ADP, showed decreased levels, with the DOR group exhibiting a reduced ADP/ATP ratio. Downregulation of genes involved in glycolysis (HK, PGK, LDH1), the TCA cycle (CS), and gluconeogenesis (PCK) was observed, along with reduced glucose content and expression of glucose transporter genes (GLUT1 and GLUT3) in DOR GCs. Additionally, decreased AMPK pathway activity and impaired mitochondrial function in DOR suggest a connection between mitochondrial dysfunction and disrupted energy metabolism. Above all, the decline in glucose-energy metabolism in DOR is closely associated with its onset and progression. Reduced glucose uptake and impaired mitochondrial function in DOR GCs lead to internal energy imbalances, hindering the AMPK signaling pathway, limiting energy production and supply, and ultimately impacting follicle development and maturation.

Luteinizing Hormone Receptor Mutation (LHRN316S) Causes Abnormal Follicular Development Revealed by Follicle Single-Cell Analysis and CRISPR/Cas9

Abnormal interaction between granulosa cells and oocytes causes disordered development of ovarian follicles. However, the interactions between oocytes and cumulus granulosa cells (CGs), oocytes and mural granulosa cells (MGs), and CGs and MGs remain to be fully explored. Using single-cell RNA-sequencing (scRNA-seq), we determined the transcriptional profiles of oocytes, CGs and MGs in antral follicles. Analysis of scRNA-seq data revealed that CGs may regulate follicular development through the BMP15-KITL-KIT-PI3K-ARF6 pathway with elevated expression of luteinizing hormone receptor (LHR). Because internalization of the LHR is regulated by Arf6, we constructed LHRN316S mice by CRISPR/Cas9 to further explore mechanisms of follicular development and novel treatment strategies for female infertility. Ovaries of LHRN316S mice exhibited reduced numbers of corpora lutea and ovulation. The LHRN316S mice had a reduced rate of oocyte maturation in vitro and decreased serum progesterone levels. Mating LHRN316S female mice with ICR wild type male mice revealed that the infertility rate of LHRN316S mice was 21.4% (3/14). Litter sizes from LHRN316S mice were smaller than those from control wild type female mice. The oocytes from LHRN316S mice had an increased rate of maturation in vitro after progesterone administration in vitro. Furthermore, progesterone treated LHRN316S mice produced offspring numbers per litter equivalent to WT mice. These findings provide key insights into cellular interactions in ovarian follicles and provide important clues for infertility treatment.

Stem cell transplantation extends the reproductive life span of naturally aging cynomolgus monkeys

The ovary is crucial for female reproduction and health, as it generates oocytes and secretes sex hormones. Transplantation of mesenchymal stem cells (MSCs) has been shown to alleviate pathological ovarian aging. However, it is unclear whether MSCs could benefit the naturally aging ovary. In this study, we first examined the dynamics of ovarian reserve of Chinese women during perimenopause. Using a naturally aging cynomolgus monkey (Macaca fascicularis) model, we found that transplanting human embryonic stem cells-derived MSC-like cells, which we called M cells, into the aging ovaries significantly decreased ovarian fibrosis and DNA damage, enhanced secretion of sex hormones and improved fertility. Encouragingly, a healthy baby monkey was born after M-cell transplantation. Moreover, single-cell RNA sequencing analysis and in vitro functional validation suggested that apoptosis, oxidative damage, inflammation, and fibrosis were mitigated in granulosa cells and stromal cells following M-cell transplantation. Altogether, these findings demonstrate the beneficial effects of M-cell transplantation on aging ovaries and expand our understanding of the molecular mechanisms underlying ovarian aging and stem cell-based alleviation of this process.

Comparison of chromatin accessibility remodeling of granulosa cells in patients with endometrioma or pelvic/tubal infertility

PURPOSE

To elucidatethe epigenetic alteration associated with impaired oogenesis in endometrioma using multi-omic approaches.

METHODS

ATAC-seq was performed on the granulosa cells (GCs) of 6 patients (3 with endometrioma and 3 without). Follicular samples from another 20 patients (10 with endometrioma and 10 without) were collected for mRNA-seq analysis of GCs and extracellular vesicles (EVs) of follicular fluid. qRT-PCR validated candidate genes in GCs from 44 newly enrolled patients (19 with endometrioma and 25 without). mRNA abundance was compared with the Mann-Whitney test. Pearson's correlation analyzed relationships between candidate genes and oocyte parameters.

RESULTS

Chromatin accessibility and gene expression profiles of GCs from endometrioma patients differed significantly from the pelvic/tubal infertility group. RNA-seq revealed most differentially expressed genes were downregulated (6216/7325) and enriched in the cellular localization pathway. Multi-omics analyses identified 22 significantly downregulated genes in the GCs of endometrioma patients, including PPIF (P < 0.0001) and VEGFA (P = 0.0148). Both genes were further confirmed by qRT-PCR. PPIF (r = 0.46, p = 0.043) and VEGFA (r = 0.45, p = 0.048) correlated with the total number of retrieved oocytes.

CONCLUSIONS

GC chromatin remodeling may disrupt GC and EV transcriptomes, interfering with somatic cell-oocyte communication and leading to compromised oogenesis in endometrioma patients.

DIA proteomic and PRM validation through human granulose cells profiles screen suitable biomarkers for polycystic ovary syndrome patients

The aim of this study is to identify differentially expressed proteins (DEPs) in granulose cells (GCs) from women with or withoutpolycystic ovary syndrome (PCOS) via data independent acquisition (DIA) proteomic analysis.A total of 63 women were recruited for this study, 34 PCOS patients as experimental group (P), and 29 women without PCOS as Normal group (NP). DIA-based proteomic analysis was performed to identify DEPs in GCs between the P and NP samples. Certain typical DEPs were further validated by Parallel reaction monitoring (PRM), and correlation analysis was performed between these DEPs and the clinical characteristics.Cell vitality was measured by CCK-8 assay. DIA analysis revealed 174 significantly DEPs, of which 7 were upregulated and 167 downregulated. Bioinformatics analysis was performed to analysis the significantly DEPs. The PRM experiment confirmed TOP2A and SPHKAP were upregulated significantly in P by comparing to NP, while GM2A, MRPS16, APOA2 and FGF2 were downregulated significantly. Most notably, Correlation analysis revealed that TOP2A, SPHKAP, MRPS16 and FGF2were positively correlated with TG, AMH and Age, but negatively correlated with Menarche age, DBIL, FT3, Basal serum FSH and LH.Meanwhile, CCK-8 assay has shown that downregulation of FGF2 could weaken cell viability. Finally, a panel of DEPs were identified in the GCs of patients with PCOS, of which certain significant DEPs might play essential roles in the pathogenesis of PCOS, could be regarded as candidate biomarkers for PCOS.

Decoding molecular features of bovine oocyte fate during antral follicle growth via single-cell multi-omics analysis†

Antral follicle size is a useful predictive marker of the competency of enclosed oocytes for yielding an embryo following in vitro maturation and fertilization. However, the molecular mechanisms underpinning oocyte developmental potential during bovine antral follicle growth are still unclear. Here, we used a modified single-cell multi-omics approach to analyze the transcriptome, DNA methylome, and chromatin accessibility in parallel for oocytes and cumulus cells collected from bovine antral follicles of different sizes. Transcriptome profiling identified three types of oocytes (small, medium, and large) that underwent different developmental trajectories, with large oocytes exhibiting the largest average follicle size and characteristics resembling metaphase-II oocytes. Differential expression analysis and real-time polymerase chain reaction assay showed that most replication-dependent histone genes were highly expressed in large oocytes. The joint analysis of multi-omics data revealed that the transcription of 20 differentially expressed genes in large oocytes was associated with both DNA methylation and chromatin accessibility. In addition, oocyte-cumulus interaction analysis showed that inflammation, DNA damage, and p53 signaling pathways were active in small oocytes, which had the smallest average follicle sizes. We further confirmed that p53 pathway inhibition in the in vitro maturation experiments using oocytes obtained from small antral follicles could improve the quality of oocytes and increased the blastocyte rate after in vitro fertilization and culture. Our work provides new insights into the intricate orchestration of bovine oocyte fate determination during antral folliculogenesis, which is instrumental for optimizing in vitro maturation techniques to optimize oocyte quality.

Evidence of apoptosis as an early event leading to cyclophosphamide-induced primordial follicle depletion in a prepubertal mouse model

Introduction

The mechanisms leading to ovarian primordial follicle depletion following gonadotoxic chemotherapy with cyclophosphamide and other cytotoxic drugs are currently understood through two main explanatory theories: apoptosis and over-activation. Discrepancies between the findings of different studies investigating these mechanisms do not allow to reach a firm conclusion. The heterogeneity of cell types in ovaries and their different degrees of sensitivity to damage, cell-cell interactions, periodical follicle profile differences, model age-dependent differences, and differences of exposure durations of tested drugs may partially explain the discrepancies among studies.

Methods

This study used intact prepubertal mice ovaries in culture as study model, in which most follicles are primordial follicles. Histological and transcriptional analyses of ovaries exposed to the active metabolite of cyclophosphamide 4-hydroperoxycyclophosphamide (4-HC) were carried out via a time-course experiment at 8, 24, 48, and 72 h.

Results

4-HC treated ovaries showed a significant decrease in primordial follicle density at 24 h, along with active DNA damage (TUNEL) and overexpressed apoptosis signals (cleaved-poly ADP ribose polymerase in immunohistochemistry and western blotting). Meanwhile 4-HC treatment significantly up-regulated H2ax, Casp 6, Casp 8, Noxa, and Bax in ovaries, and up-regulated Puma in primordial follicles (FISH).

Discussion

Our results indicated that cyclophosphamide-induced acute ovarian primordial follicle depletion was mainly related to apoptotic pathways. No evidence of follicle activation was found, neither through changes in the expression of related genes to follicle activation nor in the density of growing follicles. Further validation at protein level in 4-HC-treated prepubertal mice ovaries at 24 h confirmed these observations.

A spotlight on factors influencing the in vitro folliculogenesis of isolated preantral follicles

Female fertility preservation via complete in vitro folliculogenesis is still chimerical. Due to many factors affecting the efficiency of isolation and culture of preantral follicles, the improvement of techniques geared to fertility preservation in higher mammals seems to be at an impasse. We need an objective view of the current stand to understand how to progress further. As such, a survey was conducted to analyze the relative distribution of studies performed in ten mammalian species on preantral follicle culture available on PubMed. Using the bovine as a reference model, we explore some factors influencing data variation that contribute to the difficulty in reproducing studies. While years of research have enabled the recapitulation of folliculogenesis from as modest as the early antral follicle stage ex vivo, in vitro preantral folliculogenesis remains elusive. Herein, we revisit the classical evidence that laid the foundations for understanding preantral folliculogenesis and review the length, breadth, and depth of information that the era of big data has currently levied. Moving forward, we recognize the urgency of synthesizing the multi-disciplinary approaches to mimic folliculogenesis in vitro to achieve a translational landscape of infertility at individual and large-scale conservation levels.

TET1 overexpression affects cell proliferation and apoptosis in aging ovaries

PURPOSE

Along with the progress of society, human life expectancy has been increasing, and late marriage and late childbearing are the current trend. Since reproductive aging affects fertility, ovarian aging in women has become a major reproductive health issue in the current society. During ovarian aging, DNA methylation levels may change. The ten-eleven translocation (TET) protein family proteins TET1, TET2, and TET3 are important DNA demethylation enzymes, and differential expression of TET1, TET2, and TET3 may affect the proliferation and apoptosis of aging ovarian cells. The aim of this study was to investigate the role of TET1 in the regulation of ovarian aging.

METHODS

The expression of 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) was analyzed by immunofluorescence (IF) in young and aging ovaries of six 6-8-week-old female mice and six 6-8-month-old female mice. Then, the expression pattern of the TET protein family in young and aging ovaries of mice was investigated. To determine the impact of TET1 on ovarian development, the aging of IOSE-80, KGN, and SKOV-3 cells was induced with D-galactosidase (D-gal). Cells were then transfected using the TET1 overexpression vector or si-TET1. We assessed the proliferation and apoptosis of aging cells after transfection and analyzed the regulatory effect of TET1 expression on aging cells. Additionally, we verified the Tet1 expression in Tet1-KO mice.

RESULTS

The 5mC to 5hmC transition, oocyte maturation, and blastocyst rate were reduced in aging mice compared to young mice. In aging mice ovaries, the expression levels of Tet1, Tet2, and Tet3 were reduced significantly, with Tet1 being particularly pronounced. The overexpression of TET1 promoted proliferation and inhibited apoptosis in aging human ovarian cells. Furthermore, Tet1 expression was very low in Tet1-KO C57BL/6 J mice ovaries.

CONCLUSION

This study demonstrates that the expression levels of TET family proteins are low in aging ovaries, and the overexpression of TET1 can promote proliferation and inhibit apoptosis in aging ovarian 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.

SMARTdb: An Integrated Database for Exploring Single-cell Multi-omics Data of Reproductive Medicine

Single-cell multi-omics sequencing has greatly accelerated reproductive research in recent years, and the data are continually growing. However, utilizing these data resources is challenging for wet-lab researchers. A comprehensive platform for exploring single-cell multi-omics data related to reproduction is urgently needed. Here, we introduce the single-cell multi-omics atlas of reproduction (SMARTdb), an integrative and user-friendly platform for exploring molecular dynamics of reproductive development, aging, and disease, which covers multi-omics, multi-species, and multi-stage data. We curated and analyzed single-cell transcriptomic and epigenomic data of over 2.0 million cells from 6 species across the entire lifespan. A series of powerful functionalities are provided, such as "Query gene expression", "DIY expression plot", "DNA methylation plot", and "Epigenome browser". With SMARTdb, we found that the male germ cell-specific expression pattern of RPL39L and RPL10L is conserved between human and other model animals. Moreover, DNA hypomethylation and open chromatin may collectively regulate the specific expression pattern of RPL39L in both male and female germ cells. In summary, SMARTdb is a powerful platform for convenient data mining and gaining novel insights into reproductive development, aging, and disease. SMARTdb is publicly available at https://smart-db.cn.

Genetic links between ovarian ageing, cancer risk and de novo mutation rates

Human genetic studies of common variants have provided substantial insight into the biological mechanisms that govern ovarian ageing1. Here we report analyses of rare protein-coding variants in 106,973 women from the UK Biobank study, implicating genes with effects around five times larger than previously found for common variants (ETAA1, ZNF518A, PNPLA8, PALB2 and SAMHD1). The SAMHD1 association reinforces the link between ovarian ageing and cancer susceptibility1, with damaging germline variants being associated with extended reproductive lifespan and increased all-cause cancer risk in both men and women. Protein-truncating variants in ZNF518A are associated with shorter reproductive lifespan-that is, earlier age at menopause (by 5.61 years) and later age at menarche (by 0.56 years). Finally, using 8,089 sequenced trios from the 100,000 Genomes Project (100kGP), we observe that common genetic variants associated with earlier ovarian ageing associate with an increased rate of maternally derived de novo mutations. Although we were unable to replicate the finding in independent samples from the deCODE study, it is consistent with the expected role of DNA damage response genes in maintaining the genetic integrity of germ cells. This study provides evidence of genetic links between age of menopause and cancer risk.

The Differentiation Fate of Granulosa Cells and the Regulatory Mechanism in Ovary

Granulosa cells (GCs) are important drives of the reproductive process, not only the supporting cells for nutrition, but also cells with endocrine functions. Their differentiation and development parallel the entire menstruation period and even during pregnancy, making it tightly linked to the fate of the follicle. To elucidate the underlying mechanism is of great significance for related researches. The life course of GCs is briefly divided into five stages, from epithelial cells to pre-granulosa cells, GCs, mural and cumulus cells, lutein cells, and eventually disappear. A wide variety of genes and transcription factors participate in the regulation of different stages, and more importantly, various hormones secreted by the pituitary gland and GCs themselves play a leading role. These endogenous and exogenous signalling molecules interact to form a cross-linked communication network, promoting the development of GCs. Together with oocytes, theca cells and other functional cells in the ovary, GCs drive one of the most vital biological processes in women.

In-depth analysis of transcriptomes in ovarian cortical follicles from children and adults reveals interfollicular heterogeneity

The ovarian cortical reserve of follicles is vital for fertility. Some medical treatments are toxic to follicles, leading to premature ovarian insufficiency. Ovarian tissue cryopreservation is an established method to preserve fertility in adults and even applied in prepuberty despite unproven efficacy. Here, we analyze transcriptomes of 120 cortical follicles from children and adults for detailed comparison. We discover heterogeneity with two main types of follicles in both age groups: one with expected oocyte-granulosa profiles and another with predicted role in signaling. Transcriptional changes during growth to the secondary stage are similar overall in children and adults, but variations related to extracellular matrix, theca cells, and miRNA profiles are found. Notably, cyclophosphamide dose correlates with interferon signaling in child follicles. Additionally, morphology alone is insufficient for follicle categorization suggesting a need for additional markers. Marker genes for early follicle activation are determined. These findings will help refine follicular classification and fertility preservation techniques across critical ages.

Liraglutide improves follicle development in polycystic ovary syndrome by inhibiting CXCL10 secretion

BACKGROUND

At present, a number of clinical trials have been carried out on GLP-1 receptor agonist liraglutide in the treatment of polycystic ovary syndrome (PCOS). However, the effect of liraglutide on follicle development and its specific mechanism are still unclear.

METHODS

RNA sequencing was used to explore the molecular characteristics of granulosa cells from patients with PCOS treated with liraglutide. The levels of C-X-C motif chemokine ligand 10 (CXCL10) in follicular fluid were detected by ELISA, the expression levels of ovulation related genes and inflammatory factor genes in follicles and granulosa cells were detected by qPCR and the protein levels of connexin 43 (Cx43), Janus Kinase 2 (JAK2) and phosphorylated JAK2 were detected by Western blot. The mouse ovarian follicles culture system in vitro was used to detect the status of follicle development and ovulation.

RESULTS

In the present study, we found that liraglutide inhibited the secretion of inflammatory factors in PCOS granulosa cells, among which CXCL10 was the most significant. In addition, CXCL10 was significantly higher in granulosa cells and follicular fluid in PCOS patients than in non-PCOS patients. We applied in vitro follicle culture and other techniques to carry out the mechanism exploration which revealed that CXCL10 disrupted the homeostasis of gap junction protein alpha 1 (GJA1) between oocyte and granulosa cells before physiological ovulation, thus inhibiting follicular development and ovulation. Liraglutide inhibited CXCL10 secretion in PCOS granulosa cells by inhibiting the JAK signaling pathway and can improved dehydroepiandrosterone (DHEA)-induced follicle development disorders, which is reversed by CXCL10 supplementation.

CONCLUSIONS

The present study suggests that liraglutide inhibits CXCL10 secretion in granulosa cells through JAK signaling pathway, thereby improving the homeostasis of GJA1 between oocyte and granulosa cells before physiological ovulation and ultimately improving the follicular development and ovulation of PCOS, which provides more supportive evidence for the clinical application of liraglutide in the treatment of ovulatory disorders in PCOS.

TRIAL REGISTRATION

Not applicable.

Enhancing Oocyte Quality in Aging Mice: Insights from Mesenchymal Stem Cell Therapy and FOXO3a Signaling Pathway Activation

Ovarian aging reduced the quality of oocytes, resulting in age-related female infertility. It is reported that mesenchymal stem cells (MSCs) therapy can improve age-related ovarian function decline and the success rate of in vitro maturation (IVM) in assisted reproductive therapy. In order to investigate the effectiveness and mechanisms of MSCs to enhance oocyte quality of cumulus oocyte complexes (COCs) in advanced age, this study focus on the respective functional improvement of oocytes and granulosa cells (GCs) from aging mice and further to explore and verify the possible mechanisms. Here, we studied a popular but significant protein of follicular development, Forkhead box O-3a (FOXO3a), which is a transcription factor that mediates a variety of cellular processes, but the functions of which in regulating oocyte quality in MSCs therapy still remain inconclusive. In this study, the RNA-seq data of metaphase II (MII) oocytes and GCs isolated from COCs confirmed that, GCs of immature follicles show the most potential to be the targeted cells of bone marrow mesenchymal stem cells (BMSCs) by FOXO3a signaling pathway. Furthermore, we demonstrated the effectiveness of BMSCs co-culture with aging COCs to enhance oocyte quality and found its mechanism to function via ameliorating the biological function of GCs by alleviating FOXO3a levels. These results provide significant fundamental research on MSCs therapy on ovarian aging, as well as offering guidance for raising the success rate of assisted reproductive technology such IVM in clinical and non-clinical settings.

Transcriptome profiling reveals dysregulation of inflammatory and protein synthesis genes in PCOS

To analyze the differential expression genes of polycystic ovary syndrome (PCOS), clarify their functions and pathways, as well as the protein-protein interaction network, identify HUB genes, and explore the pathological mechanism. PCOS microarray datasets were screened from the GEO database. Common differentially expressed genes (co-DEGs) were obtained using GEO2R and Venn analysis. Enrichment and pathway analyses were conducted using the DAVID online tool, with results presented in bubble charts. Protein-protein interaction analysis was performed using the STRING tool. HUB genes were identified using Cytoscape software and further interpreted with the assistance of the GeneCards database. A total of two sets of co-DEGs (108 and 102), key proteins (15 and 55), and hub genes (10 and 10) were obtained. The co-DEGs: (1) regulated inflammatory responses and extracellular matrix, TNF, and IL-17 signaling pathways; (2) regulated ribosomes and protein translation, ribosome and immune pathways. The key proteins: (1) regulated inflammation, immunity, transcription, matrix metabolism, proliferation/differentiation, energy, and repair; (2) regulated ubiquitination, enzymes, companion proteins, respiratory chain components, and fusion proteins. The Hub genes: (1) encoded transcription factors and cytokines, playing vital roles in development and proliferation; (2) encoded ribosomes and protein synthesis, influencing hormone and protein synthesis, associated with development and infertility. The dysregulated expression of inflammation and protein synthesis genes in PCOS may be the key mechanism underlying its onset and progression.

Identification and Functional Analysis of circRNAs during Goat Follicular Development

Litter size is a crucial quantitative trait in animals, closely linked to follicular development. Circular RNA (circRNA), a type of single-stranded closed-loop endogenous RNA with stable expression, plays pivotal roles in various biological processes, yet its function in goat follicular development remains unclear. In this study, we collected large (follicle diameter > 3 mm) and small (1 mm < follicle diameter < 3 mm) follicles from black goats in the Chuanzhong region for circRNA sequencing, with the aim of elucidating the functional circRNAs that influence follicle development in goats. Differential analysis revealed that 17 circRNAs were upregulated in large follicles, and 28 circRNAs were upregulated in small follicles. Functional enrichment analysis revealed significant enrichment of pathways related to reproduction, including cellular response to follicle-stimulating hormone stimulus, the PI3K-Akt signaling pathway, the MAPK signaling pathway, and the Notch signaling pathway. Based on the ceRNA mechanism, 45 differentially expressed circRNAs were found to target and bind a total of 418 miRNAs, and an intercalation network including miR-324-3p (circRNA2497, circRNA5650), miR-202-5p (circRNA3333, circRNA5501), and miR-493-3p (circRNA4995, circRNA5508) was constructed. In addition, conservation analysis revealed that 2,239 circRNAs were conserved between goats and humans. Prediction of translation potential revealed that 154 circRNAs may potentially utilize both N6-methyladenosine (m6A) and internal ribosome entry site (IRES) translation mechanisms. Furthermore, the differential expression and circularization cleavage sites of five circRNAs were validated through RT-qPCR and DNA sequencing. Our study constructed a circRNA map in goat follicle development, offering a theoretical foundation for enhancing goat reproductive performance.

Transcriptome Dynamics and Cell Dialogs Between Oocytes and Granulosa Cells in Mouse Follicle Development

The development and maturation of follicles is a sophisticated and multistage process. The dynamic gene expression of oocytes and their surrounding somatic cells and the dialogs between these cells are critical to this process. In this study, we accurately classified the oocyte and follicle development into nine stages and profiled the gene expression of mouse oocytes and their surrounding granulosa cells and cumulus cells. The clustering of the transcriptomes showed the trajectories of two distinct development courses of oocytes and their surrounding somatic cells. Gene expression changes precipitously increased at Type 4 stage and drastically dropped afterward within both oocytes and granulosa cells. Moreover, the number of differentially expressed genes between oocytes and granulosa cells dramatically increased at Type 4 stage, most of which persistently passed on to the later stages. Strikingly, cell communications within and between oocytes and granulosa cells became active from Type 4 stage onward. Cell dialogs connected oocytes and granulosa cells in both unidirectional and bidirectional manners. TGFB2/3, TGFBR2/3, INHBA/B, and ACVR1/1B/2B of TGF-β signaling pathway functioned in the follicle development. NOTCH signaling pathway regulated the development of granulosa cells. Additionally, many maternally DNA methylation- or H3K27me3-imprinted genes remained active in granulosa cells but silent in oocytes during oogenesis. Collectively, Type 4 stage is the key turning point when significant transcription changes diverge the fate of oocytes and granulosa cells, and the cell dialogs become active to assure follicle development. These findings shed new insights on the transcriptome dynamics and cell dialogs facilitating the development and maturation of oocytes and follicles.

Deleterious variants in RNF111 impair female fertility and induce premature ovarian insufficiency in humans and mice

Premature ovarian insufficiency (POI) is a heterogeneous female disorder characterized by the loss of ovarian function before the age of 40. It represents a significant detriment to female fertility. However, the known POI-causative genes currently account for only a fraction of cases. To elucidate the genetic factors underlying POI, we conducted whole-exome sequencing on a family with three fertile POI patients and identified a deleterious missense variant in RNF111. In a subsequent replication study involving 1,030 POI patients, this variant was not only confirmed but also accompanied by the discovery of three additional predicted deleterious RNF111 variants. These variants collectively account for eight cases, representing 0.78% of the study cohort. A further study involving 500 patients with diminished ovarian reserve also identified two additional RNF111 variants. Notably, RNF111 encodes an E3 ubiquitin ligase with a regulatory role in the TGF-β/BMP signaling pathway. Our analysis revealed that RNF111/RNF111 is predominantly expressed in the oocytes of mice, monkeys, and humans. To further investigate the functional implications of RNF111 variants, we generated two mouse models: one with a heterozygous missense mutation (Rnf111+/M) and another with a heterozygous null mutation (Rnf111+/-). Both mouse models exhibited impaired female fertility, characterized by reduced litter sizes and small ovarian reserve. Additionally, RNA-seq and quantitative proteomics analysis unveiled that Rnf111 haploinsufficiency led to dysregulation in female gonad development and negative regulation of the BMP signaling pathway within mouse ovaries. In conclusion, our findings strongly suggest that monoallelic deleterious variants in RNF111 can impair female fertility and induce POI in both humans and mice.

Sheep litter size heredity basis using genome-wide selective analysis

Sheep are important herbivorous domestic animal globally, and the Chinese indigenous sheep breed has a multitude of economically significant variations due to the diverse geographical and ecological conditions. In particular, certain native breeds exhibit a visible high litter size phenotype due to the selection pressure of natural and artificial for thousands of years, offering an ideal animal model for investigating sheep's fecundity. In this study, selective signal analysis was performed on public whole-genome sequencing data from 60 sheep across eight breeds to identify candidate genes related to litter size. Results revealed that a total of 34,065,017 single-nucleotide polymorphisms (SNPs) were identified from all sheep, and 65 candidate genes (CDGs) were pinpointed from the top 1% of interacted windows and SNPs between the pairwise fixation index (FST, >0.149543) and cross-population extended haplotype homozygosity (XP-EHH, >0.701551). A total of 41 CDGs (e.g. VRTN, EYA2 and MCPH1) were annotated to 576 GO terms, of which seven terms were directly linked to follicular and embryonic development (e.g. TBXT, BMPR1B, and BMP2). In addition, 73 KEGG pathways were enriched by 21 CDGs (e.g. ENTPD5, ABCD4 and RXFP2), mainly related to Hippo (TCF4, BMPR1B and BMP2), TGF-β (BMPR1B and BMP2), PI3K-Akt (ITGB4, IL4R and PPP2R5A) and Jak-STAT signalling pathways (IL20RA and IL4R). Notably, a series of CDGs was under strong selection in sheep with high litter size traits. These findings result could improve the comprehension of the genetic underpinnings of sheep litter size. Furthermore, it provides valuable CDGS for future molecular breeding.

MCM8 interacts with DDX5 to promote R-loop resolution

MCM8 has emerged as a core gene in reproductive aging and is crucial for meiotic homologous recombination repair. It also safeguards genome stability by coordinating the replication stress response during mitosis, but its function in mitotic germ cells remains elusive. Here we found that disabling MCM8 in mice resulted in proliferation defects of primordial germ cells (PGCs) and ultimately impaired fertility. We further demonstrated that MCM8 interacted with two known helicases DDX5 and DHX9, and loss of MCM8 led to R-loop accumulation by reducing the retention of these helicases at R-loops, thus inducing genome instability. Cells expressing premature ovarian insufficiency-causative mutants of MCM8 with decreased interaction with DDX5 displayed increased R-loop levels. These results show MCM8 interacts with R-loop-resolving factors to prevent R-loop-induced DNA damage, which may contribute to the maintenance of genome integrity of PGCs and reproductive reserve establishment. Our findings thus reveal an essential role for MCM8 in PGC development and improve our understanding of reproductive aging caused by genome instability in mitotic germ cells.

The retrotransposon - derived capsid genes PNMA1 and PNMA4 maintain reproductive capacity

The human genome contains 24 gag -like capsid genes derived from deactivated retrotransposons conserved among eutherians. Although some of their encoded proteins retain the ability to form capsids and even transfer cargo, their fitness benefit has remained elusive. Here we show that the gag -like genes PNMA1 and PNMA4 support reproductive capacity. Six-week-old mice lacking either Pnma1 or Pnma4 are indistinguishable from wild-type littermates, but by six months the mutant mice become prematurely subfertile, with precipitous drops in sex hormone levels, gonadal atrophy, and abdominal obesity; overall they produce markedly fewer offspring than controls. Analysis of donated human ovaries shows that expression of both genes declines normally with aging, while several PNMA1 and PNMA4 variants identified in genome-wide association studies are causally associated with low testosterone, altered puberty onset, or obesity. These findings expand our understanding of factors that maintain human reproductive health and lend insight into the domestication of retrotransposon-derived genes.

The landscape of transcriptional profiles in human oocytes with different chromatin configurations

With increasingly used assisted reproductive technology (ART), the acquisition of high-quality oocytes and early embryos has become the focus of much attention. Studies in mice have found that the transition of chromatin conformation from non-surrounded nucleolus (NSN) to surrounded nucleolus (SN) is essential for oocyte maturation and early embryo development, and similar chromatin transition also exists in human oocytes. In this study, we collected human NSN and SN oocytes and investigated their transcriptome. The analysis of differentially expressed genes showed that epigenetic functions, cyclin-dependent kinases and transposable elements may play important roles in chromatin transition during human oocyte maturation. Our findings provide new insights into the molecular mechanism of NSN-to-SN transition of human oocyte and obtained new clues for improvement of oocyte in vitro maturation technique.

FOXL2 and NR5A1 induce human fibroblasts into steroidogenic ovarian granulosa-like cells

Human granulosa cells in different stages are essential for maintaining normal ovarian function, and granulosa cell defect is the main cause of ovarian dysfunction. To address this problem, it is necessary to induce functional granulosa cells at different stages in vitro. In this study, we established a reprogramming method to induce early- and late-stage granulosa cells with different steroidogenic abilities. We used an AMH-fluorescence-reporter system to screen candidate factors for cellular reprogramming and generated human induced granulosa-like cells (hiGC) by overexpressing FOXL2 and NR5A1. AMH-EGFP+ hiGC resembled human cumulus cells in transcriptome profiling and secreted high levels of oestrogen and progesterone, similar to late-stage granulosa cells at antral or preovulatory stage. Moreover, we identified CD55 as a cell surface marker that can be used to isolate early-stage granulosa cells. CD55+ AMH-EGFP- hiGC secreted high levels of oestrogen but low levels of progesterone, and their transcriptome profiles were more similar to early-stage granulosa cells. More importantly, CD55+ hiGC transplantation alleviated polycystic ovary syndrome (PCOS) in a mouse model. Therefore, hiGC provides a cellular model to study the developmental program of human granulosa cells and has potential to treat PCOS.

Reactive oxygen species signalling in the deterioration of quality of mammalian oocytes cultured in vitro: Protective effect of antioxidants

The in vitro fertilization (IVF) is the first choice of infertile couples worldwide to plan for conception. Besides having a significant advancement in IVF procedure, the success rate is still poor. Although several approaches have been tested to improve IVF protocol, minor changes in culture conditions, physical factors and/or drug treatment generate reactive oxygen species (ROS) in oocytes. Due to large size and huge number of mitochondria, oocyte is more susceptible towards ROS-mediated signalling under in vitro culture conditions. Elevation of ROS levels destabilize maturation promoting factor (MPF) that results in meiotic exit from diplotene as well as metaphase-II (M-II) arrest in vitro. Once meiotic exit occurs, these oocytes get further arrested at metaphase-I (M-I) stage or metaphase-III (M-III)-like stage under in vitro culture conditions. The M-I as well as M-III arrested oocytes are not fit for fertilization and limits IVF outcome. Further, the generation of excess levels of ROS cause oxidative stress (OS) that initiate downstream signalling to initiate various death pathways such as apoptosis, autophagy, necroptosis and deteriorates oocyte quality under in vitro culture conditions. The increase of cellular enzymatic antioxidants and/or supplementation of exogenous antioxidants in culture medium protect ROS-induced deterioration of oocyte quality in vitro. Although a growing body of evidence suggests the ROS and OS-mediated deterioration of oocyte quality in vitro, their downstream signalling and related mechanisms remain poorly understood. Hence, this review article summarizes the existing evidences concerning ROS and OS-mediated downstream signalling during deterioration of oocyte quality in vitro. The use of various antioxidants against ROS and OS-mediated impairment of oocyte quality in vitro has also been explored in order to increase the success rate of IVF during assisted reproductive health management.

In Vitro Growth of Mammalian Follicles and Oocytes

Mammalian ovaries contain a large number of immature follicles, most of which are destined to degenerate before ovulation [...].

GDF9His209GlnfsTer6/S428T and GDF9Q321X/S428T bi-allelic variants caused female subfertility with defective follicle enlargement

BACKGROUND

Antral follicles consist of an oocyte cumulus complex surrounding by somatic cells, including mural granulosa cells as the inner layer and theca cells as the outsider layer. The communications between oocytes and granulosa cells have been extensively explored in in vitro studies, however, the role of oocyte-derived factor GDF9 on in vivo antral follicle development remains elusive due to lack of an appropriate animal model. Clinically, the phenotype of GDF9 variants needs to be determined.

METHODS

Whole-exome sequencing (WES) was performed on two unrelated infertile women characterized by an early rise of estradiol level and defect in follicle enlargement. Besides, WES data on 1,039 women undergoing ART treatment were collected. A Gdf9Q308X/S415T mouse model was generated based on the variant found in one of the patients.

RESULTS

Two probands with bi-allelic GDF9 variants (GDF9His209GlnfsTer6/S428T, GDF9Q321X/S428T) and eight GDF9S428T heterozygotes with normal ovarian response were identified. In vitro experiments confirmed that these variants caused reduction of GDF9 secretion, and/or alleviation in BMP15 binding. Gdf9Q308X/S415T mouse model was constructed, which recapitulated the phenotypes in probands with abnormal estrogen secretion and defected follicle enlargement. Further experiments in mouse model showed an earlier expression of STAR in small antral follicles and decreased proliferative capacity in large antral follicles. In addition, RNA sequencing of granulosa cells revealed the transcriptomic profiles related to defective follicle enlargement in the Gdf9Q308X/S415T group. One of the downregulated genes, P4HA2 (a collagen related gene), was found to be stimulated by GDF9 protein, which partly explained the phenotype of defective follicle enlargement.

CONCLUSIONS

GDF9 bi-allelic variants contributed to the defect in antral follicle development. Oocyte itself participated in the regulation of follicle development through GDF9 paracrine effect, highlighting the essential role of oocyte-derived factors on ovarian response.

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.

CRL4DCAF13 E3 ubiquitin ligase targets MeCP2 for degradation to prevent DNA hypermethylation and ensure normal transcription in growing oocytes

The DNA methylation is gradually acquired during oogenesis, a process sustained by successful follicle development. However, the functional roles of methyl-CpG-binding protein 2 (MeCP2), an epigenetic regulator displaying specifical binding with methylated DNA, remains unknown in oogenesis. In this study, we found MeCP2 protein was highly expressed in primordial and primary follicle, but was almost undetectable in secondary follicles. However, in aged ovary, MeCP2 protein is significantly increased in both oocyte and granulosa cells. Overexpression of MeCP2 in growing oocyte caused transcription dysregulation, DNA hypermethylation, and genome instability, ultimately leading to follicle growth arrest and apoptosis. MeCP2 is targeted by DCAF13, a substrate recognition adaptor of the Cullin 4-RING (CRL4) E3 ligase, and polyubiquitinated for degradation in both cells and oocytes. Dcaf13-null oocyte exhibited an accumulation of MeCP2 protein, and the partial rescue of follicle growth arrest induced by Dcaf13 deletion was observed following MeCP2 knockdown. The RNA-seq results revealed that large amounts of genes were regulated by the DCAF13-MeCP2 axis in growing oocytes. Our study demonstrated that CRL4DCAF13 E3 ubiquitin ligase targets MeCP2 for degradation to ensure normal DNA methylome and transcription in growing oocytes. Moreover, in aged ovarian follicles, deceased DCAF13 and DDB1 protein were observed, indicating a potential novel mechanism that regulates ovary aging.

Validation of mitochondrial biomarkers and immune dynamics in polycystic ovary syndrome

PROBLEM

Polycystic ovary syndrome (PCOS), a prevalent endocrine-metabolic disorder, presents considerable therapeutic challenges due to its complex and elusive pathophysiology.

METHOD OF STUDY

We employed three machine learning algorithms to identify potential biomarkers within a training dataset, comprising GSE138518, GSE155489, and GSE193123. The diagnostic accuracy of these biomarkers was rigorously evaluated using a validation dataset using area under the curve (AUC) metrics. Further validation in clinical samples was conducted using PCR and immunofluorescence techniques. Additionally, we investigate the complex interplay among immune cells in PCOS using CIBERSORT to uncover the relationships between the identified biomarkers and various immune cell types.

RESULTS

Our analysis identified ACSS2, LPIN1, and NR4A1 as key mitochondria-related biomarkers associated with PCOS. A notable difference was observed in the immune microenvironment between PCOS patients and healthy controls. In particular, LPIN1 exhibited a positive correlation with resting mast cells, whereas NR4A1 demonstrated a negative correlation with monocytes in PCOS patients.

CONCLUSION

ACSS2, LPIN1, and NR4A1 emerge as PCOS-related diagnostic biomarkers and potential intervention targets, opening new avenues for the diagnosis and management of PCOS.

Spatiotemporal transcriptomic changes of human ovarian aging and the regulatory role of FOXP1

Limited understanding exists regarding how aging impacts the cellular and molecular aspects of the human ovary. This study combines single-cell RNA sequencing and spatial transcriptomics to systematically characterize human ovarian aging. Spatiotemporal molecular signatures of the eight types of ovarian cells during aging are observed. An analysis of age-associated changes in gene expression reveals that DNA damage response may be a key biological pathway in oocyte aging. Three granulosa cells subtypes and five theca and stromal cells subtypes, as well as their spatiotemporal transcriptomics changes during aging, are identified. FOXP1 emerges as a regulator of ovarian aging, declining with age and inhibiting CDKN1A transcription. Silencing FOXP1 results in premature ovarian insufficiency in mice. These findings offer a comprehensive understanding of spatiotemporal variability in human ovarian aging, aiding the prioritization of potential diagnostic biomarkers and therapeutic strategies.

A case report of a family with developmental arrest of human prokaryotic stage zygote

To study the genetic variation leading to the arrest phenotype of pronuclear (PN) zygotes. We recruited a family characterized by recurrent PN arrest during in vitro fertilization (IVF) and intracytoplasmic sperm injection cycles (ICSI) and performed whole-exome sequencing for 2 individuals. The transcriptome profiles of PN-arrest zygotes were assessed by single-cell RNA sequencing analysis. The variants were then validated by PCR amplification and Sanger sequencing in the affected individuals and other family members. A family characterized by recurrent PN arrest during IVF and ICSI cycles were enrolled after giving written informed consent. Peripheral blood samples were taken for DNA extraction. Three PN-arrest zygotes from patient III-3 were used for single-cell RNA-seq as described. This phenotype was reproduced after multiple cycles of egg retrieval and after trying different fertilization methods and multiple ovulation regimens. The mutant genes of whole exon sequencing were screened and verified. The missense variant c. C1630T (p.R544W) in RGS12 was responsible for a phenotype characterized by paternal transmission. RGS12 controls Ca2+ oscillation, which is required for oocyte activation after fertilization. Single-cell transcriptome profiling of PN-arrest zygotes revealed defective established translation, RNA processing and cell cycle, which explained the failure of complete oocyte activation. Furthermore, we identified proximal genes involved in Ca2+ oscillation-cytostatic factor-anaphase-promoting complex (Ca2+ oscillation-CSF-APC) signaling, including upregulated CaMKII, ORAI1, CDC20, and CDH1 and downregulated EMI1 and BUB3. The findings indicate abnormal spontaneous Ca2+ oscillations leading to oocytes with prolonged low CSF level and high APC level, which resulted in defective nuclear envelope breakdown and DNA replication. We have identified an RGS12 variant as the potential cause of female infertility characterized by arrest at the PN stage during multiple IVF and ICSI.

Microfluidic platform for omics analysis on single cells with diverse morphology and size: A review

BACKGROUND

Microfluidic techniques have emerged as powerful tools in single-cell research, facilitating the exploration of omics information from individual cells. Cell morphology is crucial for gene expression and physiological processes. However, there is currently a lack of integrated analysis of morphology and single-cell omics information. A critical challenge remains: what platform technologies are the best option to decode omics data of cells that are complex in morphology and size?

RESULTS

This review highlights achievements in microfluidic-based single-cell omics and isolation of cells based on morphology, along with other cell sorting methods based on physical characteristics. Various microfluidic platforms for single-cell isolation are systematically presented, showcasing their diversity and adaptability. The discussion focuses on microfluidic devices tailored to the distinct single-cell isolation requirements in plants and animals, emphasizing the significance of considering cell morphology and cell size in optimizing single-cell omics strategies. Simultaneously, it explores the application of microfluidic single-cell sorting technologies to single-cell sequencing, aiming to effectively integrate information about cell shape and size.

SIGNIFICANCE AND NOVELTY

The novelty lies in presenting a comprehensive overview of recent accomplishments in microfluidic-based single-cell omics, emphasizing the integration of different microfluidic platforms and their implications for cell morphology-based isolation. By underscoring the pivotal role of the specialized morphology of different cells in single-cell research, this review provides robust support for delving deeper into the exploration of single-cell omics data.

Unveiling the Ovarian Cell Characteristics and Molecular Mechanism of Prolificacy in Goats via Single-Nucleus Transcriptomics Data Analysis

Increases in litter size, which are influenced by ovulation, are responsible for between 74% and 96% of the economic value of genetic progress, which influences selection. For the selection and breeding of highly prolific goats, genetic mechanisms underlying variations in litter size should be elucidated. Here, we used single-nucleus RNA sequencing to analyze 44,605 single nuclei from the ovaries of polytocous and monotocous goats during the follicular phase. Utilizing known reference marker genes, we identified 10 ovarian cell types characterized by distinct gene expression profiles, transcription factor networks, and reciprocal interaction signatures. An in-depth analysis of the granulosa cells revealed three subtypes exhibiting distinct gene expression patterns and dynamic regulatory mechanisms. Further investigation of cell-type-specific prolificacy-associated transcriptional changes elucidated that "downregulation of apoptosis", "increased anabolism", and "upstream responsiveness to hormonal stimulation" are associated with prolificacy. This study provides a comprehensive understanding of the cell-type-specific mechanisms and regulatory networks in the goat ovary, providing insights into the molecular mechanisms underlying goat prolificacy. These findings establish a vital foundation for furthering understanding of the molecular mechanisms governing folliculogenesis and for improving the litter size in goats via molecular design breeding.

Integrated analysis reveals the regulatory mechanism of the neddylation inhibitor MLN4924 on the metabolic dysregulation in rabbit granulosa cells

BACKGROUND

Neddylation, an important post-translational modification (PTM) of proteins, plays a crucial role in follicular development. MLN4924 is a small-molecule inhibitor of the neddylation-activating enzyme (NAE) that regulates various biological processes. However, the regulatory mechanisms of neddylation in rabbit ovarian cells have not been emphasized. Here, the transcriptome and metabolome profiles in granulosa cells (GCs) treated with MLN4924 were utilized to identify differentially expressed genes, followed by pathway analysis to precisely define the altered metabolisms.

RESULTS

The results showed that 563 upregulated and 910 downregulated differentially expressed genes (DEGs) were mainly enriched in pathways related to cancer, cell cycle, PI3K-AKT, progesterone-mediated oocyte maturation, and PPAR signaling pathway. Furthermore, we characterized that MLN4924 inhibits PPAR-mediated lipid metabolism, and disrupts the cell cycle by promoting the apoptosis and proliferation of GCs. Importantly, we found the reduction of several metabolites in the MLN4924 treated GCs, including glycerophosphocholine, arachidic acid, and palmitic acid, which was consistent with the deregulation of PPAR signaling pathways. Furthermore, the increased metabolites included 6-Deoxy-6-sulfo-D-glucono-1,5-lactone and N-Acetyl-D-glucosaminyldiphosphodolichol. Combined with transcriptome data analyses, we identified genes that strongly correlate with metabolic dysregulation, particularly those related to glucose and lipid metabolism. Therefore, neddylation inhibition may disrupt the energy metabolism of GCs.

CONCLUSIONS

These results provide a foundation for in-depth research into the role and molecular mechanism of neddylation in ovary development.

Multiomics insights into the female reproductive aging

The human female reproductive lifespan significantly diminishes with age, leading to decreased fertility, reduced fertility quality and endocrine function disorders. While many aspects of aging in general have been extensively documented, the precise mechanisms governing programmed aging in the female reproductive system remain elusive. Recent advancements in omics technologies and computational capabilities have facilitated the emergence of multiomics deep phenotyping. Through the application and refinement of various high-throughput omics methods, a substantial volume of omics data has been generated, deepening our comprehension of the pathogenesis and molecular underpinnings of reproductive aging. This review highlights current and emerging multiomics approaches for investigating female reproductive aging, encompassing genomics, epigenomics, transcriptomics, proteomics, metabolomics, and microbiomics. We elucidate their influence on fundamental cell biology and translational research in the context of reproductive aging, address the limitations and current challenges associated with multiomics studies, and offer a glimpse into future prospects.

Multiomics Analyses Provide New Insight into Genetic Variation of Reproductive Adaptability in Tibetan Sheep

Domestication and artificial selection during production-oriented breeding have greatly shaped the level of genomic variability in sheep. However, the genetic variation associated with increased reproduction remains elusive. Here, two groups of samples from consecutively monotocous and polytocous sheep were collected for genome-wide association, transcriptomic, proteomic, and metabolomic analyses to explore the genetic variation in fecundity in Tibetan sheep. Genome-wide association study revealed strong associations between BMPR1B (p.Q249R) and litter size, as well as between PAPPA and lambing interval; these findings were validated in 1,130 individuals. Furthermore, we constructed the first single-cell atlas of Tibetan sheep ovary tissues and identified a specific mural granulosa cell subtype with PAPPA-specific expression and differential expression of BMPR1B between the two groups. Bulk RNA-seq indicated that BMPR1B and PAPPA expressions were similar between the two groups of sheep. 3D protein structure prediction and coimmunoprecipitation analysis indicated that mutation and mutually exclusive exons of BMPR1B are the main mechanisms for prolific Tibetan sheep. We propose that PAPPA is a key gene for stimulating ovarian follicular growth and development, and steroidogenesis. Our work reveals the genetic variation in reproductive performance in Tibetan sheep, providing insights and valuable genetic resources for the discovery of genes and regulatory mechanisms that improve reproductive success.

Optimization of human recombinant granulocyte-colony stimulating factor supplementation during in vitro production of porcine embryos to improve the efficiency of resource utilization of poor-quality cumulus-oocyte complexes

Granulocyte colony-stimulating factor (G-CSF), a pleiotropic cytokine, is secreted by the reproductive tract. Furthermore, our previous study indicated that human recombinant G-CSF (hrG-CSF) supplementation during porcine oocyte in vitro maturation (IVM) or during embryo in vitro culture (IVC) improved their quality and development potential when using cumulus-oocyte complexes (COCs) with more than three cumulus cell layers (CCL >3). Thus, in this study, we investigate the optimal conditions of hrG-CSF supplementation throughout the in vitro production (IVP: IVM + IVC) system to improve the embryo production efficiency of "poor-quality (CCL ≤3)" oocytes. COCs were classified into two groups according to the number of CCL (>3 and ≤3) and embryonic viability was analyzed after treatment with hrG-CSF during IVC. The mRNA transcription levels of G-CSF in COCs were compared based on their type and the period of IVM. Finally, developmental capacity and quality were evaluated after treatment with hrG-CSF for different periods of IVP. No marked effects on the developmental potential of embryos when using CCL ≤3 type COCs were observed after supplementing hrG-CSF only during IVC. Moreover, the mRNA transcription level of G-CSF increased gradually with IVM culture time and was higher in CCL ≤3 COCs than in >3. Supplementing hrG-CSF only during the IVM period resulted in the best embryo developmental potential, while supplementing hrG-CSF during the IVP period resulted in the best quality embryos, reflected in the increased total cell number and decreased apoptotic nuclei index of blastocysts. These findings indicate that "poor-quality" COCs may have a greater demand for G-CSF than "good-quality", meanwhile hrG-CSF supplementation throughout IVP improves resource utilization efficiency in poor-quality COCs.

Dysregulation of ferroptosis-related genes in granulosa cells associates with impaired oocyte quality in polycystic ovary syndrome

Background

Poor oocyte quality remains one of the major challenges for polycystic ovary syndrome (PCOS) patients during in vitro fertilization (IVF) treatment. Granulosa cells (GCs) in PCOS display altered functions and could cause an unfavorable microenvironment for oocyte growth and maturation. Ferroptosis is a new form of programmed cell death, but its role in PCOS has been largely unclarified.

Methods

Ferroptosis-related differentially expressed genes (DEGs) of GCs in women with PCOS were identified by bioinformatic analyses of GSE155489 and GSE168404 datasets. Functional enrichment analyses were conducted using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes. Core ferroptosis-related genes were further screened by random forest, and evaluated for diagnostic value by receiver operating characteristic curve analyses. Gene expression was validated by real-time quantitative polymerase chain reaction of collected GC samples, and analyzed for association with oocyte quality. In addition, gene regulatory network was constructed based on predicted RNA interactions and transcription factors, while potential therapeutic compounds were screened through molecular docking with crystallographic protein structures.

Results

A total of 14 ferroptosis-related DEGs were identified. These DEGs were mainly enriched in reactive oxygen species metabolic process, mitochondrial outer membrane, antioxidant activity as well as ferroptosis and adipocytokine signaling pathways. Eight core ferroptosis-related genes (ATF3, BNIP3, DDIT4, LPIN1, NOS2, NQO1, SLC2A1 and SLC2A6) were further selected in random forest model, which showed high diagnostic performance for PCOS. Seven of them were validated in GC samples, and five were found to be significantly and positively correlated with one or more oocyte quality parameters in PCOS patients, including oocyte retrieval rate, mature oocyte rate, normal fertilization rate, and good-quality embryo rate. Gene regulatory network revealed JUN and HMGA1 as two important transcription factors, while dicoumarol and flavin adenine dinucleotide were predicted as small molecules with therapeutic potential.

Conclusions

This is the first comprehensive report to study the differential expression of ferroptosis-related genes in GCs of PCOS and their clinical relevance with oocyte quality. Our findings could provide novel insights on the potential role of GC ferroptosis in PCOS pathogenesis, diagnosis, and targeted treatment.

miR-6881-3p contributes to diminished ovarian reserve by regulating granulosa cell apoptosis by targeting SMAD4

BACKGROUND

In our previous investigation, we revealed a significant increase in the expression of microRNA-6881-3p (miR-6881-3p) in follicular fluid granulosa cells (GCs) from women with diminished ovarian reserve (DOR) compared to those with normal ovarian reserve (NOR). However, the role of miR-6881-3p in the development of DOR remains poorly understood.

OBJECTIVE

This study aimed to elucidate the involvement of miR-6881-3p in the regulation of granulosa cells (GCs) function and the pathogenesis of DOR.

MATERIALS AND METHODS

Initially, we assessed the expression levels of miR-6881-3p in GCs obtained from human follicular fluid in both NOR and DOR cases and explored the correlation between miR-6881-3p expression and clinical outcomes in assisted reproduction technology (ART). Bioinformatic predictions and dual-luciferase reporter assays were employed to identify the target gene of miR-6881-3p. Manipulation of miR-6881-3p expression was achieved through the transfection of KGN cells with miR-6881-3p mimics, inhibitor, and miRNA negative control (NC). Following transfection, we assessed granulosa cell apoptosis and cell cycle progression via flow cytometry and quantified target gene expression through quantitative real-time polymerase chain reaction (qRT-PCR) and Western blot (WB) analysis. Finally, we examined the correlation between target gene expression levels in GCs from NOR and DOR patients and their association with ART outcomes.

RESULTS

Our findings revealed elevated miR-6881-3p levels in GCs from DOR patients, which negatively correlated with ovarian reserve function and ART outcomes. We identified a direct binding interaction between miR-6881-3p and the 3'-untranslated region of the SMAD4. Transfection with miR-6881-3p mimics induced apoptosis in KGN cell. Furthermore, miR-6881-3p expression negatively correlated with both mRNA and protein levels of the SMAD4. The mRNA and protein levels of SMAD4 were notably reduced in GCs from DOR patients, and SMAD4 mRNA expression positively correlated with ART outcomes. In addition, the mRNA levels of FSHR, CYP11A1 were notably reduced after transfection with miR-6881-3p mimics in KGN cell, while LHCGR notably increased. The mRNA and protein levels of FSHR, CYP11A1 were notably reduced in GCs from DOR patients, while LHCGR notably increased.

CONCLUSION

This study underscores the role of miR-6881-3p in directly targeting SMAD4 mRNA, subsequently diminishing granulosa cell viability and promoting apoptosis, and may affect steroid hormone regulation and gonadotropin signal reception in GCs. These findings contribute to our understanding of the pathogenesis of DOR.

Extracellular matrix-templating fibrous hydrogels promote ovarian tissue remodeling and oocyte growth

Synthetic matrices which mimic the extracellular composition of native tissue create a comprehensive model for studying development and disease. Here, we have engineered a composite material which retains cell-secreted ECM for the culture of ovarian follicles by embedding electrospun dextran fibers functionalized with basement membrane binder (BMB) peptide in PEG hydrogels. In the presence of ECM-sequestering fibers, encapsulated immature primordial follicles and ovarian stromal cells aggregated into large organoid-like structures with dense deposition of laminin, perlecan, and collagen I, leading to steroidogenesis and significantly greater rates of oocyte survival and growth. We determined that cell aggregation restored key cell-cell interactions critical for oocyte survival, whereas oocyte growth was dependent on cell-matrix interactions achieved in the presence of BMB. Here we have shown that sequestration and retention of cell-secreted ECM along synthetic fibers mimics fibrous ECM structure and restores the cell-cell and cell-matrix interactions critical for engineering an artificial ovary.

Global crotonylome identifies EP300-regulated ANXA2 crotonylation in cumulus cells as a regulator of oocyte maturation

Lysine crotonylation (Kcr), a newly discovered post-translational modification, played a crucial role in physiology and disease progression. However, the roles of crotonylation in oocyte meiotic resumption remain elusive. As abnormal cumulus cell development will cause oocyte maturation arrest and female infertility, we report that cumulus cells surrounding human meiotic arrested oocytes showed significantly lower crotonylation, which was associated with decreased EP300 expression and blocked cumulus cell expansion. In cultured human cumulus cells, exogenous crotonylation or EP300 activator promoted cell proliferation and reduced cell apoptosis, whereas EP300 knockdown induced the opposite effect. Transcriptome profiling analysis in human cumulus cells indicated that functions of crotonylation were associated with activation of epidermal growth factor receptor (EGFR) pathway. Importantly, we characterized the Kcr proteomics landscape in cumulus cells by LC-MS/MS analysis, and identified that annexin A2 (ANXA2) was crotonylated in cumulus cells in an EP300-dependent manner. Crotonylation of ANXA2 enhanced the ANXA2-EGFR binding, and then activated the EGFR pathway to affect cumulus cell proliferation and apoptosis. Using mouse oocytes IVM model and EP300 knockout mice, we further confirmed that crotonylation alteration in cumulus cells affected the oocyte maturation. Together, our results indicated that EP300-mediated crotonylation is important for cumulus cells functions and oocyte maturation.

Addition of rapamycin or co-culture with cumulus cells from younger reproductive age women does not improve rescue in vitro oocyte maturation or euploidy rates in older reproductive age women

Both spontaneously conceived pregnancies and those achieved using assisted reproduction decline with advancing maternal age. In this study, we tested if rapamycin and/or cumulus cells (CCs) from young donors could improve oocyte maturation and euploidy rates of germinal vesicle (GV) stage oocytes obtained from older women of reproductive age. A total of 498 GVs from 201 women >38 years (40.6 ± 1.8, mean ± SD) were included. GVs were randomly assigned into five groups for rescue IVM: control (with no CCs and no rapamycin); with autologous CCs; with autologous CCs and rapamycin; with CCs from young women (<35 years); and with CCs from young women and rapamycin. After 24 h of culture, the first polar body (PB) was biopsied in metaphase II oocytes, and the cytogenetic constitution was assessed using next-generation sequencing for both oocytes and PBs. Comparable maturation rates were found (56.2%, 60.0%, 46.5%, 51.7%, and 48.5% for groups 1-5, respectively; P = 0.30). Similarly, comparable euploidy rates were observed in the five groups (41.5%, 37.8%, 47.2%, 43.6%, and 47.8% for Groups 1-5, respectively; P = 0.87). Our findings indicate that rescue IVM is effective for obtaining mature euploid oocytes in older women of reproductive age, and that incubation with rapamycin or CCs obtained from young donors does not improve the maturation or euploidy rate.

Identification and molecular characterization of genes modulating progression of an oocyte from M-I to M-II in rat ovary

BACKGROUND

To achieve oocyte competence for successful fertilization, bidirectional communication between oocyte and granulosa cells is crucial. The acquisition of meiotic competency in oocyte is facilitated by various regulatory genes however, expression pattern of these genes is not well documented during meiotic transition from Metaphase-I to Metaphase-II stage. Therefore, the present research analyzed the expression pattern of regulatory genes that are involved in the transition from M-I to M-II stages in rat oocyte.

METHODS

The analysis of the data was conducted by applying an array of bioinformatic tools. The investigation of gene group interactions was carried out by employing the STRING database, which relies on co-expression information. The gene ontology (GO) analysis was performed utilizing the comparative GO database. Functional annotation for GO and pathway enrichment analysis were performed for genes involved in networking. The GO obtained through computational simulations was subsequently validated using quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis.

RESULTS

The findings of our study suggest that there is a distinct gene expression pattern in both the oocyte and granulosa cells. This pattern indicates that oocyte-secreted factors, such as BMP15 and GDF9, play a crucial role in regulating the progression of the meiotic cell cycle from the M-I to M-II stages. We have also examined the level of mRNA expression of genes including CYP11A1, CYP19A1, and STAR, which are crucial for the steroidogenesis.

CONCLUSIONS

It is fascinating to observe that the oscillatory pattern of specific key genes may hold significance in the process of in vitro oocyte maturation, specifically during the transition from the M-I to M-II stage. It might be useful for determining biomarker genes and potential pathways that play a role in attaining oocyte competency, thereby aiding in the assessment of oocyte quality for the purpose of achieving successful fertilization.

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.

Spatiotemporal Distribution and Function of Mitochondria in Oocytes

Mitochondria are energy provider organelles in eukaryotic cells that contain their own specific genome. This review addresses structural and functional properties of mitochondria, focusing on recent discoveries about the changes in quality and number of mitochondria per cell during oocyte development. We highlight how oocyte mitochondria exhibit stage-specific morphology and characteristics at different stages of development, in sharp contrast to the elongated mitochondria present in somatic cells. We then evaluate the latest transcriptomic data to elucidate the complex functions of mitochondria during oocyte maturation and the impact of mitochondria on oocyte development. Finally, we describe the methodological progress of mitochondrial replacement therapy to rescue oocytes with developmental disorders or mitochondrial diseases, hoping to provide a guiding reference to future clinical applications.

Neurotrophin-4 promotes in vitro development and maturation of human secondary follicles yielding metaphase II oocytes and successful blastocyst formation

STUDY QUESTION

Does a matrix-free culture system supplemented with neurotrophic factor 4 (NT4) improve human in vitro follicular development and meiotic maturation, ultimately resulting in fertilizable oocytes?

SUMMARY ANSWER

NT4 supplementation of in vitro culture significantly enhances the growth, steroid hormone production, and maturity potential of human secondary follicles derived from fresh ovarian medulla (from post- and pre-pubertal patients), thereby yielding fertilizable oocytes.

WHAT IS KNOWN ALREADY

Reconstituting folliculogenesis in vitro is of paramount importance in the realms of fertility preservation, reproductive biology research, and reproductive toxicity assessments. However, the efficiency of in vitro culture systems remains suboptimal, as the attainment of fertilizable oocytes from in vitro growth (IVG) of human follicles remains unachieved, with the data being particularly scant regarding follicles from prepubertal girls. We have previously found that mouse oocytes from secondary follicles derived from IVG are deficient in neuroendocrine regulation. NT4 and its corresponding receptor have been identified in human follicles. Significantly, the addition of NT4 during the IVG process markedly enhances both follicle growth and oocyte maturation rates in mice.

STUDY DESIGN SIZE DURATION

Fresh medulla tissue obtained during tissue preparation for ovarian tissue cryopreservation (OTC) were collected from 10 patients aged from 6 to 21 years old, all of whom had undergone unilateral oophorectomy as a means of fertility preservation. Isolated secondary follicles were individually cultured in vitro with or without NT4 in a matrix-free system.

PARTICIPANTS/MATERIALS SETTING METHODS

Secondary follicles, extracted via enzymatic digestion and mechanical disruption from each patient, were randomly allocated to either a control group or an NT4-supplemented group (100 ng/ml), followed by individual culture on an ultra-low attachment plate. Follicle growth and viability were assessed by microscopy. Levels of anti-Müllerian hormone (AMH), estradiol, and progesterone in the medium were quantified. An oocyte-specific marker was identified using confocal fluorescence microscopy following DEAD box polypeptide 4 (DDX4) staining. The competence of individual oocytes for maturation and fertilization were assessed after IVM and ICSI with donated sperm samples.

MAIN RESULTS AND THE ROLE OF CHANCE

Overall, isolated follicles from both groups survived up to 6 weeks with increasing diameters over the duration (P < 0.05), reaching terminal diameters of almost 1 mm with confirmed steroidogenesis and expression of oocyte marker (DDX4), and producing morphologically normal MII oocytes. When compared with the control group, the NT4 group had a similar initial follicular diameter (206 ± 61.3 vs 184 ± 93.4 μm) but exhibited a significant increase in follicular diameter from the ninth day of culture onwards (P < 0.05). From Week 3, estradiol and progesterone production were significantly increased in the NT4 group, while no significant difference was observed in AMH production between groups. The proportion of 'fast-growth' follicles in the NT4 group was significantly higher than that in the control group (13/23 vs 6/24, P < 0.05). An increased efficiency of MII oocyte maturation per live follicle in the NT4 group was also observed (control group vs NT4 group, 4/24 vs 10/23, P < 0.05). It is noteworthy that an MII oocyte obtained from the control group exhibited abnormal fertilization after ICSI. In contrast, an MII oocyte acquired from the NT4 group progressed to the blastocyst stage and showed potential for transfer.

LARGE SCALE DATA

N/A.

LIMITATIONS REASONS FOR CAUTION

The cohort examined in this study was all patients diagnosed with beta-thalassemia major. Whether this culture system is effective for patients with other diseases remains unknown. Since the chosen dose of NT4 was established based on dose finding in mice, the optimal dose for use in a human IVG system needs further confirmation. The oocytes and embryos procured from this study have not been quantified for ploidy status or epigenetic signatures.

WIDER IMPLICATIONS OF THE FINDINGS

Fresh medulla tissue obtained during tissue preparation for OTC may serve as a precious source of fertilizable oocytes for female fertility preservation, even for pre-pubertal girls, without the threat of tumour reintroduction. After further characterization and optimization of the system, this culture system holds the potential to provide a powerful future research tool, for the comprehensive exploration of human follicular development mechanisms and for conducting reproductive toxicity evaluations.

STUDY FUNDING/COMPETING INTERESTS

This work was supported by the National Key R&D Program of China (grant number 2022YFC2703000) and National Natural Science Foundation of China (grant numbers 82271651 and 81871214). The medium used in human follicle in vitro culture in this study has been applied for a national invention patent in China (No. 202211330660.7). The inventors of the patent, in order, are: Y.G., C.F., and X.L.

Core fucosylation regulates the ovarian response via FSH receptor during follicular development

INTRODUCTION

Ovarian low response to follicle-stimulating hormone (FSH) causes infertility featuring hypergonadotropic hypogonadism, ovarian failure, and/or defective ovarian response.

OBJECTIVES

N-glycosylation is essential for FSH receptor (FSHR). Core fucosylation catalyzed by fucosyltransferase 8 (FUT8) is the most common N-glycosylation. Core fucosylation level changes between individuals and plays important roles in multiple physiological and pathological conditions. This study aims to elucidate the significance of FUT8 to modulate FSHR function in female fertility.

METHODS

Samples from patients classified as poor ovary responders (PORs) were detected with lectin blot and real-time PCR. Fut8 gene knockout (Fut8-/-) mice and FUT8-knockdown human granulosa cell line (KGN-KD) were established and in vitro fertilization (IVF) assay, western blot, molecular interaction, immunofluorescence and immunoprecipitation were applied.

RESULTS

Core fucosylation is indispensable for oocyte and follicular development. FSHR is a highly core-fucosylated glycoprotein. Loss of core fucosylation suppressed binding of FSHR to FSH, and attenuated FSHR downstream signaling in granulosa cells. Transcriptomic analysis revealed the downregulation of several transcripts crucial for oocyte meiotic progression and preimplantation development in Fut8-/- mice and in POR patients. Furthermore, loss of FUT8 inhibited the interaction between granulosa cells and oocytes, reduced transzonal projection (TZP) formation and caused poor developmental competence of oocytes after fertilization in vitro. While L-fucose administration increased the core fucosylation of FSHR, and its sensitivity to FSH.

CONCLUSION

This study first reveals a significant presence of core fucosylation in female fertility control. Decreased fucosylation on FSHR reduces the interaction of FSH-FSHR and subsequent signaling, which is a feature of the POR patients. Our results suggest that core fucosylation controls oocyte and follicular development via the FSH/FSHR pathway and is essential for female fertility in mammals.

Single-Cell Transcriptome Analysis Reveals Development-Specific Networks at Distinct Synchronized Antral Follicle Sizes in Sheep Oocytes

The development of the ovarian antral follicle is a complex, highly regulated process. Oocytes orchestrate and coordinate the development of mammalian ovarian follicles, and the rate of follicular development is governed by a developmental program intrinsic to the oocyte. Characterizing oocyte signatures during this dynamic process is critical for understanding oocyte maturation and follicular development. Although the transcriptional signature of sheep oocytes matured in vitro and preovulatory oocytes have been previously described, the transcriptional changes of oocytes in antral follicles have not. Here, we used single-cell transcriptomics (SmartSeq2) to characterize sheep oocytes from small, medium, and large antral follicles. We characterized the transcriptomic landscape of sheep oocytes during antral follicle development, identifying unique features in the transcriptional atlas, stage-specific molecular signatures, oocyte-secreted factors, and transcription factor networks. Notably, we identified the specific expression of 222 genes in the LO, 8 and 6 genes that were stage-specific in the MO and SO, respectively. We also elucidated signaling pathways in each antral follicle size that may reflect oocyte quality and in vitro maturation competency. Additionally, we discovered key biological processes that drive the transition from small to large antral follicles, revealing hub genes involved in follicle recruitment and selection. Thus, our work provides a comprehensive characterization of the single-oocyte transcriptome, filling a gap in the mapping of the molecular landscape of sheep oogenesis. We also provide key insights into the transcriptional regulation of the critical sizes of antral follicular development, which is essential for understanding how the oocyte orchestrates follicular development.