Coordinated chemokine expression defines macrophage subsets across tissues.

Lung-resident macrophages, which include alveolar macrophages and interstitial macrophages (IMs), exhibit a high degree of diversity, generally attributed to different activation states, and often complicated by the influx of monocytes into the pool of tissue-resident macrophages. To gain a deeper insight into the functional diversity of IMs, here we perform comprehensive transcriptional profiling of resident IMs and reveal ten distinct chemokine-expressing IM subsets at steady state and during inflammation. Similar IM subsets that exhibited coordinated chemokine signatures and differentially expressed genes were observed across various tissues and species, indicating conserved specialized functional roles. Other macrophage types shared specific IM chemokine profiles, while also presenting their own unique chemokine signatures. Depletion of CD206hi IMs in Pf4creR26EYFP+DTR and Pf4creR26EYFPCx3cr1DTR mice led to diminished inflammatory cell recruitment, reduced tertiary lymphoid structure formation and fewer germinal center B cells in models of allergen- and infection-driven inflammation. These observations highlight the specialized roles of IMs, defined by their coordinated chemokine production, in regulating immune cell influx and organizing tertiary lymphoid tissue architecture.

Integrative Co-Registration of Elemental Imaging and Histopathology for Enhanced Spatial Multimodal Analysis of Tissue Sections through TRACE.

Summary: Elemental imaging provides detailed profiling of metal bioaccumulation, offering more precision than bulk analysis by targeting specific tissue areas. However, accurately identifying comparable tissue regions from elemental maps is challenging, requiring the integration of hematoxylin and eosin (H&E) slides for effective comparison. Facilitating the streamlined co-registration of Whole Slide Images (WSI) and elemental maps, TRACE enhances the analysis of tissue regions and elemental abundance in various pathological conditions. Through an interactive containerized web application, TRACE features real-time annotation editing, advanced statistical tools, and data export, supporting comprehensive spatial analysis. Notably, it allows for comparison of elemental abundances across annotated tissue structures and enables integration with other spatial data types through WSI co-registration. Availability and Implementation: Available on the following platforms- GitHub: jlevy44/trace_app , PyPI: trace_app , Docker: joshualevy44/trace_app , Singularity: joshualevy44/trace_app . Contact: joshua.levy@cshs.org. Supplementary information: Supplementary data are available.

Pharmacological induction of chromatin remodeling drives chemosensitization in triple-negative breast cancer.

Targeted therapies have improved outcomes for certain cancer subtypes, but cytotoxic chemotherapy remains a mainstay for triple-negative breast cancer (TNBC). The epithelial-to-mesenchymal transition (EMT) is a developmental program co-opted by cancer cells that promotes metastasis and chemoresistance. There are no therapeutic strategies specifically targeting mesenchymal-like cancer cells. We report that the US Food and Drug Administration (FDA)-approved chemotherapeutic eribulin induces ZEB1-SWI/SNF-directed chromatin remodeling to reverse EMT that curtails the metastatic propensity of TNBC preclinical models. Eribulin induces mesenchymal-to-epithelial transition (MET) in primary TNBC in patients, but conventional chemotherapy does not. In the treatment-naive setting, but not after acquired resistance to other agents, eribulin sensitizes TNBC cells to subsequent treatment with other chemotherapeutics. These findings provide an epigenetic mechanism of action of eribulin, supporting its use early in the disease process for MET induction to prevent metastatic progression and chemoresistance. These findings warrant prospective clinical evaluation of the chemosensitizing effects of eribulin in the treatment-naive setting.

Associations in cell type-specific hydroxymethylation and transcriptional alterations of pediatric central nervous system tumors.

Although intratumoral heterogeneity has been established in pediatric central nervous system tumors, epigenomic alterations at the cell type level have largely remained unresolved. To identify cell type-specific alterations to cytosine modifications in pediatric central nervous system tumors, we utilize a multi-omic approach that integrated bulk DNA cytosine modification data (methylation and hydroxymethylation) with both bulk and single-cell RNA-sequencing data. We demonstrate a large reduction in the scope of significantly differentially modified cytosines in tumors when accounting for tumor cell type composition. In the progenitor-like cell types of tumors, we identify a preponderance differential Cytosine-phosphate-Guanine site hydroxymethylation rather than methylation. Genes with differential hydroxymethylation, like histone deacetylase 4 and insulin-like growth factor 1 receptor, are associated with cell type-specific changes in gene expression in tumors. Our results highlight the importance of epigenomic alterations in the progenitor-like cell types and its role in cell type-specific transcriptional regulation in pediatric central nervous system tumors.

Hierarchal single-cell lineage tracing reveals differential fate commitment of CD8 T-cell clones in response to acute infection.

Generating balanced populations of CD8 effector and memory T cells is necessary for immediate and durable immunity to infections and cancer. Yet, a definitive understanding of CD8 differentiation remains unclear. We used CARLIN, a processive lineage recording mouse model with single-cell RNA-seq and TCR-seq to track endogenous antigen-specific CD8 T cells during acute viral infection. We identified a diverse repertoire of expanded T-cell clones represented by seven transcriptional states. TCR enrichment analysis revealed differential memory- or effector-fate biases within clonal populations. Shared Vb segments and amino acid motifs were found within biased categories despite high TCR diversity. Using single-cell CARLIN barcode-seq we tracked multi-generational clones and found that unlike unbiased or memory-biased clones, which stably retain their fate profiles, effector-biased clones could adopt memory- or effector-bias within subclones. Collectively, our study demonstrates that a heterogenous T-cell repertoire specific for a shared antigen is composed of clones with distinct TCR-intrinsic fate-biases.

Identifying tumor type and cell type-specific gene expression alterations in pediatric central nervous system tumors.

Central nervous system (CNS) tumors are the leading cause of pediatric cancer death, and these patients have an increased risk for developing secondary neoplasms. Due to the low prevalence of pediatric CNS tumors, major advances in targeted therapies have been lagging compared to other adult tumors. We collect single nuclei RNA-seq data from 84,700 nuclei of 35 pediatric CNS tumors and three non-tumoral pediatric brain tissues and characterize tumor heterogeneity and transcriptomic alterations. We distinguish cell subpopulations associated with specific tumor types including radial glial cells in ependymomas and oligodendrocyte precursor cells in astrocytomas. In tumors, we observe pathways important in neural stem cell-like populations, a cell type previously associated with therapy resistance. Lastly, we identify transcriptomic alterations among pediatric CNS tumor types compared to non-tumor tissues, while accounting for cell type effects on gene expression. Our results suggest potential tumor type and cell type-specific targets for pediatric CNS tumor treatment. Here we address current gaps in understanding single nuclei gene expression profiles of previously under-investigated tumor types and enhance current knowledge of gene expression profiles of single cells of various pediatric CNS tumors.

Potential to Enhance Large Scale Molecular Assessments of Skin Photoaging through Virtual Inference of Spatial Transcriptomics from Routine Staining.

The advent of spatial transcriptomics technologies has heralded a renaissance in research to advance our understanding of the spatial cellular and transcriptional heterogeneity within tissues. Spatial transcriptomics allows investigation of the interplay between cells, molecular pathways, and the surrounding tissue architecture and can help elucidate developmental trajectories, disease pathogenesis, and various niches in the tumor microenvironment. Photoaging is the histological and molecular skin damage resulting from chronic/acute sun exposure and is a major risk factor for skin cancer. Spatial transcriptomics technologies hold promise for improving the reliability of evaluating photoaging and developing new therapeutics. Challenges to current methods include limited focus on dermal elastosis variations and reliance on self-reported measures, which can introduce subjectivity and inconsistency. Spatial transcriptomics offers an opportunity to assess photoaging objectively and reproducibly in studies of carcinogenesis and discern the effectiveness of therapies that intervene in photoaging and preventing cancer. Evaluation of distinct histological architectures using highly-multiplexed spatial technologies can identify specific cell lineages that have been understudied due to their location beyond the depth of UV penetration. However, the cost and interpatient variability using state-of-the-art assays such as the 10x Genomics Spatial Transcriptomics assays limits the scope and scale of large-scale molecular epidemiologic studies. Here, we investigate the inference of spatial transcriptomics information from routine hematoxylin and eosin-stained (H&E) tissue slides. We employed the Visium CytAssist spatial transcriptomics assay to analyze over 18,000 genes at a 50-micron resolution for four patients from a cohort of 261 skin specimens collected adjacent to surgical resection sites for basal cell and squamous cell keratinocyte tumors. The spatial transcriptomics data was co-registered with 40x resolution whole slide imaging (WSI) information. We developed machine learning models that achieved a macro-averaged median AUC and F1 score of 0.80 and 0.61 and Spearman coefficient of 0.60 in inferring transcriptomic profiles across the slides, and accurately captured biological pathways across various tissue architectures.

Spatial Omics Driven Crossmodal Pretraining Applied to Graph-based Deep Learning for Cancer Pathology Analysis.

Graph-based deep learning has shown great promise in cancer histopathology image analysis by contextualizing complex morphology and structure across whole slide images to make high quality downstream outcome predictions (ex: prognostication). These methods rely on informative representations (i.e., embeddings) of image patches comprising larger slides, which are used as node attributes in slide graphs. Spatial omics data, including spatial transcriptomics, is a novel paradigm offering a wealth of detailed information. Pairing this data with corresponding histological imaging localized at 50-micron resolution, may facilitate the development of algorithms which better appreciate the morphological and molecular underpinnings of carcinogenesis. Here, we explore the utility of leveraging spatial transcriptomics data with a contrastive crossmodal pretraining mechanism to generate deep learning models that can extract molecular and histological information for graph-based learning tasks. Performance on cancer staging, lymph node metastasis prediction, survival prediction, and tissue clustering analyses indicate that the proposed methods bring improvement to graph based deep learning models for histopathological slides compared to leveraging histological information from existing schemes, demonstrating the promise of mining spatial omics data to enhance deep learning for pathology workflows.

Feasibility of Inferring Spatial Transcriptomics from Single-Cell Histological Patterns for Studying Colon Cancer Tumor Heterogeneity.

Background: Spatial transcriptomics involves studying the spatial organization of gene expression within tissues, offering insights into the molecular diversity of tumors. While spatial gene expression is commonly amalgamated from 1-10 cells across 50-micron spots, recent methods have demonstrated the capability to disaggregate this information at subspot resolution by leveraging both expression and histological patterns. However, elucidating such information from histology alone presents a significant challenge but if solved can better permit spatial molecular analysis at cellular resolution for instances where Visium data is not available, reducing study costs. This study explores integrating single-cell histological and transcriptomic data to infer spatial mRNA expression patterns in whole slide images collected from a cohort of stage pT3 colorectal cancer patients. A cell graph neural network algorithm was developed to align histological information extracted from detected cells with single cell RNA patterns through optimal transport methods, facilitating the analysis of cellular groupings and gene relationships. This approach leveraged spot-level expression as an intermediary to co-map histological and transcriptomic information at the single-cell level. Results: Our study demonstrated that single-cell transcriptional heterogeneity within a spot could be predicted from histological markers extracted from cells detected within a spot. Furthermore, our model exhibited proficiency in delineating overarching gene expression patterns across whole-slide images. This approach compared favorably to traditional patch-based computer vision methods as well as other methods which did not incorporate single cell expression during the model fitting procedures. Topological nuances of single-cell expression within a Visium spot were preserved using the developed methodology. Conclusion: This innovative approach augments the resolution of spatial molecular assays utilizing histology as a sole input through synergistic co-mapping of histological and transcriptomic datasets at the single-cell level, anchored by spatial transcriptomics. While initial results are promising, they warrant rigorous validation. This includes collaborating with pathologists for precise spatial identification of distinct cell types and utilizing sophisticated assays, such as Xenium, to attain deeper subcellular insights.

The Overlooked Role of Specimen Preparation in Bolstering Deep Learning-Enhanced Spatial Transcriptomics Workflows.

The application of deep learning methods to spatial transcriptomics has shown promise in unraveling the complex relationships between gene expression patterns and tissue architecture as they pertain to various pathological conditions. Deep learning methods that can infer gene expression patterns directly from tissue histomorphology can expand the capability to discern spatial molecular markers within tissue slides. However, current methods utilizing these techniques are plagued by substantial variability in tissue preparation and characteristics, which can hinder the broader adoption of these tools. Furthermore, training deep learning models using spatial transcriptomics on small study cohorts remains a costly endeavor. Necessitating novel tissue preparation processes enhance assay reliability, resolution, and scalability. This study investigated the impact of an enhanced specimen processing workflow for facilitating a deep learning-based spatial transcriptomics assessment. The enhanced workflow leveraged the flexibility of the Visium CytAssist assay to permit automated H&E staining (e.g., Leica Bond) of tissue slides, whole-slide imaging at 40x-resolution, and multiplexing of tissue sections from multiple patients within individual capture areas for spatial transcriptomics profiling. Using a cohort of thirteen pT3 stage colorectal cancer (CRC) patients, we compared the efficacy of deep learning models trained on slide prepared using an enhanced workflow as compared to the traditional workflow which leverages manual tissue staining and standard imaging of tissue slides. Leveraging Inceptionv3 neural networks, we aimed to predict gene expression patterns across matched serial tissue sections, each stemming from a distinct workflow but aligned based on persistent histological structures. Findings indicate that the enhanced workflow considerably outperformed the traditional spatial transcriptomics workflow. Gene expression profiles predicted from enhanced tissue slides also yielded expression patterns more topologically consistent with the ground truth. This led to enhanced statistical precision in pinpointing biomarkers associated with distinct spatial structures. These insights can potentially elevate diagnostic and prognostic biomarker detection by broadening the range of spatial molecular markers linked to metastasis and recurrence. Future endeavors will further explore these findings to enrich our comprehension of various diseases and uncover molecular pathways with greater nuance. Combining deep learning with spatial transcriptomics provides a compelling avenue to enrich our understanding of tumor biology and improve clinical outcomes. For results of the highest fidelity, however, effective specimen processing is crucial, and fostering collaboration between histotechnicians, pathologists, and genomics specialists is essential to herald this new era in spatial transcriptomics-driven cancer research.

Neoadjuvant Intratumoral Immunotherapy with Cowpea Mosaic Virus Induces Local and Systemic Antitumor Efficacy in Canine Mammary Cancer Patients.

The lack of optimal models to evaluate novel agents is delaying the development of effective immunotherapies against human breast cancer (BC). In this prospective open label study, we applied neoadjuvant intratumoral immunotherapy with empty cowpea mosaic virus-like particles (eCPMV) to 11 companion dogs diagnosed with canine mammary cancer (CMC), a spontaneous tumor resembling human BC. We found that two neoadjuvant intratumoral eCPMV injections resulted in tumor reduction in injected tumors in all patients and in noninjected tumors located in the ipsilateral and contralateral mammary chains of injected dogs. Tumor reduction was independent of clinical stage, tumor size, histopathologic grade, and tumor molecular subtype. RNA-seq-based analysis of injected tumors indicated a decrease in DNA replication activity and an increase in activated dendritic cell infiltration in the tumor microenvironment. Immunohistochemistry analysis demonstrated significant intratumoral increases in neutrophils, T and B lymphocytes, and plasma cells. eCPMV intratumoral immunotherapy demonstrated antitumor efficacy without any adverse effects. This novel immunotherapy has the potential for improving outcomes for human BC patients.

Potential to Enhance Large Scale Molecular Assessments of Skin Photoaging through Virtual Inference of Spatial Transcriptomics from Routine Staining.

The advent of spatial transcriptomics technologies has heralded a renaissance in research to advance our understanding of the spatial cellular and transcriptional heterogeneity within tissues. Spatial transcriptomics allows investigation of the interplay between cells, molecular pathways and the surrounding tissue architecture and can help elucidate developmental trajectories, disease pathogenesis, and various niches in the tumor microenvironment. Photoaging is the histological and molecular skin damage resulting from chronic/acute sun exposure and is a major risk factor for skin cancer. Spatial transcriptomics technologies hold promise for improving the reliability of evaluating photoaging and developing new therapeutics. Current challenges, including limited focus on dermal elastosis variations and reliance on self-reported measures, can introduce subjectivity and inconsistency. Spatial transcriptomics offer an opportunity to assess photoaging objectively and reproducibly in studies of carcinogenesis and discern the effectiveness of therapies that intervene on photoaging and prevent cancer. Evaluation of distinct histological architectures using highly-multiplexed spatial technologies can identify specific cell lineages that have been understudied due to their location beyond the depth of UV penetration. However, the cost and inter-patient variability using state-of-the-art assays such as the 10x Genomics Spatial Transcriptomics assays limits the scope and scale of large-scale molecular epidemiologic studies. Here, we investigate the inference of spatial transcriptomics information from routine hematoxylin and eosin-stained (H&E) tissue slides. We employed the Visium CytAssist spatial transcriptomics assay to analyze over 18,000 genes at a 50-micron resolution for four patients from a cohort of 261 skin specimens collected adjacent to surgical resection sites for basal and squamous keratinocyte tumors. The spatial transcriptomics data was co-registered with 40x resolution whole slide imaging (WSI) information. We developed machine learning models that achieved a macro-averaged median AUC and F1 score of 0.80 and 0.61 and Spearman coefficient of 0.60 in inferring transcriptomic profiles across the slides, and accurately captured biological pathways across various tissue architectures.

CDDO-Methyl Ester Inhibits BRAF Inhibitor Resistance and Remodels the Myeloid Compartment in BRAF-mutant Melanoma.

Approximately 50% of advanced melanomas harbor activating BRAF V600E mutations that are sensitive to BRAF inhibition. However, the duration of the response to BRAF inhibitors (BRAFi) has been limited due to the development of acquired resistance, which is preceded by recruitment of immunosuppressive myeloid cells and regulatory T cells (T regs ). While the addition of MAPK/ERK kinase 1 inhibitors (MEKi) prolongs therapeutic response to BRAF inhibition, most patients still develop resistance. Using a Braf V600E/+ /Pten -/- graft mouse model of melanoma, we now show that the addition of the methyl ester of the synthetic triterpenoid 2-cyano-3,12-dioxooleana-1,9(11)-dien-28-oic acid (C-Me) to the BRAFi vemurafenib analog PLX4720 at resistance significantly reduces tumor burden. Dual treatment remodels the BRAFi resistant-tumor microenvironment (TME), reducing infiltration of T regs and tumor associated macrophages (TAMs), and attenuates immunosuppressive cytokine production. For the first time, we characterize myeloid populations using scRNA-seq in BRAFi-resistant tumors and demonstrate that restoration of therapeutic response is associated with significant changes in immune-activated myeloid subset representation. Collectively, these studies suggest that C-Me inhibits acquired resistance to BRAFi. Use of C-Me in combination with other therapies may both inhibit melanoma growth and enhance therapeutic responsiveness more broadly.

Spatial Omics Driven Crossmodal Pretraining Applied to Graph-based Deep Learning for Cancer Pathology Analysis.

Graph-based deep learning has shown great promise in cancer histopathology image analysis by contextualizing complex morphology and structure across whole slide images to make high quality downstream outcome predictions (ex: prognostication). These methods rely on informative representations (i.e., embeddings) of image patches comprising larger slides, which are used as node attributes in slide graphs. Spatial omics data, including spatial transcriptomics, is a novel paradigm offering a wealth of detailed information. Pairing this data with corresponding histological imaging localized at 50-micron resolution, may facilitate the development of algorithms which better appreciate the morphological and molecular underpinnings of carcinogenesis. Here, we explore the utility of leveraging spatial transcriptomics data with a contrastive crossmodal pretraining mechanism to generate deep learning models that can extract molecular and histological information for graph-based learning tasks. Performance on cancer staging, lymph node metastasis prediction, survival prediction, and tissue clustering analyses indicate that the proposed methods bring improvement to graph based deep learning models for histopathological slides compared to leveraging histological information from existing schemes, demonstrating the promise of mining spatial omics data to enhance deep learning for pathology workflows.

Estrogen Therapy Induces Receptor-Dependent DNA Damage Enhanced by PARP Inhibition in ER+ Breast Cancer.

PURPOSE: Clinical evidence indicates that treatment with estrogens elicits anticancer effects in ∼30% of patients with advanced endocrine-resistant estrogen receptor α (ER)-positive breast cancer. Despite the proven efficacy of estrogen therapy, its mechanism of action is unclear and this treatment remains underused. Mechanistic understanding may offer strategies to enhance therapeutic efficacy. EXPERIMENTAL DESIGN: We performed genome-wide CRISPR/Cas9 screening and transcriptomic profiling in long-term estrogen-deprived ER+ breast cancer cells to identify pathways required for therapeutic response to the estrogen 17ß-estradiol (E2). We validated findings in cell lines, patient-derived xenografts (PDX), and patient samples, and developed a novel combination treatment through testing in cell lines and PDX models. RESULTS: Cells treated with E2 exhibited replication-dependent markers of DNA damage and the DNA damage response prior to apoptosis. Such DNA damage was partially driven by the formation of DNA:RNA hybrids (R-loops). Pharmacologic suppression of the DNA damage response via PARP inhibition with olaparib enhanced E2-induced DNA damage. PARP inhibition synergized with E2 to suppress growth and prevent tumor recurrence in BRCA1/2-mutant and BRCA1/2-wild-type cell line and PDX models. CONCLUSIONS: E2-induced ER activity drives DNA damage and growth inhibition in endocrine-resistant breast cancer cells. Inhibition of the DNA damage response using drugs such as PARP inhibitors can enhance therapeutic response to E2. These findings warrant clinical exploration of the combination of E2 with DNA damage response inhibitors in advanced ER+ breast cancer, and suggest that PARP inhibitors may synergize with therapeutics that exacerbate transcriptional stress.

Characterizing control of memory CD8 T cell differentiation by BTB-ZF transcription factor Zbtb20.

Members of the BTB-ZF transcription factor family regulate the immune system. Our laboratory identified that family member Zbtb20 contributes to the differentiation, recall responses, and metabolism of CD8 T cells. Here, we report a characterization of the transcriptional and epigenetic signatures controlled by Zbtb20 at single-cell resolution during the effector and memory phases of the CD8 T cell response. Without Zbtb20, transcriptional programs associated with memory CD8 T cell formation were up-regulated throughout the CD8 T response. A signature of open chromatin was associated with genes controlling T cell activation, consistent with the known impact on differentiation. In addition, memory CD8 T cells lacking Zbtb20 were characterized by open chromatin regions with overrepresentation of AP-1 transcription factor motifs and elevated RNA- and protein-level expressions of the corresponding AP-1 components. Finally, we describe motifs and genomic annotations from the DNA targets of Zbtb20 in CD8 T cells identified by cleavage under targets and release under nuclease (CUT&RUN). Together, these data establish the transcriptional and epigenetic networks contributing to the control of CD8 T cell responses by Zbtb20.

Alternating 17ß-Estradiol and Aromatase Inhibitor Therapies Is Efficacious in Postmenopausal Women with Advanced Endocrine-Resistant ER+ Breast Cancer.

PURPOSE: Strategies to implement estrogen therapy for advanced estrogen receptor-positive (ER+) breast cancer are underdeveloped. Preclinical data suggest that cycling treatment with 17ß-estradiol followed by estrogen deprivation can control tumor growth long-term. PATIENTS AND METHODS: Postmenopausal women with advanced ER+/HER2- breast cancer with recurrence or progression on ≥ 1 antiestrogen or aromatase inhibitor (AI)-based therapy were eligible. Patients received 17ß-estradiol (2 mg orally, three times a day) for 8 weeks followed by AI (physician's choice) for 16 weeks, alternating treatments on an 8-week/16-week schedule until disease progression. Patients then optionally received continuous single-agent treatment until a second instance of disease progression. Endpoints included 24-week clinical benefit and objective response per RECIST, and tumor genetic alterations. RESULTS: Of 19 evaluable patients, clinical benefit rate was 42.1% [95% confidence interval (CI), 23.1%-63.9%] and objective response rate (ORR) was 15.8% (95% CI, 5.7%-37.9%). One patient experienced a grade 3 adverse event related to 17ß-estradiol. Among patients who received continuous single-agent treatment until a second instance of disease progression, clinical benefit was observed in 5 of 12 (41.7%) cases. Tumor ER (ESR1) mutations were found by whole-exome profiling in 4 of 7 (57.1%) versus 2 of 9 (22.2%) patients who did versus did not experience clinical benefit from alternating 17ß-estradiol/AI therapy. The only two patients to experience objective responses to initial 17ß-estradiol had tumor ESR1 mutations. CONCLUSIONS: Alternating 17ß-estradiol/AI therapy may be a promising treatment for endocrine-refractory ER+ breast cancer, including following progression on CDK4/6 inhibitors or everolimus. Further study is warranted to determine whether the antitumor activity of 17ß-estradiol differs according to ESR1 mutation status.

Pharmacological Induction of mesenchymal-epithelial transition chemosensitizes breast cancer cells and prevents metastatic progression.

The epithelial-mesenchymal transition (EMT) is a developmental program co-opted by tumor cells that aids the initiation of the metastatic cascade. Tumor cells that undergo EMT are relatively chemoresistant, and there are currently no therapeutic avenues specifically targeting cells that have acquired mesenchymal traits. We show that treatment of mesenchymal-like triple-negative breast cancer (TNBC) cells with the microtubule-destabilizing chemotherapeutic eribulin, which is FDA-approved for the treatment of advanced breast cancer, leads to a mesenchymal-epithelial transition (MET). This MET is accompanied by loss of metastatic propensity and sensitization to subsequent treatment with other FDA-approved chemotherapeutics. We uncover a novel epigenetic mechanism of action that supports eribulin pretreatment as a path to MET induction that curtails metastatic progression and the evolution of therapy resistance.

Inferring spatial transcriptomics markers from whole slide images to characterize metastasis-related spatial heterogeneity of colorectal tumors: A pilot study.

Over 150 000 Americans are diagnosed with colorectal cancer (CRC) every year, and annually over 50 000 individuals will die from CRC, necessitating improvements in screening, prognostication, disease management, and therapeutic options. Tumor metastasis is the primary factor related to the risk of recurrence and mortality. Yet, screening for nodal and distant metastasis is costly, and invasive and incomplete resection may hamper adequate assessment. Signatures of the tumor-immune microenvironment (TIME) at the primary site can provide valuable insights into the aggressiveness of the tumor and the effectiveness of various treatment options. Spatially resolved transcriptomics technologies offer an unprecedented characterization of TIME through high multiplexing, yet their scope is constrained by cost. Meanwhile, it has long been suspected that histological, cytological, and macroarchitectural tissue characteristics correlate well with molecular information (e.g., gene expression). Thus, a method for predicting transcriptomics data through inference of RNA patterns from whole slide images (WSI) is a key step in studying metastasis at scale. In this work, we collected tissue from 4 stage-III (pT3) matched colorectal cancer patients for spatial transcriptomics profiling. The Visium spatial transcriptomics (ST) assay was used to measure transcript abundance for 17 943 genes at up to 5000 55-micron (i.e., 1-10 cells) spots per patient sampled in a honeycomb pattern, co-registered with hematoxylin and eosin (H&E) stained WSI. The Visium ST assay can measure expression at these spots through tissue permeabilization of mRNAs, which are captured through spatially (i.e., x-y positional coordinates) barcoded, gene specific oligo probes. WSI subimages were extracted around each co-registered Visium spot and were used to predict the expression at these spots using machine learning models. We prototyped and compared several convolutional, transformer, and graph convolutional neural networks to predict spatial RNA patterns at the Visium spots under the hypothesis that the transformer- and graph-based approaches better capture relevant spatial tissue architecture. We further analyzed the model's ability to recapitulate spatial autocorrelation statistics using SPARK and SpatialDE. Overall, the results indicate that the transformer- and graph-based approaches were unable to outperform the convolutional neural network architecture, though they exhibited optimal performance for relevant disease-associated genes. Initial findings suggest that different neural networks that operate on different scales are relevant for capturing distinct disease pathways (e.g., epithelial to mesenchymal transition). We add further evidence that deep learning models can accurately predict gene expression in whole slide images and comment on understudied factors which may increase its external applicability (e.g., tissue context). Our preliminary work will motivate further investigation of inference for molecular patterns from whole slide images as metastasis predictors and in other applications.

Identification of Spatial Proteomic Signatures of Colon Tumor Metastasis: A Digital Spatial Profiling Approach.

Over 150,000 Americans are diagnosed with colorectal cancer (CRC) every year, and annually >50,000 individuals are estimated to die of CRC, necessitating improvements in screening, prognostication, disease management, and therapeutic options. CRC tumors are removed en bloc with surrounding vasculature and lymphatics. Examination of regional lymph nodes at the time of surgical resection is essential for prognostication. Developing alternative approaches to indirectly assess recurrence risk would have utility in cases where lymph node yield is incomplete or inadequate. Spatially dependent, immune cell-specific (eg, tumor-infiltrating lymphocytes), proteomic, and transcriptomic expression patterns inside and around the tumor-the tumor immune microenvironment-can predict nodal/distant metastasis and probe the coordinated immune response from the primary tumor site. The comprehensive characterization of tumor-infiltrating lymphocytes and other immune infiltrates is possible using highly multiplexed spatial omics technologies, such as the GeoMX Digital Spatial Profiler. In this study, machine learning and differential co-expression analyses helped identify biomarkers from Digital Spatial Profiler-assayed protein expression patterns inside, at the invasive margin, and away from the tumor, associated with extracellular matrix remodeling (eg, granzyme B and fibronectin), immune suppression (eg, forkhead box P3), exhaustion and cytotoxicity (eg, CD8), Programmed death ligand 1-expressing dendritic cells, and neutrophil proliferation, among other concomitant alterations. Further investigation of these biomarkers may reveal independent risk factors of CRC metastasis that can be formulated into low-cost, widely available assays.