{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,7,27]],"date-time":"2026-07-27T15:00:03Z","timestamp":1785164403508,"version":"3.55.0"},"reference-count":60,"publisher":"MDPI AG","issue":"15","license":[{"start":{"date-parts":[[2026,7,25]],"date-time":"2026-07-25T00:00:00Z","timestamp":1784937600000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Int. J. Mol. Sci.","IJMS"],"abstract":"<jats:p>The leading-edge (LE) of hepatocellular carcinoma (HCC) is a critical region driving malignant progression and is closely associated with high patient mortality and marked intratumoral heterogeneity. Multi-omics integration identified elevated expression of SPARC and IGFBP7 in the LE region, which was associated with stromal remodeling-related transcriptional programs and an immune-depleted microenvironment. Cell-cell communication and pathway analyses further suggested potential links between LE-associated stromal states and pro-invasive signaling programs. Furthermore, we developed SpaPred, which demonstrated favorable performance in inferring the spatiotemporal heterogeneity of HCC at the spatial resolution. This model overcomes the limitations of existing algorithms in analyzing the composition of tissue spatial structures. Finally, integration of in silico trajectory-perturbation and pharmacogenomic drug-response analyses prioritized Oxaliplatin, Belinostat, and Temsirolimus as candidate compounds associated with LE-related transcriptional programs. These drug predictions are computational and require experimental validation. Collectively, SpaPred provides a hypothesis-generating framework for investigating spatial heterogeneity and candidate therapeutic vulnerabilities in HCC.<\/jats:p>","DOI":"10.3390\/ijms27156643","type":"journal-article","created":{"date-parts":[[2026,7,27]],"date-time":"2026-07-27T09:47:41Z","timestamp":1785145661000},"page":"6643","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":0,"title":["Deciphering the Leading-Edge Spatiotemporal Microenvironment of Hepatocellular Carcinoma for Targeted Drug Discovery Using SpaPred"],"prefix":"10.3390","volume":"27","author":[{"ORCID":"https:\/\/orcid.org\/0009-0000-3653-0503","authenticated-orcid":false,"given":"Shibo","family":"Zhang","sequence":"first","affiliation":[{"name":"College of Bioinformatics Science and Technology, Harbin Medical University, Harbin 150081, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Ziqiao","family":"Li","sequence":"additional","affiliation":[{"name":"College of Bioinformatics Science and Technology, Harbin Medical University, Harbin 150081, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Kexin","family":"Yu","sequence":"additional","affiliation":[{"name":"College of Bioinformatics Science and Technology, Harbin Medical University, Harbin 150081, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Guang","family":"Shi","sequence":"additional","affiliation":[{"name":"College of Bioinformatics Science and Technology, Harbin Medical University, Harbin 150081, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Yangguang","family":"Su","sequence":"additional","affiliation":[{"name":"College of Bioinformatics Science and Technology, Harbin Medical University, Harbin 150081, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Xin","family":"Hu","sequence":"additional","affiliation":[{"name":"College of Bioinformatics Science and Technology, Harbin Medical University, Harbin 150081, China"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-2423-8569","authenticated-orcid":false,"given":"Xiujie","family":"Chen","sequence":"additional","affiliation":[{"name":"College of Bioinformatics Science and Technology, Harbin Medical University, Harbin 150081, China"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1968","published-online":{"date-parts":[[2026,7,25]]},"reference":[{"key":"ref_1","first-page":"17","article-title":"Cancer statistics, 2023","volume":"73","author":"Siegel","year":"2023","journal-title":"CA Cancer J. Clin."},{"key":"ref_2","doi-asserted-by":"crossref","first-page":"346","DOI":"10.1038\/s41582-023-00809-y","article-title":"Single-cell and spatial transcriptomics: Deciphering brain complexity in health and disease","volume":"19","author":"Piwecka","year":"2023","journal-title":"Nat. Rev. Neurol."},{"key":"ref_3","doi-asserted-by":"crossref","first-page":"211","DOI":"10.1038\/s41586-021-03634-9","article-title":"Exploring tissue architecture using spatial transcriptomics","volume":"596","author":"Rao","year":"2021","journal-title":"Nature"},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"188720","DOI":"10.1016\/j.bbcan.2022.188720","article-title":"Understanding initiation and progression of hepatocellular carcinoma through single cell sequencing","volume":"1877","author":"Li","year":"2022","journal-title":"Biochim. Biophys. Acta Rev. Cancer"},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"4163","DOI":"10.7150\/thno.71873","article-title":"Spatial maps of hepatocellular carcinoma transcriptomes reveal spatial expression patterns in tumor immune microenvironment","volume":"12","author":"Wang","year":"2022","journal-title":"Theranostics"},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"64","DOI":"10.1038\/s41568-020-00316-y","article-title":"Publisher Correction: Mutations and mechanisms of WNT pathway tumour suppressors in cancer","volume":"21","author":"Bugter","year":"2021","journal-title":"Nat. Rev. Cancer"},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"136","DOI":"10.1016\/j.semcancer.2022.09.004","article-title":"TGF-beta in developmental and fibrogenic EMTs","volume":"86","author":"Lee","year":"2022","journal-title":"Semin. Cancer Biol."},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"425","DOI":"10.1038\/s41392-021-00828-5","article-title":"Targeting PI3K\/Akt signal transduction for cancer therapy","volume":"6","author":"He","year":"2021","journal-title":"Signal Transduct. Target. Ther."},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"1410","DOI":"10.1172\/JCI114314","article-title":"Biochemical basis of synergy between antigen and T-helper (Th) cell-mediated activation of resting human B cells","volume":"84","author":"Chartash","year":"1989","journal-title":"J. Clin. Investig."},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"442","DOI":"10.1016\/j.cell.2020.03.048","article-title":"Single-Cell Analyses Inform Mechanisms of Myeloid-Targeted Therapies in Colon Cancer","volume":"181","author":"Zhang","year":"2020","journal-title":"Cell"},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"4235","DOI":"10.1016\/j.cell.2023.07.034","article-title":"A pan-cancer single-cell panorama of human natural killer cells","volume":"186","author":"Tang","year":"2023","journal-title":"Cell"},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"104","DOI":"10.1186\/s12943-022-01569-x","article-title":"Targeting TGF-\u03b2 signal transduction for fibrosis and cancer therapy","volume":"21","author":"Peng","year":"2022","journal-title":"Mol. Cancer"},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"3","DOI":"10.1038\/s41392-021-00762-6","article-title":"Wnt\/\u03b2-catenin signalling: Function, biological mechanisms, and therapeutic opportunities","volume":"7","author":"Liu","year":"2022","journal-title":"Signal Transduct. Target. Ther."},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"eadk2091","DOI":"10.1126\/sciadv.adk2091","article-title":"IRF8-mutant B cell lymphoma evades immunity through a CD74-dependent deregulation of antigen processing and presentation in MHCII complexes","volume":"10","author":"Qiu","year":"2024","journal-title":"Sci. Adv."},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"1281","DOI":"10.1016\/j.cell.2021.01.022","article-title":"Inhibitory CD161 receptor identified in glioma-infiltrating T cells by single-cell analysis","volume":"184","author":"Mathewson","year":"2021","journal-title":"Cell"},{"key":"ref_16","doi-asserted-by":"crossref","first-page":"585","DOI":"10.1038\/s41422-023-00831-1","article-title":"An invasive zone in human liver cancer identified by Stereo-seq promotes hepatocyte-tumor cell crosstalk, local immunosuppression and tumor progression","volume":"33","author":"Wu","year":"2023","journal-title":"Cell Res."},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"218","DOI":"10.1038\/s41392-021-00641-0","article-title":"Signaling pathways in cancer-associated fibroblasts and targeted therapy for cancer","volume":"6","author":"Wu","year":"2021","journal-title":"Signal Transduct. Target. Ther."},{"key":"ref_18","doi-asserted-by":"crossref","first-page":"1872","DOI":"10.1158\/0008-5472.CAN-23-2889","article-title":"SPARC Stabilizes ApoE to Induce Cholesterol-Dependent Invasion and Sorafenib Resistance in Hepatocellular Carcinoma","volume":"84","author":"Wan","year":"2024","journal-title":"Cancer Res."},{"key":"ref_19","doi-asserted-by":"crossref","unstructured":"Yi, X., Zheng, X., Xu, H., Li, J., Zhang, T., Ge, P., Liao, D., Li, H., Lyu, X., and Ai, J. (2022). IGFBP7 and the Tumor Immune Landscape: A Novel Target for Immunotherapy in Bladder Cancer. Front. Immunol., 13.","DOI":"10.3389\/fimmu.2022.898493"},{"key":"ref_20","doi-asserted-by":"crossref","first-page":"134","DOI":"10.1016\/j.semcancer.2021.02.015","article-title":"Cellular heterogeneity and plasticity in liver cancer","volume":"82","author":"Chan","year":"2022","journal-title":"Semin. Cancer Biol."},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"1951","DOI":"10.1158\/1078-0432.CCR-14-0122","article-title":"Intratumor Heterogeneity in Hepatocellular Carcinoma","volume":"21","author":"Friemel","year":"2015","journal-title":"Clin. Cancer Res."},{"key":"ref_22","first-page":"397","article-title":"Cloud-Based GWAS Platform: An Innovative Solution for Efficient Acquisition and Analysis of Genomic Data","volume":"1","author":"Ke","year":"2025","journal-title":"Med. Res."},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"796","DOI":"10.1038\/s41746-025-02172-2","article-title":"Machine learning identifies TIME subtypes linking EGFR mutations and immune states in lung adenocarcinoma","volume":"8","author":"Gong","year":"2025","journal-title":"npj Digit. Med."},{"key":"ref_24","doi-asserted-by":"crossref","first-page":"1000","DOI":"10.21147\/j.issn.1000-9604.2025.06.11","article-title":"Metabolic reprogramming signature predicts immunotherapy efficacy in lung adenocarcinoma: Targeting SLC25A1 to overcome immune resistance","volume":"37","author":"Zhang","year":"2025","journal-title":"Chin. J. Cancer Res."},{"key":"ref_25","doi-asserted-by":"crossref","first-page":"247","DOI":"10.1158\/2326-6066.CIR-23-0440","article-title":"Engagement of CD99 Activates Distinct Programs in Ewing Sarcoma and Macrophages","volume":"12","author":"Manara","year":"2024","journal-title":"Cancer Immunol. Res."},{"key":"ref_26","doi-asserted-by":"crossref","first-page":"e2003558","DOI":"10.1002\/advs.202003558","article-title":"Multi-Arm PEG\/Peptidomimetic Conjugate Inhibitors of DR6\/APP Interaction Block Hematogenous Tumor Cell Extravasation","volume":"8","author":"Wang","year":"2021","journal-title":"Adv. Sci."},{"key":"ref_27","doi-asserted-by":"crossref","first-page":"6","DOI":"10.1186\/s13046-015-0121-1","article-title":"Characterization of genome-wide TFCP2 targets in hepatocellular carcinoma: Implication of targets FN1 and TJP1 in metastasis","volume":"34","author":"Xu","year":"2015","journal-title":"J. Exp. Clin. Cancer Res."},{"key":"ref_28","doi-asserted-by":"crossref","first-page":"2757","DOI":"10.7150\/thno.94537","article-title":"Targeting tumor cell-to-macrophage communication by blocking Vtn-C1qbp interaction inhibits tumor progression via enhancing macrophage phagocytosis","volume":"14","author":"Zhang","year":"2024","journal-title":"Theranostics"},{"key":"ref_29","doi-asserted-by":"crossref","first-page":"5029","DOI":"10.1038\/s41467-023-40271-4","article-title":"Spatial transcriptomics reveals distinct and conserved tumor core and edge architectures that predict survival and targeted therapy response","volume":"14","author":"Arora","year":"2023","journal-title":"Nat. Commun."},{"key":"ref_30","doi-asserted-by":"crossref","first-page":"1053","DOI":"10.1023\/A:1008213732429","article-title":"Oxaliplatin: A review of preclinical and clinical studies","volume":"9","author":"Raymond","year":"1998","journal-title":"Ann. Oncol."},{"key":"ref_31","doi-asserted-by":"crossref","first-page":"3361","DOI":"10.1200\/JCO.2011.41.2395","article-title":"Epigenetic therapy using belinostat for patients with unresectable hepatocellular carcinoma: A multicenter phase I\/II study with biomarker and pharmacokinetic analysis of tumors from patients in the Mayo Phase II Consortium and the Cancer Therapeutics Research Group","volume":"30","author":"Yeo","year":"2012","journal-title":"J. Clin. Oncol."},{"key":"ref_32","doi-asserted-by":"crossref","first-page":"561","DOI":"10.1159\/000518297","article-title":"Phase II Trial of the Combination of Temsirolimus and Sorafenib in Advanced Hepatocellular Carcinoma with Tumor Mutation Profiling","volume":"10","author":"Kelley","year":"2021","journal-title":"Liver Cancer"},{"key":"ref_33","doi-asserted-by":"crossref","first-page":"93","DOI":"10.1007\/978-1-4939-3578-9_5","article-title":"The Gene Expression Omnibus Database","volume":"1418","author":"Clough","year":"2016","journal-title":"Methods Mol. Biol."},{"key":"ref_34","doi-asserted-by":"crossref","first-page":"293","DOI":"10.1038\/s41587-023-01767-y","article-title":"Dictionary learning for integrative, multimodal and scalable single-cell analysis","volume":"42","author":"Hao","year":"2024","journal-title":"Nat. Biotechnol."},{"key":"ref_35","doi-asserted-by":"crossref","first-page":"1289","DOI":"10.1038\/s41592-019-0619-0","article-title":"Fast, sensitive and accurate integration of single-cell data with Harmony","volume":"16","author":"Korsunsky","year":"2019","journal-title":"Nat. Methods"},{"key":"ref_36","doi-asserted-by":"crossref","first-page":"1396","DOI":"10.1126\/science.1254257","article-title":"Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma","volume":"344","author":"Patel","year":"2014","journal-title":"Science"},{"key":"ref_37","doi-asserted-by":"crossref","first-page":"494","DOI":"10.1038\/s41586-018-0414-6","article-title":"RNA velocity of single cells","volume":"560","author":"Soldatov","year":"2018","journal-title":"Nature"},{"key":"ref_38","doi-asserted-by":"crossref","unstructured":"Gene Ontology Consortium, Aleksander, S.A., Balhoff, J., Carbon, S., Cherry, J.M., Drabkin, H.J., Ebert, D., Feuermann, M., Gaudet, P., and Harris, N.L. (2023). The Gene Ontology knowledgebase in 2023. Genetics, 224, iyad031.","DOI":"10.1093\/genetics\/iyad031"},{"key":"ref_39","doi-asserted-by":"crossref","first-page":"27","DOI":"10.1093\/nar\/28.1.27","article-title":"KEGG: Kyoto Encyclopedia of Genes and Genomes","volume":"28","author":"Kanehisa","year":"1999","journal-title":"Nucleic Acids Res."},{"key":"ref_40","doi-asserted-by":"crossref","first-page":"1349","DOI":"10.1038\/s41587-022-01273-7","article-title":"Spatially informed cell-type deconvolution for spatial transcriptomics","volume":"40","author":"Ma","year":"2022","journal-title":"Nat. Biotechnol."},{"key":"ref_41","doi-asserted-by":"crossref","unstructured":"Gaujoux, R., and Seoighe, C. (2010). A flexible R package for nonnegative matrix factorization. BMC Bioinform., 11.","DOI":"10.1186\/1471-2105-11-367"},{"key":"ref_42","doi-asserted-by":"crossref","first-page":"e50","DOI":"10.1093\/nar\/gkab043","article-title":"SPOTlight: Seeded NMF regression to deconvolute spatial transcriptomics spots with single-cell transcriptomes","volume":"49","author":"Nieto","year":"2021","journal-title":"Nucleic Acids Res."},{"key":"ref_43","doi-asserted-by":"crossref","first-page":"517","DOI":"10.1038\/s41587-021-00830-w","article-title":"Robust decomposition of cell type mixtures in spatial transcriptomics","volume":"40","author":"Cable","year":"2022","journal-title":"Nat. Biotechnol."},{"key":"ref_44","doi-asserted-by":"crossref","first-page":"661","DOI":"10.1038\/s41587-021-01139-4","article-title":"Cell2location maps fine-grained cell types in spatial transcriptomics","volume":"40","author":"Kleshchevnikov","year":"2022","journal-title":"Nat. Biotechnol."},{"key":"ref_45","doi-asserted-by":"crossref","first-page":"1352","DOI":"10.1038\/s41592-021-01264-7","article-title":"Deep learning and alignment of spatially resolved single-cell transcriptomes with Tangram","volume":"18","author":"Biancalani","year":"2021","journal-title":"Nat. Methods"},{"key":"ref_46","doi-asserted-by":"crossref","first-page":"78","DOI":"10.1186\/s13059-021-02286-2","article-title":"Giotto: A toolbox for integrative analysis and visualization of spatial expression data","volume":"22","author":"Dries","year":"2021","journal-title":"Genome Biol."},{"key":"ref_47","doi-asserted-by":"crossref","first-page":"1408","DOI":"10.1038\/s41587-020-0591-3","article-title":"Generalizing RNA velocity to transient cell states through dynamical modeling","volume":"38","author":"Bergen","year":"2020","journal-title":"Nat. Biotechnol."},{"key":"ref_48","doi-asserted-by":"crossref","first-page":"690","DOI":"10.1016\/j.cell.2021.12.045","article-title":"Mapping transcriptomic vector fields of single cells","volume":"185","author":"Qiu","year":"2022","journal-title":"Cell"},{"key":"ref_49","doi-asserted-by":"crossref","first-page":"1244","DOI":"10.1093\/bioinformatics\/btv723","article-title":"PharmacoGx: An R package for analysis of large pharmacogenomic datasets","volume":"32","author":"Smirnov","year":"2016","journal-title":"Bioinformatics"},{"key":"ref_50","doi-asserted-by":"crossref","first-page":"D1227","DOI":"10.1093\/nar\/gkad1040","article-title":"DGIdb 5.0: Rebuilding the drug\u2013gene interaction database for precision medicine and drug discovery platforms","volume":"52","author":"Cannon","year":"2023","journal-title":"Nucleic Acids Res."},{"key":"ref_51","doi-asserted-by":"crossref","first-page":"bbab260","DOI":"10.1093\/bib\/bbab260","article-title":"oncoPredict: An R package for predicting in vivo or cancer patient drug response and biomarkers from cell line screening data","volume":"22","author":"Maeser","year":"2021","journal-title":"Brief. Bioinform."},{"key":"ref_52","doi-asserted-by":"crossref","first-page":"D1328","DOI":"10.1093\/nar\/gkae883","article-title":"Comparative toxicogenomics database\u2019s 20th anniversary: Update 2025","volume":"53","author":"Davis","year":"2024","journal-title":"Nucleic Acids Res."},{"key":"ref_53","doi-asserted-by":"crossref","first-page":"D1075","DOI":"10.1093\/nar\/gkv1075","article-title":"The SIDER database of drugs and side effects","volume":"44","author":"Kuhn","year":"2016","journal-title":"Nucleic Acids Res."},{"key":"ref_54","doi-asserted-by":"crossref","first-page":"D1265","DOI":"10.1093\/nar\/gkad976","article-title":"DrugBank 6.0: The DrugBank Knowledgebase for 2024","volume":"52","author":"Knox","year":"2024","journal-title":"Nucleic Acids Res."},{"key":"ref_55","doi-asserted-by":"crossref","first-page":"D1180","DOI":"10.1093\/nar\/gkad1004","article-title":"The ChEMBL Database in 2023: A drug discovery platform spanning multiple bioactivity data types and time periods","volume":"52","author":"Zdrazil","year":"2024","journal-title":"Nucleic Acids Res."},{"key":"ref_56","doi-asserted-by":"crossref","first-page":"D1373","DOI":"10.1093\/nar\/gkac956","article-title":"PubChem 2023 update","volume":"51","author":"Kim","year":"2022","journal-title":"Nucleic Acids Res."},{"key":"ref_57","doi-asserted-by":"crossref","first-page":"59","DOI":"10.1038\/s41568-024-00763-x","article-title":"The present and future of the Cancer Dependency Map","volume":"25","author":"Arafeh","year":"2025","journal-title":"Nat. Rev. Cancer"},{"key":"ref_58","first-page":"100141","article-title":"clusterProfiler 4.0: A universal enrichment tool for interpreting omics data","volume":"2","author":"Wu","year":"2021","journal-title":"Innovation"},{"key":"ref_59","doi-asserted-by":"crossref","first-page":"e113","DOI":"10.1093\/nar\/gkaa802","article-title":"Stable gene expression for normalisation and single-sample scoring","volume":"48","author":"Bhuva","year":"2020","journal-title":"Nucleic Acids Res."},{"key":"ref_60","doi-asserted-by":"crossref","unstructured":"Jin, S., Plikus, M.V., and Nie, Q. (2023). CellChat for systematic analysis of cell-cell communication from single-cell and spatially resolved transcriptomics. bioRxiv.","DOI":"10.1101\/2023.11.05.565674"}],"container-title":["International Journal of Molecular Sciences"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/1422-0067\/27\/15\/6643\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,7,27]],"date-time":"2026-07-27T14:08:45Z","timestamp":1785161325000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/1422-0067\/27\/15\/6643"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,7,25]]},"references-count":60,"journal-issue":{"issue":"15","published-online":{"date-parts":[[2026,8]]}},"alternative-id":["ijms27156643"],"URL":"https:\/\/doi.org\/10.3390\/ijms27156643","relation":{},"ISSN":["1422-0067"],"issn-type":[{"value":"1422-0067","type":"electronic"}],"subject":[],"published":{"date-parts":[[2026,7,25]]}}}