A comparative cross-platform analysis of cuproptosis-related genes in human nonobstructive azoospermia: An observational study.

This study aimed to identify novel biomarkers associated with cuproptosis in human nonobstructive azoospermia (NOA). We obtained 4 NOA microarray datasets (GSE145467, GSE9210, GSE108886, and GSE45885) from the NCBI Gene Expression Omnibus database and merged them into training set. Another NOA dataset (GSE45887) was used as validation set. Differentially expressed cuproptosis-related genes were identified from training set. Gene Ontology function and Kyoto Encyclopedia of Genes and Genomes pathway analyses were conducted. Least absolute shrinkage and selection operator regression and support vector machine-recursive feature elimination were used to identify hub cuproptosis-related genes. We calculated the expression of the hub cuproptosis-related genes in both validation set and patients with NOA. Gene set variation analysis was used to explore their potential biological functions. The risk prediction model was built by logistic regression analysis and was evaluated in the validation set. Finally, we constructed a competing endogenous RNA network. The training set included 29 patents in the control group and 92 in the NOA group, and 10 cuproptosis-related differentially expressed genes were identified. Subsequently, we screened 6 hub cuproptosis-related genes (DBT, GCSH, NFE2L2, NLRP3, PDHA1, and SLC31A1) by least absolute shrinkage and selection operator regression and support vector machine-recursive feature elimination. GCSH, NFE2L2, NLRP3, and SLC31A1 expressed higher in NOA group than in control group (P < .05) in the validation set (4 patients in control and 16 in NOA groups), while the expression levels of GCSH, NFE2L2, NLRP3, PDHA1, and SLC31A1 were higher in NOA group than in control group (P < .05) in our patients (3 patients in control and 4 in NOA groups). The model based on the 6-gene signature showed superior performance with an AUC value of 0.970 in training set, while 1.0 in validation set. Gene set variation analysis revealed a higher enrichment score of "homologous recombination" in the high expression groups of the 6 hub genes. Finally, we constructed a competing endogenous RNA network and found hsa-miR-335-3p and hsa-miR-1-3p were the most frequently related to the 6 hub genes. DBT, GCSH, NFE2L2, NLRP3, PDHA1, and SLC31A1 may serve as predictors of cuproptosis and play important roles in the NOA pathogenesis.

Jiang S ，Wei Y ，Li Y ，Liu W ，Wang Z ，Meng X ，Zhu Q ，Shen L ... -  《-》

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》

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 ... -  《-》

Integrative bioinformatics approaches for identifying potential biomarkers and pathways involved in non-obstructive azoospermia.

Non-obstructive azoospermia (NOA) is a disease related to spermatogenic disorders. Currently, the specific etiological mechanism of NOA is unclear. This study aimed to use integrated bioinformatics to screen biomarkers and pathways involved in NOA and reveal their potential molecular mechanisms. GSE145467 and GSE108886 gene expression profiles were obtained from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) between NOA tissues and matched obstructive azoospermia (OA) tissues were identified using the GEO2R tool. Common DEGs in the two datasets were screened out by the VennDiagram package. For the functional annotation of common DEGs, DAVID v.6.8 was used to perform Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. In accordance with data collected from the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database, a protein-protein interaction (PPI) network was constructed by Cytoscape. Cytohubba in Cytoscape was used to screen the hub genes. Furthermore, the hub genes were validated based on a separate dataset, GSE9210. Finally, potential micro RNAs (miRNAs) of hub genes were predicted by miRWalk 3.0. A total of 816 common DEGs, including 52 common upregulated and 764 common downregulated genes in two datasets, were screened out. Some of the more important of these pathways, including focal adhesion, PI3K-Akt signaling pathway, cell cycle, oocyte meiosis, AMP-activated protein kinase (AMPK) signaling pathway, FoxO signaling pathway, and Huntington disease, were involved in spermatogenesis. We further identified the top 20 hub genes from the PPI network, including CCNB2, DYNLL2, HMMR, NEK2, KIF15, DLGAP5, NUF2, TTK, PLK4, PTTG1, PBK, CEP55, CDKN3, CDC25C, MCM4, DNAI1, TYMS, PPP2R1B, DNAI2, and DYNLRB2, which were all downregulated genes. In addition, potential miRNAs of hub genes, including hsa-miR-3666, hsa-miR-130b-3p, hsa-miR-15b-5p, hsa-miR-6838-5p, and hsa-miR-195-5p, were screened out. Taken together, the identification of the above hub genes, miRNAs and pathways will help us better understand the mechanisms associated with NOA, and provide potential biomarkers and therapeutic targets for NOA.

Hu T ，Luo S ，Xi Y ，Tu X ，Yang X ，Zhang H ，Feng J ，Wang C ，Zhang Y ... -  《-》

Identification of biomarkers associated with macrophage infiltration in non-obstructive azoospermia using single-cell transcriptomic and microarray data.

Non-obstructive azoospermia (NOA) is a common clinical cause of male infertility. Research suggests that macrophages are linked to testicular function; however, their involvement in NOA remains unknown. To evaluate the importance of macrophages infiltration in NOA and identify the macrophage-related biomarkers, the gene-expression microarray data GSE45885 and the single-cell transcriptomic data GSE149512 were utilized from the Gene Expression Omnibus (GEO). A single-sample gene set enrichment analysis (ssGSEA) was conducted to investigate immune cell proliferation. The Seurat package was used for the single-cell data analysis, and the limma package was used to identify the differentially expressed genes between the NOA and normal samples. Moreover, we conducted a weighted gene co-expression network analysis (WGCNA) to identify the macrophage-related key modules and genes, and conducted Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses for the functional exploration. To identify the macrophage-related biomarkers, we conducted least absolute shrinkage and selection operator (LASSO) and support vector machine-recursive feature elimination (SVM-RFE) analyses. Real-time quantitative polymerase chain reaction (RT-qPCR) was used to verify the marker genes present in NOA. We confirmed that open reading frame 72 gene on chromosome 9 (C9orf72) [area under the curve (AUC) =0.861] and cartilage-associated protein (CRTAP) (AUC =0.917) were the hub genes of NOA, and the RT-qPCR analysis revealed the critical expression of both genes in NOA. Through the combination of tissue transcriptomic and single-cell RNA-sequencing analyses, we concluded that macrophage infiltration is significant in different subtypes of NOA, and we hypothesized that C9orf72 and CRTAP play critical roles in NOA due to their high expression in macrophages.

Luo X ，Zheng H ，Nai Z ，Li M ，Li Y ，Lin N ，Li Y ，Wu Z ... -  《-》