Biomarker Identification in Pancreatic Cancer Through Concordant Differential Expression and Interpretable Machine Learning Analyses - Maciá-Escalante, Sonia; López-Aladid, Rubén; Muñoz-Tovar, Rebeca; López-Herrero, Manuel; Navarro-Sellés, Ana; Garmendia, Leonor; Puerto, Carolina; Fossati-Vázquez, Mariana; Parente-Muñoz, Pablo (medRxiv, 2026)
Background Pancreatic ductal adenocarcinoma is one of the most aggressive and lethal malignancies of the gastrointestinal tract. The poor prognosis is largely attributed to late-stage diagnosis, pronounced tumor heterogeneity, and limited therapeutic efficacy. These challenges underscore the urgent need for the identification of robust molecular biomarkers and novel therapeutic targets.Methods Gene expression data from a total of 146 pancreatic tissue samples, comprising 72 normal and 74 tumor specimens obtained from the Pan-Cancer Atlas(TCGA) were analyzed. Differential gene expression analysis was conducted using the DESeq2 package, followed by functional enrichment analysis based on GO and KEGG. A classification model was developed using the XGBoost algorithm and evaluated through 500 bootstrapping iterations and 5-fold cross-validation to ensure robustness and generalizability. Model interpretability was assessed using SHAP (SHapley Additive exPlanations) values to identify genes with the highest predictive contribution.Results A comprehensive transcriptomic analysis revealed significant dysregulation of multiple genes between normal and tumor pancreatic tissues. Genes such as GJB3, S100A2, MSLN, and SLC2A1 were notably overexpressed, whereas DEFA6, APOB, and RBP2 exhibited marked downregulation, indicative of impaired exocrine function and aberrant epithelial reprogramming. The XGBoost classification model achieved an average area under the curve (AUC) of 0.9868 and an overall accuracy of 98.6%. SHAP (SHapley Additive exPlanations) analysis identified GJB3, LINC02086, and TSPAN1 as key predictive features. Six genes were concurrently identified as differentially expressed and highly influential within the model, supporting their potential utility as robust biomarkers for pancreatic tumor characterization.Conclusions Pancreatic ductal adenocarcinoma is marked by extensive transcriptomic reprogramming. The integration of differential gene expression analysis with interpretable machine learning enabled the identification of a molecular signature with potential diagnostic and therapeutic relevance.Competing Interest StatementThe authors have declared no competing interest.Funding StatementThis study did not receive any fundingAuthor DeclarationsI confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.YesI confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.YesI understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).YesI have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable.YesDe-identified gene-level resources and analytic artifacts generated for this study will be made publicly available upon publication at a persistent repository (DOI to be provided upon acceptance). Released files will include, subject to data-use permissions, gene-level normalized expression matrices; the full differential expression table (deseq2_results_annotated.csv); Gene Ontology and KEGG over-representation results (downregulated_GO_Biological_Process_2021_results.csv) and upregulated_GO_Biological_Process_2021_results.csv); SHAP summaries aggregated across evaluation runs (shap_values_bootstrap_summary.csv) and per-bootstrap classification metrics (Classification Metrics Across 500 Bootstraps.csv). Reproducible code will be deposited alongside the data, comprising R scripts for DESeq2 and clusterProfiler workflows and Python scripts for model training and explainability with XGBoost and SHAP.