Altered Gene Expression in the Testis of Infertile Patients with Nonobstructive Azoospermia

The molecular mechanism of nonobstructive azoospermia (NOA) remains unclear. The aim of this study was to identify gene expression changes in NOA patients and to explore potential biomarkers and therapeutic targets.

The gene expression profiles of GSE45885 and GSE145467 were collected from the Gene Expression Omnibus (GEO) database, and the differences between NOA and normal spermatogenesis were analyzed. Enrichment analysis was performed to explore biological functions for common differentially expressed genes (DEGs) in GSE45885 and GSE145467. Coexpression analysis of DEGs in GSE45885 was performed, and two modules with the highest correlation with NOA were screened. Key genes were then screened from the intersection genes of the two modules and common DEGs and PPI network. The expression of key genes was validated by quantitative real-time polymerase chain reaction (qRT-PCR) experiments. Finally, through miRTarBase, miRDB, and RAID, the miRNAs were predicted to regulate key genes, respectively.

A total of 345 common DEGs were identified and they were mainly related to spermatogenesis, insulin signaling pathway. Coexpression analysis of DEGs in GSE45885 yielded eight modules; MEblack and MEturquoise had the highest correlation with NOA. Six genes in MEturquoise and RNF141 in MEblack were identified as key genes. qRT-PCR experiments validated the differential expression of key genes between NOA and control. Furthermore, RNF141 was regulated by the largest number of miRNAs.

Our findings suggest that the significant change expression of key genes may be potential markers and therapeutic targets of NOA and may have some impact on the development of NOA.