Abstract
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.
1 Introduction
During the past two decades, an alarming rise in infertility particularly in developed/developing countries has become a matter of great concern (1). Around 15% of couples suffering from infertility globally, out of which 50%, is exclusively caused by the male partner (2). Male infertility is a multifaceted pathological condition with diverse manifestations, ranging from total absence of testicular sperm to specific changes in sperm quality (3). Chromosomal aberrations (like Klinefelter syndrome 47XXY) or variance (like 9qh+ heteromorphism) and other genetic mutations contribute upto 15-20% of total male infertility; however around 30-50% of cases remain untreatable due to idiopathic in origin/nature (4–6). Therefore, the mechanistic details of testicular spermatogenesis need to be reexamined at the cellular and molecular level (7, 8).
Spermatogenesis is gonadotropin regulated, highly synchronized, developmental programme involving multi-step events like stem cell renewal/differentiation, epigenomic remodeling, and meiotic divisions for generating millions of spermatozoa inside testes (9–11). The continuity of this process throughout adult life is supported by an active spermatogonial stem cell (SSC) population that critically maintains the key balance between the capacity of self-renewal (to sustain the stem cell population) and differentiation [to produce spermatogonial progenitor cells (SPCs), which subsequently develop into mature sperm via meiotic divisions] (12–14). Therefore, the cellular and molecular landscape associated with the SSC micro-environment turns out to be critical for regulating male fertility. The rapid intrinsic speed (around 1000-1500 sperm per heartbeat in adult men) of spermatogenic progression significantly highlights the degree of orchestration governing the transcriptional dynamics of testicular cells (15–17). Figure 1 represents the histological architecture of the testis showing arrays of multiple developing germ cells (Gc) (from SSC to mature motile sperms) alongside with visual representation of the different stages and cellular interactions essential for spermatogenesis.
Figure 1
Functional and/or anatomical defects in testicular somatic and/or Gc lead to spermatogenic impairment (with a qualitative and quantitative decline in sperm count clinically defined as azoospermia, oligozoospermia, teratozoospermia, and asthenospermia, etc.) causing male sub-infertility/infertility (18). Azoospermia, (the most severe form where semen without sperm) are of two types: obstructive azoospermia (OA) and non-obstructive azoospermia (NOA). OA stems from normal sperm production but obstructed sperm delivery, while NOA, accounting for around 10% of male infertility, results from lack of spermatogenesis in the testis (19). Three specific clinical categories are described in NOA as primary testicular failure like Sertoli cell only syndrome (SCOS)/Gc aplasia, Gc maturational arrest and hypo-spermatogenesis (20, 21). Supplementary Table 1 highlights (A) clinical categories of semen quality, (B) etiologies and frequencies of male infertility and (C) the hallmarks of primary testicular failure.
Notably, multiple genetic complexities involving over 2,000 spermatogenic genes observed in different forms of male infertility having substantial variations in semen quality and testicular histology (5, 22). Over recent decades assisted reproductive techniques (ART) e.g.- in vitro fertilization (IVF) and intra-cytoplasmic sperm injection (ICSI) have successfully facilitated biological fatherhood for men with extremely low sperm counts (3). However, addressing the fundamental basis of poor sperm quality and quantity remains challenging due to a limited understanding of the intricate molecular and cellular interactions driving human sperm production (18, 23).
An impressive progress in deploying high-throughput sequencing technologies in biomedical research has been observed during the past decade (24). Within this scenario, both bulk and single-cell RNA-sequencing (Sc-RNA-seq) collectively have emerged as revolutionary tools, revealing the molecular complexities governing testicular physiology and elucidating mechanisms underlying testicular dysregulation/disorders (25). However, analyzing bulk RNAs from diverse cell types faces adverse challenges in detecting expression data in rare cell populations, possibly diluting the presence of uncommon transcripts. In contrast, Sc-RNA-seq generates distinct RNA fingerprints for each cell (at single-cell resolution) type, transforming our understanding of differential transcriptomic profiles from complex tissues like the testis and therefore found to be critically informative in the field of clinical endocrinology, especially in deciphering infertility (26, 27). Table 1 summarizes the advantages of the Sc-RNA-seq strategy over the conventional bulk transcriptomic approach.
Table 1
| Approach/ Methods | Aim | Protocol | Popular platforms | Merits | Limitations | Direct Application(s) | Key References |
|---|---|---|---|---|---|---|---|
| Bulk RNA- Seq | To investigate the global gene expression profile of an organ/tissue without discriminating different cell populations | Total RNAs extracted from the tissues, depletion of rRNAs, cDNA preparation, adaptor ligation, library preparation, amplification, sequencing and generation of FASTQ files | Illumina (HiSeq/ NovaSeq etc), Ion torrent, Pacific Bio- Sciences & Oxford Nanopore etc | Well-developed for both short and long reads, cost-effective high throughput technique | Single cell specific transcriptomic data are unachievable, noisy gene expression profile with limited accuracy and detection rate, no spatial information available | Typically used for evaluating tissue and age specific differential gene expression profiles in model (reference guided genome assembly) and non-model (Trinity based de novo assembly) organisms, Identifications of splice variants/ transcript isoforms etc | (24) |
| Single Cell RNA-seq | To investigate the gene transcription at the resolution of a single cell level obtained from a complex tissue/ organ | Single-cell suspension is prepared [either by crude enzymatic digestion /mechanical dissociation, or sorting cells labeled by transgenic reporters (GFP/RFP/YFP etc) or based on antibodies etc.] followed by RNA extraction and subsequent tagging of the target transcripts with cell-specific unique molecular identifiers (UMIs), next cDNA preparation and other protocol same as performed for Bulk RNA-Seq | Two different capture systems- aqueous droplet and micro-bead based 10X Genomics & Microfluidic chip-based Fluigidm C1 Sequencing- Same as for bulk RNA-seq Illumina (HiSeq/ NovaSeq etc), | Most efficient at the level of individual cells (capable of estimating more than 10,000 cells) | High cost, limited number of transcripts can be detected/ identified, no spatial information available | Characterization of multifaceted transcriptional dynamics in heterozygous cell population in complex organ/ tissue network | (28, 29) |
| Spatial transcriptiomics | To explore the global imaging of transcriptional profile of an organ/tissue in situ | Ex situ Sequencing: Freshly frozen tissue are cryo-sectioned using OTC and oligo-bead (each having a linker sequence, a spatial barcode, a UMI sequence, and a polyT tail) arrays are placed on the tissue-section slide and slide then transferred to a tube for cDNA synthesis and library preparation and subsequent sequencing. Data are aligned for gene expression matrix and visualized via cluster mapping by barcode location matrix. | 10x Visium, Slide-Seq, NanoString GeoMx, HDST etc. | Discover de novo RNA sequence, capable to detect splice-isoforms , SNVs, whole transcriptional data in intact tissue architecture and exact location , 3D visualization possible with multiple sections | Low (10 µm) resolution, cell-cell boundaries not defined, dead spaces between two beads, low detection efficiency particularly for transcripts with low abundance and lack of single cell specific data [however, individual spot (single measurement site) based data can be visualized at single cell resolution by deconvolution tool] | To investigate cellular heterogeneity with spatial information, critical for developmental gene expression histopathology, disease progression, tumour and stem cell micro- environment etc. | (30, 31) |
| Ex situ Sequencing: Involves direct read out (via in situ cDNA synthesis, rolling amplification and sequencing) the transcript sequence within the tissue. Data are aligned and visualized via reconstitution. | |||||||
| Hybridization based: Target transcripts are detected by hybridization of complementary fluorescent probes flowed by imaging and decoding the data for visualization. | MERFISH, SeqFISH+ , oligo FISSEQ, Vizgen, 10X Xenium, Nanostring CosMx etc | Higher sub- celluar resolution (400nm) , well defined cell-cell boundaries with proper segmentation, 3D imaging possible without slicing and single RNA molecule is also detectable | Restricted to only known targeted transcripts (on that basis probes are synthesized), Chances of Optical crowding, High cost |
Summarizes the essential aspects of different RNA-Seq strategies, including their aims, key features, popular platforms, merits, limitations, and direct applications.
The table is divided into three sections: Bulk RNA-Seq, Single Cell RNA-Seq, and Spatial Transcriptomics.
Despite being critically relevant, comprehensive reviews of recent literature revealing the dynamic landscape of single-cell testicular transcriptomics are currently limited (25, 32). The current study dissects the spermatogenic landscapes viz., SSC/SPC niche, and cellular transcriptional dynamics during the meiotic onset and progression. Finally, this review attempts to reveal the intricate tapestry of testicular development and disorders through next-generation sequencing technologies.
2 Methodologies adopted for Sc-RNA-seq
The Sc-RNA-seq generates high-resolution transcriptomic data that precisely identifies minor transcriptomic variations at a single-cell level from a complex organ (33, 34). This tool has significant implications for gonadal development and disorders [viz., infertility, Disorders/Differences of Sex Development (DSDs), and cancer] (26, 27). Sc-RNA-seq produces transcriptomic data employing methods like conventional cellular purification via enzyme-based crude mechanical isolation or transgenic reporter gene/antibody-based sorting of target cells by sophisticated magnetic-activated cell sorter (MACS) or fluorescence-activated cell sorter (FACS) columns (35). Figure 2 represents the fundamental working principles of bulk and Sc-RNA-seq.
Figure 2
Single-cell isolation and sequencing can be achieved at two different platforms, first droplet-based 10× Genomics Chromium where individual cells get encapsulated in aqueous droplets along with micro-beads loaded with barcoded primers (36) and second Fluigidm C1, where individual cells are physically captured on microfluidic chips (37). After generating Sc-RNA-seq libraries, 10x Genomics utilizes 3′ end-counting unique molecular identifiers (UMIs), a genetic barcode that ensures a distinct cell population is traced and categorized as a specific lineage within a particular tissue/organ (32). In contrast, Fluidigm C1 libraries are automatically generated having full-length transcriptomes with all possible mRNA isoforms/variants (due to alternative splicing/transcription start site, etc.) (37). Notably, a new method Smart-seq2 has been developed recently having enhanced yield and length of cDNA libraries with better reverse transcription, template switching, and pre-amplification, improving the accuracy of transcript detection at a cheaper cost (38). Figure 3 schematically illustrates the available Sc-RNA-seq platforms and downstream workflow for data analyses.
Figure 3
Sc-RNA-seq data are plotted in bi-dimensional space and interpreted by methods like t-Distributed Stochastic Neighbor Embedding (t-SNE) and/or Uniform Manifold Approximation and Projection (UMAP) (33, 34). Although, these methods clearly depict cellular heterogeneity, Sc-RNA-seq data demonstrate cellular condition in a non-functional state as the tissue architecture gets dissociated/disaggregated (39). Therefore, new technology has been developed that captures the spatial transcriptomic data from the tissues via individual-cell laser-capture micro-dissection or multiplexed in situ hybridization (39). However, this method also shows limited sensitivity and depends on prior information on types of cells and/or genes to be targeted/analyzed. Recent Slide-seq technology overcomes these limitations and is considered to be the most suitable for generating high-throughput spatial transcriptomic information at 10-μm resolution (40). Figure 4 depicts the workflow of four popular methodologies being deployed for spatial transcriptomics.
Figure 4
3 Fetal testicular morphogenesis
The ontology of testicular germ and somatic cell lineages are developmentally independent (15, 41, 42). The bipotential fetal gonads are formed at the gonadal ridge (GR) during embryonic life [embryonic days (E) 8.0-9.5 in mice and weeks of gestation (W) 3.5-4 in humans] and appear distinctly (by E9.5-10 and W5) with thickening on the ventral side of the mesonephros (43). However, during E6.5- 7.5 and W3 murine and human primordial gem cells (PGCs) emerged from posterior-epiblast and amnion respectively. After getting specified (at E7.5 in mice and W4 in humans), these PGCs start migrating towards the developing GR and finally colonize the nascent gonads by E10.5 in mice and W8 in humans (44). The gonadal PGCs stop proliferation, lose pluripotency to gain competence for meiotic commitment via upregulating RNA binding DAZL protein (45). During E11.5 in mice and W6-7 in humans these nascent gonads undergo sex determination (SD), triggered by Sry expression in pre-Sertoli cells of XY gonads (46). The SD event involves step-wise transcriptional cascades generating male (Sry, Sox9, Fgf9, Amh, Dmrt1 etc.) and female (Wnt4, Rspo1, Foxl2, Runx1, etc.) specific robust genetic programming (47). The hallmark of fetal testicular morphogenesis is completed with the formation of the testicular cord and associated vasculature pattern due to selective (XY specific) migration of the mesonephric endothelial cells by E12.5 in mice and W6.5-7 in humans (48, 49). Notably, any developmental defect in the testicular cord formation or abnormality in Sertoli cells (Sc)/Leydig cells (Lc) differentiation directly leads to spermatogenic impairment and infertility during adulthood. Male Gc fail to enter meiosis due to active degradation of Retinoic acid (RA) by CYP26B1 expressed in fetal murine testes from E12.5-13.5 and transformed into mitotically quiescent (G0 arrested) gonocytes (50, 51). During fetal life, E15.5 in mice and W8.5-9 in humans, tubular supportive Sc, PTc, and interstitial fetal Lc (FLc) all are well demarcated (43).
4 Exploring cellular heterogeneity in the testis
Testis is a highly complex organ showing remarkable cellular heterogeneity constituted by both germline and somatic cells. Multiple somatic cell types [like nourishing Sc within seminiferous tubules (ST), supportive PTc covering the ST, interstitial (in inter-tubular space) steroidogenic Lc, interstitial macrophages, etc.] and diverse arrays of developing Gc [including pre-meiotic SSCs/SPCs (2n), different stages of meiotic Gc like primary spermatocytes (2n), secondary spermatocytes (n), round and elongated spermatids (n) and post-meiotic spermatozoa (n) etc] within ST contribute the tissue–organization of testis (Figure 1). Notably, almost all possible cellular events [like mitosis, meiosis, cellular differentiation, apoptosis, epigenetic remodeling, histone–protamine exchange and meiotic sex chromosome inactivation (MSCI) etc.] are evident in testis (52). Table 2 summarizes distinct/key cell markers of individual testicular cells with respect to the appropriate functional significance regulating male fertility. We here have discussed individual major cell–types contributing to the tissue organization of adult testis.
Table 2
| S. No. | Specific Marker/ Enriched Transcript | Full Name | Specific Cell type | Key Function/ Role |
|---|---|---|---|---|
| 1. | ID4 | Inhibitor of Differentiation 4 | SSCs / SPCs and Undifferentiated Spermatogonia (Spg) | Typical SSC marker |
| 2. | GFRA1 | GDNF family receptor α-1 | SSCs/ SPCs and Undifferentiated Spg | Critical for SSC self-renewal and differentiation. |
| 3. | MAGEA4 | Melanoma Antigen Gene Family A, Member 4 | Undifferentiated and differentiated Spg | Essential for normal testicular development and function, with its expression pattern and function tightly regulated to ensure proper spermatogenesis and male fertility. |
| 4. | SYCP3 | Synaptonemal Complex Protein 3 | First meiotic prophase, including preleptotene, leptotene, and zygotene spermatocytes (Sct) | Assembly and stability of the synaptonemal complex, as well as for the progression of meiotic recombination and synapsis between homologous chromosomes. |
| 5. | DMC1 | DNA Meiotic Recombination Protein 1 | First meiotic prophase | Key component of the meiotic recombination machinery and is specifically involved in the repair of DNA double-strand breaks (DSBs) that occur during meiosis |
| 6. | PIWIL1 | Piwi Like RNA-Mediated Gene Silencing 1 | Leptotene, pachytene,and zygotene Sct | Maintenance of genomic integrity through transposon silencing |
| 7. | PGK2 | phosphoglycerate kinase | Spermatids | Critical for meiotic germ cell metabolism |
| 8. | ACR | Acrosin | Pachytene Sct to postmeiotic spermatid stage | Also known as the Acrosin gene, plays a crucial role in testicular development by encoding the Acrosin protein |
| 9. | GAPDHS | Glyceraldehyde-3-Phosphate Dehydrogenase, Testis-Specific | Primary spermatocytes | Essential for regulating energy metabolism and sperm maturation processes |
| 10. | PRM1 | Protamine 1 | Spermatid | Ensures the proper packaging and condensation of DNA within sperm cells, which is essential for male fertility and successful reproduction. |
| 11. | SRY | Sex-Determining Region Y | Sertoli cells | Master regulator of male sex determination |
| 12. | SOX9 | SRY (Sex Determining Region Y)-Box 9 | Sertoli cells | Critical role in promoting testicular morphogenesis Sertoli cell differentiation |
| 13. | WT1 | Wilms Tumor 1 | Sertoli cells, germ cells, and interstitial cells (Leydig cells) | Differentiation and maintenance of Sertoli cells. |
| 14. | DMRT1 | Doublesex and Mab-3 Related Transcription Factor 1 | Sertoli cells, Germ cells | Maintains testicular differentiation by suppressing Foxl2 |
| 15. | FGF9 | Fibroblast Growth Factor 9 | Sertoli cells | Critical role in promoting testicular morphogenesis Sertoli cell differentiation |
| 16. | ALDH1A1 | Retinaldehyde dehydrogenase 1A1 | Leydig cells | Catalyzes the conversion of retinaldehyde to retinoic acid, a form of vitamin A. Retinoic acid serves as a signaling molecule that regulates the expression of genes involved in germ cell differentiation and meiosis initiation. Additionally, ALDH1A1 has been implicated in protecting germ cells from oxidative stress, thereby promoting their survival and proper development within the testes. |
| 17. | AMH | Anti-Müllerian Hormone | Sertoli cells | Regression of the Müllerian ducts, preventing the development of female reproductive structures. |
| 18. | IGFBP5 | Insulin-Like Growth Factor Binding Protein 5 | Sertoli and Leydig cells | Bioavailability and activity of insulin-like growth factors (IGFs), which are involved in testicular development and spermatogenesis |
| 19. | INSL3 | Insulin-Like 3 | Sertoli cells | Descent of the testes during fetal development |
| 20. | RSPO1 | R-Spondin 1 | Leydig cells | Regulation of Wnt signaling pathways, which play important roles ovarian differentiation |
| 21. | KITL (KIT) | Kit Ligand (Stem Cell Factor) | Sertoli and Leydig cells | Induction of proliferation/ differentiation / survival of SSC ‘ Spermatogonial cells |
| 22. | NR5A1 (SF1) | Nuclear Receptor Subfamily 5 Group A Member 1 (Steroidogenic Factor 1) | Sertoli and Leydig cells | Upstream regulator of Sry, induces steroidogenesis, sexual differentiation. |
| 23. | BMPR2 | Bone Morphogenetic Protein Receptor Type 2 | Sertoli cells, Leydig cells, and germ cells | Regulation of cell growth, differentiation, and apoptosis. |
| 24. | AR | Androgen Receptor | Sertoli cells, Leydig cells, and germ cells | Binds with testosterone, critical for testicular development, spermatogenesis, and male sexual differentiation. |
| 25. | CYP17A1 | Cytochrome P450 Family 17 Subfamily A Member 1 | Leydig cells | It is involved in the biosynthesis of androgens, such as testosterone and dehydroepiandrosterone (DHEA). |
| 26. | HSD3B1 | Hydroxy-Delta-5-Steroid Dehydrogenase, 3 Beta- and Steroid Delta-Isomerase 1 | Leydig cells | Catalyzes the conversion of pregnenolone and 17-hydroxypregnenolone to dehydroepiandrosterone (DHEA) and androstenedione, respectively, which are precursors for testosterone biosynthesis. |
| 27. | LHCGR | Luteinizing Hormone/Chorionic Gonadotropin Receptor | Leydig cells and germ cells | Mediates the effects of luteinizing hormone (LH) on Leydig cells, stimulating testosterone production. |
| 28. | ACVR2A | Activin A Receptor Type 2A | Sertoli cells, Leydig cells, and germ cells | Mediates the effects of activin and other ligands on cell growth, differentiation, and apoptosis. |
| 29. | DHH | Desert Hedgehog Homolog | Sertoli cells | Critical role in testicular development, Sertoli cell differentiation, germ cell proliferation and spermatogenesis. |
| 30. | EFNB2 | Ephrin B2 | Sertoli cells | Involved in cell-cell signaling and adhesion processes during testicular development and spermatogenesis. |
| 31. | NOTCH4 | Notch Receptor 4 | Sertoli cells, Leydig cells, and germ cells | Critical in Leydig cell fate determination, differentiation, and proliferation. |
| 32. | PDGFB | Platelet-Derived Growth Factor Subunit B | Leydig cells | Involved in cell growth, differentiation, and angiogenesis. |
| 33. | CXCR4 | C-X-C Motif Chemokine Receptor 4 | Sertoli cells, Leydig cells, and germ cells | Homing of gonocytes/ SSC to basement membrane to establish SSC pool/ niche |
| 34. | GDNF | Glial Cell Line-Derived Neurotrophic Factor | Sertoli cells, and peritubular cells | Regulation of SSC self-renewal and differentiation. |
| 35. | GJA1 | Gap Junction Protein Alpha 1 | Sertoli cells and germ cells | Forms gap junctions between adjacent Sertoli cells to form BTB. |
| 36. | FGFR1/R2/R3 | Fibroblast Growth Factor Receptor 1/2/3 | Sertoli cells, Leydig cells, and germ cells | Mediates the effects of fibroblast growth factors (FGFs) on cell growth, differentiation, and survival. |
| 37. | BMP4/8b | Bone Morphogenetic Protein 4/8b | Sertoli cells | Critical for testicular development, PGC specification, germ cell differentiation, and spermatogenesis. |
| 38. | HMGA1 | High Mobility Group AT-Hook 1 | Sertoli cells, Leydig cells, and germ cells | Involved in chromatin remodeling, gene regulation, and cell differentiation. |
| 39. | SMAD1/5/7 | SMAD Family Member 1/5/7 | Sertoli cells, Leydig cells, and germ cells | Regulation of transforming growth factor-beta (TGF-β) signaling pathways, which play important roles in testicular development, germ cell differentiation, and spermatogenesis. |
| 40. | STRA8 | Stimulated by Retinoic Acid Gene 8 | Meiotic germ cells | Meiotic gatekeeper |
| 41. | DMC1 | DNA Meiotic Recombination Protein 1 | Meiotic germ cells | Involved in the repair of DNA double-strand breaks during meiotic recombination, ensuring proper chromosome segregation and genetic diversity in sperm cells. |
| 42. | RAD51 | RAD51 Recombinase | Germ cells | Involved in homologous recombination repair of DNA double-strand breaks during meiosis and DNA repair processes. |
| 43. | ADAMTS2 | ADAM Metallopeptidase with Thrombospondin Type 1 Motif 2 | Sertoli cells, Leydig cells, and germ cells | Involved in extracellular matrix remodeling and tissue homeostasis. |
| 44. | CSF1 | Colony Stimulating Factor 1 | Sertoli cells | Regulation of SSC and macrophage differentiation and function |
| 45. | CDKN2D | Cyclin-Dependent Kinase Inhibitor 2D | Sertoli cells, Leydig cells, and germ cells | Involved in the regulation of cell cycle progression and proliferation |
| 46. | DAZL | Deleted in Azoospermia Like | germ cells | Germ cell licensing or commitment towards meiosis |
| 47. | RHOX5 | Reproductive Homeobox 5 | Sertoli cells and germ cells | Critical Androgen responsive gene supports male fertility. |
| 48. | RHOX 10 | Reproductive Homeobox 10 | SSCs / SPCs and Undifferentiated Spg | Maker of SSC , important for self-renewal |
| 49. | SLC25A4 | Solute Carrier Family 25 Member 4 | Sertoli cells, Leydig cells, and germ cells | Involved in the transport of metabolites across mitochondrial membranes and energy metabolism. |
| 50. | FSHR | Follicle-Stimulating Hormone Receptor | Sertoli cells | Promote Proliferation of Sertoli cells and support pre-meiotic germ cell differentiation. |
Essential Genes Influencing Spermatogenesis Identified Across Multiple Single-Cell RNA Sequencing Research Studies.
4.1 The germ-line cell types
The Gc residing within ST passes through a series of developmental stages [spermatogonia (2n), primary and secondary spermatocytes (n), round and elongated spermatids (n) etc] culminating in the production of mature sperm. Here we discuss different stages of developing Gc.
4.1.1 Primordial germ cells
Primordial germ cells (PGCs) represent the earliest identifiable precursors of Gc. PGCs originate via direct induction of amnion and/or posterior-epiblast (PE) cells from adjacent extra-embryonic ectoderm (EEE) and visceral endoderm (VE) by morphogens like bone morphogenetic proteins (Bmp4 and Bmp8b), Activin A and transcription factors like Blimp1 (Transcriptional repressor B lymphocyte induced maturation protein 1, also known as PR domain-containing protein 1, Prdm1), Prdm14 (PR domain zinc-finger protein 14) and AP2γ (activating enhancer-binding protein 2γ or Tfap2c) and Sox17 (only in humans) etc. (44). The PGCs migrate to colonize the developing GR and subsequently gain the meiotic commitment (lose pluripotency by upregulating DAZL) by E12.5 in mice and W8-9 in humans (50).
4.1.2 SSCs
SSCs are established from gonocytes post-birth (postnatal age 2-3 days in mice) which guarantee the continuity of sperm production throughout adult life (53). In rodents and primates, the duration of spermatogenesis varies among species, (in days; man 64, rhesus monkey 42, mouse 35, and rat 52) typically spanning 35 to 70 days (10, 54). In rodents As (A single spermatogonia) and in primates Apale (proliferative stem/progenitor pool) or Adark (slow diving reserve stem pool) are considered to be the functional SSC populations respectively (13, 14). The key balance between self-renewal and differentiation of these cells remains critical for spermatogenic maintenance throughout adult life (13). The Sc and PTc contribute to establishing the micro-environment of SSC that governs such regulated balance (55).
4.1.3 SPCs
SPCs are constituted with undifferentiated Spermatogonia (Spg), type Apaired or Aaligned connected with common cytoplasmic bridges. In rodents, these pre-meiotic (mitotic) Gc exhibit clonal fragmentation and get differentiated into a transit amplifying (TA) population (A1 to A4) induced by RA and then further transformed to differentiated Spg type B stage having meiotic fate (56). In primates, no such TA population is observed however, in rhesus monkeys (not human) 4 categories of Spg B (B1,B2,B3,B4) are found. In rodents, RA further drives the meiotic entry of Spg B by inducing Stra8, Rec8 and Meiosin expressions in these cells (50, 57).
4.1.4 Meiotic Gc
Meiotic Gc represent the most critical stage of spermatogenesis (58) with two successive divisions. In the first division, diploid (2n) primary spermatocytes transform into haploid (n) secondary spermatocytes followed by in the second division, these spermatocytes further divide to produce round spermatids. Meiotic Gc are pivotal in reducing the chromosome number by half, ensuring that the resulting sperm cells possess the correct genetic material for fertilization upon union with an egg cell. The juxtacrine testosterone (T) signaling through Sc regulates the meiotic progression (59).
4.1.5 Post-meiotic Gc
Post-meiotic Gc maturation mainly involves the transformation of a round spermatid to an elongated structure and such elongated spermatids get condensed and further mature into spermatozoa. In this step, the histone–protamine transition is observed leading to massive hetero-chromatinization leading to transcriptional silencing. However, the functional maturation of male gamete occurs with spermiogenesis, which involves complex cellular changes, viz., development of acrosome and flagellum, etc. Unraveling this metamorphosis is crucial for addressing male fertility concerns and developing strategies for male reproductive health and contraception (60). Figure 5 elucidates the testicular transcriptional dynamics concerning complex cellular heterogeneity. Figure 5A exhibits the differential degrees of heterogeneous/diverse gene expression profiles observed in multiple developmental stages of human (i) spermatogonia, (ii) spermatocytes, (iii) spermatids, and (iv) the entire adult human testis. Figure 5B is the graphical representation of developing Gc-specific marker expression during the entire spermatogenic progression. Figure 6 illustrates the comparative heat maps showing the expression profiles of critical genes involved/associated with different developmental/maturational stages of (A) spermatogonia (B) spermatocytes (C) spermatids and (D) whole adult testis.
Figure 5
Figure 6
Notably, Dazl (essential for gametogenic competence to Gc) expression was found to be consistent in all stages of spermatogonial development [like SSCs, progenitor spermatogonia (Spg)-1,2, early differentiating Spg-1,2 and late differentiating Spg] (Figure 6A). However, Dazl expression gets upregulated in high leptotene-zygotene spermatocytes (Sc1-1); declined in subsequent pachytene stage. Essential meiotic genes like Dmc1, Scml1 and Rad51 were found to be elevated in pre-leptotene stage, whereas, RA synthesizing genes like Aldh1a1 was highly enriched in somatic Sc/PTc. Finally, Dmrt1, a critical regulator testicular development was found to get augmented in early differentiating Spg-2, pre-leptotene Sc1-1 and late round spermatid (Std1) (Figure 6).
4.2 The somatic cells
Five distinct somatic cell types (that undergo mitosis but not meiosis) provide the critical structural and nutritional support are discussed in the following section.
4.2.1 Sertoli cells
Sertoli cells (Sc) are the major somatic cells that are the target of both FSH and T, known to dictate spermatogenesis by providing direct structural and nutritional support to all stages of developing Gc within the ST. Sc nurture and regulate the meiotic progression and completion of Gc. Sc-Sc tight junctions establish the blood-testis barrier (BTB) essential for Gc development and their directional movement towards the tubular lumen (61). FSH is critical for Sc proliferation and expansion of the pre-meiotic Gc population, whereas T governs the BTB formation, dynamics and Gc meiosis (4, 62). Sc derived glial cell line-derived neurotrophic factor (GDNF) and stem cell factor (SCF) play crucial role in determining the SSC fate decision from self-renewal to differentiation. Fetal Sc originate from the coelomic epithelium (63) and get specified (at E11.5 in mice and W6-7 in humans) by Sry driven Sox9 expression (43).The development of Sc is distinctly divided into three phases, immature (during fetal/neonatal/perinatal phase), maturing (during juvenile/pre-pubertal life) and mature (at puberty and adulthood) (64–66). Hormonal (FSH & T) stimulations to fetal/neonatal/infant Sc, show limited responses like self-proliferation and local expansion of pre-meiotic Gc but are found to be insufficient to induce robust onset of spermatogonial differentiation. The functional immaturity in terms of hormonal responsiveness of younger Sc is considered the underlying cause. Restricted plasma membrane localization of FSH-Receptor (Fsh-r), inadequate expression of associated Gαs sub-units and poor binding of T with Androgen Receptor (Ar) have been demonstrated in immature Sc during infancy (67–70). Therefore, the functional maturation of pubertal Sc is found to be developmentally critical for inducing Gc differentiation. Such pubertal maturational event in Sc become a prerequisite for spermatogenic onset and male fertility (15, 65). This is further evident by the presence of persistent immature Sc in adult testes that results in severe oligozoospermic conditions with infertility in rodents (71–75) and in humans (64, 76–79).
4.2.2 Peritubular cells
Peritubular cells (PTc) constitute the outermost covering of the ST and further provides structural support by creating peristaltic waves through their contractile elements, potentially prompting fluid movement within the tubule lumen to facilitate the expulsion of spermatozoa (80). PTc are the targets of T signaling and shown to promote Sc maturation. Furthermore, crosstalk between PTc with adjacent Lc play a vital role in spermatogenic development. Along with Sc, these cells govern the maintenance of the SSC micro-environment (81, 82).
4.2.3 Leydig cells
Leydig cells (Lc) are the key steroidogenic cells located within the inter-tubular interstitium and the sole target of LH for T biosynthesis. Lc are critical for male reproductive development viz., differentiation of the genital tract, fetal virilization, pubertal maturation and regulation of fertility during adulthood (82, 83). Two [fetal (FLc) and adult (ALc)] and four (fetal, neonatal, pubertal, and adult) distinct populations of Lc are reported in mice and primates with differential morphology and functions respectively (41). In mice, FLc exhibit dual origins (both coelomic epithelium and notch-active Nestin-positive perivascular cells located at the gonad–mesonephros borders) and first detected as Nr5a1 (also known as Ad4BP/SF-1) positive cells on E12.5 (W9-10 in humans). Murine FLc lack HSD17β3 enzyme and thereby produce only androstenedione (precursor of T) and promote the initial virilization and differentiation of male genitalia. Post-birth by P (post-natal age in days) 7-12, FLc undergo periodic regression and replaced by T-producing ALc during P18-25. Nestin-positive perivascular cells and FLc are the progenitor populations for ALc (84). Recent research indicates that 5-20% of FLc are maintained in adult testis (37). T acts on ALc, PTc and Sc to promote spermatogenesis, however, Gc do not express Ar. Sc shown to be the most critical target of T signaling for regulating male fertility (85). Although the extent of LH responsiveness of ALc remain unaltered with testicular aging, overall T production declines with testicular aging and testicular metabolic complications (86–88).
4.2.4 Macrophages
Macrophages, (being CSF1R+, CD206+ MHCII-) are present in the testicular interstitium and regulate the steroidogenic activity of ALc (84, 89). Another distinct population of adult macrophages (CSF1R-, CD206-, MHCII+) has been observed at various patches of the peritubular region which actively contributes to the SSC niche and promotes spermatogonial differentiation (80, 90). Intriguingly in fetal life, yolk-sac-derived murine macrophages (at E7.5) play a critical role in directing the selective migration of mesonephric endothelial cells towards developing XY gonads (not XX) and directs testicular cord (at E12.5) formation (49). Recent findings indicate that monocytes derived from hematopoietic stem cells (HSCs) of fetal liver get colonized in embryonic mouse testes in Sc dependent manner and further differentiate into testis resident macrophages during postnatal life (84). During adulthood, CD206+ interstitial macrophages promote ALc proliferation and steroidogenesis and thereby indirectly regulate male fertility. Furthermore, these macrophages also maintain the immunosuppressive environment within the testis (91).
4.2.5 Endothelial cells
Endothelial cells (ECs) form the border walls of blood vessels, facilitating nutrient exchange for developing Gc/sperm. Dysfunction in these cells can affect sperm development and overall testicular health, impacting fertility. During murine testicular morphogenesis, a key feature involves sex-specific vascularization where endothelial cells migrate from the neighboring mesonephros into the fetal XY gonads (during E11-12.5), encircling Sc-Gc clusters and prompting the formation of seminiferous cords. These cells are potent sources of vascular endothelial growth factor A (Vegf A) critical for tissue remodeling/homeostasis (92). ECs also contribute to the SSC niche by producing GDNF that sustains SSCs during extended culture periods (93).
Notably, moderate expressions of Dmrt1, Sox9, Wt1, Amh, Dhh and Aldh1a1 were found in murine Sc from E18.5 to P2D upto P7D. However, the transcript levels of the genes essential for induction of Gc differentiation like Scf/Kit-ligand, Gdnf and Csf-1 etc, remained low during this time of testicular development. In consistent with previous reports, Lc specific Insl3, Cyp17a1, Hsd3b1, Lh/cg-receptor and PTc restricted Ar expressions were observed from E18.5 to P7D; whereas meiotic gatekeeper Stra8 expression was enriched at Gc on P7D (Figure 7C). Figures 7A, B show age-dependent dynamic gene expression profiles observed in different cell types during three [(i) E18.5, (ii) P2D (post-natal 2 days of age) & (iii) P7D (post-natal 7 days of age)] distinct developmental phases of murine testicular maturation. Figure 7C depicts the comparative heat maps revealing the expression profiles of key spermatogenic genes in different cell types during (i) fetal (E18.5) testicular developmental and neonatal [(ii) P2D, (iii) P7D] testis maturation.
Figure 7
5 A comprehensive review of literature on testicular Sc-cell-RNA-seq revealed the molecular basis of testicular dysfunctions/disorders
The literature search was conducted during December – 2023- February 2024. Initially, systematic inquiries were performed in the PubMed database. Given the review’s primary focus on “single-cell RNA sequencing” and “Testis”, a tailored search strategy was applied to connect these two aspects. In detail, concerning the former aspect, the terms “single-cell” (term A) or “Sc-RNA-seq” (term B) were designated. Each of these terms was paired with one of the subsequent testis-related (the latter aspect-related) or animal model-related search terms, which included “testis” (term 1), “spermatogenesis” (term 2), “testicular development” (term 3), “azoospermia” (term 4), “Klinefelter Syndrome” (term 5), “cryptozoospermia” (term 6), “oligozoospermia” (term 7), “asthenospermia” (term 8), “teratospermia” (term 9), “orchitis” (term 10), “cryptorchidism” (term 11), “male and gonad” (term 12), “Humans” (term 13), “Monkey” (term 14), “Primates” (term 15), “mouse” (term 16), “Rat” (term 17), “Testicular diseases” (term 18), and “Testicular atlas” (term 19).
To reduce the likelihood of overlooking pertinent studies, particularly preprint studies, an additional search is conducted on the GEO dataset (https://www.https://www.ncbi.nlm.nih.gov/geo/) employing combinations of refined above search terms. The searches yielded approximately ~132 records from the last 10 years (2015–2024). Among these, 60 articles contained original Sc-RNA-seq data (from at least one donor) related to human testicular samples, while the remaining 68 articles involved re-analyzed data from prior publications or the datasets from other primates, rodents, livestock animals, or non-mammalian species. Starting in 2015, Guo et al. pioneered Sc-RNA-seq on 233 individual male and female human PGCs from 15 embryos between 4 and 19 weeks of gestation (94), followed by two studies in 2017 that sequenced adult testicular cells (95, 96). Interestingly, despite being from different teams, both adult studies focused on spermatogonia. Since 2018, numerous studies have conducted Sc-RNA-seq on human testicular samples, encompassing both prenatal and postnatal stages, normal and abnormal testes. Figure 8 represents the timeline of literature published/available in PubMed. Table 3 is the chronological summary of original research articles published on Sc-RNA-seq data for normal testicular samples collected from (A) primates (human and monkeys) and (B) rodents. Table 4 discusses country-wise Sc-RNA-seq data on clinical human samples with multiple forms of impaired spermatogenesis/infertility reported so far. Notably, Tables 3, 4 will serve as a critical repositories/valuable resource for Sc-RNA-seq datasets specific to distinct testicular cell types generated to date.
Figure 8
Table 3
| S. No | Year | Species | Num ber of Donar(s) | Sample Target Group (specific cell types) | Number of Cells examined | Age of tissue | GEO Accession Number | Tool and Sequencing Platform | Major Findings/ Critical Outcome | References |
|---|---|---|---|---|---|---|---|---|---|---|
| 1 | 2015 | Human | 15 embryos from 9 donars | Primordial Germ Cells (PGCs) | 197-319 cells | Fetal gonads at 4 to 19 weeks of gestation | Tang method followed by Illumina HiSeq2000/25 00 sequencer for 100 bp or 150 bp paired- end sequencing. | Inactivated X chromosomes get reactivated in PGCs by 4 weeks of age, gene expression pattern found to be homogeneous during 4 and 11 weeks of embryonic age in PGCs and global erasure of DNA methylation gets completed by 10 to 11 weeks. | (94) | |
| 2 | 2017 | Human | 5 | SSEA4+ SSC and differentiating c- KIT+ spermatog onia | Single cell preparation 15,000- 20,000 out of which 175(92 filtered) analysed | Adult Unknow n | GSE92280 | Fluidigm C1 followed by Illumina HiSeq2500 | Sc-RNA-seq clustering analysis reveals four cellular/developmental states during SSC differentiation, e.g.- Cluster A for state 1 key markers Id4, Etv5 (major transcription factors) and Txnip (for inhibiting glucose update); Cluster B for state 2 having Fgfr1/2 , Bmpr etc; Cluster C for state 3 with Ndufa/b (metabolic mitochondria), Sohlh2, Nr5a1, Thgbr1 (differentiation) , lsm3/4/6 (RNA splicing) Cluster D for state 4 with Cdc45, Rec8, Fanca (cell cycle, DNA repair) etc | (95) |
| 3 | 2017 | Human | Spermatogonial arrest (n = 1) and normal spermatogenesis (n = 7) | Spermatogonia | 20 and 105 | Unknown | GSE91063 | Tang method/Shallow RNA-seq followed by Ion Torrent PGM | Prominent heterogeneous transcriptomic profiles of Oct4, Utf1, Magea4, Boll, Prm2, Ddx5, Tspy1, Eef1a1 and Ngn3 in spermatogonial cells either with spermatogenic arrest (20 cells) or in normal spermatogenesis (105 cells). | (96) |
| 4 | 2017 | Human | 29 (17 female and 12 male) embryos | Human fetal germ cells (FGCs) | 2,167 single FGCs for analysis 1204 (1068 filtered) | 21 weeks of gestation | GSE86146 | Modified smart-seq2 followed by Illumina HiSeq 2500 | Reciprocal BMP and Notch signaling network found between FGCs and supportive gonadal cells, male FGCs develop through stages of migration, mitosis, and cell-cycle arrest (without RA signaling and meiosis) | (97) |
| 5 | 2017 | Human | 6 (all are biological fathers) | Germ cell subtypes | 500 | 38-52 years | SRP069329 | Laser capture microdissection, 5-µm-thickness of fixed testis; Illumina HiSeq 2000 | This study has identified over 4000 genes including the expression of 110 RNA-binding proteins and 137 long non-coding RNAs, previously unknown in spermatogenesis. Minimal transcriptional changes found between quiescent Adark to active Apale spermatogonia. | (98) |
| 6 | 2018 | Human | 4 adults and 2 infants | Human infant and adult testis | ~650,000 out of which 7830 filtered | 17-25 years in young adults, 13 months for infants | GSE120508 | 10× Genomics Chromium Controller followed by Illumina HiSeq 2500 | This study reveals > 8000 differentially regulated genes and uniquely explored multiple transposable elements (TE) like LTR12C/D/E, SVA_D, AluYa5, LTR10A, and LTR40C elements, and long non-coding RNAs during multiple steps of spermatogenesis. Five distinct states of spermatogonia (including a novel early SSC state, termed State 0) with discrete transcriptional status and developmental plasticity has been reported. | (99) |
| 7 | 2018 | Human and mice | 1610 (2 adult normal men, 7 OA and 1 NOA) | Spermatogonia, meiotic germ cells and haploid spermatids Spermatogonia, meiotic germ cells and haploid spermatids | 4,651 and 7,134 spermatogenic cells from mice and men for analysis 36451 for 10x and 635 for C1 | Organ donor 34.3 ±7.2 years) than biopsy patients (42.4 ± 1.2years) | GSE108970, GSE108974, GSE108977, GSE109049, GSE109033, GSE109037 | Both 10x Genomics and Fluidigm C1 followed by Illumina NextSeq 500; Illumina HiSeq 3000 | This study reports Sc-RNA-seq of >62,000 spermatogenic cells from mice and humans showing phylogenetic resemblance in SSC fate regulation. The hepatic stellate cell activation pathway is shown to be associated with SSC fate. Unique 3-gene identifiers distinguish 11 spermatogenic cell types in mice and humans | (100) |
| 8 | 2018 | Human | 10 (2 Norma l men, 7 OA and 1 NOA) | Spermatogonia, meiotic germ cells, and haploid spermatids | Out of 3,243 testicular cells, 3,028 were filtered. Additionally, 2,854 cells with normal spermatogenesis and 174 cells from a non-obstructive azoospermic (NOA) donor were analyzed. | 2 normal adults 30 & 60 years; 7 OA (27,29, 34,39,41,43,44 years) 1 NOA (24 years) | GSE106487 | Modified Smart-seq2 followed by Illumina HiSeq 4000. | Total 7,378 expressed genes and 122,443 mRNA molecules in each individual cell are detected. Total 17 clusters have been reported out of which, 3 spermatogonia subtypes, 7 spermatocyte subtypes, and 4 spermatid subtypes. Stage dependent expressions of Hmga1, Piwil4, Tex29, Scml1 and Ccdc112 have been found with critical contribution of FGF and BMP pathways in SSC development. | (101) |
| 9 | 2019 | Human | 4 | Spermatogonial cells | 33,585 | 2 neonatal (day 2 and day 7 of age) and 2 fertile adults aged 37 and 42 years) | GSE124263 | 10x Genomics followed by Illumina HiSeq4000 | Neonatal and adult testicular “niche” cells and factors are defined. In adult testes, 4 clusters/ subsets of undifferentiated spermatogonia (SPG) with distinct markers [slow proliferative 2 populations of SPG -A (having Utf1+ subset -1, Nanos3+ subset - 2) and 2 populations mitotic SPG -B (Nanos3, L1td1, Asb9, Dmrt1, Tuba3d, Dn mt1, Calr + early diff-SPG B and Sohlh2 + Diff-SPG)] have been identified. | (102) |
| 10 | 2020 | Human | 4 | Somatic and Germ cell and pubertal testis | 19223(12854 filtered) | Juvenile (7-11 years) and pre-pubertal (13-14 years) testes. 2 Testoster one-suppress ed transgen der female testis (26,60 years) | GSE134144 | 10 x genomics followed by Illumina HiSeq2500 (Homo sapiens) Illumina NovaSeq 6000 (Homo sapiens) | This transcriptional atlas of the developing human testis provides major insights into developmental changes and key factors/pathways that accompany male puberty. Data show distinctive phases of germ cell differentiation, common progenitor for Leydig and peritubular myoid cells, two distinct transcriptional states of pre-pubertal Sertoli cells. | (103) |
| 11 | 2020 | Human/ Monkey/Mouse | 4 Human; 5 Monkey | Somatic and Germ cells of testis | Human (13,837); Monkey (21,574) | 20-40 years; 4-13 years | GSE142585 | Drop-seq followed by Illumina HiSeq 2500 (Homo sapiens) Illumina HiSeq 2500 (Macaca mulatta) | Multi-species analyses reveal conserved and divergent testicular transcriptional dynamics. Xenotransplantation showsTspan33 as a marker of SSC, 6 different sub-populations of spermatogonial cells are defined based on specific markers. Phylogenetic differences in ligand-receptor interaction in testis are defined. | (104) |
| 12 | 2020 | Cynomolgus | 4 (1 infant, 1 juveni le 2 pubertal) | Spermatogonial cells, meiotic germ cells | filtered data with 17,792 cells | 1 year old infant, 2 year old juvenile, and 4 years old pubertal | E-MTAB-8979 | 10x Genomic followed by Illumina HiSeq 4000 or NovaSeq 6000 using paired- end sequencing | This data reveals critical markers outlining SSC differentiation (Utf1, Uchl1 and Dmrt1). Genes associated with DNA damage response pathway has been examined throughout spermatogenic progression. Phylogenetic conserved features along with divergent multicopy and ampliconic gene content found in meiotic sex chromosome inactive (MSCI) for primates and rodents. | (105) |
| 13 | 2020 | Human | 2 | Somatic and germ cells of testis | 853 PGCs, 2,854 normal testicular cells and 228 Sertoli cells | Male embryo and testis from (NOA) man | Unknown | Modified STRT-seq | All testicular cells express ACE2 and Sertoli cells show highest expression level, with age ACE2 expression gets decreased | (106) |
| 14 | 2020 | Human | 29 | Somatic and germ cell of pubertal testis | Spermatogo nial cells | fertile men aged between 30 and 50 y | SE144085 | 10 x Genomics | Gdnf and BMP8b broadly support long-term spermatogonial (Spg) culture, while activin A selectively promotes differentiated Kit+ Spg cells. Inhibition of AKT pathway specifically support primitive human undifferentiated Spg. | (107) |
| 15 | 2020 | Human and mice | 02 | Germ cells | 2554 | Adult men with OA 40 and 45 years of age | Drop microfluidics system flowed by Illumina NextSeq 500/550 | Spermatogenic genes display lower rats of mutations with low diversity in the population. Genes, remain silent show diverge and rapid evolutionary rates.Germ;ine specific mutational signmature is generated by TCR following 3’-pyrimidine rule. | (108) | |
| 16 | 2020 | Human, Mice | 3 | PGCs/ Gonocytes | 16429 | Fetal testes 17–18 weeks of gestation | GSE153819, GSE86146, GSE124263, GSE117101. | 10x Genomics, followed by 2 × 150 paired- end sequencing on Illumina HiSeq 4000 or NovaSeq 6000 | This study presents in vitro protocol for germ-line specification. The h-iPS generated PGC-like cells further get reconstituted as M-prospermatogonia-like cells and T1 prospermatogonia-like cells and closely resemble with human T1- prospermatogonia in vivo exhibiting diminished proliferation. Dynamic and stage-specific regulation of transposable elements during prospermatogonial specification have been found. | (109) |
| 17 | 2021 | Human | 2 | PGCs and somatic cells | ∼32,500 cells | Fetal first trimester and neonatal 5 months old | GSE143356 GSE161617 | 10X genomics | A transcriptional cell atlas of the fetal and postnatal human testes. Fetal sematic cells originate from common progenitor by 7 weeks post fertilization and by 14 weeks, PGCs exit mitosis, downregulate pluripotent transcription factors, and strongly resemble the state 0 spermatogonia. | (110) |
| 18 | 2021 | Human, Mice | 2 | Spermatogonial cells | Slide-Seq | 3-10 months in mice and 25 and 32 years for human | PRJNA668433 | Slide-seq tools | Combining Slide-seq with targeted in situ RNA sequencing, this study demonstrates specific specific differences in the cellular compositions of spermatogonial microenvironment between mouse and human testes. Furthermore, In developing Gc Habp4 (hyaluronan binding protein 4) has been identified as a potential regulator of chromatin remodelling | (111) |
| 19 | 2022 | Human Mouse | 22 | Human gonadal and adjacent extra-gonadal tissues (female n = 33, male n = 22) Mouse fetal gonads | 347,709, 96,174 and 40,742 cells, 63,929 cells | First and second trimester s of gestation (6– 21 PCW) Embryon ic 10.5, 11.5 and 12.5 days | E- MTAB -10551 | 10× Genomics followed by Illumina HiSeq4000; Novaseq 6000 | In both species, sex determination involves upregulated Sry and sPax8s, a gonadal lineage located at the gonadal–mesonephric interface. In human testes, Siglec 15+ and Trem2+ fetal testicular macrophage populations are identified. This study provides a comprehensive spatiotemporal map of human and mouse gonadal differentiation, which can guide ex vivo spermatogenesis. They identified specific human regulatory programs governing the development of germ-line and somatic cell lineages, comparing these stages with equivalent ones in mice | (27) |
| 20 | 2022 | Human | 9 male embryos/ fetuses and 6 female embryos/ fetuses | PGCs, gonocytes and fetal somatic cells | 31006 for 10x 709 for modified STRT-seq | 6 - 23 weeks of gestation | HRA000344 | 10× Genomics/mod ified STRT-seq Followed by Illumina HiSeq 4000 | Clustering analyses of fetal gonadal tissues reveal several novel cell subsets, viz., Pou5f1+Sparc+ FGCs and Krt19+ somatic cells. BMP signalling found to be critical in cell and stage-manner and promotes the gonocyte-to-spermatogonium transition. Biosynthesis of T also get transferred from fetal Sc to ALc. | (112) |
| 21 | 2022 | Human/ Mouse | 12 (4 Young and /8 Old) | All testicular cells | 44,657 | (17-22 years); and (62-76 years) | GSE182786; PRJNA757777 | 10x Genomics, Followed by Illumina NovaSeq 6000 | Age dependent spermatogenic dysregulation found to be more prominent as compare to that found in SSC. Altered pathways includes common inflammatory and hedgehog signalling selective metabolic signalling in Sc, T production in ALC, apoptosis in PTc etc . The extent of such dysregulation found to be associated with obesity in older but not in younger men. | (113) |
| 22 | 2022 | Monkey | 3 | Whole Testicular tissues | 7,500 nuclei | Young (5-6 years); Old (18-21 years) | GSA (CRA007812) | 10× Genomics followed by NovaSeq 6000 (Illumina, 20012866) | This work depicts in-depth transcriptomic traits of testicular aging at single-cell resolution, Sc found to be the most susceptible to aging, with significant decline in Wt 1 expression. | (114) |
| 23 | 2023 | Human/ All primates/Mice | 2 (Human), 3 (Chimpanze e), 2(Bon obo), 2 (Gorilla), 1 (Gibb on), 2 (maca que), 2 (marmoset) | Whole adult testes | 97,521 single- nuclei | 28, 32 years (Human) 16,21,45 years (Chimpanzee), 15,36 years (Bonobo), 43,51 | s E- MTAB -11063; E- MTAB -11064 | 10x Genomics followed by NextSeq 550 (Illumina) | This study elucidates the phylogenetic analyses of spermatogenic transcriptional dynamics based on single-nucleus transcriptomic data obtained from testes across 11 mammalian species. Results indicate critical aspects of molecular evolution of spermatogenesis like conserved transcriptome differences between X- and Y-bearing spermatids and meiotic sex-chromosome inactivation (MSCI) process etc. | (115) |
| 24 | 2023 | Human/ Mouse | 4 norma l and 9 NOA patients | Meiotic germ cells | Normal (1,097); NOA1 (176 cells) and NOA2 (130 cells) | unknown | GSE235324 | Illumina HiSeq 4000 (humans); Illumina NovaSeq 6000 (mice) | This study provides a multi-omics-based spermatogenic landscape at single-cell resolution. Abnormal DNA hyper-methylation has been detected in leptotene spermatocytes of NOA patients. Functionally, functional interference in DNA demethylation affects meiotic recombination and fertility. | (116) |
| 25 | 2023 | Human | – | Parthenocytes (PGs) and Andandrocytes (AGs) | – | – | NA | HiSeq2500 and NovaSeq platforms | They showed transcriptomic profiles comparison between different stages of early development in AGs and PGs. PGs activate primitive genes until the 4-cell stage with increased methylation, challenging the belief that only early stages are totipotent and hinting at totipotency in cleavage-stage AGs, emphasizing paternal transcript significance. | (117) |
| 26 | 2023 | Baboon and Rhesus Macaqu es | 2 baboo ns, 2 macaq ues | SSC | 5000 cells | baboons (new- born and 26 months) macaque s (15 and 20 months) | GSE222105; | 10x Genomics followed by Illumina HiSeq2500 | Human spermatogonia form discrete groups, whereas baboon and rhesus spermatogonia are found to be less heterogeneous. A cross-species analysis reveals significant differences between primate and mouse SSCs. For example, genes coding for components and regulators of the actin cytoskeleton cortical for cell- adhesion are enriched in primate-specific SSC. In humans, both SSC and progenitor spermatogonia were found to be Adark, while Apale spermatogonia appears biased towards differentiation. | (118) |
| 27 | 2023 | Human/ Mouse | 4 | SSC | Not Applicable | 7–18 weeks old mice and Human males 21- 48 years | PRJNA668433) | Slide-Seq | PTN and EPHA1 are potential SSC niche factors niche and spatial alteration in ligand-receptor paring leads to diabetes-induced infertility | (119) |
| S.No | Year | Species | Sample Target Group or specific cell types | Number of cells examined | Age of the samples | GEO Accession Number | Platform | Major Findings and Critical Outcome | References |
|---|---|---|---|---|---|---|---|---|---|
| 1 | 2018 | Mouse/ Human | Adult spermatogon ia, spermatocyte s, spermatids, Steady-state Spermatogenic cells | 62,000 cells | 6-day postpartum (P6) and adult | GSE108970; GSE108977; GSE109033 | Illumina NextSeq 500; Illumina HiSeq 3000 | The Mammalian Spermatogenesis Single- Cell Transcriptome, from Spermatogonial Stem Cells to Spermatids. | (100) |
| 2 | 2018 | Mouse | 2,500 cells from 8-week adult mice testis | 2550 cells | 2 adult (8 weeks) testis | GSE104556 | Illumina HiSeq 2500 | They covered all stages comprehensively, offering an ideal resource for marker discovery and analysis of differentiation processes. This dataset serves as a reference for future studies in disrupted spermatogenesis using single-cell RNA sequencing. | (120, 121) |
| 3 | 2018 | Mouse | Adult testis | 35,000 cells | 7,8,9,20 weeks | GSE112393 | Illumina HiSeq 2500 | They utilize single-cell RNA sequencing to create a resource on mouse spermatogenesis which aims to unravel the diverse cell types within the adult testis, pinpoint factors influencing germ cell differentiation, and uncover various types of somatic cells involved in this process. | (122) |
| 4 | 2018 | Mouse | Embryonic stage somatic cells | 435 cells | Embryo day E10.5, E11.5, E12.5, E13.5, E16.5 | GSE97519 | Illumina HiSeq2000 | They have identified a progenitor cell population expressing single Nr5a1 before sex determination. These cells undergo temporal fate specification, exhibiting competence windows to differentiate initially into Sertoli cells or later into fetal Leydig cells. | (123) |
| 5 | 2019 | Mouse | Spermatogonia and Spermatocyte | – | Postnatal (6d, 14d, 18d, 25d, 30d); Adult (8 weeks) | GSE121904 | Illumina NextSeq; 500 | Dynamic transcriptome profiles within spermatogonial and spermatocyte populations during postnatal testis maturation revealed by single-cell sequencing | (124) |
| 6 | 2020 | Mouse | Germ and somatic cell development during the perinatal period | 8,916 cells | Embryonic stage 18.5, P2, P7 | GSE130593 | Illumina HiSeq 4000 | This study pinpoints specific signaling pathways between somatic and germ cells during the perinatal phase, serving as a valuable resource for understanding testicular cell development during this time. | (125) |
| 7 | 2020 | Mouse | 3-week mice injected with busulfan and rescued with AOS | 27,000 cells | Treatment 3 weeks of age to 8 weeks of age | GSE131629 | Illumina HiSeq 2000 | They concluded that alginate oligosaccharides improved blood and testis metabolomes as well as the gut microbiota to support the recovery of spermatogenesis | (126) |
| 8 | 2020 | Murine | ECs from 11 mouse tissues including testis | >32,000 single Endothelia l Cells (ECs) | 8-weeks male mice | Array Express E-MTAB- 8077 | Illumina HiSeq 4000 | Reported is a comprehensive murine atlas consisting of over 32,000 single ECs transcriptomes obtained from 11 mouse tissues including testis. Within these subclusters, diverse classical and tissue- specialized endothelial cell subtypes have been delineated. | (127) |
| 9 | 2021 | Mouse/ Human | – | – | 3 and 10 months testis; Two healthy men (25 and 32 years) | PRJNA668433 | Nextera XT; NextSeq 500/550 High and Mid Output Kit v2.5 (75 Cycles) | They generate a spatial transcriptome atlas of the mammalian testis, which is used to reveal the spatial organization of the testicular microenvironment and profile its changes under diabetic conditions. | (111) |
| 10 | 2021 | Mouse | PGCs | 14,750 | E10.5, W12.5, and E16.5 | GSE136220 | Illumina HiSeq 4000 (Mus musculus) | This study provides a molecular roadmap of GC sex determination at single-cell resolution and will serve as a valuable resource for future studies of gonad development, function and disease | (128) |
| 11 | 2021 | Mouse | cultured THY-1+ GCs w/wo FGF9 | 4,143 control and 4,482 FGF9-treated cells | – | – | NextSeq 500 (Illumina) | They concluded that FGF9 is an important regulator of SSC proliferation, operating through p38 MAPK phosphorylation and upregulating Etv5 and Bcl6b in turn. | (129) |
| 12 | 2022 | Mouse | Young and Old | 37,571 cells | 2 and 24 months | PRJCA00856 | 10X Genomics Chromium; NovaSeq 6000 | This study reveals the first detailed map of aging in mouse testes at the single-cell level. It identifies changes in gene profiles, disruptions in stem cell balance, and the emergence of specific aging-related cell types, like macrophages. The findings offer a vast resource for studying age-related subfertility and suggest potential paths for diagnostics and targeted treatments against testicular aging by focusing on oxidative stress and inflammation. | (130) |
| 13 | 2022 | Rat | Adult LCs before and after EDS treatment | 10,000 | 12 weeks | PRJCA006440 | Illumina NovaSeq 6000 | Identification of Rat Testicular Leydig Precursor Cells by Single-Cell-RNA- Sequence Analysis | (85) |
| 14 | 2022 | Rat | Testes from 3 animals were collected at each of 1, 3, and 7 weeks post-EDS treatment | 4,000 | 12 weeks | PRJCA006139 | Illumina HiSeq PE150 | The study explored the effects of the androgen environment on the regulation of spermatogenesis. As this is the first single- cell RNA-Seq dataset for rat testes, it can also serve as a reference for future studies | (131) |
| 15 | 2022 | Mouse | WT and Alkbh5 KO testis | WT- 5,596 cells; KO- 6,816 cells | 12 week | GSE190396 | HiSeq X Ten (Mus musculus) | This study presents the inaugural single-cell RNA sequencing profile of ALKBH5- deficient mice's testes. It underscores ALKBH5's significance in germ cell development and spermatogenesis, providing fresh molecular insights. These discoveries might form the groundwork for further investigations into the origins and therapies for male infertility. | (132) |
| 16 | 2023 | Mouse | Control and Sertoli cell specific Scf- cKO mice testis | – | 8-week-old | GSE161040 | HiSeq X Ten (Mus musculus) | Here authors utilised fluorescent reporter mice and revealed widespread expression of the stem cell factor (Scf) across different testicular stromal cells, boosting Scf specifically in Sertoli cells significantly improved spermatogenesis, underscoring their crucial role in regulating this process. | (133) |
| 17 | 2023 | Mouse | fetal to neonatal transition | 25,613 | E18.5; Postnatal (day1; day3; day6) | GSE164439 | Illumina HiSeq 2000 (Mus musculus) | The single-cell chromatin accessibility landscape in mouse perinatal testis development | (134) |
| 18 | 2023 | Mouse/H uman | Adult male mice; Mouse SSCs; 4 healthy men; Human SSCs | – | (7–18 weeks old) | PRJNA668433 | Illumina NovaSeq S2 flow cell 100 cycle kit | Their findings suggested that inferring ligand-receptor (LR) pairs within the spermatogonial stem cell niche identified PTN and EPHA1 as potential factors within this niche. They also indicated that the spatial changes in LR pair expression are the root cause of infertility induced by diabetes | (118) |
| 19 | 2023 | Mouse | Adult WT, GRTH- Knockout, and GRTH knock-in mutant mice | 7,000 cell for each sample | 8 weeks old | GSE221226 | Illumina NovaSeq 6000 | These studies highlight the significance of pGRTH in acrosome biogenesis and the progression of round spermatids RS into ES during spermiogenesis. | (135) |
| 20 | 2023 | Mouse/H uman | WT mice treated with busulfan | 10,000 cells | 8 weeks old | GSE164787 Reanalysed GSE124263 | Illumina NovaSeq 6000 | This study provides a blueprint for understanding the development of the male germline and supporting somatic cells in humans. The germ cell subset markers identified are candidates to be used for clinical applications, including SSC therapy for treating infertility. | (136) |
| 21 | 2023 | Mouse | adult Ddx43+/+, Ddx43KI/KI and Ddx43– /–testes | 17,133 | adult | PRJNA650016; PRJNA838233 | Illumina NovaSeq 6000 | Their identification of RNA helicase DDX43 as a crucial regulator of chromatin remodeling in spermiogenesis showcases its pivotal role in this process. These findings underscore the significance of DDX43 in spermiogenesis and emphasize the efficacy of the single-cell-based approach in delineating cell-state-specific control in male germline development. | (137) |
| 22 | 2023 | Mouse | control and Kdm6a conditional knockout (“cKO”) Whole testis | WT- 19,378 cells; KO- 16,740 cells | – | GSE215112 | Illumina NovaSeq 6000 | The histone demethylase KDM6A, exhibits transient expression just before and during meiotic entry in spermatogenesis, playing a crucial role in maintaining epigenetic states across generations in the male germline. | (138) |
Summarizing the key findings of mammalian testicular single-cell RNA-seq datasets. (A). Normal healthy Human/ Monkeys. (B). Exclusively Rodents.
Table 4
| S.No | Year | Population | Number of Donar(s) | Age | Pathology/ Etiology reported | Number of Cells examined | GEO Accession Number | Platform | Major Findings/ Critical Outcome | References |
|---|---|---|---|---|---|---|---|---|---|---|
| 1 | 2018 | Chinese | 2 Healthy donors, 7 Obstructiv e azoosperm ia (OA) donors, and 1 Non- obstructiv e azoosperm ia (NOA) donor | Healthy (30, 60 years); OA (39, 43, 27, 34, 44, 41, 29); and NOA (24 years) | SCOS found in one NOA donor | 2854 (normal) 174 (NOA) | GSE106487 | Modified Smart- seq2 followed by Illumina HiSeq 4000 | Enrichment of γH2AX signal in NOA somatic cells, suggests their activation of DNA damage response mechanisms. | (101) |
| 2 | 2019 | German | 1 donor | Unknown | 47,XXY Klinefelter Syndrome (KS) | 3289 | GSE130151 | 10x Genomic s; Illumina NovaSeq 6000 | On the transcriptional level, Klinefelter Gc (both spermatogonia and sperm) exhibit normal DNA methylation, however, Klinefelter patients show variations in the DNA methylation of imprinted regions. | (139) |
| 3 | 2020 | Danish | 1 KS donor | Unknown | KS | 3,289 | Reanalyzed GSE130151 | Illumina NovaSeq6000 | KS likely causes early testicular damage, leading to fewer Gc, loss in Sc, and ALc changes. Genetic studies pinpoint specific gene alterations, especially in X-escapee genes, driving these issues. Sc dysfunction is emphasized as pivotal in Gc loss, urging wider research across ages and more KS patients to understand this condition better. | (140) |
| 4 | 2020 | Chinese | 4 (2 OA and 2NOA) | Unknown | NOA | 1212 | OEP000778 | Singleron GEXSC OPE foll owed by Illumina HiSeq X Ten | Transcriptional difference in OA and NOA, Sc-Sc junction genes Syne2, Atp2b1, Mtdh and elevated Fate1. | (141) |
| 5 | 2020 | Chinese | 10 Healthy donors and 7 NOA donor | Infant to Adult | Testicular torsion , benign testicle mass , or contralateral testis to cryptorchidism NOA (SCOS) | 88,723 | GSE149512 | 10X Genomic s followed by Illumina NovaSeq 6000 | This study has compared testicular cells (from 10 normal and 7 NOA) and found prominent maturational defects. Blocking Wnt signaling helps immature cells to support germ cell survival, offering new diagnostic /treatment strategies for NOA. | (64) |
| 6 | 2021 | Chinese | 7 OA, 1 NOA and 2 normal persons | 34, 36, 50 years | Hypospermatoge nesis | 480 | GSE157421 | (STRT- seq) follo wed by Illumina HiSeq XTEN | Changes in autophagy-related genes, such as the upregulation of Sqstm1 and the down-regulation of Lc3a in spermatids, were identified in NOA Gc , while Cst3- mediated autophagy potentially aids in maintaining SSC through Oct4, Id1, and Nanos3; additionally, inhibition of Cst3 disrupted meiosis and spermatid formation. | (142) |
| 7 | 2021 | Italian | 3 iNOA donors and 1 OA donor | iNOA (32/41/37 years) and OA (37 years) | idiopathic SCOS/ germ cell aplasia (OA and iNOA) | 3880 | GSE154535 | 10X Genomic s followed by Illumina NextSeq5 00/Nova Seq 6000 | idiopathic Gc aplasia viz., immaturity of ALc, chronic tissue inflammation, fibrosis, and senescence phenotype of the testicular somatic cells, ratio of Lamin A/C transcripts and an active DNA damage response in Lc and PTc. | (143) |
| 8 | 2021 | German | 3 healthy and 3 cryptozoo spermic donor | 31, 33, 59 years healthy and 36, 39, 25 years Diseased | cryptozoospermi a | 30,000 | GSE153947 | 10X Genomic s followed by Illumina NovaSeq 6000 | Alterations in the crypto-spermatogonial compartment having more of Piwil4+ undifferentiated spermatogonia, prolonged expression of Egr4 and fewer SSC/ Adark spermatogonia. . | (144) |
| 9 | 2021 | Americ an | 12 donors | Mix age | KS and idiopathic azoospermic (iNOA) male infertility | 26300 | GSE169062 | 10x Genomic s followed by Illumina HiSeq 2500 | This study identifies a subpopulation of Sc (within multiple individuals having KS) lacking transcription from the XIST locus resulting elevated X-linked gene expressions. Furthermore, 72 pathways upregulated in KS suggesting changes in interstitial cells due to loss of X inactivation in Sc. | (145) |
| 10 | 2022 | Mix | – | -Do- | NOA | – | GSE149512; GSM4504195; GSM4504196; GSM4504197; GSM4504195; GSM4504196; GSM4504197 (Data Reanalyzed) | – | Elevated ALc and macrophages in iNOA patients testes Lc specific markers Lhx9, Klf8, Klf4, Arid5b and Rxrg) in NOA, and macrophages-specific TFs (such as Pou2f2, Spib, Irf5, Cebpa, Elk4 and Klf6) in NOA, | (146) |
List of Human Participants Exhibiting Spermatogenic Dysfunction or Infertility: A Detailed Examination.
6 New progress
This comprehensive review of literature on testicular Sc-RNA-seq highlights the complexity of testicular tissue composition and the uniqueness of cellular transcriptional dynamics. The cutting-edge approach has unveiled a heterogeneous landscape within the testis, providing insights into spermatogenic Gc, somatic PTc and interstitial Lc and supportive Sc associated with SSC microenvironments. Studies employing Sc-RNA-seq have examined the transcriptional profiles of various Gc populations at distinct developmental stages elucidating the molecular regulation of spermatogonial differentiation, meiosis, and sperm maturation. The application of Sc-RNA-seq has not only enriched our understanding of testicular cell heterogeneity but has also identified novel markers and regulatory pathways critical for male infertility and reproductive health (147). Integration of multi-omics data and advanced in-silico analyses offers a comprehensive view of the testicular transcriptomic landscape, fostering discoveries with profound implications for both basic biology and clinical applications in reproductive medicine.
A significant development in the field of testicular transcriptomics requires an inclusive assimilation/integration of data obtained from Sc-RNA-seq analyses with spatial transcriptomic information (32). For example, Chen and colleagues employed the “spatial transcriptomic analysis” technology, creating a spatial atlas that delineates testicular gene expression at near-single-cell resolution in the human testis (111). Similarly, Garcia-Alonso et al. conducted integrated analyses of male gonadal development by fusing Sc-RNA-seq and spatial transcriptomics, offering a more comprehensive and in-depth perspective (27).
7 Current challenges
The integration of these methodologies is anticipated to augment future investigations on the human testis and male infertility, affording a comprehensive understanding of functional/impaired spermatogenesis. Although Sc-RNA-seq analyses have corroborated the notable presence of immune cells (e.g., macrophages - both M1 and M2, T cells, mast cells, and B cells) in the testicular microenvironment, the functional roles (and also transcriptional patterns) of these cells in regulating male infertility development remain elusive (85). However, despite the current high cost of Sc-RNA-seq, several of the studies sincerely lack a direct clinical orientation Therefore, diverse testicular pathologies, viz., cryptorchidism, AZFc deletions, azoospermia post chemo/radiotherapy, and teratospermia, demand focused exploration of Sc-RNA-seq analyses. The ongoing refinement and cost mitigation of this technology holds future promise for clinical applications, such as diagnostics, classification, prediction of sperm retrieval, and evaluation of hormonal treatment efficacy. Figure 9 illustrates the probable applications of Sc-RNA-Seq data towards framing appropriate preclinical studies (both in vitro and in vivo approach) leading to potential future diagnostic (assays/tools/kits) in male reproductive health care.
Figure 9
8 Concluding remarks and future directions
Exploring the intricate landscape of testicular function and pathology requires a comprehensive understanding spanning from the fundamental processes of development to the complexities of the disorders. This review has examined the uncharted territory of the testicular single-cell transcriptional atlas, illuminating the diverse cellular populations and their orchestrated molecular dynamics. From the genesis of Gc to the intricate orchestration of somatic cell interactions, this exploration navigates the transcriptional intricacies shaping testicular development. Furthermore, it delves into the disruptive signatures underlying various testicular disorders, shedding light on potential avenues for diagnostic and therapeutic advancements. By unraveling the cellular intricacies from a single-cell perspective, this review aims to provide comprehensive insights bridging the developmental trajectory to the manifestations of disorders within the testicular microenvironment. We also summarize the previous attempts deploying Sc-RNA-Seq on clinical testicular biopsies (Table 4) that have significantly enriched our understanding of testicular physiology e.g.- regulation of SSC micro-environment and testicular somatic cell maturation, etc. Finally, the increasing availability of publicly accessible testicular Sc-RNA-seq datasets provides researchers with greater possibilities for conducting personalized data analyses to address diverse research requirements. We anticipate that this review will not only facilitate the comprehensive understanding of spermatogenic transcription during normal and impaired (sub-fertile/infertile) conditions but also serve as a valuable reference for critical diagnosis and cure for some forms of unexplained/idiopathic male infertility.
Statements
Author contributions
MT: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Writing – review & editing. IB: Visualization, Formal Analysis, Writing – original draft, Writing – review & editing. MM: Data curation, Formal analysis, Writing – review & editing. VB: Data curation, Formal analysis, Writing – review & editing. MC: Conceptualization, Data curation, Investigation, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
IB sincerely thanks the past and present Vice-Chancellor, Central University of Kerala, Kasaragod, Kerala, India for necessary administrative support.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2024.1394812/full#supplementary-material
References
1
InhornMCPatrizioP. Infertility around the globe: new thinking on gender, reproductive technologies and global movements in the 21st century. Hum Reprod Update. (2015) 21:411–26. doi: 10.1093/humupd/dmv016
2
EisenbergMLEstevesSCLambDJHotalingJMGiwercmanAHwangKet al. Male infertility. Nat Rev Dis Primers. (2023) 9:49. doi: 10.1038/s41572-023-00459-w
3
AgarwalABaskaranSParekhNChoCLHenkelRVijSet al. Male infertility. Lancet. (2021) 397:319–33. doi: 10.1016/S0140-6736(20)32667-2
4
BhattacharyaIDeySBanerjeeA. Revisiting the gonadotropic regulation of mammalian spermatogenesis: evolving lessons during the past decade. Front Endocrinol (Lausanne). (2023) 14:1110572. doi: 10.3389/fendo.2023.1110572
5
KrauszCRiera-EscamillaA. Genetics of male infertility. Nat Rev Urol. (2018) 15:369–84. doi: 10.1038/s41585-018-0003-3
6
TournayeHKrauszCOatesRD. Novel concepts in the aetiology of male reproductive impairment. Lancet Diabetes Endocrinol. (2017) 5:544–53. doi: 10.1016/S2213-8587(16)30040-7
7
CannarellaRCondorelliRAMongioìLMLa VigneraSCalogeroAE. Molecular biology of spermatogenesis: novel targets of apparently idiopathic male infertility. Int J Mol Sci. (2020) 21:1728. doi: 10.3390/ijms21051728
8
MajumdarSSBhattacharyaI. Genomic and post-genomic leads toward regulation of spermatogenesis. Prog Biophys Mol Biol. (2013) 113:409–22. doi: 10.1016/j.pbiomolbio.2013.01.002
9
RussellLDEttlinRAHikimAPCleggED. Mammalian spermatogenesis. In: RussellLDEttlinRSinha HikimAPCleggED, editors. Histological and histopathological evaluation of the testis. Cache River, Florida (1990). p. 1–40.
10
SharpeRM. Regulation of spermatogenesis. The physiology of reproduction Vol. 1. New York: Raven Press (1994) p. 1363–434.
11
BhattacharyaISharmaSSMajumdarSS. Etiology of male infertility: an update. Reprod Sci. (2023) 1(4):942–65. doi: 10.1007/s43032-023-01401-x
12
SongHWWilkinsonMF. In vitro spermatogenesis: A long journey to get tails. Spermatogenesis. (2012) 2:238–44. doi: 10.4161/spmg.22069
13
FayomiAPOrwigKE. Spermatogonial stem cells and spermatogenesis in mice, monkeys, and men. Stem Cell Res. (2018) 29:207–14. doi: 10.1016/j.scr.2018.04.009
14
OatleyJMBrinsterRL. The germline stem cell niche unit in mammalian testes. Physiol Rev. (2012) 92:577–95. doi: 10.1152/physrev.00025.2011
15
BhattacharyaISen SharmaSMajumdarSS. Pubertal orchestration of hormones and testis in primates. Mol Reprod Dev. (2019) 86:1505–30. doi: 10.1002/mrd.23246
16
RabbaniMZhengXManskeGLVargoAShamiANLiJZet al. Decoding the spermatogenesis program: new insights from transcriptomic analyses. Annu Rev Genet. (2022) 56:339–68. doi: 10.1146/annurev-genet-080320-040045
17
GriswoldMD. Spermatogenesis: the commitment to meiosis. Physiol Rev. (2016) 96:1–17. doi: 10.1152/physrev.00013.2015
18
TüttelmannFRuckertCRöpkeA. Disorders of spermatogenesis: Perspectives for novel genetic diagnostics after 20 years of unchanged routine. Med Genet. (2018) 30:12–20. doi: 10.1007/s11825-018-0181-7
19
ZhankinaRBaghbanNAskarovMSaipiyevaDIbragimovAKadirovaBet al. Mesenchymal stromal/stem cells and their exosomes for restoration of spermatogenesis in non-obstructive azoospermia: a systemic review. Stem Cell Res Ther. (2021) 12:229. doi: 10.1186/s13287-021-02295-9
20
CioppiFRostaVKrauszC. Genetics of azoospermia. Int J Mol Sci. (2021) 22:3264. doi: 10.3390/ijms22063264
21
BehreHMBergmannMSimoniMTüttelmannF. Primary testicular failure. 2015 aug 30. In: FeingoldKRAnawaltBBlackmanMRBoyceAChrousosGCorpasEde HerderWWDhatariyaKDunganKHoflandJKalraSKaltsasGKapoorNKochCKoppPKorbonitsMKovacsCSKuohungWLaferrèreBLevyMMcGeeEAMcLachlanRNewMPurnellJSahayRShahASSingerFSperlingMAStratakisCATrenceDLWilsonDP, editors. MDText.com, Inc, South Dartmouth (MA (2000).
22
MatzukMMLambDJ. The biology of infertility: research advances and clinical challenges. Nat Med. (2008) 14:1197–213. doi: 10.1038/nm.f.1895
23
NetoFTLBachPVNajariBBLiPSGoldsteinM. Spermatogenesis in humans and its affecting factors. Semin Cell Dev Biol. (2016) 59:10–26. doi: 10.1016/j.semcdb.2016.04.009
24
StarkRGrzelakMHadfieldJ. RNA sequencing: the teenage years. Nat Rev Genet. (2019) 20:631–56. doi: 10.1038/s41576-019-0150-2
25
SuzukiSDiazVDHermannBP. What has single-cell RNA-seq taught us about mammalian spermatogenesis? Biol Reprod. (2019) 101:617–34. doi: 10.1093/biolre/ioz088
26
EstermannMASmithCA. Applying single-cell analysis to gonadogenesis and DSDs (Disorders/differences of sex development). Int J Mol Sci. (2020) 21:6614. doi: 10.3390/ijms21186614
27
Garcia-AlonsoLLorenziVMazzeoCIAlves-LopesJPRobertsKSancho-SerraCet al. Single-cell roadmap of human gonadal development. Nature. (2022) 607:540–7. doi: 10.1038/s41586-022-04918-4
28
JovicDLiangXZengHLinLXuFLuoY. Single-cell RNA sequencing technologies and applications: A brief overview. Clin Transl Med. (2022) 12:e694. doi: 10.1002/ctm2.694
29
ZhangSLiXLinJLinQWongKC. Review of single-cell RNA-seq data clustering for cell-type identification and characterization. RNA. (2023) 29:517–30. doi: 10.1261/rna.078965.121
30
RaoABarkleyDFrançaGSYanaiI. Exploring tissue architecture using spatial transcriptomics. Nature. (2021) 596:211–20. doi: 10.1038/s41586-021-03634-9
31
ZhangYLinXYaoZSunDLinXWangXet al. Deconvolution algorithms for inference of the cell-type composition of the spatial transcriptome. Comput Struct Biotechnol J. (2022) 21:176–84. doi: 10.1016/j.csbj.2022.12.001
32
DongFPingPMaYChenXF. Application of single-cell RNA sequencing on human testicular samples: a comprehensive review. Int J Biol Sci. (2023) 19:2167–97. doi: 10.7150/ijbs.82191
33
ChengCChenWJinHChenX. A review of single-cell RNA-seq annotation, integration, and cell-cell communication. Cells. (2023) 12:1970. doi: 10.3390/cells12151970
34
SlovinSCarissimoAPanarielloFGrimaldiABouchéVGambardellaGet al. Single-cell RNA sequencing analysis: A step-by-step overview. Methods Mol Biol. (2021) 2284:343–65. doi: 10.1007/978-1-0716-1307-8_19
35
HwangBLeeJHBangD. Single-cell RNA sequencing technologies and bioinformatics pipelines. Exp Mol Med. (2018) 50:1–14. doi: 10.1038/s12276-018-0071-8
36
ZhengGXTerryJMBelgraderPRyvkinPBentZWWilsonRet al. Massively parallel digital transcriptional profiling of single cells. Nat Commun. (2017) 8:14049. doi: 10.1038/ncomms14049
37
PengXWuJBrunmeirRKimSYZhangQDingCet al. TELP, a sensitive and versatile library construction method for next-generation sequencing. Nucleic Acids Res. (2015) 43:e35. doi: 10.1093/nar/gku818
38
PicelliSBjörklundÅKFaridaniORSagasserSWinbergGSandbergR. Smart-seq2 for sensitive full-length transcriptome profiling in single cells. Nat Methods. (2013) 10:1096–8. doi: 10.1038/nmeth.2639
39
AspMBergenstråhleJLundebergJ. Spatially resolved transcriptomes-next generation tools for tissue exploration. Bioessays. (2020) 42:e1900221. doi: 10.1002/bies.201900221
40
RodriquesSGStickelsRRGoevaAMartinCAMurrayEVanderburgCRet al. Slide-seq: A scalable technology for measuring genome-wide expression at high spatial resolution. Science. (2019) 363:1463–7. doi: 10.1126/science.aaw1219
41
BhattacharyaIDeyS. Emerging concepts on Leydig cell development in fetal and adult testis. Front Endocrinol (Lausanne). (2023) 13:1086276. doi: 10.3389/fendo.2022.1086276
42
ReyRJossoNRacineC. Sexual differentiation. In: FeingoldKRAnawaltBBlackmanMRet al, editors. Endotext. MDText.com, Inc, South Dartmouth (MA (2000). Available at: https://www.ncbi.nlm.nih.gov/books/NBK279001/.
43
RotgersEJørgensenAYaoHH. At the crossroads of fate-somatic cell lineage specification in the fetal gonad. Endocr Rev. (2018) 39:739–59. doi: 10.1210/er.2018-00010
44
HancockGVWamaithaSEPeretzLClarkAT. Mammalian primordial germ cell specification. Development. (2021) 148:dev189217. doi: 10.1242/dev.189217
45
NichollsPKSchorleHNaqviSHuYCFanYCarmellMAet al. Mammalian germ cells are determined after PGC colonization of the nascent gonad. Proc Natl Acad Sci U.S.A. (2019) 116:25677–87. doi: 10.1073/pnas.1910733116
46
SpillerCKoopmanPBowlesJ. Sex determination in the mammalian germline. Annu Rev Genet. (2017) 51:265–85. doi: 10.1146/annurev-genet-120215-035449
47
CapelB. Vertebrate sex determination: evolutionary plasticity of a fundamental switch. Nat Rev Genet. (2017) 18:675–89. doi: 10.1038/nrg.2017.60
48
CoolJDeFalcoTCapelB. Testis formation in the fetal mouse: dynamic and complex de novo tubulogenesis. Wiley Interdiscip Rev Dev Biol. (2012) 1:847–59. doi: 10.1002/wdev.62
49
DeFalcoTBhattacharyaIWilliamsAVSamsDMCapelB. Yolk-sac-derived macrophages regulate fetal testis vascularization and morphogenesis. Proc Natl Acad Sci U.S.A. (2014) 111:E2384–93. doi: 10.1073/pnas.1400057111
50
BhattacharyaISharmaPPurohitSKothiyalSDasMBanerjeeA. Recent update on retinoic acid-driven initiation of spermatogonial differentiation. Front Cell Dev Biol. (2022) 10:833759. doi: 10.3389/fcell.2022.833759
51
BowlesJKoopmanP. Sex determination in mammalian germ cells: extrinsic versus intrinsic factors. Reproduction. (2010) 139:943–58. doi: 10.1530/REP-10-0075
52
SekidoRLovell-BadgeR. Genetic control of testis development. Sex Dev. (2013) 7:21–32. doi: 10.1159/000342221
53
KubotaHBrinsterRL. Spermatogonial stem cells. Biol Reprod. (2018) 99:52–74. doi: 10.1093/biolre/ioy077
54
HermannBPSukhwaniMHanselMCOrwigKE. Spermatogonial stem cells in higher primates: are there differences from those in rodents? Reproduction. (2010) 139:479–93. doi: 10.1530/REP-09-0255
55
de RooijDGReppingSvan PeltAM. Role for adhesion molecules in the spermatogonial stem cell niche. Cell Stem Cell. (2008) 3:467–8. doi: 10.1016/j.stem.2008.10.010
56
de RooijDGGrootegoedJA. Spermatogonial stem cells. Curr Opin Cell Biol. (1998) 10:694–701. doi: 10.1016/s0955-0674(98)80109-9
57
de RooijDG. The nature and dynamics of spermatogonial stem cells. Development. (2017) 144:3022–30. doi: 10.1242/dev.146571
58
BillmyreKK. Chromosome-specific behaviors during early meiosis. Curr Top Dev Biol. (2023) 151:127–54. doi: 10.1016/bs.ctdb.2022.05.002
59
LaroseHKentTMaQShamiANHarerimanaNLiJZet al. Regulation of meiotic progression by Sertoli-cell androgen signaling. Mol Biol Cell. (2020) 31:2841–62. doi: 10.1091/mbc.E20-05-0334
60
HaoSLNiFDYangWX. The dynamics and regulation of chromatin remodeling during spermiogenesis. Gene. (2019) 706:201–10. doi: 10.1016/j.gene.2019.05.027
61
GriswoldMD. 50 years of spermatogenesis: Sertoli cells and their interactions with germ cells. Biol Reprod. (2018) 99:87–100. doi: 10.1093/biolre/ioy027
62
WalkerWHChengJ. FSH and testosterone signaling in Sertoli cells. Reproduction. (2005) 130:15–28. doi: 10.1530/rep.1.00358
63
KarlJCapelB. Sertoli cells of the mouse testis originate from the coelomic epithelium. Dev Biol. (1998) 203:323–33. doi: 10.1006/dbio.1998.9068
64
ZhaoLYaoCXingXJingTLiPZhuZet al. Single-cell analysis of developing and azoospermia human testicles reveals central role of Sertoli cells. Nat Commun. (2020) 11:5683. doi: 10.1038/s41467-020-19414-4
65
BhattacharyaIBasuSPradhanBSSarkarHNagarajanPMajumdarSS. Testosterone augments FSH signaling by upregulating the expression and activity of FSH-Receptor in Pubertal Primate Sertoli cells. Mol Cell Endocrinol. (2019) 482:70–80. doi: 10.1016/j.mce.2018.12.012
66
SharpeRMMcKinnellCKivlinCFisherJS. Proliferation and functional maturation of Sertoli cells, and their relevance to disorders of testis function in adulthood. Reproduction. (2003) 125:769–84. doi: 10.1530/rep.0.1250769
67
BhattacharyaISharmaSSSarkarHGuptaAPradhanBSMajumdarSS. FSH mediated cAMP signalling upregulates the expression of Gα subunits in pubertal rat Sertoli cells. Biochem Biophys Res Commun. (2021) 569:100–5. doi: 10.1016/j.bbrc.2021.06.094
68
BhattacharyaIBasuSSardaKGautamMNagarajanPPradhanBSet al. Low levels of Gαs and Ric8b in testicular sertoli cells may underlie restricted FSH action during infancy in primates. Endocrinology. (2015) 156:1143–55. doi: 10.1210/en.2014-1746
69
BhattacharyaIPradhanBSSardaKGautamMBasuSMajumdarSS. A switch in Sertoli cell responsiveness to FSH may be responsible for robust onset of germ cell differentiation during prepubartal testicular maturation in rats. Am J Physiol Endocrinol Metab. (2012) 303:E886–98. doi: 10.1152/ajpendo.00293.2012
70
MajumdarSSSardaKBhattacharyaIPlantTM. Insufficient androgen and FSH signaling may be responsible for the azoospermia of the infantile primate testes despite exposure to an adult-like hormonal milieu. Hum Reprod. (2012) 27:2515–25. doi: 10.1093/humrep/des184
71
PradhanBSBhattacharyaISarkarRMajumdarSS. Downregulation of sostdc1 in testicular sertoli cells is prerequisite for onset of robust spermatogenesis at puberty. Sci Rep. (2019) 9:11458. doi: 10.1038/s41598-019-47930-x
72
PradhanBSBhattacharyaISarkarRMajumdarSS. Pubertal down-regulation of Tetraspanin 8 in testicular Sertoli cells is crucial for male fertility. Mol Hum Reprod. (2020) 26:760–72. doi: 10.1093/molehr/gaaa055
73
GuptaAVatsAGhosalAMandalKSarkarRBhattacharyaIet al. Follicle-stimulating hormone-mediated decline in miR-92a-3p expression in pubertal mice Sertoli cells is crucial for germ cell differentiation and fertility. Cell Mol Life Sci. (2022) 79:136. doi: 10.1007/s00018-022-04174-9
74
HeinrichABhandaryBPotterSJRatnerNDeFalcoT. Cdc42 activity in Sertoli cells is essential for maintenance of spermatogenesis. Cell Rep. (2021) 37:109885. doi: 10.1016/j.celrep.2021.109885
75
HeinrichAPotterSJGuoLRatnerNDeFalcoT. Distinct roles for rac1 in sertoli cell function during testicular development and spermatogenesis. Cell Rep. (2020) 31:107513. doi: 10.1016/j.celrep.2020.03.077
76
MaMYangSZhangZLiPGongYLiuLet al. Sertoli cells from non-obstructive azoospermia and obstructive azoospermia patients show distinct morphology, Raman spectrum and biochemical phenotype. Hum Reprod. (2013) 28:1863–73. doi: 10.1093/humrep/det068
77
PlottonISanchezPDurandPLejeuneH. Decrease of both stem cell factor and clusterin mRNA levels in testicular biopsies of azoospermic patients with constitutive or idiopathic but not acquired spermatogenic failure. Hum Reprod. (2006) 21:2340–5. doi: 10.1093/humrep/del158
78
StegerKReyRKlieschSLouisFSchleicherGBergmannM. Immunohistochemical detection of immature Sertoli cell markers in testicular tissue of infertile adult men: a preliminary study. Int J Androl. (1996) 19:122–8. doi: 10.1111/j.1365-2605.1996.tb00448.x
79
PaduchDAHilzSGrimsonASchlegelPNJedlickaAEWrightWW. Aberrant gene expression by Sertoli cells in infertile men with Sertoli cell-only syndrome [published correction appears in PLoS One. 2019 Aug 20;14(8):e0221648]. PloS One. (2019) 14:e0216586. doi: 10.1371/journal.pone.0216586
80
PotterSJDeFalcoT. Role of the testis interstitial compartment in spermatogonial stem cell function. Reproduction. (2017) 153:R151–62. doi: 10.1530/REP-16-0588
81
MishraJGautamMDadhichRKowtharapuBSMajumdarSS. Peritubular cells may modulate Leydig cell-mediated testosterone production through a nonclassic pathway. Fertil Steril. (2012) 98:1308–17.e1. doi: 10.1016/j.fertnstert.2012.07.1124
82
ZhangTSunPGengQFanHGongYHuYet al. Disrupted spermatogenesis in a metabolic syndrome model: the role of vitamin A metabolism in the gut-testis axis. Gut. (2022) 71:78–87. doi: 10.1136/gutjnl-2020-323347
83
ZirkinBRPapadopoulosV. Leydig cells: formation, function, and regulation. Biol Reprod. (2018) 99:101–11. doi: 10.1093/biolre/ioy059
84
GuXHeinrichALiSYDeFalcoT. Testicular macrophages are recruited during a narrow fetal time window and promote organ-specific developmental functions. Nat Commun. (2023) 14:1439. doi: 10.1038/s41467-023-37199-0
85
GuanXChenPJiMWenXChenDZhaoXet al. Identification of rat testicular leydig precursor cells by single-cell-RNA-sequence analysis. Front Cell Dev Biol. (2022) 10:805249. doi: 10.3389/fcell.2022.805249
86
ChoubeyMRanjanABoraPSBaltazarFMartinLJKrishnaA. Role of adiponectin as a modulator of testicular function during aging in mice. Biochim Biophys Acta Mol Basis Dis. (2019) 1865:413–27. doi: 10.1016/j.bbadis.2018.11.019
87
ChoubeyMRanjanABoraPSKrishnaA. Protective role of adiponectin against testicular impairment in high-fat diet/streptozotocin-induced type 2 diabetic mice. Biochimie. (2020) 168:41–52. doi: 10.1016/j.biochi.2019.10.014
88
ChoubeyMRanjanAKrishnaA. Adiponectin/AdipoRs signaling as a key player in testicular aging and associated metabolic disorders. Vitam Horm. (2021) 115:611–34. doi: 10.1016/bs.vh.2020.12.024
89
BhushanSMeinhardtA. The macrophages in testis function. J Reprod Immunol. (2017) 119:107–12. doi: 10.1016/j.jri.2016.06.008
90
DeFalcoTPotterSJWilliamsAVWallerBKanMJCapelB. Macrophages contribute to the spermatogonial niche in the adult testis. Cell Rep. (2015) 12:1107–19. doi: 10.1016/j.celrep.2015.07.015
91
GuXLiSYDeFalcoT. Immune and vascular contributions to organogenesis of the testis and ovary. FEBS J. (2022) 289:2386–408. doi: 10.1111/febs.15848
92
SargentKMMcFeeRMSpuri GomesRCuppAS. Vascular endothelial growth factor A: just one of multiple mechanisms for sex-specific vascular development within the testis? J Endocrinol. (2015) 227:R31–50. doi: 10.1530/JOE-15-0342
93
BhangDHKimBJKimBGSchadlerKBaekKHKimYHet al. Testicular endothelial cells are a critical population in the germline stem cell niche. Nat Commun. (2018) 9:4379. doi: 10.1038/s41467-018-06881-z
94
GuoFYanLGuoHLiLHuBZhaoYet al. The transcriptome and DNA methylome landscapes of human primordial germ cells. Cell. (2015) 161:1437–52. doi: 10.1016/j.cell.2015.05.015
95
GuoJGrowEJYiCMlcochovaHMaherGJLindskogCet al. Chromatin and single-cell RNA-seq profiling reveal dynamic signaling and metabolic transitions during human spermatogonial stem cell development. Cell Stem Cell. (2017) 21:533–46. e6. doi: 10.1016/j.stem.2017.09.003
96
NeuhausNYoonJTerwortNKlieschSSeggewissJHugeAet al. Single-cell gene expression analysis reveals diversity among human spermatogonia. Mol Hum Reprod. (2017) 23:79–90. doi: 10.1093/molehr/gaw079
97
WangMLiuXChangGChenYAnGYanLet al. Single-cell RNA-seq analysis maps development of human germline cells and gonadal niche interactions. Cell Stem Cell. (2017) 20:858–873.e4. doi: 10.1016/j.stem.2017.03.007
98
JanSZVormerTLJongejanARölingMDSilberSJde RooijDGet al. Unraveling transcriptome dynamics in human spermatogenesis. Development. (2017) 144:3659–73. doi: 10.1242/dev.152413
99
GuoJGrowEJMlcochovaHMaherGJLindskogCNieXet al. The adult human testis transcriptional cell atlas. Cell Res. (2018) 28:1141–57. doi: 10.1038/s41422-018-0099-2
100
HermannBPChengKSinghARoa-De La CruzLMutojiKNChenICet al. The mammalian spermatogenesis single-cell transcriptome, from spermatogonial stem cells to spermatids. Cell Rep. (2018) 25:1650–1667.e8. doi: 10.1016/j.celrep.2018.10.026
101
WangMLiuXChangGChenYAnGYanLet al. Single-cell RNA sequencing analysis reveals sequential cell fate transition during human spermatogenesis. Cell Stem Cell. (2018) 23:599–614.e4. doi: 10.1016/j.stem.2018.08.007
102
SohniATanKSongHWBurowDde RooijDGLaurentLet al. The neonatal and adult human testis defined at the single-cell level. Cell Rep. (2019) 26:1501–1517.e4. doi: 10.1016/j.celrep.2019.01.045
103
GuoJNieXGieblerMMlcochovaHWangYGrowEJet al. The dynamic transcriptional cell atlas of testis development during human puberty. Cell Stem Cell. (2020) 26:262–276.e4. doi: 10.1016/j.stem.2019.12.005
104
ShamiANZhengXMunyokiSKMaQManskeGLGreenCDet al. Single-cell RNA sequencing of human, macaque, and mouse testes uncovers conserved and divergent features of mammalian spermatogenesis. Dev Cell. (2020) 54:529–547.e12. doi: 10.1016/j.devcel.2020.05.010
105
LauXMunusamyPNgMJSangrithiM. Single-cell RNA sequencing of the cynomolgus macaque testis reveals conserved transcriptional profiles during mammalian spermatogenesis. Dev Cell. (2020) 54:548–566.e7. doi: 10.1016/j.devcel.2020.07.018
106
LiuXChenYTangWZhangLChenWYanZet al. Single-cell transcriptome analysis of the novel coronavirus (SARS-CoV-2) associated gene ACE2 expression in normal and non-obstructive azoospermia (NOA) human male testes. Sci China Life Sci. (2020) 63:1006–15. doi: 10.1007/s11427-020-1705-0
107
TanKSongHWThompsonMMunyokiSSukhwaniMHsiehTCet al. Transcriptome profiling reveals signaling conditions dictating human spermatogonia fate in vitro. Proc Natl Acad Sci U.S.A. (2020) 117:17832–41. doi: 10.1073/pnas.2000362117
108
XiaBYanYBaronMWagnerFBarkleyDChiodinMet al. Widespread transcriptional scanning in the testis modulates gene evolution rates. Cell. (2020) 180:248–262.e21. doi: 10.1016/j.cell.2019.12.015
109
HwangYSSuzukiSSeitaYItoJSakataYAsoHet al. Reconstitution of prospermatogonial specification in vitro from human induced pluripotent stem cells. Nat Commun. (2020) 11:5656. doi: 10.1038/s41467-020-19350-3
110
GuoJSosaEChitiashviliTNieXRojasEJOliverEet al. Single-cell analysis of the developing human testis reveals somatic niche cell specification and fetal germline stem cell establishment. Cell Stem Cell. (2021) 28:764–778.e4. doi: 10.1016/j.stem.2020.12.004
111
ChenHMurrayESinhaALaumasALiJLesmanDet al. Dissecting mammalian spermatogenesis using spatial transcriptomics. Cell Rep. (2021) 37:109915. doi: 10.1016/j.celrep.2021.109915
112
WangRLiuXLiLYangMYongJZhaiFet al. Dissecting human gonadal cell lineage specification and sex determination using A single-cell RNA-seq approach. Genomics Proteomics Bioinf. (2022) 20:223–45. doi: 10.1016/j.gpb.2022.04.002
113
NieXMunyokiSKSukhwaniMSchmidNMisselAEmeryBRet al. Single-cell analysis of human testis aging and correlation with elevated body mass index. Dev Cell. (2022) 57:1160–1176.e5. doi: 10.1016/j.devcel.2022.04.004
114
HuangDZuoYZhangCSunGJingYLeiJet al. A single-nucleus transcriptomic atlas of primate testicular aging reveals exhaustion of the spermatogonial stem cell reservoir and loss of Sertoli cell homeostasis. Protein Cell. (2023) 14:888–907. doi: 10.1093/procel/pwac057
115
MuratFMbengueNWingeSBTrefzerTLeushkinESeppMet al. The molecular evolution of spermatogenesis across mammals. Nature. (2023) 613:308–16. doi: 10.1038/s41586-022-05547-7
116
HuangYLiLAnGYangXCuiMSongXet al. Single-cell multi-omics sequencing of human spermatogenesis reveals a DNA demethylation event associated with male meiotic recombination. Nat Cell Biol. (2023) 25:1520–34. doi: 10.1038/s41556-023-01232-7
117
de CastroPVendrellXEscrichLGrauNGonzalez-MartinRQuiñoneroAet al. Comparative single-cell transcriptomic profiles of human androgenotes and parthenogenotes during early development. Fertil Steril. (2023) 119:675–87. doi: 10.1016/j.fertnstert.2022.12.027
118
SinghAHermannBP. Conserved transcriptome features define prepubertal primate spermatogonial stem cells as adark spermatogonia and identify unique regulators. Int J Mol Sci. (2023) 24:4755. doi: 10.3390/ijms24054755
119
RajachandranSZhangXCaoQCaldeira-BrantALZhangXSongYet al. Dissecting the spermatogonial stem cell niche using spatial transcriptomics. Cell Rep. (2023) 42:112737. doi: 10.1016/j.celrep.2023.112737
120
LukassenSBoschEEkiciABWinterpachtA. Single-cell RNA sequencing of adult mouse testes. Sci Data. (2018) 5:180192. doi: 10.1038/sdata.2018.192
121
LukassenSBoschEEkiciABWinterpachtA. Characterization of germ cell differentiation in the male mouse through single-cell RNA sequencing. Sci Rep. (2018) 8:6521. doi: 10.1038/s41598-018-24725-0
122
GreenCDMaQManskeGLShamiANZhengXMariniSet al. A comprehensive roadmap of murine spermatogenesis defined by single-cell RNA-seq. Dev Cell. (2018) 46:651–667.e10. doi: 10.1016/j.devcel.2018.07.025
123
StévantINeirijnckYBorelCEscoffierJSmithLBAntonarakisSEet al. Deciphering cell lineage specification during male sex determination with single-cell RNA sequencing. Cell Rep. (2018) 22:1589–99. doi: 10.1016/j.celrep.2018.01.043
124
GriveKJHuYShuEGrimsonAElementoOGrenierJKet al. Dynamic transcriptome profiles within spermatogonial and spermatocyte populations during postnatal testis maturation revealed by single-cell sequencing. PloS Genet. (2019) 15:e1007810. doi: 10.1371/journal.pgen.1007810
125
TanKSongHWWilkinsonMF. Single-cell RNAseq analysis of testicular germ and somatic cell development during the perinatal period. Development. (2020) 147:dev183251. doi: 10.1242/dev.183251
126
ZhaoYZhangPGeWFengYLiLSunZet al. Alginate oligosaccharides improve germ cell development and testicular microenvironment to rescue busulfan disrupted spermatogenesis. Theranostics. (2020) 10:3308–24. doi: 10.7150/thno.43189
127
KaluckaJde RooijLPMHGoveiaJRohlenovaKDumasSJMetaEet al. Single-cell transcriptome atlas of murine endothelial cells. Cell. (2020) 180:764–779.e20. doi: 10.1016/j.cell.2020.01.015
128
MayèreCNeirijnckYSararolsPRandsCMStévantIKühneFet al. Single-cell transcriptomics reveals temporal dynamics of critical regulators of germ cell fate during mouse sex determination. FASEB J. (2021) 35:e21452. doi: 10.1096/fj.202002420R
129
YangFWhelanECGuanXDengBWangSSunJet al. FGF9 promotes mouse spermatogonial stem cell proliferation mediated by p38 MAPK signalling. Cell Prolif. (2021) 54:e12933. doi: 10.1111/cpr.12933
130
ZhangWXiaSXiaoWSongYTangLCaoMet al. A single-cell transcriptomic landscape of mouse testicular aging. J Adv Res. (2023) 53:219–34. doi: 10.1016/j.jare.2022.12.007
131
GuanXJiMWenXHuangFZhaoXChenDet al. Single-cell RNA sequencing of adult rat testes after Leydig cell elimination and restoration. Sci Data. (2022) 9:106. doi: 10.1038/s41597-022-01225-5
132
HongSShenXLuoCSunF. Comparative analysis of the testes from wild-type and Alkbh5-knockout mice using single-cell RNA sequencing. G3 (Bethesda). (2022) 12:jkac130. doi: 10.1093/g3journal/jkac130
133
PengYJTangXTShuHSDongWShaoHZhouBO. Sertoli cells are the source of stem cell factor for spermatogenesis. Development. (2023) 150:dev200706. doi: 10.1242/dev.200706
134
SuenHCRaoSLukACSZhangRYangLQiHet al. The single-cell chromatin accessibility landscape in mouse perinatal testis development. Elife. (2023) 12:e75624. doi: 10.7554/eLife.75624
135
KavarthapuRAnbazhaganRPalSDufauML. Single-cell transcriptomic profiling of the mouse testicular germ cells reveals important role of phosphorylated GRTH/DDX25 in round spermatid differentiation and acrosome biogenesis during spermiogenesis. Int J Mol Sci. (2023) 24:3127. doi: 10.3390/ijms24043127
136
JinCWangZLiPTangJJiaoTLiYet al. Decoding the spermatogonial stem cell niche under physiological and recovery conditions in adult mice and humans. Sci Adv. (2023) 9:eabq3173. doi: 10.1126/sciadv.abq3173
137
TanHWangWZhouCWangYZhangSYangPet al. Single-cell RNA-seq uncovers dynamic processes orchestrated by RNA-binding protein DDX43 in chromatin remodeling during spermiogenesis. Nat Commun. (2023) 14:2499. doi: 10.1038/s41467-023-38199-w
138
WaltersBWRainsfordSRHeuerRADiasNHuangXde RooijDet al. KDM6A/UTX promotes spermatogenic gene expression across generations and is not required for male fertility†. Biol Reprod. (2024) 110:391–407. doi: 10.1093/biolre/ioad141
139
LaurentinoSHeckmannLDi PersioSLiXMeyer Zu HörsteG. High-resolution analysis of germ cells from men with sex chromosomal aneuploidies reveals normal transcriptome but impaired imprinting. Clin Epigenet. (2019) 11:127. doi: 10.1186/s13148-019-0720-3
140
WingeSBSoraggiSSchierupMHRajpert-De MeytsEAlmstrupK. Integration and reanalysis of transcriptomics and methylomics data derived from blood and testis tissue of men with 47,XXY Klinefelter syndrome indicates the primary involvement of Sertoli cells in the testicular pathogenesis. Am J Med Genet C Semin Med Genet. (2020) 184:239–55. doi: 10.1002/ajmg.c.31793
141
ChenSAnGWangHWuXPingPHuLet al. Human obstructive (postvasectomy) and nonobstructive azoospermia - Insights from scRNA-Seq and transcriptome analysis. Genes Dis. (2020) 9:766–76. doi: 10.1016/j.gendis.2020.09.004
142
WangMXuYZhangYChenYChangGAnGet al. Deciphering the autophagy regulatory network via single-cell transcriptome analysis reveals a requirement for autophagy homeostasis in spermatogenesis. Theranostics. (2021) 11:5010–27. doi: 10.7150/thno.55645
143
AlfanoMTasciniASPederzoliFLocatelliINebuloniMGianneseFet al. Aging, inflammation and DNA damage in the somatic testicular niche with idiopathic germ cell aplasia. Nat Commun. (2021) 12:5205. doi: 10.1038/s41467-021-25544-0
144
Di PersioSTekathTSiebert-KussLMCremersJFWistubaJLiXet al. Single-cell RNA-seq unravels alterations of the human spermatogonial stem cell compartment in patients with impaired spermatogenesis. Cell Rep Med. (2021) 2:100395. doi: 10.1016/j.xcrm.2021.100395
145
MahyariEGuoJLimaACLewinsohnDPStendahlAMVigh-ConradKAet al. Comparative single-cell analysis of biopsies clarifies pathogenic mechanisms in Klinefelter syndrome. Am J Hum Genet. (2021) 108:1924–45. doi: 10.1016/j.ajhg.2021.09.001
146
TangXJXiaoQHWangXLHeYTianYNXiaBTet al. Single-cell transcriptomics-based study of transcriptional regulatory features in the non-obstructive azoospermia testis. Front Genet. (2022) 13:875762. doi: 10.3389/fgene.2022.875762
147
Di PersioSNeuhausN. Human spermatogonial stem cells and their niche in male (in)fertility: novel concepts from single-cell RNA-sequencing. Hum Reprod. (2023) 38:1–13. doi: 10.1093/humrep/deac245
Summary
Keywords
single-cell RNA-sequencing (Sc-RNA-seq), testis, spatial transcriptomics, spermatogenesis, male infertility
Citation
Tirumalasetty MB, Bhattacharya I, Mohiuddin MS, Baki VB and Choubey M (2024) Understanding testicular single cell transcriptional atlas: from developmental complications to male infertility. Front. Endocrinol. 15:1394812. doi: 10.3389/fendo.2024.1394812
Received
02 March 2024
Accepted
14 June 2024
Published
11 July 2024
Volume
15 - 2024
Edited by
Maria Eugenia Teves, Virginia Commonwealth University, United States
Reviewed by
Carlos Córdova-Fletes, Universidad Autonoma de Nuevo León, Mexico
Shivani Srivastava, Yale University, United States
Updates
Copyright
© 2024 Tirumalasetty, Bhattacharya, Mohiuddin, Baki and Choubey.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mayank Choubey, choubeymayank48@gmail.com; Mayank.Choubey@nyulangone.org
†These authors have contributed equally to this work
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.