Single-cell multi-omics sequencing of human spermatogenesis reveals a DNA demethylation event associated with male meiotic recombination

SETD1B-mediated broad H3K4me3 controls proper temporal patterns of gene expression critical for spermatid development

Epigenetic programming governs cell fate determination during development through intricately controlling sequential gene activation and repression. Although H3K4me3 is widely recognized as a hallmark of gene activation, its role in modulating transcription output and timing within a continuously developing system remains poorly understood. In this study, we provide a detailed characterization of the epigenomic landscapes in developing male germ cells. We identified thousands of spermatid-specific broad H3K4me3 domains regulated by the SETD1B-RFX2 axis, representing a previously underappreciated form of H3K4me3. These domains, overlapping with H3K27ac-marked enhancers and promoters, play critical roles in orchestrating robust transcription and accurate temporal control of gene expression. Mechanistically, these broad H3K4me3 compete effectively with regular H3K4me3 for transcriptional machinery, thereby ensuring robust levels and precise timing of master gene expression in mouse spermiogenesis. Disruption of this mechanism compromises the accuracy of transcription dosage and timing, ultimately impairing spermiogenesis. Additionally, we unveil remarkable changes in the distribution of heterochromatin marks, including H3K27me3 and H3K9me2, during the mitosis-to-meiosis transition and completion of meiotic recombination, which closely correlates with gene silencing. This work underscores the highly orchestrated epigenetic regulation in spermatogenesis, highlighting the previously unrecognized role of Setd1b in the formation of broad H3K4me3 domains and transcriptional control, and provides an invaluable resource for future studies toward the elucidation of spermatogenesis.

Single-cell sequencing in non-obstructive azoospermia: insights from primary and re-analysis studies

Non-obstructive azoospermia (NOA) constitutes one of the most severe forms of male infertility. Recent advancements in single-cell sequencing have significantly contributed to understanding the molecular landscape of NOA in human testicular tissues, elucidating the factors that underpin spermatogenic dysfunction. This technology has improved our understanding of the condition at a cellular level. Concurrently, bioinformatics developments have facilitated the re-analysis of publicly available single-cell datasets, offering novel insights into the disorder. Nevertheless, a comprehensive review integrating primary and re-analysis studies of single-cell sequencing in NOA is lacking. This review systematically evaluates 10 primary studies reporting original single-cell sequencing data of human NOA testicular samples and 22 secondary studies that re-analyzed these published data. We explore single-cell sequencing applications in germ cells, Sertoli cells, and Leydig cells, offering a comprehensive overview of molecular insights into spermatogenic dysfunction. Our review highlights novel findings in secondary studies, including the roles of transcriptional regulators, RNA transcription, endocrine disruptors, and microtubular cytoskeleton, thereby bridging primary studies and re-analysis studies. Additionally, we discussed future research directions and the challenges of translating single-cell research findings into clinical applications. In summary, single-cell sequencing offers a high-resolution, single-cell perspective of NOA testicular tissue, paving the way for innovative therapeutic strategies in male infertility.

Genome-wide transcriptional silencing and mRNA stabilization allow the coordinated expression of the meiotic program in mice

The transcriptional dynamic of mammalian cells when these transit from the ubiquitous mitotic to a meiotic-specific program is key to understand this switch central to sexual reproduction. By quantifying active RNA polymerase II and nascent transcripts using single cell dataset and ethynyl-uridine pool-down with sorted cells from synchronized testes, we detailed the transcriptional activity of murine male germ cells. When spermatogonia differentiate, transcription slows down, reaching minimal activity at meiotic entry and resumes during pachytene stage. This event, we termed EMLT (for early meiotic low transcription), is distinct from the silencing of sex chromosomes as it is independent of Setdb1, though it is accompanied by the same chromatin mark, H3K9me3. EMLT is delayed in Stra8KO but occurs in mutants altering meiotic chromosome structure or double-strand break formation or repair. By comparing transcript abundance and nascent transcription we unveil a massive event of messenger RNA stabilization that parallels EMLT. Altogether our data indicate that meiosis is initiated with a nearly silent genome, and we propose that the stabilization of transcripts at that time facilitates the meiotic entry by synchronizing the expression of several meiotic subprograms.

Characterization of PRDM9 Multifunctionality in Yak Testes Through Protein Interaction Mapping

Meiotic recombination is initiated by the formation of programmed DNA double-strand breaks during spermatogenesis. PRDM9 determines the localization of recombination hotspots by interacting with several protein complexes in mammals. The function of PRDM9 is not well understood during spermatogenesis in mice or yaks. In this study, we applied yeast two-hybrid assays combined with next-generation sequencing techniques to screen the complete set of PRDM9-interacting proteins and explore its novel functions in yak spermatogenesis. Our results showed that 267 PRDM9-interacting proteins were identified. The gene ontology (GO) analysis of the interacting proteins revealed that the GO terms were primarily associated with spermatogenesis, positive regulation of double-strand break repair via homologous recombination, RNA splicing, the ubiquitin-dependent ERAD pathway, and other biological processes. MKX and PDCD5 were verified to be strongly interacting with PRDM9 and expressed in prophase I of meiosis in both mouse and yak testes. The localizations of RNA splicing genes including THOC5, DDX5, and XRCC6 were expressed in spermatocytes. Cattleyak is the hybrid offspring of a yak and a domestic cow, and the male offspring are sterile. The gene expression of the interacting proteins was also examined in the sterile male hybrid of yak and cattle. Among the 58 detected genes, 55 were downregulated in cattleyak. In conclusion, we established a complete PRDM9 interaction network, and a novel function of PRDM9 was identified, which will further promote our understanding of spermatogenesis. It also provides new insights for the study of hybrid male sterility.

Site-specific DNA demethylation during spermatogenesis presets the sites of nucleosome retention in mouse sperm

DNA methylation patterns are inherited from the parental germline to the embryo. In mature sperm, the sites of unmethylated DNA are tightly coupled to sites of histone retention at gene regulatory elements that are implicated in paternal epigenetic inheritance. The timing and mechanism of site-specific DNA demethylation in the male germline currently remains unknown. Here, we perform genome-wide profiling of DNA methylation during spermatogenesis by capturing methylated DNA through interaction with a methyl-DNA binding protein domain (MBD). Our data demonstrate that there is a site-specific change in DNA methylation during the mitosis-to-meiosis transition. Importantly, the genomic sites that are demethylated during this transition predetermine nucleosome retention sites in spermatozoa. These results suggest that site-specific DNA demethylation during the mitosis-to-meiosis transition of spermatogenesis prepares embryonic gene expression after fertilization. We therefore propose DNA demethylation during spermatogenesis as a novel phase of epigenetic reprogramming that contributes to embryonic gene regulation.

Epigenetic characterization of adult rhesus monkey spermatogonial stem cells identifies key regulators of stem cell homeostasis

Spermatogonial stem cells (SSCs) play crucial roles in the preservation of male fertility. However, successful ex vivo expansion of authentic human SSCs remains elusive due to the inadequate understanding of SSC homeostasis regulation. Using rhesus monkeys (Macaca mulatta) as a representative model, we characterized SSCs and progenitor subsets through single-cell RNA sequencing using a cell-specific network approach. We also profiled chromatin status and major histone modifications (H3K4me1, H3K4me3, H3K27ac, H3K27me3 and H3K9me3), and subsequently mapped promoters and active enhancers in TSPAN33+ putative SSCs. Comparing the epigenetic changes between fresh TSPAN33+ cells and cultured TSPAN33+ cells (resembling progenitors), we identified the regulatory elements with higher activity in SSCs, and the potential transcription factors and signaling pathways implicated in SSC regulation. Specifically, TGF-β signaling is activated in monkey putative SSCs. We provided evidence supporting its role in promoting self-renewal of monkey SSCs in culture. Overall, this study outlines the epigenetic landscapes of monkey SSCs and provides clues for optimization of the culture condition for primate SSCs expansion.

Chromatin accessibility: biological functions, molecular mechanisms and therapeutic application

The dynamic regulation of chromatin accessibility is one of the prominent characteristics of eukaryotic genome. The inaccessible regions are mainly located in heterochromatin, which is multilevel compressed and access restricted. The remaining accessible loci are generally located in the euchromatin, which have less nucleosome occupancy and higher regulatory activity. The opening of chromatin is the most important prerequisite for DNA transcription, replication, and damage repair, which is regulated by genetic, epigenetic, environmental, and other factors, playing a vital role in multiple biological progresses. Currently, based on the susceptibility difference of occupied or free DNA to enzymatic cleavage, solubility, methylation, and transposition, there are many methods to detect chromatin accessibility both in bulk and single-cell level. Through combining with high-throughput sequencing, the genome-wide chromatin accessibility landscape of many tissues and cells types also have been constructed. The chromatin accessibility feature is distinct in different tissues and biological states. Research on the regulation network of chromatin accessibility is crucial for uncovering the secret of various biological processes. In this review, we comprehensively introduced the major functions and mechanisms of chromatin accessibility variation in different physiological and pathological processes, meanwhile, the targeted therapies based on chromatin dynamics regulation are also summarized.

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.

Single-cell sequencing to multi-omics: technologies and applications

Cells, as the fundamental units of life, contain multidimensional spatiotemporal information. Single-cell RNA sequencing (scRNA-seq) is revolutionizing biomedical science by analyzing cellular state and intercellular heterogeneity. Undoubtedly, single-cell transcriptomics has emerged as one of the most vibrant research fields today. With the optimization and innovation of single-cell sequencing technologies, the intricate multidimensional details concealed within cells are gradually unveiled. The combination of scRNA-seq and other multi-omics is at the forefront of the single-cell field. This involves simultaneously measuring various omics data within individual cells, expanding our understanding across a broader spectrum of dimensions. Single-cell multi-omics precisely captures the multidimensional aspects of single-cell transcriptomes, immune repertoire, spatial information, temporal information, epitopes, and other omics in diverse spatiotemporal contexts. In addition to depicting the cell atlas of normal or diseased tissues, it also provides a cornerstone for studying cell differentiation and development patterns, disease heterogeneity, drug resistance mechanisms, and treatment strategies. Herein, we review traditional single-cell sequencing technologies and outline the latest advancements in single-cell multi-omics. We summarize the current status and challenges of applying single-cell multi-omics technologies to biological research and clinical applications. Finally, we discuss the limitations and challenges of single-cell multi-omics and potential strategies to address them.

The role of epigenetics in women's reproductive health: the impact of environmental factors

This paper explores the significant role of epigenetics in women's reproductive health, focusing on the impact of environmental factors. It highlights the crucial link between epigenetic modifications-such as DNA methylation and histones post-translational modifications-and reproductive health issues, including infertility and pregnancy complications. The paper reviews the influence of pollutants like PM2.5, heavy metals, and endocrine disruptors on gene expression through epigenetic mechanisms, emphasizing the need for understanding how dietary, lifestyle choices, and exposure to chemicals affect gene expression and reproductive health. Future research directions include deeper investigation into epigenetics in female reproductive health and leveraging gene editing to mitigate epigenetic changes for improving IVF success rates and managing reproductive disorders.

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.

The mechanism of INO80D involved in chromatin remodeling regulating spermatogenesis in Chinese mitten crab (Eriocheir sinensis)

The INO80D protein, a component of the INO80 chromatin remodeling complex, plays a pivotal role in chromatin remodeling, gene expression, and DNA repair within mammalian sperm. In contrast to the condensed nuclear structure of mammalian sperm, Chinese mitten crab, Eriocheir sinensis, exhibits a distinctively decondensed sperm nucleus. The distribution and function of INO80D during the E. sinensis spermatogenesis were previously enigmatic. Our research endeavored to elucidate the distribution and function of INO80D, thereby enhancing our comprehension of sperm decondensation and the process of spermatogenesis in this species. Employing transcriptome sequencing, RT-qPCR, western blot analysis, and immunofluorescence techniques, we observed a pronounced upregulation of INO80D in the adult E. sinensis in comparison to the juvenile. The protein predominantly resides in the cellular nucleus, with high levels in spermatogonia and spermatocytes, less in stage I and III spermatids, and lowest in mature sperm. The results indicated that INO80D is initially instrumental in chromatin decondensation to facilitate gene accessibility and DNA repair during the early phases of spermatogenesis. Its role subsequently shifts to maintaining decondensed chromatin stability and genetic integrity during spermiogenesis. The sustained presence of INO80D during spermiogenesis is essential for the ultimate maturation of the decondensed sperm nucleus, imperative for preserving the unique decondensed state and the protection of genetic material in E. sinensis. Our study concludes that INO80D exerts a multifaceted influence on the spermatogenesis of E. sinensis, impacting chromatin decondensation, genetic integrity, and the regulation of early gene expression. This understanding could potentially improve crab breeding in aquaculture.

Understanding testicular single cell transcriptional atlas: from developmental complications to male infertility

Spermatogenesis is a multi-step biological process where mitotically active diploid (2n) spermatogonia differentiate into haploid (n) spermatozoa via regulated meiotic programming. The alarming rise in male infertility has become a global concern during the past decade thereby demanding an extensive profiling of testicular gene expression. Advancements in Next-Generation Sequencing (NGS) technologies have revolutionized our empathy towards complex biological events including spermatogenesis. However, despite multiple attempts made in the past to reveal the testicular transcriptional signature(s) either with bulk tissues or at the single-cell, level, comprehensive reviews on testicular transcriptomics and associated disorders are limited. Notably, technologies explicating the genome-wide gene expression patterns during various stages of spermatogenic progression provide the dynamic molecular landscape of testicular transcription. Our review discusses the advantages of single-cell RNA-sequencing (Sc-RNA-seq) over bulk RNA-seq concerning testicular tissues. Additionally, we highlight the cellular heterogeneity, spatial transcriptomics, dynamic gene expression and cell-to-cell interactions with distinct cell populations within the testes including germ cells (Gc), Sertoli cells (Sc), Peritubular cells (PTc), Leydig cells (Lc), etc. Furthermore, we provide a summary of key finding of single-cell transcriptomic studies that have shed light on developmental mechanisms implicated in testicular disorders and male infertility. These insights emphasize the pivotal roles of Sc-RNA-seq in advancing our knowledge regarding testicular transcriptional landscape and may serve as a potential resource to formulate future clinical interventions for male reproductive health.

Genome-wide DNA methylation changes in human spermatogenesis

Sperm production and function require the correct establishment of DNA methylation patterns in the germline. Here, we examined the genome-wide DNA methylation changes during human spermatogenesis and its alterations in disturbed spermatogenesis. We found that spermatogenesis is associated with remodeling of the methylome, comprising a global decline in DNA methylation in primary spermatocytes followed by selective remethylation, resulting in a spermatids/sperm-specific methylome. Hypomethylated regions in spermatids/sperm were enriched in specific transcription factor binding sites for DMRT and SOX family members and spermatid-specific genes. Intriguingly, while SINEs displayed differential methylation throughout spermatogenesis, LINEs appeared to be protected from changes in DNA methylation. In disturbed spermatogenesis, germ cells exhibited considerable DNA methylation changes, which were significantly enriched at transposable elements and genes involved in spermatogenesis. We detected hypomethylation in SVA and L1HS in disturbed spermatogenesis, suggesting an association between the abnormal programming of these regions and failure of germ cells progressing beyond meiosis.

The Expression and Epigenetic Characteristics of the HSF2 Gene in Cattle-Yak and the Correlation with Its Male Sterility

Aberrant expression of the heat shock proteins and factors was revealed to be closely associated with male reproduction. Heat shock factor 2 (HSF2) is a transcription factor that is involved in the regulation of diverse developmental pathways. However, the role and the corresponding molecular mechanism of HSF2 in male cattle-yak sterility are still poorly understood. Therefore, the aim of this study was to obtain the sequence and the biological information of the cattle-yak HSF2 gene and to investigate the spatiotemporal expression profiles of the locus during the development of cattle-yak testes. Additionally, the differential expression was analyzed between the cattle-yak and the yak, and the methylation of corresponding promoter regions was compared. Our results showed an additional 54 bp fragment and a missense mutation (lysine to glutamic acid) were presented in the cattle-yak HSF2 gene, which correlated with enriched expression in testicular tissue. In addition, the expression of the HSF2 gene showed dynamic changes during the growth of the testes, reaching a peak in adulthood. The IHC indicated that HSF2 protein was primarily located in spermatocytes (PS), spermatogonia (SP), and Sertoli cells (SC) in cattle-yak testes, compared with the corresponding cells of cattle and the yak. Furthermore, bisulfite-sequencing PCR (BSP) revealed that the methylated CpG sites in the promoter region of the cattle-yak HSF2 were more numerous than in the yak counterpart, which suggests hypermethylation of this region in the cattle-yak. Taken together, the low expression abundance and hypermethylation of HSF2 may underpin the obstruction of spermatogenesis, which leads to male cattle-yak infertility. Our study provided a basic guideline for the HSF2 gene in male reproduction and a new insight into the mechanisms of male cattle-yak sterility.

Single-cell multi-omics profiling of human preimplantation embryos identifies cytoskeletal defects during embryonic arrest

Human in vitro fertilized embryos exhibit low developmental capabilities, and the mechanisms that underlie embryonic arrest remain unclear. Here using a single-cell multi-omics sequencing approach, we simultaneously analysed alterations in the transcriptome, chromatin accessibility and the DNA methylome in human embryonic arrest due to unexplained reasons. Arrested embryos displayed transcriptome disorders, including a distorted microtubule cytoskeleton, increased genomic instability and impaired glycolysis, which were coordinated with multiple epigenetic reprogramming defects. We identified Aurora A kinase (AURKA) repression as a cause of embryonic arrest. Mechanistically, arrested embryos induced through AURKA inhibition resembled the reprogramming abnormalities of natural embryonic arrest in terms of the transcriptome, the DNA methylome, chromatin accessibility and H3K4me3 modifications. Mitosis-independent sequential activation of the zygotic genome in arrested embryos showed that YY1 contributed to human major zygotic genome activation. Collectively, our study decodes the reprogramming abnormalities and mechanisms of human embryonic arrest and the key regulators of zygotic genome activation.