Search

Supplemental Movie 2: Animated EVAC streamlines (grey lines), showing the path taken by interstitial fluid, overlaid on vascular perfusion measurements (colour scale, measured in vivo using arterial spin labelling) and the location of blood vessels (yellow structures, measured ex vivo using micro-CT). The data were acquired in an LS174T colorectal carcinoma tumour xenograft.

Additional data describing the reproducibility of convectionMRI data, and a figure showing the convectionMRI acquisition sequence.

Several distinct fluid flow phenomena occur in solid tumors, including intravascular blood flow and interstitial convection. Interstitial fluid pressure is often raised in solid tumors, which can limit drug delivery. To probe low-velocity flow in tumors resulting from raised interstitial fluid pressure, we developed a novel MRI technique named convection-MRI, which uses a phase-contrast acquisition with a dual-inversion vascular nulling preparation to separate intra- and extravascular flow. Here, we report the results of experiments in flow phantoms, numerical simulations, and tumor xenograft models to investigate the technical feasibility of convection-MRI. We observed a significant correlation between estimates of effective fluid pressure from convection-MRI with gold-standard, invasive measurements of interstitial fluid pressure in mouse models of human colorectal carcinoma. Our results show how convection-MRI can provide insights into the growth and responsiveness to vascular-targeting therapy in colorectal cancers.Significance: A noninvasive method for measuring low-velocity fluid flow caused by raised fluid pressure can be used to assess changes caused by therapy. Cancer Res; 78(7); 1859–72. ©2018 AACR.

Supplemental Movie 4: Animated EVAC streamlines (grey lines), showing the path taken by interstitial fluid, overlaid on vascular perfusion measurements (colour scale, measured in vivo using arterial spin labelling) and the location of blood vessels (yellow structures, measured ex vivo using micro-CT). The data were acquired in an SW1222 colorectal carcinoma tumour xenograft.

Supplemental Movie 1: Schematic movie illustrating the physical principles underpinning the convection-MRI sequence

Supplemental Movie 5: Interstitial convection, measured using EVAC-MRI and animated using a particle simulation (left), compared with MRI signal enhancement with time, following injection with a contrast agent (Gd-DTPA), in an LS174T tumour.

High-resolution copy of Figure 6b. Fluid flow measured in vivo with convectionMRI is shown as grey streamlines and perfusion is shown as a colorscale. The location of blood vessels is represented by yellow volume renderings, acquired using ex vivo microvascular casting and imaged with micro-CT.

Supplemental Movie 3: Animated EVAC streamlines (grey lines), showing the path taken by interstitial fluid, overlaid on vascular perfusion measurements (colour scale, measured in vivo using arterial spin labelling) and the location of blood vessels (yellow structures, measured ex vivo using micro-CT). The data were acquired in an LS174T colorectal carcinoma tumour xenograft.

PDF file - 1763KB, This supplemental material includes S1), a figure with the fit of the VERDICT model to the diffusion data from all the samples, S2), a figure demonstrating the quality of the diffusion images and the signal attenuation from an example sample of each cell line, TS1) a table with the orientation angle estimates of the VERDICT model, S3) a figure to demonstrate that VERDICT performs better than a tensor model, a two-tensor model and an isotropic version of the VERDICT model. We also include details of all the histological procedures as well as extra information about the choice of the vascular compartment in the VERDICT model.

We introduce a hybrid, discrete-continuum method to model blood pressure across large microvascular networks with incomplete architectural information. Our multiscale approach couples a 1D Poiseuille flow description in large, discrete arteriolar and venular networks, via point sources of flux, to a continuum-based Darcy model for transport in the capillary bed. We test our hybrid framework using a vascular network imaged from the mouse brain medulla / pons using multi-fluorescence high-resolution episcopic microscopy. The fully-resolved vascular network is used to predicted the hydraulic conductivity of the capillary network and to generate a fully-discrete pressure solution to benchmark our novel method against. We show that the discrete-continuum methodology is a computationally tractable and effective tool for predicting blood pressure in real-world microvascular tissues when capillary microvessels are not well-defined.

Studies of the structural and molecular features of the lymphatic vasculature, which clears fluid, macromolecules and leukocytes from the tissue microenvironment, have largely relied on animal models, with limited information in human organs beyond traditional immunohistochemical assessment. Here, we use three-dimensional imaging and single-cell RNA-sequencing to study lymphatics in the human kidney. We found a hierarchical arrangement of lymphatic vessels within human kidneys, initiating along specialised nephron epithelium in the renal cortex and displaying a distinct, kidney-specific transcriptional profile. In chronic transplant rejection we found kidney allograft lymphatic expansion alongside a loss of structural hierarchy, with human leukocyte antigen-expressing lymphatic vessels infiltrating the medulla, presenting a putative target for alloreactive antibodies. This occurred concurrently with lymphatic vessels invading and interconnecting tertiary lymphoid structures at early stages of lymphocyte colonisation. Analysis of intercellular signalling revealed upregulation of co-inhibitory molecule-mediated CD4+T cell-lymphatic crosstalk in rejecting kidneys, potentially acting to limit local alloimmune responses. Overall, we delineate novel structural and molecular features of human kidney lymphatics and reveal perturbations to their phenotype and transcriptome in the context of alloimmunity.SUMMARYLymphatics regulate fluid balance and immune cell accumulation but are under-studied in human organs such as the kidney. Jafree and colleagues profiled human kidney lymphatics using three-dimensional imaging and single-cell RNA-sequencing, revealing structural and transcriptional perturbations in rejecting kidney transplants.

AIM: The adapter protein p130Cas, encoded by the Bcar1 gene, is a key regulator of cell movement, adhesion, and cell cycle control in diverse cell types. Bcar1 constitutive knockout mice are embryonic lethal by embryonic days (E) 11.5-12.5, but the role of Bcar1 in embryonic development remains unclear. Here, we investigated the role of Bcar1 specifically in cardiovascular development and defined the cellular and molecular mechanisms disrupted following targeted Bcar1 deletions. METHODS AND RESULTS: We crossed Bcar1 floxed mice with Cre transgenic lines allowing for cell-specific knockout either in smooth muscle and early cardiac tissues (SM22-Cre), mature smooth muscle cells (smMHC-Cre), endothelial cells (Tie2-Cre), second heart field cells (Mef2c-Cre), or neural crest cells (NCC) (Pax3-Cre) and characterised these conditional knock outs using a combination of histological and molecular biology techniques.Conditional knockout of Bcar1 in SM22-expressing smooth muscle cells and cardiac tissues (Bcar1SM22KO) was embryonically lethal from E14.5-15.5 due to severe cardiovascular defects, including abnormal ventricular development and failure of outflow tract (OFT) septation leading to a single outflow vessel reminiscent of persistent truncus arteriosus. SM22-restricted loss of Bcar1 was associated with failure of OFT cushion cells to undergo differentiation to septal mesenchymal cells positive for SMC-specific α-actin, and disrupted expression of proteins and transcription factors involved in epithelial-to-mesenchymal transformation (EMT). Furthermore, knockout of Bcar1 specifically in NCC (Bcar1PAX3KO) recapitulated part of the OFT septation and aortic sac defects seen in the Bcar1SM22KO mutants, indicating a cell-specific requirement for Bcar1 in NCC essential for OFT septation. In contrast, conditional knockouts of Bcar1 in differentiated smooth muscle, endothelial cells, and second heart field cells survived to term and were phenotypically normal at birth and post-natally. CONCLUSIONS: Our work reveals a cell-specific requirement for Bcar1 in NCC, early myogenic and cardiac cells, essential for OFT septation, myocardialisation and EMT/cell cycle regulation and differentiation to myogenic lineages. TRANSLATIONAL PERSPECTIVE: The molecular pathways coordinating cardiogenesis and the remodelling of the OFT are complex, and dysregulation of these pathways causes human heart defects. Our findings highlight a specific requirement for Bcar1 essential for cardiogenesis. Furthermore, the failure of OFT septation in Bcar1SM22KO mice resembles persistent truncus arteriosus (PTA), a feature of several human congenital heart diseases, including DiGeorge Syndrome. Our findings have implications for the mechanisms underlying the pathogenesis of congenital heart disease, and suggest that mice with conditional Bcar1 deletions may be useful models for dissecting mechanisms involved in the pathogenesis of human heart defects.

Objectives Effects of liver disease on portal venous (PV), hepatic arterial (HA), total liver blood flow (TLBF), and cardiac function are poorly understood. Terlipressin modulates PV flow but effects on HA, TLBF, and sepsis/acute-on-chronic liver failure (ACLF)-induced haemodynamic changes are poorly characterised. In this study, we investigated the effects of terlipressin and sepsis/ACLF on hepatic haemodynamics and cardiac function in a rodent cirrhosis model using caval subtraction phase-contrast (PC) MRI and cardiac cine MRI. Methods Sprague-Dawley rats (n = 18 bile duct–ligated (BDL), n = 16 sham surgery controls) underwent caval subtraction PCMRI to estimate TLBF and HA flow and short-axis cardiac cine MRI for systolic function at baseline, following terlipressin and lipopolysaccharide (LPS) infusion, to model ACLF. Results All baseline hepatic haemodynamic/cardiac systolic function parameters (except heart rate and LV mass) were significantly different in BDL rats. Following terlipressin, baseline PV flow (sham 181.4 ± 12.1 ml/min/100 g; BDL 68.5 ± 10.1 ml/min/100 g) reduced (sham − 90.3 ± 11.1 ml/min/100 g, p < 0.0001; BDL − 31.0 ± 8.0 ml/min/100 g, p = 0.02), sham baseline HA flow (33.0 ± 11.3 ml/min/100 g) increased (+ 92.8 ± 21.3 ml/min/100 g, p = 0.0003), but BDL baseline HA flow (83.8 ml/min/100 g) decreased (− 34.4 ± 7.5 ml/min/100 g, p = 0.11). Sham baseline TLBF (214.3 ± 16.7 ml/min/100 g) was maintained (+ 2.5 ± 14.0 ml/min/100 g, p > 0.99) but BDL baseline TLBF (152.3 ± 18.7 ml/min/100 g) declined (− 65.5 ± 8.5 ml/min/100 g, p = 0.0004). Following LPS, there were significant differences between cohort and change in HA fraction (p = 0.03) and TLBF (p = 0.01) with BDL baseline HA fraction (46.2 ± 4.6%) reducing (− 20.9 ± 7.5%, p = 0.03) but sham baseline HA fraction (38.2 ± 2.0%) remaining unchanged (+ 2.9 ± 6.1%, p > 0.99). Animal cohort and change in systolic function interactions were significant only for heart rate (p = 0.01) and end-diastolic volume (p = 0.03). Conclusions Caval subtraction PCMRI and cardiac MRI in a rodent model of cirrhosis demonstrate significant baseline hepatic haemodynamic/cardiac differences, failure of the HA buffer response post-terlipressin and an altered HA fraction response in sepsis, informing potential translation to ACLF patients. Key Points Caval subtraction phase-contrast and cardiac MRI demonstrate: • Significant differences between cirrhotic/non-cirrhotic rodent hepatic blood flow and cardiac systolic function at baseline. • Failure of the hepatic arterial buffer response in cirrhotic rodents in response to terlipressin. • Reductions in hepatic arterial flow fraction in the setting of acute-on-chronic liver failure.

Imaging intact human organs from the organ to the cellular scale in three dimensions is a goal of biomedical imaging. To meet this challenge, we developed hierarchical phase-contrast tomography (HiP-CT), an X-ray phase propagation technique using the European Synchrotron Radiation Facility (ESRF)’s Extremely Brilliant Source (EBS). The spatial coherence of the ESRF-EBS combined with our beamline equipment, sample preparation and scanning developments enabled us to perform non-destructive, three-dimensional (3D) scans with hierarchically increasing resolution at any location in whole human organs. We applied HiP-CT to image five intact human organ types: brain, lung, heart, kidney and spleen. HiP-CT provided a structural overview of each whole organ followed by multiple higher-resolution volumes of interest, capturing organotypic functional units and certain individual specialized cells within intact human organs. We demonstrate the potential applications of HiP-CT through quantification and morphometry of glomeruli in an intact human kidney and identification of regional changes in the tissue architecture in a lung from a deceased donor with coronavirus disease 2019 (COVID-19).

3D microscopy of large biological samples (>0.5 cm3) is transforming biological research. Many existing techniques require trade‐offs between image resolution, sample size, and method complexity. A simple robust instrument with the potential to conduct large‐volume 3D imaging currently exists in the form of the optical high‐resolution episcopic microscopy (HREM). However, the development of the instrument to date is limited to single‐fluorescent wavelength imaging with nonspecific eosin staining. Herein, developments to realize the potential of the HREM to become multifluorescent high‐resolution episcopic microscopy (MF‐HREM) are presented. MF‐HREM is a serial‐sectioning and block‐facing wide‐field fluorescence imaging technique, which does not require tissue clearing or optical sectioning. Multiple developments are detailed in sample preparation and image postprocessing to enable multiple specific stains in large samples and show how these enable segmentation and quantification of the data. The application of MF‐HREM is demonstrated in a variety of biological contexts: 3D imaging of whole tumor vascular networks and tumor cell invasion in xenograft tumors up to 7.5 mm3 at resolutions of 2.75 μm, quantification of glomeruli volume in the adult mouse kidney, and quantification of vascular networks and white‐matter track orientation in adult mouse brain.

Introduction To combine numerical simulations, in vitro and in vivo experiments to evaluate the feasibility of measuring diffusion exchange across the cell membrane with diffusion exchange spectroscopy (DEXSY). Methods DEXSY acquisitions were simulated over a range of permeabilities in nerve tissue and yeast substrates. In vitro measurements were performed in a yeast substrate and in vivo measurements in mouse tumor xenograft models, all at 9.4 T. Results Diffusion exchange was observed in simulations over a physiologically relevant range of cell permeability values. In vitro and in vivo measures also provided evidence of diffusion exchange, which was quantified with the Diffusion Exchange Index (DEI). Conclusions Our findings provide preliminary evidence that DEXSY can be used to make in vivo measurements of diffusion exchange and cell membrane permeability.

Nanoparticles have the potential to increase the efficacy of anticancer drugs whilst reducing off-target side effects. However, there remain uncertainties regarding the cellular uptake kinetics of nanoparticles which could have implications for nanoparticle design and delivery. Polymersomes are nanoparticle candidates for cancer therapy which encapsulate chemotherapy drugs. Here we develop a mathematical model to simulate the uptake of polymersomes via endocytosis, a process by which polymersomes bind to the cell surface before becoming internalised by the cell where they then break down, releasing their contents which could include chemotherapy drugs. We focus on two in vitro configurations relevant to the testing and development of cancer therapies: a well-mixed culture model and a tumour spheroid setup. Our mathematical model of the well-mixed culture model comprises a set of coupled ordinary differential equations for the unbound and bound polymersomes and associated binding dynamics. Using a singular perturbation analysis we identify an optimal number of ligands on the polymersome surface which maximises internalised polymersomes and thus intracellular chemotherapy drug concentration. In our mathematical model of the spheroid, a multiphase system of partial differential equations is developed to describe the spatial and temporal distribution of bound and unbound polymersomes via advection and diffusion, alongside oxygen, tumour growth, cell proliferation and viability. Consistent with experimental observations, the model predicts the evolution of oxygen gradients leading to a necrotic core. We investigate the impact of two different internalisation functions on spheroid growth, a constant and a bond dependent function. It was found that the constant function yields faster uptake and therefore chemotherapy delivery. We also show how various parameters, such as spheroid permeability, lead to travelling wave or steady-state solutions.

Human organs are complex, three-dimensional and multiscale systems. Spatially mapping the human body down through its hierarchy, from entire organs to their individual functional units and specialised cells, is a major obstacle to fully understanding health and disease. To meet this challenge, we developed hierarchical phase-contrast tomography (HiP-CT), an X-ray phase propagation technique utilising the European Synchrotron Radiation Facility’s Extremely Brilliant Source: the world’s first high-energy 4th generation X-ray source. HiP-CT enabled three-dimensional and non-destructive imaging at near-micron resolution in soft tissues at one hundred thousand times the voxel size whilst maintaining the organ’s structure. We applied HiP-CT to image five intact human parenchymal organs: brain, lung, heart, kidney and spleen. These were hierarchically assessed with HiP-CT, providing a structural overview of the whole organ alongside detail of the organ’s individual functional units and cells. The potential applications of HiP-CT were demonstrated through quantification and morphometry of glomeruli in an intact human kidney, and identification of regional changes to the architecture of the air-tissue interface and alveolar morphology in the lung of a deceased COVID-19 patient. Overall, we show that HiP-CT is a powerful tool which can provide a comprehensive picture of structural information for whole intact human organs, encompassing precise details on functional units and their constituent cells to better understand human health and disease.

Noninvasive measurements of liver perfusion and fibrosis in cirrhotic small animals can help develop treatments for haemodynamic complications of liver disease. Here, we measure liver perfusion in cirrhotic rodents using flow‐sensitive alternating inversion recovery arterial spin labelling (FAIR ASL), evaluating agreement with previously validated caval subtraction phase‐contrast magnetic resonance imaging (PCMRI) total liver blood flow (TLBF). Baseline differences in cirrhotic rodents and the haemodynamic effects of acute inflammation were investigated using FAIR ASL and tissue T1. Sprague–Dawley rats (nine bile duct ligated [BDL] and ten sham surgery controls) underwent baseline hepatic FAIR ASL with T1 measurement and caval subtraction PCMRI (with two‐dimensional infra‐/supra‐hepatic inferior vena caval studies), induction of inflammation with intravenous lipopolysaccharide (LPS) and repeat liver FAIR ASL with T1 measurement after ~90 minutes. The mean difference between FAIR ASL hepatic perfusion and caval subtraction PCMRI TLBF was −51 ± 30 ml/min/100 g (Bland–Altman 95% limits‐of‐agreement ±258 ml/min/100 g). The FAIR ASL coefficient of variation was smaller than for caval subtraction PCMRI (29.3% vs 50.1%; P = .03). At baseline, FAIR ASL liver perfusion was lower in BDL rats (199 ± 32 ml/min/100 g vs sham 316 ± 24 ml/min/100 g; P = .01) but liver T1 was higher (BDL 1533 ± 50 vs sham 1256 ± 18 ms; P = .0004). Post‐LPS FAIR ASL liver perfusion response differences were observed between sham/BDL rats (P = .02), approaching significance in sham (+78 ± 33 ml/min/100 g; P = .06) but not BDL rats (−49 ± 40 ml/min/100 g; P = .47). Post‐LPS differences in liver tissue T1 were nonsignificant (P = .35). FAIR ASL hepatic perfusion and caval subtraction PCMRI TLBF agreement was modest, with significant baseline FAIR ASL liver perfusion and tissue T1 differences in rodents with advanced cirrhosis compared with controls. Following inflammatory stress, differences in hepatic perfusion response were detected between cirrhotic/control animals, but liver T1 was unaffected. Findings underline the potential of FAIR ASL in the assessment of vasoactive treatments for patients with chronic liver disease and inflammation., Flow‐sensitive alternating inversion recovery arterial spin labelling (FAIR ASL) was used to measure liver perfusion and T1 at baseline and following inflammatory stress in cirrhotic and control rats. Baseline agreement of FAIR ASL liver perfusion with previously validated caval subtraction phase‐contrast MRI was modest. Following inflammatory stress, significant differences in hepatic perfusion response were detected between cirrhotic and control animals but liver T1 was unaffected, highlighting the potential of FAIR ASL in the assessment of chronic liver disease and inflammation.

Three-dimensional microscopy of large biological samples (>0.5 cm3) is transforming biological research. Many existing techniques require trade-offs between image resolution, sample size and method complexity. A simple robust instrument with the potential to perform large volume 3D imaging currently exists in the form of the Optical HREM, however the development of the instrument to date is limited to single fluorescent wavelength imaging with non-specific eosin staining. This work presents developments to realize the potential of the HREM to become Multi-fluorescent High Resolution Episcopic Microscopy (MF-HREM).MF-HREM is a serial-sectioning and block-facing wide-field fluorescence imaging technique, which does not require tissue clearing or optical sectioning. Multiple developments are detailed in sample preparation and image post-processing to enable multiple specific stains in large samples, and show how these enable segmentation and quantification of the data. The application of MF-HREM is demonstrated in a variety of biological contexts: 3D imaging of whole tumor vascular networks and tumor cell invasion in xenograft tumors up to 7.5 mm3 at resolutions of 2.75 μm, quantification of glomeruli volume in the adult mouse kidney, and quantification of vascular networks and white matter track orientation in adult mouse brain.

Advances in biological imaging have accelerated our understanding of human physiology in both health and disease. As these advances have developed, the opportunities gained by integrating with cutting-edge mathematical models have become apparent yet remain challenging. Combined imaging-modelling approaches provide unprecedented opportunity to correlate data on tissue architecture and function, across length and time scales, to better understand the mechanisms that underpin fundamental biology and also to inform clinical decisions. Here we discuss the opportunities and challenges of such approaches, providing literature examples across a range of organ systems. Given the breadth of the field we focus on the intersection of continuum modelling and in vivo imaging applied to the vasculature and blood flow, though our rationale and conclusions extend widely. We propose three key research pillars (image acquisition, image processing, mathematical modelling) and present their respective advances as well as future opportunity via better integration. Multidisciplinary efforts that develop imaging and modelling tools concurrently, and share them open-source with the research community, provide exciting opportunity for advancing these fields.

Understanding the uptake of a drug by diseased tissue, and the drug's subsequent spatiotemporal distribution, are central factors in the development of effective targeted therapies. However, the interaction between the pathophysiology of diseased tissue and individual therapeutic agents can be complex, and can vary across tissue types and across subjects. Here, we show that the combination of mathematical modelling, high-resolution optical imaging of intact and optically cleared tumour tissue from animal models, and in vivo imaging of vascular perfusion predicts the heterogeneous uptake, by large tissue samples, of specific therapeutic agents, as well as their spatiotemporal distribution. In particular, by using murine models of colorectal cancer and glioma, we report and validate predictions of steady-state blood flow and intravascular and interstitial fluid pressure in tumours, of the spatially heterogeneous uptake of chelated gadolinium by tumours, and of the effect of a vascular disrupting agent on tumour vasculature.

Several distinct fluid flow phenomena occur in solid tumors, including intravascular blood flow and interstitial convection. Interstitial fluid pressure is often raised in solid tumors, which can limit drug delivery. To probe low-velocity flow in tumors resulting from raised interstitial fluid pressure, we developed a novel MRI technique named convection-MRI, which uses a phase-contrast acquisition with a dual-inversion vascular nulling preparation to separate intra- and extravascular flow. Here, we report the results of experiments in flow phantoms, numerical simulations, and tumor xenograft models to investigate the technical feasibility of convection-MRI. We observed a significant correlation between estimates of effective fluid pressure from convection-MRI with gold-standard, invasive measurements of interstitial fluid pressure in mouse models of human colorectal carcinoma. Our results show how convection-MRI can provide insights into the growth and responsiveness to vascular-targeting therapy in colorectal cancers. Significance: A noninvasive method for measuring low-velocity fluid flow caused by raised fluid pressure can be used to assess changes caused by therapy. Cancer Res; 78(7); 1859–72. ©2018 AACR.

The neurovascular mechanisms underpinning the local regulation of cerebral blood flow (CBF) and oxygen transport remain elusive. In this study we have combined novel in vivo imaging of cortical microvascular and mural cell architecture with mathematical modelling of blood flow and oxygen transport, to provide new insights into CBF regulation that would be inaccessible in a conventional experimental context. Our study implicates vasomotion of smooth muscle actin-covered vessels, rather than pericyte-covered capillaries, as the main mechanism for modulating tissue oxygenation. We also resolve seemingly paradoxical observations in the literature around reduced blood velocity in response to arteriolar constrictions and deduce the cause to be propagation of constrictions to upstream penetrating arterioles. We provide support for pericytes acting as signalling conduits for upstream smooth muscle activation, and erythrocyte deformation as a complementary regulatory mechanism. Finally, we caution against the use of blood velocity as a proxy measurement for flow. Our combined imaging-modelling platform complements conventional experimentation allowing cerebrovascular physiology to be probed in unprecedented detail.

Heterogeneity of lymphatic vessels during embryogenesis is critical for organ-specific lymphatic function. Little is known about lymphatics in the developing kidney, despite their established roles in pathology of the mature organ. We performed three-dimensional imaging to characterize lymphatic vessel formation in the mammalian embryonic kidney at single-cell resolution. In mouse, we visually and quantitatively assessed the development of kidney lymphatic vessels, remodeling from a ring-like anastomosis under the nascent renal pelvis; a site of VEGF-C expression, to form a patent vascular plexus. We identified a heterogenous population of lymphatic endothelial cell clusters in mouse and human embryonic kidneys. Exogenous VEGF-C expanded the lymphatic population in explanted mouse embryonic kidneys. Finally, we characterized complex kidney lymphatic abnormalities in a genetic mouse model of polycystic kidney disease. Our study provides novel insights into the development of kidney lymphatic vasculature; a system which likely has fundamental roles in renal development, physiology and disease., eLife digest In most organs in the body, fluid tends to build up in the spaces between cells, especially if the organs become inflamed. Each organ has a ‘waste disposal system’; a set of specialized tubes called lymphatic vessels, to clear away this excess fluid and keep a check on inflammation. Defects in these tubes have been linked to a wide range of diseases including heart attacks, obesity, dementia and cancer. The kidneys are responsible for filtering blood and balancing many of the body’s chemical processes. Polycystic kidney disease (PKD) is the most common genetic kidney disorder and it results in cysts filled with fluid building up in the kidney. The growth of cysts in PKD may be due to a problem with the lymphatic vessels. However, compared to other organs, how lymphatic vessels first form within the kidney and what they do is not well understood. Now, Jafree et al. have used three-dimensional imaging to study how lymphatic vessels form in the kidneys of mice and humans. The experiments showed that lymphatic vessels first appear when mouse kidneys are about half developed, and start to grow rapidly when the kidneys are thought to begin filtering blood. Clusters of cells that may help lymphatic vessels to grow were also found hidden deep within the kidneys of mouse embryos. Treating the kidneys with a factor that stimulates the growth of lymphatic vessels increased the numbers of these clusters. Jafree et al. found similar clusters of cells in human kidneys, suggesting that lymphatic vessels in the kidneys of different mammals may develop in the same way. Further experiments showed that the lymphatic vessels of kidneys in mice with PKD become distorted early on in the disease, when cysts are still small and before the mice develop symptoms. In the future, identifying drugs that target kidney lymphatic vessels may lead to more effective treatments for patients with PKD and other kidney diseases.

Cancers exhibit spatially heterogeneous, unique vascular architectures across individual samples, cell-lines and patients. This inherently disorganised collection of leaky blood vessels contribute significantly to suboptimal treatment efficacy. Preclinical tools are urgently required which incorporate the inherent variability and heterogeneity of tumours to optimise and engineer anti-cancer therapies. In this study, we present a novel computational framework which incorporates whole, realistic tumours extracted ex vivo to efficiently simulate vascular blood flow and interstitial fluid transport in silico for validation against in vivo biomedical imaging. Our model couples Poiseuille and Darcy descriptions of vascular and interstitial flow, respectively, and incorporates spatially heterogeneous blood vessel lumen and interstitial permeabilities to generate accurate predictions of tumour fluid dynamics. Our platform enables highly-controlled experiments to be performed which provide insight into how tumour vascular heterogeneity contributes to tumour fluid transport. We detail the application of our framework to an orthotopic murine glioma (GL261) and a human colorectal carcinoma (LS147T), and perform sensitivity analysis to gain an understanding of the key biological mechanisms which determine tumour fluid transport. Finally we mimic vascular normalization by modifying parameters, such as vascular and interstitial permeabilities, and show that incorporating realistic vasculatures is key to modelling the contrasting fluid dynamic response between tumour samples. Contrary to literature, we show that reducing tumour interstitial fluid pressure is not essential to increase interstitial perfusion and that therapies should seek to develop an interstitial fluid pressure gradient. We also hypothesise that stabilising vessel diameters and permeabilities are not key responses following vascular normalization and that therapy may alter interstitial hydraulic conductivity. Consequently, we suggest that normalizing the interstitial microenvironment may provide a more effective means to increase interstitial perfusion within tumours., Author summary The structure of tumours varies widely, with dense and chaotically-formed networks of blood vessels that differ between each individual tumour and even between different regions of the same tumour. This atypical environment can inhibit the delivery of anti-cancer therapies. Computational tools are urgently required which facilitate a deeper understanding of the relationship between blood vessel architectures and therapeutic response. We have developed a computational framework which integrates the complex tumour vascular architecture to predict fluid transport across all lengths scales in whole tumours. We apply our model to two tumour cell-lines and show that differences in their inherent vascular structures influence flow through cancerous tissue. We also use our platform to predict the fluid dynamic response following vascular normalization therapy in realistic, static tumour networks and show that the response is dependent on tumour vascular architecture. We hypothesise that therapy may alter the permeability of interstitial tissue to fluid transport and show that lowering interstitial fluid pressure is not a necessary therapeutic outcome to increase tumour perfusion.