Junctional Adhesion Molecule (JAM)-C recruitment of Pard3 and drebrin to cell contacts initiates neuron-glia recognition and layer-specific cell sorting in developing cerebella.

Sorting maturing neurons into distinct layers is critical for brain development, with disruptions leading to neurological disorders and pediatric cancers. Lamination coordinates where, when, and how cells interact, facilitating events that direct migrating neurons to their destined positions within emerging neural networks and control the wiring of connections in functional circuits. While the role of adhesion molecule expression and presentation in driving adhesive recognition during neuronal migration along glial fibers is recognized, the mechanisms by which the spatial arrangement of these molecules on the cell surface dictates adhesive specificity and translates contact-based external cues into intracellular responses like polarization and cytoskeletal organization remain largely unexplored. We used the cerebellar granule neuron (CGN) system to demonstrate that JAM-C receptor cis-binding on the same cell and trans-binding to neighboring cells controls the recruitment of the Pard3 polarity protein and drebrin microtubule-actin crosslinker at CGN to glial adhesion sites, complementing previous studies that showed Pard3 controls JAM-C exocytic surface presentation. Leveraging advanced imaging techniques, specific probes for cell recognition, and analytical methods to dissect adhesion dynamics, our findings reveal: 1) JAM-C cis or trans mutants result in reduced adhesion formation between CGNs and cerebellar glia, 2) these mutants exhibit delayed recruitment of Pard3 at the adhesion sites, and 3) CGNs with JAM-C mutations experience postponed sorting and entry into the cerebellar molecular layer (ML). By developing a conditional system to image adhesion components from two different cells simultaneously, we made it possible to investigate the dynamics of cell recognition on both sides of neuron-glial contacts and the subsequent recruitment of proteins required for CGN migration. This system and an approach that calculates local correlation based on convolution kernels at the cell adhesions site revealed that CGN to CGN JAM recognition preferentially recruits higher levels of Pard3 and drebrin than CGN to glia JAM recognition. The long latency time of CGNs in the inner external germinal layer (EGL) can be attributed to the combined strength of CGN-CGN contacts and the less efficient Pard3 recruitment by CGN-BG contacts, acting as gatekeepers to ML entry. As CGNs eventually transition to glia binding for radial migration, our research demonstrates that establishing permissive JAM-recognition sites on glia via cis and trans interactions of CGN JAM-C serves as a critical temporal checkpoint for sorting at the EGL to ML boundary. This mechanism integrates intrinsic and extrinsic cellular signals, facilitating heterotypic cell sorting into the ML and dictating the precise spatial organization within the cerebellar architecture.

A spatial cell atlas of neuroblastoma reveals developmental, epigenetic and spatial axis of tumor heterogeneity.

Neuroblastoma is a pediatric cancer arising from the developing sympathoadrenal lineage with complex inter- and intra-tumoral heterogeneity. To chart this complexity, we generated a comprehensive cell atlas of 55 neuroblastoma patient tumors, collected from two pediatric cancer institutions, spanning a range of clinical, genetic, and histologic features. Our atlas combines single-cell/nucleus RNA-seq (sc/scRNA-seq), bulk RNA-seq, whole exome sequencing, DNA methylation profiling, spatial transcriptomics, and two spatial proteomic methods. Sc/snRNA-seq revealed three malignant cell states with features of sympathoadrenal lineage development. All of the neuroblastomas had malignant cells that resembled sympathoblasts and the more differentiated adrenergic cells. A subset of tumors had malignant cells in a mesenchymal cell state with molecular features of Schwann cell precursors. DNA methylation profiles defined four groupings of patients, which differ in the degree of malignant cell heterogeneity and clinical outcomes. Using spatial proteomics, we found that neuroblastomas are spatially compartmentalized, with malignant tumor cells sequestered away from immune cells. Finally, we identify spatially restricted signaling patterns in immune cells from spatial transcriptomics. To facilitate the visualization and analysis of our atlas as a resource for further research in neuroblastoma, single cell, and spatial-omics, all data are shared through the Human Tumor Atlas Network Data Commons at www.humantumoratlas.org.

Progressive development of melanoma-induced cachexia differentially impacts organ systems in mice.

Cachexia is a systemic wasting syndrome that increases cancer-associated mortality. How cachexia progressively and differentially impacts distinct tissues is largely unknown. Here, we find that the heart and skeletal muscle undergo wasting at early stages and are the tissues transcriptionally most impacted by cachexia. We also identify general and organ-specific transcriptional changes that indicate functional derangement by cachexia even in tissues that do not undergo wasting, such as the brain. Secreted factors constitute a top category of cancer-regulated genes in host tissues, and these changes include upregulation of the angiotensin-converting enzyme (ACE). ACE inhibition with the drug lisinopril improves muscle force and partially impedes cachexia-induced transcriptional changes, although wasting is not prevented, suggesting that cancer-induced host-secreted factors can regulate tissue function during cachexia. Altogether, by defining prevalent and temporal and tissue-specific responses to cachexia, this resource highlights biomarkers and possible targets for general and tissue-tailored anti-cachexia therapies.

Matricellular Protein WISP2 Is an Endogenous Inhibitor of Collagen Linearization and Cancer Metastasis.

Collagen remodeling contributes to many physiologic and pathologic processes. In primary tumors, the linearization of collagen fibers promotes cancer cell invasion and metastasis and is indicative of poor prognosis. However, it remains unknown whether there are endogenous inhibitors of collagen linearization that could be exploited therapeutically. Here, we show that collagen linearization is controlled by two secreted matricellular proteins with antagonistic functions. Specifically, WISP1 was secreted by cancer cells, bound to type I collagen (Col I), and linearized Col I via its cysteine-rich C-terminal (CT) domain. In contrast, WISP2, which lacks a CT domain, inhibited Col I linearization by preventing WISP1-Col I binding. Analysis of patient data revealed that WISP2 expression is lower in most solid tumors, in comparison with normal tissues. Consequently, genetic or pharmacologic restoration of higher WISP2 levels impaired collagen linearization and prevented tumor cell invasion and metastasis in vivo in models of human and murine breast cancer. Thus, this study uncovers WISP2 as the first inhibitor of collagen linearization ever identified and reveals that collagen architecture can be normalized and metastasis inhibited by therapeutically restoring a high WISP2:WISP1 ratio. SIGNIFICANCE: Two secreted factors, WISP1 and WISP2, antagonistically regulate collagen linearization, and therapeutically increasing the WISP2:WISP1 ratio in tumors limits collagen linearization and inhibits metastasis.See related commentary by Barcus and Longmore, p. 5611.

Oxygen Tension and the VHL-Hif1α Pathway Determine Onset of Neuronal Polarization and Cerebellar Germinal Zone Exit.

Postnatal brain circuit assembly is driven by temporally regulated intrinsic and cell-extrinsic cues that organize neurogenesis, migration, and axo-dendritic specification in post-mitotic neurons. While cell polarity is an intrinsic organizer of morphogenic events, environmental cues in the germinal zone (GZ) instructing neuron polarization and their coupling during postnatal development are unclear. We report that oxygen tension, which rises at birth, and the von Hippel-Lindau (VHL)-hypoxia-inducible factor 1α (Hif1α) pathway regulate polarization and maturation of post-mitotic cerebellar granule neurons (CGNs). At early postnatal stages with low GZ vascularization, Hif1α restrains CGN-progenitor cell-cycle exit. Unexpectedly, cell-intrinsic VHL-Hif1α pathway activation also delays the timing of CGN differentiation, germinal zone exit, and migration initiation through transcriptional repression of the partitioning-defective (Pard) complex. As vascularization proceeds, these inhibitory mechanisms are downregulated, implicating increasing oxygen tension as a critical switch for neuronal polarization and cerebellar GZ exit.

Correlative three-dimensional super-resolution and block-face electron microscopy of whole vitreously frozen cells.

Within cells, the spatial compartmentalization of thousands of distinct proteins serves a multitude of diverse biochemical needs. Correlative super-resolution (SR) fluorescence and electron microscopy (EM) can elucidate protein spatial relationships to global ultrastructure, but has suffered from tradeoffs of structure preservation, fluorescence retention, resolution, and field of view. We developed a platform for three-dimensional cryogenic SR and focused ion beam-milled block-face EM across entire vitreously frozen cells. The approach preserves ultrastructure while enabling independent SR and EM workflow optimization. We discovered unexpected protein-ultrastructure relationships in mammalian cells including intranuclear vesicles containing endoplasmic reticulum-associated proteins, web-like adhesions between cultured neurons, and chromatin domains subclassified on the basis of transcriptional activity. Our findings illustrate the value of a comprehensive multimodal view of ultrastructural variability across whole cells.

Deep multiomics profiling of brain tumors identifies signaling networks downstream of cancer driver genes.

High throughput omics approaches provide an unprecedented opportunity for dissecting molecular mechanisms in cancer biology. Here we present deep profiling of whole proteome, phosphoproteome and transcriptome in two high-grade glioma (HGG) mouse models driven by mutated RTK oncogenes, PDGFRA and NTRK1, analyzing 13,860 proteins and 30,431 phosphosites by mass spectrometry. Systems biology approaches identify numerous master regulators, including 41 kinases and 23 transcription factors. Pathway activity computation and mouse survival indicate the NTRK1 mutation induces a higher activation of AKT downstream targets including MYC and JUN, drives a positive feedback loop to up-regulate multiple other RTKs, and confers higher oncogenic potency than the PDGFRA mutation. A mini-gRNA library CRISPR-Cas9 validation screening shows 56% of tested master regulators are important for the viability of NTRK-driven HGG cells, including TFs (Myc and Jun) and metabolic kinases (AMPKa1 and AMPKa2), confirming the validity of the multiomics integrative approaches, and providing novel tumor vulnerabilities.

Interleukin-15 Enhances Anti-GD2 Antibody-Mediated Cytotoxicity in an Orthotopic PDX Model of Neuroblastoma.

PURPOSE: Immunotherapy with IL2, GM-CSF, and an anti-disialoganglioside (GD2) antibody significantly increases event-free survival in children with high-risk neuroblastoma. However, therapy failure in one third of these patients and IL2-related toxicities pose a major challenge. We compared the immunoadjuvant effects of IL15 with those of IL2 for enhancing antibody-dependent cell-mediated cytotoxicity (ADCC) in neuroblastoma. EXPERIMENTAL DESIGN: We tested ADCC against neuroblastoma patient-derived xenografts (PDX) in vitro and in vivo and examined the functional and migratory properties of NK cells activated with IL2 and IL15. RESULTS: In cell culture, IL15-activated NK cells induced higher ADCC against two GD+ neuroblastoma PDXs than did IL2-activated NK cells (P < 0.001). This effect was dose-dependent (P < 0.001) and was maintained across several effector-to-tumor ratios. As compared with IL2, IL15 also improved chemotaxis of NK cells, leading to higher numbers of tumorsphere-infiltrating NK cells in vitro (P = 0.002). In an orthotopic PDX model, animals receiving chemoimmunotherapy with an anti-GD2 antibody, GM-CSF, and a soluble IL15/IL15Rα complex had greater tumor regression than did those receiving chemotherapy alone (P = 0.012) or combined with anti-GD2 antibody and GM-CSF with (P = 0.016) or without IL2 (P = 0.035). This was most likely due to lower numbers of immature tumor-infiltrating NK cells (DX5+CD27+) after IL15/IL15Rα administration (P = 0.029) and transcriptional upregulation of Gzmd. CONCLUSIONS: The substitution of IL15 for IL2 leads to significant tumor regression in vitro and in vivo and supports clinical testing of IL15 for immunotherapy in pediatric neuroblastoma.

The tumor cell-secreted matricellular protein WISP1 drives pro-metastatic collagen linearization.

Collagen linearization is a hallmark of aggressive tumors and a key pathogenic event that promotes cancer cell invasion and metastasis. Cell-generated mechanical tension has been proposed to contribute to collagen linearization in tumors, but it is unknown whether other mechanisms play prominent roles in this process. Here, we show that the secretome of cancer cells is by itself able to induce collagen linearization independently of cell-generated mechanical forces. Among the tumor cell-secreted factors, we find a key role in this process for the matricellular protein WISP1 (CCN4). Specifically, WISP1 directly binds to type I collagen to promote its linearization in vitro (in the absence of cells) and in vivo in tumors. Consequently, WISP1-induced type I collagen linearization facilitates tumor cell invasion and promotes spontaneous breast cancer metastasis, without significantly affecting gene expression. Furthermore, higher WISP1 expression in tumors from cancer patients correlates with faster progression to metastatic disease and poor prognosis. Altogether, these findings reveal a conceptually novel mechanism whereby pro-metastatic collagen linearization critically depends on a cancer cell-secreted factor.

The Dynamic Epigenetic Landscape of the Retina During Development, Reprogramming, and Tumorigenesis.

In the developing retina, multipotent neural progenitors undergo unidirectional differentiation in a precise spatiotemporal order. Here we profile the epigenetic and transcriptional changes that occur during retinogenesis in mice and humans. Although some progenitor genes and cell cycle genes were epigenetically silenced during retinogenesis, the most dramatic change was derepression of cell-type-specific differentiation programs. We identified developmental-stage-specific super-enhancers and showed that most epigenetic changes are conserved in humans and mice. To determine how the epigenome changes during tumorigenesis and reprogramming, we performed integrated epigenetic analysis of murine and human retinoblastomas and induced pluripotent stem cells (iPSCs) derived from murine rod photoreceptors. The retinoblastoma epigenome mapped to the developmental stage when retinal progenitors switch from neurogenic to terminal patterns of cell division. The epigenome of retinoblastomas was more similar to that of the normal retina than that of retina-derived iPSCs, and we identified retina-specific epigenetic memory.

Measuring Absolute Blood Perfusion in Mice Using Dynamic Contrast-Enhanced Ultrasound.

We investigated the feasibility of estimating absolute tissue blood perfusion using dynamic contrast-enhanced ultrasound (CEUS) imaging in mice. We developed a novel method of microbubble administration and a model-free approach to estimate absolute kidney perfusion, and explored the kidney as a reference organ to estimate absolute perfusion of a neuroblastoma tumor. We performed CEUS on the kidneys of CD1 nude mice using the VisualSonics VEVO 2100 imaging system. We estimated individual kidney blood perfusion using the burst-replenishment (BR) technique. We repeated the kidney imaging on the mice after a week. We performed CEUS imaging of a neuroblastoma mouse xenograft tumor along with its right kidney using two sets of microbubble administration parameters to estimate absolute tumor blood perfusion. We performed statistical tests at a significance level of 0.05. Our estimated absolute kidney perfusion (425 ± 123 mL/min/100 g) was within the range of previously reported values. There was no statistical difference between the estimated absolute kidney blood perfusions from the 2 wk of imaging (paired t-test, p = 0.09). We estimated the absolute blood perfusion in the neuroblastoma tumor to be 16.49 and 16.9 mL/min/100 g for the two sets of microbubble administration parameters (Wilcoxon rank-sum test, p = 0.6). We have established the kidney as a reliable reference organ in which to estimate absolute perfusion of other tissues. Using a neuroblastoma tumor, we have determined the feasibility of estimating absolute blood perfusion in tissues using contrast-enhanced ultrasound imaging.

Detection of Phenotypic Alterations Using High-Content Analysis of Whole-Slide Images.

Tumors exhibit spatial heterogeneity, as manifested in immunohistochemistry (IHC) staining patterns. Current IHC quantification methods lose information by reducing this heterogeneity in each whole-slide image (WSI) or in selective fields of view to a single staining index. The aim of this study was to investigate the sensitivity of an IHC quantification method that uses this heterogeneity to reliably compare IHC staining patterns. We virtually partitioned WSIs by a grid of square tiles, and computed the staining index distributions to quantify heterogeneities. We used samples from these distributions as inputs to non-parametric statistical comparisons. We applied our grid method to fixed tumor samples from 26 tumors obtained from a double-blind preclinical study of a patient-derived orthotopic xenograft model of pediatric neuroblastoma in CD1 nude mice. We compared the results of our grid method to the results based on whole-slide indices, the current practice. We show that our grid method reliably detects phenotypic alterations that other tests based on whole-slide indices fail to detect. Based on robustness and increased sensitivity of statistical inference, we conclude that our method of whole-slide grid quantification is superior to existing whole-slide quantification techniques.

Emergent Stratification in Solid Tumors Selects for Reduced Cohesion of Tumor Cells: A Multi-Cell, Virtual-Tissue Model of Tumor Evolution Using CompuCell3D.

Tumor cells and structure both evolve due to heritable variation of cell behaviors and selection over periods of weeks to years (somatic evolution). Micro-environmental factors exert selection pressures on tumor-cell behaviors, which influence both the rate and direction of evolution of specific behaviors, especially the development of tumor-cell aggression and resistance to chemotherapies. In this paper, we present, step-by-step, the development of a multi-cell, virtual-tissue model of tumor somatic evolution, simulated using the open-source CompuCell3D modeling environment. Our model includes essential cell behaviors, microenvironmental components and their interactions. Our model provides a platform for exploring selection pressures leading to the evolution of tumor-cell aggression, showing that emergent stratification into regions with different cell survival rates drives the evolution of less cohesive cells with lower levels of cadherins and higher levels of integrins. Such reduced cohesivity is a key hallmark in the progression of many types of solid tumors.

Development and characterization of a human orthotopic neuroblastoma xenograft.

Neuroblastoma is a pediatric cancer of the developing sympathoadrenal lineage. The tumors are known to develop from the adrenal gland or paraspinal ganglia and have molecular and cellular features of sympathetic neurons such as dense core vesicles and catecholamine production. Here we present the detailed molecular, cellular, genetic and epigenetic characterization of an orthotopic xenograft derived from a high-risk stage 4 neuroblastoma patient. Overall, the xenografted tumor retained the high risk features of the primary tumor and showed aggressive growth and metastasis in the mouse. Also, the genome was preserved with no additional copy number variations, structural variations or aneuploidy. There were 13 missense mutations identified in the xenograft that were not present in the patient's primary tumor and there were no new nonsense mutations. None of the missense mutations acquired in the xenograft were in known cancer genes. We also demonstrate the feasibility of using the orthotopic neuroblastoma xenograft to test standard of care chemotherapy and molecular targeted therapeutics. Finally, we optimized a new approach to produce primary cultures of the neuroblastoma xenografts for high-throughput drug screening which can be used to test new combinations of therapeutic agents for neuroblastoma.

Targeting the DNA repair pathway in Ewing sarcoma.

Ewing sarcoma (EWS) is a tumor of the bone and soft tissue that primarily affects adolescents and young adults. With current therapies, 70% of patients with localized disease survive, but patients with metastatic or recurrent disease have a poor outcome. We found that EWS cell lines are defective in DNA break repair and are sensitive to PARP inhibitors (PARPis). PARPi-induced cytotoxicity in EWS cells was 10- to 1,000-fold higher after administration of the DNA-damaging agents irinotecan or temozolomide. We developed an orthotopic EWS mouse model and performed pharmacokinetic and pharmacodynamic studies using three different PARPis that are in clinical development for pediatric cancer. Irinotecan administered on a low-dose, protracted schedule previously optimized for pediatric patients was an effective DNA-damaging agent when combined with PARPis; it was also better tolerated than combinations with temozolomide. Combining PARPis with irinotecan and temozolomide gave complete and durable responses in more than 80% of the mice.

Ovarian tumor attachment, invasion, and vascularization reflect unique microenvironments in the peritoneum: insights from xenograft and mathematical models.

Ovarian cancer relapse is often characterized by metastatic spread throughout the peritoneal cavity with tumors attached to multiple organs. In this study, interaction of ovarian cancer cells with the peritoneal tumor microenvironment was evaluated in a xenograft model based on intraperitoneal injection of fluorescent SKOV3.ip1 ovarian cancer cells. Intra-vital microscopy of mixed GFP-red fluorescent protein (RFP) cell populations injected into the peritoneum demonstrated that cancer cells aggregate and attach as mixed spheroids, emphasizing the importance of homotypic adhesion in tumor formation. Electron microscopy provided high resolution structural information about local attachment sites. Experimental measurements from the mouse model were used to build a three-dimensional cellular Potts ovarian tumor model (OvTM) that examines ovarian cancer cell attachment, chemotaxis, growth, and vascularization. OvTM simulations provide insight into the relative influence of cancer cell-cell adhesion, oxygen availability, and local architecture on tumor growth and morphology. Notably, tumors on the mesentery, omentum, or spleen readily invade the "open" architecture, while tumors attached to the gut encounter barriers that restrict invasion and instead rapidly expand into the peritoneal space. Simulations suggest that rapid neovascularization of SKOV3.ip1 tumors is triggered by constitutive release of angiogenic factors in the absence of hypoxia. This research highlights the importance of cellular adhesion and tumor microenvironment in the seeding of secondary ovarian tumors on diverse organs within the peritoneal cavity. Results of the OvTM simulations indicate that invasion is strongly influenced by features underlying the mesothelial lining at different sites, but is also affected by local production of chemotactic factors. The integrated in vivo mouse model and computer simulations provide a unique platform for evaluating targeted therapies for ovarian cancer relapse.

Multi-scale modeling of tissues using CompuCell3D.

The study of how cells interact to produce tissue development, homeostasis, or diseases was, until recently, almost purely experimental. Now, multi-cell computer simulation methods, ranging from relatively simple cellular automata to complex immersed-boundary and finite-element mechanistic models, allow in silico study of multi-cell phenomena at the tissue scale based on biologically observed cell behaviors and interactions such as movement, adhesion, growth, death, mitosis, secretion of chemicals, chemotaxis, etc. This tutorial introduces the lattice-based Glazier-Graner-Hogeweg (GGH) Monte Carlo multi-cell modeling and the open-source GGH-based CompuCell3D simulation environment that allows rapid and intuitive modeling and simulation of cellular and multi-cellular behaviors in the context of tissue formation and subsequent dynamics. We also present a walkthrough of four biological models and their associated simulations that demonstrate the capabilities of the GGH and CompuCell3D.

Front instabilities and invasiveness of simulated 3D avascular tumors.

We use the Glazier-Graner-Hogeweg model to simulate three-dimensional (3D), single-phenotype, avascular tumors growing in an homogeneous tissue matrix (TM) supplying a single limiting nutrient. We study the effects of two parameters on tumor morphology: a diffusion-limitation parameter defined as the ratio of the tumor-substrate consumption rate to the substrate-transport rate, and the tumor-TM surface tension. This initial model omits necrosis and oxidative/hypoxic metabolism effects, which can further influence tumor morphology, but our simplified model still shows significant parameter dependencies. The diffusion-limitation parameter determines whether the growing solid tumor develops a smooth (noninvasive) or fingered (invasive) interface, as in our earlier two-dimensional (2D) simulations. The sensitivity of 3D tumor morphology to tumor-TM surface tension increases with the size of the diffusion-limitation parameter, as in 2D. The 3D results are unexpectedly close to those in 2D. Our results therefore may justify using simpler 2D simulations of tumor growth, instead of more realistic but more computationally expensive 3D simulations. While geometrical artifacts mean that 2D sections of connected 3D tumors may be disconnected, the morphologies of 3D simulated tumors nevertheless correlate with the morphologies of their 2D sections, especially for low-surface-tension tumors, allowing the use of 2D sections to partially reconstruct medically-important 3D-tumor structures.

3D multi-cell simulation of tumor growth and angiogenesis.

We present a 3D multi-cell simulation of a generic simplification of vascular tumor growth which can be easily extended and adapted to describe more specific vascular tumor types and host tissues. Initially, tumor cells proliferate as they take up the oxygen which the pre-existing vasculature supplies. The tumor grows exponentially. When the oxygen level drops below a threshold, the tumor cells become hypoxic and start secreting pro-angiogenic factors. At this stage, the tumor reaches a maximum diameter characteristic of an avascular tumor spheroid. The endothelial cells in the pre-existing vasculature respond to the pro-angiogenic factors both by chemotaxing towards higher concentrations of pro-angiogenic factors and by forming new blood vessels via angiogenesis. The tumor-induced vasculature increases the growth rate of the resulting vascularized solid tumor compared to an avascular tumor, allowing the tumor to grow beyond the spheroid in these linear-growth phases. First, in the linear-spherical phase of growth, the tumor remains spherical while its volume increases. Second, in the linear-cylindrical phase of growth the tumor elongates into a cylinder. Finally, in the linear-sheet phase of growth, tumor growth accelerates as the tumor changes from cylindrical to paddle-shaped. Substantial periods during which the tumor grows slowly or not at all separate the exponential from the linear-spherical and the linear-spherical from the linear-cylindrical growth phases. In contrast to other simulations in which avascular tumors remain spherical, our simulated avascular tumors form cylinders following the blood vessels, leading to a different distribution of hypoxic cells within the tumor. Our simulations cover time periods which are long enough to produce a range of biologically reasonable complex morphologies, allowing us to study how tumor-induced angiogenesis affects the growth rate, size and morphology of simulated tumors.

Contact-inhibited chemotaxis in de novo and sprouting blood-vessel growth.

Blood vessels form either when dispersed endothelial cells (the cells lining the inner walls of fully formed blood vessels) organize into a vessel network (vasculogenesis), or by sprouting or splitting of existing blood vessels (angiogenesis). Although they are closely related biologically, no current model explains both phenomena with a single biophysical mechanism. Most computational models describe sprouting at the level of the blood vessel, ignoring how cell behavior drives branch splitting during sprouting. We present a cell-based, Glazier-Graner-Hogeweg model (also called Cellular Potts Model) simulation of the initial patterning before the vascular cords form lumens, based on plausible behaviors of endothelial cells. The endothelial cells secrete a chemoattractant, which attracts other endothelial cells. As in the classic Keller-Segel model, chemotaxis by itself causes cells to aggregate into isolated clusters. However, including experimentally observed VE-cadherin-mediated contact inhibition of chemotaxis in the simulation causes randomly distributed cells to organize into networks and cell aggregates to sprout, reproducing aspects of both de novo and sprouting blood-vessel growth. We discuss two branching instabilities responsible for our results. Cells at the surfaces of cell clusters attempting to migrate to the centers of the clusters produce a buckling instability. In a model variant that eliminates the surface-normal force, a dissipative mechanism drives sprouting, with the secreted chemical acting both as a chemoattractant and as an inhibitor of pseudopod extension. Both mechanisms would also apply if force transmission through the extracellular matrix rather than chemical signaling mediated cell-cell interactions. The branching instabilities responsible for our results, which result from contact inhibition of chemotaxis, are both generic developmental mechanisms and interesting examples of unusual patterning instabilities.