Predicting the Risk of Lumbar Prolapsed Disc: A Gene Signature-Based Machine Learning Analysis.

{"sci_version": "1.2", "schema": "https://corpus.bntrg.co.uk/sci/1.0/sci.schema.json", "id": "pmid:40319430", "kind": "science-record", "revision": 3, "content_hash": "sha256:d2e8c07f9dbb940cca31ecf58edf490e6a8cf428df16a4664f2874d84381b315", "created_at": "2026-11-01T20:30:28Z", "updated_at": "2026-11-01T22:14:01Z", "generator": {"name": "corpus-assay-engine", "run_id": "reading-sci-2026-11-01"}, "paper": {"state": "complete", "title": "Predicting the Risk of Lumbar Prolapsed Disc: A Gene Signature-Based Machine Learning Analysis.", "doi": "10.1007/s40122-025-00744-4", "pmid": "40319430", "authors": ["Fengfeng Wang", "Fei Meng", "Stanley Sau Ching Wong"], "journal": "Pain and therapy", "year": 2025, "abstract": "INTRODUCTION: Lumbar prolapsed disc (LPD) is a leading cause of low back pain, contributing significantly to global disability and healthcare burden. This study aimed to develop machine learning models to predict the risk of LPD by analysing gene expression profiles for early detection.\nMETHODS: Transcriptomic data from peripheral blood samples were obtained from the Gene Expression Omnibus (GEO) database, with dataset GSE150408 used for training and GSE124272 for testing. The training dataset included 17 patients with sciatica resulting from LPD, all of whom had magnetic resonance imaging confirmation of single-level LPD at either the L4/5 or L5/S1 levels. Data from 17 healthy volunteers were used as controls. Recursive feature elimination (RFE) was employed to identify the most relevant gene signatures among 23 pain-related genes. Machine learning models, including support vector machine (SVM), random forest, k-nearest neighbours (KNN), logistic regression, and Extreme Gradient Boosting (XGBoost), were trained and evaluated. Model performance was assessed using accuracy, area under the curve (AUC), F1 score, and Matthews correlation coefficient (MCC).\nRESULTS: Eight key gene signatures were identified as significant predictors of LPD, with MMP9 exhibiting the highest importance score. Most of these genes were differentially expressed between patients with LPD and healthy controls (p < 0.05). Among the models, random forest demonstrated the highest accuracy (0.80, 95% CI 0.73-0.85) and MCC (0.64, 95% CI 0.53-0.76), followed by KNN, XGBoost, and SVM. Overall, the random forest model exhibited the most robust performance in predicting the risk of LPD.\nCONCLUSION: The results of our study suggest that machine learning models based on pain-related gene signatures may identify patients at high risk of developing LPD with reasonably high accuracy. These prediction models could perhaps be integrated into clinical diagnostic tools to enhance early diagnosis and prevention.", "claims": [{"id": "claim-method", "statement": "We used t test or Mann–Whitney U test to compare the expression levels of selected gene signatures between two groups, with significance determined by a p value less than 0.05.", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:40319430", "locator": "sections[10] \"Statistical Analysis\" paragraph 2", "quote": "We used t test or Mann–Whitney U test to compare the expression levels of selected gene signatures between two groups, with significance determined by a p value less than 0.05."}], "state": "complete"}, {"id": "claim-data-access", "statement": "The datasets used and/or analysed in the current study were obtained from the publicly available Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE150408 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE124272).", "kind": "methodological", "target": null, "evidence": [{"asset": "source:pmid:40319430", "locator": "sections[21] \"notes\" paragraph 1", "quote": "The datasets used and/or analysed in the current study were obtained from the publicly available Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE150408 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE124272)."}], "state": "complete"}, {"id": "claim-target-result", "statement": "Among the models, random forest demonstrated the highest accuracy (0.80, 95% CI 0.73–0.85) and MCC (0.64, 95% CI 0.53–0.76), followed by KNN, XGBoost, and SVM.", "kind": "quantitative", "target": {"value": 95.0, "unit": "per cent", "tolerance": null, "source": "Among the models, random forest demonstrated the highest accuracy (0.80, 95% CI 0.73–0.85) and MCC (0.64, 95% CI 0.53–0.76), followed by KNN, XGBoost, and SVM."}, "evidence": [{"asset": "source:pmid:40319430", "locator": "sections[0] \"Abstract\" paragraph 3", "quote": "Among the models, random forest demonstrated the highest accuracy (0.80, 95% CI 0.73–0.85) and MCC (0.64, 95% CI 0.53–0.76), followed by KNN, XGBoost, and SVM."}], "state": "complete"}], "gaps": []}, "experiment": {"state": "partial", "objective": "Among the models, random forest demonstrated the highest accuracy (0.80, 95% CI 0.73–0.85) and MCC (0.64, 95% CI 0.53–0.76), followed by KNN, XGBoost, and SVM.", "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:40319430"], "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:GSE124272"], "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:machine_learning_random_forest", "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:GSE124272"], "outputs": ["the replicated result and its comparison"], "instruments": [], "reagents": []}], "materials": [], "instruments": [], "reagents": [], "conditions": {}, "gaps": ["step:run_analysis blocked on method:machine_learning_random_forest"]}, "data": {"state": "partial", "assets": [{"id": "source:pmid:40319430", "role": "manuscript", "format": "xml", "media_type": "application/xml", "uri": "https://doi.org/10.1007/s40122-025-00744-4", "sha256": "sha256:2717b409446d1ed76f3f0f5353ea608a740cb4786a4badf9b32549d56d4cacb9", "bytes": 31518, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "10.1007/s40122-025-00744-4", "version": null, "note": "Full article fetched from PMC; the checksum covers the rendered section body, not the original XML bytes.", "structure": {"pmcid": "PMC12085505", "source": "europe_pmc_fulltext", "coverage": "whole_body_without_references", "text_hash_scope": "rendered_section_body"}}, {"id": "data:geo:GSE124272", "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=GSE124272", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE124272", "version": null, "note": "The datasets used and/or analysed in the current study were obtained from the publicly available Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE150408 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE124272).", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The datasets used and/or analysed in the current study were obtained from the publicly available Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE150408 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE124272).", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM3526782", "GSM3526783", "GSM3526784", "GSM3526785", "GSM3526786", "GSM3526787", "GSM3526788", "GSM3526789", "GSM3526790", "GSM3526791", "GSM3526792", "GSM3526793", "GSM3526794", "GSM3526795", "GSM3526796", "GSM3526797"], "sample_count": 16, "group_sizes": null, "group_sizes_basis": null, "platform": "21185", "taxon": "Homo sapiens", "title": "Transcriptome signatures reveal candidate key genes in the whole blood of patients with lumbar disc prolapse", "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:GSE150408", "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=GSE150408", "sha256": null, "bytes": null, "access": "open", "state": "resolved", "reader": null, "reader_status": "reader_required", "identifier": "GSE150408", "version": null, "note": "The datasets used and/or analysed in the current study were obtained from the publicly available Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE150408 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE124272).", "structure": {"repository": "GEO", "access_route": "GEO accession URL", "access_statement_verbatim": "The datasets used and/or analysed in the current study were obtained from the publicly available Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE150408 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE124272).", "access_statement_kind": "data_availability_statement", "sample_accessions": ["GSM4548691", "GSM4548692", "GSM4548693", "GSM4548694", "GSM4548695", "GSM4548696", "GSM4548697", "GSM4548698", "GSM4548699", "GSM4548700", "GSM4548701", "GSM4548702", "GSM4548703", "GSM4548704", "GSM4548705", "GSM4548706", "GSM4548707", "GSM4548708", "GSM4548709", "GSM4548710", "GSM4548711", "GSM4548712", "GSM4548713", "GSM4548714", "GSM4548715", "GSM4548716", "GSM4548717", "GSM4548718", "GSM4548719", "GSM4548720", "GSM4548721", "GSM4548722", "GSM4548723", "GSM4548724", "GSM4548725", "GSM4548726", "GSM4548727", "GSM4548728", "GSM4548729", "GSM4548730", "GSM4548731", "GSM4548732", "GSM4548733", "GSM4548734", "GSM4548735", "GSM4548736", "GSM4548737", "GSM4548738", "GSM4548739", "GSM4548740", "GSM4548741", "GSM4548742", "GSM4548743", "GSM4548744", "GSM4548745", "GSM4548746", "GSM4548747", "GSM4548748", "GSM4548749"], "sample_count": 59, "group_sizes": null, "group_sizes_basis": null, "platform": "21185", "taxon": "Homo sapiens", "title": "Identification of key genes potentially related to intervertebral disc degeneration by microarray analysis", "resolution_checked_at": "2026-11-01T17:30:42Z", "resolution_check": "identifier_resolution_only_no_download_no_licence_read", "downloaded": false, "licence_read": false}}], "gaps": ["asset:data:geo:GSE124272 not downloaded", "asset:data:geo:GSE124272 licence not read", "asset:data:geo:GSE150408 not downloaded", "asset:data:geo:GSE150408 licence not read"]}, "analysis": {"state": "partial", "methods": [{"id": "method-1", "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-2", "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-3", "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": 0.7, "band": "medium", "quote": "Gradient Boosting"}}, {"id": "method-4", "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-5", "family": null, "specific": null, "label": "t test or Mann–Whitney U test", "dsl_operations": [], "missing_operations": [], "expressible": null, "parameters": {}, "data_shapes": [], "evidence": {"score": 1.0, "band": "verbatim", "quote": "We used t test or Mann–Whitney U test to compare the expression levels of selected gene signatures between two groups, with significance determined by a p value less than 0.05."}}], "targets": [{"id": "target-published-estimate", "claim_id": "claim-target-result", "operation": "cross_validate", "expected_value": 95.0, "unit": "per cent", "tolerance": null, "state": "complete", "note": "The published estimate as the reader quoted it from the paper. No tolerance is published because the paper states none."}], "gaps": ["method.reader_named_unmapped", "method.machine_learning_random_forest"]}, "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": false, "required_operations": [], "missing_operations": ["machine_learning_random_forest"], "basis": "The paper names the method it used and the corpus's method extractor maps it onto the DSL catalogue: random_forest."}, "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": ["GSE150408", "GSE124272"], "resolved_accessions": ["GSE150408", "GSE124272"], "statements": ["The datasets used and/or analysed in the current study were obtained from the publicly available Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE150408 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE124272)."], "basis": "The paper names GSE150408, GSE124272 in Gene Expression Omnibus (GEO) database and the accession resolves in the repository. No dataset was downloaded and no licence was read."}}, "can_run": {"dsl_operations": [], "assets": ["source:pmid:40319430", "data:geo:GSE124272"], "steps": []}, "must_build": {"dsl_operations": [{"operation": "machine_learning_random_forest", "unlocks_papers": null, "priority": null}], "converters": [], "readers": [], "capabilities": [], "data_gaps": ["data:geo:GSE124272"]}, "gaps": []}, "review": {"state": "partial", "machine_critique": {"kind": "extraction_review", "reader": {"model": "MiniMax-M3.1-Flash", "provider": "minimax", "provider_family": "minimax", "read_at": "2026-11-01T19:51:41Z", "independence": "single_reader"}, "evidence_check": {"method": "verbatim", "availability": "verbatim"}, "calls": 5, "usage": {"prompt_tokens": 8991, "completion_tokens": 1733, "cached_tokens": 264}}, "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:40319430", "type": "paper", "label": "pmid:40319430"}, {"id": "source:pmid:40319430", "type": "manuscript", "label": "PMID 40319430 full text"}, {"id": "data:geo:GSE124272", "type": "dataset", "label": "GEO GSE124272"}], "activities": [{"id": "activity-reading-sci-2026-11-01", "type": "classification", "used": [{"id": "pmid:40319430", "type": "paper", "label": "pmid:40319430"}, {"id": "source:pmid:40319430", "type": "manuscript", "label": "PMID 40319430 full text"}, {"id": "data:geo:GSE124272", "type": "dataset", "label": "GEO GSE124272"}], "agent": "agent:corpus-assay-engine"}, {"id": "activity-read-40319430", "type": "extraction", "used": ["source:pmid:40319430"], "agent": "agent:minimax-reader"}, {"id": "activity-resolve-GSE124272", "type": "accession_resolution", "used": ["data:geo:GSE124272"], "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:40319430", "parent_hash": "sha256:4b785bc375104eaca96765e50aa698262975c1c2b269a0224d3f59aaa73fda55", "fork_of": null, "generation": 2}, "gaps": ["asset:data:geo:GSE124272 licence not read", "asset:data:geo:GSE124272 not downloaded", "asset:data:geo:GSE150408 licence not read", "asset:data:geo:GSE150408 not downloaded", "execution.not_run", "method.machine_learning_random_forest", "method.reader_named_unmapped", "review.human_critique", "review.independent_reader", "review.machine_critique", "step:run_analysis blocked on method:machine_learning_random_forest"], "assays": [{"role": "analysis", "definition_id": "machine_learning", "definition_version": "f2132e96ae6c7a4e", "definition_digest": null, "required_operations": [], "support_state": "partial", "source_hashes": {"dsl_catalogue": "f2132e96ae6c7a4e"}, "confidence": 0.7, "evidence": ["random forest"], "order_basis": "analysis_only", "id": "assay:1"}]}