Single-cell RNA sequencing reveals regulation of fetal ovary development in the monkey (Macaca fascicularis)

Temporal metabolomics analysis reveals the metabolic patterns in goat cumulus cells during oocyte maturation

Cumulus cells play a crucial role in the oocyte growth and maturation processes through providing necessary nutrients and growth signals by gap junction communication. However, a global overview of metabolic events in goat cumulus cells is still lacking. In the present study, we collected cumulus cells from goat cumulus-oocyte complexes (COCs) at different developmental stages. Metabolomics analysis was performed to investigate the global metabolic patterns in cumulus cells during oocyte in vitro maturation. In particular, we revealed the several significantly altered metabolic pathways and metaboliccharacteristics in goat cumulus cells, including the accumulation of fatty acids, steroid hormones metabolism, active catabolism of arginine during meiotic resumption, and a progressive decline in nucleotide metabolism. In conclusion, the dataset generated by our metabolomic profiling will provide valuable information to understand the key metabolic pathways and metabolites involved in COCs development.

PICALM Regulating the Generation of Amyloid β-Peptide to Promote Anthracycline-Induced Cardiotoxicity

Anthracyclines are chemotherapeutic drugs used to treat solid and hematologic malignancies. However, life-threatening cardiotoxicity, with cardiac dilation and heart failure, is a drawback. A combination of in vivo for single cell/nucleus RNA sequencing and in vitro approaches is used to elucidate the underlying mechanism. Genetic depletion and pharmacological blocking peptides on phosphatidylinositol binding clathrin assembly (PICALM) are used to evaluate the role of PICALM in doxorubicin-induced cardiotoxicity in vivo. Human heart tissue samples are used for verification. Patients with end-stage heart failure and chemotherapy-induced cardiotoxicity have thinner cell membranes compared to healthy controls do. Using the doxorubicin-induced cardiotoxicity mice model, it is possible to replicate the corresponding phenotype in patients. Cellular changes in doxorubicin-induced cardiotoxicity in mice, especially in cardiomyocytes, are identified using single cell/nucleus RNA sequencing. Picalm expression is upregulated only in cardiomyocytes with doxorubicin-induced cardiotoxicity. Amyloid β-peptide production is also increased after doxorubicin treatment, which leads to a greater increase in the membrane permeability of cardiomyocytes. Genetic depletion and pharmacological blocking peptides on Picalm reduce the generation of amyloid β-peptide. This alleviates the doxorubicin-induced cardiotoxicity in vitro and in vivo. In human heart tissue samples of patients with chemotherapy-induced cardiotoxicity, PICALM, and amyloid β-peptide are elevated as well.

Comparative genomics of macaques and integrated insights into genetic variation and population history

The crab-eating macaques ( Macaca fascicularis ) and rhesus macaques ( M. mulatta ) are widely studied nonhuman primates in biomedical and evolutionary research. Despite their significance, the current understanding of the complex genomic structure in macaques and the differences between species requires substantial improvement. Here, we present a complete genome assembly of a crab-eating macaque and 20 haplotype-resolved macaque assemblies to investigate the complex regions and major genomic differences between species. Segmental duplication in macaques is ∼42% lower, while centromeres are ∼3.7 times longer than those in humans. The characterization of ∼2 Mbp fixed genetic variants and ∼240 Mbp complex loci highlights potential associations with metabolic differences between the two macaque species (e.g., CYP2C76 and EHBP1L1 ). Additionally, hundreds of alternative splicing differences show post-transcriptional regulation divergence between these two species (e.g., PNPO ). We also characterize 91 large-scale genomic differences between macaques and humans at a single-base-pair resolution and highlight their impact on gene regulation in primate evolution (e.g., FOLH1 and PIEZO2 ). Finally, population genetics recapitulates macaque speciation and selective sweeps, highlighting potential genetic basis of reproduction and tail phenotype differences (e.g., STAB1 , SEMA3F , and HOXD13 ). In summary, the integrated analysis of genetic variation and population genetics in macaques greatly enhances our comprehension of lineage-specific phenotypes, adaptation, and primate evolution, thereby improving their biomedical applications in human diseases.

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.

Roles of immune microenvironment in the female reproductive maintenance and regulation: novel insights into the crosstalk of immune cells

Female fertility decline is an accumulative consequence caused by complex factors, among them, the disruption of the immune profile in female reproduction stands out as a crucial contributor. Presently, the effects of immune microenvironment (IME) on the female reproductive process have attracted increasing attentions for their dynamic but precisive roles. Immunocytes including macrophages, dendritic cells, T cells, B cells and neutrophils, with diverse subpopulations as well as high plasticity functioned dynamically in the process of female reproduction through indirect intercellular communication via specific cytokine release transduced by molecular signal networks or direct cell-cell contact to maintain the stability of the reproductive process have been unveiled. The immune profile of female reproduction in each stage has also been meticulously unveiled. Especially, the application of single-cell sequencing (scRNA-seq) technology in this process reveals the distribution map of immune cells, which gives a novel insight for the homeostasis of IME and provides a research direction for better exploring the role of immune cells in female reproduction. Here, we provide an all-encompassing overview of the latest advancements in immune modulation within the context of the female reproductive process. Our approach involves structuring our summary in accordance with the physiological sequence encompassing gonadogenesis, folliculogenesis within the ovaries, ovulation through the fallopian tubes, and the subsequent stages of embryo implantation and development within the uterus. Our overarching objective is to construct a comprehensive portrayal of the immune microenvironment (IME), thereby accentuating the pivotal role played by immune cells in governing the intricate female reproductive journey. Additionally, we emphasize the pressing need for heightened attention directed towards strategies that focus on immune interventions within the female reproductive process, with the ultimate aim of enhancing female fertility.

Single-cell transcriptomic atlas of goat ovarian aging

BACKGROUND

The ovaries are one of the first organs that undergo degenerative changes earlier in the aging process, and ovarian aging is shown by a decrease in the number and quality of oocytes. However, little is known about the molecular mechanisms of female age-related fertility decline in different types of ovarian cells during aging, especially in goats. Therefore, the aim of this study was to reveal the mechanisms driving ovarian aging in goats at single-cell resolution.

RESULTS

For the first time, we surveyed the single-cell transcriptomic landscape of over 27,000 ovarian cells from newborn, young and aging goats, and identified nine ovarian cell types with distinct gene-expression signatures. Functional enrichment analysis showed that ovarian cell types were involved in their own unique biological processes, such as Wnt beta-catenin signalling was enriched in germ cells, whereas ovarian steroidogenesis was enriched in granulosa cells (GCs). Further analysis showed that ovarian aging was linked to GCs-specific changes in the antioxidant system, oxidative phosphorylation, and apoptosis. Subsequently, we identified a series of dynamic genes, such as AMH, CRABP2, THBS1 and TIMP1, which determined the fate of GCs. Additionally, FOXO1, SOX4, and HIF1A were identified as significant regulons that instructed the differentiation of GCs in a distinct manner during ovarian aging.

CONCLUSIONS

This study revealed a comprehensive aging-associated transcriptomic atlas characterizing the cell type-specific mechanisms during ovarian aging at the single-cell level and offers new diagnostic biomarkers and potential therapeutic targets for age-related goat ovarian diseases.

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

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

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

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

A single-cell transcriptomic atlas characterizes cell types and their molecular features in yak ovarian cortex

The ovary as one of the most dynamic organs produces steroids to orchestrate female secondary sexual characteristics, harbors ovarian reserve for oocytes, releases mature oocytes for fertilization, and maintains pregnancy. Yak (Bos grunniens) is the only bovid animal that can adapt to the harsh climatic conditions on the Qinghai-Tibetan Plateau (altitudes of over 3000 m above sea level). However, the cellular atlas is composed of oocytes and other somatic cells, and their individual molecular characteristics remain to be elucidated in the yak ovary. Here, single-cell RNA sequencing (scRNA-seq) was performed to delineate the molecular signature of various cell types in the yak ovarian cortex. A cellular atlas of yak ovarian cortex was constructed successfully on the basis of the differentially expressed genes (DEGs) from the distinct cell types and their functional enrichment analysis, comprising endothelial cells, nature kill cells, stromal cells, smooth muscle cells, oocytes, macrophages, epithelial cells, and granulosa cells. Meanwhile, the signature genes were determined based on their expression specificity in each cell type. A cell-to-cell communication network was built in light of the differentially overexpressed ligand and receptor genes from each cell type. Further, the oocytes were subdivided into four subtypes based on their individual DEGs and the functional enrichment of the DEGs. FST and TOP2A were identified as maker genes for oocytes by immunostaining in the yak ovarian cortex. The cellular atlas reveals the biological characteristics of the ovarian cortex at the cellular molecular level and provides insights into female reproductive biology via cellular communications in the yak.

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

BACKGROUND

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

SUMMARY

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

KEY MESSAGES

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

Single-cell roadmap of human gonadal development

Gonadal development is a complex process that involves sex determination followed by divergent maturation into either testes or ovaries¹. Historically, limited tissue accessibility, a lack of reliable in vitro models and critical differences between humans and mice have hampered our knowledge of human gonadogenesis, despite its importance in gonadal conditions and infertility. Here, we generated a comprehensive map of first- and second-trimester human gonads using a combination of single-cell and spatial transcriptomics, chromatin accessibility assays and fluorescent microscopy. We extracted human-specific regulatory programmes that control the development of germline and somatic cell lineages by profiling equivalent developmental stages in mice. In both species, we define the somatic cell states present at the time of sex specification, including the bipotent early supporting population that, in males, upregulates the testis-determining factor SRY and sPAX8s, a gonadal lineage located at the gonadal-mesonephric interface. In females, we resolve the cellular and molecular events that give rise to the first and second waves of granulosa cells that compartmentalize the developing ovary to modulate germ cell differentiation. In males, we identify human SIGLEC15⁺ and TREM2⁺ fetal testicular macrophages, which signal to somatic cells outside and inside the developing testis cords, respectively. This study provides a comprehensive spatiotemporal map of human and mouse gonadal differentiation, which can guide in vitro gonadogenesis.

Single-Cell Atlas of the Chinese Tongue Sole (Cynoglossus semilaevis) Ovary Reveals Transcriptional Programs of Oogenesis in Fish

Oogenesis is a highly orchestrated process that depends on regulation by autocrine/paracrine hormones and growth factors. However, many details of the molecular mechanisms that regulate fish oogenesis remain elusive. Here, we performed a single-cell RNA sequencing (scRNA-seq) analysis of the molecular signatures of distinct ovarian cell categories in adult Chinese tongue sole (Cynoglossus semilaevis). We characterized the successive stepwise development of three germ cell subtypes. Notably, we identified the cellular composition of fish follicle walls, including four granulosa cell types and one theca cell type, and we proposed important transcription factors (TFs) showing high activity in the regulation of cell identity. Moreover, we found that the extensive niche–germline bidirectional communications regulate fish oogenesis, whereas ovulation in fish is accompanied by the coordination of simultaneous and tightly sequential processes across different granulosa cells. Additionally, a systems biology analysis of the homologous genes shared by Chinese tongue sole and macaques revealed remarkably conserved biological processes in germ cells and granulosa cells across vertebrates. Our results provide key insights into the cell-type-specific mechanisms underlying fish oogenesis at a single-cell resolution, which offers important clues for exploring fish breeding mechanisms and the evolution of vertebrate reproductive systems.

Single cell RNA sequencing techniques and applications in research of ovary development and related diseases

The ovary is a highly organized composite of germ cells and various types of somatic cells, whose communications dictate ovary development to generate functional oocytes. The differences between individual cells might have profound effects on ovary functions. Single cell RNA sequencing techniques are promising approaches to explore the cell type composition of organisms, the dynamics of transcriptomes, the regulatory network between genes and the signaling pathways between cell types at the single cell resolution. In this review, we provide an overview of the currently available single cell RNA sequencing techniques including Smart-seq2 and Drop-seq, as well as their applications in biological and clinical research to provide a better understanding on the molecular mechanisms underlying ovary development and associated diseases.

Exposures during the prepuberty period and future offspring's health: evidence from human cohort studies†

Emerging evidence suggests that exposures in prepuberty, particularly in fathers-to-be, may impact the phenotype of future offspring. Analyses of the RHINESSA cohort find that offspring of father’s exposed to tobacco smoking or overweight that started in prepuberty demonstrate poorer respiratory health in terms of more asthma and lower lung function. A role of prepuberty onset smoking for offspring fat mass is suggested in the RHINESSA and ALSPAC cohorts, and historic studies suggest that ancestral nutrition during prepuberty plays a role for grand-offspring’s health and morbidity. Support for causal relationships between ancestral exposures and (grand-)offspring’s health in humans has been enhanced by advancements in statistical analyses that optimize the gain while accounting for the many complexities and deficiencies in human multigeneration data. The biological mechanisms underlying such observations have been explored in experimental models. A role of sperm small RNA in the transmission of paternal exposures to offspring phenotypes has been established, and chemical exposures and overweight have been shown to influence epigenetic programming in germ cells. For example, exposure of adolescent male mice to smoking led to differences in offspring weight and alterations in small RNAs in the spermatozoa of the exposed fathers. It is plausible that male prepuberty may be a time window of particular susceptibility, given the extensive epigenetic reprogramming taking place in the spermatocyte precursors at this age. In conclusion, epidemiological studies in humans, mechanistic research, and biological plausibility, all support the notion that exposures in the prepuberty of males may influence the phenotype of future offspring.

Immune and vascular contributions to organogenesis of the testis and ovary

Gonad development is a highly regulated process that coordinates cell specification and morphogenesis to produce sex-specific organ structures that are required for fertility, such as testicular seminiferous tubules and ovarian follicles. While sex determination occurs within specialized gonadal supporting cells, sexual differentiation is evident throughout the entire organ, including within the interstitial compartment, which contains immune cells and vasculature. While immune and vascular cells have been traditionally appreciated for their supporting roles during tissue growth and homeostasis, an increasing body of evidence supports the idea that these cell types are critical drivers of sexually dimorphic morphogenesis of the gonad. Myeloid immune cells, such as macrophages, are essential for multiple aspects of gonadogenesis and fertility, including for forming and maintaining gonadal vasculature in both sexes at varying stages of life. While vasculature is long known for supporting organ growth and serving as an export mechanism for gonadal sex steroids in utero, it is also an important component of fetal testicular morphogenesis and differentiation; additionally, it is vital for ovarian corpus luteal function and maintenance of pregnancy. These findings point toward a new paradigm in which immune cells and blood vessels are integral components of sexual differentiation and organogenesis. In this review, we discuss the state of the field regarding the diverse roles of immune and vascular cells during organogenesis of the testis and ovary and highlight outstanding questions in the field that could stimulate new research into these previously underappreciated constituents of the gonad.

Single-cell RNA sequencing reveals species-specific time spans of cell cycle transitions in early oogenesis

Oogenesis is a highly regulated process and its basic cellular events are evolutionarily conserved. However, the time spans of oogenesis differ substantially among species. To explore these interspecies differences in oogenesis, we performed single-cell RNA-sequencing on mouse and monkey female germ cells and downloaded the single-cell RNA-sequencing data of human female germ cells. The cell cycle analyses indicate that the period and extent of cell cycle transitions are significantly different between the species. Moreover, hierarchical clustering of critical cell cycle genes and the interacting network of cell cycle regulators also exhibit distinguished patterns across species. We propose that differences in the regulation of cell cycle transitions may underlie female germ cell developmental allochrony between species. A better understanding of the cell cycle transition machinery will provide new insights into the interspecies differences in female germ cell developmental time spans.

The Application of Single-Cell RNA Sequencing in Mammalian Meiosis Studies

Meiosis is a cellular division process that produces gametes for sexual reproduction. Disruption of complex events throughout meiosis, such as synapsis and homologous recombination, can lead to infertility and aneuploidy. To reveal the molecular mechanisms of these events, transcriptome studies of specific substages must be conducted. However, conventional methods, such as bulk RNA-seq and RT-qPCR, are not able to detect the transcriptional variations effectively and precisely, especially for identifying cell types and stages with subtle differences. In recent years, mammalian meiotic transcriptomes have been intensively studied at the single-cell level by using single-cell RNA-seq (scRNA-seq) approaches, especially through two widely used platforms, Smart-seq2 and Drop-seq. The scRNA-seq protocols along with their downstream analysis enable researchers to accurately identify cell heterogeneities and investigate meiotic transcriptomes at a higher resolution. In this review, we compared bulk RNA-seq and scRNA-seq to show the advantages of the scRNA-seq in meiosis studies; meanwhile, we also pointed out the challenges and limitations of the scRNA-seq. We listed recent findings from mammalian meiosis (male and female) studies where scRNA-seq applied. Next, we summarized the scRNA-seq analysis methods and the meiotic marker genes from spermatocytes and oocytes. Specifically, we emphasized the different features of the two scRNA-seq protocols (Smart-seq2 and Drop-seq) in the context of meiosis studies and discussed their strengths and weaknesses in terms of different research purposes. Finally, we discussed the future applications of scRNA-seq in the meiosis field.