3D histology reveals that immune response to pancreatic precancers is heterogeneous and depends on global pancreas structure.

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal cancer for which few effective therapies exist. Immunotherapies specifically are ineffective in pancreatic cancer, in part due to its unique stromal and immune microenvironment. Pancreatic intraepithelial neoplasia, or PanIN, is the main precursor lesion to PDAC. Recently it was discovered that PanINs are remarkably abundant in the grossly normal pancreas, suggesting that the vast majority will never progress to cancer. Here, through construction of 48 samples of cm3-sized human pancreas tissue, we profiled the immune microenvironment of 1,476 PanINs in 3D and at single-cell resolution to better understand the early evolution of the pancreatic tumor microenvironment and to determine how inflammation may play a role in cancer progression. We found that bulk pancreatic inflammation strongly correlates to PanIN cell fraction. We found that the immune response around PanINs is highly heterogeneous, with distinct immune hotspots and cold spots that appear and disappear in a span of tens of microns. Immune hotspots generally mark locations of higher grade of dysplasia or locations near acinar atrophy. The immune composition at these hotspots is dominated by naïve, cytotoxic, and regulatory T cells, cancer associated fibroblasts, and tumor associated macrophages, with little similarity to the immune composition around less-inflamed PanINs. By mapping FOXP3+ cells in 3D, we found that regulatory T cells are present at higher density in larger PanIN lesions compared to smaller PanINs, suggesting that the early initiation of PanINs may not exhibit an immunosuppressive response. This analysis demonstrates that while PanINs are common in the pancreases of most individuals, inflammation may play a pivotal role, both at the bulk and the microscopic scale, in demarcating regions of significance in cancer progression.

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.

Age-Related Increases in IGFBP2 Increase Melanoma Cell Invasion and Lipid Synthesis.

Aged patients with melanoma (>65 years old) have more aggressive disease relative to young patients (<55 years old) for reasons that are not completely understood. Analysis of the young and aged secretome from human dermal fibroblasts identified >5-fold levels of IGF-binding protein 2 (IGFBP2) in the aged fibroblast secretome. IGFBP2 functionally triggers upregulation of the PI3K-dependent fatty acid biosynthesis program in melanoma cells. Melanoma cells co-cultured with aged dermal fibroblasts have higher levels of lipids relative to those co-cultured with young dermal fibroblasts, which can be lowered by silencing IGFBP2 expression in fibroblasts prior to treating with conditioned media. Conversely, ectopically treating melanoma cells with recombinant IGFBP2 in the presence of conditioned media from young fibroblasts or overexpressing IGFBP2 in melanoma cells promoted lipid synthesis and accumulation in melanoma cells. Treatment of young mice with rIGFBP2 increases tumor growth. Neutralizing IGFBP2 in vitro reduces migration and invasion in melanoma cells, and in vivo studies demonstrate that neutralizing IGFBP2 in syngeneic aged mice reduces tumor growth and metastasis. Our results suggest that aged dermal fibroblasts increase melanoma cell aggressiveness through increased secretion of IGFBP2, stressing the importance of considering age when designing studies and treatment. SIGNIFICANCE: The aged microenvironment drives metastasis in melanoma cells. This study reports that IGFBP2 secretion by aged fibroblasts induces lipid accumulation in melanoma cells, driving an increase in tumor invasiveness. Neutralizing IGFBP2 decreases melanoma tumor growth and metastasis.

E-cadherin Induces Serine Synthesis to Support Progression and Metastasis of Breast Cancer.

The loss of E-cadherin, an epithelial cell adhesion molecule, has been implicated in metastasis by mediating the epithelial-mesenchymal transition (EMT), which promotes invasion and migration of cancer cells. However, recent studies have demonstrated that E-cadherin supports the survival and proliferation of metastatic cancer cells. Here, we identified a metabolic role for E-cadherin in breast cancer by upregulating the de novo serine synthesis pathway (SSP). The upregulated SSP provided metabolic precursors for biosynthesis and resistance to oxidative stress, enabling E-cadherin+ breast cancer cells to achieve faster tumor growth and enhanced metastases. Inhibition of PHGDH, a rate-limiting enzyme in the SSP, significantly and specifically hampered proliferation of E-cadherin+ breast cancer cells and rendered them vulnerable to oxidative stress, inhibiting their metastatic potential. These findings reveal that E-cadherin reprograms cellular metabolism, promoting tumor growth and metastasis of breast cancers.

Power-law growth models explain incidences and sizes of pancreatic cancer precursor lesions and confirm spatial genomic findings.

Pancreatic ductal adenocarcinoma is a rare but lethal cancer. Recent evidence suggests that pancreatic intraepithelial neoplasia (PanIN), a microscopic precursor lesion that gives rise to pancreatic cancer, is larger and more prevalent than previously believed. Better understanding of the growth-law dynamics of PanINs may improve our ability to understand how a miniscule fraction makes the transition to invasive cancer. Here, using three-dimensional tissue mapping, we analyzed >1000 PanINs and found that lesion size is distributed according to a power law. Our data suggest that in bulk, PanIN size can be predicted by general growth behavior without consideration for the heterogeneity of the pancreatic microenvironment or an individual's age, history, or lifestyle. Our models suggest that intraductal spread and fusing of lesions drive our observed size distribution. This analysis lays the groundwork for future mathematical modeling efforts integrating PanIN incidence, morphology, and molecular features to understand tumorigenesis and demonstrates the utility of combining experimental measurement with dynamic modeling in understanding tumorigenesis.

High-motility pro-tumorigenic monocytes drive macrophage enrichment in the tumor microenvironment.

Enrichment of tumor-associated macrophages (TAMΦs) in the tumor microenvironment correlates with worse clinical outcomes in triple-negative breast cancer (TNBC) patients, prompting the development of therapies to inhibit TAMΦ infiltration. However, the lackluster efficacy of CCL2-based chemotaxis blockade in clinical trials suggests that a new understanding of monocyte/macrophage infiltration may be necessary. Here we demonstrate that random migration, and not only chemotaxis, drives macrophage tumor infiltration. We identified tumor- associated monocytes (TAMos) that display a dramatically enhanced migration capability, induced rapidly by the tumor microenvironment, that drives effective tumor infiltration, in contrast to low-motility differentiated macrophages. TAMo, not TAMΦ, promotes cancer cell proliferation through activation of the MAPK pathway. IL-6 secreted both by cancer cells and TAMo themselves enhances TAMo migration by increasing dendritic protrusion dynamics and myosin- based contractility via the JAK2/STAT3 signaling pathway. Independent from CCL2 mediated chemotaxis, IL-6 driven enhanced migration and pro-proliferative effect of TAMo were validated in a syngeneic TNBC mouse model. Depletion of IL-6 in cancer cells significantly attenuated monocyte infiltration and reversed TAMo-induced cancer cell proliferation. This work reveals the critical role random migration plays in monocyte driven TAMΦ enrichment in a tumor and pinpoints IL-6 as a potential therapeutic target in combination with CCL2 to ameliorate current strategies against TAMΦ infiltration.

qMAP enabled microanatomical mapping of human skin aging.

Aging is a major driver of diseases in humans. Identifying features associated with aging is essential for designing robust intervention strategies and discovering novel biomarkers of aging. Extensive studies at both the molecular and organ/whole-body physiological scales have helped determined features associated with aging. However, the lack of meso-scale studies, particularly at the tissue level, limits the ability to translate findings made at molecular scale to impaired tissue functions associated with aging. In this work, we established a tissue image analysis workflow - quantitative micro-anatomical phenotyping (qMAP) - that leverages deep learning and machine vision to fully label tissue and cellular compartments in tissue sections. The fully mapped tissue images address the challenges of finding an interpretable feature set to quantitatively profile age-related microanatomic changes. We optimized qMAP for skin tissues and applied it to a cohort of 99 donors aged 14 to 92. We extracted 914 microanatomic features and found that a broad spectrum of these features, represented by 10 cores processes, are strongly associated with aging. Our analysis shows that microanatomical features of the skin can predict aging with a mean absolute error (MAE) of 7.7 years, comparable to state-of-the-art epigenetic clocks. Our study demonstrates that tissue-level architectural changes are strongly associated with aging and represent a novel category of aging biomarkers that complement molecular markers. Our results highlight the complex and underexplored multi-scale relationship between molecular and tissue microanatomic scales.

Senescent cell heterogeneity and responses to senolytic treatment are related to cell cycle status during cell growth arrest.

Cellular senescence has been strongly linked to aging and age-related diseases. It is well established that the phenotype of senescent cells is highly heterogeneous and influenced by their cell type and senescence-inducing stimulus. Recent single-cell RNA-sequencing studies identified heterogeneity within senescent cell populations. However, proof of functional differences between such subpopulations is lacking. To identify functionally distinct senescent cell subpopulations, we employed high-content image analysis to measure senescence marker expression in primary human endothelial cells and fibroblasts. We found that G2-arrested senescent cells feature higher senescence marker expression than G1-arrested senescent cells. To investigate functional differences, we compared IL-6 secretion and response to ABT263 senolytic treatment in G1 and G2 senescent cells. We determined that G2-arrested senescent cells secrete more IL-6 and are more sensitive to ABT263 than G1-arrested cells. We hypothesize that cell cycle dependent DNA content is a key contributor to the heterogeneity within senescent cell populations. This study demonstrates the existence of functionally distinct senescent subpopulations even in culture. This data provides the first evidence of selective cell response to senolytic treatment among senescent cell subpopulations. Overall, this study emphasizes the importance of considering the senescent cell heterogeneity in the development of future senolytic therapies.

PanIN or IPMN? Redefining Lesion Size in 3 Dimensions.

Pancreatic ductal adenocarcinoma (PDAC) develops from 2 known precursor lesions: a majority (∼85%) develops from pancreatic intraepithelial neoplasia (PanIN), and a minority develops from intraductal papillary mucinous neoplasms (IPMNs). Clinical classification of PanIN and IPMN relies on a combination of low-resolution, 3-dimensional (D) imaging (computed tomography, CT), and high-resolution, 2D imaging (histology). The definitions of PanIN and IPMN currently rely heavily on size. IPMNs are defined as macroscopic: generally >1.0 cm and visible in CT, and PanINs are defined as microscopic: generally <0.5 cm and not identifiable in CT. As 2D evaluation fails to take into account 3D structures, we hypothesized that this classification would fail in evaluation of high-resolution, 3D images. To characterize the size and prevalence of PanINs in 3D, 47 thick slabs of pancreas were harvested from grossly normal areas of pancreatic resections, excluding samples from individuals with a diagnosis of an IPMN. All patients but one underwent preoperative CT scans. Through construction of cellular resolution 3D maps, we identified >1400 ductal precursor lesions that met the 2D histologic size criteria of PanINs. We show that, when 3D space is considered, 25 of these lesions can be digitally sectioned to meet the 2D histologic size criterion of IPMN. Re-evaluation of the preoperative CT images of individuals found to possess these large precursor lesions showed that nearly half are visible on imaging. These findings demonstrate that the clinical classification of PanIN and IPMN fails in evaluation of high-resolution, 3D images, emphasizing the need for re-evaluation of classification guidelines that place significant weight on 2D assessment of 3D structures.

3D genomic mapping reveals multifocality of human pancreatic precancers.

Pancreatic intraepithelial neoplasias (PanINs) are the most common precursors of pancreatic cancer, but their small size and inaccessibility in humans make them challenging to study1. Critically, the number, dimensions and connectivity of human PanINs remain largely unknown, precluding important insights into early cancer development. Here, we provide a microanatomical survey of human PanINs by analysing 46 large samples of grossly normal human pancreas with a machine-learning pipeline for quantitative 3D histological reconstruction at single-cell resolution. To elucidate genetic relationships between and within PanINs, we developed a workflow in which 3D modelling guides multi-region microdissection and targeted and whole-exome sequencing. From these samples, we calculated a mean burden of 13 PanINs per cm3 and extrapolated that the normal intact adult pancreas harbours hundreds of PanINs, almost all with oncogenic KRAS hotspot mutations. We found that most PanINs originate as independent clones with distinct somatic mutation profiles. Some spatially continuous PanINs were found to contain multiple KRAS mutations; computational and in situ analyses demonstrated that different KRAS mutations localize to distinct cell subpopulations within these neoplasms, indicating their polyclonal origins. The extensive multifocality and genetic heterogeneity of PanINs raises important questions about mechanisms that drive precancer initiation and confer differential progression risk in the human pancreas. This detailed 3D genomic mapping of molecular alterations in human PanINs provides an empirical foundation for early detection and rational interception of pancreatic cancer.

A 3D in vitro assay to study combined immune cell infiltration and cytotoxicity.

Immune cell-mediated killing of cancer cells in a solid tumor is prefaced by a multi-step infiltration cascade of invasion, directed migration, and cytotoxic activities. In particular, immune cells must invade and migrate through a series of different extracellular matrix (ECM) boundaries and domains before reaching and killing their target tumor cells. These infiltration events are a central challenge to the clinical success of CAR T cells against solid tumors. The current standard in vitro cell killing assays measure cell cytotoxicity in an obstacle-free, two-dimensional (2D) microenvironment, which precludes the study of 3D immune cell-ECM interactions. Here, we present a 3D combined infiltration/cytotoxicity assay based on an oil-in-water microtechnology. This assay measures stromal invasion following extravasation, migration through the stromal matrix, and invasion of the solid tumor in addition to cell killing. We compare this 3D cytotoxicity assay to the benchmark 2D assay through tumor assembloid cocultures with immune cells and engineered immune cells. This assay is amenable to an array of imaging techniques, which allows direct observation and quantification of each stage of infiltration in different immune and oncological contexts. We establish the 3D infiltration/cytotoxicity assay as an important tool for the mechanistic study of immune cell interactions with the tumor microenvironment.

Small Extracellular Vesicles Promote Stiffness-mediated Metastasis.

Tissue stiffness is a critical prognostic factor in breast cancer and is associated with metastatic progression. Here we show an alternative and complementary hypothesis of tumor progression whereby physiologic matrix stiffness affects the quantity and protein cargo of small extracellular vesicles (EV) produced by cancer cells, which in turn aid cancer cell dissemination. Primary patient breast tissue released by cancer cells on matrices that model human breast tumors (25 kPa; stiff EVs) feature increased adhesion molecule presentation (ITGα2ß1, ITGα6ß4, ITGα6ß1, CD44) compared with EVs from softer normal tissue (0.5 kPa; soft EVs), which facilitates their binding to extracellular matrix proteins including collagen IV, and a 3-fold increase in homing ability to distant organs in mice. In a zebrafish xenograft model, stiff EVs aid cancer cell dissemination. Moreover, normal, resident lung fibroblasts treated with stiff and soft EVs change their gene expression profiles to adopt a cancer-associated fibroblast phenotype. These findings show that EV quantity, cargo, and function depend heavily on the mechanical properties of the extracellular microenvironment. SIGNIFICANCE: Here we show that the quantity, cargo, and function of breast cancer-derived EVs vary with mechanical properties of the extracellular microenvironment.

MLL1 regulates cytokine-driven cell migration and metastasis.

Cell migration is a critical contributor to metastasis. Cytokine production and its role in cancer cell migration have been traditionally associated with immune cells. We find that the histone methyltransferase Mixed-Lineage Leukemia 1 (MLL1) controls 3D cell migration via cytokines, IL-6, IL-8, and TGF-ß1, secreted by the cancer cells themselves. MLL1, with its scaffold protein Menin, controls actin filament assembly via the IL-6/8/pSTAT3/Arp3 axis and myosin contractility via the TGF-ß1/Gli2/ROCK1/2/pMLC2 axis, which together regulate dynamic protrusion generation and 3D cell migration. MLL1 also regulates cell proliferation via mitosis-based and cell cycle-related pathways. Mice bearing orthotopic MLL1-depleted tumors exhibit decreased lung metastatic burden and longer survival. MLL1 depletion leads to lower metastatic burden even when controlling for the difference in primary tumor growth rates. Combining MLL1-Menin inhibitor with paclitaxel abrogates tumor growth and metastasis, including preexistent metastasis. These results establish MLL1 as a potent regulator of cell migration and highlight the potential of targeting MLL1 in patients with metastatic disease.

E-cadherin interacts with EGFR resulting in hyper-activation of ERK in multiple models of breast cancer.

The loss of intercellular adhesion molecule E-cadherin is a hallmark of the epithelial-mesenchymal transition (EMT), during which tumor cells transition into an invasive phenotype. Accordingly, E-cadherin has long been considered a tumor suppressor gene; however, E-cadherin expression is paradoxically correlated with breast cancer survival rates. Using novel multi-compartment organoids and multiple in vivo models, we show that E-cadherin promotes a hyper-proliferative phenotype in breast cancer cells via interaction with the transmembrane receptor EGFR. The E-cad and EGFR interaction results in activation of the MEK/ERK signaling pathway, leading to a significant increase in proliferation via activation of transcription factors, including c-Fos. Pharmacological inhibition of MEK activity in E-cadherin positive breast cancer significantly decreases both tumor growth and macro-metastasis in vivo. This work provides evidence for a novel role of E-cadherin in breast tumor progression and identifies a new target to treat hyper-proliferative E-cadherin-positive breast tumors, thus providing the foundation to utilize E-cadherin as a biomarker for specific therapeutic success.

Generative interpolation and restoration of images using deep learning for improved 3D tissue mapping.

The development of novel imaging platforms has improved our ability to collect and analyze large three-dimensional (3D) biological imaging datasets. Advances in computing have led to an ability to extract complex spatial information from these data, such as the composition, morphology, and interactions of multi-cellular structures, rare events, and integration of multi-modal features combining anatomical, molecular, and transcriptomic (among other) information. Yet, the accuracy of these quantitative results is intrinsically limited by the quality of the input images, which can contain missing or damaged regions, or can be of poor resolution due to mechanical, temporal, or financial constraints. In applications ranging from intact imaging (e.g. light-sheet microscopy and magnetic resonance imaging) to sectioning based platforms (e.g. serial histology and serial section transmission electron microscopy), the quality and resolution of imaging data has become paramount. Here, we address these challenges by leveraging frame interpolation for large image motion (FILM), a generative AI model originally developed for temporal interpolation, for spatial interpolation of a range of 3D image types. Comparative analysis demonstrates the superiority of FILM over traditional linear interpolation to produce functional synthetic images, due to its ability to better preserve biological information including microanatomical features and cell counts, as well as image quality, such as contrast, variance, and luminance. FILM repairs tissue damages in images and reduces stitching artifacts. We show that FILM can decrease imaging time by synthesizing skipped images. We demonstrate the versatility of our method with a wide range of imaging modalities (histology, tissue-clearing/light-sheet microscopy, magnetic resonance imaging, serial section transmission electron microscopy), species (human, mouse), healthy and diseased tissues (pancreas, lung, brain), staining techniques (IHC, H&E), and pixel resolutions (8 nm, 2 µm, 1mm). Overall, we demonstrate the potential of generative AI in improving the resolution, throughput, and quality of biological image datasets, enabling improved 3D imaging.

Substrate stiffness modulates the emergence and magnitude of senescence phenotypes in dermal fibroblasts.

Cellular senescence is a major driver of aging and disease. Here we show that substrate stiffness modulates the emergence and magnitude of senescence phenotypes after exposure to senescence inducers. Using a primary dermal fibroblast model, we show that decreased substrate stiffness accelerates senescence-associated cell-cycle arrest and regulates the expression of conventional protein-based biomarkers of senescence. We found that the expression of these senescence biomarkers, namely p21WAF1/CIP1 and p16INK4a are mechanosensitive and are in-part regulated by myosin contractility through focal adhesion kinase (FAK)-ROCK signaling. Interestingly, at the protein level senescence-induced dermal fibroblasts on soft substrates (0.5 kPa) do not express p21WAF1/CIP1 and p16INK4a at comparable levels to induced cells on stiff substrates (4GPa). However, cells express CDKN1a, CDKN2a, and IL6 at the RNA level across both stiff and soft substrates. Moreover, when cells are transferred from soft to stiff substrates, senescent cells recover an elevated expression of p21WAF1/CIP1 and p16INK4a at levels comparable to senescence cells on stiff substrates, pointing to a mechanosensitive regulation of the senescence phenotype. Together, our results indicate that the emergent senescence phenotype depends critically on the local mechanical environments of cells and that senescent cells actively respond to changing mechanical cues.

A disrupted compartment boundary underlies abnormal cardiac patterning and congenital heart defects.

Failure of septation of the interventricular septum (IVS) is the most common congenital heart defect (CHD), but mechanisms for patterning the IVS are largely unknown. We show that a Tbx5+/Mef2cAHF+ progenitor lineage forms a compartment boundary bisecting the IVS. This coordinated population originates at a first- and second heart field interface, subsequently forming a morphogenetic nexus. Ablation of Tbx5+/Mef2cAHF+ progenitors cause IVS disorganization, right ventricular hypoplasia and mixing of IVS lineages. Reduced dosage of the CHD transcription factor TBX5 disrupts boundary position and integrity, resulting in ventricular septation defects (VSDs) and patterning defects, including Slit2 and Ntn1 misexpression. Reducing NTN1 dosage partly rescues cardiac defects in Tbx5 mutant embryos. Loss of Slit2 or Ntn1 causes VSDs and perturbed septal lineage distributions. Thus, we identify essential cues that direct progenitors to pattern a compartment boundary for proper cardiac septation, revealing new mechanisms for cardiac birth defects.

Magnetic Resonance Imaging-Based Assessment of Pancreatic Fat Strongly Correlates With Histology-Based Assessment of Pancreas Composition.

OBJECTIVE: The aim of the study is to assess the relationship between magnetic resonance imaging (MRI)-based estimation of pancreatic fat and histology-based measurement of pancreatic composition. MATERIALS AND METHODS: In this retrospective study, MRI was used to noninvasively estimate pancreatic fat content in preoperative images from high-risk individuals and disease controls having normal pancreata. A deep learning algorithm was used to label 11 tissue components at micron resolution in subsequent pancreatectomy histology. A linear model was used to determine correlation between histologic tissue composition and MRI fat estimation. RESULTS: Twenty-seven patients (mean age 64.0 ± 12.0 years [standard deviation], 15 women) were evaluated. The fat content measured by MRI ranged from 0% to 36.9%. Intrapancreatic histologic tissue fat content ranged from 0.8% to 38.3%. MRI pancreatic fat estimation positively correlated with microanatomical composition of fat (r = 0.90, 0.83 to 0.95], P < 0.001); as well as with pancreatic cancer precursor ( r = 0.65, P < 0.001); and collagen ( r = 0.46, P < 0.001) content, and negatively correlated with pancreatic acinar ( r = -0.85, P < 0.001) content. CONCLUSIONS: Pancreatic fat content, measurable by MRI, correlates to acinar content, stromal content (fibrosis), and presence of neoplastic precursors of cancer.

Engineering self-propelled tumor-infiltrating CAR T cells using synthetic velocity receptors.

Chimeric antigen receptor (CAR) T cells express antigen-specific synthetic receptors, which upon binding to cancer cells, elicit T cell anti-tumor responses. CAR T cell therapy has enjoyed success in the clinic for hematological cancer indications, giving rise to decade-long remissions in some cases. However, CAR T therapy for patients with solid tumors has not seen similar success. Solid tumors constitute 90% of adult human cancers, representing an enormous unmet clinical need. Current approaches do not solve the central problem of limited ability of therapeutic cells to migrate through the stromal matrix. We discover that T cells at low and high density display low- and high-migration phenotypes, respectively. The highly migratory phenotype is mediated by a paracrine pathway from a group of self-produced cytokines that include IL5, TNFα, IFNγ, and IL8. We exploit this finding to "lock-in" a highly migratory phenotype by developing and expressing receptors, which we call velocity receptors (VRs). VRs target these cytokines and signal through these cytokines' cognate receptors to increase T cell motility and infiltrate lung, ovarian, and pancreatic tumors in large numbers and at doses for which control CAR T cells remain confined to the tumor periphery. In contrast to CAR therapy alone, VR-CAR T cells significantly attenuate tumor growth and extend overall survival. This work suggests that approaches to the design of immune cell receptors that focus on migration signaling will help current and future CAR cellular therapies to infiltrate deep into solid tumors.