Integrative analyses of potential biomarkers and pathways for non-obstructive azoospermia.

Background: Non-obstructive azoospermia (NOA) is the most severe form of male infertility. Currently, the molecular mechanisms underlying NOA pathology have not yet been elucidated. Hence, elucidation of the mechanisms of NOA and exploration of potential biomarkers are essential for accurate diagnosis and treatment of this disease. In the present study, we aimed to screen for biomarkers and pathways involved in NOA and reveal their potential molecular mechanisms using integrated bioinformatics. Methods: We downloaded two gene expression datasets from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) in NOA and matched the control group tissues were identified using the limma package in R software. Subsequently, Gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), gene set enrichment analysis (GSEA), protein-protein interaction (PPI) network, gene-microRNAs network, and transcription factor (TF)-hub genes regulatory network analyses were performed to identify hub genes and associated pathways. Finally, we conducted immune infiltration analysis using CIBERSORT to evaluate the relationship between the hub genes and the NOA immune infiltration levels. Results: We identified 698 common DEGs, including 87 commonly upregulated and 611 commonly downregulated genes in the two datasets. GO analysis indicated that the most significantly enriched gene was protein polyglycylation, and KEGG pathway analysis revealed that the DEGs were most significantly enriched in taste transduction and pancreatic secretion signaling pathways. GSEA showed that DEGs affected the biological functions of the ribosome, focaladhesion, and protein_expor. We further identified the top 31 hub genes from the PPI network, and friends analysis of hub genes in the PPI network showed that NR4A2 had the highest score. In addition, immune infiltration analysis found that CD8+ T cells and plasma cells were significantly correlated with ODF3 expression, whereas naive B cells, plasma cells, monocytes, M2 macrophages, and resting mast cells showed significant variation in the NR4A2 gene expression group, and there were differences in T cell regulatory immune cell infiltration in the FOS gene expression groups. Conclusion: The present study successfully constructed a regulatory network of DEGs between NOA and normal controls and screened three hub genes using integrative bioinformatics analysis. In addition, our results suggest that functional changes in several immune cells in the immune microenvironment may play an important role in spermatogenesis. Our results provide a novel understanding of the molecular mechanisms of NOA and offer potential biomarkers for its diagnosis and treatment.

Zhong Y ，Chen X ，Zhao J ，Deng H ，Li X ，Xie Z ，Zhou B ，Xian Z ，Li X ，Luo G ，Li H ... -  《Frontiers in Genetics》

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

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

Alagundagi DB ，Ghate SD ，Shetty P ，Gollapalli P ，Shetty P ，Patil P ... -  《-》

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》

Six potential biomarkers in septic shock: a deep bioinformatics and prospective observational study.

Septic shock occurs when sepsis is related to severe hypotension and leads to a remarkable high number of deaths. The early diagnosis of septic shock is essential to reduce mortality. High-quality biomarkers can be objectively measured and evaluated as indicators to accurately predict disease diagnosis. However, single-gene prediction efficiency is inadequate; therefore, we identified a risk-score model based on gene signature to elevate predictive efficiency. The gene expression profiles of GSE33118 and GSE26440 were downloaded from the Gene Expression Omnibus (GEO) database. These two datasets were merged, and the differentially expressed genes (DEGs) were identified using the limma package in R software. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichments of DEGs were performed. Subsequently, Lasso regression and Boruta feature selection algorithm were combined to identify the hub genes of septic shock. GSE9692 was then subjected to weighted gene co-expression network analysis (WGCNA) to identify the septic shock-related gene modules. Subsequently, the genes within such modules that matched with septic shock-related DEGs were identified as the hub genes of septic shock. To further understand the function and signaling pathways of hub genes, we performed gene set variation analysis (GSVA) and then used the CIBERSORT tool to analyze the immune cell infiltration pattern of diseases. The diagnostic value of hub genes in septic shock was determined using receiver operating characteristic (ROC) analysis and verified using quantitative PCR (qPCR) and Western blotting in our hospital patients with septic shock. A total of 975 DEGs in the GSE33118 and GSE26440 databases were obtained, of which 30 DEGs were remarkably upregulated. With the use of Lasso regression and Boruta feature selection algorithm, six hub genes (CD177, CLEC5A, CYSTM1, MCEMP1, MMP8, and RGL4) with expression differences in septic shock were screened as potential diagnostic markers for septic shock among the significant DEGs and were further validated in the GSE9692 dataset. WGCNA was used to identify the co-expression modules and module-trait correlation. Enrichment analysis showed significant enrichment in the reactive oxygen species pathway, hypoxia, phosphatidylinositol 3-kinases (PI3K)/Protein Kinase B (AKT)/mammalian target of rapamycin (mTOR) signaling, nuclear factor-κβ/tumor necrosis factor alpha (NF-κβ/TNF-α), and interleukin-6 (IL-6)/Janus Kinase (JAK)/Signal Transducers and Activators of Transcription 3 (STAT3) signaling pathways. The receiver operating characteristic curve (ROC) of these signature genes was 0.938, 0.914, 0.939, 0.956, 0.932, and 0.914, respectively. In the immune cell infiltration analysis, the infiltration of M0 macrophages, activated mast cells, neutrophils, CD8 T cells, and naive B cells was more significant in the septic shock group. In addition, higher expression levels of CD177, CLEC5A, CYSTM1, MCEMP1, MMP8, and RGL4 messenger RNA (mRNA) were observed in peripheral blood mononuclear cells (PBMCs) isolated from septic shock patients than from healthy donors. Higher expression levels of CD177 and MMP8 proteins were also observed in the PBMCs isolated from septic shock patients than from control participants. CD177, CLEC5A, CYSTM1, MCEMP1, MMP8, and RGL4 were identified as hub genes, which were of considerable value in the early diagnosis of septic shock patients. These preliminary findings are of great significance for studying immune cell infiltration in the pathogenesis of septic shock, which should be further validated in clinical studies and basic studies.

Kong C ，Zhu Y ，Xie X ，Wu J ，Qian M ... -  《Frontiers in Immunology》