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

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

Follicular Hyperstimulation Dysgenesis: New Explanation for Adverse Effects of Excessive FSH in Ovarian Stimulation

High follicle-stimulating hormone (FSH) doses during ovarian stimulation protocols for assisted reproductive technologies (ART) are detrimental to ovulatory follicle function and oocyte quality. However, the mechanisms are unclear. In a small ovarian reserve heifer model, excessive FSH doses lead to phenotypic heterogeneity of ovulatory size follicles, with most follicles displaying signs of premature luteinization and a range in severity of abnormalities. By performing whole transcriptome analyses of granulosa cells, cumulus cells, and oocytes from individual follicles of animals given standard or excessive FSH doses, we identified progressive changes in the transcriptomes of the 3 cell types, with increasing severity of follicular abnormality with the excessive doses. The granulosa and cumulus cells each diverged progressively from their normal phenotypes and became highly similar to each other in the more severely affected follicles. Pathway analysis indicates a possible dysregulation of the final stages of folliculogenesis, with processes characteristic of ovulation and luteinization occurring concurrently rather than sequentially in the most severely affected follicles. These changes were associated with disruptions in key pathways in granulosa and cumulus cells, which may account for previously reported reduced estradiol production, enhanced progesterone and oxytocin production and diminished ovulation rates. Predicted deficiencies in oocyte survival, stress response, and fertilization suggest likely reductions in oocyte health, which could further compromise oocyte quality and ART outcomes.

IGF2R, KCNQ1, PLAGL1, and SNRPN DNA methylation is completed in bovine by the early antral follicle stage

Imprinted genes are inherited with different DNA methylation patterns depending on the maternal or paternal origin of the allele. In cattle (Bos taurus), abnormal methylation of these genes is linked to the large offspring syndrome, a neonatal overgrowth phenotype analogous to the human Beckwith-Wiedemann syndrome. We hypothesized that in bovine oocytes, some of the methylation patterns on maternally imprinted genes are acquired in the last phase of folliculogenesis. The pyrosequencing analysis of IGF2R, KCNQ1, PLAGL1, and SNRPN imprinted genes showed no clear progression of methylation in oocytes from follicles 1-2 mm (late pre antral/early antral) and up. Instead, these oocytes displayed complete methylation at the imprinted differentially methylated regions (>80%). Other mechanisms related to imprint maintenance should be investigated to explain the hypomethylation at IGF2R, KCNQ1, PLAGL1, and SNRPN maternally imprinted sites observed in some bovine embryos.

Integration of transcriptomics and non-targeted metabolomics reveals the underlying mechanism of follicular atresia in Chinese buffalo

Follicular atresia is a complex physiological process, which results in the waste of follicles and oocytes from the ovary. Elucidating the physiological mechanism of follicular atresia will hopefully reverse the fate of follicles, thereby improve the reproductive efficiency of female animals. However, there are still many gaps to be filled during the follicular atresia process. In this study, we first comprehensively summarized and compared a variety of methods to classify Chinese buffalo follicles with different extent of atresia. Then follicular fluid and granulosa cells from the corresponding follicles with different extent of atresia were collected for non-targeted metabolomics and transcriptomics analysis, respectively. After the detection and analysis of 129 follicles, a reasonable classification standard was formed: on the basis of morphological classification, the relative concentrations of estradiol (E2) and progesterone (PROG) in the follicular fluid were determined, follicles with an estradiol-to-progesterone (E2/PROG) ratio >5 were classified as healthy follicles (HF), 1≤ E2/PROG ≤5 as early atretic follicles (EF) and E2/PROG <1 as late atretic follicles (LF). Correspondingly, follicles with granulosa cells apoptosis rate less than 15 % were divided into HF, 15%-25% were classified as EF and more than 25 % were classified as LF. The integration analysis of non-targeted metabolomics and transcriptomics highlights the following three aspects: (1) Atresia seriously damaged the lipid metabolism homeostasis of follicle, in which PPARγ play important roles. (2) Energy metabolism and nucleotide metabolism of atretic follicles were inhibited. (3) Bilirubin is involved in follicular atresia, and it may be the main force to prevent lipid peroxidation in follicular cells. In summary, results of this study provide new understanding of the molecular mechanisms of Chinese buffalo follicular atresia.

Specific imprinted genes demethylation in association with oocyte donor's age and culture conditions in bovine embryos assessed at day 7 and 12 post insemination

The production of bovine embryos through in vitro maturation and fertilization is an important tool of the genomic revolution in dairy cattle. Gene expression analysis of these embryos revealed differences according to the culture conditions or oocyte donor's pubertal status compared to in vivo derived embryos. We hypothesized that some of the methylation patterns in oocytes are acquired in the last step of folliculogenesis and could be influenced by the environment created in the follicles containing these oocytes. These altered patterns may not be erased during the first week of embryonic development in culture or may be sensitive to the conditions during that time. To quantify the changes related to culture conditions, an in vivo control group consisting of embryos (Day 12 post fertilization for all groups) obtained from superovulated and artificially inseminated cows was compared to in vitro produced (IVP) embryos cultured with or without Fetal Bovine Serum (FBS). To measure the effect of the oocytes donor's age, we also compared a fourth group consisting of IVP embryos produced with oocytes collected following ovarian stimulation of pre-pubertal animals. Embryonic disk and trophoblast cells were processed separately and the methylation status of ten imprinted genes (H19, MEST, KCNQ1, SNRPN, PEG3, NNAT, GNASXL, IGF2R, PEG10, and PLAGL1) was assessed by pyrosequencing. Next, ten Day 7 blastocysts were produced following the same methodology as for the D12 embryos (four groups) to observe the most interesting genes (KCNQ1, SNRPN, IGF2R and PLAGL1) at an earlier developmental stage. For all samples, we observed overall lower methylation levels and greater variability in the three in vitro groups compared to the in vivo group. The individual embryo analysis indicated that some embryos were deviant from the others and some were not affected. We concluded that IGF2R, SNRPN, and PEG10 were particularly sensitive to culture conditions and the presence of FBS, while KCNQ1 and PLAGL1 were more affected in embryos derived from pre-pubertal donors. This work provides markers at the single imprinted control region (ICR) resolution to assess the culture environment required to minimize epigenetic perturbations in bovine embryos generated by assisted reproduction techniques, thus laying the groundwork for a better comprehension of the complex interplay between in vitro conditions and imprinted genes.

Effect of follicle size and atresia grade on mitochondrial membrane potential and steroidogenic acute regulatory protein expression in bovine granulosa cells

During follicular development, granulosa cells undergo functional and structural changes affecting their steroidogenic activity. Oestrogen synthesis mainly occurs in the endoplasmic reticulum and relies on aromatase activity to convert androgens that arise from theca cells. In the present study, indicators of mitochondria-related steroidogenic capacity, as steroidogenic acute regulatory (StAR) protein expression and mitochondrial membrane potential (MMP), have been evaluated in bovine granulosa cells (GCs) and related to follicle growth and atresia. Atresia was estimated by morphological examination of follicle walls and cumulus-oocyte complexes (COC) and assessed by terminal deoxynucleotidyl transferase (TdT)-mediated dUTP nick end labeling (TUNEL) assay for apoptosis detection. Bovine ovarian follicles were macroscopically classified according to their atresia grade and grouped into small, medium or large follicles. After follicle opening, the COCs were morphologically classified for follicle atresia and the GCs were collected. Granulosa cells were fixed for immunofluorescence (IF) and TUNEL assay, frozen for western blotting (WB) or freshly maintained for MMP analyses. StAR protein expression was assessed using both IF and WB analyses. The follicle atresia grade could be efficiently discriminated based on either follicle wall or COC morphological evaluations. Granulosa cells collected from small non-atretic follicles showed a higher (P <0.01) MMP and WB-based StAR protein expression than small atretic follicles. For IF analysis, StAR protein expression in large atretic follicles was higher (P <0.05) than that in large non-atretic follicles. These results suggest a role played by mitochondria in GC steroidogenic activity, which declines in healthy follicles along with their growth. In large follicles, steroidogenic activity increases with atresia and is possibly associated with progesterone production.