Blood vessel wall shear stress determines regions of liposome accumulation in angiogenic vasculature.

Nanoparticles used for drug delivery often require intravenous administration exposing them to fluid forces within the vasculature, yet the impact of blood flow on nanoparticle delivery remains incompletely understood. Here, we utilized transgenic zebrafish embryos to investigate the relationship between the accumulation of fluorescently labeled PEGylated liposomes and various hemodynamic factors (such as flow velocity, wall shear stress (WSS), and flow pattern) across a wide range of angiogenic blood vessels. We reconstructed 3D models of vascular structures from confocal images and used computational fluid dynamics to calculate local WSS, velocities, and define flow patterns. The spatial distribution of fluorescently labeled liposomes was subsequently mapped within the same 3D space and correlated with local hemodynamic parameters. Through the integration of computational fluid dynamics and in vivo experimentation, we show that liposomes accumulated in vessel regions with WSS between 0.1-0.8 Pa, displaying an inverse linear correlation (R2 > 0.85) between time-averaged wall shear stress and liposome localization in vivo. Interestingly, flow pattern did not appear to impact liposome accumulation. Collectively, our findings suggest the potential of stealth liposomes for passive targeting of low-flow vasculature, including capillaries and intricate angiogenic vasculature resembling that of tumor vessel networks.

The development of brain pericytes requires expression of the transcription factor nkx3.1 in intermediate precursors.

Brain pericytes are one of the critical cell types that regulate endothelial barrier function and activity, thus ensuring adequate blood flow to the brain. The genetic pathways guiding undifferentiated cells into mature pericytes are not well understood. We show here that pericyte precursor populations from both neural crest and head mesoderm of zebrafish express the transcription factor nkx3.1 develop into brain pericytes. We identify the gene signature of these precursors and show that an nkx3.1-, foxf2a-, and cxcl12b-expressing pericyte precursor population is present around the basilar artery prior to artery formation and pericyte recruitment. The precursors later spread throughout the brain and differentiate to express canonical pericyte markers. Cxcl12b-Cxcr4 signaling is required for pericyte attachment and differentiation. Further, both nkx3.1 and cxcl12b are necessary and sufficient in regulating pericyte number as loss inhibits and gain increases pericyte number. Through genetic experiments, we have defined a precursor population for brain pericytes and identified genes critical for their differentiation.

SAM domain variants of EPHB4 associated with aberrant signaling are linked to lymphatic-related fetal hydrops and facial dysmorphology.

Variants in EPHB4 (Ephrin type B receptor 4), a transmembrane tyrosine kinase receptor, have been identified in individuals with various vascular anomalies including Capillary Malformation-Arteriovenous Malformation syndrome 2 and lymphatic-related (non-immune) fetal hydrops (LRHF). Here, we identify two novel variants in EPHB4 that disrupt the SAM domain in two unrelated individuals. Proband 1 presented within the LRHF phenotypic spectrum with hydrops, and proband 2 presented with large nuchal translucency prenatally that spontaneously resolved in addition to dysmorphic features on exam postnatally. These are the first disease associated variants identified that do not disrupt EPHB4 protein expression or tyrosine-kinase activity. We identify that EPHB4 SAM domain disruptions can lead to aberrant downstream signaling, with a loss of the SAM domain resulting in elevated MAPK signaling in proband 1, and a missense variant within the SAM domain resulting in increased cell proliferation in proband 2. This data highlights that a functional SAM domain is required for proper EPHB4 function and vascular development.

Vessel Metrics: A software tool for automated analysis of vascular structure in confocal imaging.

Images contain a wealth of information that is often under analyzed in biological studies. Developmental models of vascular disease are a powerful way to quantify developmentally regulated vessel phenotypes to identify the roots of the disease process. We present vessel Metrics, a software tool specifically designed to analyze developmental vascular microscopy images that will expedite the analysis of vascular images and provide consistency between research groups. We developed a segmentation algorithm that robustly quantifies different image types, developmental stages, organisms, and disease models at a similar accuracy level to a human observer. We validate the algorithm on confocal, lightsheet, and two photon microscopy data in a zebrafish model expressing fluorescent protein in the endothelial nuclei. The tool accurately segments data taken by multiple scientists on varying microscopes. We validate vascular parameters such as vessel density, network length, and diameter, across developmental stages, genetic mutations, and drug treatments, and show a favorable comparison to other freely available software tools. Additionally, we validate the tool in a mouse model. Vessel Metrics reduces the time to analyze experimental results, improves repeatability within and between institutions, and expands the percentage of a given vascular network analyzable in experiments.

rasa1-related arteriovenous malformation is driven by aberrant venous signalling.

Arteriovenous malformations (AVMs) develop where abnormal endothelial signalling allows direct connections between arteries and veins. Mutations in RASA1, a Ras GTPase activating protein, lead to AVMs in humans and, as we show, in zebrafish rasa1 mutants. rasa1 mutants develop cavernous AVMs that subsume part of the dorsal aorta and multiple veins in the caudal venous plexus (CVP) - a venous vascular bed. The AVMs progressively enlarge and fill with slow-flowing blood. We show that the AVM results in both higher minimum and maximum flow velocities, resulting in increased pulsatility in the aorta and decreased pulsatility in the vein. These hemodynamic changes correlate with reduced expression of the flow-responsive transcription factor klf2a. Remodelling of the CVP is impaired with an excess of intraluminal pillars, which is a sign of incomplete intussusceptive angiogenesis. Mechanistically, we show that the AVM arises from ectopic activation of MEK/ERK in the vein of rasa1 mutants, and that cell size is also increased in the vein. Blocking MEK/ERK signalling prevents AVM initiation in mutants. Alterations in venous MEK/ERK therefore drive the initiation of rasa1 AVMs.

Altered Expression of TMEM43 Causes Abnormal Cardiac Structure and Function in Zebrafish.

Arrhythmogenic cardiomyopathy (ACM) is an inherited heart muscle disease caused by heterozygous missense mutations within the gene encoding for the nuclear envelope protein transmembrane protein 43 (TMEM43). The disease is characterized by myocyte loss and fibro-fatty replacement, leading to life-threatening ventricular arrhythmias and sudden cardiac death. However, the role of TMEM43 in the pathogenesis of ACM remains poorly understood. In this study, we generated cardiomyocyte-restricted transgenic zebrafish lines that overexpress eGFP-linked full-length human wild-type (WT) TMEM43 and two genetic variants (c.1073C>T, p.S358L; c.332C>T, p.P111L) using the Tol2-system. Overexpression of WT and p.P111L-mutant TMEM43 was associated with transcriptional activation of the mTOR pathway and ribosome biogenesis, and resulted in enlarged hearts with cardiomyocyte hypertrophy. Intriguingly, mutant p.S358L TMEM43 was found to be unstable and partially redistributed into the cytoplasm in embryonic and adult hearts. Moreover, both TMEM43 variants displayed cardiac morphological defects at juvenile stages and ultrastructural changes within the myocardium, accompanied by dysregulated gene expression profiles in adulthood. Finally, CRISPR/Cas9 mutants demonstrated an age-dependent cardiac phenotype characterized by heart enlargement in adulthood. In conclusion, our findings suggest ultrastructural remodeling and transcriptomic alterations underlying the development of structural and functional cardiac defects in TMEM43-associated cardiomyopathy.

Stroke-associated intergenic variants modulate a human FOXF2 transcriptional enhancer.

SNPs associated with human stroke risk have been identified in the intergenic region between Forkhead family transcription factors FOXF2 and FOXQ1, but we lack a mechanism for the association. FoxF2 is expressed in vascular mural pericytes and is important for maintaining pericyte number and stabilizing small vessels in zebrafish. The stroke-associated SNPs are located in a previously unknown transcriptional enhancer for FOXF2, functional in human cells and zebrafish. We identify critical enhancer regions for FOXF2 gene expression, including binding sites occupied by transcription factors ETS1, RBPJ, and CTCF. rs74564934, a stroke-associated SNP adjacent to the ETS1 binding site, decreases enhancer function, as does mutation of RPBJ sites. rs74564934 is significantly associated with the increased risk of any stroke, ischemic stroke, small vessel stroke, and elevated white matter hyperintensity burden in humans. Foxf2 has a conserved function cross-species and is expressed in vascular mural pericytes of the vessel wall. Thus, stroke-associated SNPs modulate enhancer activity and expression of a regulator of vascular stabilization, FOXF2, thereby modulating stroke risk.

Hnrnpul1 controls transcription, splicing, and modulates skeletal and limb development in vivo.

Mutations in RNA-binding proteins can lead to pleiotropic phenotypes including craniofacial, skeletal, limb, and neurological symptoms. Heterogeneous nuclear ribonucleoproteins (hnRNPs) are involved in nucleic acid binding, transcription, and splicing through direct binding to DNA and RNA, or through interaction with other proteins in the spliceosome. We show a developmental role for Hnrnpul1 in zebrafish, resulting in reduced body and fin growth and missing bones. Defects in craniofacial tendon growth and adult-onset caudal scoliosis are also seen. We demonstrate a role for Hnrnpul1 in alternative splicing and transcriptional regulation using RNA-sequencing, particularly of genes involved in translation, ubiquitination, and DNA damage. Given its cross-species conservation and role in splicing, it would not be surprising if it had a role in human development. Whole-exome sequencing detected a homozygous frameshift variant in HNRNPUL1 in 2 siblings with congenital limb malformations, which is a candidate gene for their limb malformations. Zebrafish Hnrnpul1 mutants suggest an important developmental role of hnRNPUL1 and provide motivation for exploring the potential conservation of ancient regulatory circuits involving hnRNPUL1 in human development.

Endothelial Semaphorin 3fb regulates Vegf pathway-mediated angiogenic sprouting.

Vessel growth integrates diverse extrinsic signals with intrinsic signaling cascades to coordinate cell migration and sprouting morphogenesis. The pro-angiogenic effects of Vascular Endothelial Growth Factor (VEGF) are carefully controlled during sprouting to generate an efficiently patterned vascular network. We identify crosstalk between VEGF signaling and that of the secreted ligand Semaphorin 3fb (Sema3fb), one of two zebrafish paralogs of mammalian Sema3F. The sema3fb gene is expressed by endothelial cells in actively sprouting vessels. Loss of sema3fb results in abnormally wide and stunted intersegmental vessel artery sprouts. Although the sprouts initiate at the correct developmental time, they have a reduced migration speed. These sprouts have persistent filopodia and abnormally spaced nuclei suggesting dysregulated control of actin assembly. sema3fb mutants show simultaneously higher expression of pro-angiogenic (VEGF receptor 2 (vegfr2) and delta-like 4 (dll4)) and anti-angiogenic (soluble VEGF receptor 1 (svegfr1)/ soluble Fms Related Receptor Tyrosine Kinase 1 (sflt1)) pathway components. We show increased phospho-ERK staining in migrating angioblasts, consistent with enhanced Vegf activity. Reducing Vegfr2 kinase activity in sema3fb mutants rescues angiogenic sprouting. Our data suggest that Sema3fb plays a critical role in promoting endothelial sprouting through modulating the VEGF signaling pathway, acting as an autocrine cue that modulates intrinsic growth factor signaling.

Semaphorin 3fa Controls Ocular Vascularization From the Embryo Through to the Adult.

Purpose: Pathological blood vessel growth in the eye is implicated in several diseases that result in vision loss, including age-related macular degeneration and diabetic retinopathy. The limits of current disease therapies have created the need to identify and characterize new antiangiogenic drugs. Here, we identify the secreted chemorepellent semaphorin-3fa (Sema3fa) as an endogenous anti-angiogenic in the eye. Methods: We generated a CRISPR/Cas9 sema3fa zebrafish mutant line, sema3faca304/304. We assessed the retinal and choroidal vasculature in both larval and adult wild-type and sema3fa mutant zebrafish. Results: We find sema3fa mRNA is expressed by the ciliary marginal zone, neural retina, and retinal pigment epithelium of zebrafish larvae as choroidal vascularization emerges and the hyaloid/retinal vasculature is remodeled. The hyaloid vessels of sema3fa mutants develop appropriately but fail to remodel during the larval period, with adult mutants exhibiting a denser network of capillaries in the retinal periphery than seen in wild-type. The choroid vasculature is also defective in that it develops precociously, and aberrant, leaky sprouts are present in the normally avascular outer retina of both sema3faca304/304 larvae and adult fish. Conclusions: Sema3fa is a key endogenous signal for maintaining an avascular retina and preventing pathologic vascularization. Furthermore, we provide a new experimentally accessible model for studying choroid neovascularization (CNV) resulting from primary changes in the retinal environment that lead to downstream vessel infiltration.

Development of vascular regulation in the zebrafish embryo.

The thin endothelial wall of a newly formed vessel is under enormous stress at the onset of blood flow, rapidly acquiring support from mural cells (pericytes and vascular smooth muscle cells; vSMCs) during development. Mural cells then develop vasoactivity (contraction and relaxation) but we have little information as to when this first develops or the extent to which pericytes and vSMCs contribute. For the first time, we determine the dynamic developmental acquisition of vasoactivity in vivo in the cerebral vasculature of zebrafish. We show that pericyte-covered vessels constrict in response to α1-adrenergic receptor agonists and dilate in response to nitric oxide donors at 4 days postfertilization (dpf) but have heterogeneous responses later, at 6 dpf. In contrast, vSMC-covered vessels constrict at 6 dpf, and dilate at both stages. Using genetic ablation, we demonstrate that vascular constriction and dilation is an active response. Our data suggest that both pericyte- and vSMC-covered vessels regulate their diameter in early development, and that their relative contributions change over developmental time.

foxc1 is required for embryonic head vascular smooth muscle differentiation in zebrafish.

Vascular smooth muscle of the head derives from neural crest, but developmental mechanisms and early transcriptional drivers of the vSMC lineage are not well characterized. We find that in early development, the transcription factor foxc1b is expressed in mesenchymal cells that associate with the vascular endothelium. Using timelapse imaging, we observe that foxc1b expressing mesenchymal cells differentiate into acta2 expressing vascular mural cells. We show that in zebrafish, while foxc1b is co-expressed in acta2 positive smooth muscle cells that associate with large diameter vessels, it is not co-expressed in capillaries where pdgfrß positive pericytes are located. In addition to being an early marker of the lineage, foxc1 is essential for vSMC differentiation; we find that foxc1 loss of function mutants have defective vSMC differentiation and that early genetic ablation of foxc1b or acta2 expressing populations blocks vSMC differentiation. Furthermore, foxc1 is expressed upstream of acta2 and is required for acta2 expression in vSMCs. Using RNA-Seq we determine an enriched intersectional gene expression profile using dual expression of foxc1b and acta2 to identify novel vSMC markers. Taken together, our data suggests that foxc1 is a marker of vSMCs and plays a critical functional role in promoting their differentiation.

MicroRNA26 attenuates vascular smooth muscle maturation via endothelial BMP signalling.

As small regulatory transcripts, microRNAs (miRs) act as genetic 'fine tuners' of posttranscriptional events, and as genetic switches to promote phenotypic switching. The miR miR26a targets the BMP signalling effector, smad1. We show that loss of miR26a leads to hemorrhage (a loss of vascular stability) in vivo, suggesting altered vascular differentiation. Reduction in miR26a levels increases smad1 mRNA and phospho-Smad1 (pSmad1) levels. We show that increasing BMP signalling by overexpression of smad1 also leads to hemorrhage. Normalization of Smad1 levels through double knockdown of miR26a and smad1 rescues hemorrhage, suggesting a direct relationship between miR26a, smad1 and vascular stability. Using an in vivo BMP genetic reporter and pSmad1 staining, we show that the effect of miR26a on smooth muscle differentiation is non-autonomous; BMP signalling is active in embryonic endothelial cells, but not in smooth muscle cells. Nonetheless, increased BMP signalling due to loss of miR26a results in an increase in acta2-expressing smooth muscle cell numbers and promotes a differentiated smooth muscle morphology. Similarly, forced expression of smad1 in endothelial cells leads to an increase in smooth muscle cell number and coverage. Furthermore, smooth muscle phenotypes caused by inhibition of the BMP pathway are rescued by loss of miR26a. Taken together, our data suggest that miR26a modulates BMP signalling in endothelial cells and indirectly promotes a differentiated smooth muscle phenotype. Our data highlights how crosstalk from BMP-responsive endothelium to smooth muscle is important for smooth muscle differentiation.

Mutations in ILK, encoding integrin-linked kinase, are associated with arrhythmogenic cardiomyopathy.

Arrhythmogenic cardiomyopathy is a genetic heart muscle disorder characterized by fibro-fatty replacement of cardiomyocytes leading to life-threatening ventricular arrhythmias, heart failure, and sudden cardiac death. Mutations in genes encoding cardiac junctional proteins are known to cause about half of cases, while remaining genetic causes are unknown. Using exome sequencing, we identified 2 missense variants (p.H33N and p.H77Y) that were predicted to be damaging in the integrin-linked kinase (ILK) gene in 2 unrelated families. The p.H33N variant was found to be de novo. ILK links integrins and the actin cytoskeleton, and is essential for the maintenance of normal cardiac function. Both of the new variants are located in the ILK ankyrin repeat domain, which binds to the first LIM domain of the adaptor proteins PINCH1 and PINCH2. In silico binding studies proposed that the human variants disrupt the ILK-PINCH complex. Recombinant mutant ILK expressed in H9c2 rat myoblast cells shows aberrant prominent cytoplasmic localization compared to the wild-type. Expression of human wild-type and mutant ILK under the control of the cardiac-specific cmlc2 promotor in zebrafish shows that p.H77Y and p.P70L, a variant previously reported in a dilated cardiomyopathy family, cause cardiac dysfunction and death by about 2-3 weeks of age. Our findings provide genetic and functional evidence that ILK is a cardiomyopathy disease gene and highlight its relevance for diagnosis and genetic counseling of inherited cardiomyopathies.

Pericyte Biology in Zebrafish.

The zebrafish is an outstanding model for studying vascular biology in vivo. Pericytes and vascular smooth muscle cells can be imaged as they associate with vessels and provide stability and integrity to the vasculature. In zebrafish, pericytes associate with the cerebral and trunk vasculature on the second day of development, as assayed by pdgfrß and notch3 markers. In the head, cerebral pericytes are neural crest derived, except for the pericytes of the hindbrain vasculature, which are mesoderm derived. Similar to the hindbrain, pericytes on the trunk vasculature are also mesoderm derived. Regardless of their location, pericyte development depends on a complex interaction between blood flow and signalling pathways, such as Notch, SONIC HEDGEHOG and BMP signalling, all of which positively regulate pericyte numbers.Pericyte numbers rapidly increase as development proceeds in order to stabilize both the blood-brain barrier and the vasculature and hence, prevent haemorrhage. Consequently, compromised pericyte development results in compromised vascular integrity, which then evolves into detrimental pathologies. Some of these pathologies have been modelled in zebrafish by inducing mutations in the notch3, foxc1 and foxf2 genes. These zebrafish models provide insights into the mechanisms of disease as associated with pericyte biology. Going forward, these models may be key contributors in elucidating the role of vascular mural cells in regulating vessel diameter and hence, blood flow.

Nanoparticle localization in blood vessels: dependence on fluid shear stress, flow disturbances, and flow-induced changes in endothelial physiology.

Nanoparticles in the bloodstream are subjected to complex fluid forces as they move through the curves and branches of healthy or tumor vasculature. While nanoparticles are known to preferentially accumulate in angiogenic vessels, little is known about the flow conditions in these vessels and how these conditions may influence localization. Here, we report a methodology which combines confocal imaging of nanoparticle-injected transgenic zebrafish embryos, 3D modeling of the vasculature, particle mapping, and computational fluid dynamics, to quantitatively assess the effects of fluid forces on nanoparticle distribution in vivo. Six-fold lower accumulation was found in zebrafish arteries compared to the lower velocity veins. Nanoparticle localization varied inversely with shear stress. Highest accumulation was present in regions of disturbed flow found at branch points and curvatures in the vasculature. To further investigate cell-particle association under flow, human endothelial cells were exposed to nanoparticles under hemodynamic conditions typically found in human vessels. Physiological adaptations of endothelial cells to 20 hours of flow enhanced nanoparticle accumulation in regions of disturbed flow. Overall our results suggest that fluid shear stress magnitude, flow disturbances, and flow-induced changes in endothelial physiology modulate nanoparticle localization in angiogenic vessels.

Quantum dot interactions and flow effects in angiogenic zebrafish (Danio rerio) vessels and human endothelial cells.

Nanoparticle (NP) interactions with biological tissues are affected by the size, shape and surface chemistry of the NPs. Here we use in vivo (zebrafish) and in vitro (HUVEC) models to investigate association of quantum dots (QDs) with endothelial cells and the effect of fluid flow. After injection into the developing zebrafish, circulating QDs associate with endothelium and penetrate surrounding tissue parenchyma over time. Amino-functionalized QDs cluster, interact with cells, and clear more rapidly than carboxy-functionalized QDs in vivo, highlighting charge influences. QDs show stronger accumulation in slow-flowing, small caliber venous vessels than in fast-flowing high caliber arterial vessels. Parallel-plate flow experiments with HUVEC support these findings, showing reduced QD-EC association with increasing flow. In vivo, flow arrest after nanoparticle injection still results in venous accumulation at 18 h. Overall our results suggest that both QD charge and blood flow modulate particle-endothelial cell interactions.

The LIM-homeodomain transcription factor Islet2a promotes angioblast migration.

Angioblasts of the developing vascular system require many signaling inputs to initiate their migration, proliferation and differentiation into endothelial cells. What is less studied is which intrinsic cell factors interpret these extrinsic signals. Here, we show the Lim homeodomain transcription factor islet2a (isl2a) is expressed in the lateral posterior mesoderm prior to angioblast migration. isl2a deficient angioblasts show disorganized migration to the midline to form axial vessels and fail to spread around the tailbud of the embryo. Isl2a morphants have fewer vein cells and decreased vein marker expression. We demonstrate that isl2a is required cell autonomously in angioblasts to promote their incorporation into the vein, and is permissive for vein identity. Knockout of isl2a results in decreased migration and proliferation of angioblasts during intersegmental artery growth. Since Notch signaling controls both artery-vein identity and tip-stalk cell formation, we explored the interaction of isl2a and Notch. We find that isl2a expression is negatively regulated by Notch activity, and that isl2a positively regulates flt4, a VEGF-C receptor repressed by Notch during angiogenesis. Thus Isl2a may act as an intermediate between Notch signaling and genetic programs controlling angioblast number and migration, placing it as a novel transcriptional regulator of early angiogenesis.

Patterning mechanisms of the sub-intestinal venous plexus in zebrafish.

Despite considerable interest in angiogenesis, organ-specific angiogenesis remains less well characterized. The vessels that absorb nutrients from the yolk and later provide blood supply to the developing digestive system are primarily venous in origin. In zebrafish, these are the vessels of the Sub-intestinal venous plexus (SIVP) and they represent a new candidate model to gain an insight into the mechanisms of venous angiogenesis. Unlike other vessel beds in zebrafish, the SIVP is not stereotypically patterned and lacks obvious sources of patterning information. However, by examining the area of vessel coverage, number of compartments, proliferation and migration speed we have identified common developmental steps in SIVP formation. We applied our analysis of SIVP development to obd mutants that have a mutation in the guidance receptor PlexinD1. obd mutants show dysregulation of nearly all parameters of SIVP formation. We show that the SIVP responds to a unique combination of pathways that control both arterial and venous growth in other systems. Blocking Shh, Notch and Pdgf signaling has no effect on SIVP growth. However Vegf promotes sprouting of the predominantly venous plexus and Bmp promotes outgrowth of the structure. We propose that the SIVP is a unique model to understand novel mechanisms utilized in organ-specific angiogenesis.

Hematopoietic stem cell arrival triggers dynamic remodeling of the perivascular niche.

Hematopoietic stem and progenitor cells (HSPCs) can reconstitute and sustain the entire blood system. We generated a highly specific transgenic reporter of HSPCs in zebrafish. This allowed us to perform high-resolution live imaging on endogenous HSPCs not currently possible in mammalian bone marrow. Using this system, we have uncovered distinct interactions between single HSPCs and their niche. When an HSPC arrives in the perivascular niche, a group of endothelial cells remodel to form a surrounding pocket. This structure appears conserved in mouse fetal liver. Correlative light and electron microscopy revealed that endothelial cells surround a single HSPC attached to a single mesenchymal stromal cell. Live imaging showed that mesenchymal stromal cells anchor HSPCs and orient their divisions. A chemical genetic screen found that the compound lycorine promotes HSPC-niche interactions during development and ultimately expands the stem cell pool into adulthood. Our studies provide evidence for dynamic niche interactions upon stem cell colonization. PAPERFLICK: