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.

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.

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.

Cytoplasmic accumulation and plasma membrane association of anillin and Ect2 promote confined migration and invasion.

Cells migrating in confinement experience mechanical challenges whose consequences on cell migration machinery remain only partially understood. Here, we demonstrate that a pool of the cytokinesis regulatory protein anillin is retained during interphase in the cytoplasm of different cell types. Confinement induces recruitment of cytoplasmic anillin to plasma membrane at the poles of migrating cells, which is further enhanced upon nuclear envelope (NE) rupture(s). Rupture events also enable the cytoplasmic egress of predominantly nuclear RhoGEF Ect2. Anillin and Ect2 redistributions scale with microenvironmental stiffness and confinement, and are observed in confined cells in vitro and in invading tumor cells in vivo. Anillin, which binds actomyosin at the cell poles, and Ect2, which activates RhoA, cooperate additively to promote myosin II contractility, and promote efficient invasion and extravasation. Overall, our work provides a mechanistic understanding of how cytokinesis regulators mediate RhoA/ROCK/myosin II-dependent mechanoadaptation during confined migration and invasive cancer progression.

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.

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.

Precision-engineered biomimetics: the human fallopian tube.

The fallopian tube has an essential role in several physiological and pathological processes from pregnancy to ovarian cancer. However, there are no biologically relevant models to study its pathophysiology. The state-of-the-art organoid model has been compared to two-dimensional tissue sections and molecularly assessed providing only cursory analyses of the model's accuracy. We developed a novel multi-compartment organoid model of the human fallopian tube that was meticulously tuned to reflect the compartmentalization and heterogeneity of the tissue's composition. We validated this organoid's molecular expression patterns, cilia-driven transport function, and structural accuracy through a highly iterative platform wherein organoids are compared to a three-dimensional, single-cell resolution reference map of a healthy, transplantation-quality human fallopian tube. This organoid model was precision-engineered to match the human microanatomy. One sentence summary: Tunable organoid modeling and CODA architectural quantification in tandem help design a tissue-validated organoid model.

Uncovering the spatial landscape of molecular interactions within the tumor microenvironment through latent spaces.

Recent advances in spatial transcriptomics (STs) enable gene expression measurements from a tissue sample while retaining its spatial context. This technology enables unprecedented in situ resolution of the regulatory pathways that underlie the heterogeneity in the tumor as well as the tumor microenvironment (TME). The direct characterization of cellular co-localization with spatial technologies facilities quantification of the molecular changes resulting from direct cell-cell interaction, as it occurs in tumor-immune interactions. We present SpaceMarkers, a bioinformatics algorithm to infer molecular changes from cell-cell interactions from latent space analysis of ST data. We apply this approach to infer the molecular changes from tumor-immune interactions in Visium spatial transcriptomics data of metastasis, invasive and precursor lesions, and immunotherapy treatment. Further transfer learning in matched scRNA-seq data enabled further quantification of the specific cell types in which SpaceMarkers are enriched. Altogether, SpaceMarkers can identify the location and context-specific molecular interactions within the TME from ST data.

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.

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.

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.

Lung Cancer Recurrence Risk Prediction through Integrated Deep Learning Evaluation.

Background: Prognostic risk factors for completely resected stage IA non-small-cell lung cancers (NSCLCs) have advanced minimally over recent decades. Although several biomarkers have been found to be associated with cancer recurrence, their added value to TNM staging and tumor grade are unclear. Methods: Features of preoperative low-dose CT image and histologic findings of hematoxylin- and eosin-stained tissue sections of resected lung tumor specimens were extracted from 182 stage IA NSCLC patients in the National Lung Screening Trial. These features were combined to predict the risk of tumor recurrence or progression through integrated deep learning evaluation (IDLE). Added values of IDLE to TNM staging and tumor grade in progression risk prediction and risk stratification were evaluated. Results: The 5-year AUC of IDLE was 0.817 ± 0.037 as compared to the AUC = 0.561 ± 0.042 and 0.573 ± 0.044 from the TNM stage and tumor grade, respectively. The IDLE score was significantly associated with cancer recurrence (p < 0.0001) even after adjusting for TNM staging and tumor grade. Synergy between chest CT image markers and histological markers was the driving force of the deep learning algorithm to produce a stronger prognostic predictor. Conclusions: Integrating markers from preoperative CT images and pathologist's readings of resected lung specimens through deep learning can improve risk stratification of stage 1A NSCLC patients over TNM staging and tumor grade alone. Our study suggests that combining markers from nonoverlapping platforms can increase the cancer risk prediction accuracy.

Deep learning identification of stiffness markers in breast cancer.

While essential to our understanding of solid tumor progression, the study of cell and tissue mechanics has yet to find traction in the clinic. Determining tissue stiffness, a mechanical property known to promote a malignant phenotype in vitro and in vivo, is not part of the standard algorithm for the diagnosis and treatment of breast cancer. Instead, clinicians routinely use mammograms to identify malignant lesions and radiographically dense breast tissue is associated with an increased risk of developing cancer. Whether breast density is related to tumor tissue stiffness, and what cellular and non-cellular components of the tumor contribute the most to its stiffness are not well understood. Through training of a deep learning network and mechanical measurements of fresh patient tissue, we create a bridge in understanding between clinical and mechanical markers. The automatic identification of cellular and extracellular features from hematoxylin and eosin (H&E)-stained slides reveals that global and local breast tissue stiffness best correlate with the percentage of straight collagen. Importantly, the percentage of dense breast tissue does not directly correlate with tissue stiffness or straight collagen content.

Cell Trafficking at the Intersection of the Tumor-Immune Compartments.

Migration is an essential cellular process that regulates human organ development and homeostasis as well as disease initiation and progression. In cancer, immune and tumor cell migration is strongly associated with immune cell infiltration, immune escape, and tumor cell metastasis, which ultimately account for more than 90% of cancer deaths. The biophysics and molecular regulation of the migration of cancer and immune cells have been extensively studied separately. However, accumulating evidence indicates that, in the tumor microenvironment, the motilities of immune and cancer cells are highly interdependent via secreted factors such as cytokines and chemokines. Tumor and immune cells constantly express these soluble factors, which produce a tightly intertwined regulatory network for these cells' respective migration. A mechanistic understanding of the reciprocal regulation of soluble factor-mediated cell migration can provide critical information for the development of new biomarkers of tumor progression and of tumor response to immuno-oncological treatments. We review the biophysical andbiomolecular basis for the migration of immune and tumor cells and their associated reciprocal regulatory network. We also describe ongoing attempts to translate this knowledge into the clinic.

Algal extracellular organic matter mediated photocatalytic degradation of estrogens.

Estrogens are among the most concerned emerging contaminants in the wastewater treatment effluent due to their sexual disruption in aquatic wildlife. The use of microalgae for secondary wastewater effluent polishing is a promising approach due to the economic benefit and value-added products. In this study, three microalgae species, including Selenastrum capricornutum, Scenedesmus quadricauda and Chlorella vulgaris were selected to conduct batch experiments to examine important mechanisms, especially the role of algal extracellular organic matter (AEOM) on two selected estrogens (17ß-estradiol, E2 and 17α-ethynylestradiol, EE2) removal. Results showed that estrogens could not be significantly degraded under visible light irradiation and adsorption of estrogens by microalgae was negligible. All three living microalgae cultures have ability to remove E2 and EE2, and Selenastrum capricornutum showed the highest E2 and EE2 removal efficiency of 91% and 83%, respectively, corresponding to the reduction of predicted estrogenic activity of 86%. AEOM from three microalgae cultures could induce photodegradation of estrogens, and AEOM from Selenastrum capricornutum and Chlorella vulgaris achieved 100% of E2 and EE2 removal under visible light irradiation. Fluorescence excitation-emission matrix spectroscopy identified humic/fulvic-like substances in AEOM from three microalgae cultures, which might be responsible for inducing the indirect photolysis of E2 and EE2. Therefore, in the living microalgae cultures, the major estrogens removal mechanisms should include biotransformation as well as AEOM meditated photocatalytic degradation. Since removal rates through photodegradation could be faster than biotransformation, the AEOM mediated photocatalytic degradation can play a potential role to remove emerging contaminants when using microalgae technology for wastewater effluent treatment.

A microfluidic cell-migration assay for the prediction of progression-free survival and recurrence time of patients with glioblastoma.

Clinical scores, molecular markers and cellular phenotypes have been used to predict the clinical outcomes of patients with glioblastoma. However, their clinical use has been hampered by confounders such as patient co-morbidities, by the tumoral heterogeneity of molecular and cellular markers, and by the complexity and cost of high-throughput single-cell analysis. Here, we show that a microfluidic assay for the quantification of cell migration and proliferation can categorize patients with glioblastoma according to progression-free survival. We quantified with a composite score the ability of primary glioblastoma cells to proliferate (via the protein biomarker Ki-67) and to squeeze through microfluidic channels, mimicking aspects of the tight perivascular conduits and white-matter tracts in brain parenchyma. The assay retrospectively categorized 28 patients according to progression-free survival (short-term or long-term) with an accuracy of 86%, predicted time to recurrence and correctly categorized five additional patients on the basis of survival prospectively. RNA sequencing of the highly motile cells revealed differentially expressed genes that correlated with poor prognosis. Our findings suggest that cell-migration and proliferation levels can predict patient-specific clinical outcomes.