{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,10,13]],"date-time":"2026-10-13T06:22:32Z","timestamp":1791872552676,"version":"4.3.4"},"reference-count":48,"publisher":"Frontiers Media SA","license":[{"start":{"date-parts":[[2022,11,24]],"date-time":"2022-11-24T00:00:00Z","timestamp":1669248000000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":["frontiersin.org"],"crossmark-restriction":true},"short-container-title":["Front. Genet."],"abstract":"<jats:p><jats:bold>Background:<\/jats:bold> 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.<\/jats:p><jats:p><jats:bold>Methods:<\/jats:bold> 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.<\/jats:p><jats:p><jats:bold>Results:<\/jats:bold> 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<jats:sup>+<\/jats:sup> 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.<\/jats:p><jats:p><jats:bold>Conclusion:<\/jats:bold> 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.<\/jats:p>","DOI":"10.3389\/fgene.2022.988047","type":"journal-article","created":{"date-parts":[[2022,11,24]],"date-time":"2022-11-24T09:46:10Z","timestamp":1669283170000},"update-policy":"https:\/\/doi.org\/10.3389\/crossmark-policy","source":"Crossref","is-referenced-by-count":8,"title":["Integrative analyses of potential biomarkers and pathways for non-obstructive azoospermia"],"prefix":"10.3389","volume":"13","author":[{"given":"Yucheng","family":"Zhong","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Xiaoqing","family":"Chen","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Jun","family":"Zhao","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Hao","family":"Deng","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Xiaohang","family":"Li","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Zhongju","family":"Xie","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Bingyu","family":"Zhou","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Zhuojie","family":"Xian","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Xiaoqin","family":"Li","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Guoqun","family":"Luo","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Huan","family":"Li","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1965","published-online":{"date-parts":[[2022,11,24]]},"reference":[{"key":"B1","doi-asserted-by":"publisher","first-page":"1515","DOI":"10.14670\/HH-24.1515","article-title":"Expression of the proto-oncogene c-fos and the immunolocalization of c-fos, phosphorylated c-fos and estrogen receptor beta in the human testis","volume":"24","author":"Ara\u00fajo","year":"2009","journal-title":"Histol. Histopathol."},{"key":"B2","doi-asserted-by":"publisher","first-page":"2079","DOI":"10.1016\/j.molimm.2011.07.001","article-title":"Dietary intake of 17\u03b1-ethinylestradiol promotes leukocytes infiltration in the gonad of the hermaphrodite gilthead seabream","volume":"48","author":"Cabas","year":"2011","journal-title":"Mol. Immunol."},{"key":"B3","doi-asserted-by":"publisher","first-page":"2451","DOI":"10.1093\/hmg\/ddaa101","article-title":"Whole-exome sequencing of a large Chinese azoospermia and severe oligospermia cohort identifies novel infertility causative variants and genes","volume":"29","author":"Chen","year":"2020","journal-title":"Hum. Mol. Genet."},{"key":"B4","doi-asserted-by":"publisher","first-page":"187","DOI":"10.1007\/978-1-0716-1839-4_12","article-title":"geneExpressionFromGEO: An R package to facilitate data reading from gene expression Omnibus (GEO)","volume":"2401","author":"Chicco","year":"2022","journal-title":"Methods Mol. Biol."},{"key":"B5","doi-asserted-by":"publisher","first-page":"265","DOI":"10.1093\/humupd\/dmx008","article-title":"Sperm recovery and ICSI outcomes in klinefelter syndrome: A systematic review and meta-analysis","volume":"23","author":"Corona","year":"2017","journal-title":"Hum. Reprod. Update"},{"key":"B6","doi-asserted-by":"publisher","first-page":"3117","DOI":"10.3390\/ani11113117","article-title":"The distribution, expression patterns and functional analysis of NR1D1 and NR4A2 in the reproductive axis tissues of the male Tianzhu White Yak","volume":"11","author":"Dai","year":"2021","journal-title":"Animals."},{"key":"B7","doi-asserted-by":"publisher","first-page":"167149","DOI":"10.1016\/j.jmb.2021.167149","article-title":"Protein interaction network-based deep learning framework for identifying disease-associated human proteins","volume":"433","author":"Das","year":"2021","journal-title":"J. Mol. Biol."},{"key":"B8","doi-asserted-by":"publisher","first-page":"561","DOI":"10.1007\/s00418-021-02023-7","article-title":"Impaired macrophages and failure of steroidogenesis and spermatogenesis in rat testes with cytokines deficiency induced by diacerein","volume":"156","author":"de Oliveira","year":"2021","journal-title":"Histochem. Cell Biol."},{"key":"B9","doi-asserted-by":"publisher","first-page":"617133","DOI":"10.3389\/fgene.2021.617133","article-title":"Weighted correlation gene network analysis reveals new potential mechanisms and biomarkers in non-obstructive azoospermia","volume":"12","author":"Dong","year":"2021","journal-title":"Front. Genet."},{"key":"B10","doi-asserted-by":"publisher","first-page":"470","DOI":"10.1016\/s0015-0282(16)59421-7","article-title":"Immunocompetent cells in human testis in health and disease","volume":"48","author":"el-Demiry","year":"1987","journal-title":"Fertil. Steril."},{"key":"B11","doi-asserted-by":"publisher","first-page":"249","DOI":"10.1159\/000316604","article-title":"Klinefelter's syndrome: A clinical and therapeutical update","volume":"4","author":"Forti","year":"2010","journal-title":"Sex. Dev."},{"key":"B12","doi-asserted-by":"publisher","first-page":"2628","DOI":"10.1093\/bioinformatics\/btz931","article-title":"ShinyGO: A graphical gene-set enrichment tool for animals and plants","volume":"36","author":"Ge","year":"2020","journal-title":"Bioinformatics"},{"key":"B13","doi-asserted-by":"publisher","first-page":"57","DOI":"10.3390\/ijms22010057","article-title":"T lymphocytes and testicular immunity: A new insight into immune regulation in testes","volume":"22","author":"Gong","year":"2020","journal-title":"Int. J. Mol. Sci."},{"key":"B14","doi-asserted-by":"publisher","first-page":"3844","DOI":"10.1093\/bioinformatics\/btac409","article-title":"ggtranscript: an R package for the visualization and interpretation of transcript isoforms using ggplot2","volume":"38","author":"Gustavsson","year":"2022","journal-title":"Bioinformatics"},{"key":"B15","doi-asserted-by":"publisher","first-page":"130","DOI":"10.1016\/j.jri.2011.01.010","article-title":"Toll-like receptors and signalling in spermatogenesis and testicular responses to inflammation--a perspective","volume":"88","author":"Hedger","year":"2011","journal-title":"J. Reprod. Immunol."},{"key":"B16","doi-asserted-by":"publisher","first-page":"91","DOI":"10.1186\/s12958-018-0401-7","article-title":"Spermatogenesis improved by suppressing the high level of endogenous gonadotropins in idiopathic non-obstructive azoospermia: A case control pilot study","volume":"16","author":"Hu","year":"2018","journal-title":"Reprod. Biol. Endocrinol."},{"key":"B17","doi-asserted-by":"publisher","first-page":"W367","DOI":"10.1093\/nar\/gkac364","article-title":"GEOexplorer: A webserver for gene expression analysis and visualisation","volume":"50","author":"Hunt","year":"2022","journal-title":"Nucleic Acids Res."},{"key":"B18","doi-asserted-by":"publisher","first-page":"62","DOI":"10.1016\/s0022-5347(17)38662-7","article-title":"Evaluation of the azoospermic patient","volume":"142","author":"Jarow","year":"1989","journal-title":"J. Urol."},{"key":"B19","doi-asserted-by":"publisher","first-page":"3818","DOI":"10.1038\/s41598-022-07624-3","article-title":"GOntoSim: A semantic similarity measure based on LCA and common descendants","volume":"2","author":"Kamran","year":"2022","journal-title":"Sci. Rep."},{"key":"B20","doi-asserted-by":"publisher","first-page":"D545","DOI":"10.1093\/nar\/gkaa970","article-title":"Kegg: Integrating viruses and cellular organisms","volume":"49","author":"Kanehisa","year":"2021","journal-title":"Nucleic Acids Res."},{"key":"B21","doi-asserted-by":"publisher","first-page":"8208","DOI":"10.1158\/0008-5472.CAN-04-2134","article-title":"Nuclear hormone receptor NR4A2 is involved in cell transformation and apoptosis","volume":"64","author":"Ke","year":"2004","journal-title":"Cancer Res."},{"key":"B22","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1002\/mrd.1124","article-title":"Keratins: Unraveling the coordinated construction of scaffolds in spermatogenic cells","volume":"61","author":"Kierszenbaum","year":"2002","journal-title":"Mol. Reprod. Dev."},{"key":"B23","doi-asserted-by":"publisher","first-page":"271","DOI":"10.1016\/j.beem.2010.08.006","article-title":"Male infertility: Pathogenesis and clinical diagnosis","volume":"25","author":"Krausz","year":"2011","journal-title":"Best. Pract. Res. Clin. Endocrinol. Metab."},{"key":"B24","doi-asserted-by":"publisher","first-page":"bbaa219","DOI":"10.1093\/bib\/bbaa219","article-title":"A review of digital cytometry methods: Estimating the relative abundance of cell types in a bulk of cells","volume":"22","author":"Le","year":"2021","journal-title":"Brief. Bioinform."},{"key":"B25","doi-asserted-by":"publisher","first-page":"417","DOI":"10.1016\/j.cels.2015.12.004","article-title":"The Molecular Signatures Database (MSigDB) hallmark gene set collection","volume":"1","author":"Liberzon","year":"2015","journal-title":"Cell Syst."},{"key":"B26","doi-asserted-by":"publisher","first-page":"175","DOI":"10.3791\/62528","article-title":"Three differential expression analysis methods for RNA sequencing: Limma, EdgeR, DESeq2","volume":"2021","author":"Liu","year":"2021","journal-title":"J. Vis. Exp."},{"key":"B27","doi-asserted-by":"publisher","first-page":"30","DOI":"10.1111\/j.2047-2927.2013.00156.x","article-title":"Seminal, clinical and colour-Doppler ultrasound correlations of prostatitis-like symptoms in males of infertile couples","volume":"2","author":"Lotti","year":"2014","journal-title":"Andrology"},{"key":"B28","doi-asserted-by":"publisher","first-page":"2021","DOI":"10.1210\/me.2007-0370","article-title":"The orphan nuclear receptor NUR77 regulates hormone-induced StAR transcription in Leydig cells through cooperation with Ca2+\/calmodulin-dependent protein kinase I","volume":"22","author":"Martin","year":"2008","journal-title":"Mol. Endocrinol."},{"key":"B29","doi-asserted-by":"publisher","first-page":"3","DOI":"10.1095\/biolreprod.110.083824","article-title":"Nuclear receptors in Leydig cell gene expression and function","volume":"83","author":"Martin","year":"2010","journal-title":"Biol. Reprod."},{"key":"B30","doi-asserted-by":"publisher","first-page":"10","DOI":"10.1016\/j.semcdb.2016.04.009","article-title":"Spermatogenesis in humans and its affecting factors","volume":"59","author":"Neto","year":"2016","journal-title":"Semin. Cell Dev. Biol."},{"key":"B31","doi-asserted-by":"publisher","first-page":"725276","DOI":"10.3389\/fdata.2021.725276","article-title":"Statistical enrichment analysis of samples: A general-purpose tool to annotate metadata neighborhoods of biological samples","volume":"4","author":"Nguyen","year":"2021","journal-title":"Front. Big Data"},{"key":"B32","doi-asserted-by":"publisher","first-page":"W405","DOI":"10.1093\/nar\/gkac480","article-title":"HemI 2.0: An online service for heatmap illustration","volume":"50","author":"Ning","year":"2022","journal-title":"Nucleic Acids Res."},{"key":"B33","doi-asserted-by":"publisher","first-page":"214","DOI":"10.1126\/science.abg3067","article-title":"Olfactory receptor 2 in vascular macrophages drives atherosclerosis by NLRP3-dependent IL-1 production","volume":"375","author":"Orecchioni","year":"2022","journal-title":"Science"},{"key":"B34","doi-asserted-by":"publisher","first-page":"337","DOI":"10.1016\/j.mcna.2017.10.008","article-title":"Male infertility diagnosis and treatment in the era of in vitro fertilization and intracytoplasmic sperm injection","volume":"102","author":"Pan","year":"2018","journal-title":"Med. Clin. North Am."},{"key":"B35","doi-asserted-by":"publisher","first-page":"101479","DOI":"10.1016\/j.beem.2020.101479","article-title":"Genetic mutations contributing to non-obstructive azoospermia","volume":"34","author":"Pe\u00f1a","year":"2020","journal-title":"Best. Pract. Res. Clin. Endocrinol. Metab."},{"key":"B36","doi-asserted-by":"publisher","first-page":"102","DOI":"10.1002\/mrd.1136","article-title":"Molecular cloning of Odf3 encoding a novel coiled-coil protein of sperm tail outer dense fibers","volume":"61","author":"Petersen","year":"2002","journal-title":"Mol. Reprod. Dev."},{"key":"B37","first-page":"336","article-title":"Genetics of human male infertility","volume":"50","author":"Poongothai","year":"2009","journal-title":"Singap. Med. J."},{"key":"B38","doi-asserted-by":"publisher","first-page":"e106","DOI":"10.1002\/cpbi.106","article-title":"Network building with the Cytoscape BioGateway app explained in five use cases","volume":"72","author":"Puig","year":"2020","journal-title":"Curr. Protoc. Bioinforma."},{"key":"B39","doi-asserted-by":"publisher","first-page":"48","DOI":"10.1016\/j.jsbmb.2015.04.016","article-title":"Nuclear receptor 4A (NR4A) family - orphans no more","volume":"157","author":"Safe","year":"2016","journal-title":"J. Steroid Biochem. Mol. Biol."},{"key":"B40","doi-asserted-by":"publisher","first-page":"857","DOI":"10.1007\/s43032-021-00837-3","article-title":"Histone methylation regulates gene expression in the round spermatids to set the RNA payloads of sperm","volume":"29","author":"Sarkar","year":"2022","journal-title":"Reprod. Sci."},{"key":"B41","doi-asserted-by":"publisher","first-page":"D605","DOI":"10.1093\/nar\/gkaa1074","article-title":"The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene\/measurement sets","volume":"49","author":"Szklarczyk","year":"2021","journal-title":"Nucleic Acids Res."},{"key":"B42","doi-asserted-by":"publisher","first-page":"554","DOI":"10.1016\/S2213-8587(16)30043-2","article-title":"Concepts in diagnosis and therapy for male reproductive impairment","volume":"5","author":"Tournaye","year":"2017","journal-title":"Lancet. Diabetes Endocrinol."},{"key":"B43","doi-asserted-by":"publisher","first-page":"602","DOI":"10.1016\/j.steroids.2008.11.017","article-title":"Molecular mechanisms of testosterone action in spermatogenesis","volume":"74","author":"Walker","year":"2009","journal-title":"Steroids"},{"key":"B44","doi-asserted-by":"publisher","first-page":"599","DOI":"10.1016\/j.stem.2018.08.007","article-title":"Single-cell RNA sequencing analysis reveals sequential cell fate transition during human spermatogenesis","volume":"23","author":"Wang","year":"2018","journal-title":"Cell Stem Cell"},{"key":"B45","doi-asserted-by":"publisher","first-page":"100141","DOI":"10.1016\/j.xinn.2021.100141","article-title":"clusterProfiler 4.0: A universal enrichment tool for interpreting omics data","volume":"2","author":"Wu","year":"2021","journal-title":"Innovation."},{"key":"B46","doi-asserted-by":"publisher","first-page":"5683","DOI":"10.1038\/s41467-020-19414-4","article-title":"Single-cell analysis of developing and azoospermia human testicles reveals central role of Sertoli cells","volume":"11","author":"Zhao","year":"2020","journal-title":"Nat. Commun."},{"key":"B47","doi-asserted-by":"publisher","first-page":"458","DOI":"10.1186\/s12864-021-07793-x","article-title":"Online database for brain cancer-implicated genes: Exploring the subtype-specific mechanisms of brain cancer","volume":"22","author":"Zhao","year":"2021","journal-title":"BMC Genomics"},{"key":"B48","doi-asserted-by":"publisher","first-page":"3934","DOI":"10.1002\/jcb.27677","article-title":"Identification and functional analysis of spermatogenesis-associated gene modules in azoospermia by weighted gene coexpression network analysis","volume":"120","author":"Zheng","year":"2019","journal-title":"J. Cell. Biochem."}],"container-title":["Frontiers in Genetics"],"original-title":[],"link":[{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fgene.2022.988047\/full","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2022,11,24]],"date-time":"2022-11-24T09:46:14Z","timestamp":1669283174000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.frontiersin.org\/articles\/10.3389\/fgene.2022.988047\/full"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2022,11,24]]},"references-count":48,"alternative-id":["10.3389\/fgene.2022.988047"],"URL":"https:\/\/doi.org\/10.3389\/fgene.2022.988047","relation":{},"ISSN":["1664-8021"],"issn-type":[{"value":"1664-8021","type":"electronic"}],"subject":[],"published":{"date-parts":[[2022,11,24]]},"article-number":"988047"}}