The Novel Key Genes of Non-obstructive Azoospermia Affect Spermatogenesis: Transcriptomic Analysis Based on RNA-Seq and scRNA-Seq Data

Abstract

Non-obstructive azoospermia (NOA) is one of the most important causes of male infertility. It is mainly characterized by the absence of sperm in semen repeatedly or the number of sperm is small and not fully developed. At present, its pathogenesis remains largely unknown. The goal of this study is to identify hub genes that might affect biomarkers related to spermatogenesis. Using the clinically significant transcriptome and single-cell sequencing data sets on the Gene Expression Omnibus (GEO) database, we identified candidate hub genes related to spermatogenesis. Based on them, we performed Gene Ontology (GO) functional enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment pathway analyses, protein-protein interaction (PPI) network analysis, principal component analysis (PCA), cell cluster analysis, and pseudo-chronological analysis. We identified a total of 430 differentially expressed genes, of which three have not been reported related to spermatogenesis (C22orf23, TSACC, and TTC25), and the expression of these three hub genes was different in each type of sperm cells. The results of the pseudo-chronological analysis of the three hub genes indicated that TTC25 was in a low expression state during the whole process of sperm development, while the expression of C22orf23 had two fluctuations in the differentiating spermatogonia and late primary spermatocyte stages, and TSACC showed an upward trend from the spermatogonial stem cell stage to the spermatogenesis stage. Our research found that the three hub genes were different in the trajectory of sperm development, indicating that they might play important roles in different sperm cells. This result is of great significance for revealing the pathogenic mechanism of NOA and further research.