Integrated multiple microarray studies by robust rank aggregation to identify immune-associated biomarkers in Crohn's disease based on three machine learning methods.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:36792688", "kind": "science-record", "revision": 3, "content_hash": "sha256:1e2c30b1f63cfba451aabc6ebaec62ddfe8c9c3f7032a936bd96ec0f956b1f7f", "created_at": "2026-11-01T19:01:23Z", "updated_at": "2026-11-01T22:08:49Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Integrated multiple microarray studies by robust rank aggregation to identify immune-associated biomarkers in Crohn's disease based on three machine learning methods.", "doi": "10.1038/s41598-022-26345-1", "pmid": "36792688", "authors": ["Zi-An Chen", "Hui-Hui Ma", "Yan Wang", "Hui Tian", "Jian-Wei Mi", "Dong-Mei Yao", "Chuan-Jie Yang"], "journal": "Scientific reports", "year": 2023, "abstract": "Crohn's disease (CD) is a complex autoimmune disorder presumed to be driven by complex interactions of genetic, immune, microbial and even environmental factors. Intrinsic molecular mechanisms in CD, however, remain poorly understood. The identification of novel biomarkers in CD cases based on larger samples through machine learning approaches may inform the diagnosis and treatment of diseases. A comprehensive analysis was conducted on all CD datasets of Gene Expression Omnibus (GEO); our team then used the robust rank aggregation (RRA) method to identify differentially expressed genes (DEGs) between controls and CD patients. PPI (protein‒protein interaction) network and functional enrichment analyses were performed to investigate the potential functions of the DEGs, with molecular complex detection (MCODE) identifying some important functional modules from the PPI network. Three machine learning algorithms, support vector machine-recursive feature elimination (SVM-RFE), random forest (RF), and least absolute shrinkage and selection operator (LASSO), were applied to determine characteristic genes, which were verified by ROC curve analysis and immunohistochemistry (IHC) using clinical samples. Univariable and multivariable logistic regression were used to establish a machine learning score for diagnosis. Single-sample GSEA (ssGSEA) was performed to examine the correlation between immune infiltration and biomarkers. In total, 5 datasets met the inclusion criteria: GSE75214, GSE95095, GSE126124, GSE179285, and GSE186582. Based on RRA integrated analysis, 203 significant DEGs were identified (120 upregulated genes and 83 downregulated genes), and MCODE revealed some important functional modules in the PPI network. Machine learning identified LCN2, REG1A, AQP9, CCL2, GIP, PROK2, DEFA5, CXCL9, and NAMPT; AQP9, PROK2, LCN2, and NAMPT were further verified by ROC curves and IHC in the external cohort. The final machine learning score was defined as [Expression level of AQP9 × (2.644)] + [Expression level of LCN2 × (0.958)] + [Expression level of NAMPT × (1.115)]. ssGSEA showed markedly elevated levels of dendritic cells and innate immune cells, such as macrophages and NK cells, in CD, consistent with the gene enrichment results that the DEGs are mainly involved in the IL-17 signaling pathway and humoral immune response. The selected biomarkers analyzed by the RRA method and machine learning are highly reliable. These findings improve our understanding of the molecular mechanisms of CD pathogenesis.", "claims": [{"id": "claim-method", "statement": "Using RRA, all genes were sorted for each dataset and ranked based on their logFC with the limma package.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:36792688", "locator": "sections[5] \"RRA analysis and identification of DEGs in the integrated cohort\" paragraph 2", "quote": "Using RRA, all genes were sorted for each dataset and ranked based on their logFC with the limma package."}], "state": "complete"}, {"id": "claim-data-access", "statement": "The datasets analyzed during the current research are all available in Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/). The data used to support the findings of this study are available from the corresponding author upon request.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:36792688", "locator": "sections[34] \"notes\" paragraph 1", "quote": "The datasets analyzed during the current research are all available in Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/). The data used to support the findings of this study are available from the corresponding author upon request."}], "state": "complete"}, {"id": "claim-target-result", "statement": "Based on RRA integrated analysis, 203 significant DEGs were identified (120 upregulated genes and 83 downregulated genes), and MCODE revealed some important functional modules in the PPI network.", "kind": "qualitative", "target": null, "evidence": [{"asset": "source:pmid:36792688", "locator": "sections[0] \"Abstract\" paragraph 1", "quote": "Based on RRA integrated analysis, 203 significant DEGs were identified (120 upregulated genes and 83 downregulated genes), and MCODE revealed some important functional modules in the PPI network."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "partial", "objective": "Based on RRA integrated analysis, 203 significant DEGs were identified (120 upregulated genes and 83 downregulated genes), and MCODE revealed some important functional modules in the PPI network.", "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:36792688"], "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": "blocked", "blocked_on": "data:accession_absent", "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": [], "outputs": ["the source data with its repository provenance"], "instruments": [], "reagents": []}, {"id": "run_analysis", "order": 3, "name": "Run the published analysis", "kind": "analysis", "state": "blocked", "blocked_on": "method:group_comparison_wilcoxon", "description": "Apply the paper's statistical method to the resolved series and compare the result with the published target.", "duration_hours": null, "parameters": {}, "inputs": [], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": ["step:run_analysis blocked on method:group_comparison_wilcoxon", "step:obtain_source_data blocked on data:accession_absent"]}, "data": {"state": "partial", "assets": [{"id": "source:pmid:36792688", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1038/s41598-022-26345-1", "sha256": "sha256:753c17d757ac4fc6e2643ec95c1113b006acd9f03c5e1ff86157da28dd998382", "bytes": 39683, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1038/s41598-022-26345-1", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC9931764", "source": "europe_pmc_fulltext", "coverage": "whole_body_without_references", "text_hash_scope": "rendered_section_body"}}], "gaps": []}, "analysis": {"state": "partial", "methods": [{"id": "method-1", "family": "group_comparison", "specific": "wilcoxon", "label": "Group comparison", "dsl_operations": [], "missing_operations": ["group_comparison_wilcoxon"], "expressible": false, "parameters": {}, "data_shapes": ["two_or_more_groups_continuous"], "evidence": {"score": 0.8, "band": "high", "quote": "Wilcoxon"}}, {"id": "method-2", "family": "differential_expression", "specific": "limma_voom", "label": "Differential expression analysis", "dsl_operations": ["normalise_counts", "differential_expression_fit", "fdr_bh"], "missing_operations": [], "expressible": true, "parameters": {"log2fc_threshold": "d ranked based on their logFC with the limma package"}, "data_shapes": ["count_matrix_with_condition"], "evidence": {"score": 0.8, "band": "high", "quote": "limma"}}, {"id": "method-3", "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.7, "band": "medium", "quote": "random forest"}}, {"id": "method-4", "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": "support vector machine"}}, {"id": "method-5", "family": "machine_learning", "specific": "decision_tree", "label": "Machine learning", "dsl_operations": ["train_classifier", "cross_validate", "roc_auc"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["labelled_feature_matrix"], "evidence": {"score": 0.8, "band": "high", "quote": "decision tree"}}, {"id": "method-6", "family": "regression_modelling", "specific": "logistic", "label": "Regression modelling", "dsl_operations": ["glm_logit_fit"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["predictors_and_outcome"], "evidence": {"score": 0.7, "band": "medium", "quote": "logistic regression"}}, {"id": "method-7", "family": "correlation_analysis", "specific": "spearman", "label": "Correlation analysis", "dsl_operations": ["spearman_rho"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["paired_numeric_columns"], "evidence": {"score": 0.8, "band": "high", "quote": "Spearman's correlation"}}, {"id": "method-8", "family": "network_analysis", "specific": "ppi_network", "label": "Network / pathway analysis", "dsl_operations": ["build_graph", "centrality_scores", "pathway_enrich"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["edge_list"], "evidence": {"score": 0.7, "band": "medium", "quote": "PPI"}}, {"id": "method-9", "family": "diagnostic_evaluation", "specific": "roc", "label": "Diagnostic test evaluation", "dsl_operations": ["confusion_counts", "roc_curve", "sensitivity_specificity", "solve_confusion_matrix", "classification_metrics", "published_precision_envelope"], "missing_operations": [], "expressible": true, "parameters": {}, "data_shapes": ["reference_standard_and_test"], "evidence": {"score": 0.7, "band": "medium", "quote": "ROC curve"}}, {"id": "method-10", "family": null, "specific": null, "label": "robust rank aggregation (RRA)", "dsl_operations": [], "missing_operations": [], "expressible": null, "parameters": {}, "data_shapes": [], "evidence": {"score": 1.0, "band": "verbatim", "quote": "Using RRA, all genes were sorted for each dataset and ranked based on their logFC with the limma package."}}], "targets": [], "gaps": ["method.reader_named_unmapped", "method.group_comparison_wilcoxon"]}, "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": "blocked", "tier": "R4", "blocking_axis": "data", "blocked_on": "data:accession_absent", "axes": {"method": {"tier": "R1", "expressible": false, "required_operations": [], "missing_operations": ["group_comparison_wilcoxon"], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: wilcoxon."}, "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": "R4", "reachable": false, "accessions": [], "resolved_accessions": [], "statements": ["The datasets analyzed during the current research are all available in Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/). The data used to support the findings of this study are available from the corresponding author upon request."], "basis": "The paper names no accession to its source data, so the data has to come from the authors."}}, "can_run": {"dsl_operations": [], "assets": ["source:pmid:36792688"], "steps": []}, "must_build": {"dsl_operations": [{"operation": "group_comparison_wilcoxon", "unlocks_papers": null, "priority": null}], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:None"]}, "gaps": ["data:accession_absent"]}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-11-01T18:27:32Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 11737, "completion_tokens": 1911, "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:36792688", "type": "paper", "label": "pmid:36792688"}, {"id": "source:pmid:36792688", "type": "manuscript", "label": "PMID 36792688 full text"}, {"id": "data:geo:None", "type": "dataset", "label": "GEO None"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:36792688", "type": "paper", "label": "pmid:36792688"}, {"id": "source:pmid:36792688", "type": "manuscript", "label": "PMID 36792688 full text"}, {"id": "data:geo:None", "type": "dataset", "label": "GEO None"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-36792688", "type": "extraction", "used": ["source:pmid:36792688"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-None", "type": "accession_resolution", "used": ["data:geo:None"], "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:36792688", "parent_hash": "sha256:cfebe952d047cdba270229bb8a19228fe8f631721b7b9c8925f1dce6e2d12bf6", "fork_of": null, "generation": 2}, "gaps": ["execution.not_run", "method.group_comparison_wilcoxon", "method.reader_named_unmapped", "review.human_critique", "review.independent_reader", "review.machine_critique", "step:obtain_source_data blocked on data:accession_absent", "step:run_analysis blocked on method:group_comparison_wilcoxon"], "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": ["Integrated multiple microarray studies by robust rank aggregation to identi"], "order_basis": "paper_text_order:0", "id": "assay:1"}, {"role": "data_generation", "definition_id": "ihc", "definition_version": "99e1195c7d6f", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"assay_lexicon": "99e1195c7d6f"}, "confidence": 0.99, "evidence": ["hich were verified by ROC curve analysis and immunohistochemistry (IHC) using clinical samples. Univariable an", "ROC curve analysis and immunohistochemistry (IHC) using clinical samples. Univariable and mul", "AMPT were further verified by ROC curves and IHC in the external cohort. The final machine le", "tis and CD tissue specimens were examined by IHC for the expression levels of AQP9 and PROK2."], "order_basis": "paper_text_order:1283", "id": "assay:2"}, {"role": "analysis", "definition_id": "group_comparison", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 0.8, "evidence": ["Wilcoxon"], "order_basis": "consumes_earlier_output", "id": "assay:3"}]}