Integrated molecular-network analysis reveals infertility-associated key genes and transcription factors in the non-obstructive azoospermia

Background Male infertility is a multifactorial reproductive health problem with complex causes. Non-obstructive azoospermia (NOA) is characterized by failure of spermatogenesis, leading to the absence of spermatozoa in ejaculates. The molecular mechanism underlying the NOA is still not well understood. Objectives This study aims to identify the key genes involved in male infertility that could be a potential biomarker in the diagnosis and prognosis of azoospermia. Study design The microarray expression profiles dataset GSE45885 and GSE45887 were downloaded from the NCBI's Gene Expression Omnibus (GEO) database and analyzed for male infertility-associated differentially expressed genes (DEGs) using the GEO2R tool. The common DEGs between the two datasets were combined and their protein–protein interaction (PPI) network was constructed using Cytoscape to reveal the hub genes by topology and module analysis. In addition, transcription factors (TFs) and protein kinases regulating the hub genes were identified using the X2K tool. Then, the expression of the hub genes was validated by analyzing the GSE190752 microarray dataset. Further, the PPI network was screened for biological roles and enriched pathways using DAVID software. Results About 256 DEGs associated with NOA were identified and constructed the PPI network to find the infertility-associated proteins. The biological processes linked with these proteins were spermatogenesis, cell differentiation, flagellated sperm motility, and spermatid development. The topology and module analysis of the infertility-associated protein network identified the hub genes TEX38, FAM71F, PRR30, FAM166A, LYZL6, TPPP2, ARMC12, SPACA4, and FAM205A, which were found to be upregulated in the non-obstructive azoospermia. In addition, a total of 23 transcription factors and 3 protein kinases that are regulating these key hub genes were identified. Further these hub genes expression was validated using the microarray data and found that their expression was increased in the testicular biopsies obtained from NOA subjects, compared to healthy individuals. Conclusion The identified key genes and its associated transcription factors are known to regulate the infertility-related processes in the non-obstructive azoospermia. Also, the clinical sample-based microarray data validation for the expression of these key hub genes indicates their potentiality to develop them as diagnostic or prognostic biomarkers for NOA.