Single-Cell RNA Sequencing of Human, Macaque, and Mouse Testes Uncovers Conserved and Divergent Features of Mammalian Spermatogenesis

An organism-wide atlas of hormonal signaling based on the mouse lemur single-cell transcriptome

Hormones mediate long-range cell communication and play vital roles in physiology, metabolism, and health. Traditionally, endocrinologists have focused on one hormone or organ system at a time. Yet, hormone signaling by its very nature connects cells of different organs and involves crosstalk of different hormones. Here, we leverage the organism-wide single cell transcriptional atlas of a non-human primate, the mouse lemur (Microcebus murinus), to systematically map source and target cells for 84 classes of hormones. This work uncovers previously-uncharacterized sites of hormone regulation, and shows that the hormonal signaling network is densely connected, decentralized, and rich in feedback loops. Evolutionary comparisons of hormonal genes and their expression patterns show that mouse lemur better models human hormonal signaling than mouse, at both the genomic and transcriptomic levels, and reveal primate-specific rewiring of hormone-producing/target cells. This work complements the scale and resolution of classical endocrine studies and sheds light on primate hormone regulation.

Construction and analysis of a joint diagnostic model of machine learning for cryptorchidism based on single-cell sequencing

BACKGROUND

Cryptorchidism is a condition in which one or both of a baby's testicles do not fully descend into the bottom of the scrotum. Newborns with cryptorchidism are at increased risk of developing infertility later in life. The aim of this study was to develop a novel diagnostic model for cryptorchidism and to identify new biomarkers associated with cryptorchidism.

METHODS

The study data were obtained from RNA sequencing data of cryptorchid patients from Nantong University Hospital and the Gene Expression Omnibus (GEO) database. Differential expression analysis was used to obtain differentially expressed genes (DEGs) between the control and cryptorchid groups. These DEGs were analyzed for their functions by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment using GSEA software. Random Forest algorithm was used to screen central genes based on these DEGs. Neuralnet software package was used to develop artificial neural network models. Based on clinical data, receiver operating characteristic (ROC) was used to validate the models. Single-cell sequencing analysis was used for the pathogenesis of cryptorchidism.

RESULTS

We obtained a total of 525 important DEGs related to cryptorchidism, which are mainly associated with biological functions such as supramolecular complexes and microtubule cytoskeleton. Random forest approach screening obtained eight hub genes. A neural network based on the hub genes showed a 100% success rate of the model. Finally, single-cell sequencing analysis validated the hub genes.

CONCLUSION

We developed a novel diagnostic model for cryptorchidism using artificial neural networks and validated its utility as an effective diagnostic tool.

Characterization of the human fetal gonad and reproductive tract by single-cell transcriptomics

During human fetal development, sex differentiation occurs not only in the gonads but also in the adjacent developing reproductive tract. However, while the cellular composition of male and female human fetal gonads is well described, that of the adjacent developing reproductive tract remains poorly characterized. Here, we performed single-cell transcriptomics on male and female human fetal gonads together with the adjacent developing reproductive tract from first and second trimesters, highlighting the morphological and molecular changes during sex differentiation. We validated different cell populations of the developing reproductive tract and gonads and compared the molecular signatures between the first and second trimesters, as well as between sexes, to identify conserved and sex-specific features. Together, our study provides insights into human fetal sex-specific gonadogenesis and development of the reproductive tract beyond the gonads.

The dynamic genetic determinants of increased transcriptional divergence in spermatids

Cis-genetic effects are key determinants of transcriptional divergence in discrete tissues and cell types. However, how cis- and trans-effects act across continuous trajectories of cellular differentiation in vivo is poorly understood. Here, we quantify allele-specific expression during spermatogenic differentiation at single-cell resolution in an F1 hybrid mouse system, allowing for the comprehensive characterisation of cis- and trans-genetic effects, including their dynamics across cellular differentiation. Collectively, almost half of the genes subject to genetic regulation show evidence for dynamic cis-effects that vary during differentiation. Our system also allows us to robustly identify dynamic trans-effects, which are less pervasive than cis-effects. In aggregate, genetic effects were strongest in round spermatids, which parallels their increased transcriptional divergence we identified between species. Our approach provides a comprehensive quantification of the variability of genetic effects in vivo, and demonstrates a widely applicable strategy to dissect the impact of regulatory variants on gene regulation in dynamic systems.

Single-cell transcriptome analyses reveal critical regulators of spermatogonial stem cell fate transitions

BACKGROUND

Spermatogonial stem cells (SSCs) are the foundation cells for continual spermatogenesis and germline regeneration in mammals. SSC activities reside in the undifferentiated spermatogonial population, and currently, the molecular identities of SSCs and their committed progenitors remain unclear.

RESULTS

We performed single-cell transcriptome analysis on isolated undifferentiated spermatogonia from mice to decipher the molecular signatures of SSC fate transitions. Through comprehensive analysis, we delineated the developmental trajectory and identified candidate transcription factors (TFs) involved in the fate transitions of SSCs and their progenitors in distinct states. Specifically, we characterized the Asingle spermatogonial subtype marked by the expression of Eomes. Eomes+ cells contained enriched transplantable SSCs, and more than 90% of the cells remained in the quiescent state. Conditional deletion of Eomes in the germline did not impact steady-state spermatogenesis but enhanced SSC regeneration. Forced expression of Eomes in spermatogenic cells disrupted spermatogenesis mainly by affecting the cell cycle progression of undifferentiated spermatogonia. After injury, Eomes+ cells re-enter the cell cycle and divide to expand the SSC pool. Eomes+ cells consisted of 7 different subsets of cells at single-cell resolution, and genes enriched in glycolysis/gluconeogenesis and the PI3/Akt signaling pathway participated in the SSC regeneration process.

CONCLUSIONS

In this study, we explored the molecular characteristics and critical regulators of subpopulations of undifferentiated spermatogonia. The findings of the present study described a quiescent SSC subpopulation, Eomes+ spermatogonia, and provided a dynamic transcriptional map of SSC fate determination.

scKWARN: Kernel-weighted-average robust normalization for single-cell RNA-seq data

MOTIVATION

Single-cell RNA-seq normalization is an essential step to correct unwanted biases caused by sequencing depth, capture efficiency, dropout, and other technical factors. Existing normalization methods primarily reduce biases arising from sequencing depth by modeling count-depth relationship and/or assuming a specific distribution for read counts. However, these methods may lead to over or under-correction due to presence of technical biases beyond sequencing depth and the restrictive assumption on models and distributions.

RESULTS

We present scKWARN, a Kernel Weighted Average Robust Normalization designed to correct known or hidden technical confounders without assuming specific data distributions or count-depth relationships. scKWARN generates a pseudo expression profile for each cell by borrowing information from its fuzzy technical neighbors through a kernel smoother. It then compares this profile against the reference derived from cells with the same bimodality patterns to determine the normalization factor. As demonstrated in both simulated and real datasets, scKWARN outperforms existing methods in removing a variety of technical biases while preserving true biological heterogeneity.

AVAILABILITY AND IMPLEMENTATION

scKWARN is freely available at https://github.com/cyhsuTN/scKWARN.

Favorable culture conditions for spermatogonial propagation in human and non-human primate primary testicular cell cultures: a systematic review and meta-analysis

Introduction: Autologous transplantation of spermatogonial stem cells (SSCs) isolated from cryopreserved testicular biopsies obtained before oncological treatment could restore fertility in male childhood cancer survivors. There is a clear necessity for in vitro propagation of the limited SSCs from the testicular biopsy prior to transplantation due to limited numbers of spermatogonia in a cryopreserved testicular biopsy. Still, there is no consensus regarding their optimal culture method. Methods: We performed a systematic review and meta-analysis of studies reporting primary testicular cell cultures of human and non-human primate origin through use of Pubmed, EMBASE, and Web of Science core collection databases. Of 760 records, we included 42 articles for qualitative and quantitative analysis. To quantify in vitro spermatogonial propagation, spermatogonial colony doubling time (CDT) was calculated, which measures the increase in the number of spermatogonial colonies over time. A generalized linear mixed model analysis was used to assess the statistical effect of various culture conditions on CDT. Results: Our analysis indicates decreased CDTs, indicating faster spermatogonial propagation in cultures with a low culture temperature (32°C); with use of non-cellular matrices; use of StemPro-34 medium instead of DMEM; use of Knockout Serum Replacement; and when omitting additional growth factors in the culture medium. Discussion: The use of various methods and markers to detect the presence of spermatogonia within the reported cultures could result in detection bias, thereby potentially influencing comparability between studies. However, through use of CDT in the quantitative analysis this bias was reduced. Our results provide insight into critical culture conditions to further optimize human spermatogonial propagation in vitro, and effectively propagate and utilize these cells in a future fertility restoration therapy and restore hope of biological fatherhood for childhood cancer survivors.

TAF7L regulates early stages of male germ cell development in the rat

Male germ cell development is dependent on the orchestrated regulation of gene networks. TATA-box binding protein associated factors (TAFs) facilitate interactions of TATA-binding protein with the TATA element, which is known to coordinate gene transcription during organogenesis. TAF7 like (Taf7l) is situated on the X chromosome and has been implicated in testis development. We examined the biology of TAF7L in testis development using the rat. Taf7l was prominently expressed in preleptotene to leptotene spermatocytes. To study the impact of TAF7L on the testis we generated a global loss-of-function rat model using CRISPR/Cas9 genome editing. Exon 3 of the Taf7l gene was targeted. A founder was generated possessing a 110 bp deletion within the Taf7l locus, which resulted in a frameshift and the premature appearance of a stop codon. The mutation was effectively transmitted through the germline. Deficits in TAF7L did not adversely affect pregnancy or postnatal survival. However, the Taf7l disruption resulted in male infertility due to compromised testis development and failed sperm production. Mutant germ cells suffer meiotic arrest at late zygotene/early pachynema stages, with defects in sex body formation. This testis phenotype was more pronounced than previously described for the subfertile Taf7l null mouse. We conclude that TAF7L is essential for male germ cell development in the rat.

scRNA-seq reveals that origin recognition complex subunit 6 regulates mouse spermatogonial cell proliferation and apoptosis via activation of Wnt/β-catenin signaling

The regulation of spermatogonial proliferation and apoptosis is of great significance for maintaining spermatogenesis. The single-cell RNA sequencing (scRNA-seq) analysis of the testis was performed to identify genes upregulated in spermatogonia. Using scRNA-seq analysis, we identified the spermatogonia upregulated gene origin recognition complex subunit 6 ( Orc6 ), which is involved in DNA replication and cell cycle regulation; its protein expression in the human and mouse testis was detected by western blot and immunofluorescence. To explore the potential function of Orc6 in spermatogonia, the C18-4 cell line was transfected with control or Orc6 siRNA. Subsequently, 5-ethynyl-2-deoxyuridine (EdU) and terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL) assays, flow cytometry, and western blot were used to evaluate its effects on proliferation and apoptosis. It was revealed that ORC6 could promote proliferation and inhibit apoptosis of C18-4 cells. Bulk RNA sequencing and bioinformatics analysis indicated that Orc6 was involved in the activation of wingless/integrated (Wnt)/β-catenin signaling. Western blot revealed that the expression of β-catenin protein and its phosphorylation (Ser675) were significantly decreased when silencing the expression of ORC6. Our findings indicated that Orc6 was upregulated in spermatogonia, whereby it regulated proliferation and apoptosis by activating Wnt/β-catenin signaling.

Integrated single cell transcriptome sequencing analysis reveals species-specific genes and molecular pathways for pig spermiogenesis

Mammalian spermatogenesis is a highly complicated and intricately organized process involving spermatogonia propagation (mitosis) and meiotic differentiation into mature sperm cells (spermiogenesis). In pigs, spermatogonia development and the role of somatic cells in spermatogenesis were previously investigated in detail. However, the characterization of key molecules fundamental to pig spermiogenesis remains less explored. Here we compared spermatogenesis between humans and pigs, focusing on spermiogenesis, by integrative testicular single-cell RNA sequencing (scRNA-seq) analysis. Human and pig testicular cells were clustered into 26 different groups, with cell-type-specific markers and signalling pathways. For spermiogenesis, pseudo-time analysis classified the lineage differentiation routes for round, elongated spermatids and spermatozoa. Moreover, markers and molecular pathways specific to each type of spermatids were examined for humans and pigs, respectively. Furthermore, high-dimensional weighted gene co-expression network analysis (hdWGCNA) identified gene modules specific for each type of human and pig spermatids. Hub genes (pig: SNRPD2.1 related to alternative splicing; human: CATSPERZ, Ca[2+] ion channel) potentially involved in spermiogenesis were also revealed. Taken together, our integrative analysis found that human and pig spermiogeneses involve specific genes and molecular pathways and provided resources and insights for further functional investigation on spermatid maturation and male reproductive ability.

A postmeiotically bifurcated roadmap of honeybee spermatogenesis marked by phylogenetically restricted genes

Haploid males of hymenopteran species produce gametes through an abortive meiosis I followed by meiosis II that can either be symmetric or asymmetric in different species. Thus, one spermatocyte could give rise to two spermatids with either equal or unequal amounts of cytoplasm. It is currently unknown what molecular features accompany these postmeiotic sperm cells especially in species with asymmetric meiosis II such as bees. Here we present testis single-cell RNA sequencing datasets from the honeybee (Apis mellifera) drones of 3 and 14 days after emergence (3d and 14d). We show that, while 3d testes exhibit active, ongoing spermatogenesis, 14d testes only have late-stage spermatids. We identify a postmeiotic bifurcation in the transcriptional roadmap during spermatogenesis, with cells progressing toward the annotated spermatids (SPT) and small spermatids (sSPT), respectively. Despite an overall similarity in their transcriptomic profiles, sSPTs express the fewest genes and the least RNA content among all the sperm cell types. Intriguingly, sSPTs exhibit a relatively high expression level for Hymenoptera-restricted genes and a high mutation load, suggesting that the special meiosis II during spermatogenesis in the honeybee is accompanied by phylogenetically young gene activities.

Characterization of long-term ex vivo expansion of tree shrew spermatogonial stem cells

Tree shrews ( Tupaia belangeri chinensis) share a close relationship to primates and have been widely used in biomedical research. We previously established a spermatogonial stem cell (SSC)-based gene editing platform to generate transgenic tree shrews. However, the influences of long-term expansion on tree shrew SSC spermatogenesis potential remain unclear. Here, we examined the in vivo spermatogenesis potential of tree shrew SSCs cultured across different passages. We found that SSCs lost spermatogenesis ability after long-term expansion (>50 passages), as indicated by the failure to colonize the seminiferous epithelium and generate donor spermatogonia (SPG)-derived spermatocytes or spermatids marking spermatogenesis. RNA sequencing (RNA-seq) analysis of undifferentiated SPGs across different passages revealed significant gene expression changes after sub-culturing primary SPG lines for more than 40 passages on feeder layers. Specifically, DNA damage response and repair genes (e.g., MRE11, SMC3, BLM, and GEN1) were down-regulated, whereas genes associated with mitochondrial function (e.g., NDUFA9, NDUFA8, NDUFA13, and NDUFB8) were up-regulated after expansion. The DNA damage accumulation and mitochondrial dysfunction were experimentally validated in high-passage cells. Supplementation with nicotinamide adenine dinucleotide (NAD +) precursor nicotinamide riboside (NR) exhibited beneficial effects by reducing DNA damage accumulation and mitochondrial dysfunction in SPG elicited by long-term culture. Our research presents a comprehensive analysis of the genetic and physiological attributes critical for the sustained expansion of undifferentiated SSCs in tree shrews and proposes an effective strategy for extended in vitro maintenance.

A single-cell transcriptomic landscape of mouse testicular aging

INTRODUCTION

Advanced paternal age of reproduction is an increasing trend, especially in developed countries and areas. This trend results in elevated risks of adverse reproductive outcomes such as reduced fertility rates, increased pregnancy loss, and poor childhood health. Yet, a systematic profiling of aging-associated molecular and cellular alterations in testicular tissue is still missing.

OBJECTIVES

We aimed to dissect aging-associated molecular characteristics in testes of mice.

METHODS

Single-cell transcriptomic sequencing and analysis were conducted in testes of young (2 months old) and old mice (24 months old). Immunofluorescences and immunochemistry were used to characterize aging-associated phenotypes and verify single cell sequence results.

RESULTS

Here, we constructed the first single-cell transcriptomic atlases of testes of young and old mice. In-depth dissection of aging-dependent transcriptional alterations in specific cell types revealed multiple dysregulated biological processes such as increased 'senescence-associated secretory phenotype' and 'inflammation', which were major aging-associated characteristics. Further analysis of aging-related differentially expressed genes uncovered a disrupted balance of undifferentiated and differentiated spermatogonia stem cells in spermatogonia, indicative of a potential role of spermatogonia stem cells in aging-associated subfertility. Importantly, for the first time, our results identified an increased subtype of aging-specific macrophages, which may contribute to a hostile proinflammatory microenvironment during testicular aging.

CONCLUSION

Taken together, our findings depict the distinct single-cell transcriptional features of the aged mouse testes and provide enormous resources for a comprehensive understanding of the cell-type-specific molecular mechanisms underlying mouse testicular aging, which may shed light on developing novel potential diagnostic biomarkers and therapeutic targets for age-associated male subfertility.

Transcription and chromatin regulation by TAF4b during cellular quiescence of developing prospermatogonia

Prospermatogonia (ProSpg) link the embryonic development of male primordial germ cells to the healthy establishment of postnatal spermatogonia and spermatogonial stem cells. While these spermatogenic precursor cells undergo the characteristic transitions of cycling and quiescence, the transcriptional events underlying these developmental hallmarks remain unknown. Here, we investigated the expression and function of TBP-associated factor 4b (Taf4b) in the timely development of quiescent mouse ProSpg using an integration of gene expression profiling and chromatin mapping. We find that Taf4b mRNA expression is elevated during the transition of mitotic-to-quiescent ProSpg and Taf4b-deficient ProSpg are delayed in their entry into quiescence. Gene ontology, protein network analysis, and chromatin mapping demonstrate that TAF4b is a direct and indirect regulator of chromatin and cell cycle-related gene expression programs during ProSpg quiescence. Further validation of these cell cycle mRNA changes due to the loss of TAF4b was accomplished via immunostaining for proliferating cell nuclear antigen (PCNA). Together, these data indicate that TAF4b is a key transcriptional regulator of the chromatin and quiescent state of the developing mammalian spermatogenic precursor lineage.

TAF7L REGULATES EARLY STAGES OF MALE GERM CELL DEVELOPMENT

Male germ cell development is dependent on the orchestrated regulation of gene networks. TATA-box binding protein associated factors (TAFs) facilitate interactions of TATA-binding protein with the TATA element, which is known to coordinate gene transcription during organogenesis. TAF7 like (Taf7l) is situated on the X chromosome and has been implicated in testis development. We examined the biology of TAF7L in testis development using the rat. Taf7l was prominently expressed in preleptotene to leptotene spermatocytes. To study the impact of TAF7L on the testis we generated a global loss-of-function rat model using CRISPR/Cas9 genome editing. Exon 3 of the Taf7l gene was targeted. A founder was generated possessing a 110 bp deletion within the Taf7l locus, which resulted in a frameshift and the premature appearance of a stop codon. The mutation was effectively transmitted through the germline. Deficits in TAF7L did not adversely affect pregnancy or postnatal survival. However, the Taf7l disruption resulted in male infertility due to compromised testis development and failed sperm production. Mutant germ cells suffer meiotic arrest at the zygotene stage, with defects in sex body formation and meiotic sex chromosome inactivation. This testis phenotype was more pronounced than previously described for the subfertile Taf7l null mouse. We conclude that TAF7L is essential for male germ cell development in the rat.

Integrated analysis of the ubiquitination mechanism reveals the specific signatures of tissue and cancer

BACKGROUND

Ubiquitination controls almost all cellular processes. The dysregulation of ubiquitination signals is closely associated with the initiation and progression of multiple diseases. However, there is little comprehensive research on the interaction and potential function of ubiquitination regulators (UBRs) in spermatogenesis and cancer.

METHODS

We systematically characterized the mRNA and protein expression of UBRs across tissues and further evaluated their roles in testicular development and spermatogenesis. Subsequently, we explored the genetic alterations, expression perturbations, cancer hallmark-related pathways, and clinical relevance of UBRs in pan-cancer.

RESULTS

This work reveals heterogeneity in the expression patterns of UBRs across tissues, and the expression pattern in testis is the most distinct. UBRs are dynamically expressed during testis development, which are critical for normal spermatogenesis. Furthermore, UBRs have widespread genetic alterations and expression perturbations in pan-cancer. The expression of 79 UBRs was identified to be closely correlated with the activity of 32 cancer hallmark-related pathways, and ten hub genes were screened for further clinical relevance analysis by a network-based method. More than 90% of UBRs can affect the survival of cancer patients, and hub genes have an excellent prognostic classification for specific cancer types.

CONCLUSIONS

Our study provides a comprehensive analysis of UBRs in spermatogenesis and pan-cancer, which can build a foundation for understanding male infertility and developing cancer drugs in the aspect of ubiquitination.

The molecular evolution of mammalian spermatogenesis

The testis is a key male reproductive organ that produces gametes through the process of spermatogenesis. Testis morphologies, sperm phenotypes, and the process of spermatogenesis evolve rapidly in mammals, presumably due to the evolutionary pressure on males to give rise to their own offspring. Here, we review studies illuminating the molecular evolution of the testis, in particular large-scale transcriptomic studies, which were based on bulk tissue samples and, more recently, individual cells. Together with various genomic and epigenomic data, these studies have unveiled the cellular source, molecular mechanisms, and evolutionary forces that underlie the rapid phenotypic evolution of the testis. They also revealed shared (ancestral) and species-specific spermatogenic gene expression programs. The insights and available data that have accumulated also provide a valuable resource for the investigation and treatment of male fertility disorders - a dramatically increasing problem in modern industrial societies.

A single-cell landscape of triptolide-associated testicular toxicity in mice

Triptolide is a key active component of the widely used traditional Chinese herb medicine Tripterygium wilfordii Hook. F. Although triptolide exerts multiple biological activities and shows promising efficacy in treating inflammatory-related diseases, its well-known safety issues, especially reproductive toxicity has aroused concerns. However, a comprehensive dissection of triptolide-associated testicular toxicity at single cell resolution is still lacking. Here, we observed testicular toxicity after 14 days of triptolide exposure, and then constructed a single-cell transcriptome map of 59,127 cells in mouse testes upon triptolide-treatment. We identified triptolide-associated shared and cell-type specific differentially expressed genes, enriched pathways, and ligand-receptor pairs in different cell types of mouse testes. In addition to the loss of germ cells, our results revealed increased macrophages and the inflammatory response in triptolide-treated mouse testes, suggesting a critical role of inflammation in triptolide-induced testicular injury. We also found increased reactive oxygen species (ROS) signaling and downregulated pathways associated with spermatid development in somatic cells, especially Leydig and Sertoli cells, in triptolide-treated mice, indicating that dysregulation of these signaling pathways may contribute to triptolide-induced testicular toxicity. Overall, our high-resolution single-cell landscape offers comprehensive information regarding triptolide-associated gene expression profiles in major cell types of mouse testes at single cell resolution, providing an invaluable resource for understanding the underlying mechanism of triptolide-associated testicular injury and additional discoveries of therapeutic targets of triptolide-induced male reproductive toxicity.

Dissecting the spermatogonial stem cell niche using spatial transcriptomics

Spermatogonial stem cells (SSCs) in the testis support the lifelong production of sperm. SSCs reside within specialized microenvironments called "niches," which are essential for SSC self-renewal and differentiation. However, our understanding of the molecular and cellular interactions between SSCs and niches remains incomplete. Here, we combine spatial transcriptomics, computational analyses, and functional assays to systematically dissect the molecular, cellular, and spatial composition of SSC niches. This allows us to spatially map the ligand-receptor (LR) interaction landscape in both mouse and human testes. Our data demonstrate that pleiotrophin regulates mouse SSC functions through syndecan receptors. We also identify ephrin-A1 as a potential niche factor that influences human SSC functions. Furthermore, we show that the spatial re-distribution of inflammation-related LR interactions underlies diabetes-induced testicular injury. Together, our study demonstrates a systems approach to dissect the complex organization of the stem cell microenvironment in health and disease.

A pathogenic variant in the uncharacterized RNF212B gene results in severe aneuploidy male infertility and repeated IVF failure

Quantitative and qualitative spermatogenic impairments are major causes of men's infertility. Although in vitro fertilization (IVF) is effective, some couples persistently fail to conceive. To identify causal variants in patients with severe male infertility factor and repeated IVF failures, we sequenced the exome of two consanguineous family members who underwent several failed IVF cycles and were diagnosed with low sperm count and motility. We identified a rare homozygous nonsense mutation in a previously uncharacterized gene, RNF212B, as the causative variant. Recurrence was identified in another unrelated, infertile patient who also faced repeated failed IVF treatments. scRNA-seq demonstrated meiosis-specific expression of RNF212B. Sequence analysis located a protein domain known to be associated with aneuploidy, which can explain multiple IVF failures. Accordingly, FISH analysis revealed a high aneuploidy rate in the patients' sperm cells and their IVF embryos. Finally, inactivation of the Drosophila orthologs significantly reduced male fertility. Given that members of the evolutionary conserved RNF212 gene family are involved in meiotic recombination and crossover maturation, our findings indicate a critical role of RNF212B in meiosis, genome stability, and in human fertility. Since recombination is completely absent in Drosophila males, our findings may indicate an additional unrelated role for the RNF212-like paralogs in spermatogenesis.

Identification of follicle-stimulating hormone-responsive genes in Sertoli cells during early postnatal mouse testis development

BACKGROUND

In the mouse testis, Sertoli cells rapidly divide during a narrow window of time pre-pubertally and differentiate thereafter. The number of Sertoli cells determines the testis size and germ cell-carrying capacity. Follicle-stimulating hormone (FSH) binds its cognate FSH-receptors expressed on Sertoli cells and acts as a mitogen to regulate their proliferation. Fshb-/- mutant adult male mice have reduced Sertoli cell number and testis size and reduced sperm number and motility. However, FSH-responsive genes in early postnatal mouse Sertoli cells are unknown.

OBJECTIVES

To identify FSH-responsive genes in early postnatal mouse Sertoli cells.

MATERIALS AND METHODS

A fluorescence-activated cell sorting method was developed to rapidly purify Sertoli cells from control and Fshb-/- mice carrying a Sox9 GfpKI allele. These pure Sertoli cells were used for large-scale gene expression analyses.

RESULTS

We show that mouse Sertoli cells rarely divide beyond postnatal day 7. Our in vivo BrdU labeling studies indicate loss of FSH results in a 30% reduction in Sertoli cell proliferation in mice at 5 days of age. Flowsorted GFP+ Sertoli cells with maximal Fshr expression were 97%-98% pure and mostly devoid of Leydig and germ cells as assessed by Taqman qPCR quantification of gene expression and immunolabeling of the corresponding cell-specific markers. Large-scale gene expression analysis identified several differentially regulated genes in flow-sorted GFP+ Sertoli cells obtained from testis of control and Fshb-/- mice at 5 days of age. The top 25 networks identified by pathway analysis include those related to the cell cycle, cell survival and most importantly, carbohydrate and lipid metabolism and molecular transport.

DISCUSSION

Several of the FSH-responsive genes identified in this study could serve as useful markers for Sertoli cell proliferation in normal physiology, toxicant-induced Sertoli cell/testis injury, and other pathological conditions.

CONCLUSION

Our studies reveal that FSH-regulates macromolecular metabolism and molecular transport networks of genes in early postnatal Sertoli cells most likely in preparation for establishment of functional associations with germ cells to successfully coordinate spermatogenesis.

Loss Of Chromodomain of Male-Specific Lethal 3 (MSL3) Does Not Affect Spermatogenesis In Rodents

Msl3 is a member of the chromatin-associated male-specific lethal MSL complex which is responsible for the transcriptional upregulation of genes on the X chromosome in males Drosophila. Although the dosage complex operates differently in mammals, the Msl3 gene is conserved from flies to humans. Msl3 is required for meiotic entry during Drosophila oogenesis. Recent reports indicate that also in primates, Msl3 is expressed in undifferentiated germline cells before meiotic entry. However, if Msl3 plays a role in the meiotic entry of mammals has yet to be explored. To study this, we used mouse spermatogenesis as a study model. Analyses of single cells RNA-seq data revealed that, in mice, Msl3 is mostly expressed in meiotic cells. To test the role of Msl3 in meiosis, we used a male germline-specific Stra8-iCre driver and a newly generated Msl3flox conditional knock-out mouse line. Msl3 conditional loss-of-function in spermatogonia did not cause spermatogenesis defects or changes in the expression of genes related to meiosis. Our data suggest that, in mice, Msl3 exhibits delayed expression compared to Drosophila and primates, and loss-of-function mutations disrupting the chromodomain of Msl3 alone do not impede meiotic entry in rodents.

Analysis of Notch1 signaling in mammalian sperm development

OBJECTIVE

A mammalian Delta-Notch signaling component, Notch1, has been suggested for its expression during the normal sperm development although its conditional deletion caused no apparent abnormalities. Since we established our original transgenic mouse system that enabled labeling of past and ongoing Notch1 signaling at a cellular level, we tried to validate that observation in vivo. Our transgenic mouse system used Cre/loxP system to induce tandem dsRed expression upon Notch1 signaling.

RESULTS

To our surprise, we were unable to observe tandem dsRed expression in the seminiferous tubules where the sperms developed. In addition, tandem dsRed expression was lacking in the somatic cells of the next generation in our transgenic mouse system, suggesting that sperms received no Notch1 signaling during their development. To validate this result, we conducted re-analysis of four single-cell RNA-seq datasets from mouse and human testes and showed that Notch1 expression was little in the sperm cell lineage. Collectively, our results posed a question into the involvement of Notch1 in the normal sperm development although this observation may help the interpretation of the previous result that Notch1 conditional deletion caused no apparent abnormalities in murine spermatogenesis.

Single-cell transcriptomic profiling of microvascular endothelial cell heterogeneity in congenital diaphragmatic hernia

Congenital diaphragmatic hernia (CDH) is a neonatal anomaly that includes pulmonary hypoplasia and hypertension. We hypothesized that microvascular endothelial cell (EC) heterogeneity is different in CDH lungs and related to lung underdevelopment and remodeling. To test this, we evaluated rat fetuses at E21.5 in a nitrofen model of CDH to compare lung transcriptomes among healthy controls (2HC), nitrofen-exposed controls (NC) and nitrofen-exposed subjects with CDH. Single-cell RNA sequencing with unbiased clustering revealed 3 distinct microvascular EC clusters: a general population (mvEC), a proliferative population and a population high in hemoglobin. Only the CDH mvEC cluster had a distinct inflammatory transcriptomic signature as compared to the 2HC and NC endothelial cells, e.g. greater activation and adhesion of inflammatory cells and production of reactive oxygen species. Furthermore, CDH mvECs had downregulated Ca4, Apln and Ednrb gene expression. Those genes are markers for ECs important to lung development, gas exchange and alveolar repair (mvCa4+). mvCa4+ ECs were reduced in CDH (2HC [22.6%], NC [13.1%] and CDH [5.3%], p < 0.0001). Overall, these findings identify transcriptionally distinct microvascular endothelial cell clusters in CDH, including the distinctly inflammatory mvEC cluster and the depleted group of mvCa4+ ECs, which together may contribute to pathogenesis.

Developmental origins of mammalian spermatogonial stem cells: New perspectives on epigenetic regulation and sex chromosome function

Male and female germ cells undergo genome-wide reprogramming during their development, and execute sex-specific programs to complete meiosis and successfully generate healthy gametes. While sexually dimorphic germ cell development is fundamental, similarities and differences exist in the basic processes governing normal gametogenesis. At the simplest level, male gamete generation in mammals is centred on the activity of spermatogonial stem cells (SSCs), and an equivalent cell state is not present in females. Maintaining this unique SSC epigenetic state, while keeping to germ cell-intrinsic developmental programs, poses challenges for the correct completion of spermatogenesis. In this review, we highlight the origins of spermatogonia, comparing and contrasting them with female germline development to emphasize specific developmental processes that are required for their function as germline stem cells. We identify gaps in our current knowledge about human SSCs and further discuss the impact of the unique regulation of the sex chromosomes during spermatogenesis, and the roles of X-linked genes in SSCs.

Genomic study of TEX15 variants: prevalence and allelic heterogeneity in men with spermatogenic failure

Introduction: Human spermatogenesis is a highly intricate process that requires the input of thousands of testis-specific genes. Defects in any of them at any stage of the process can have detrimental effects on sperm production and/or viability. In particular, the function of many meiotic proteins encoded by germ cell specific genes is critical for maturation of haploid spermatids and viable spermatozoa, necessary for fertilization, and is also extremely sensitive to even the slightest change in coding DNA. Methods: Here, using whole exome and genome approaches, we identified and reported novel, clinically significant variants in testis-expressed gene 15 (TEX15), in unrelated men with spermatogenic failure (SPGF). Results: TEX15 mediates double strand break repair during meiosis. Recessive loss-of-function (LOF) TEX15 mutations are associated with SPGF in humans and knockout male mice are infertile. We expand earlier reports documenting heterogeneous allelic pathogenic TEX15 variants that cause a range of SPGF phenotypes from oligozoospermia (low sperm) to nonobstructive azoospermia (no sperm) with meiotic arrest and report the prevalence of 0.6% of TEX15 variants in our patient cohort. Among identified possible LOF variants, one homozygous missense substitution c.6835G>A (p.Ala2279Thr) co-segregated with cryptozoospermia in a family with SPGF. Additionally, we observed numerous cases of inferred in trans compound heterozygous variants in TEX15 among unrelated individuals with varying degrees of SPGF. Variants included splice site, insertions/deletions (indels), and missense substitutions, many of which resulted in LOF effects (i.e., frameshift, premature stop, alternative splicing, or potentially altered posttranslational modification sites). Conclusion: In conclusion, we performed an extensive genomic study of familial and sporadic SPGF and identified potentially damaging TEX15 variants in 7 of 1097 individuals of our combined cohorts. We hypothesize that SPGF phenotype severity is dictated by individual TEX15 variant's impact on structure and function. Resultant LOFs likely have deleterious effects on crossover/recombination in meiosis. Our findings support the notion of increased gene variant frequency in SPGF and its genetic and allelic heterogeneity as it relates to complex disease such as male infertility.

Recent progress in understanding the genomic architecture of sexual conflict

Genomic conflict between the sexes over shared traits is widely assumed to be resolved through the evolution of sex-biased expression and the subsequent emergence of sexually dimorphic phenotypes. However, while there is support for a broad relationship between genome-wide patterns of expression level and sexual conflict, recent studies suggest that sex differences in the nature and strength of interactions between loci are instead key to conflict resolution. Furthermore, the advent of new technologies for measuring and perturbing expression means we now have much more power to detect genomic signatures of sexual conflict. Here, we review our current understanding of the genomic architecture of sexual conflict in the light of these new studies and highlight the potential for novel approaches to address outstanding knowledge gaps.

Single-Cell RNA Sequencing Reveals Atlas of Yak Testis Cells

Spermatogenesis is a complex process that involves proliferation and differentiation of diploid male germ cells into haploid flagellated sperm and requires intricate interactions between testicular somatic cells and germ cells. The cellular heterogeneity of this process presents a challenge in analyzing the different cell types at various developmental stages. Single-cell RNA sequencing (scRNA-seq) provides a useful tool for exploring cellular heterogeneity. In this study, we performed a comprehensive and unbiased single-cell transcriptomic study of spermatogenesis in sexually mature 4-year-old yak using 10× Genomics scRNA-seq. Our scRNA-seq analysis identified six somatic cell types and various germ cells, including spermatogonial stem cells, spermatogonia, early-spermatocytes, late-spermatocytes, and spermatids in yak testis. Pseudo-timing analysis showed that Leydig and myoid cells originated from common progenitor cells in yaks. Moreover, functional enrichment analysis demonstrated that the top expressed genes in yak testicular somatic cells were significantly enriched in the cAMP signaling pathway, PI3K-Akt signaling pathway, MAPK signaling pathway, and ECM receptor interactions. Throughout the spermatogenesis process, genes related to spermatogenesis, cell differentiation, DNA binding, and ATP binding were expressed. Using immunohistochemical techniques, we identified candidate marker genes for spermatogonial stem cells and Sertoli cells. Our research provides new insights into yak spermatogenesis and the development of various types of cells in the testis, and presents more reliable marker proteins for in vitro culture and identification of yak spermatogonial stem cells in the later stage.

Once upon a Testis: The Tale of Cyclic Nucleotide Phosphodiesterase in Testicular Cancers

Phosphodiesterases are key regulators that fine tune the intracellular levels of cyclic nucleotides, given their ability to hydrolyze cAMP and cGMP. They are critical regulators of cAMP/cGMP-mediated signaling pathways, modulating their downstream biological effects such as gene expression, cell proliferation, cell-cycle regulation but also inflammation and metabolic function. Recently, mutations in PDE genes have been identified and linked to human genetic diseases and PDEs have been demonstrated to play a potential role in predisposition to several tumors, especially in cAMP-sensitive tissues. This review summarizes the current knowledge and most relevant findings regarding the expression and regulation of PDE families in the testis focusing on PDEs role in testicular cancer development.

Novel Candidate Regulators and Developmental Trajectory of Pituitary Thyrotropes

The pituitary gland regulates growth, metabolism, reproduction, the stress response, uterine contractions, lactation, and water retention. It secretes hormones in response to hypothalamic input, end organ feedback, and diurnal cues. The mechanisms by which pituitary stem cells are recruited to proliferate, maintain quiescence, or differentiate into specific cell types, especially thyrotropes, are not well understood. We used single-cell RNA sequencing in juvenile P7 mouse pituitary cells to identify novel factors in pituitary cell populations, with a focus on thyrotropes and rare subtypes. We first observed cells coexpressing markers of both thyrotropes and gonadotropes, such as Pou1f1 and Nr5a1. This was validated in vivo by both immunohistochemistry and lineage tracing of thyrotropes derived from Nr5a1-Cre; mTmG mice and demonstrates that Nr5a1-progenitors give rise to a proportion of thyrotropes during development. Our data set also identifies novel factors expressed in pars distalis and pars tuberalis thyrotropes, including the Shox2b isoform in all thyrotropes and Sox14 specifically in Pou1f1-negative pars tuberalis thyrotropes. We have therefore used single-cell transcriptomics to determine a novel developmental trajectory for thyrotropes and potential novel regulators of thyrotrope populations.

Single-cell transcriptomic dissection of the toxic impact of di(2-ethylhexyl) phthalate on immature testicular development at the neonatal stage

Di(2-ethylhexyl) phthalate (DEHP) early exposure leads to immature testicular injury, and we aimed to utilize single-cell RNA (scRNA) sequencing to comprehensively assess the toxic effect of DEHP on testicular development. Therefore, we gavaged pregnant C57BL/6 mice with 750 mg/kg body weight DEHP from gestational day 13.5 to delivery and performed scRNA sequencing of neonatal testes at postnatal day 5.5. The results revealed the gene expression dynamics in testicular cells. DEHP disrupted the developmental trajectory of germ cells and the balance between the self-renewal and differentiation of spermatogonial stem cells. Additionally, DEHP caused an abnormal developmental trajectory, cytoskeletal damage and cell cycle arrest in Sertoli cells; disrupted the metabolism of testosterone in Leydig cells; and disturbed the developmental trajectory in peritubular myoid cells. Elevated oxidative stress and excessive apoptosis mediated by p53 were observed in almost all testicular cells. The intercellular interactions among four cell types were altered, and biological processes related to glial cell line-derived neurotrophic factor (GDNF), transforming growth factor-β (TGF-β), NOTCH, platelet-derived growth factor (PDGF) and WNT signaling pathways were enriched after DEHP treatment. These findings systematically describe the damaging effects of DEHP on the immature testes and provide substantial novel insights into the reproductive toxicity of DEHP.

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

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

Transcriptomic analysis of the Non-Obstructive Azoospermia (NOA) to address gene expression regulation in human testis

There is a need to understand the molecular basis of testes under Non-Obstructive Azoospermia (NOA), a state of failed spermatogenesis. There has been a lack of attention to the transcriptome at the level of alternatively spliced mRNAs (iso-mRNAs) and the mechanism of gene expression regulation. Hence, we aimed to establish a reliable iso-mRNA profile of NOA-testes, and explore molecular mechanisms - especially those related to gene expression regulation. We sequenced mRNAs from testicular samples of donors with complete spermatogenesis (control samples) and a failure of spermatogenesis (NOA samples). We identified differentially expressed genes and their iso-mRNAs via standard NGS data analyses. We then listed these iso-mRNAs hierarchically based on the extent of consistency of differential quantities across samples and groups, and validated the lists via RT-qPCRs (for 80 iso-mRNAs). In addition, we performed extensive bioinformatic analysis of the splicing features, domains, interactions, and functions of differentially expressed genes and iso-mRNAs. Many top-ranking down-regulated genes and iso-mRNAs, i.e., those down-regulated more consistently across the NOA samples, are associated with mitosis, replication, meiosis, cilium, RNA regulation, and post-translational modifications such as ubiquitination and phosphorylation. Most down-regulated iso-mRNAs correspond to full-length proteins that include all expected domains. The predominance of alternative promoters and termination sites in these iso-mRNAs indicate their gene expression regulation via promoters and UTRs. We compiled a new, comprehensive list of human transcription factors (TFs) and used it to identify TF-'TF gene' interactions with potential significance in down-regulating genes under the NOA condition. The results indicate that RAD51 suppression by HSF4 prevents SP1-activation, and SP1, in turn, could regulate multiple TF genes. This potential regulatory axis and other TF interactions identified in this study could explain the down-regulation of multiple genes in NOA-testes. Such molecular interactions may also have key regulatory roles during normal human spermatogenesis.

Retinoic acid-induced differentiation of porcine prospermatogonia in vitro

Spermatogenesis is an intricate developmental process occurring in testes by which spermatogonial stem cells (SSCs) self-renew and differentiate into mature sperm. The molecular mechanisms for SSC self-renewal and differentiation, while have been well studied in mice, may differ between mice and domestic animals including pigs. To gain knowledge about the molecular mechanisms for porcine SSC self-renewal and differentiation that have so far been poorly understood, here we isolated and enriched prospermatogonia from neonatal porcine testes, and exposed the cells to retinoic acid, a direct inducer for spermatogonial differentiation. We then identified that retinoic acid could induce porcine prospermatogonial differentiation, which was accompanied by a clear transcriptomic alteration, as revealed by the RNA-sequencing analysis. We also compared retinoic acid-induced in vitro porcine spermatogonial differentiation with the in vivo process, and compared retinoic acid-induced in vitro spermatogonial differentiation between pigs and mice. Furthermore, we analyzed retinoic acid-induced differentially expressed long non-coding RNAs (lncRNAs), and demonstrated that a pig-specific lncRNA, lncRNA-106504875, positively regulated porcine spermatogonial proliferation by targeting the core transcription factor ZBTB16. Taken together, these results would help to elucidate the roles of retinoic acid in porcine spermatogonial differentiation, thereby contributing to further knowledge about the molecular mechanisms underlying porcine SSC development and, in the long run, to optimization of both long-term culture and induced differentiation systems for porcine SSCs.

Dynamic cellular and molecular characteristics of spermatogenesis in the viviparous marine teleost Sebastes schlegelii†

Spermatogenesis is a dynamic cell developmental process that is essential for reproductive success. Vertebrates utilize a variety of reproductive strategies, including sperm diversity, and internal and external fertilization. Research on the cellular and molecular dynamic changes involved in viviparous teleost spermatogenesis, however, is currently lacking. Here, we combined cytohistology, 10 × genomic single-cell RNA-seq, and transcriptome technology to determine the dynamic development characteristics of the spermatogenesis of Sebastes schlegelii. The expressions of lhcgr (Luteinizing hormone/Choriogonadotropin receptor), fshr (follicle-stimulating hormone receptor), ar (androgen receptor), pgr (progesterone receptor), and cox (cyclo-oxygen-ase), as well as the prostaglandin E and F levels peaked during the maturation period, indicating that they were important for sperm maturation and mating. Fifteen clusters were identified based on the 10 × genomic single-cell results. The cell markers of the sub-cluster were identified by their upregulation; piwil, dazl, and dmrt1 were upregulated and identified as spermatogonium markers, and sycp1/3 and spo11 were identified as spermatocyte markers. For S. schlegelii, the sperm head nucleus was elongated (spherical to streamlined in shape), which is a typical characteristic for sperm involved in internal fertilization. We also identified a series of crucial genes associated with spermiogenesis, such as spata6, spag16, kif20a, trip10, and klf10, while kif2c, kifap3, fez2, and spaca6 were found to be involved in nucleus elongation. The results of this study will enrich our cellular and molecular knowledge of spermatogenesis and spermiogenesis in fish that undergo internal fertilization.

Conditional deletion of Hspa5 leads to spermatogenesis failure and male infertility in mice

Heat shock proteins (HSPs) have important roles in different developmental stages of spermatogenesis. The heat shock 70 kDa protein 5 (HSPA5) is an important component of the unfolded protein response that promotes cell survival under endoplasmic reticulum (ER) stress conditions. In this study, we explored the function of HSPA5 in spermatogenesis, by generating a germ cell-specific deletion mutant of the Hspa5 gene (conditional knockout of the Hspa5 gene, Hspa5-cKO) using CRISPR/Cas9 technology and the Cre/Loxp system. Hspa5 knockout resulted in severe germ cell loss and vacuolar degeneration of seminiferous tubules, leading to complete arrest of spermatogenesis, testicular atrophy, and male infertility in adult mice. Furthermore, defects occurred in the spermatogenic epithelium of Hspa5-cKO mice as early as Cre recombinase expression. Germ cell ablation of Hspa5 impaired spermatogonia proliferation and differentiation from post-natal day 7 (P7) to P10, which led to a dramatic reduction of differentiated spermatogonia, compromised meiosis, and led to impairment of testis development and the disruption of the first wave of spermatogenesis. Consistent with these results, single-cell RNA sequencing (scRNA-seq) analysis showed that germ cells, especially differentiated spermatogonia, were dramatically reduced in Hspa5-cKO testes compared with controls at P10, further confirming that HSPA5 is crucial for germ cell development. These results suggest that HSPA5 is indispensable for normal spermatogenesis and male reproduction in mice.

Germ Cell Maintenance and Sustained Testosterone and Precursor Hormone Production in Human Prepubertal Testis Organ Culture with Tissues from Boys 7 Years+ under Conditions from Adult Testicular Tissue

Human prepubertal testicular tissues are rare, but organ culture conditions to develop a system for human in vitro-spermatogenesis are an essential option for fertility preservation in prepubertal boys subjected to gonadotoxic therapy. To avoid animal testing in line with the 3Rs principle, organ culture conditions initially tested on human adult testis tissue were applied to prepubertal samples (n = 3; patient ages 7, 9, and 12 years). Tissues were investigated by immunostaining and transmission electron microscopy (TEM), and the collected culture medium was profiled for steroid hormones by liquid chromatography-tandem mass spectrometry (LC-MS/MS). Culture conditions proved suitable for prepubertal organ culture since SSCs and germ cell proliferation could be maintained until the end of the 3-week-culture. Leydig cells (LCs) were shown to be competent for steroid hormone production. Three additional testis tissues from boys of the same age were examined for the number of germ cells and undifferentiated spermatogonia (SPG). Using TEM micrographs, eight tissues from patients aged 1.5 to 13 years were examined, with respect to the sizes of mitochondria (MT) in undifferentiated SPG and compared with those from two adult testicular tissues. Mitochondrial sizes were shown to be comparable between adults and prepubertal boys from approximately 7 years of age, which suggests the transition of SSCs from normoxic to hypoxic metabolism at about or before this time period.

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

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

scONE-seq: A single-cell multi-omics method enables simultaneous dissection of phenotype and genotype heterogeneity from frozen tumors

Single-cell multi-omics can provide a unique perspective on tumor cellular heterogeneity. Most previous single-cell whole-genome RNA sequencing (scWGS-RNA-seq) methods demonstrate utility with intact cells from fresh samples. Among them, many are not applicable to frozen samples that cannot produce intact single-cell suspensions. We have developed scONE-seq, a versatile scWGS-RNA-seq method that amplifies single-cell DNA and RNA without separating them from each other and hence is compatible with frozen biobanked samples. We benchmarked scONE-seq against existing methods using fresh and frozen samples to demonstrate its performance in various aspects. We identified a unique transcriptionally normal-like tumor clone by analyzing a 2-year frozen astrocytoma sample, demonstrating that performing single-cell multi-omics interrogation on biobanked tissue by scONE-seq could enable previously unidentified discoveries in tumor biology.

Cross-species cell-type assignment from single-cell RNA-seq data by a heterogeneous graph neural network

Cross-species comparative analyses of single-cell RNA sequencing (scRNA-seq) data allow us to explore, at single-cell resolution, the origins of the cellular diversity and evolutionary mechanisms that shape cellular form and function. Cell-type assignment is a crucial step to achieve that. However, the poorly annotated genome and limited known biomarkers hinder us from assigning cell identities for nonmodel species. Here, we design a heterogeneous graph neural network model, CAME, to learn aligned and interpretable cell and gene embeddings for cross-species cell-type assignment and gene module extraction from scRNA-seq data. CAME achieves significant improvements in cell-type characterization across distant species owing to the utilization of non-one-to-one homologous gene mapping ignored by early methods. Our large-scale benchmarking study shows that CAME significantly outperforms five classical methods in terms of cell-type assignment and model robustness to insufficiency and inconsistency of sequencing depths. CAME can transfer the major cell types and interneuron subtypes of human brains to mouse and discover shared cell-type-specific functions in homologous gene modules. CAME can align the trajectories of human and macaque spermatogenesis and reveal their conservative expression dynamics. In short, CAME can make accurate cross-species cell-type assignments even for nonmodel species and uncover shared and divergent characteristics between two species from scRNA-seq data.

Perspectives: Approaches for Studying Livestock Spermatogonia

At present, the knowledge base on characteristics and biology of spermatogonia in livestock is limited in comparison to rodents, yet the importance of studying these cells for comparative species analysis and enhancing reproductive capacity in food animals is high. Previous studies have established that although many core attributes of organ physiology and mechanisms governing essential cellular functions are conserved across eutherians, significant differences exist between mice and higher order mammals. In this chapter, we briefly discuss distinguishing aspects of testicular anatomy and the spermatogenic lineage in livestock and critical considerations for studying spermatogonial stem cell biology in these species.

Comparative single-cell transcriptomic analysis reveals differences in signaling pathways in gonadal primordial germ cells between chicken (Gallus gallus) and zebra finch (Taeniopygia guttata)

Primordial germ cells (PGCs) have been used in avian genetic resource conservation and transgenic animal production. Despite their potential applications to numerous avian taxa facing extinction due to habitat loss and degradation, research has largely focused on poultry, such as chickens, in part owing to the difficulty in obtaining intact PGCs from other species. Recently, phenotypic differences between PGCs of chicken and zebra finch, a wild bird with vocal learning, in early embryonic development have been reported. In this study, we used advanced single-cell RNA sequencing (scRNA-seq) technology to evaluate zebra finch and chicken PGCs and surrounding cells, and to identify species-specific characteristics. We constructed single-cell transcriptome landscapes of chicken gonadal PGCs for a comparison with previously reported scRNA-seq data for zebra finch. We identified interspecific differences in several signaling pathways in gonadal PGCs and somatic cells. In particular, NODAL and insulin signaling pathway activity levels were higher in zebra finch than in chickens, whereas activity levels of the downstream FGF signaling pathway, involved in the proliferation of chicken PGCs, were higher in chickens. This study is the first cross-species single-cell transcriptomic analysis targeting birds, revealing differences in germ cell development between phylogenetically distant Galliformes and Passeriformes. Our results provide a basis for understanding the reproductive physiology of avian germ cells and for utilizing PGCs in the restoration of endangered birds and the production of transgenic birds.

Optimization of Testicular Fixation-Embedding Techniques for Improved Evaluation of Mammalian Spermatogonial Morphology and Function

Numerous methods have been successfully used to evaluate mammalian spermatogonial biology However, the conventional light microscopy assays present a challenge in precisely identifying spermatogonial phenotypes, which can result in discrepancies between molecular and morphological findings. Such precise association could lead to a more robust interpretation of spermatogonial activity in steady-state spermatogenesis, which may facilitate the translation from basic research to clinical applications. In this chapter, we present two histological processing methods that enable a comprehensive analysis of spermatogonial morphology and function, involving fixation of mammalian testicular tissue in glutaraldehyde and embedding in plastic resin. These techniques have proven to be effective in light microscopy studies.

The molecular evolution of spermatogenesis across mammals

The testis produces gametes through spermatogenesis and evolves rapidly at both the morphological and molecular level in mammals^1-6, probably owing to the evolutionary pressure on males to be reproductively successful⁷. However, the molecular evolution of individual spermatogenic cell types across mammals remains largely uncharacterized. Here we report evolutionary analyses of single-nucleus transcriptome data for testes from 11 species that cover the three main mammalian lineages (eutherians, marsupials and monotremes) and birds (the evolutionary outgroup), and include seven primates. We find that the rapid evolution of the testis was driven by accelerated fixation rates of gene expression changes, amino acid substitutions and new genes in late spermatogenic stages, probably facilitated by reduced pleiotropic constraints, haploid selection and transcriptionally permissive chromatin. We identify temporal expression changes of individual genes across species and conserved expression programs controlling ancestral spermatogenic processes. Genes predominantly expressed in spermatogonia (germ cells fuelling spermatogenesis) and Sertoli (somatic support) cells accumulated on X chromosomes during evolution, presumably owing to male-beneficial selective forces. Further work identified transcriptomal differences between X- and Y-bearing spermatids and uncovered that meiotic sex-chromosome inactivation (MSCI) also occurs in monotremes and hence is common to mammalian sex-chromosome systems. Thus, the mechanism of meiotic silencing of unsynapsed chromatin, which underlies MSCI, is an ancestral mammalian feature. Our study illuminates the molecular evolution of spermatogenesis and associated selective forces, and provides a resource for investigating the biology of the testis across mammals.

Perspectives: Methods for Evaluating Primate Spermatogonial Stem Cells

Mammalian spermatogenesis is a complex, highly productive process generating millions of sperm per day. Spermatogonial stem cells (SSCs) are at the foundation of spermatogenesis and can either self-renew, producing more SSCs, or differentiate to initiate spermatogenesis and produce sperm. The biological potential of SSCs to produce and maintain spermatogenesis makes them a promising tool for the treatment of male infertility. However, translating knowledge from rodents to higher primates (monkeys and humans) is challenged by different vocabularies that are used to describe stem cells and spermatogenic lineage development in those species. Furthermore, while rodent SSCs are defined by their biological potential to produce and maintain spermatogenesis in a transplant assay, there is no equivalent routine and accessible bioassay to test monkey and human SSCs or replicate their functions in vitro. This chapter describes progress characterizing, isolating, culturing, and transplanting SSCs in higher primates.

Single-cell RNA sequencing uncovers dynamic roadmap and cell-cell communication during buffalo spermatogenesis

Spermatogenesis carries the task of precise intergenerational transmission of genetic information from the paternal genome and involves complex developmental processes regulated by the testicular microenvironment. Studies performed mainly in mouse models have established the theoretical basis for spermatogenesis, yet the wide interspecies differences preclude direct translation of the findings, and farm animal studies are progressing slowly. More than 32,000 cells from prepubertal (3-month-old) and pubertal (24-month-old) buffalo testes were analyzed by using single-cell RNA sequencing (scRNA-seq), and dynamic gene expression roadmaps of germ and somatic cell development were generated. In addition to identifying the dynamic processes of sequential cell fate transitions, the global cell-cell communication essential to maintain regular spermatogenesis in the buffalo testicular microenvironment was uncovered. The findings provide the theoretical basis for establishing buffalo germline stem cells in vitro or culturing organoids and facilitating the expansion of superior livestock breeding.

The impact of paternal age on new mutations and disease in the next generation

Advanced paternal age is associated with an increased risk of fathering children with genetic disorders and other adverse reproductive consequences. However, the mechanisms underlying this phenomenon remain largely unexplored. In this review, we focus on the impact of paternal age on de novo mutations that are an important contributor to genetic disease and can be studied both indirectly through large-scale sequencing studies and directly in the tissue in which they predominantly arise-the aging testis. We discuss the recent data that have helped establish the origins and frequency of de novo mutations, and highlight experimental evidence about the close link between new mutations, parental age, and genetic disease. We then focus on a small group of rare genetic conditions, the so-called "paternal age effect" disorders that show a strong association between paternal age and disease prevalence, and discuss the underlying mechanism ("selfish selection") and implications of this process in more detail. More broadly, understanding the causes and consequences of paternal age on genetic risk has important implications both for individual couples and for public health advice given that the average age of fatherhood is steadily increasing in many developed nations.

Determination of the nanoparticle- and cell-specific toxicological mechanisms in 3D liver spheroids using scRNAseq analysis

Engineered nanomaterials (ENMs) are commonly used in consumer products, allowing exposure to target organs such as the lung, liver, and skin that could lead to adverse health effects in humans. To better reflect on toxicological effects in liver cells, it is important to consider the contribution of hepatocyte morphology, function, and intercellular interactions in a dynamic 3D microenvironment. Herein, we used a 3D liver spheroid model containing hepatocyte and Kupffer cells (KCs) to study the effects of three different material compositions, namely vanadium pentoxide (V2O5), titanium dioxide (TiO2), or graphene oxide (GO). Additionally, we used single-cell RNA sequencing (scRNAseq) to determine the nanoparticle (NP) and cell-specific toxicological responses. A general finding was that hepatocytes exhibit more variation in gene expression and adaptation of signaling pathways than KCs. TNF-α production tied to the NF-κB pathway was a commonly affected pathway by all NPs while impacts on the metabolic function of hepatocytes were unique to V2O5. V2O5 NPs also showed the largest number of differentially expressed genes in both cell types, many of which are related to pro-inflammatory and apoptotic response pathways. There was also evidence of mitochondrial ROS generation and caspase-1 activation after GO and V2O5 treatment, in association with cytokine production. All considered, this study provides insight into the impact of nanoparticles on gene responses in key liver cell types, providing us with a scRNAseq platform that can be used for high-content screening of nanomaterial impact on the liver, for use in biosafety and biomedical applications.

Decoding the Spermatogenesis Program: New Insights from Transcriptomic Analyses

Spermatogenesis is a complex differentiation process coordinated spatiotemporally across and along seminiferous tubules. Cellular heterogeneity has made it challenging to obtain stage-specific molecular profiles of germ and somatic cells using bulk transcriptomic analyses. This has limited our ability to understand regulation of spermatogenesis and to integrate knowledge from model organisms to humans. The recent advancement of single-cell RNA-sequencing (scRNA-seq) technologies provides insights into the cell type diversity and molecular signatures in the testis. Fine-grained cell atlases of the testis contain both known and novel cell types and define the functional states along the germ cell developmental trajectory in many species. These atlases provide a reference system for integrated interspecies comparisons to discover mechanistic parallels and to enable future studies. Despite recent advances, we currently lack high-resolution data to probe germ cell-somatic cell interactions in the tissue environment, but the use of highly multiplexed spatial analysis technologies has begun to resolve this problem. Taken together, recent single-cell studies provide an improvedunderstanding of gametogenesis to examine underlying causes of infertility and enable the development of new therapeutic interventions.