A New Era of Data-Driven Cancer Research and Care: Opportunities and Challenges.

People diagnosed with cancer and their formal and informal caregivers are increasingly faced with a deluge of complex information, thanks to rapid advancements in the type and volume of diagnostic, prognostic, and treatment data. This commentary discusses the opportunities and challenges that the society faces as we integrate large volumes of data into regular cancer care.

Spatial multi-omics reveal intratumoral humoral immunity niches associated with tertiary lymphoid structures in pancreatic cancer immunotherapy pathologic responders.

Pancreatic adenocarcinoma (PDAC) is a rapidly progressing cancer that responds poorly to immunotherapies. Intratumoral tertiary lymphoid structures (TLS) have been associated with rare long-term PDAC survivors, but the role of TLS in PDAC and their spatial relationships within the context of the broader tumor microenvironment remain unknown. We generated a spatial multi-omics atlas encompassing 26 PDAC tumors from patients treated with combination immunotherapies. Using machine learning-enabled H&E image classification models and unsupervised gene expression matrix factorization methods for spatial transcriptomics, we characterized cellular states within TLS niches spanning across distinct morphologies and immunotherapies. Unsupervised learning generated a TLS-specific spatial gene expression signature that significantly associates with improved survival in PDAC patients. These analyses demonstrate TLS-associated intratumoral B cell maturation in pathological responders, confirmed with spatial proteomics and BCR profiling. Our study also identifies spatial features of pathologic immune responses, revealing TLS maturation colocalizing with IgG/IgA distribution and extracellular matrix remodeling. HIGHLIGHTS: Integrated multi-modal spatial profiling of human PDAC tumors from neoadjuvant immunotherapy clinical trials reveal diverse spatial niches enriched in TLS.TLS maturity is influenced by tumor location and the cellular neighborhoods in which TLS immune cells are recruited.Unsupervised machine learning of genome-wide signatures on spatial transcriptomics data characterizes the TLS-enriched TME and associates TLS transcriptomes with survival outcomes in PDAC.Interactions of spatially variable gene expression patterns showed TLS maturation is coupled with immunoglobulin distribution and ECM remodeling in pathologic responders.Intratumoral plasma cell and immunoglobin gene expression spatial dynamics demonstrate trafficking of TLS-driven humoral immunity in the PDAC TME. Significance: We report a spatial multi-omics atlas of PDAC tumors from a series of immunotherapy neoadjuvant clinical trials. Intratumorally, pathologic responders exhibit mature TLS that propagate plasma cells into malignant niches. Our findings offer insights on the role of TLS-associated humoral immunity and stromal remodeling during immunotherapy treatment.

Patient-derived Organoid Pharmacotyping As A Predictive Tool for Therapeutic Selection in Pancreatic Ductal Adenocarcinoma.

OBJECTIVE: We integrate a new approach to chemosensitivity data for clinically-relevant regimen matching, and demonstrate the relationship with clinical outcomes in a large PDO biobank. SUMMARY BACKGROUND DATA: Pancreatic ductal adenocarcinoma (PDAC) usually recurs following potentially curative resection. Prior studies related patient-derived organoid (PDO) chemosensitivity with clinical responses. METHODS: PDOs were established from pre-treatment biopsies in a multi-institution clinical trial (n=21) and clinical specimens at a high-volume pancreatectomy center (n=74, of which 48 were pre-treated). PDO in vitro chemosensitivities to standard-of-care chemotherapeutics (pharmacotypes) were matched to potential clinically-relevant regimens by a weighted nearest-neighbors analysis. Clinical outcomes were then compared for patients who had well-matched versus poorly-matched treatment according to this metric. RESULTS: Our function matched 91% of PDOs to a standard-of-care regimen (9% pan-resistant). PDOs poorly-matched to the neoadjuvant regimen received would have matched to an alternative in 34% of cases. Patients receiving neoadjuvant chemotherapy well-matched to their pharmacotype experienced improved CA 19-9 response (60% decreased to normal when well-matched, 29% when poorly-matched, P<0.05) and lymph node down-staging (33% N0 after poorly-matched, 69% after well-matched, P<0.05). Patients receiving both well-matched neoadjuvant and adjuvant chemotherapy experienced improved recurrence-free- and overall survival (median RFS 8.5 mo poorly-matched, 15.9 mo well-matched, P<0.05; median OS 19.5 vs. 30.3 mo, P<0.05). CONCLUSION: In vitro PDO pharmacotyping can inform PDAC therapy selection. We demonstrate improved outcomes including survival for patients treated with regimens well-matched to their PDO chemosensitivities. A subsequent prospective study using PDO pharmacotype matching could improve oncologic outcomes and improve quality of life by avoiding therapies not expected to be effective.

Hypoxia induces ROS-resistant memory upon reoxygenation in vivo promoting metastasis in part via MUC1-C.

Hypoxia occurs in 90% of solid tumors and is associated with metastasis and mortality. Breast cancer cells that experience intratumoral hypoxia are 5x more likely to develop lung metastasis in animal models. Using spatial transcriptomics, we determine that hypoxic cells localized in more oxygenated tumor regions (termed 'post-hypoxic') retain expression of hypoxia-inducible and NF-kB-regulated genes, even in the oxygen-rich bloodstream. This cellular response is reproduced in vitro under chronic hypoxic conditions followed by reoxygenation. A subset of genes remains increased in reoxygenated cells. MUC1/MUC1-C is upregulated by both HIF-1α and NF-kB-p65 during chronic hypoxia. Abrogating MUC1 decreases the expression of superoxide dismutase enzymes, causing reactive oxygen species (ROS) production and cell death. A hypoxia-dependent genetic deletion of MUC1, or MUC1-C inhibition by GO-203, increases ROS levels in circulating tumor cells (CTCs), reducing the extent of metastasis. High MUC1 expression in tumor biopsies is associated with recurrence, and MUC1+ CTCs have lower ROS levels than MUC1- CTCs in patient-derived xenograft models. This study demonstrates that therapeutically targeting MUC1-C reduces hypoxia-driven metastasis.

Claudin 7 suppresses invasion and metastasis through repression of a smooth muscle actin program.

Metastasis initiates when cancer cells escape from the primary tumor, which requires changes to intercellular junctions. Claudins are transmembrane proteins that form the tight junction, and their expression is reduced in aggressive breast tumors. However, claudins' roles during breast cancer metastasis remain unclear. We used gain- and loss-of-function genetics in organoids isolated from murine breast cancer models to establish that Cldn7 suppresses invasion and metastasis. Transcriptomic analysis revealed that Cldn7 knockdown induced smooth muscle actin (SMA)-related genes and a broader mesenchymal phenotype. We validated our results in human cell lines, fresh human tumor tissue, bulk RNA-seq, and public single-cell RNA-seq data. We consistently observed an inverse relationship between Cldn7 expression and expression of SMA-related genes. Furthermore, knockdown and overexpression of SMA-related genes demonstrated that they promote breast cancer invasion. Our data reveal that Cldn7 suppresses breast cancer invasion and metastasis through negative regulation of SMA-related and mesenchymal gene expression.

Computational methods and biomarker discovery strategies for spatial proteomics: a review in immuno-oncology.

Advancements in imaging technologies have revolutionized our ability to deeply profile pathological tissue architectures, generating large volumes of imaging data with unparalleled spatial resolution. This type of data collection, namely, spatial proteomics, offers invaluable insights into various human diseases. Simultaneously, computational algorithms have evolved to manage the increasing dimensionality of spatial proteomics inherent in this progress. Numerous imaging-based computational frameworks, such as computational pathology, have been proposed for research and clinical applications. However, the development of these fields demands diverse domain expertise, creating barriers to their integration and further application. This review seeks to bridge this divide by presenting a comprehensive guideline. We consolidate prevailing computational methods and outline a roadmap from image processing to data-driven, statistics-informed biomarker discovery. Additionally, we explore future perspectives as the field moves toward interfacing with other quantitative domains, holding significant promise for precision care in immuno-oncology.

Identifying a gene signature of metastatic potential by linking pre-metastatic state to ultimate metastatic fate.

Identifying the key molecular pathways that enable metastasis by analyzing the eventual metastatic tumor is challenging because the state of the founder subclone likely changes following metastatic colonization. To address this challenge, we labeled primary mouse pancreatic ductal adenocarcinoma (PDAC) subclones with DNA barcodes to characterize their pre-metastatic state using ATAC-seq and RNA-seq and determine their relative in vivo metastatic potential prospectively. We identified a gene signature separating metastasis-high and metastasis-low subclones orthogonal to the normal-to-PDAC and classical-to-basal axes. The metastasis-high subclones feature activation of IL-1 pathway genes and high NF-κB and Zeb/Snail family activity and the metastasis-low subclones feature activation of neuroendocrine, motility, and Wnt pathway genes and high CDX2 and HOXA13 activity. In a functional screen, we validated novel mediators of PDAC metastasis in the IL-1 pathway, including the NF-κB targets Fos and Il23a, and beyond the IL-1 pathway including Myo1b and Tmem40. We scored human PDAC tumors for our signature of metastatic potential from mouse and found that metastases have higher scores than primary tumors. Moreover, primary tumors with higher scores are associated with worse prognosis. We also found that our metastatic potential signature is enriched in other human carcinomas, suggesting that it is conserved across epithelial malignancies. This work establishes a strategy for linking cancer cell state to future behavior, reveals novel functional regulators of PDAC metastasis, and establishes a method for scoring human carcinomas based on metastatic potential.

PanIN and CAF transitions in pancreatic carcinogenesis revealed with spatial data integration.

This study introduces a new imaging, spatial transcriptomics (ST), and single-cell RNA-sequencing integration pipeline to characterize neoplastic cell state transitions during tumorigenesis. We applied a semi-supervised analysis pipeline to examine premalignant pancreatic intraepithelial neoplasias (PanINs) that can develop into pancreatic ductal adenocarcinoma (PDAC). Their strict diagnosis on formalin-fixed and paraffin-embedded (FFPE) samples limited the single-cell characterization of human PanINs within their microenvironment. We leverage whole transcriptome FFPE ST to enable the study of a rare cohort of matched low-grade (LG) and high-grade (HG) PanIN lesions to track progression and map cellular phenotypes relative to single-cell PDAC datasets. We demonstrate that cancer-associated fibroblasts (CAFs), including antigen-presenting CAFs, are located close to PanINs. We further observed a transition from CAF-related inflammatory signaling to cellular proliferation during PanIN progression. We validate these findings with single-cell high-dimensional imaging proteomics and transcriptomics technologies. Altogether, our semi-supervised learning framework for spatial multi-omics has broad applicability across cancer types to decipher the spatiotemporal dynamics of carcinogenesis.

Individual variation in the emergence of anterior-to-posterior neural fates from human pluripotent stem cells.

Variability between human pluripotent stem cell (hPSC) lines remains a challenge and opportunity in biomedicine. In this study, hPSC lines from multiple donors were differentiated toward neuroectoderm and mesendoderm lineages. We revealed dynamic transcriptomic patterns that delineate the emergence of these lineages, which were conserved across lines, along with individual line-specific transcriptional signatures that were invariant throughout differentiation. These transcriptomic signatures predicted an antagonism between SOX21-driven forebrain fates and retinoic acid-induced hindbrain fates. Replicate lines and paired adult tissue demonstrated the stability of these line-specific transcriptomic traits. We show that this transcriptomic variation in lineage bias had both genetic and epigenetic origins, aligned with the anterior-to-posterior structure of early mammalian development, and was present across a large collection of hPSC lines. These findings contribute to developing systematic analyses of PSCs to define the origin and consequences of variation in the early events orchestrating individual human development.

Phosphodiesterase-5 inhibition collaborates with vaccine-based immunotherapy to reprogram myeloid cells in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is highly lethal and resistant to immunotherapy. Although immune recognition can be enhanced with immunomodulatory agents including checkpoint inhibitors and vaccines, few patients experience clinical efficacy because the tumor immune microenvironment (TiME) is dominated by immunosuppressive myeloid cells that impose T cell inhibition. Inhibition of phosphodiesterase-5 (PDE5) was reported to downregulate metabolic regulators arginase and inducible NOS in immunosuppressive myeloid cells and enhance immunity against immune-sensitive tumors, including head and neck cancers. We show for the first time to our knowledge that combining a PDE5 inhibitor, tadalafil, with a mesothelin-specific vaccine, anti-programmed cell death protein 1, and anti-cytotoxic T lymphocyte-associated protein 4 yields antitumor efficacy even against immune-resistant PDAC. To determine immunologic advantages conferred by tadalafil, we profiled the TiME using mass cytometry and single-cell RNA-sequencing analysis with Domino to infer intercellular signaling. Our analyses demonstrated that tadalafil reprograms myeloid cells to be less immunosuppressive. Moreover, tadalafil synergized with the vaccine, enhancing T cell activation including mesothelin-specific T cells. Tadalafil treatment was also associated with myeloid/T cell signaling axes important for antitumor responses (e.g., Cxcr3, Il12). Our study shows that PDE5 inhibition combined with vaccine-based immunotherapy promotes pro-inflammatory states of myeloid cells, activation of T cells, and enhanced myeloid/T cell crosstalk to yield antitumor efficacy against immune-resistant PDAC.

Is There a Mathematician on Board?

As convergence science is emerging as a theme in cancer research, scientists from diverse backgrounds, including mathematics, are increasingly entering our research community. The captain of a Southwest Airlines flight recently paged in jest for a mathematician to help support fuel calculations that would enable a flight to arrive safely at the American Association for Cancer Research (AACR) Annual Meeting, epitomizing the need for well-trained mathematicians to address pressing problems. Here, we summarize the roles mathematicians can play in cancer research and the support needed to facilitate their entry into the cancer research field. The inclusion of scientific diversity across quantitative and engineering disciplines is critical for advancing the understanding and treatment of cancer.

Cytokines drive the formation of memory-like NK cell subsets via epigenetic rewiring and transcriptional regulation.

Activation of natural killer (NK) cells with the cytokines interleukin-12 (IL-12), IL-15, and IL-18 induces their differentiation into memory-like (ML) NK cells; however, the underlying epigenetic and transcriptional mechanisms are unclear. By combining ATAC-seq, CITE-seq, and functional analyses, we discovered that IL-12/15/18 activation results in two main human NK fates: reprogramming into enriched memory-like (eML) NK cells or priming into effector conventional NK (effcNK) cells. eML NK cells had distinct transcriptional and epigenetic profiles and enhanced function, whereas effcNK cells resembled cytokine-primed cNK cells. Two transcriptionally discrete subsets of eML NK cells were also identified, eML-1 and eML-2, primarily arising from CD56bright or CD56dim mature NK cell subsets, respectively. Furthermore, these eML subsets were evident weeks after transfer of IL-12/15/18-activated NK cells into patients with cancer. Our findings demonstrate that NK cell activation with IL-12/15/18 results in previously unappreciated diverse cellular fates and identifies new strategies to enhance NK therapies.

Integration of Clinical Trial Spatial Multiomics Analysis and Virtual Clinical Trials Enables Immunotherapy Response Prediction and Biomarker Discovery.

Due to the lack of treatment options, there remains a need to advance new therapeutics in hepatocellular carcinoma (HCC). The traditional approach moves from initial molecular discovery through animal models to human trials to advance novel systemic therapies that improve treatment outcomes for patients with cancer. Computational methods that simulate tumors mathematically to describe cellular and molecular interactions are emerging as promising tools to simulate the impact of therapy entirely in silico, potentially greatly accelerating delivery of new therapeutics to patients. To facilitate the design of dosing regimens and identification of potential biomarkers for immunotherapy, we developed a new computational model to track tumor progression at the organ scale while capturing the spatial heterogeneity of the tumor in HCC. This computational model of spatial quantitative systems pharmacology was designed to simulate the effects of combination immunotherapy. The model was initiated using literature-derived parameter values and fitted to the specifics of HCC. Model validation was done through comparison with spatial multiomics data from a neoadjuvant HCC clinical trial combining anti-PD1 immunotherapy and a multitargeted tyrosine kinase inhibitor cabozantinib. Validation using spatial proteomics data from imaging mass cytometry demonstrated that closer proximity between CD8 T cells and macrophages correlated with nonresponse. We also compared the model output with Visium spatial transcriptomics profiling of samples from posttreatment tumor resections in the clinical trial and from another independent study of anti-PD1 monotherapy. Spatial transcriptomics data confirmed simulation results, suggesting the importance of spatial patterns of tumor vasculature and TGFß in tumor and immune cell interactions. Our findings demonstrate that incorporating mathematical modeling and computer simulations with high-throughput spatial multiomics data provides a novel approach for patient outcome prediction and biomarker discovery. Significance: Incorporating mathematical modeling and computer simulations with high-throughput spatial multiomics data provides an effective approach for patient outcome prediction and biomarker discovery.

Natural killer cells associate with epithelial cells in the pancreatic ductal adenocarcinoma tumor microenvironment.

Background: Pancreatic ductal adenocarcinoma (PDAC) is the most common form of pancreatic cancer. PDAC's poor prognosis and resistance to immunotherapy are attributed in part to its dense, fibrotic tumor microenvironment (TME), which is known to inhibit immune cell infiltration. We recently demonstrated that PDAC patients with higher natural killer (NK) cell content and activation have better survival rates. However, NK cell interactions in the PDAC TME have yet to be deeply studied. We show here that NK cells are present and active in the human PDAC TME. Methods: We used imaging mass cytometry (IMC) to assess NK cell content, function, and spatial localization in human PDAC samples. Then, we used CellChat, a tool to infer ligand-receptor interactions, on a human PDAC scRNAseq dataset to further define NK cell interactions in PDAC. Results: Spatial analyses showed for the first time that active NK cells are present in the PDAC TME, and both associate and interact with malignant epithelial cell ducts. We also found that fibroblast-rich, desmoplastic regions limit NK cell infiltration in the PDAC TME. CellChat analysis identified that the CD44 receptor on NK cells interacts with PDAC extracellular matrix (ECM) components such as collagen, fibronectin and laminin expressed by fibroblasts and malignant epithelial cells. This led us to hypothesize that these interactions play roles in regulating NK cell motility in desmoplastic PDAC TMEs. Using 2D and 3D in vitro assays, we found that CD44 neutralization significantly increased NK cell invasion through matrix. Conclusions: Targeting ECM-immune cell interactions may increase NK cell invasion into the PDAC TME.

Virtual patient analysis identifies strategies to improve the performance of predictive biomarkers for PD-1 blockade.

Patients with metastatic triple-negative breast cancer (TNBC) show variable responses to PD-1 inhibition. Efficient patient selection by predictive biomarkers would be desirable, but is hindered by the limited performance of existing biomarkers. Here, we leveraged in-silico patient cohorts generated using a quantitative systems pharmacology model of metastatic TNBC, informed by transcriptomic and clinical data, to explore potential ways to improve patient selection. We tested 90 biomarker candidates, including various cellular and molecular species, by a cutoff-based biomarker testing algorithm combined with machine learning-based feature selection. Combinations of pre-treatment biomarkers improved the specificity compared to single biomarkers at the cost of reduced sensitivity. On the other hand, early on-treatment biomarkers, such as the relative change in tumor diameter from baseline measured at two weeks after treatment initiation, achieved remarkably higher sensitivity and specificity. Further, blood-based biomarkers had a comparable ability to tumor- or lymph node-based biomarkers in identifying a subset of responders, potentially suggesting a less invasive way for patient selection.

Career trajectories of MD-PhD physician scientists: The loss of women investigators.

Advances in biomedical research require a robust physician scientist workforce. Despite being equally successful at securing early career awards from the NIH as men, women MD-PhD physician scientists are less likely to serve as principal investigators on mid- and later careers awards. Here, we discuss the causes of gender disparities in academic medicine, the implications of losing highly trained women physician scientists, and the institutional and systemic changes needed to sustain this pool of talented investigators.

Leveraging multi-omics data to empower quantitative systems pharmacology in immuno-oncology.

Understanding the intricate interactions of cancer cells with the tumor microenvironment (TME) is a pre-requisite for the optimization of immunotherapy. Mechanistic models such as quantitative systems pharmacology (QSP) provide insights into the TME dynamics and predict the efficacy of immunotherapy in virtual patient populations/digital twins but require vast amounts of multimodal data for parameterization. Large-scale datasets characterizing the TME are available due to recent advances in bioinformatics for multi-omics data. Here, we discuss the perspectives of leveraging omics-derived bioinformatics estimates to inform QSP models and circumvent the challenges of model calibration and validation in immuno-oncology.

Transfer Learning Reveals Cancer-Associated Fibroblasts Are Associated with Epithelial-Mesenchymal Transition and Inflammation in Cancer Cells in Pancreatic Ductal Adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive malignancy characterized by an immunosuppressive tumor microenvironment enriched with cancer-associated fibroblasts (CAF). This study used a convergence approach to identify tumor cell and CAF interactions through the integration of single-cell data from human tumors with human organoid coculture experiments. Analysis of a comprehensive atlas of PDAC single-cell RNA sequencing data indicated that CAF density is associated with increased inflammation and epithelial-mesenchymal transition (EMT) in epithelial cells. Transfer learning using transcriptional data from patient-derived organoid and CAF cocultures provided in silico validation of CAF induction of inflammatory and EMT epithelial cell states. Further experimental validation in cocultures demonstrated integrin beta 1 (ITGB1) and vascular endothelial factor A (VEGFA) interactions with neuropilin-1 mediating CAF-epithelial cell cross-talk. Together, this study introduces transfer learning from human single-cell data to organoid coculture analyses for experimental validation of discoveries of cell-cell cross-talk and identifies fibroblast-mediated regulation of EMT and inflammation. SIGNIFICANCE: Adaptation of transfer learning to relate human single-cell RNA sequencing data to organoid-CAF cocultures facilitates discovery of human pancreatic cancer intercellular interactions and uncovers cross-talk between CAFs and tumor cells through VEGFA and ITGB1.

Personalized neoantigen vaccine and pembrolizumab in advanced hepatocellular carcinoma: a phase 1/2 trial.

Programmed cell death protein 1 (PD-1) inhibitors have modest efficacy as a monotherapy in hepatocellular carcinoma (HCC). A personalized therapeutic cancer vaccine (PTCV) may enhance responses to PD-1 inhibitors through the induction of tumor-specific immunity. We present results from a single-arm, open-label, phase 1/2 study of a DNA plasmid PTCV (GNOS-PV02) encoding up to 40 neoantigens coadministered with plasmid-encoded interleukin-12 plus pembrolizumab in patients with advanced HCC previously treated with a multityrosine kinase inhibitor. Safety and immunogenicity were assessed as primary endpoints, and treatment efficacy and feasibility were evaluated as secondary endpoints. The most common treatment-related adverse events were injection-site reactions, observed in 15 of 36 (41.6%) patients. No dose-limiting toxicities or treatment-related grade ≥3 events were observed. The objective response rate (modified intention-to-treat) per Response Evaluation Criteria in Solid Tumors 1.1 was 30.6% (11 of 36 patients), with 8.3% (3 of 36) of patients achieving a complete response. Clinical responses were associated with the number of neoantigens encoded in the vaccine. Neoantigen-specific T cell responses were confirmed in 19 of 22 (86.4%) evaluable patients by enzyme-linked immunosorbent spot assays. Multiparametric cellular profiling revealed active, proliferative and cytolytic vaccine-specific CD4+ and CD8+ effector T cells. T cell receptor ß-chain (TCRß) bulk sequencing results demonstrated vaccination-enriched T cell clone expansion and tumor infiltration. Single-cell analysis revealed posttreatment T cell clonal expansion of cytotoxic T cell phenotypes. TCR complementarity-determining region cloning of expanded T cell clones in the tumors following vaccination confirmed reactivity against vaccine-encoded neoantigens. Our results support the PTCV's mechanism of action based on the induction of antitumor T cells and show that a PTCV plus pembrolizumab has clinical activity in advanced HCC. ClinicalTrials.gov identifier: NCT04251117 .