Angiogenesis requires co-ordination of multiple signalling inputs to regulate the behaviour of endothelial cells (ECs) as they form vascular networks. Vascular endothelial growth factor (VEGF) is essential for angiogenesis and induces downstream signalling pathways including increased cytosolic calcium levels. Here we show that transmembrane protein 33 (tmem33), which has no known function in multicellular organisms, is essential to mediate effects of VEGF in both zebrafish and human ECs. We find that tmem33 localises to the endoplasmic reticulum in zebrafish ECs and is required for cytosolic calcium oscillations in response to Vegfa. tmem33-mediated endothelial calcium oscillations are critical for formation of endothelial tip cell filopodia and EC migration. Global or endothelial-cell-specific knockdown of tmem33 impairs multiple downstream effects of VEGF including ERK phosphorylation, Notch signalling and embryonic vascular development. These studies reveal a hitherto unsuspected role for tmem33 and calcium oscillations in the regulation of vascular development.

Neurovascular coupling (through which local cerebral blood flow changes in response to neural activation are mediated) is impaired in many diseases including diabetes. Current preclinical rodent models of neurovascular coupling rely on invasive surgery and instrumentation, but transgenic zebrafish coupled with advances in imaging techniques allow non-invasive quantification of cerebrovascular anatomy, neural activation, and cerebral vessel haemodynamics. We therefore established a novel non-invasive, non-anaesthetised zebrafish larval model of neurovascular coupling, in which visual stimulus evokes neuronal activation in the optic tectum that is associated with a specific increase in red blood cell speed in tectal blood vessels. We applied this model to the examination of the effect of glucose exposure on cerebrovascular patterning and neurovascular coupling. We found that chronic exposure of zebrafish to glucose impaired tectal blood vessel patterning and neurovascular coupling. The nitric oxide donor sodium nitroprusside rescued all these adverse effects of glucose exposure on cerebrovascular patterning and function. Our results establish the first non-mammalian model of neurovascular coupling, offering the potential to perform more rapid genetic modifications and high-throughput screening than is currently possible using rodents. Furthermore, using this zebrafish model, we reveal a potential strategy to ameliorate the effects of hyperglycemia on cerebrovascular function.

The motivation of developing simple minimal models for neuro-glio-vascular (NGV) system arises from a recent modeling study elucidating the bidirectional information flow within the NGV system having 89 dynamic equations (1). While this was one of the first attempts at formulating a comprehensive model for neuro-glio-vascular system, it poses severe restrictions in scaling up to network levels. On the contrary, low-­dimensional models are convenient devices in simulating large networks that also provide an intuitive understanding of the complex interactions occurring within the NGV system. The key idea underlying the proposed models is to describe the glio-vascular system as a lumped system, which takes neural firing rate as input and returns an “energy” variable (analogous to ATP) as output. To this end, we present two models: biophysical neuro-energy (Model 1 with five variables), comprising KATP channel activity governed by neuronal ATP dynamics, and the dynamic threshold (Model 2 with three variables), depicting the dependence of neural firing threshold on the ATP dynamics. Both the models show different firing regimes, such as continuous spiking, phasic, and tonic bursting depending on the ATP production coefficient, ɛp, and external current. We then demonstrate that in a network comprising such energy-dependent neuron units, ɛp could modulate the local field potential (LFP) frequency and amplitude. Interestingly, low-frequency LFP dominates under low ɛp conditions, which is thought to be reminiscent of seizure-like activity observed in epilepsy. The proposed “neuron-energy” unit may be implemented in building models of NGV networks to simulate data obtained from multimodal neuroimaging systems, such as functional near infrared spectroscopy coupled to electroencephalogram and functional magnetic resonance imaging coupled to electroencephalogram. Such models could also provide a theoretical basis for devising optimal neurorehabilitation strategies, such as non-invasive brain stimulation for stroke patients.

We identify a novel endothelial membrane behaviour in transgenic zebrafish. Cerebral blood vessels extrude large transient spherical structures that persist for an average of 23 min before regressing into the parent vessel. We term these structures “kugeln”, after the German for sphere. Kugeln are only observed arising from the cerebral vessels and are present as late as 28 days post fertilization. Kugeln do not communicate with the vessel lumen and can form in the absence of blood flow. They contain little or no cytoplasm, but the majority are highly positive for nitric oxide reactivity. Kugeln do not interact with brain lymphatic endothelial cells (BLECs) and can form in their absence, nor do they perform a scavenging role or interact with macrophages. Inhibition of actin polymerization, Myosin II, or Notch signalling reduces kugel formation, while inhibition of VEGF or Wnt dysregulation (either inhibition or activation) increases kugel formation. Kugeln represent a novel Notch‐dependent NO‐containing endothelial organelle restricted to the cerebral vessels, of currently unknown function.

Zebrafish transgenic lines and light sheet fluorescence microscopy allow in-depth insights into three-dimensional vascular development in vivo. However, quantification of the zebrafish cerebral vasculature in 3D remains highly challenging. Here, we describe and test an image analysis workflow for 3D quantification of the total or regional zebrafish brain vasculature, called zebrafish vasculature quantification “ZVQ”. It provides the first landmark- or object-based vascular inter-sample registration of the zebrafish cerebral vasculature, producing Population Average Maps allowing rapid assessment of intra- and inter-group vascular anatomy. ZVQ also extracts a range of quantitative vascular parameters from a user-specified Region of Interest including volume, surface area, density, branching points, length, radius, and complexity. Application of ZVQ to thirteen experimental conditions, including embryonic development, pharmacological manipulations and morpholino induced gene knockdown, shows ZVQ is robust, allows extraction of biologically relevant information and quantification of vascular alteration, and can provide novel insights into vascular biology. To allow dissemination, the code for quantification, a graphical user interface, and workflow documentation are provided. Together, ZVQ provides the first open-source quantitative approach to assess the 3D cerebrovascular architecture in zebrafish.