Identifying potential biomarkers for non-obstructive azoospermia using WGCNA and machine learning algorithms.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:37854191", "kind": "science-record", "revision": 4, "content_hash": "sha256:b1ba362119a0d3daa09e0f1612871940dac36adf66d92d1a8094accd9b9204c3", "created_at": "2026-11-01T19:02:42Z", "updated_at": "2026-11-01T22:52:54Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Identifying potential biomarkers for non-obstructive azoospermia using WGCNA and machine learning algorithms.", "doi": "10.3389/fendo.2023.1108616", "pmid": "37854191", "authors": ["Qizhen Tang", "Quanxin Su", "Letian Wei", "Kenan Wang", "Tao Jiang"], "journal": "Frontiers in endocrinology", "year": 2023, "abstract": "Objective: 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.\nMethods: 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.\nResults: 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.\nConclusions: 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.", "claims": [{"id": "claim-method", "statement": "Finally, we defined the 5 most significant genes (IL20RB, C9orf117, HILS1, PAOX, and DZIP1), which were identified, by XGB, as potential biomarkers", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:37854191", "locator": "sections[17] \"Identification of potential biomarkers\" paragraph 2", "quote": "Finally, we defined the 5 most significant genes (IL20RB, C9orf117, HILS1, PAOX, and DZIP1), which were identified, by XGB, as potential biomarkers"}], "state": "complete"}, {"id": "claim-data-access", "statement": "The original contributions presented in the study are included in the article/ Supplementary Material . Further inquiries can be directed to the corresponding author.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:37854191", "locator": "sections[24] \"Data availability statement\" paragraph 2", "quote": "The original contributions presented in the study are included in the article/ Supplementary Material . Further inquiries can be directed to the corresponding author."}], "state": "complete"}, {"id": "claim-target-result", "statement": "IL20RB, C9orf117, HILS1, PAOX, and DZIP1 were identified as potential NOA biomarkers.", "kind": "qualitative", "target": null, "evidence": [{"asset": "source:pmid:37854191", "locator": "sections[0] \"Abstract\" paragraph 3", "quote": "IL20RB, C9orf117, HILS1, PAOX, and DZIP1 were identified as potential NOA biomarkers."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "complete", "objective": "IL20RB, C9orf117, HILS1, PAOX, and DZIP1 were identified as potential NOA biomarkers.", "steps": [{"id": "read_published_estimate", "order": 1, "name": "Read the published estimate from the source", "kind": "other", "state": "complete", "blocked_on": null, "description": "The claim, method and evidence span published with the source article.", "duration_hours": null, "parameters": {}, "inputs": ["source:pmid:37854191"], "outputs": ["the published estimate and its evidence"], "instruments": [], "reagents": []}, {"id": "obtain_source_data", "order": 2, "name": "Obtain the source data from the named repository", "kind": "data-preparation", "state": "complete", "blocked_on": null, "description": "The GEO series accession resolves and its sample accessions are listed. This lane did not download the dataset or read its licence.", "duration_hours": null, "parameters": {}, "inputs": ["data:geo:GSE108886"], "outputs": ["the source data with its repository provenance"], "instruments": [], "reagents": []}, {"id": "run_analysis", "order": 3, "name": "Run the published analysis", "kind": "analysis", "state": "complete", "blocked_on": null, "description": "Apply the paper's statistical method to the resolved series and compare the result with the published target.", "duration_hours": null, "parameters": {}, "inputs": ["data:geo:GSE108886"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": []}, "data": {"state": "partial", "assets": [{"id": "source:pmid:37854191", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.3389/fendo.2023.1108616", "sha256": "sha256:f341f6c428b0d6ed16ba94829c757bcf3457fe2bf80793dc05f55dcab7bcd828", "bytes": 33277, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.3389/fendo.2023.1108616", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC10579891", "source": "europe_pmc_fulltext", "coverage": "whole_body_without_references", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE108886", "role": "raw", "format": "geo-series", "media_type": "application/vnd.ncbi.geo-series+json", "uri": "https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE108886", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE108886", "version": null, "note": "The original contributions presented in the study are included in the article/ Supplementary Material . Further inquiries can be directed to the corresponding author.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The original contributions presented in the study are included in the article/ Supplementary Material . Further inquiries can be directed to the corresponding author.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM2915444", "GSM2915445", "GSM2915446", "GSM2915447", "GSM2915448", "GSM2915449", "GSM2915450", "GSM2915451", "GSM2915452", "GSM2915453", "GSM2915454", "GSM2915455"], "sample_count": 12, "group_sizes": null, "group_sizes_basis": null, "platform": "10558", "taxon": "Homo sapiens", "title": "Spermatogenomics: correlating testicular gene expression to human male infertility", "resolution_checked_at": "2026-11-01T15:39:54Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}, {"id": "data:geo:GSE145467", "role": "raw", "format": "geo-series", "media_type": "application/vnd.ncbi.geo-series+json", "uri": "https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE145467", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE145467", "version": null, "note": "The original contributions presented in the study are included in the article/ Supplementary Material . Further inquiries can be directed to the corresponding author.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The original contributions presented in the study are included in the article/ Supplementary Material . Further inquiries can be directed to the corresponding author.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM4318655", "GSM4318656", "GSM4318657", "GSM4318658", "GSM4318660", "GSM4318662", "GSM4318665", "GSM4318668", "GSM4318670", "GSM4318672", "GSM4318673", "GSM4318674", "GSM4318675", "GSM4318676", "GSM4318677", "GSM4318678", "GSM4318679", "GSM4318680", "GSM4318681", "GSM4318682"], "sample_count": 20, "group_sizes": null, "group_sizes_basis": null, "platform": "4133", "taxon": "Homo sapiens", "title": "Transcriptome changes in patients with severely impaired spermatogenesis", "resolution_checked_at": "2026-11-01T15:39:54Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}, {"id": "data:geo:GSE45885", "role": "raw", "format": "geo-series", "media_type": "application/vnd.ncbi.geo-series+json", "uri": "https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE45885", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE45885", "version": null, "note": "The original contributions presented in the study are included in the article/ Supplementary Material . Further inquiries can be directed to the corresponding author.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The original contributions presented in the study are included in the article/ Supplementary Material . Further inquiries can be directed to the corresponding author.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM1118241", "GSM1118242", "GSM1118243", "GSM1118244", "GSM1118245", "GSM1118246", "GSM1118247", "GSM1118248", "GSM1118249", "GSM1118250", "GSM1118251", "GSM1118252", "GSM1118253", "GSM1118254", "GSM1118255", "GSM1118256", "GSM1118257", "GSM1118258", "GSM1118259", "GSM1118260", "GSM1118261", "GSM1118262", "GSM1118263", "GSM1118264", "GSM1118265", "GSM1118266", "GSM1118267", "GSM1118268", "GSM1118269", "GSM1118270", "GSM1118271"], "sample_count": 31, "group_sizes": null, "group_sizes_basis": null, "platform": "6244", "taxon": "Homo sapiens", "title": "Potential biomarkers of non-obstructive azoospermia identified in microarray gene expression analysis", "resolution_checked_at": "2026-11-01T15:39:54Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}, {"id": "data:geo:GSE45887", "role": "raw", "format": "geo-series", "media_type": "application/vnd.ncbi.geo-series+json", "uri": "https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE45887", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE45887", "version": null, "note": "The original contributions presented in the study are included in the article/ Supplementary Material . Further inquiries can be directed to the corresponding author.", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The original contributions presented in the study are included in the article/ Supplementary Material . Further inquiries can be directed to the corresponding author.", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM1118281", "GSM1118282", "GSM1118283", "GSM1118284", "GSM1118285", "GSM1118286", "GSM1118287", "GSM1118288", "GSM1118289", "GSM1118290", "GSM1118291", "GSM1118292", "GSM1118293", "GSM1118294", "GSM1118295", "GSM1118296", "GSM1118297", "GSM1118298", "GSM1118299", "GSM1118300"], "sample_count": 20, "group_sizes": null, "group_sizes_basis": null, "platform": "6244", "taxon": "Homo sapiens", "title": "The gene expression analysis of paracrine/autocrine factors in patients with spermatogenetic failure compared to normal spermatogenesis", "resolution_checked_at": "2026-11-01T18:06:44Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}], "gaps": ["asset:data:geo:GSE108886 not downloaded", "asset:data:geo:GSE108886 licence not read", "asset:data:geo:GSE145467 not downloaded", "asset:data:geo:GSE145467 licence not read", "asset:data:geo:GSE45885 not downloaded", "asset:data:geo:GSE45885 licence not read", "asset:data:geo:GSE45887 not downloaded", "asset:data:geo:GSE45887 licence not read"]}, "analysis": {"state": "complete", "methods": [{"id": "method-1", "family": "machine_learning", "specific": "gradient_boosting", "label": "Machine learning", "dsl_operations": ["train_classifier", "cross_validate", "roc_auc"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["labelled_feature_matrix"], "evidence": {"score": 1.0, "band": "verbatim", "quote": "Finally, we defined the 5 most significant genes (IL20RB, C9orf117, HILS1, PAOX, and DZIP1), which were identified, by XGB, as potential biomarkers"}}, {"id": "method-2", "family": "machine_learning", "specific": "random_forest", "label": "Machine learning", "dsl_operations": [], "missing_operations": ["machine_learning_random_forest"], "expressible": false, "parameters": {}, "data_shapes": ["labelled_feature_matrix"], "evidence": {"score": 0.8, "band": "high", "quote": "Random Forest"}}, {"id": "method-3", "family": "machine_learning", "specific": "support_vector_machine", "label": "Machine learning", "dsl_operations": ["train_classifier", "cross_validate", "roc_auc"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["labelled_feature_matrix"], "evidence": {"score": 0.7, "band": "medium", "quote": "SVM"}}], "targets": [], "gaps": []}, "capability": {"state": "complete", "requirements": [{"capability_id": "capability:computational_analysis", "label": "Computational analysis platform", "assay_families": ["meta_analysis", "machine_learning", "survival_analysis", "statistical_regression", "pca", "molecular_docking", "network_analysis", "differential_expression", "gwas"], "availability_state": "owned", "tier_unlocked": "R1", "candidate": "Automated Assays Laravel application and simulated adapters"}], "gaps": []}, "readiness": {"state": "complete", "tier": "R1", "blocking_axis": "none", "blocked_on": "nothing:replicable_now", "axes": {"method": {"tier": "R1", "expressible": true, "required_operations": ["train_classifier", "cross_validate", "roc_auc"], "missing_operations": [], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: gradient_boosting."}, "capability": {"tier": "R1", "capability_id": "capability:computational_analysis", "availability_state": "owned", "basis": "The computational analysis is owned. The physical bench the paper used is not, so a physical validation cannot be run here."}, "data": {"tier": "R1", "reachable": true, "accessions": ["GSE45885", "GSE45887", "GSE108886", "GSE145467"], "resolved_accessions": ["GSE45885", "GSE45887", "GSE108886", "GSE145467"], "statements": ["The original contributions presented in the study are included in the article/ Supplementary Material . Further inquiries can be directed to the corresponding author."], "basis": "The paper names GSE45885, GSE45887, GSE108886, GSE145467 in GEO and the accession resolves in the repository. No dataset was downloaded and no licence was read."}}, "can_run": {"dsl_operations": ["train_classifier", "cross_validate", "roc_auc"], "assets": ["source:pmid:37854191", "data:geo:GSE108886"], "steps": ["run_analysis"]}, "must_build": {"dsl_operations": [], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE108886"]}, "gaps": []}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-11-01T18:56:26Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 9326, "completion_tokens": 1609, "cached_tokens": 256}}, "human_critique": null, "stages": [{"id": "extraction", "state": "complete", "reader": "MiniMax-M3.1-Flash", "provider": "minimax"}], "gaps": ["review.independent_reader", "review.machine_critique", "review.human_critique"]}, "execution": {"state": "absent", "status": "not_run", "mode": null, "runs": [], "gaps": ["execution.not_run"]}, "provenance": {"entities": [{"id": "pmid:37854191", "type": "paper", "label": "pmid:37854191"}, {"id": "source:pmid:37854191", "type": "manuscript", "label": "PMID 37854191 full text"}, {"id": "data:geo:GSE108886", "type": "dataset", "label": "GEO GSE108886"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:37854191", "type": "paper", "label": "pmid:37854191"}, {"id": "source:pmid:37854191", "type": "manuscript", "label": "PMID 37854191 full text"}, {"id": "data:geo:GSE108886", "type": "dataset", "label": "GEO GSE108886"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-37854191", "type": "extraction", "used": ["source:pmid:37854191"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE108886", "type": "accession_resolution", "used": ["data:geo:GSE108886"], "agent": "agent:corpus-assay-engine"}], "agents": [{"id": "agent:corpus-assay-engine", "type": "software", "label": "Corpus assay engine"}, {"id": "agent:minimax-reader", "type": "software", "label": "MiniMax literature reader"}], "source_run_id": "reading-sci-2026-11-01", "dsl_catalogue_fingerprint": "343fb30756ba6ec3", "capability_table_fingerprint": "2cb3e42f047d2ecf", "snapshot_id": "ba34b49490d9bfaf"}, "lineage": {"parent_id": "pmid:37854191", "parent_hash": "sha256:b2ab2239373d03b81f903eae938d78403406acf7f27496583d8d0a05eda2b342", "fork_of": null, "generation": 3}, "gaps": ["asset:data:geo:GSE108886 licence not read", "asset:data:geo:GSE108886 not downloaded", "asset:data:geo:GSE145467 licence not read", "asset:data:geo:GSE145467 not downloaded", "asset:data:geo:GSE45885 licence not read", "asset:data:geo:GSE45885 not downloaded", "asset:data:geo:GSE45887 licence not read", "asset:data:geo:GSE45887 not downloaded", "execution.not_run", "review.human_critique", "review.independent_reader", "review.machine_critique"], "assays": [{"role": "data_generation", "definition_id": "microarray", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.93, "evidence": ["as potential biomarkers. Methods: Three NOA microarray datasets (GSE45885, GSE108886, and GSE145467"], "order_basis": "paper_text_order:416", "id": "assay:1"}, {"role": "analysis", "definition_id": "machine_learning", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": ["train_classifier", "cross_validate", "roc_auc"], "support_state": "supported", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 1.0, "evidence": ["Finally, we defined the 5 most significant genes (IL20RB, C9orf117, HILS1, PAOX, and DZIP1), which were identified, by XGB, as potential biomarkers"], "order_basis": "consumes_earlier_output", "id": "assay:2"}]}