Glioblastoma (GBM) remains the most malignant primary brain tumor, with a median survival rarely exceeding 2 years. Tumor heterogeneity and an immunosuppressive microenvironment are key factors contributing to the poor response rates of current therapeutic approaches. GBM-associated macrophages (GAMs) often exhibit immunosuppressive features that promote tumor progression. However, their dynamic interactions with GBM tumor cells remain poorly understood. Here, we used patient-derived GBM stem cell cultures and combined single-cell RNA sequencing of GAM-GBM co-cultures and real-time in vivo monitoring of GAM-GBM interactions in orthotopic zebrafish xenograft models to provide insight into the cellular, molecular, and spatial heterogeneity. Our analyses revealed substantial heterogeneity across GBM patients in GBM-induced GAM polarization and the ability to attract and activate GAMs – features that correlated with patient survival. Differential gene expression analysis, immunohistochemistry on original tumor samples, and knock-out experiments in zebrafish subsequently identified LGALS1 as a primary regulator of immunosuppression. Overall, our work highlights that GAM-GBM interactions can be studied in a clinically relevant way using co-cultures and avatar models, while offering new opportunities to identify promising immune-modulating targets.

A ‘GGGGCC’ repeat expansion in the first intron of the C9orf72 gene is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). The exact mechanism resulting in these neurodegenerative diseases remains elusive, but C9 repeat RNA toxicity has been implicated as a gain-of-function mechanism. Our aim was to use a zebrafish model for C9orf72 RNA toxicity to identify modifiers of the ALS-linked phenotype. We discovered that the RNA-binding protein heterogeneous nuclear ribonucleoprotein K (HNRNPK) reverses the toxicity of both sense and antisense repeat RNA, which is dependent on its subcellular localization and RNA recognition, and not on C9orf72 repeat RNA binding. We observed HNRNPK cytoplasmic mislocalization in C9orf72 ALS patient fibroblasts, induced pluripotent stem cell (iPSC)-derived motor neurons and post-mortem motor cortex and spinal cord, in line with a disrupted HNRNPK function in C9orf72 ALS. In C9orf72 ALS/FTD patient tissue, we discovered an increased nuclear translocation, but reduced expression of ribonucleotide reductase regulatory subunit M2 (RRM2), a downstream target of HNRNPK involved in the DNA damage response. Last but not least, we showed that increasing the expression of HNRNPK or RRM2 was sufficient to mitigate DNA damage in our C9orf72 RNA toxicity zebrafish model. Overall, our study strengthens the relevance of RNA toxicity as a pathogenic mechanism in C9orf72 ALS and demonstrates its link with an aberrant DNA damage response, opening novel therapeutic avenues for C9orf72 ALS/FTD.

Previous studies have shown that Vasohibin 1 (Vash1) is stimulated by VEGFs in endothelial cells and that its overexpression interferes with angiogenesis in vivo Recently, Vash1 was found to mediate tubulin detyrosination, a post-translational modification that is implicated in many cell functions, such as cell division. Here, we used the zebrafish embryo to investigate the cellular and subcellular mechanisms of Vash1 on endothelial microtubules during formation of the trunk vasculature. We show that microtubules within venous-derived secondary sprouts are strongly and selectively detyrosinated in comparison with other endothelial cells, and that this difference is lost upon vash1 knockdown. Vash1 depletion in zebrafish specifically affected secondary sprouting from the posterior cardinal vein, increasing endothelial cell divisions and cell number in the sprouts. We show that altering secondary sprout numbers and structure upon Vash1 depletion leads to defective lymphatic vessel formation and ectopic lymphatic progenitor specification in the zebrafish trunk.

Previous studies have shown that Vasohibin-1 (Vash-1) is stimulated by VEGFs in endothelial cells and that its overexpression interferes with angiogenesis in vivo . Recently, Vasohibin-1 was found to mediate tubulin detyrosination, a post-translational modification that is implicated in many cell functions, such as cell division. Here we used the zebrafish embryo to investigate the cellular and subcellular mechanisms of Vash-1 on endothelial microtubules during formation of the trunk vasculature. We show that microtubules within venous-derived secondary sprouts are strongly and selectively detyrosinated in comparison with other endothelial cells, and that this difference is lost upon vash-1 knockdown. Vasohibin-1 depletion in zebrafish specifically affected secondary sprouting from the posterior cardinal vein, increasing both the number of sprouts and endothelial cell divisions. We show that altering secondary sprout numbers and structure upon vash-1 depletion leads to a failure in the development and specification of lymphatic vessels of the zebrafish trunk. SUMMARY Vasohibin-1 mediated detyrosination of endothelial microtubules is selectively required for adequate behaviour of venous secondary sprouting and subsequent formation of functional lymphatics in the zebrafish trunk. ### Competing Interest Statement The authors have declared no competing interest.

How developing vascular networks acquire the right balance of arteries, veins and lymphatic vessels to efficiently supply and drain tissues is poorly understood. In zebrafish embryos, the robust and regular 50:50 global balance of intersegmental veins and arteries that form along the trunk prompts the intriguing question of how does the organism keep 'count'? Previous studies have suggested that the ultimate fate of an intersegmental vessel (ISV) is determined by the identity of the approaching secondary sprout emerging from the posterior cardinal vein. Here, we show that the formation of a balanced trunk vasculature involves an early heterogeneity in endothelial cell behaviour and Notch signalling activity in the seemingly identical primary ISVs that is independent of secondary sprouting and flow. We show that Notch signalling mediates the local patterning of ISVs, and an adaptive flow-mediated mechanism subsequently fine-tunes the global balance of arteries and veins along the trunk. We propose that this dual mechanism provides the adaptability required to establish a balanced network of arteries, veins and lymphatic vessels.

How developing vascular networks acquire the right balance of arteries, veins and lymphatics to efficiently supply and drain tissues is poorly understood. In zebrafish embryos, the robust and regular 50:50 global balance of intersegmental veins and arteries that form along the trunk, prompts the intriguing question how the organism keeps count. Previous studies suggest that the ultimate fate of an intersegmental vessel (ISV) is determined by the identity of the approaching secondary sprout emerging from the posterior cardinal vein (PCV). Here, using high time-resolution imaging, advanced cell tracking and computational analysis, we show that the formation of a balanced trunk vasculature involves an early heterogeneity in endothelial cell (EC) behavior in the seemingly identical primary ISVs and an adaptive flow-mediated mechanism that fine-tunes the balance of arteries and veins along the trunk. Detailed examination of the trunk vasculature dynamics throughout development reveals the frequent formation of three-way vascular connections between primary ISVs, the dorsal aorta (DA) and the PCV. Differential resolution of these connections into arteries or veins is mediated by polarized cell movement of the ECs within the ISV. Quantitative analysis of the cellular organization, polarity and directional movement of ECs in primary ISVs identifies an early differential behavior between future arteries and veins that is largely specified in the ECs of the individual ISVs, is dependent on Dll4/Notch, and occurs even in the absence of secondary sprouting. Notch signaling is involved in a local patterning mechanism normally favoring the formation of alternating arteries and veins. The global artery-vein balance is however maintained through a flow-dependent mechanism that can overwrite the local patterning. We propose that this dual mechanism driving arterio-venous identity during developmental angiogenesis in the zebrafish trunk provides the adaptability required to establish a balanced network of arteries, veins and lymphatic vessels.

Formation of a regularly branched blood vessel network is crucial in development and physiology. Here we show that the expression of the Notch ligand Dll4 fluctuates in individual endothelial cells within sprouting vessels in the mouse retina in vivo and in correlation with dynamic cell movement in mouse embryonic stem cell-derived sprouting assays. We also find that sprout elongation and branching associates with a highly differential phase pattern of Dll4 between endothelial cells. Stimulation with pathologically high levels of Vegf, or overexpression of Dll4, leads to Notch dependent synchronization of Dll4 fluctuations within clusters, both in vitro and in vivo. Our results demonstrate that the Vegf-Dll4/Notch feedback system normally operates to generate heterogeneity between endothelial cells driving branching, whilst synchronization drives vessel expansion. We propose that this sensitive phase transition in the behaviour of the Vegf-Dll4/Notch feedback loop underlies the morphogen function of Vegfa in vascular patterning.

29 Formation of a regularly branched blood vessel network is crucial in development and 30 physiology. Here we show that the expression of the Notch ligand Dll4 fluctuates in 31 individual endothelial cells within sprouting vessels in the mouse retina in vivo and in 32 correlation with dynamic cell movement in mouse embryonic stem cell-derived 33 sprouting assays. We also find that sprout elongation and branching associates with a 34 highly differential phase pattern of Dll4 between endothelial cells. Stimulation with 35 pathologically high levels of Vegf, or overexpression of Dll4, leads to Notch 36 dependent synchronization of Dll4 fluctuations within clusters, both in vitro and in 37

How vascular tubes build, maintain and adapt continuously perfused lumens to meet local metabolic needs remains poorly understood. Recent studies showed that blood flow itself plays a critical role in the remodelling of vascular networks, and suggested it is also required for the lumenization of new vascular connections. However, it is still unknown how haemodynamic forces contribute to the formation of new vascular lumens during blood vessel morphogenesis. Here we report that blood flow drives lumen expansion during sprouting angiogenesis in vivo by inducing spherical deformations of the apical membrane of endothelial cells, in a process that we have termed inverse blebbing. We show that endothelial cells react to these membrane intrusions by local and transient recruitment and contraction of actomyosin, and that this mechanism is required for single, unidirectional lumen expansion in angiogenic sprouts. Our work identifies inverse membrane blebbing as a cellular response to high external pressure. We show that in the case of blood vessels such membrane dynamics can drive local cell shape changes required for global tissue morphogenesis, shedding light on a pressure-driven mechanism of lumen formation in vertebrates.

Patterning of functional blood vessel networks is achieved by pruning of superfluous connections. The cellular and molecular principles of vessel regression are poorly understood. Here we show that regression is mediated by dynamic and polarized migration of endothelial cells, representing anastomosis in reverse. Establishing and analyzing the first axial polarity map of all endothelial cells in a remodeling vascular network, we propose that balanced movement of cells maintains the primitive plexus under low shear conditions in a metastable dynamic state. We predict that flow-induced polarized migration of endothelial cells breaks symmetry and leads to stabilization of high flow/shear segments and regression of adjacent low flow/shear segments.

Sprouting angiogenesis drives blood vessel growth in healthy and diseased tissues. Vegf and Dll4/Notch signalling cooperate in a negative feedback loop that specifies endothelial tip and stalk cells to ensure adequate vessel branching and function. Current concepts posit that endothelial cells default to the tip-cell phenotype when Notch is inactive. Here we identify instead that the stalk-cell phenotype needs to be actively repressed to allow tip-cell formation. We show this is a key endothelial function of neuropilin-1 (Nrp1), which suppresses the stalk-cell phenotype by limiting Smad2/3 activation through Alk1 and Alk5. Notch downregulates Nrp1, thus relieving the inhibition of Alk1 and Alk5, thereby driving stalk-cell behaviour. Conceptually, our work shows that the heterogeneity between neighbouring endothelial cells established by the lateral feedback loop of Dll4/Notch utilizes Nrp1 levels as the pivot, which in turn establishes differential responsiveness to TGF-β/BMP signalling.

[This corrects the article DOI: 10.1371/journal.pbio.1002125.].

Patterning of functional blood vessel networks is achieved by pruning of superfluous connections. The cellular and molecular principles of vessel regression are poorly understood. Here we show that regression is mediated by dynamic and polarized migration of endothelial cells, representing anastomosis in reverse. Establishing and analyzing the first axial polarity map of all endothelial cells in a remodeling vascular network, we propose that balanced movement of cells maintains the primitive plexus under low shear conditions in a metastable dynamic state. We predict that flow-induced polarized migration of endothelial cells breaks symmetry and leads to stabilization of high flow/shear segments and regression of adjacent low flow/shear segments.

Vessel sprouting by migrating tip and proliferating stalk endothelial cells (ECs) is controlled by genetic signals (such as Notch), but it is unknown whether metabolism also regulates this process. Here, we show that ECs relied on glycolysis rather than on oxidative phosphorylation for ATP production and that loss of the glycolytic activator PFKFB3 in ECs impaired vessel formation. Mechanistically, PFKFB3 not only regulated EC proliferation but also controlled the formation of filopodia/lamellipodia and directional migration, in part by compartmentalizing with F-actin in motile protrusions. Mosaic in vitro and in vivo sprouting assays further revealed that PFKFB3 overexpression overruled the pro-stalk activity of Notch, whereas PFKFB3 deficiency impaired tip cell formation upon Notch blockade, implying that glycolysis regulates vessel branching.

Summary The importance of the blood- and lymph vessels in the transport of essential fluids, gases, macromolecules and cells in vertebrates warrants optimal insight into the regulatory mechanisms underlying their development. Mouse and zebrafish models of lymphatic development are instrumental for gene discovery and gene characterization but are challenging for certain aspects, e.g. no direct accessibility of embryonic stages, or non-straightforward visualization of early lymphatic sprouting, respectively. We previously demonstrated that the Xenopus tadpole is a valuable model to study the processes of lymphatic development. However, a fluorescent Xenopus reporter directly visualizing the lymph vessels was lacking. Here, we created transgenic Tg(Flk1:eGFP) Xenopus laevis reporter lines expressing green fluorescent protein (GFP) in blood- and lymph vessels driven by the Flk1 (VEGFR-2) promoter. We also established a high-resolution fluorescent dye labeling technique selectively and persistently visualizing lymphatic endothelial cells, even in conditions of impaired lymph vessel formation or drainage function upon silencing of lymphangiogenic factors. Next, we applied the model to dynamically document blood and lymphatic sprouting and patterning of the initially avascular tadpole fin. Furthermore, quantifiable models of spontaneous or induced lymphatic sprouting into the tadpole fin were developed for dynamic analysis of loss-of-function and gain-of-function phenotypes using pharmacologic or genetic manipulation. Together with angiography and lymphangiography to assess functionality, Tg(Flk1:eGFP) reporter tadpoles readily allowed detailed lymphatic phenotyping of live tadpoles by fluorescence microscopy. The Tg(Flk1:eGFP) tadpoles represent a versatile model for functional lymph/angiogenomics and drug screening.

To supply tissues with nutrients and oxygen, the cardiovascular system forms a seamless, hierarchically branched, network of lumenized tubes. Here, we show that maintenance of patent vessel lumens requires the Bα regulatory subunit of protein phosphatase 2A (PP2A). Deficiency of Bα in zebrafish precludes vascular lumen stabilization resulting in perfusion defects. Similarly, inactivation of PP2A-Bα in cultured ECs induces tubulogenesis failure due to alteration of cytoskeleton dynamics, actomyosin contractility and maturation of cell-extracellular matrix (ECM) contacts. Mechanistically, we show that PP2A-Bα controls the activity of HDAC7, an essential transcriptional regulator of vascular stability. In the absence of PP2A-Bα, transcriptional repression by HDAC7 is abrogated leading to enhanced expression of the cytoskeleton adaptor protein ArgBP2. ArgBP2 hyperactivates RhoA causing inadequate rearrangements of the EC actomyosin cytoskeleton. This study unravels the first specific role for a PP2A holoenzyme in development: the PP2A-Bα/HDAC7/ArgBP2 axis maintains vascular lumens by balancing endothelial cytoskeletal dynamics and cell-matrix adhesion.