Single-cell RNA-seq unravels alterations of the human spermatogonial stem cell compartment in patients with impaired spermatogenesis

Identification of two hidden clinical subgroups among men with idiopathic cryptozoospermia

STUDY QUESTION

Are there subgroups among patients with cryptozoospermia pointing to distinct etiologies?

SUMMARY ANSWER

We reveal two distinct subgroups of cryptozoospermic (Crypto) patients based on testicular tissue composition, testicular volume, and FSH levels.

WHAT IS KNOWN ALREADY

Cryptozoospermic patients present with a sperm concentration below 0.1 million/ml. While the etiology of the severely impaired spermatogenesis remains largely unknown, alterations of the spermatogonial compartment have been reported including a reduction of the reserve stem cells in these patients.

STUDY DESIGN, SIZE, DURATION

To assess whether there are distinct subgroups among cryptozoospermic patients, we applied the statistical method of cluster analysis. For this, we retrospectively selected 132 cryptozoospermic patients from a clinical database who underwent a testicular biopsy in the frame of fertility treatment at a university hospital. As controls (Control), we selected 160 patients with obstructive azoospermia and full spermatogenesis. All 292 patients underwent routine evaluation for endocrine, semen, and histological parameters (i.e. the percentage of tubules with elongated spermatids). Moreover, outcome of medically assisted reproduction (MAR) was assessed for cryptozoospermic (n = 73) and Control patients (n = 87), respectively. For in-depth immunohistochemical and histomorphometrical analyses, representative tissue samples from cryptozoospermic (n = 27) and Control patients (n = 12) were selected based on cluster analysis results and histological parameters.

PARTICIPANTS/MATERIALS, SETTING, METHODS

This study included two parts: firstly using clinical parameters of the entire cohort of 292 patients, we performed principal component analysis (PCA) followed by hierarchical clustering on principal components (i.e. considering hormonal values, ejaculate parameters, and histological information). Secondly, for histological analyses seminiferous tubules were categorized according to the most advanced germ cell type present in sections stained with Periodic acid Schif. On the selected cohort of 39 patients (12 Control, 27 cryptozoospermic), we performed immunohistochemistry for spermatogonial markers melanoma-associated antigen 4 (MAGEA4) and piwi like RNA-mediated gene silencing 4 (PIWIL4) followed by quantitative analyses. Moreover, the morphologically defined Adark spermatogonia, which are considered to be the reserve stem cells, were quantified.

MAIN RESULTS AND THE ROLE OF CHANCE

The PCA and hierarchical clustering revealed three different clusters, one of them containing all Control samples. The main factors driving the sorting of patients to the clusters were the percentage of tubules with elongated spermatids (Cluster 1, all Control patients and two cryptozoospermic patients), the percentage of tubules with spermatocytes (Cluster 2, cryptozoospermic patients), and tubules showing a Sertoli cells only phenotype (Cluster 3, cryptozoospermic patients). Importantly, the percentage of tubules containing elongated spermatids was comparable between Clusters 2 and 3. Additional differences were higher FSH levels (P < 0.001) and lower testicular volumes (P < 0.001) in Cluster 3 compared to Cluster 2. In the spermatogonial compartment of both cryptozoospermic Clusters, we found lower numbers of MAGEA4+ and Adark spermatogonia but higher proportions of PIWIL4+ spermatogonia, which were significantly correlated with a lower percentage of tubules containing elongated spermatids. In line with this common alteration, the outcome of MAR was comparable between Controls as well as both cryptozoospermic Clusters.

LIMITATIONS, REASONS FOR CAUTION

While we have uncovered the existence of subgroups within the cohort of cryptozoospermic patients, comprehensive genetic analyses remain to be performed to unravel potentially distinct etiologies.

WIDER IMPLICATIONS OF THE FINDINGS

The novel insight that cryptozoospermic patients can be divided into two subgroups will facilitate the strategic search for underlying genetic etiologies. Moreover, the shared alterations of the spermatogonial stem cell compartment between the two cryptozoospermic subgroups could represent a general response mechanism to the reduced output of sperm, which may be associated with a progressive phenotype. This study therefore offers novel approaches towards the understanding of the etiology underlying the reduced sperm formation in cryptozoospermic patients.

STUDY FUNDING/COMPETING INTEREST(S)

German research foundation CRU 326 (grants to: SDP, NN). Moreover, we thank the Faculty of Medicine of the University of Münster for the financial support of Lena Charlotte Schülke through the MedK-program. We acknowledge support from the Open Access Publication Fund of the University of Münster. The authors have no potential conflicts of interest.

TRIAL REGISTRATION NUMBER

N/A.

Adult Human, but Not Rodent, Spermatogonial Stem Cells Retain States with a Foetal-like Signature

Spermatogenesis involves a complex process of cellular differentiation maintained by spermatogonial stem cells (SSCs). Being critical to male reproduction, it is generally assumed that spermatogenesis starts and ends in equivalent transcriptional states in related species. Based on single-cell gene expression profiling, it has been proposed that undifferentiated human spermatogonia can be subclassified into four heterogenous subtypes, termed states 0, 0A, 0B, and 1. To increase the resolution of the undifferentiated compartment and trace the origin of the spermatogenic trajectory, we re-analysed the single-cell (sc) RNA-sequencing libraries of 34 post-pubescent human testes to generate an integrated atlas of germ cell differentiation. We then used this atlas to perform comparative analyses of the putative SSC transcriptome both across human development (using 28 foetal and pre-pubertal scRNA-seq libraries) and across species (including data from sheep, pig, buffalo, rhesus and cynomolgus macaque, rat, and mouse). Alongside its detailed characterisation, we show that the transcriptional heterogeneity of the undifferentiated spermatogonial cell compartment varies not only between species but across development. Our findings associate 'state 0B' with a suppressive transcriptomic programme that, in adult humans, acts to functionally oppose proliferation and maintain cells in a ready-to-react state. Consistent with this conclusion, we show that human foetal germ cells-which are mitotically arrested-can be characterised solely as state 0B. While germ cells with a state 0B signature are also present in foetal mice (and are likely conserved at this stage throughout mammals), they are not maintained into adulthood. We conjecture that in rodents, the foetal-like state 0B differentiates at birth into the renewing SSC population, whereas in humans it is maintained as a reserve population, supporting testicular homeostasis over a longer reproductive lifespan while reducing mutagenic load. Together, these results suggest that SSCs adopt differing evolutionary strategies across species to ensure fertility and genome integrity over vastly differing life histories and reproductive timeframes.

DDX20 is required for cell-cycle reentry of prospermatogonia and establishment of spermatogonial stem cell pool during testicular development in mice

RNA-binding proteins (RBPs), as key regulators of mRNA fate, are abundantly expressed in the testis. However, RBPs associated with human male infertility remain largely unknown. Through bioinformatic analyses, we identified 62 such RBPs, including an evolutionarily conserved RBP, DEAD-box helicase 20 (DDX20). Male germ-cell-specific inactivation of Ddx20 at E15.5 caused T1-propsermatogonia (T1-ProSG) to fail to reenter cell cycle during the first week of testicular development in mice. Consequently, neither the foundational spermatogonial stem cell (SSC) pool nor progenitor spermatogonia were ever formed in the knockout testes. Mechanistically, DDX20 functions to control the translation of its target mRNAs, many of which encode cell-cycle-related regulators, by interacting with key components of the translational machinery in prospermatogonia. Our data demonstrate a previously unreported function of DDX20 as a translational regulator of critical cell-cycle-related genes, which is essential for cell-cycle reentry of T1-ProSG and formation of the SSC pool.

Identification of EGR4 as a prospective target for inhibiting tumor cell proliferation and a novel biomarker in colorectal cancer

EGR4 (Early Growth Response 4) is a member of the EGR family, involving in tumorigenesis. However, the function and action mechanism of EGR4 in the pathogenesis of colorectal cancer (CRC) remain unclear. To address this, we assessed the prognosis of CRC based on EGR4 using the Kaplan-Meier plotter tool and tissue microarray. The abundance of immunoinfiltration was evaluated through ssGSEA, TISIDB, and TIMER. In vitro experiments involving knockdown or overexpression of EGR4 were performed, and RNA-sequencing was conducted to explore potential mechanisms. Furthermore, we used oxaliplatin and 5-fluorouracil to validate the impact of EGR4 on chemo-resistance. Pan-cancer analysis and tissue microarray showed that EGR4 was highly expressed in CRC and significantly correlated with an unfavorable prognosis. Moreover, EGR4 expression was associated with immunoinfiltration and cancer-associated fibroblasts in the CRC microenvironment. Functional enrichment demonstrated that high-expressional EGR4 were involved in chromatin and nucleosome assembly. Additionally, EGR4 promoted the proliferation of CRC cells. Mechanistically, EGR4 upregulated TNFα to activate the NF-κB signaling pathway, and its knockdown reduced p65 nuclear translocation. Importantly, combining shEGR4 with oxaliplatin and 5-fluorouracil significantly inhibited CRC proliferation. Taken together, these findings provide new insights into the potential prognosis and therapeutic targets of EGR4 in CRC.

Scrutinizing the human TEX genes in the context of human male infertility

BACKGROUND

Infertility affects around 15% of all couples worldwide and is increasingly linked to variants in genes specifically expressed in the testis. Well-established causes of male infertility include pathogenic variants in the genes TEX11, TEX14, and TEX15, while few studies have recently reported variants in TEX13B, TEX13C, FAM9A (TEX39A), and FAM9B (TEX39B).

OBJECTIVES

We aimed at screening for novel potential candidate genes among the human TEX ("testis expressed") genes as well as verifying previously described disease associations in this set of genes.

MATERIALS AND METHODS

To this end, we screened the exome sequencing data of 1305 men, including 1056 crypto- and azoospermic individuals, and determined cell-specific expression by analyzing testis-specific single-cell RNA sequencing data for genes with identified variants. To investigate the overarching role in male fertility, we generated testis-specific knockdown (KD) models of all 10 orthologous TEX genes in Drosophila melanogaster.

RESULTS

We detected rare potential disease-causing variants in TEX10, TEX13A, TEX13B, TEX13C, TEX13D, ZFAND3 (TEX27), TEX33, FAM9A (TEX39A), and FAM9B (TEX39B), in 28 infertile men, of which 15 men carried variants in TEX10, TEX27, and TEX33. The KD of TEX2, TEX9, TEX10, TEX13, ZFAND3 (TEX27), TEX28, TEX30, NFX1 (TEX42), TEX261, and UTP4 (TEX292) in Drosophila resulted in normal fertility.

DISCUSSION

Based on our findings, the autosomal dominant predicted genes TEX10 and ZFAND3 (TEX27) and the autosomal recessive predicted gene TEX33, which all three are conceivably required for germ cell maturation, were identified as novel potential candidate genes for human non-obstructive azoospermia. We additionally identified hemizygous loss-of-function (LoF) variants in TEX13B, TEX13C, and FAM9A (TEX39A) as unlikely monogenic culprits of male infertility as LoF variants were also found in control men.

CONCLUSION

Our findings concerning the X-linked genes TEX13B, TEX13C, and FAM9A (TEX39A) contradict previous reports and will decrease false-positive reports in genetic diagnostics of azoospermic men.

Molecular mechanisms of cellular dysfunction in testes from men with non-obstructive azoospermia

Male factor infertility affects 50% of infertile couples worldwide; the most severe form, non-obstructive azoospermia (NOA), affects 10-15% of infertile males. Treatment for individuals with NOA is limited to microsurgical sperm extraction paired with in vitro fertilization intracytoplasmic sperm injection. Unfortunately, spermatozoa are only retrieved in ~50% of patients, resulting in live birth rates of 21-46%. Regenerative therapies could provide a solution; however, understanding the cell-type-specific mechanisms of cellular dysfunction is a fundamental necessity to develop precision medicine strategies that could overcome these abnormalities and promote regeneration of spermatogenesis. A number of mechanisms of cellular dysfunction have been elucidated in NOA testicular cells. These mechanisms include abnormalities in both somatic cells and germ cells in NOA testes, such as somatic cell immaturity, aberrant growth factor signalling, increased inflammation, increased apoptosis and abnormal extracellular matrix regulation. Future cell-type-specific investigations in identifying modulators of cellular transcription and translation will be key to understanding upstream dysregulation, and these studies will require development of in vitro models to functionally interrogate spermatogenic niche dysfunction in both somatic and germ cells.

The testicular microvasculature in Klinefelter syndrome is immature with compromised integrity and characterized by excessive inflammatory cross-talk

STUDY QUESTION

Does Klinefelter syndrome (KS) lead to a distinct gene expression pattern at single-cell level in the testes that could provide insight into the reported microvascular dysfunction in the testes?

SUMMARY ANSWER

A distinct gene expression pattern within microvascular-associated cells of males with KS suggests excessive endothelial cell (EC) activation, disorganized vessel formation, and the presence of immature vessels with compromised integrity.

WHAT IS KNOWN ALREADY

Recent studies show that males with KS exhibit microvascular dysfunction in their testes, which affects blood flow and is associated with lower circulating levels of testosterone.

STUDY DESIGN, SIZE, DURATION

A comparative cross-sectional study of males with KS (n = 6), non-obstructive azoospermia (NOA) (n = 5), cryptozoospermia (n = 3), and controls (n = 15) was carried out.

PARTICIPANTS/MATERIALS, SETTING, METHODS

We analyzed publicly available single-cell RNA sequencing data of testicular cells from males with KS, males with NOA, males with cryptozoospermia, and controls. The integration of these datasets allowed us to analyze gene expression profiles and communication patterns among the cell types within the testis and to identify capillary ECs to investigate changes at the microvascular level.

MAIN RESULTS AND THE ROLE OF CHANCE

Rooted in changes at the single-cell level, our study demonstrates a shift in gene expression forming the foundation for altered cellular communication, microvascular remodeling, and pro-inflammatory responses within the testes of males with KS. We identified genes that were dysregulated in capillary ECs from males with KS (Padj < 0.05). Specifically, the unique microvascular gene expression in males with KS indicated enhanced capillary EC activation and increased inflammatory cross-talk, leading to impaired vessel maturation and increased EC barrier permeability.

LIMITATIONS, REASONS FOR CAUTION

Our study is constrained by an unbalanced design, with varying sample sizes and number of cells within each group. We acknowledge the restricted access to clinical information. In addition, our findings were deduced from changes in gene expression, which limits us to infer potential biological consequences arising from these alterations. Furthermore, the absence of a pre-pubertal age group limits the generalizability of our findings and warrants further investigation.

WIDER IMPLICATIONS OF THE FINDINGS

This study offers novel insights into the testicular pathophysiology in KS and underscores the potential contribution of microvascular dysfunction to the hypogonadism and infertility observed in males with KS. While this study aims to better understand the microvascular dysfunction in KS, the precise connections to testosterone deficiency and testicular atrophy remain to be fully elucidated.

STUDY FUNDING/COMPETING INTEREST(S)

A.S. was supported by the Independent Research Fund Denmark (0134-00130B). C.H.G. was supported by Novo Nordisk Foundation (NNF15OC0016474, NNF20OC0060610), 'Fonden til lægevidenskabens fremme', the Familien Hede Nielsen foundation and the Independent Research Fund Denmark (0134-00406A). E.B.J. was supported by Aarhus University and E.B.J. and C.H.G by the Independent Research Fund Denmark (2096-00165A). J.M.K. was supported by Lundbeckfonden (R307-2018-3667), Carlsberg Fonden (CF19-0687), Novo Nordisk Fonden (0073440) and Steno Diabetes Center Aarhus (SDCA). The authors declare no conflicts of interest.

TRIAL REGISTRATION NUMBER

N/A.

The construction of a testis transcriptional cell atlas from embryo to adult reveals various somatic cells and their molecular roles

BACKGROUND

The testis is a complex organ that undergoes extensive developmental changes from the embryonic stage to adulthood. The development of germ cells, which give rise to spermatozoa, is tightly regulated by the surrounding somatic cells.

METHODS

To better understand the dynamics of these changes, we constructed a transcriptional cell atlas of the testis, integrating single-cell RNA sequencing data from over 26,000 cells across five developmental stages: fetal germ cells, infants, childhood, peri-puberty, and adults. We employed various analytical techniques, including clustering, cell type assignments, identification of differentially expressed genes, pseudotime analysis, weighted gene co-expression network analysis, and evaluation of paracrine cell-cell communication, to comprehensively analyze this transcriptional cell atlas of the testis.

RESULTS

Our analysis revealed remarkable heterogeneity in both somatic and germ cell populations, with the highest diversity observed in Sertoli and Myoid somatic cells, as well as in spermatogonia, spermatocyte, and spermatid germ cells. We also identified key somatic cell genes, including RPL39, RPL10, RPL13A, FTH1, RPS2, and RPL18A, which were highly influential in the weighted gene co-expression network of the testis transcriptional cell atlas and have been previously implicated in male infertility. Additionally, our analysis of paracrine cell-cell communication supported specific ligand-receptor interactions involved in neuroactive, cAMP, and estrogen signaling pathways, which support the crucial role of somatic cells in regulating germ cell development.

CONCLUSIONS

Overall, our transcriptional atlas provides a comprehensive view of the cell-to-cell heterogeneity in the testis and identifies key somatic cell genes and pathways that play a central role in male fertility across developmental stages.

GEFAAR: a generic framework for the analysis of antimicrobial resistance providing statistics and cluster analyses

Easy access to antimicrobial resistance data and meaningful visualization is essential to guide the empirical antimicrobial treatment and to promote the rational use of antimicrobial agents. Currently available solutions are commonly externally hosted, centralized systems. However, there is a need for close monitoring by local analysis tools. To fill this gap, we developed GEFAAR-a generic framework for the analysis of antimicrobial resistance data. Following the example of the German Robert Koch Institute (RKI), an interactive web-application is provided to determine basic pathogen and resistance statistics. In addition to the RKI's externally maintained database, our application provides a generic framework to import tabular data and to analyze them safely in a local environment. Moreover, our application offers an intuitive web-based user interface to visualize resistance trend analysis as well as advanced cluster analyses on species- or clinic/unit level to generate alerts of potential transmission events.

CSNK1G2-AS1 promotes metastasis, colony formation and serves as a biomarker in testicular germ cell tumor cells

Background/Aim: Some long non-coding RNAs (lncRNAs) have been found to significantly participate in the progression of TGCTs. In comparison to the normal testis, the TGCT tissues showed significantly decreased CSNK1G2-AS1 expression, however, its effect on TGCTs and its mechanism are still unclear. The aim of this study is to investigate the effect of CSNK1G2-AS1 on TGCTs and explore the mechanism underlying its effect on TGCTs. Materials and Methods: In this study, to evaluate the expression of CSNK1G2-AS1 in tissue samples from TGCTs, the UCSC and GEPIA databases were applied and qRT-PCR was conducted. The Kaplan-Meier Plotter was applied to analyze the correlation between CSNK1G2-AS1 methylation levels and the prognosis of TGCTs patients. The assays of MTS, clone formation, transwell, and flow cytometry were performed to investigate the effect of CSNK1G2-AS1 overexpression on the proliferation, metastasis, and apoptosis of TGCT cells, respectively. Finally, western blotting was conducted to determine the expressions of the proteins associated with EMT and AKT. Results: Our study first found that, compared to the normal testis, TGCTs tissue showed significantly decreased CSNK1G2-AS1 expression, and hypomethylation of CSNK1G2-AS1 was significantly correlated with a better prognosis with TGCTs patients. In vitro, we found that overexpression of CSNK1G2-AS1 dramatically promoted the clone formation, invasion, and migration of TGCT cells, but inhibited apoptosis. And CSNK1G2-AS1 overexpression significantly decreased the expression of EMT-associated proteins ZO-1 but increased the expression and phosphorylation of AKT. Conclusions: CSNK1G2-AS1 may play an essential role in the pathogenesis and metastasis of TGCTs through the EMT- and AKT-mediated signal pathways.

OGRE: calculate, visualize, and analyze overlap between genomic input regions and public annotations

BACKGROUND

Modern genome sequencing leads to an ever-growing collection of genomic annotations. Combining these elements with a set of input regions (e.g. genes) would yield new insights in genomic associations, such as those involved in gene regulation. The required data are scattered across different databases making a manual approach tiresome, unpractical, and prone to error. Semi-automatic approaches require programming skills in data parsing, processing, overlap calculation, and visualization, which most biomedical researchers lack. Our aim was to develop an automated tool providing all necessary algorithms, benefiting both bioinformaticians and researchers without bioinformatic training.

RESULTS

We developed overlapping annotated genomic regions (OGRE) as a comprehensive tool to associate and visualize input regions with genomic annotations. It does so by parsing regions of interest, mining publicly available annotations, and calculating possible overlaps between them. The user can thus identify location, type, and number of associated regulatory elements. Results are presented as easy to understand visualizations and result tables. We applied OGRE to recent studies and could show high reproducibility and potential new insights. To demonstrate OGRE's performance in terms of running time and output, we have conducted a benchmark and compared its features with similar tools.

CONCLUSIONS

OGRE's functions and built-in annotations can be applied as a downstream overlap association step, which is compatible with most genomic sequencing outputs, and can thus enrich pre-existing analyses pipelines. Compared to similar tools, OGRE shows competitive performance, offers additional features, and has been successfully applied to two recent studies. Overall, OGRE addresses the lack of tools for automatic analysis, local genomic overlap calculation, and visualization by providing an easy to use, end-to-end solution for both biologists and computational scientists.

Potency Assays for Mesenchymal Stromal Cell Secretome-Based Products for Tissue Regeneration

Adult stem cells maintaining tissue homeostasis and regeneration are tightly regulated by their specific microenvironments or stem cell niches. The dysfunction of niche components may alter the activity of stem cells and ultimately lead to intractable chronic or acute disorders. To overcome this dysfunction, niche-targeting regenerative medicine treatments such as gene, cell, and tissue therapy are actively investigated. Here, multipotent mesenchymal stromal cells (MSCs), and particularly their secretomes, are of high interest due to their potency to recover and reactivate damaged or lost stem cell niches. However, a workflow for the development of MSC secretome-based products is not fully covered by regulatory authorities, and and this issue significantly complicates their clinical translation and has possibly been expressed in a huge number of failed clinical trials. One of the most critical issues in this regard relates to the development of potency assays. In this review, guidelines for biologicals and cell therapies are considered to be applied for the development of potency assays for the MSC secretome-based products that aim for tissue regeneration. Specific attention is paid to their possible effects on stem cell niches and to a spermatogonial stem cell niche in particular.

Machine learning and integrative analysis identify the common pathogenesis of azoospermia complicated with COVID-19

Background

Although more recent evidence has indicated COVID-19 is prone to azoospermia, the common molecular mechanism of its occurrence remains to be elucidated. The aim of the present study is to further investigate the mechanism of this complication.

Methods

To discover the common differentially expressed genes (DEGs) and pathways of azoospermia and COVID-19, integrated weighted co-expression network (WGCNA), multiple machine learning analyses, and single-cell RNA-sequencing (scRNA-seq) were performed.

Results

Therefore, we screened two key network modules in the obstructive azoospermia (OA) and non-obstructive azoospermia (NOA) samples. The differentially expressed genes were mainly related to the immune system and infectious virus diseases. We then used multiple machine learning methods to detect biomarkers that differentiated OA from NOA. Enrichment analysis showed that azoospermia patients and COVID-19 patients shared a common IL-17 signaling pathway. In addition, GLO1, GPR135, DYNLL2, and EPB41L3 were identified as significant hub genes in these two diseases. Screening of two different molecular subtypes revealed that azoospermia-related genes were associated with clinicopathological characteristics of age, hospital-free-days, ventilator-free-days, charlson score, and d-dimer of patients with COVID-19 (P < 0.05). Finally, we used the Xsum method to predict potential drugs and single-cell sequencing data to further characterize whether azoospermia-related genes could validate the biological patterns of impaired spermatogenesis in cryptozoospermia patients.

Conclusion

Our study performs a comprehensive and integrated bioinformatics analysis of azoospermia and COVID-19. These hub genes and common pathways may provide new insights for further mechanism research.

WWC2 expression in the testis: Implications for spermatogenesis and male fertility

The family of WWC proteins is known to regulate cell proliferation and organ growth control via the Hippo signaling pathway. As WWC proteins share a similar domain structure and a common set of interacting proteins, they are supposed to fulfill compensatory functions in cells and tissues. While all three WWC family members WWC1, WWC2, and WWC3 are found co-expressed in most human organs including lung, brain, kidney, and liver, in the testis only WWC2 displays a relatively high expression. In this study, we investigated the testicular WWC2 expression in spermatogenesis and male fertility. We show that the Wwc2 mRNA expression level in mouse testes is increased during development in parallel with germ cell proliferation and differentiation. The cellular expression of each individual WWC family member was evaluated in published single-cell mRNA datasets of murine and human testes demonstrating a high WWC2 expression predominantly in early spermatocytes. In line with this, immunohistochemistry revealed cytosolic WWC2 protein expression in primary spermatocytes from human testes displaying full spermatogenesis. In accordance with these findings, markedly lower WWC2 expression levels were detected in testicular tissues from mice and men lacking germ cells. Finally, analysis of whole-exome sequencing data of male patients affected by infertility and unexplained severe spermatogenic failure revealed several heterozygous, rare WWC2 gene variants with a proposed damaging function and putative impact on WWC2 protein structure. Taken together, our findings provide novel insights into the testicular expression of WWC2 and show its cell-specific expression in spermatocytes. As rare WWC2 variants were identified in the background of disturbed spermatogenesis, WWC2 may be a novel candidate gene for male infertility.

Deciphering the Molecular Characteristics of Human Idiopathic Nonobstructive Azoospermia from the Perspective of Germ Cells

Nonobstructive azoospermia (NOA) is one of the most important causes of male infertility, accounting for 10-15% of infertile men worldwide. Among these, more than 70% of cases are idiopathic NOA (iNOA), whose pathogenesis and molecular basis remain unknown. This work profiles 3696 human testicular single-cell transcriptomes from 17 iNOA patients, which are classified into four classes with different arrest periods and variable cell proportions based on the gene expression patterns and pathological features. Genes related to the cell cycle, energy production, and gamete generation show obvious abnormalities in iNOA germ cells. This work identifies several candidate causal genes for iNOA, including CD164, LELP1, and TEX38, which are significantly downregulated in iNOA germ cells. Notably, CD164 knockdown promotes apoptosis in spermatogonia. Cellular communications between spermatogonial stem cells and Sertoli cells are disturbed in iNOA patients. Moreover, BOD1L2, C1orf194, and KRTCAP2 are found to indicate testicular spermatogenic capacity in a variety of testicular diseases, such as Y-chromosome microdeletions and Klinefelter syndrome. In general, this study analyzes the pathogenesis of iNOA from the perspective of germ cell development, transcription factor (TF) regulatory networks, as well as germ cell and somatic cell interactions, which provides new ideas for clinical diagnosis.

Conserved Transcriptome Features Define Prepubertal Primate Spermatogonial Stem Cells as Adark Spermatogonia and Identify Unique Regulators

Antineoplastic treatments for cancer and other non-malignant disorders can result in long-term or permanent male infertility by ablating spermatogonial stem cells (SSCs). SSC transplantation using testicular tissue harvested before a sterilizing treatment is a promising approach for restoring male fertility in these cases, but a lack of exclusive biomarkers to unequivocally identify prepubertal SSCs limits their therapeutic potential. To address this, we performed single-cell RNA-seq on testis cells from immature baboons and macaques and compared these cells with published data from prepubertal human testis cells and functionally-defined mouse SSCs. While we found discrete groups of human spermatogonia, baboon and rhesus spermatogonia appeared less heterogenous. A cross-species analysis revealed cell types analogous to human SSCs in baboon and rhesus germ cells, but a comparison with mouse SSCs revealed significant differences with primate SSCs. Primate-specific SSC genes were enriched for components and regulators of the actin cytoskeleton and participate in cell-adhesion, which may explain why the culture conditions for rodent SSCs are not appropriate for primate SSCs. Furthermore, correlating the molecular definitions of human SSC, progenitor and differentiating spermatogonia with the histological definitions of Adark/Apale spermatogonia indicates that both SSCs and progenitor spermatogonia are Adark, while Apale spermatogonia appear biased towards differentiation. These results resolve the molecular identity of prepubertal human SSCs, define novel pathways that could be leveraged for advancing their selection and propagation in vitro, and confirm that the human SSC pool resides entirely within Adark spermatogonia.

Transcriptome analyses in infertile men reveal germ cell-specific expression and splicing patterns

The process of spermatogenesis-when germ cells differentiate into sperm-is tightly regulated, and misregulation in gene expression is likely to be involved in the physiopathology of male infertility. The testis is one of the most transcriptionally rich tissues; nevertheless, the specific gene expression changes occurring during spermatogenesis are not fully understood. To better understand gene expression during spermatogenesis, we generated germ cell-specific whole transcriptome profiles by systematically comparing testicular transcriptomes from tissues in which spermatogenesis is arrested at successive steps of germ cell differentiation. In these comparisons, we found thousands of differentially expressed genes between successive germ cell types of infertility patients. We demonstrate our analyses' potential to identify novel highly germ cell-specific markers (TSPY4 and LUZP4 for spermatogonia; HMGB4 for round spermatids) and identified putatively misregulated genes in male infertility (RWDD2A, CCDC183, CNNM1, SERF1B). Apart from these, we found thousands of genes showing germ cell-specific isoforms (including SOX15, SPATA4, SYCP3, MKI67). Our approach and dataset can help elucidate genetic and transcriptional causes for male infertility.

Human spermatogonial stem cells and their niche in male (in)fertility: novel concepts from single-cell RNA-sequencing

The amount of single-cell RNA-sequencing (scRNA-seq) data produced in the field of human male reproduction has steadily increased. Transcriptional profiles of thousands of testicular cells have been generated covering the human neonatal, prepubertal, pubertal and adult period as well as different types of male infertility; the latter include non-obstructive azoospermia, cryptozoospermia, Klinefelter syndrome and azoospermia factor deletions. In this review, we provide an overview of transcriptional changes in different testicular subpopulations during postnatal development and in cases of male infertility. Moreover, we review novel concepts regarding the existence of spermatogonial and somatic cell subtypes as well as their crosstalk and provide corresponding marker genes to facilitate their identification. We discuss the potential clinical implications of scRNA-seq findings, the need for spatial information and the necessity to corroborate findings by exploring other levels of regulation, including at the epigenetic or protein level.

Application of single-cell RNA sequencing on human testicular samples: a comprehensive review

So far there has been no comprehensive review using systematic literature search strategies to show the application of single-cell RNA sequencing (scRNA-seq) in the human testis of the whole life cycle (from embryos to aging males). Here, we summarized the application of scRNA-seq analyses on various human testicular biological samples. A systematic search was conducted in PubMed and Gene Expression Omnibus (GEO), focusing on English researches published after 2009. Articles related to GEO data-series were also retrieved in PubMed or BioRxiv. 81 full-length studies were finally included in the review. ScRNA-seq has been widely used on different human testicular samples with various library strategies, and new cell subtypes such as State 0 spermatogonial stem cells (SSC) and stage_a/b/c Sertoli cells (SC) were identified. For the development of normal testes, scRNA-seq-based evidence showed dynamic transcriptional changes of both germ cells and somatic cells from embryos to adults. And dysregulated metabolic signaling or hedgehog signaling were revealed by scRNA-seq in aged SC or Leydig cells (LC), respectively. For infertile males, scRNA-seq studies revealed profound changes of testes, such as the increased proportion of immature SC/LC of Klinefelter syndrome, the somatic immaturity and altered germline autophagy of patients with non-obstructive azoospermia, and the repressed differentiation of SSC in trans-females receiving testosterone inhibition therapy. Besides, the re-analyzing of public scRNA-seq data made further discoveries such as the potential vulnerability of testicular SARS-CoV-2 infection, and both evolutionary conservatism and divergence among species. ScRNA-seq analyses would unveil mechanisms of testes' development and changes so as to help developing novel treatments for male infertility.

The Molecular Signature of Human Testicular Peritubular Cells Revealed by Single-Cell Analysis

Peritubular cells of the human testis form a small compartment surrounding the seminiferous tubules. They are crucial for sperm transport, and they emerge as contributors to the spermatogonial stem cell niche. They are among the least known cell types of the human body. We employed single-cell RNA sequencing of cultured human testicular peritubular cells (HTPCs), which had been isolated from testicular samples of donors with normal spermatogenesis. The significant overlap between our results and recently published ex vivo data indicates that HTPCs are a highly adequate cellular model to define and study these cells. Thus, based on the expression of several markers, HTPCs can be classified as testicular smooth muscle cells. Small differences between the in vivo/in vitro expressed genes may be due to cellular plasticity. Plasticity was also shown upon addition of FCS to the culture medium. Based on transcriptome similarities, four cellular states were identified. Further analyses confirmed the presence of known stem cell niche-relevant factors (e.g., GDNF) and identified unknown functions, e.g., the ability to produce retinoic acid. Therefore, HTPCs allow us to define the signature(s) and delineate the functions of human testicular peritubular cells. The data may also serve as a resource for future studies to better understand male (in)fertility.

Mapping Human Reproduction with Single-Cell Genomics

The trillions of cells in the human body develop as a result of the fusion of two extremely specialized cells: an oocyte and a sperm. This process is essential for the continuation of our species, as it ensures that parental genetic information is mixed and passed on from generation to generation. In addition to producing oocytes, the female reproductive system must provide the environment for the appropriate development of the fetus until birth. New genomic and computational tools offer unique opportunities to study the tight spatiotemporal regulatory mechanisms that are required for the cycle of human reproduction. This review explores how single-cell technologies have been used to build cellular atlases of the human reproductive system across the life span and how these maps have proven useful to better understand reproductive pathologies and dissect the heterogeneity of in vitro model systems. Expected final online publication date for the Annual Review of Genomics and Human Genetics, Volume 23 is October 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Single-cell analysis of human testis aging and correlation with elevated body mass index

Aging men display reduced reproductive health; however, testis aging is poorly understood at the molecular and genomic levels. Here, we utilized single-cell RNA-seq to profile over 44,000 cells from both young and older men and examined age-related changes in germline development and in the testicular somatic cells. Age-related changes in spermatogonial stem cells appeared modest, whereas age-related dysregulation of spermatogenesis and somatic cells ranged from moderate to severe. Altered pathways included signaling and inflammation in multiple cell types, metabolic signaling in Sertoli cells, hedgehog signaling and testosterone production in Leydig cells, cell death and growth in testicular peritubular cells, and possible developmental regression in both Leydig and peritubular cells. Remarkably, the extent of dysregulation correlated with body mass index in older but not in younger men. Collectively, we reveal candidate molecular mechanisms underlying the complex testicular changes conferred by aging and their possible exacerbation by concurrent chronic conditions such as obesity.

In vitro propagation of XXY human Klinefelter spermatogonial stem cells: A step towards new fertility opportunities

Klinefelter Syndrome (KS) is characterized by a masculine phenotype, supernumerary sex chromosomes (47, XXY), and impaired fertility due to loss of spermatogonial stem cells (SSCs). Early testicular cryopreservation could be an option for future fertility treatments in these patients, including SSCs transplantation or in vitro spermatogenesis. It is critically essential to adapt current in vitro SSCs propagation systems as a fertility option for KS patients. KS human testicular samples (13,15- and 17-year-old non-mosaic KS boys) were donated by patients enrolled in an experimental testicular tissue banking program. Testicular cells were isolated from cryopreserved tissue and propagated in long-term culture for 110 days. Cell-specific gene expression confirmed the presence of all four main cell types found in testes: Spermatogonia, Sertoli, Leydig, and Peritubular cells. A population of ZBTB16⁺ undifferentiated spermatogonia was identified throughout the culture using digital PCR. Flow cytometric analysis also detected an HLA^-/CD9⁺/CD49f⁺ population, indicating maintenance of a stem cell subpopulation among the spermatogonial cells. FISH staining for chromosomes X and Y showed most cells containing an XXY karyotype with a smaller number containing either XY or XX. Both XY and XX populations were able to be enriched by magnetic sorting for CD9 as a spermatogonia marker. Molecular karyotyping demonstrated genomic stability of the cultured cells, over time. Finally, single-cell RNAseq analysis confirmed transcription of ID4, TCN2, and NANOS 3 within a population of putative SSCs population. This is the first study showing successful isolation and long-term in vitro propagation of human KS testicular cells. These findings could inform the development of therapeutic fertility options for KS patients, either through in vitro spermatogenesis or transplantation of SSC, in vivo.