Identifying potential biomarkers for non-obstructive azoospermia using WGCNA and machine learning algorithms.
The cause and mechanism of non-obstructive azoospermia (NOA) is complicated; therefore, an effective therapy strategy is yet to be developed. This study aimed to analyse the pathogenesis of NOA at the molecular biological level and to identify the core regulatory genes, which could be utilised as potential biomarkers.
Three NOA microarray datasets (GSE45885, GSE108886, and GSE145467) were collected from the GEO database and merged into training sets; a further dataset (GSE45887) was then defined as the validation set. Differential gene analysis, consensus cluster analysis, and WGCNA were used to identify preliminary signature genes; then, enrichment analysis was applied to these previously screened signature genes. Next, 4 machine learning algorithms (RF, SVM, GLM, and XGB) were used to detect potential biomarkers that are most closely associated with NOA. Finally, a diagnostic model was constructed from these potential biomarkers and visualised as a nomogram. The differential expression and predictive reliability of the biomarkers were confirmed using the validation set. Furthermore, the competing endogenous RNA network was constructed to identify the regulatory mechanisms of potential biomarkers; further, the CIBERSORT algorithm was used to calculate immune infiltration status among the samples.
A total of 215 differentially expressed genes (DEGs) were identified between NOA and control groups (27 upregulated and 188 downregulated genes). The WGCNA results identified 1123 genes in the MEblue module as target genes that are highly correlated with NOA positivity. The NOA samples were divided into 2 clusters using consensus clustering; further, 1027 genes in the MEblue module, which were screened by WGCNA, were considered to be target genes that are highly correlated with NOA classification. The 129 overlapping genes were then established as signature genes. The XGB algorithm that had the maximum AUC value (AUC=0.946) and the minimum residual value was used to further screen the signature genes. IL20RB, C9orf117, HILS1, PAOX, and DZIP1 were identified as potential NOA biomarkers. This 5 biomarker model had the highest AUC value, of up to 0.982, compared to other single biomarker models; additionally, the results of this biomarker model were verified in the validation set.
As IL20RB, C9orf117, HILS1, PAOX, and DZIP1 have been determined to possess the strongest association with NOA, these five genes could be used as potential therapeutic targets for NOA patients. Furthermore, the model constructed using these five genes, which possessed the highest diagnostic accuracy, may be an effective biomarker model that warrants further experimental validation.
Tang Q
,Su Q
,Wei L
,Wang K
,Jiang T
... -
《Frontiers in Endocrinology》
Identification and validation of SHC1 and FGFR1 as novel immune-related oxidative stress biomarkers of non-obstructive azoospermia.
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.
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.
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.
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.
Pan Y
,Chen X
,Zhou H
,Xu M
,Li Y
,Wang Q
,Xu Z
,Ren C
,Liu L
,Liu X
... -
《Frontiers in Endocrinology》
Evaluation of immune status in testis and macrophage polarization associated with testicular damage in patients with nonobstructive azoospermia.
Immune cells residing in the testicular interstitial space form the immunological microenvironment of the testis. They are assumed to play a role in maintaining testicular homeostasis and immune privilege. However, the immune status and related cell polarization in patients with nonobstructive azoospermia (NOA) remains poorly characterized. System evaluation of the testis immunological microenvironment in NOA patients may help to reveal the mechanisms of idiopathic azoospermia.
The gene expression patterns of immune cells in normal human testes were systematically analyzed by single-cell RNA sequencing (scRNA-seq) and preliminarily verification by the human protein atlas (HPA) online database. The immune cell infiltration profiles and immune status of patients with NOA was analyzed by single-sample gene set enrichment analysis (ssGSEA) and gene set variation analysis (GSVA) based on four independent public microarray datasets (GSE45885, GSE45887, GSE9210, and GSE145467), obtained from Gene Expression Omnibus (GEO) online database. The relationship between immune cells and spermatogenesis score was further analyzed by Spearman correlation analysis. Finally, immunohistochemistry (IHC) staining was performed to identify the main immune cell types and their polarization status in patients with NOA.
Both scRNA-seq and HPA analysis showed that testicular macrophages represent the largest pool of immune cells in the normal testis, and also exhibit an attenuated inflammatory response by expressing high levels of tolerance proteins (CD163, IL-10, TGF-β, and VEGF) and reduced expression of TLR signaling pathway-related genes. Correlation analysis revealed that the testicular immune score and macrophages including M1 and M2 macrophages were significantly negatively correlated with spermatogenesis score in patients with NOA (GSE45885 and GSE45887). In addition, the number of M1 and M2 macrophages was significantly higher in patients with NOA (GSE9210 and GSE145467) than in normal testis. GSVA analysis indicated that the immunological microenvironment in NOA tissues was manifested by activated immune system and pro-inflammatory status. IHC staining results showed that the number of M1 and M2 macrophages was significantly higher in NOA tissues than in normal testis and negatively correlated with the Johnson score.
Testicular macrophage polarization may play a vital role in NOA development and is a promising potential therapeutic target.
Zheng W
,Zhang S
,Jiang S
,Huang Z
,Chen X
,Guo H
,Li M
,Zheng S
... -
《-》