Spermatogenesis is a highly regulated process that produces sperm to transmit genetic information to the next generation. Although extensively studied in mice, our current understanding of primate spermatogenesis is limited to populations defined by state-specific markers defined from rodent data. As between-species differences have been reported in the process duration and cellular differentiation hierarchy, it remains unclear how molecular markers and cell states are conserved or have diverged from mice to man. To address this challenge, we employ single-cell RNA-sequencing to identify transcriptional signatures of major germ and somatic cell-types of the testes in human, macaque and mice. This approach reveals differences in expression throughout spermatogenesis, including the stem/progenitor pool of spermatogonia, classical markers of differentiation, potential regulators of meiosis, the kinetics of RNA turnover during spermatid differentiation, and germ cell-soma communication. These datasets provide a rich foundation for future targeted mechanistic studies of primate germ cell development and in vitro gametogenesis. potential communications with the germline. Overall, this study provides the first single-cell comparative analysis of the spermatogenesis program between primates and rodents. Such a new resource is expected to improve our knowledge base for future studies of germ cell development in primates, and ultimately improve our understanding of the intrinsic and extrinsic evolutionary changes of the gametogenesis program. Knowledge gained from these data will inform fertility restoration efforts, including SSC culture and in vitro gametogenesis.
ResultsSingle-cell sequencing identifies major germ and somatic cell types of adult human and macaque testes. Using the Drop-seq platform, we generated single cell transcriptome data from ~14K and ~22K adult human and rhesus macaque testes, respectively. Using these datasets, our overall strategy was to first identify major cell types and states in the adult human and rhesus macaque testes separately, then combine them, along with our previous reported mouse data 12 , to define evolutionarily conserved or diverged programs within germ cells and the somatic cells of the testis (Figure 1A). The doublet rates for the human and macaque datasets is estimated to be <2%, therefore, allowing reliable analysis of individual cells (Figure S1A). After QC filtering (see Methods), 13,837 and 21,574 cells were retained from human and macaque, respectively, with an average of 3,210 unique molecular identifiers (UMIs) per cell, and 1,304 genes detected per cell. This sequencing depth and the number of detected genes was sufficient to define major cell types 13,14 . Systematic comparisons of technical batches (such as Human 1.1-1.5, Figure S1B-C) and biological replicates (4 humans and 5 macaques, Figure S1D-E) confirmed high batch-to-batch or individual-to-individual concordance, despite the variable proportions of cell types in different samples, likely due to differences i...