SUMMARY Spermatogenesis is a complex and dynamic cellular differentiation process critical to male reproduction and sustained by spermatogonial stem cells (SSCs). Although patterns of gene expression have been described for aggregates of certain spermato- genic cell types, the full continuum of gene expression patterns underlying ongoing spermatogenesis in steady state was previously unclear. Here, we catalog single-cell transcriptomes for >62,000 individual spermatogenic cells from immature (postnatal day 6) and adult male mice and adult men. This allowed us to resolve SSC and progenitor spermatogonia, elucidate the full range of gene expression changes during male meiosis and spermiogenesis, and derive unique gene expression signatures for multiple mouse and human spermatogenic cell types and/or subtypes. These transcriptome datasets provide an information-rich resource for studies of SSCs, male meiosis, testicular cancer, male infertility, or contraceptive development, as well as a gene expression roadmap to be emulated in efforts to achieve spermatogenesis in vitro.

Spermatogonial stem cells (SSCs) are a subset of undifferentiated spermatogonia responsible for ongoing spermatogenesis in mammalian testes. Spermatogonial stem cells arise from morphologically homogeneous prospermatogonia, but growing evidence suggests that only a subset of prospermatogonia develops into the foundational SSC pool. This predicts that subtypes of undifferentiated spermatogonia with discrete mRNA and protein signatures should be distinguishable in neonatal testes. We used single-cell quantitative RT-PCR to examine mRNA levels of 172 genes in individual spermatogonia from 6-day postnatal (P6) mouse testes. Cells enriched from P6 testes using the StaPut or THY1(+) magnetic cell sorting methods exhibited considerable heterogeneity in the abundance of specific germ cell and stem cell mRNAs, segregating into one somatic and three distinct spermatogonial clusters. However, P6 Id4-eGFP(+) transgenic spermatogonia, which are known to be enriched for SSCs, were more homogeneous in their mRNA levels, exhibiting uniform levels for the majority of genes examined (122 of 172). Interestingly, these cells displayed nonuniform (50 of 172) expression of a smaller cohort of these genes, suggesting there is substantial heterogeneity even within the Id4-eGFP(+) population. Further, although immunofluorescence staining largely demonstrated conformity between mRNA and protein levels, some proteins were observed in patterns that were disparate from those detected for the corresponding mRNAs in Id4-eGFP(+) spermatogonia (e.g., Kit, Sohlh2, Stra8), suggesting additional heterogeneity is introduced at the posttranscriptional level. Taken together, these data demonstrate the existence of multiple spermatogonial subtypes in P6 mouse testes and raise the intriguing possibility that these subpopulations may correlate with the development of functionally distinct spermatogenic cell types.

Summary Spermatogonial stem cells (SSCs) must balance self-renewal with production of transit-amplifying progenitors that differentiate in response to retinoic acid (RA) before entering meiosis. This self-renewal vs. differentiation spermatogonial fate decision is critical for maintaining tissue homeostasis, as imbalances cause spermatogenesis defects that can lead to human testicular cancer or infertility. A great deal of effort has been exerted to understand how the SSC population is maintained. In contrast, little is known about the essential program of differentiation initiated by retinoic acid (RA) that precedes meiosis, and the pathways and proteins involved are poorly defined. We recently reported a novel role for RA in stimulating the PI3/AKT/mTOR kinase signaling pathway to activate translation of repressed mRNAs such as Kit. Here, we examined the requirement for mTOR complex 1 (mTORC1) in mediating the RA signal to direct spermatogonial differentiation in the neonatal testis. We found that in vivo inhibition of mTORC1 by rapamycin blocked spermatogonial differentiation, which led to an accumulation of undifferentiated spermatogonia. In addition, rapamycin also blocked the RA-induced translational activation of mRNAs encoding KIT, SOHLH1, and SOHLH2 without affecting expression of STRA8. These findings highlight dual roles for RA in germ cell development – transcriptional activation of genes, and kinase signaling to stimulate translation of repressed messages required for spermatogonial differentiation.

Precise separation of spermatogonial stem cells (SSCs) from progenitor spermatogonia that lack stem cell activity and are committed to differentiation remains a challenge. To distinguish between these spermatogonial subtypes, we identified genes that exhibited bimodal mRNA levels at the single-cell level among undifferentiated spermatogonia from Postnatal Day 6 mouse testes, including Tspan8, Epha2, and Pvr, each of which encode cell surface proteins useful for cell selection. Transplantation studies provided definitive evidence that a TSPAN8-high subpopulation is enriched for SSCs. RNA-seq analyses identified genes differentially expressed between TSPAN8-high and -low subpopulations that clustered into multiple biological pathways potentially involved in SSC renewal or differentiation, respectively. Methyl-seq analysis identified hypomethylated domains in the promoters of these genes in both subpopulations that colocalized with peaks of histone modifications defined by ChIP-seq analysis. Taken together, these results demonstrate functional heterogeneity among mouse undifferentiated spermatogonia and point to key biological characteristics that distinguish SSCs from progenitor spermatogonia.

In the mammalian testis, sustained spermatogenesis relies on spermatogonial stem cells (SSCs); their progeny either remain as stem cells (self-renewal) or proliferate and differentiate to enter meiosis in response to retinoic acid (RA). Here, we sought to uncover elusive mechanisms regulating a key switch fundamental to spermatogonial fate: the capacity of spermatogonia to respond to RA. Using the developing mouse testis as a model, we found that spermatogonia and precursor prospermatogonia exhibit a heterogeneous capacity to respond to RA with at least two underlying causes. First, progenitor spermatogonia are prevented from responding to RA by catabolic activity of cytochrome P450 family 26 enzymes. Second, a smaller subset of undifferentiated spermatogonia enriched for SSCs exhibit catabolism-independent RA insensitivity. Moreover, for the first time, we observed that precursor prospermatogonia are heterogeneous and comprise subpopulations that exhibit the same differential RA responsiveness found in neonatal spermatogonia. We propose a novel model by which mammalian prospermatogonial and spermatogonial fates are regulated by their intrinsic capacity to respond (or not) to the differentiation signal provided by RA before, and concurrent with, the initiation of spermatogenesis.