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.

Spermatogonial stem cells (SSCs) maintain spermatogenesis throughout adulthood through balanced self-renewal and differentiation, yet the regulatory logic of these fate decisions is poorly understood. The transcription factors Sal-like 4 (SALL4) and promyelocytic leukemia zinc finger (PLZF; also known as ZBTB16) are known to be required for normal SSC function, but their targets are largely unknown. ChIP-seq in mouse THY1+ spermatogonia identified 4176 PLZF-bound and 2696 SALL4-bound genes, including 1149 and 515 that were unique to each factor, respectively, and 1295 that were bound by both factors. PLZF and SALL4 preferentially bound gene promoters and introns, respectively. Motif analyses identified putative PLZF and SALL4 binding sequences, but rarely both at shared sites, indicating significant non-autonomous binding in any given cell. Indeed, the majority of PLZF/SALL4 shared sites contained only PLZF motifs. SALL4 also bound gene introns at sites containing motifs for the differentiation factor DMRT1. Moreover, mRNA levels for both unique and shared target genes involved in both SSC self-renewal and differentiation were suppressed following SALL4 or PLZF knockdown. Together, these data reveal the full profile of PLZF and SALL4 regulatory targets in undifferentiated spermatogonia, including SSCs, which will help elucidate mechanisms controlling the earliest cell fate decisions in spermatogenesis.

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.

Mycobacterial infections in laboratory zebrafish (Danio rerio) are common and widespread in research colonies. Mycobacteria within free living amoebae have been shown to be transmission vectors for mycobacteriosis. Paramecium caudatum are commonly used as a first food for zebrafish, and we investigated this ciliate’s potential to serve as a vector of Mycobacterium marinum and M. chelonae. The ability of live P. caudatum to transmit these mycobacteria to larval, juvenile and adult zebrafish was evaluated. Infections were defined by histologic observation of granulomas containing acid-fast bacteria in extraintestinal locations. In both experiments, fish fed paramecia containing mycobacteria became infected at a higher incidence than controls. Larvae (exposed at 4 days post hatch) fed paramecia with M. marinum exhibited an incidence of 30% (24/80) and juveniles (exposed at 21 days post hatch) showed 31% incidence (14/45). Adult fish fed a gelatin food matrix containing mycobacteria within paramecia or mycobacteria alone for 2 wk resulted in infections when examined 8 wk after exposure as follows: M. marinum OSU 214 47% (21/45), M. marinum CH 47% (9/19), M. chelonae 38% (5/13). In contrast, fish feed mycobacteria alone in this diet did not become infected, except for 2 fish (5%) in the M. marinum OSU 214 low dose group. These results demonstrate that Paramecium caudatum can act as a vector for mycobacteria. This provides a useful animal model for evaluation of natural mycobacterial infections and demonstrates the possibility of mycobacterial transmission in zebrafish facilities via contaminated paramecia cultures.