Integration and gene co-expression network analysis of scRNA-seq transcriptomes reveal heterogeneity and key functional genes in human spermatogenesis

• Published in: Scientific reports Vol. 11; no. 1; p. 19089
• Main Authors: Salehi, Najmeh, Karimi-Jafari, Mohammad Hossein, Totonchi, Mehdi, Amiri-Yekta, Amir
• Format: Journal Article
• Language: English
• Published: England Nature Publishing Group 27.09.2021; Nature Publishing Group UK; Nature Portfolio
• Subjects: Cell differentiation; Datasets as Topic; Gene Regulatory Networks; Genes; Genetic Heterogeneity; Genotypes; Heterogeneity; Histocompatibility antigen HLA; Humans; Infertility; Infertility, Male - genetics; Male; Molecular modelling; Nucleotide sequence; Phenotypes; RNA-Seq; Single-Cell Analysis; Spermatogenesis; Spermatogenesis - genetics; Spermatogonia; Stem cells; Transcriptomes
• ISSN: 2045-2322; 2045-2322
• DOI: 10.1038/s41598-021-98267-3

Spermatogenesis is a complex process of cellular division and differentiation that begins with spermatogonia stem cells and leads to functional spermatozoa production. However, many of the molecular mechanisms underlying this process remain unclear. Single-cell RNA sequencing (scRNA-seq) is used to sequence the entire transcriptome at the single-cell level to assess cell-to-cell variability. In this study, more than 33,000 testicular cells from different scRNA-seq datasets with normal spermatogenesis were integrated to identify single-cell heterogeneity on a more comprehensive scale. Clustering, cell type assignments, differential expressed genes and pseudotime analysis characterized 5 spermatogonia, 4 spermatocyte, and 4 spermatid cell types during the spermatogenesis process. The UTF1 and ID4 genes were introduced as the most specific markers that can differentiate two undifferentiated spermatogonia stem cell sub-cellules. The C7orf61 and TNP can differentiate two round spermatid sub-cellules. The topological analysis of the weighted gene co-expression network along with the integrated scRNA-seq data revealed some bridge genes between spermatogenesis's main stages such as DNAJC5B, C1orf194, HSP90AB1, BST2, EEF1A1, CRISP2, PTMS, NFKBIA, CDKN3, and HLA-DRA. The importance of these key genes is confirmed by their role in male infertility in previous studies. It can be stated that, this integrated scRNA-seq of spermatogenic cells offers novel insights into cell-to-cell heterogeneity and suggests a list of key players with a pivotal role in male infertility from the fertile spermatogenesis datasets. These key functional genes can be introduced as candidates for filtering and prioritizing genotype-to-phenotype association in male infertility.