Intravoxel incoherent motion as a surrogate marker of perfused vascular density in rat brain.

Intravoxel incoherent motion (IVIM) MRI has emerged as a valuable technique for the assessment of tissue characteristics and perfusion. However, there is limited knowledge about the relationship between IVIM-derived measures and changes at the level of the vascular network. In this study, we investigated the potential use of IVIM MRI as a noninvasive tool for measuring changes in cerebral vascular density. Variations in quantitative immunohistochemical measurements of the vascular density across different regions in the rat brain (cortex, corpus callosum, hippocampus, thalamus, and hypothalamus) were related to the pseudo-diffusion coefficient D* and the flowing blood fraction f in healthy Wistar rats. We assessed whether region-wise differences in the vascular density are reflected by variations in the IVIM measurements and found a significant positive relationship with the pseudo-diffusion coefficient (p < 0.05, ß = 0.24). The effect of cerebrovascular alterations, such as blood-brain barrier (BBB) disruption on the perfusion-related IVIM parameters, is not well understood. Therefore, we investigated the effect of BBB disruption on the IVIM measures in a rat model of metabolic and vascular comorbidities (ZSF1 obese rat) and assessed whether this affects the relationship between the cerebral vascular density and the noninvasive IVIM measurements. We observed increased vascular permeability without detecting any differences in diffusivity, suggesting that BBB leakage is present before changes in the tissue integrity. We observed no significant difference in the relationship between cerebral vascular density and the IVIM measurements in our model of comorbidities compared with healthy normotensive rats.

Fluid flow as a driver of embryonic morphogenesis.

Fluid flow is a powerful morphogenic force during embryonic development. The physical forces created by flowing fluids can either create morphogen gradients or be translated by mechanosensitive cells into biological changes in gene expression. In this Primer, we describe how fluid flow is created in different systems and highlight the important mechanosensitive signalling pathways involved for sensing and transducing flow during embryogenesis. Specifically, we describe how fluid flow helps establish left-right asymmetry in the early embryo and discuss the role of flow of blood, lymph and cerebrospinal fluid in sculpting the embryonic cardiovascular and nervous system.

Prdm16 Supports Arterial Flow Recovery by Maintaining Endothelial Function.

[Figure: see text].

Personalised virtual gene panels reduce interpretation workload and maintain diagnostic rates of proband-only clinical exome sequencing for rare disorders.

PURPOSE: The increased adoption of genomic strategies in the clinic makes it imperative for diagnostic laboratories to improve the efficiency of variant interpretation. Clinical exome sequencing (CES) is becoming a valuable diagnostic tool, capable of meeting the diagnostic demand imposed by the vast array of different rare monogenic disorders. We have assessed a clinician-led and phenotype-based approach for virtual gene panel generation for analysis of targeted CES in patients with rare disease in a single institution. METHODS: Retrospective survey of 400 consecutive cases presumed by clinicians to have rare monogenic disorders, referred on singleton basis for targeted CES. We evaluated diagnostic yield and variant workload to characterise the usefulness of a clinician-led approach for generation of virtual gene panels that can incorporate up to three different phenotype-driven gene selection methods. RESULTS: Abnormalities of the nervous system (54.5%), including intellectual disability, head and neck (19%), skeletal system (16%), ear (15%) and eye (15%) were the most common clinical features reported in referrals. Combined phenotype-driven strategies for virtual gene panel generation were used in 57% of cases. On average, 7.3 variants (median=5) per case were retained for clinical interpretation. The overall diagnostic rate of proband-only CES using personalised phenotype-driven virtual gene panels was 24%. CONCLUSIONS: Our results show that personalised virtual gene panels are a cost-effective approach for variant analysis of CES, maintaining diagnostic yield and optimising the use of resources for clinical genomic sequencing in the clinic.

Mechanism of COVID-19-Induced Cardiac Damage from Patient, In Vitro and Animal Studies.

PURPOSE OF REVIEW: Though patient studies have been important for understanding the disease, research done in animals and cell culture complement our knowledge from patient data and provide insight into the mechanism of the disease. Understanding how COVID causes damage to the heart is essential to understanding possible long-term consequences. RECENT FINDINGS: COVID-19 is primarily a disease that attacks the lungs; however, it is known to have important consequences in many other tissues including the heart. Though myocarditis does occur in some patients, for most cases of cardiac damage, the injury arises from scarring either due to myocardial infarction or micro-infarction. The main focus is on how COVID affects blood flow through the coronaries. We review how endothelial activation leads to a hypercoagulative state in COVID-19. We also emphasize the effects that the cytokine storm can directly have on the regulation of coronary blood flow. Since the main two cell types that can be infected in the heart are pericytes and cardiomyocytes, we further describe the known effects on pericyte function and how that can further lead to microinfarcts within the heart. Though many of these effects are systemic, this review focuses on the consequences on cardiac tissue of this dysregulation and the role that it has in the formation of myocardial scarring.

B55α/PP2A Limits Endothelial Cell Apoptosis During Vascular Remodeling: A Complementary Approach To Disrupt Pathological Vessels?

RATIONALE: How endothelial cells (ECs) migrate and form an immature vascular plexus has been extensively studied. Yet, mechanisms underlying vascular remodeling remain poorly established. A better understanding of these processes may lead to the design of novel therapeutic strategies complementary to current angiogenesis inhibitors. OBJECTIVE: Starting from our previous observations that PP2A (protein phosphatase 2) regulates the HIF (hypoxia-inducible factor)/PHD-2 (prolyl hydroxylase 2)-constituted oxygen machinery, we hypothesized that this axis could play an important role during blood vessel formation, tissue perfusion, and oxygen restoration. METHODS AND RESULTS: We show that the PP2A regulatory subunit B55α is at the crossroad between vessel pruning and vessel maturation. Blood vessels with high B55α counter cell stress conditions and thrive for stabilization and maturation. When B55α is inhibited, ECs cannot cope with cell stress and undergo apoptosis, leading to massive pruning of nascent blood vessels. Mechanistically, we found that the B55α/PP2A complex restrains PHD-2 activity, promoting EC survival in a HIF-dependent manner, and furthermore dephosphorylates p38, altogether protecting ECs against cell stress occurring, for example, during the onset of blood flow. In tumors, EC-specific B55α deficiency induces pruning of immature-like tumor blood vessels resulting in delayed tumor growth and metastasis, without affecting nonpathological vessels. Consistently, systemic administration of a pan-PP2A inhibitor disrupts vascular network formation and tumor progression in vivo without additional effects on B55α-deficient vessels. CONCLUSIONS: Our data underline a unique role of the B55α/PP2A phosphatase complex in vessel remodeling and suggest the use of PP2A-inhibitors as potent antiangiogenic drugs targeting specifically nascent blood vessels with a mode-of-action complementary to VEGF-R (vascular endothelial growth factor receptor)-targeted therapies. Graphical Abstract: A graphical abstract is available for this article.

Organ-Specific Endothelial Cell Differentiation and Impact of Microenvironmental Cues on Endothelial Heterogeneity.

Endothelial cells throughout the body are heterogeneous, and this is tightly linked to the specific functions of organs and tissues. Heterogeneity is already determined from development onwards and ranges from arterial/venous specification to microvascular fate determination in organ-specific differentiation. Acknowledging the different phenotypes of endothelial cells and the implications of this diversity is key for the development of more specialized tissue engineering and vascular repair approaches. However, although novel technologies in transcriptomics and proteomics are facilitating the unraveling of vascular bed-specific endothelial cell signatures, still much research is based on the use of insufficiently specialized endothelial cells. Endothelial cells are not only heterogeneous, but their specialized phenotypes are also dynamic and adapt to changes in their microenvironment. During the last decades, strong collaborations between molecular biology, mechanobiology, and computational disciplines have led to a better understanding of how endothelial cells are modulated by their mechanical and biochemical contexts. Yet, because of the use of insufficiently specialized endothelial cells, there is still a huge lack of knowledge in how tissue-specific biomechanical factors determine organ-specific phenotypes. With this review, we want to put the focus on how organ-specific endothelial cell signatures are determined from development onwards and conditioned by their microenvironments during adulthood. We discuss the latest research performed on endothelial cells, pointing out the important implications of mimicking tissue-specific biomechanical cues in culture.

Linagliptin prevents left ventricular stiffening by reducing titin cleavage and hypophosphorylation.

The metabolic syndrome (MetS) is an escalating problem worldwide, causing left ventricular stiffening, an early characteristic of diastolic dysfunction for which no treatment exists. As diastolic dysfunction and stiffening in MetS patients are associated with increased circulating dipeptidyl peptidase-4 (DPP-4) levels, we investigated whether the clinically approved DPP-4 inhibitor linagliptin reduces left ventricular stiffness in MetS-induced cardiac disease. Sixteen-week-old obese ZSF1 rats, displaying the MetS and left ventricular stiffness, received linagliptin-supplemented or placebo diet for four weeks. Linagliptin significantly reduced obesity, hyperlipidaemia, and hyperglycaemia and improved left ventricular relaxation. This improved relaxation was related to decreased cardiac fibrosis and cardiomyocyte passive stiffness (Fpassive ). The reduced Fpassive was the result of titin isoform switching from the stiff N2B to the more flexible N2BA and increased phosphorylation of total titin and specifically its N2Bus region (S4080 and S3391). Importantly, DPP-4 directly cleaved titin in vitro, resulting in an increased Fpassive , which was prevented by simultaneous administration of linagliptin. In conclusion, linagliptin improves left ventricular stiffness in obese ZSF1 rats by preventing direct DPP4-mediated titin cleavage, as well as by modulating both titin isoform levels and phosphorylation. Reducing left ventricular stiffness by administering linagliptin might prevent MetS-induced early diastolic dysfunction in human.

Leukoencephalopathy with calcifications and cysts: Genetic and phenotypic spectrum.

Biallelic mutations in SNORD118, encoding the small nucleolar RNA U8, cause leukoencephalopathy with calcifications and cysts (LCC). Given the difficulty in interpreting the functional consequences of variants in nonprotein encoding genes, and the high allelic polymorphism across SNORD118 in controls, we set out to provide a description of the molecular pathology and clinical spectrum observed in a cohort of patients with LCC. We identified 64 affected individuals from 56 families. Age at presentation varied from 3 weeks to 67 years, with disease onset after age 40 years in eight patients. Ten patients had died. We recorded 44 distinct, likely pathogenic, variants in SNORD118. Fifty two of 56 probands were compound heterozygotes, with parental consanguinity reported in only three families. Forty nine of 56 probands were either heterozygous (46) or homozygous (three) for a mutation involving one of seven nucleotides that facilitate a novel intramolecular interaction between the 5' end and 3' extension of precursor-U8. There was no obvious genotype-phenotype correlation to explain the marked variability in age at onset. Complementing recently published functional analyses in a zebrafish model, these data suggest that LCC most often occurs due to combinatorial severe and milder mutations, with the latter mostly affecting 3' end processing of precursor-U8.

Simulating flow induced migration in vascular remodelling.

Shear stress induces directed endothelial cell (EC) migration in blood vessels leading to vessel diameter increase and induction of vascular maturation. Other factors, such as EC elongation and interaction between ECs and non-vascular areas are also important. Computational models have previously been used to study collective cell migration. These models can be used to predict EC migration and its effect on vascular remodelling during embryogenesis. We combined live time-lapse imaging of the remodelling vasculature of the quail embryo yolk sac with flow quantification using a combination of micro-Particle Image Velocimetry and computational fluid dynamics. We then used the flow and remodelling data to inform a model of EC migration during remodelling. To obtain the relation between shear stress and velocity in vitro for EC cells, we developed a flow chamber to assess how confluent sheets of ECs migrate in response to shear stress. Using these data as an input, we developed a multiphase, self-propelled particles (SPP) model where individual agents are driven to migrate based on the level of shear stress while maintaining appropriate spatial relationship to nearby agents. These agents elongate, interact with each other, and with avascular agents at each time-step of the model. We compared predicted vascular shape to real vascular shape after 4 hours from our time-lapse movies and performed sensitivity analysis on the various model parameters. Our model shows that shear stress has the largest effect on the remodelling process. Importantly, however, elongation played an especially important part in remodelling. This model provides a powerful tool to study the input of different biological processes on remodelling.

ALG13 X-linked intellectual disability: New variants, glycosylation analysis, and expanded phenotypes.

Pathogenic variants in ALG13 (ALG13 UDP-N-acetylglucosaminyltransferase subunit) cause an X-linked congenital disorder of glycosylation (ALG13-CDG) where individuals have variable clinical phenotypes that include developmental delay, intellectual disability, infantile spasms, and epileptic encephalopathy. Girls with a recurrent de novo c.3013C>T; p.(Asn107Ser) variant have normal transferrin glycosylation. Using a highly sensitive, semi-quantitative flow injection-electrospray ionization-quadrupole time-of-flight mass spectrometry (ESI-QTOF/MS) N-glycan assay, we report subtle abnormalities in N-glycans that normally account for <0.3% of the total plasma glycans that may increase up to 0.5% in females with the p.(Asn107Ser) variant. Among our 11 unrelated ALG13-CDG individuals, one male had abnormal serum transferrin glycosylation. We describe seven previously unreported subjects including three novel variants in ALG13 and report a milder neurodevelopmental course. We also summarize the molecular, biochemical, and clinical data for the 53 previously reported ALG13-CDG individuals. We provide evidence that ALG13 pathogenic variants may mildly alter N-linked protein glycosylation in both female and male subjects, but the underlying mechanism remains unclear.

Impaired SMAD1/5 Mechanotransduction and Cx37 (Connexin37) Expression Enable Pathological Vessel Enlargement and Shunting.

OBJECTIVE: Impaired ALK1 (activin receptor-like kinase-1)/Endoglin/BMP9 (bone morphogenetic protein 9) signaling predisposes to arteriovenous malformations (AVMs). Activation of SMAD1/5 signaling can be enhanced by shear stress. In the genetic disease hereditary hemorrhagic telangiectasia, which is characterized by arteriovenous malformations, the affected receptors are those involved in the activation of mechanosensitive SMAD1/5 signaling. To elucidate how genetic and mechanical signals interact in AVM development, we sought to identify targets differentially regulated by BMP9 and shear stress. Approach and Results: We identify Cx37 (Connexin37) as a differentially regulated target of ligand-induced and mechanotransduced SMAD1/5 signaling. We show that stimulation of endothelial cells with BMP9 upregulated Cx37, whereas shear stress inhibited this expression. This signaling was SMAD1/5-dependent, and in the absence of SMAD1/5, there was an inversion of the expression pattern. Ablated SMAD1/5 signaling alone caused AVM-like vascular malformations directly connecting the dorsal aorta to the inlet of the heart. In yolk sacs of mouse embryos with an endothelial-specific compound heterozygosity for SMAD1/5, addition of TNFα (tumor necrosis factor-α), which downregulates Cx37, induced development of these direct connections bypassing the yolk sac capillary bed. In wild-type embryos undergoing vascular remodeling, Cx37 was globally expressed by endothelial cells but was absent in regions of enlarging vessels. TNFα and endothelial-specific compound heterozygosity for SMAD1/5 caused ectopic regions lacking Cx37 expression, which correlated to areas of vascular malformations. Mechanistically, loss of Cx37 impairs correct directional migration under flow conditions. CONCLUSIONS: Our data demonstrate that Cx37 expression is differentially regulated by shear stress and SMAD1/5 signaling, and that reduced Cx37 expression is permissive for capillary enlargement into shunts.

Microvascular and lymphatic dysfunction in HFpEF and its associated comorbidities.

Heart failure with preserved ejection fraction (HFpEF) is a complex heterogeneous disease for which our pathophysiological understanding is still limited and specific prevention and treatment strategies are lacking. HFpEF is characterised by diastolic dysfunction and cardiac remodelling (fibrosis, inflammation, and hypertrophy). Recently, microvascular dysfunction and chronic low-grade inflammation have been proposed to participate in HFpEF development. Furthermore, several recent studies demonstrated the occurrence of generalized lymphatic dysfunction in experimental models of risk factors for HFpEF, including obesity, hypercholesterolaemia, type 2 diabetes mellitus (T2DM), hypertension, and aging. Here, we review the evidence for a combined role of coronary (micro)vascular dysfunction and lymphatic vessel alterations in mediating key pathological steps in HFpEF, including reduced cardiac perfusion, chronic low-grade inflammation, and myocardial oedema, and their impact on cardiac metabolic alterations (oxygen and nutrient supply/demand imbalance), fibrosis, and cardiomyocyte stiffness. We focus primarily on HFpEF caused by metabolic risk factors, such as obesity, T2DM, hypertension, and aging.

Shear Stress and VE-Cadherin.

Objective- Vascular fusion represents an important mechanism of vessel enlargement during development; however, its significance in postnatal vessel enlargement is still unknown. During fusion, 2 adjoining vessels merge to share 1 larger lumen. The aim of this research was to identify the molecular mechanism responsible for vascular fusion. Approach and Results- We previously showed that both low shear stress and DAPT ( N-[ N-(3,5-difluorophenacetyl)-L-alanyl]- S-phenylglycine t-butyl ester) treatment in the embryo result in a hyperfused vascular plexus and that increasing shear stress levels could prevent DAPT-induced fusion. We, therefore, investigated vascular endothelial-cadherin (VEC) phosphorylation because this is a common downstream target of low shear stress and DAPT treatment. VEC phosphorylation increases after DAPT treatment and decreased shear stress. The increased phosphorylation occurred independent of the cleavage of the Notch intracellular domain. Increasing shear stress rescues hyperfusion by DAPT treatment by causing the association of the phosphatase vascular endothelial-protein tyrosine phosphatase with VEC, counteracting VEC phosphorylation. Finally, Src (proto-oncogene tyrosine-protein kinase Src) inhibition prevents VEC phosphorylation in endothelial cells and can rescue hyperfusion induced by low shear stress and DAPT treatment. Moesin, a VEC target that was previously reported to mediate endothelial cell rearrangement during lumenization, relocalizes to cell membranes in vascular beds undergoing hyperfusion. Conclusions- This study provides the first evidence that VEC phosphorylation, induced by DAPT treatment and low shear stress, is involved in the process of fusion during vascular remodeling.

Characterisation of MHC class I genes in the koala.

Koala (Phascolarctos cinereus) populations are on the decline across the majority of Australia's mainland. Two major diseases threatening the long-term survival of affected koala populations are caused by obligate intracellular pathogens: Chlamydia and koala retrovirus (KoRV). To improve our understanding of the koala immune system, we characterised their major histocompatibility complex (MHC) class I genes, which are centrally involved in presenting foreign peptides derived from intracellular pathogens to cytotoxic T cells. A total of 11 class I genes were identified in the koala genome. Three genes, Phci-UA, UB and UC, showed relatively high genetic variability and were expressed in all 12 examined tissues, whereas the other eight genes had tissue-specific expression and limited polymorphism. Evidence of diversifying selection was detected in Phci-UA and UC, while gene conversion may have played a role in creating new alleles at Phci-UB. We propose that Phci-UA, UB and UC are likely classical MHC genes of koalas, and further research is needed to understand their role in koala chlamydial and KoRV infections.

Simultaneous imaging of blood flow dynamics and vascular remodelling during development.

Normal vascular development requires blood flow. Time-lapse imaging techniques have revolutionised our understanding of developmental biology, but measuring changes in blood flow dynamics has met with limited success. Ultrasound biomicroscopy and optical coherence tomography can concurrently image vascular structure and blood flow velocity, but these techniques lack the resolution to accurately calculate fluid forces such as shear stress. This is important because hemodynamic forces are biologically active and induce changes in the expression of genes important for vascular development. Regional variations in shear stress, rather than the overall level, control processes such as vessel enlargement and regression during vascular remodelling. We present a technique to concurrently visualise vascular remodelling and blood flow dynamics. We use an avian embryonic model and inject an endothelial-specific dye and fluorescent microspheres. The motion of the microspheres is captured with a high-speed camera and the velocity of the blood flow in and out of the region of interest is quantified by micro-particle image velocitymetry (µPIV). The vessel geometry and flow are used to numerically solve the flow physics with computational fluid dynamics (CFD). Using this technique, we can analyse changes in shear stress, pressure drops and blood flow velocities over a period of 10 to 16 h. We apply this to study the relationship between shear stress and chronic changes in vessel diameter during embryonic development, both in normal development and after TGFß stimulation. This technique allows us to study the interaction of biomolecular and biomechanical signals during vascular remodelling using an in vivo developmental model.

Flow dynamics control the location of sprouting and direct elongation during developmental angiogenesis.

Angiogenesis is tightly controlled by a number of signalling pathways. Although our understanding of the molecular mechanisms involved in angiogenesis has rapidly increased, the role that biomechanical signals play in this process is understudied. We recently developed a technique to simultaneously analyse flow dynamics and vascular remodelling by time-lapse microscopy in the capillary plexus of avian embryos and used this to study the hemodynamic environment present during angiogenic sprouting. We found that sprouts always form from a vessel at lower pressure towards a vessel at higher pressure, and that sprouts form at the location of a shear stress minimum, but avoid locations where two blood streams merge even if this point is at a lower level of shear stress than the sprouting location. Using these parameters, we were able to successfully predict sprout location in quail embryos. We also found that the pressure difference between two vessels is permissive to elongation, and that sprouts will either change direction or regress if the pressure difference becomes negative. Furthermore, the sprout elongation rate is proportional to the pressure difference between the two vessels. Our results show that flow dynamics are predictive of the location of sprout formation in perfused vascular networks and that pressure differences across the interstitium can guide sprout elongation.

TWIST1 Integrates Endothelial Responses to Flow in Vascular Dysfunction and Atherosclerosis.

RATIONALE: Blood flow-induced shear stress controls endothelial cell (EC) physiology during atherosclerosis via transcriptional mechanisms that are incompletely understood. The mechanosensitive transcription factor TWIST is expressed during embryogenesis, but its role in EC responses to shear stress and focal atherosclerosis is unknown. OBJECTIVE: To investigate whether TWIST regulates endothelial responses to shear stress during vascular dysfunction and atherosclerosis and compare TWIST function in vascular development and disease. METHODS AND RESULTS: The expression and function of TWIST1 was studied in EC in both developing vasculature and during the initiation of atherosclerosis. In zebrafish, twist was expressed in early embryonic vasculature where it promoted angiogenesis by inducing EC proliferation and migration. In adult porcine and murine arteries, TWIST1 was expressed preferentially at low shear stress regions as evidenced by quantitative polymerase chain reaction and en face staining. Moreover, studies of experimental murine carotid arteries and cultured EC revealed that TWIST1 was induced by low shear stress via a GATA4-dependent transcriptional mechanism. Gene silencing in cultured EC and EC-specific genetic deletion in mice demonstrated that TWIST1 promoted atherosclerosis by inducing inflammation and enhancing EC proliferation associated with vascular leakiness. CONCLUSIONS: TWIST expression promotes developmental angiogenesis by inducing EC proliferation and migration. In addition to its role in development, TWIST is expressed preferentially at low shear stress regions of adult arteries where it promotes atherosclerosis by inducing EC proliferation and inflammation. Thus, pleiotropic functions of TWIST control vascular disease and development.

Mosaic structural variation in children with developmental disorders.

Delineating the genetic causes of developmental disorders is an area of active investigation. Mosaic structural abnormalities, defined as copy number or loss of heterozygosity events that are large and present in only a subset of cells, have been detected in 0.2-1.0% of children ascertained for clinical genetic testing. However, the frequency among healthy children in the community is not well characterized, which, if known, could inform better interpretation of the pathogenic burden of this mutational category in children with developmental disorders. In a case-control analysis, we compared the rate of large-scale mosaicism between 1303 children with developmental disorders and 5094 children lacking developmental disorders, using an analytical pipeline we developed, and identified a substantial enrichment in cases (odds ratio = 39.4, P-value 1.073e - 6). A meta-analysis that included frequency estimates among an additional 7000 children with congenital diseases yielded an even stronger statistical enrichment (P-value 1.784e - 11). In addition, to maximize the detection of low-clonality events in probands, we applied a trio-based mosaic detection algorithm, which detected two additional events in probands, including an individual with genome-wide suspected chimerism. In total, we detected 12 structural mosaic abnormalities among 1303 children (0.9%). Given the burden of mosaicism detected in cases, we suspected that many of the events detected in probands were pathogenic. Scrutiny of the genotypic-phenotypic relationship of each detected variant assessed that the majority of events are very likely pathogenic. This work quantifies the burden of structural mosaicism as a cause of developmental disorders.

Blood flow can signal during angiogenesis not only through mechanotransduction, but also by affecting growth factor distribution.

Growth factors, such as VEGF, promote the sprouting of new blood vessels. Growth factors are generally produced far from the endothelium, and the transport of these proteins is often assumed to occur through diffusion. When sprouting occurs in a perfused vascular bed, however, interstitial flow is present that can modify protein transport. We recently developed a technique to analyze flow dynamics and vascular remodeling simultaneously in avian embryos. In this study, we extend our technique to model interstitial flow through the porous matrix of the mesenchymal tissue and use this to investigate how flow in the blood vessels affects the distribution of growth factors in the mesenchyme, using VEGF as a prototypical angiogenic molecule. We find that flow controls sprouting location and elongation, both through the direct action of mechanical force and through indirect effects on growth factor distribution. Most importantly, we find that the distribution of VEGF is regulated by interstitial flow, and the effect of diffusion of VEGF is negligible.