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.

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 cm 3 -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.

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.

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.

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.

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.

Three-dimensional assessments are necessary to determine the true, spatially-resolved composition of tissues.

Methods for spatially resolved cellular profiling using thinly cut sections have enabled in-depth quantitative tissue mapping to study inter-sample and intra-sample differences in normal human anatomy and disease onset and progression. These methods often profile extremely limited regions, which may impact the evaluation of heterogeneity due to tissue sub-sampling. Here, we applied CODA, a deep learning-based tissue mapping platform, to reconstruct the three-dimensional (3D) microanatomy of grossly normal and cancer-containing human pancreas biospecimens obtained from individuals who underwent pancreatic resection. To compare inter- and intra-sample heterogeneity, we assessed bulk and spatially resolved tissue composition in a cohort of two-dimensional (2D) whole slide images (WSIs) and a cohort of thick slabs of pancreas tissue that were digitally reconstructed in 3D from serial sections. To demonstrate the marked under sampling of 2D assessments, we simulated the number of WSIs and tissue microarrays (TMAs) necessary to represent the compositional heterogeneity of 3D data within 10% error to reveal that tens of WSIs and hundreds of TMA cores are sometimes needed. We show that spatial correlation of different pancreatic structures decay significantly within a span of microns, demonstrating that 2D histological sections may not be representative of their neighboring tissues. In sum, we demonstrate that 3D assessments are necessary to accurately assess tissue composition in normal and abnormal specimens and in order to accurately determine neoplastic content. These results emphasize the importance of intra-sample heterogeneity in tissue mapping efforts.

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 reveals that pancreatic intraepithelial neoplasms (PanINs), the microscopic precursor lesions in the pancreatic ducts that can give rise to invasive pancreatic cancer, are significantly 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 of these lesions makes the transition to invasive cancer. Here, using artificial intelligence (AI)-based three-dimensional (3D) tissue mapping method, we measured the volumes of >1,000 PanIN and found that lesion size is distributed according to a power law with a fitted exponent of -1.7 over > 3 orders of magnitude. Our data also suggest that PanIN growth is not very sensitive to the pancreatic microenvironment or an individual's age, family history, and lifestyle, and is rather shaped by general growth behavior. We analyze several models of PanIN growth and fit the predicted size distributions to the observed data. The best fitting models suggest that both intraductal spread of PanIN lesions and fusing of multiple lesions into large, highly branched structures drive PanIN growth patterns. This work lays the groundwork for future mathematical modeling efforts integrating PanIN incidence, morphology, genomic, and transcriptomic features to understand pancreas tumorigenesis, and demonstrates the utility of combining experimental measurement of human tissues with dynamic modeling for understanding cancer tumorigenesis.

Digitize your Biology! Modeling multicellular systems through interpretable cell behavior.

Cells are fundamental units of life, constantly interacting and evolving as dynamical systems. While recent spatial multi-omics can quantitate individual cells' characteristics and regulatory programs, forecasting their evolution ultimately requires mathematical modeling. We develop a conceptual framework-a cell behavior hypothesis grammar-that uses natural language statements (cell rules) to create mathematical models. This allows us to systematically integrate biological knowledge and multi-omics data to make them computable. We can then perform virtual "thought experiments" that challenge and extend our understanding of multicellular systems, and ultimately generate new testable hypotheses. In this paper, we motivate and describe the grammar, provide a reference implementation, and demonstrate its potential through a series of examples in tumor biology and immunotherapy. Altogether, this approach provides a bridge between biological, clinical, and systems biology researchers for mathematical modeling of biological systems at scale, allowing the community to extrapolate from single-cell characterization to emergent multicellular behavior.

Tissue clearing and 3D reconstruction of digitized, serially sectioned slides provide novel insights into pancreatic cancer.

Pancreatic cancer is currently the third leading cause of cancer death in the United States. The clinical hallmarks of this disease include abdominal pain that radiates to the back, the presence of a hypoenhancing intrapancreatic lesion on imaging, and widespread liver metastases. Technologies such as tissue clearing and three-dimensional (3D) reconstruction of digitized serially sectioned hematoxylin and eosin-stained slides can be used to visualize large (up to 2- to 3-centimeter cube) tissues at cellular resolution. When applied to human pancreatic cancers, these 3D visualization techniques have provided novel insights into the basis of a number of the clinical characteristics of this disease. Here, we describe the clinical features of pancreatic cancer, review techniques for clearing and the 3D reconstruction of digitized microscope slides, and provide examples that illustrate how 3D visualization of human pancreatic cancer at the microscopic level has revealed features not apparent in 2D microscopy and, in so doing, has closed the gap between bench and bedside. Compared with animal models and 2D microscopy, studies of human tissues in 3D can reveal the difference between what can happen and what does happen in human cancers.

Three-dimensional genomic mapping of human pancreatic tissue reveals striking multifocality and genetic heterogeneity in precancerous lesions.

Pancreatic intraepithelial neoplasia (PanIN) is a precursor to pancreatic cancer and represents a critical opportunity for cancer interception. However, the number, size, shape, and connectivity of PanINs in human pancreatic tissue samples are largely unknown. In this study, we quantitatively assessed human PanINs using CODA, a novel machine-learning pipeline for 3D image analysis that generates quantifiable models of large pieces of human pancreas with single-cell resolution. Using a cohort of 38 large slabs of grossly normal human pancreas from surgical resection specimens, we identified striking multifocality of PanINs, with a mean burden of 13 spatially separate PanINs per cm3 of sampled tissue. Extrapolating this burden to the entire pancreas suggested a median of approximately 1000 PanINs in an entire pancreas. In order to better understand the clonal relationships within and between PanINs, we developed a pipeline for CODA-guided multi-region genomic analysis of PanINs, including targeted and whole exome sequencing. Multi-region assessment of 37 PanINs from eight additional human pancreatic tissue slabs revealed that almost all PanINs contained hotspot mutations in the oncogene KRAS, but no gene other than KRAS was altered in more than 20% of the analyzed PanINs. PanINs contained a mean of 13 somatic mutations per region when analyzed by whole exome sequencing. The majority of analyzed PanINs originated from independent clonal events, with distinct somatic mutation profiles between PanINs in the same tissue slab. A subset of the analyzed PanINs contained multiple KRAS mutations, suggesting a polyclonal origin even in PanINs that are contiguous by rigorous 3D assessment. This study leverages a novel 3D genomic mapping approach to describe, for the first time, the spatial and genetic multifocality of human PanINs, providing important insights into the initiation and progression of pancreatic neoplasia.

Intraparenchymal metastases as a cause for local recurrence of pancreatic cancer.

Using CODA, a technique for three-dimensional reconstruction of large tissues, Kiemen et al. report observation of a microscopic focus of pancreatic cancer found in the vasculature of grossly normal human pancreas tissue resected adjacent to a large tumour. They use TP53 and SMAD4 staining to relate the small focus to the primary tumour. This report describes a represents a probable case of intraparenchymal metastasis of pancreatic cancer, revealing a probable cause of local recurrence.

CODA: quantitative 3D reconstruction of large tissues at cellular resolution.

A central challenge in biology is obtaining high-content, high-resolution information while analyzing tissue samples at volumes relevant to disease progression. We address this here with CODA, a method to reconstruct exceptionally large (up to multicentimeter cubed) tissues at subcellular resolution using serially sectioned hematoxylin and eosin-stained tissue sections. Here we demonstrate CODA's ability to reconstruct three-dimensional (3D) distinct microanatomical structures in pancreas, skin, lung and liver tissues. CODA allows creation of readily quantifiable tissue volumes amenable to biological research. As a testbed, we assess the microanatomy of the human pancreas during tumorigenesis within the branching pancreatic ductal system, labeling ten distinct structures to examine heterogeneity and structural transformation during neoplastic progression. We show that pancreatic precancerous lesions develop into distinct 3D morphological phenotypes and that pancreatic cancer tends to spread far from the bulk tumor along collagen fibers that are highly aligned to the 3D curves of ductal, lobular, vascular and neural structures. Thus, CODA establishes a means to transform broadly the structural study of human diseases through exploration of exhaustively labeled 3D microarchitecture.

Characterization of tumor-associated macrophages in prostate cancer transgenic mouse models.

BACKGROUND: Tumor-associated macrophages (TAMs) are critical components of the tumor microenvironment (TME) in prostate cancer. Commonly used orthotopic models do not accurately reflect the complete TME of a human patient or the natural initiation and progression of a tumor. Therefore, genetically engineered mouse models are essential for studying the TME as well as advancing TAM-targeted therapies. Two common transgenic (TG) models of prostate cancer are Hi-Myc and transgenic adenocarcinoma of the mouse prostate (TRAMP), but the TME and TAM characteristics of these models have not been well characterized. METHODS: To advance the Hi-Myc and TRAMP models as tools for TAM studies, macrophage infiltration and characteristics were assessed using histopathologic, flow cytometric, and expression analyses in these models at various timepoints during tumor development and progression. RESULTS: In both Hi-Myc and TRAMP models, macrophages adopt a more pro-tumor phenotype in higher histological grade tumors and in older prostate tissue. However, the Hi-Myc and TRAMP prostates differ in their macrophage density, with Hi-Myc tumors exhibiting increased macrophage density and TRAMP tumors exhibiting decreased macrophage density compared to age-matched wild type mice. CONCLUSIONS: The macrophage density and the adenocarcinoma cancer subtype of Hi-Myc appear to better mirror patient tumors, suggesting that the Hi-Myc model is the more appropriate in vivo TG model for studying TAMs and TME-targeted therapies.