Abstract:Background: Non-obstructive azoospermia (NOA) is a major contributor of male infertility. Herein, we used existing datasets to identify novel biomarkers for the diagnosis and prognosis of NOA, which could have great significance in the field of male infertility. Methods: NOA datasets were obtained from the Gene Expression Omnibus (GEO) database. CIBERSORT was utilized to analyze the distributions of 22 immune cell populations. Hub genes were identified by applying weighted gene co-expression network analysis (WGCNA), machine learning methods, and protein-protein interaction (PPI) network analysis. The expression of hub genes was verified in external datasets and was assessed by receiver operating characteristic (ROC) curve analysis. Gene set enrichment analysis (GSEA) was applied to explore the important functions and pathways of hub genes. The mRNA-microRNA (miRNA)-transcription factors (TFs) regulatory network and potential drugs were predicted based on hub genes. Single-cell RNA sequencing data from the testes of patients with NOA were applied for analyzing the distribution of hub genes in single-cell clusters. Furthermore, testis tissue samples were obtained from patients with NOA and obstructive azoospermia (OA) who underwent testicular biopsy. RT-PCR and Western blot were used to validate hub gene expression. Results: Two immune-related oxidative stress hub genes (SHC1 and FGFR1) were identified. Both hub genes were highly expressed in NOA samples compared to control samples. ROC curve analysis showed a remarkable prediction ability (AUCs > 0.8). GSEA revealed that hub genes were predominantly enriched in toll-like receptor and Wnt signaling pathways. A total of 24 TFs, 82 miRNAs, and 111 potential drugs were predicted based on two hub genes. Single-cell RNA sequencing data in NOA patients indicated that SHC1 and FGFR1 were highly expressed in endothelial cells and Leydig cells, respectively. RT-PCR and Western blot results showed that mRNA and protein levels of both hub genes were significantly upregulated in NOA testis tissue samples, which agree with the findings from analysis of the microarray data. Conclusion: It appears that SHC1 and FGFR1 could be significant immune-related oxidative stress biomarkers for detecting and managing patients with NOA. Our findings provide a novel viewpoint for illustrating potential pathogenesis in men suffering from infertility.

Identification of seven genes essential for male fertility through a genome-wide association study of non-obstructive azoospermia and RNA interference-mediated large-scale functional screening in Drosophila.

Heterozygous Loss-of-function Variants in LHX8 Cause Female Infertility Characterized by Oocyte Maturation Arrest.

Association analysis identifies new risk loci for non-obstructive azoospermia in Chinese men

Excess of Rare Variants in Genes That Are Key Epigenetic Regulators of Spermatogenesis in the Patients with Non-Obstructive Azoospermia.

Identification of Risk Genes in Chinese Nonobstructive Azoospermia Patients Based on Whole-Exome Sequencing

Weighted Correlation Gene Network Analysis Reveals New Potential Mechanisms and Biomarkers in Non-obstructive Azoospermia

CHES-1-like, the Ortholog of a Non-Obstructive Azoospermia-Associated Gene, Blocks Germline Stem Cell Differentiation by Upregulating Dpp Expression Indrosophilatestis

Genetic Insights into Azoospermia and Severe Oligozoospermia: Discovering Seven SNPs through GWAS and In Silico Analysis

Single-Cell Transcriptomics-Based Study of Transcriptional Regulatory Features in the Non-Obstructive Azoospermia Testis

Identification and Functional Analysis of Spermatogenesis-Associated Gene Modules in Azoospermia by Weighted Gene Coexpression Network Analysis

Human obstructive (postvasectomy) and nonobstructive azoospermia – Insights from scRNA-Seq and transcriptome analysis

Whole Exome Sequencing Identifies Genes Associated With Non-Obstructive Azoospermia

Single-cell multi-omics analysis reveals dysfunctional Wnt signaling of spermatogonia in non-obstructive azoospermia

Low-Frequency Germline Variants Across 6p22.2-6p21.33 Are Associated With Non-Obstructive Azoospermia In Han Chinese Men

Deciphering the Molecular Characteristics of Human Idiopathic Nonobstructive Azoospermia from the Perspective of Germ Cells

Drosophila Pif1A is essential for spermatogenesis and is the homolog of human CCDC157, a gene associated with idiopathic NOA

A Genome-Wide Association Study Reveals That Variants Within the HLA Region Are Associated with Risk for Nonobstructive Azoospermia.

Construction and external validation of a 5-gene random forest model to diagnose non-obstructive azoospermia based on the single-cell RNA sequencing of testicular tissue

Defective piRNA Processing and Azoospermia

Identification and validation of SHC1 and FGFR1 as novel immune-related oxidative stress biomarkers of non-obstructive azoospermia

Critical Roles of Long Noncoding RNAs in Drosophila Spermatogenesis.