Organ-wide 3D-imaging and topological analysis of the continuous microvascular network in a murine lymph node

Multiscale and multimodal imaging for three-dimensional vascular and histomorphological organ structure analysis of the pancreas

Exocrine and endocrine pancreas are interconnected anatomically and functionally, with vasculature facilitating bidirectional communication. Our understanding of this network remains limited, largely due to two-dimensional histology and missing combination with three-dimensional imaging. In this study, a multiscale 3D-imaging process was used to analyze a porcine pancreas. Clinical computed tomography, digital volume tomography, micro-computed tomography and Synchrotron-based propagation-based imaging were applied consecutively. Fields of view correlated inversely with attainable resolution from a whole organism level down to capillary structures with a voxel edge length of 2.0 µm. Segmented vascular networks from 3D-imaging data were correlated with tissue sections stained by immunohistochemistry and revealed highly vascularized regions to be intra-islet capillaries of islets of Langerhans. Generated 3D-datasets allowed for three-dimensional qualitative and quantitative organ and vessel structure analysis. Beyond this study, the method shows potential for application across a wide range of patho-morphology analyses and might possibly provide microstructural blueprints for biotissue engineering.

Subcapsular sinus macrophages maximize germinal center development in non-draining lymph nodes during blood-borne viral infection

Lymph node (LN) germinal centers (GCs) are critical sites for B cell activation and differentiation. GCs develop after specialized CD169+ macrophages residing in LN sinuses filter antigens (Ags) from the lymph and relay these Ags into proximal B cell follicles. Many viruses, however, first reach LNs through the blood during viremia (virus in the blood), rather than through lymph drainage from infected tissue. How LNs capture viral Ag from the blood to allow GC development is not known. Here, we followed Zika virus (ZIKV) dissemination in mice and subsequent GC formation in both infected tissue-draining and non-draining LNs. From the footpad, ZIKV initially disseminated through two LN chains, infecting LN macrophages and leading to GC formation. Despite rapid ZIKV viremia, non-draining LNs were not infected for several days. Non-draining LN infection correlated with virus-induced vascular leakage and neutralization of permeability reduced LN macrophage attrition. Depletion of non-draining LN macrophages significantly decreased GC B cells in these nodes. Thus, although LNs inefficiently captured viral Ag directly from the blood, GC formation in non-draining LNs proceeded similarly to draining LNs through LN sinus CD169+ macrophages. Together, our findings reveal a conserved pathway allowing LN macrophages to activate antiviral B cells in LNs distal from infected tissue after blood-borne viral infection.

Computational fluid dynamic modeling of the lymphatic system: a review of existing models and future directions

Historically, research into the lymphatic system has been overlooked due to both a lack of knowledge and limited recognition of its importance. In the last decade however, lymphatic research has gained substantial momentum and has included the development of a variety of computational models to aid understanding of this complex system. This article reviews existing computational fluid dynamic models of the lymphatics covering each structural component including the initial lymphatics, pre-collecting and collecting vessels, and lymph nodes. This is followed by a summary of limitations and gaps in existing computational models and reasons that development in this field has been hindered to date. Over the next decade, efforts to further characterize lymphatic anatomy and physiology are anticipated to provide key data to further inform and validate lymphatic fluid dynamic models. Development of more comprehensive multiscale- and multi-physics computational models has the potential to significantly enhance the understanding of lymphatic function in both health and disease.

Three-dimensional Imaging Reveals Immune-driven Tumor-associated High Endothelial Venules as a Key Correlate of Tumor Rejection Following Depletion of Regulatory T Cells

High endothelial venules (HEV) are specialized post capillary venules that recruit naïve T cells and B cells into secondary lymphoid organs (SLO) such as lymph nodes (LN). Expansion of HEV networks in SLOs occurs following immune activation to support development of an effective immune response. In this study, we used a carcinogen-induced model of fibrosarcoma to examine HEV remodeling after depletion of regulatory T cells (Treg). We used light sheet fluorescence microscopy imaging to visualize entire HEV networks, subsequently applying computational tools to enable topological mapping and extraction of numerical descriptors of the networks. While these analyses revealed profound cancer- and immune-driven alterations to HEV networks within LNs, these changes did not identify successful responses to treatment. The presence of HEV networks within tumors did however clearly distinguish responders from nonresponders. Finally, we show that a successful treatment response is dependent on coupling tumor-associated HEV (TA-HEV) development to T-cell activation implying that T-cell activation acts as the trigger for development of TA-HEVs which subsequently serve to amplify the immune response by facilitating extravasation of T cells into the tumor mass.

Three-dimensional analyses of vascular network morphology in a murine lymph node by X-ray phase-contrast tomography with a 2D Talbot array

With growing molecular evidence for correlations between spatial arrangement of blood vasculature and fundamental immunological functions, carried out in distinct compartments of the subdivided lymph node, there is an urgent need for three-dimensional models that can link these aspects. We reconstructed such models at a 1.84 µm resolution by the means of X-ray phase-contrast imaging with a 2D Talbot array in a short time without any staining. In addition reconstructions are verified in immunohistochemistry staining as well as in ultrastructural analyses. While conventional illustrations of mammalian lymph nodes depict the hilus as a definite point of blood and lymphatic vessel entry and exit, our method revealed that multiple branches enter and emerge from an area that extends up to one third of the organ's surface. This could be a prerequisite for the drastic and location-dependent remodeling of vascularization, which is necessary for lymph node expansion during inflammation. Contrary to corrosion cast studies we identified B-cell follicles exhibiting a two times denser capillary network than the deep cortical units of the T-cell zone. In addition to our observation of high endothelial venules spatially surrounding the follicles, this suggests a direct connection between morphology and B-cell homing. Our findings will deepen the understanding of functional lymph node composition and lymphocyte migration on a fundamental basis.

The potential role of ischaemia-reperfusion injury in chronic, relapsing diseases such as rheumatoid arthritis, Long COVID, and ME/CFS: evidence, mechanisms, and therapeutic implications

Ischaemia–reperfusion (I–R) injury, initiated via bursts of reactive oxygen species produced during the reoxygenation phase following hypoxia, is well known in a variety of acute circumstances. We argue here that I–R injury also underpins elements of the pathology of a variety of chronic, inflammatory diseases, including rheumatoid arthritis, ME/CFS and, our chief focus and most proximally, Long COVID. Ischaemia may be initiated via fibrin amyloid microclot blockage of capillaries, for instance as exercise is started; reperfusion is a necessary corollary when it finishes. We rehearse the mechanistic evidence for these occurrences here, in terms of their manifestation as oxidative stress, hyperinflammation, mast cell activation, the production of marker metabolites and related activities. Such microclot-based phenomena can explain both the breathlessness/fatigue and the post-exertional malaise that may be observed in these conditions, as well as many other observables. The recognition of these processes implies, mechanistically, that therapeutic benefit is potentially to be had from antioxidants, from anti-inflammatories, from iron chelators, and via suitable, safe fibrinolytics, and/or anti-clotting agents. We review the considerable existing evidence that is consistent with this, and with the biochemical mechanisms involved.

Quantification of vascular networks in photoacoustic mesoscopy

Mesoscopic photoacoustic imaging (PAI) enables non-invasive visualisation of tumour vasculature. The visual or semi-quantitative 2D measurements typically applied to mesoscopic PAI data fail to capture the 3D vessel network complexity and lack robust ground truths for assessment of accuracy. Here, we developed a pipeline for quantifying 3D vascular networks captured using mesoscopic PAI and tested the preservation of blood volume and network structure with topological data analysis. Ground truth data of in silico synthetic vasculatures and a string phantom indicated that learning-based segmentation best preserves vessel diameter and blood volume at depth, while rule-based segmentation with vesselness image filtering accurately preserved network structure in superficial vessels. Segmentation of vessels in breast cancer patient-derived xenografts (PDXs) compared favourably to ex vivo immunohistochemistry. Furthermore, our findings underscore the importance of validating segmentation methods when applying mesoscopic PAI as a tool to evaluate vascular networks in vivo.

Open-source analysis and visualization of segmented vasculature datasets with VesselVio

Vascular networks are fundamental components of biological systems. Quantitative analysis and observation of the features of these networks can improve our understanding of their roles in health and disease. Recent advancements in imaging technologies have enabled the generation of large-scale vasculature datasets, but barriers to analyzing these datasets remain. Modern analysis options are mainly limited to paid applications or open-source terminal-based software that requires programming knowledge with high learning curves. Here, we describe VesselVio, an open-source application developed to analyze and visualize pre-binarized vasculature datasets and pre-constructed vascular graphs. Vasculature datasets and graphs can be loaded with annotations and processed with custom parameters. Here, the program is tested on ground-truth datasets and is compared with current pipelines. The utility of VesselVio is demonstrated by the analysis of multiple formats of 2D and 3D datasets acquired with several imaging modalities, including annotated mouse whole-brain vasculature volumes.

Integrative lymph node-mimicking models created with biomaterials and computational tools to study the immune system

The lymph node (LN) is a vital organ of the lymphatic and immune system that enables timely detection, response, and clearance of harmful substances from the body. Each LN comprises of distinct substructures, which host a plethora of immune cell types working in tandem to coordinate complex innate and adaptive immune responses. An improved understanding of LN biology could facilitate treatment in LN-associated pathologies and immunotherapeutic interventions, yet at present, animal models, which often have poor physiological relevance, are the most popular experimental platforms. Emerging biomaterial engineering offers powerful alternatives, with the potential to circumvent limitations of animal models, for in-depth characterization and engineering of the lymphatic and adaptive immune system. In addition, mathematical and computational approaches, particularly in the current age of big data research, are reliable tools to verify and complement biomaterial works. In this review, we first discuss the importance of lymph node in immunity protection followed by recent advances using biomaterials to create in vitro/vivo LN-mimicking models to recreate the lymphoid tissue microstructure and microenvironment, as well as to describe the related immuno-functionality for biological investigation. We also explore the great potential of mathematical and computational models to serve as in silico supports. Furthermore, we suggest how both in vitro/vivo and in silico approaches can be integrated to strengthen basic patho-biological research, translational drug screening and clinical personalized therapies. We hope that this review will promote synergistic collaborations to accelerate progress of LN-mimicking systems to enhance understanding of immuno-complexity.

Intranodal delivery of modified docetaxel: Innovative therapeutic method to inhibit tumor cell growth in lymph nodes

Delivery of chemotherapeutic agents into metastatic lymph nodes (LNs) is challenging as they are unevenly distributed in the body. They are difficult to access via traditional systemic routes of drug administration, which produce significant adverse effects and result in low accumulation of drugs into the cancerous LN. To improve the survival rate of patients with LN metastasis, a lymphatic drug delivery system (LDDS) has been developed to target metastatic LN by delivering chemotherapy agents into sentinel LN (SLN) under ultrasound guidance. The LDDS is an advanced method that can be applied in the early stage of the progression of tumor cells in the SLN before tumor mass formation has occurred. Here we investigated the optimal physicochemical ranges of chemotherapeutic agents’ solvents with the aim of increasing treatment efficacy using the LDDS. We found that an appropriate osmotic pressure range for drug administration was 700–3,000 kPa, with a viscosity < 40 mPa⋅s. In these physicochemical ranges, expansion of lymphatic vessels and sinuses, drug retention, and subsequent antitumor effects could be more precisely controlled. Furthermore, the antitumor effects depended on the tumor progression stage in the SLN, the injection rate, and the volumes of administered drugs. We anticipate these optimal ranges to be a starting point for developing more effective drug regimens to treat metastatic LN with the LDDS.

Development of Stereo NIR-II Fluorescence Imaging System for 3D Tumor Vasculature in Small Animals

Near-infrared-II (NIR-II, 1000–1700 nm) fluorescence imaging boasts high spatial resolution and deep tissue penetration due to low light scattering, reduced photon absorption, and low tissue autofluorescence. NIR-II biological imaging is applied mainly in the noninvasive visualization of blood vessels and tumors in deep tissue. In the study, a stereo NIR-II fluorescence imaging system was developed for acquiring three-dimension (3D) images on tumor vasculature in real-time, on top of the development of fluorescent semiconducting polymer dots (IR-TPE Pdots) with ultra-bright NIR-II fluorescence (1000–1400 nm) and high stability to perform long-term fluorescence imaging. The NIR-II imaging system only consists of one InGaAs camera and a moving stage to simulate left-eye view and right-eye view for the construction of 3D in-depth blood vessel images. The system was validated with blood vessel phantom of tumor-bearing mice and was applied successfully in obtaining 3D blood vessel images with 0.6 mm- and 5 mm-depth resolution and 0.15 mm spatial resolution. The NIR-II stereo vision provides precise 3D information on the tumor microenvironment and blood vessel path.

Mathematical Modeling of Lymph Node Drainage Function by Neural Network

The lymph node (LN) represents a key structural component of the lymphatic system network responsible for the fluid balance in tissues and the immune system functioning. Playing an important role in providing the immune defense of the host organism, LNs can also contribute to the progression of pathological processes, e.g., the spreading of cancer cells. To gain a deeper understanding of the transport function of LNs, experimental approaches are used. Mathematical modeling of the fluid transport through the LN represents a complementary tool for studying the LN functioning under broadly varying physiological conditions. We developed an artificial neural network (NN) model to describe the lymph node drainage function. The NN model predicts the flow characteristics through the LN, including the exchange with the blood vascular systems in relation to the boundary and lymphodynamic conditions, such as the afferent lymph flow, Darcy’s law constants and Starling’s equation parameters. The model is formulated as a feedforward NN with one hidden layer. The NN complements the computational physics-based model of a stationary fluid flow through the LN and the fluid transport across the blood vessel system of the LN. The physical model is specified as a system of boundary integral equations (IEs) equivalent to the original partial differential equations (PDEs; Darcy’s Law and Starling’s equation) formulations. The IE model has been used to generate the training dataset for identifying the NN model architecture and parameters. The computation of the output LN drainage function characteristics (the fluid flow parameters and the exchange with blood) with the trained NN model required about 1000-fold less central processing unit (CPU) time than computationally tracing the flow characteristics of interest with the physics-based IE model. The use of the presented computational models will allow for a more realistic description and prediction of the immune cell circulation, cytokine distribution and drug pharmacokinetics in humans under various health and disease states as well as assisting in the development of artificial LN-on-a-chip technologies.

Characterizing perfusion defects in metastatic lymph nodes at an early stage using high-frequency ultrasound and micro-CT imaging

A perfusion defect in a metastatic lymph node (LN) can be visualized as a localized area of low contrast on contrast-enhanced CT, MRI or ultrasound images. Hypotheses for perfusion defects include abnormal hemodynamics in neovascular vessels or a decrease in blood flow in pre-existing blood vessels in the parenchyma due to compression by LN tumor growth. However, the mechanisms underlying perfusion defects in LNs during the early stage of LN metastasis have not been investigated. We show that tumor mass formation with very few microvessels was associated with a perfusion defect in a non-enlarged LN at the early stage of LN metastasis in a LN adenopathy mouse (LN size circa 10 mm). We found in a mouse model of LN metastasis, induced using non-keratinizing tumor cells, that during the formation of the perfusion defect in a non-enlarged LN, the number of blood vessels ≤ 50 μm in diameter decreased, while those of > 50 μm in diameter increased. The methods used were contrast-enhanced high-frequency ultrasound and contrast-enhanced micro-CT imaging systems, with a maximum spatial resolution of > 30 μm. Furthermore, we found no tumor angiogenesis or oxygen partial pressure (pO₂) changes in the metastatic LN. Our results demonstrate that the perfusion defect appears to be a specific form of tumorigenesis in the LN, which is a vascular-rich organ. We anticipate that a perfusion defect on ultrasound, CT or MRI images will be used as an indicator of a non-enlarged metastatic LN at an early stage.

Hierarchical imaging and computational analysis of three-dimensional vascular network architecture in the entire postnatal and adult mouse brain

The formation of new blood vessels and the establishment of vascular networks are crucial during brain development, in the adult healthy brain, as well as in various diseases of the central nervous system. Here, we describe a step-by-step protocol for our recently developed method that enables hierarchical imaging and computational analysis of vascular networks in postnatal and adult mouse brains. The different stages of the procedure include resin-based vascular corrosion casting, scanning electron microscopy, synchrotron radiation and desktop microcomputed tomography imaging, and computational network analysis. Combining these methods enables detailed visualization and quantification of the 3D brain vasculature. Network features such as vascular volume fraction, branch point density, vessel diameter, length, tortuosity and directionality as well as extravascular distance can be obtained at any developmental stage from the early postnatal to the adult brain. This approach can be used to provide a detailed morphological atlas of the entire mouse brain vasculature at both the postnatal and the adult stage of development. Our protocol allows the characterization of brain vascular networks separately for capillaries and noncapillaries. The entire protocol, from mouse perfusion to vessel network analysis, takes ~10 d.

A Straightforward Method for 3D Visualization of B Cell Clusters and High Endothelial Venules in Lymph Nodes Highlights Differential Roles of TNFRI and -II

Whole mount tissue immunolabeling and imaging of complete organs has tremendous benefits in characterizing organ morphology. Here, we present a straightforward method for immunostaining, clearing and imaging of whole murine peripheral lymph nodes (PLNs) for detailed analysis of their architecture and discuss all procedures in detail in a step-by-step approach. Given the importance of tumor necrosis factor receptor (TNFR) signaling in development of PLNs we used TNFRI^-/- and TNFRII^-/- mice models as proof-of-concept for this technique by visualizing and analyzing structural changes in PLN B cell clusters and high endothelial venules (HEVs). Samples were subjected to de- and rehydration with methanol, labeled with antibodies for B cells, T cells and high endothelial venules (HEVs) and optically cleared using benzyl alcohol-benzyl benzoate. Imaging was done using LaVision light sheet microscope and analysis with Imaris software. Using these techniques, we confirmed previous findings that TNFRI signaling is essential for formation of individual B cell clusters. In addition, Our data suggest that TNFRII signaling is also to some extent involved in this process as TNFRII-/- PLNs had a B cell cluster morphology reminiscent of TNFRI^-/- PLNs. Moreover, visualization and objective quantification of the complete PLN high endothelial vasculature unveiled reduced volume, length and branching points of HEVs in TNFRI^-/- PLNs, revealing an earlier unrecognized contribution of TNFRI signaling in HEV morphology. Together, these results underline the potential of whole mount tissue staining and advanced imaging techniques to unravel even subtle changes in lymphoid tissue architecture.

Study of the physicochemical properties of drugs suitable for administration using a lymphatic drug delivery system

Lymph node (LN) metastasis is thought to account for 20‐30% of deaths from head and neck cancer. The lymphatic drug delivery system (LDDS) is a new technology that enables the injection of drugs into a sentinel LN (SLN) during the early stage of tumor metastasis to treat the SLN and secondary metastatic LNs. However, the optimal physicochemical properties of the solvent used to carry the drug have not been determined. Here, we show that the osmotic pressure and viscosity of the solvent influenced the antitumor effect of cisplatin (CDDP) in a mouse model of LN metastasis. Tumor cells were inoculated into the proper axillary LN (PALN), and the LDDS was used to inject CDDP solution into the subiliac LN (SiLN) to treat the tumor cells in the downstream PALN. CDDP dissolved in saline had no therapeutic effects in the PALN after it was injected into the SiLN using the LDDS or into the tail vein (as a control). However, CDDP solution with an osmotic pressure of ~ 1,900 kPa and a viscosity of ~ 12 mPa⋅s suppressed tumor growth in the PALN after it was injected into the SiLN using the LDDS. The high osmotic pressure dilated the lymphatic vessels and sinuses to enhance drug flow in the PALN, and the high viscosity increased the retention of CDDP in the PALN. Our results demonstrate that optimizing the osmotic pressure and viscosity of the solvent can enhance the effects of CDDP, and possibly other anticancer drugs, after administration using the LDDS.

Lymph node stromal cell subsets-Emerging specialists for tailored tissue-specific immune responses

The effective priming of adaptive immune responses depends on the precise dispatching of lymphocytes and antigens into and within lymph nodes (LNs), which are strategically dispersed throughout the body. Over the past decade, a growing body of evidence has advanced our understanding of lymph node stromal cells (LNSCs) from viewing them as mere accessory cells to seeing them as critical cellular players for the modulation of adaptive immune responses. In this review, we summarize current advances on the pivotal roles that LNSCs play in orchestrating adaptive immune responses during homeostasis and infection, and highlight the imprinting of location-specific information by micro-environmental cues into LNSCs, thereby tailoring tissue-specific immune responses.

Angiogenesis in Lymph Nodes Is a Critical Regulator of Immune Response and Lymphoma Growth

Tumor-induced remodeling of the microenvironment in lymph nodes (LNs) includes the formation of blood vessels, which goes beyond the regulation of metabolism, and shaping a survival niche for tumor cells. In contrast to solid tumors, which primarily rely on neo-angiogenesis, hematopoietic malignancies usually grow within pre-vascularized autochthonous niches in secondary lymphatic organs or the bone marrow. The mechanisms of vascular remodeling in expanding LNs during infection-induced responses have been studied in more detail; in contrast, insights into the conditions of lymphoma growth and lodging remain enigmatic. Based on previous murine studies and clinical trials in human, we conclude that there is not a universal LN-specific angiogenic program applicable. Instead, signaling pathways that are tightly connected to autochthonous and infiltrating cell types contribute variably to LN vascular expansion. Inflammation related angiogenesis within LNs relies on dendritic cell derived pro-inflammatory cytokines stimulating vascular endothelial growth factor-A (VEGF-A) expression in fibroblastic reticular cells, which in turn triggers vessel growth. In high-grade B cell lymphoma, angiogenesis correlates with poor prognosis. Lymphoma cells immigrate and grow in LNs and provide pro-angiogenic growth factors themselves. In contrast to infectious stimuli that impact on LN vasculature, they do not trigger the typical inflammatory and hypoxia-related stroma-remodeling cascade. Blood vessels in LNs are unique in selective recruitment of lymphocytes via high endothelial venules (HEVs). The dissemination routes of neoplastic lymphocytes are usually disease stage dependent. Early seeding via the blood stream requires the expression of the homeostatic chemokine receptor CCR7 and of L-selectin, both cooperate to facilitate transmigration of tumor and also of protective tumor-reactive lymphocytes via HEV structures. In this view, the HEV route is not only relevant for lymphoma cell homing, but also for a continuous immunosurveillance. We envision that HEV functional and structural alterations during lymphomagenesis are not only key to vascular remodeling, but also impact on tumor cell accessibility when targeted by T cell-mediated immunotherapies.

Chemotaxing neutrophils enter alternate branches at capillary bifurcations

Upon tissue injury or microbial invasion, a large number of neutrophils converge from blood to the sites of injury or infection in a short time. The migration through a limited number of paths through tissues and capillary networks seems efficient and 'traffic jams' are generally avoided. However, the mechanisms that guide efficient trafficking of large numbers of neutrophils through capillary networks are not well understood. Here we show that pairs of neutrophils arriving closely one after another at capillary bifurcations migrate to alternating branches in vivo and in vitro. Perturbation of chemoattractant gradients and the increased hydraulic resistance induced by the first neutrophil in one branch biases the migration of the following neutrophil towards the other branch. These mechanisms guide neutrophils to efficiently navigate through capillary networks and outline the effect of inter-neutrophil interactions during migration on overall lymphocyte trafficking patterns in confined environments.

A hybrid discrete-continuum approach for modelling microcirculatory blood flow

In recent years, biological imaging techniques have advanced significantly and it is now possible to digitally reconstruct microvascular network structures in detail, identifying the smallest capillaries at sub-micron resolution and generating large 3D structural data sets of size >106 vessel segments. However, this relies on ex vivo imaging; corresponding in vivo measures of microvascular structure and flow are limited to larger branching vessels and are not achievable in three dimensions for the smallest vessels. This suggests the use of computational modelling to combine in vivo measures of branching vessel architecture and flows with ex vivo data on complete microvascular structures to predict effective flow and pressures distributions. In this paper, a hybrid discrete-continuum model to predict microcirculatory blood flow based on structural information is developed and compared with existing models for flow and pressure in individual vessels. A continuum-based Darcy model for transport in the capillary bed is coupled via point sources of flux to flows in individual arteriolar vessels, which are described explicitly using Poiseuille's law. The venular drainage is represented as a spatially uniform flow sink. The resulting discrete-continuum framework is parameterized using structural data from the capillary network and compared with a fully discrete flow and pressure solution in three networks derived from observations of the rat mesentery. The discrete-continuum approach is feasible and effective, providing a promising tool for extracting functional transport properties in situations where vascular branching structures are well defined.

From whole-organ imaging to in-silico blood flow modeling: A new multi-scale network analysis for revisiting tissue functional anatomy

We present a multi-disciplinary image-based blood flow perfusion modeling of a whole organ vascular network for analyzing both its structural and functional properties. We show how the use of Light-Sheet Fluorescence Microscopy (LSFM) permits whole-organ micro-vascular imaging, analysis and modelling. By using adapted image post-treatment workflow, we could segment, vectorize and reconstruct the entire micro-vascular network composed of 1.7 million vessels, from the tissue-scale, inside a ∼ 25 × 5 × 1 = 125mm3 volume of the mouse fat pad, hundreds of times larger than previous studies, down to the cellular scale at micron resolution, with the entire blood perfusion modeled. Adapted network analysis revealed the structural and functional organization of meso-scale tissue as strongly connected communities of vessels. These communities share a distinct heterogeneous core region and a more homogeneous peripheral region, consistently with known biological functions of fat tissue. Graph clustering analysis also revealed two distinct robust meso-scale typical sizes (from 10 to several hundred times the cellular size), revealing, for the first time, strongly connected functional vascular communities. These community networks support heterogeneous micro-environments. This work provides the proof of concept that in-silico all-tissue perfusion modeling can reveal new structural and functional exchanges between micro-regions in tissues, found from community clusters in the vascular graph.

High-resolution 3D imaging and topological mapping of the lymph node conduit system

The conduit network is a hallmark of lymph node microanatomy, but lack of suitable imaging technology has prevented comprehensive investigation of its topology. We employed an extended-volume imaging system to capture the conduit network of an entire murine lymph node (comprising over 280,000 segments). The extensive 3D images provide a comprehensive overview of the regions supplied by conduits, including perivascular sleeves and distinctive “follicular reservoirs” within B cell follicles, surrounding follicular dendritic cells. A 3D topology map of conduits within the T-cell zone showed homogeneous branching, but conduit density was significantly higher in the superficial T-cell zone compared with the deep zone, where distances between segments are sufficient for T cells to lose contact with fibroblastic reticular cells. This topological mapping of the conduit anatomy can now aid modeling of its roles in lymph node function, as we demonstrate by simulating T-cell motility in the different T-cell zones.

Extended-volume confocal imaging allowed 3D visualisation of the fine network of conduits within lymph nodes; the resulting map of conduit topology underscores structural differences between the deep and superficial T cell zone and identifies "follicular reservoirs" within B cell follicles that concentrate lymphoid fluid around follicular dendritic cells.

Quantification of the Whole Lymph Node Vasculature Based on Tomography of the Vessel Corrosion Casts

Lymph nodes (LN) are crucial for immune function, and comprise an important interface between the blood and lymphatic systems. Blood vessels (BV) in LN are highly specialized, featuring high endothelial venules across which most of the resident lymphocytes crossed. Previous measurements of overall lymph and BV flow rates demonstrated that fluid also crosses BV walls, and that this is important for immune function. However, the spatial distribution of the BV in LN has not been quantified to the degree necessary to analyse the distribution of transmural fluid movement. In this study, we seek to quantify the spatial localization of LNBV, and to predict fluid movement across BV walls. MicroCT imaging of murine popliteal LN showed that capillaries were responsible for approximately 75% of the BV wall surface area, and that this was mostly distributed around the periphery of the node. We then modelled blood flow through the BV to obtain spatially resolved hydrostatic pressures, which were then combined with Starling’s law to predict transmural flow. Much of the total 10 nL/min transmural flow (under normal conditions) was concentrated in the periphery, corresponding closely with surface area distribution. These results provide important insights into the inner workings of LN, and provide a basis for further exploration of the role of LN flow patterns in normal and pathological functions.

Engineering the vasculature for islet transplantation

The microvasculature in the pancreatic islet is highly specialized for glucose sensing and insulin secretion. Although pancreatic islet transplantation is a potentially life-changing treatment for patients with insulin-dependent diabetes, a lack of blood perfusion reduces viability and function of newly transplanted tissues. Functional vasculature around an implant is not only necessary for the supply of oxygen and nutrients but also required for rapid insulin release kinetics and removal of metabolic waste. Inadequate vascularization is particularly a challenge in islet encapsulation. Selectively permeable membranes increase the barrier to diffusion and often elicit a foreign body reaction including a fibrotic capsule that is not well vascularized. Therefore, approaches that aid in the rapid formation of a mature and robust vasculature in close proximity to the transplanted cells are crucial for successful islet transplantation or other cellular therapies. In this paper, we review various strategies to engineer vasculature for islet transplantation. We consider properties of materials (both synthetic and naturally derived), prevascularization, local release of proangiogenic factors, and co-transplantation of vascular cells that have all been harnessed to increase vasculature. We then discuss the various other challenges in engineering mature, long-term functional and clinically viable vasculature as well as some emerging technologies developed to address them. The benefits of physiological glucose control for patients and the healthcare system demand vigorous pursuit of solutions to cell transplant challenges. STATEMENT OF SIGNIFICANCE: Insulin-dependent diabetes affects more than 1.25 million people in the United States alone. Pancreatic islets secrete insulin and other endocrine hormones that control glucose to normal levels. During preparation for transplantation, the specialized islet blood vessel supply is lost. Furthermore, in the case of cell encapsulation, cells are protected within a device, further limiting delivery of nutrients and absorption of hormones. To overcome these issues, this review considers methods to rapidly vascularize sites and implants through material properties, pre-vascularization, delivery of growth factors, or co-transplantation of vessel supporting cells. Other challenges and emerging technologies are also discussed. Proper vascular growth is a significant component of successful islet transplantation, a treatment that can provide life-changing benefits to patients.

Glioblastoma multiforme restructures the topological connectivity of cerebrovascular networks

Glioblastoma multiforme alters healthy tissue vasculature by inducing angiogenesis and vascular remodeling. To fully comprehend the structural and functional properties of the resulting vascular network, it needs to be studied collectively by considering both geometric and topological properties. Utilizing Single Plane Illumination Microscopy (SPIM), the detailed capillary structure in entire healthy and tumor-bearing mouse brains could be resolved in three dimensions. At the scale of the smallest capillaries, the entire vascular systems of bulk U87- and GL261-glioblastoma xenografts, their respective cores, and healthy brain hemispheres were modeled as complex networks and quantified with fundamental topological measures. All individual vessel segments were further quantified geometrically and modular clusters were uncovered and characterized as meta-networks, facilitating an analysis of large-scale connectivity. An inclusive comparison of large tissue sections revealed that geometric properties of individual vessels were altered in glioblastoma in a relatively subtle way, with high intra- and inter-tumor heterogeneity, compared to the impact on the vessel connectivity. A network topology analysis revealed a clear decomposition of large modular structures and hierarchical network organization, while preserving most fundamental topological classifications, in both tumor models with distinct growth patterns. These results augment our understanding of cerebrovascular networks and offer a topological assessment of glioma-induced vascular remodeling. The findings may help understand the emergence of hypoxia and necrosis, and prove valuable for therapeutic interventions such as radiation or antiangiogenic therapy.

Microfluidic blood vasculature replicas using backside lithography

Blood vessels in living tissues are an organized and hierarchical network of arteries, arterioles, capillaries, veinules and veins. Their sizes, lengths, shapes and connectivity are set up for an optimum perfusion of the tissues in which they deploy. In order to study the hemodynamics and hemophysics of blood flows and also to investigate artificial vasculature for organs on a chip, it is essential to reproduce most of these geometric features. Common microfluidic techniques produce channels with a uniform height and a rectangular cross section that do not capture the size hierarchy observed in vivo. This paper presents a new single-mask photolithography process using an optical diffuser to produce a backside exposure leading to microchannels with both a rounded cross section and a direct proportionality between local height and local width, allowing a one-step design of intrinsically hierarchical networks.

3D Super-Resolution US Imaging of Rabbit Lymph Node Vasculature in Vivo by Using Microbubbles

Background Variations in lymph node (LN) microcirculation can be indicative of metastasis. The identification and quantification of metastatic LNs remains essential for prognosis and treatment planning, but a reliable noninvasive imaging technique is lacking. Three-dimensional super-resolution (SR) US has shown potential to noninvasively visualize microvascular networks in vivo. Purpose To study the feasibility of three-dimensional SR US imaging of rabbit LN microvascular structure and blood flow by using microbubbles. Materials and Methods In vivo studies were carried out to image popliteal LNs of two healthy male New Zealand white rabbits aged 6-8 weeks. Three-dimensional, high-frame-rate, contrast material-enhanced US was achieved by mechanically scanning with a linear imaging probe. Individual microbubbles were identified, localized, and tracked to form three-dimensional SR images and super-resolved velocity maps. Acoustic subaperture processing was used to improve image contrast and to generate enhanced power Doppler and color Doppler images. Vessel size and blood flow velocity distributions were evaluated and assessed by using Student paired t test. Results SR images revealed microvessels in the rabbit LN, with branches clearly resolved when separated by 30 µm, which is less than half of the acoustic wavelength and not resolvable by using power or color Doppler. The apparent size distribution of most vessels in the SR images was below 80 µm and agrees with micro-CT data, whereas most of those detected with Doppler techniques were larger than 80 µm in the images. The blood flow velocity distribution indicated that most of the blood flow in rabbit popliteal LN was at velocities lower than 5 mm/sec. Conclusion Three-dimensional super-resolution US imaging using microbubbles allows noninvasive nonionizing visualization and quantification of lymph node microvascular structures and blood flow dynamics with resolution below the wave diffraction limit. This technology has potential for studying the physiologic functions of the lymph system and for clinical detection of lymph node metastasis. Published under a CC BY 4.0 license. Online supplemental material is available for this article.

Call for a Multidisciplinary Effort to Map the Lymphatic System with Advanced Medical Imaging Techniques: A Review of the Literature and Suggestions for Future Anatomical Research

This review intends to rekindle efforts to map the lymphatic system by using a more modern approach, based on medical imaging. The structure, function, and pathologies associated with the lymphatic system are first discussed to highlight the need for more accurately mapping the lymphatic system. Next, the need for an interdisciplinary approach, with a central role for the anatomist, to come up with better maps of the lymphatic system is emphasized. The current approaches on lymphatic system research involving medical imaging will be discussed and suggestions will be made for an all-encompassing effort to thoroughly map the entire lymphatic system. A first-hand account of our integration as anatomists in the radiotherapy department is given as an example of interdisciplinary collaboration. From this account, it will become clear that the interdisciplinary collaboration of anatomists in the clinical disciplines involved in lymphatic system research/treatment still holds great promise in terms of improving clinical regimens that are currently being employed. As such, we hope that our fellow anatomists will join us in an interdisciplinary effort to map the lymphatic system, because this could, in a relatively short timeframe, provide improved treatment options for patients with cancer or lymphatic pathologies all over the world. Anat Rec, 302:1681-1695, 2019. © 2019 American Association for Anatomy.

Computer-assembled cross-species/cross-modalities two-pore physiologically based pharmacokinetic model for biologics in mice and rats

Two-pore physiologically-based pharmacokinetic (PBPK) models can be expected to describe the tissue distribution and elimination kinetics of soluble proteins, endogenous or dosed, as function of their size. In this work, we amalgamated our previous two-pore PBPK model for an inert domain antibody (dAb) in mice with the cross-species platform PBPK model for monoclonal antibodies described in literature into a unified two-pore platform that describes protein modalities of different sizes and includes neonatal Fc receptor (FcRn) mediated recycling. This unified PBPK model was parametrized for organ-specific lymph flow rates and the endosomal recycling rate constant using an extended tissue distribution time-course dataset that included an inert dAb, albumin and IgG in rats and mice. The model was evaluated by comparing the ab initio predictions for the tissue distribution and elimination properties of albumin-binding dAbs (AlbudAbs^TM) in mice and rats with the experimental observations. Due to the large number of molecular species and reactions involved in large-scale PBPK models, we have also developed and deployed a Matlab^TM script for automating the assembly of SimBiology^TM-based two-pore biologics PBPK models which drastically cuts the time and effort required for model building.

Treatment of false-negative metastatic lymph nodes by a lymphatic drug delivery system with 5-fluorouracil

Metastatic lymph nodes (LNs) may be the origin of systemic metastases. It will be important to develop a strategy that prevents systemic metastasis by treating these LNs at an early stage. False‐negative metastatic LNs, which are found during the early stage of metastasis development, are those that contain tumor cells but have a size and shape similar to LNs that do not host tumor cells. Here, we show that 5‐fluorouracil (5‐FU), delivered by means of a novel lymphatic drug delivery system (LDDS), can treat LNs with false‐negative metastases in a mouse model. The effects of 5‐FU on four cell lines were investigated using in vitro cytotoxicity and cell survival assays. The therapeutic effects of LDDS‐administered 5‐FU on false‐negative metastatic LNs were evaluated using bioluminescence imaging, high‐frequency ultrasound (US), and histology in MHX10/Mo‐lpr/lpr mice. These experimental animals develop LNs that are similar in size to human LNs. We found that all cell lines showed sensitivity to 5‐FU in the in vitro assays. Furthermore, a concentration‐dependent effect of 5‐FU to inhibit tumor growth was observed in tumor cells with low invasive growth characteristics, although a significant reduction in metastatic LN volume was not detected in MHX10/Mo‐lpr/lpr mice. Adverse effects of 5‐FU were not detected. 5‐Fluorouracil administration with a LDDS is an effective treatment method for false‐negative metastatic LNs. We anticipate that the delivery of anticancer drugs by a LDDS will be of great benefit in the prevention and treatment of cancer metastasis via LNs.

Integrative Computational Modeling of the Lymph Node Stromal Cell Landscape

Adaptive immune responses develop in secondary lymphoid organs such as lymph nodes (LNs) in a well-coordinated series of interactions between migrating immune cells and resident stromal cells. Although many processes that occur in LNs are well understood from an immunological point of view, our understanding of the fundamental organization and mechanisms that drive these processes is still incomplete. The aim of systems biology approaches is to unravel the complexity of biological systems and describe emergent properties that arise from interactions between individual constituents of the system. The immune system is greater than the sum of its parts, as is the case with any sufficiently complex system. Here, we review recent work and developments of computational LN models with focus on the structure and organization of the stromal cells. We explore various mathematical studies of intranodal T cell motility and migration, their interactions with the LN-resident stromal cells, and computational models of functional chemokine gradient fields and lymph flow dynamics. Lastly, we discuss briefly the importance of hybrid and multi-scale modeling approaches in immunology and the technical challenges involved.

Probe-based three-dimensional confocal laser endomicroscopy of brain tumors: technical note

Background

Confocal laser endomicroscopy (CLE) is used during fluorescence-guided brain tumor surgery for intraoperative microscopy of tumor tissue with cellular resolution. CLE could augment and expedite intraoperative decision-making and potentially aid in diagnosis and removal of tumor tissue.

Objective

To describe an extension of CLE imaging modality that produces Z-stack images and three-dimensional (3D) pseudocolored volumetric images.

Materials and methods

Hand-held probe-based CLE was used to collect images from GL261-luc2 gliomas in C57BL/6 mice and from human brain tumor biopsies. The mice were injected with fluorescein sodium (FNa) before imaging. Patients received FNa intraoperatively, and biopsies were imaged immediately in the operating room. Some specimens were counterstained with acridine orange, acriflavine, or Hoechst and imaged on a benchtop confocal microscope. CLE images at various depths were acquired automatically, compiled, rendered into 3D volumes using Fiji software and reviewed by a neuropathologist and neurosurgeons.

Results

CLE imaging, Z-stack acquisition, and 3D image rendering were performed using 19 mouse gliomas and 31 human tumors, including meningiomas, gliomas, and pituitary adenomas. Volumetric images and Z-stacks provided additional information about fluorescence signal distribution, cytoarchitecture, and the course of abnormal vasculature.

Conclusion

3D and Z-stack CLE imaging is a unique new option for live intraoperative endomicroscopy of brain tumors. The 3D images afford an increased spatial understanding of tumor cellular architecture and visualization of related structures compared with two-dimensional images. Future application of specific fluorescent probes could benefit from this rapid in vivo imaging technology for interrogation of brain tumor tissue.

Differential DARC/ACKR1 expression distinguishes venular from non-venular endothelial cells in murine tissues

BACKGROUND

Intravascular leukocyte recruitment in most vertebrate tissues is restricted to postcapillary and collecting venules, whereas capillaries and arterioles usually support little or no leukocyte adhesion. This segmental restriction is thought to be mediated by endothelial, rather than hemodynamic, differences. The underlying mechanisms are largely unknown, in part because effective tools to distinguish, isolate, and analyze venular endothelial cells (V-ECs) and non-venular endothelial cells (NV-ECs) have been unavailable. We hypothesized that the atypical chemokine receptor DARC (Duffy Antigen Receptor for Chemokines, a.k.a. ACKR1 or CD234) may distinguish V-ECs versus NV-ECs in mice.

METHODS

We generated a rat-anti-mouse monoclonal antibody (MAb) that specifically recognizes the erythroid and endothelial forms of native, surface-expressed DARC. Using this reagent, we characterized DARC expression and distribution in the microvasculature of murine tissues.

RESULTS

DARC was exquisitely restricted to post-capillary and small collecting venules and completely absent from arteries, arterioles, capillaries, veins, and most lymphatics in every tissue analyzed. Accordingly, intravital microscopy showed that adhesive leukocyte-endothelial interactions were restricted to DARC+ venules. DARC was detectable over the entire circumference of V-ECs, but was more concentrated at cell-cell junctions. Analysis of single-cell suspensions suggested that the frequency of V-ECs among the total microvascular EC pool varies considerably between different tissues.

CONCLUSIONS

Immunostaining of endothelial DARC allows the identification and isolation of intact V-ECs from multiple murine tissues. This strategy may be useful to dissect the mechanisms underlying segmental microvascular specialization in healthy and diseased tissues and to characterize the role of EC subsets in tissue-homeostasis, immune surveillance, infection, inflammation, and malignancies.

Topological Small-World Organization of the Fibroblastic Reticular Cell Network Determines Lymph Node Functionality

Fibroblastic reticular cells (FRCs) form the cellular scaffold of lymph nodes (LNs) and establish distinct microenvironmental niches to provide key molecules that drive innate and adaptive immune responses and control immune regulatory processes. Here, we have used a graph theory-based systems biology approach to determine topological properties and robustness of the LN FRC network in mice. We found that the FRC network exhibits an imprinted small-world topology that is fully regenerated within 4 wk after complete FRC ablation. Moreover, in silico perturbation analysis and in vivo validation revealed that LNs can tolerate a loss of approximately 50% of their FRCs without substantial impairment of immune cell recruitment, intranodal T cell migration, and dendritic cell-mediated activation of antiviral CD8+ T cells. Overall, our study reveals the high topological robustness of the FRC network and the critical role of the network integrity for the activation of adaptive immune responses.