Image-based modelling of skeletal muscle oxygenation

Histological analysis of the medial gastrocnemius muscle in young healthy children

Introduction: Histological data on muscle fiber size and proportion in (very) young typically developing (TD) children is not well documented and data on capillarization and satellite cell content are also lacking. Aims: This study investigated the microscopic properties of the medial gastrocnemius muscle in growing TD children, grouped according to age and gender to provide normal reference values in healthy children. Methods: Microbiopsies of the medial gastrocnemius (MG) muscle were collected in 46 TD boys and girls aged 2-10 years subdivided into 4 age groups (2-4, 4-6, 6-8 and 8-10 years). Sections were immunostained to assess fiber type cross-sectional area (fCSA) and proportion, the number of satellite cells (SC), capillary to fiber ratio (C/F), capillary density for type I and II fiber (CFD), capillary domain, capillary-to-fiber perimeter exchange index (CFPE) and heterogeneity index. fCSA was normalized to fibula length2 and the coefficient of variation (CV) was calculated to reflect fCSA intrasubject variability. Results: Absolute fCSA of all fibers increased with age (r = 0.72, p < 0.001) but more in boys (+112%, p < 0.05) than in girls (+48%, p > 0.05) Normalized fCSA, CV and fiber proportion did not differ between age groups and gender. C/F was strongly correlated with age in boys (r = 0.83, p < 0.001), and to a lesser extent in girls (r = 0.37, p = 0.115), while other capillary parameters as well as the number of SC remained stable with increasing age in boys and girls. Discussion: This study provides reference values of histological measures in MG according to age in normally growing boys and girls. These data may be used as a reference to determine disease impact and efficacy of therapeutic approach on the muscle.

Synchronized research on endothelial dysfunction and microcirculation structure in dorsal skin of rats with type 2 diabetes mellitus

The aim of this study was to explore changes in the microvascular tone as measured by laser Doppler flowmetry (LDF) and the microcirculation structure of the dorsal skin of rats with type 2 diabetes mellitus. The diabetic rat model was induced by a diet of high-sugar and high-lipid fodder combined with the injection of streptozotocin into the abdominal cavity. Depending on the interval between the development of diabetes and the experiments, the diabetic rats were subdivided into three groups. The evaluation of microvascular tone was based on the amplitude responses of the LDF signal fluctuations in the appropriate frequency range in the dorsal skin of the rats during a thermal test (at 42 °C). The nitric oxide (NO) level in plasma was also used as a marker of endothelial dysfunction. Changes in the microcirculation structure in the diabetic rats were estimated by measuring the microvascular density in the choke vessels of the dorsal skin of the rats. The experimental results with respect to red blood cell (RBC)-related parameters showed decreased hematocrit and hemoglobin levels and increased standard deviation of the width of the RBC distribution in three diabetic rats. The increasing fluctuation amplitudes diminished in the endothelial frequency range in response to the thermal test and this was accompanied by abnormal NO levels in plasma of the diabetic groups as compared with healthy rats. A significant reduction in the microvascular density of the choke vessels of the dorsal skin was found only in the diabetic group at the most advanced stage of diabetes in this experiment. Thus, we suggest that endothelial dysfunction occurs in diabetic rats and changes in the microcirculation structure of the dorsal skin occur in a later stage of diabetes development. A. Photograph of measurement method by using a LDF probe and heating device in the dorsal skin of the rat. B. Dorsal skin LDF signals of a healthy rat during the thermal stimuli test. (a) Blood flow signal record for the test. Wavelet filtration of blood flow signal in (b) myogenic range, (c) neurogenic range, and (d) endothelial range.

The capillary fascicle in skeletal muscle: Structural and functional physiology of RBC distribution in capillary networks

KEY POINTS

The capillary module, consisting of parallel capillaries from arteriole to venule, is classically considered as the building block of complex capillary networks. In skeletal muscle, this structure fails to address how blood flow is regulated along the entire length of the synchronously contracting muscle fibres. Using intravital video microscopy of resting extensor digitorum longus muscle in rats, we demonstrated the capillary fascicle as a series of interconnected modules forming continuous columns that align naturally with the dimensions of the muscle fascicle. We observed structural heterogeneity for module topology, and functional heterogeneity in space and time for capillary-red blood cell (RBC) haemodynamics within a module and between modules. We found that module RBC haemodynamics were independent of module resistance, providing direct evidence for microvascular flow regulation at the level of the capillary module. The capillary fascicle is an updated paradigm for characterizing blood flow and RBC distribution in skeletal muscle capillary networks.

ABSTRACT

Capillary networks are the fundamental site of oxygen exchange in the microcirculation. The capillary module (CM), consisting of parallel capillaries from terminal arteriole (TA) to post-capillary venule (PCV), is classically considered as the building block of complex capillary networks. In skeletal muscle, this structure fails to address how blood flow is regulated along the entire length of the synchronously contracting muscle fibres, requiring co-ordination from numerous modules. It has previously been recognized that TAs and PCVs interact with multiple CMs, creating interconnected networks. Using label-free intravital video microscopy of resting extensor digitorum longus muscle in rats, we found that these networks form continuous columns of linked CMs spanning thousands of microns, herein denoted as the capillary fascicle (CF); this structure aligns naturally with the dimensions of the muscle fascicle. We measured capillary-red blood cell (RBC) haemodynamics and module topology (n = 9 networks, 327 modules, 1491 capillary segments). The average module had length 481 μm, width 157 μm and 9.51 parallel capillaries. We observed structural heterogeneity for CM topology, and functional heterogeneity in space and time for capillary-RBC haemodynamics within a module and between modules. There was no correlation between capillary RBC velocity and lineal density. A passive inverse relationship between module length and haemodynamics was remarkably absent, providing direct evidence for microvascular flow regulation at the level of the CM. In summary, the CF is an updated paradigm for characterizing RBC distribution in skeletal muscle, and strengthens the theory of capillary networks as major contributors to the signal that regulates capillary perfusion.

Effect of adipose tissue thickness and tissue optical properties on the differential pathlength factor estimation for NIRS studies on human skeletal muscle

We propose a quantitative and systematic investigation of the differential pathlength factor (DPF) behavior for skeletal muscles and its dependence on different factors, such as the subcutaneous adipose tissue thickness (ATT), the variations of the tissue absorption (µ_a ) and reduced scattering (µ'_s ) coefficients, and the source-detector distance. A time domain (TD) NIRS simulation study is performed in a two-layer geometry mimicking a human skeletal muscle with an overlying adipose tissue layer. The DPF decreases when µ_a increases, while it increases when µ'_s increases. Moreover, a positive correlation between DPF and ATT is found. These results are supported by an in-vivo TD NIRS study on vastus lateralis and biceps brachii muscles of eleven subjects at rest, showing a high inter-subject and inter-muscle variability.

Tracing the footsteps of autophagy in computational biology

Autophagy plays a crucial role in maintaining cellular homeostasis through the degradation of unwanted materials like damaged mitochondria and misfolded proteins. However, the contribution of autophagy toward a healthy cell environment is not only limited to the cleaning process. It also assists in protein synthesis when the system lacks the amino acids’ inflow from the extracellular environment due to diet consumptions. Reduction in the autophagy process is associated with diseases like cancer, diabetes, non-alcoholic steatohepatitis, etc., while uncontrolled autophagy may facilitate cell death. We need a better understanding of the autophagy processes and their regulatory mechanisms at various levels (molecules, cells, tissues). This demands a thorough understanding of the system with the help of mathematical and computational tools. The present review illuminates how systems biology approaches are being used for the study of the autophagy process. A comprehensive insight is provided on the application of computational methods involving mathematical modeling and network analysis in the autophagy process. Various mathematical models based on the system of differential equations for studying autophagy are covered here. We have also highlighted the significance of network analysis and machine learning in capturing the core regulatory machinery governing the autophagy process. We explored the available autophagic databases and related resources along with their attributes that are useful in investigating autophagy through computational methods. We conclude the article addressing the potential future perspective in this area, which might provide a more in-depth insight into the dynamics of autophagy.

Correlative 3D Imaging and Microfluidic Modelling of Human Pulmonary Lymphatics using Immunohistochemistry and High-resolution μCT

Lung lymphatics maintain fluid homoeostasis by providing a drainage system that returns fluid, cells and metabolites to the circulatory system. The 3D structure of the human pulmonary lymphatic network is essential to lung function, but it is poorly characterised. Image-based 3D mathematical modelling of pulmonary lymphatic microfluidics has been limited by the lack of accurate and representative image geometries. This is due to the microstructural similarity of the lymphatics to the blood vessel network, the lack of lymphatic-specific biomarkers, the technical limitations associated with image resolution in 3D, and sectioning artefacts present in 2D techniques. We present a method that combines lymphatic specific (D240 antibody) immunohistochemistry (IHC), optimised high-resolution X-ray microfocus computed tomography (μCT) and finite-element mathematical modelling to assess the function of human peripheral lung tissue. The initial results identify lymphatic heterogeneity within and between lung tissue. Lymphatic vessel volume fraction and fractal dimension significantly decreases away from the lung pleural surface (p < 0.001, n = 25 and p < 0.01, n = 20, respectively). Microfluidic modelling successfully shows that in lung tissue the fluid derived from the blood vessels drains through the interstitium into the lymphatic vessel network and this drainage is different in the subpleural space compared to the intralobular space. When comparing lung tissue from health and disease, human pulmonary lymphatics were significantly different across five morphometric measures used in this study (p ≤ 0.0001). This proof of principle study establishes a new engineering technology and workflow for further studies of pulmonary lymphatics and demonstrates for the first time the combination of correlative μCT and IHC to enable 3D mathematical modelling of human lung microfluidics at micrometre resolution.

Tumor Ensemble-Based Modeling and Visualization of Emergent Angiogenic Heterogeneity in Breast Cancer

There is a critical need for new tools to investigate the spatio-temporal heterogeneity and phenotypic alterations that arise in the tumor microenvironment. However, computational investigations of emergent inter- and intra-tumor angiogenic heterogeneity necessitate 3D microvascular data from 'whole-tumors' as well as "ensembles" of tumors. Until recently, technical limitations such as 3D imaging capabilities, computational power and cost precluded the incorporation of whole-tumor microvascular data in computational models. Here, we describe a novel computational approach based on multimodality, 3D whole-tumor imaging data acquired from eight orthotopic breast tumor xenografts (i.e. a tumor 'ensemble'). We assessed the heterogeneous angiogenic landscape from the microvascular to tumor ensemble scale in terms of vascular morphology, emergent hemodynamics and intravascular oxygenation. We demonstrate how the abnormal organization and hemodynamics of the tumor microvasculature give rise to unique microvascular niches within the tumor and contribute to inter- and intra-tumor heterogeneity. These tumor ensemble-based simulations together with unique data visualization approaches establish the foundation of a novel 'cancer atlas' for investigators to develop their own in silico systems biology applications. We expect this hybrid image-based modeling framework to be adaptable for the study of other tissues (e.g. brain, heart) and other vasculature-dependent diseases (e.g. stroke, myocardial infarction).

Integrated method for quantitative morphometry and oxygen transport modeling in striated muscle

Identifying structural limitations in O2 transport is primarily restricted by current methods employed to characterize the nature of physiological remodeling. Inadequate resolution or breadth of available data has impaired development of routine diagnostic protocols and effective therapeutic strategies. Understanding O2 transport within striated muscle faces major challenges, most notably in quantifying how well individual fibers are supplied by the microcirculation, which has necessitated exploring tissue O2 supply using theoretical modeling of diffusive exchange. With capillary domains identified as a suitable model for the description of local O2 supply and requiring less computation than numerically calculating the trapping regions that are supplied by each capillary via biophysical transport models, we sought to design a high-throughput method for histological analysis. We present an integrated package that identifies optimal protocols for identification of important input elements, processing of digitized images with semiautomated routines, and incorporation of these data into a mathematical modeling framework with computed output visualized as the tissue partial pressure of O2 (Po2) distribution across a biopsy sample. Worked examples are provided using muscle samples from experiments involving rats and humans.

NEW & NOTEWORTHY Progress in quantitative morphometry and analytical modeling has tended to develop independently. Real diagnostic power lies in harnessing both disciplines within one user-friendly package. We present a semiautomated, high-throughput tool for determining muscle phenotype from biopsy material, which also provides anatomically relevant input to quantify tissue oxygenation, in a coherent package not previously available to nonspecialist investigators.

Methods to label, image, and analyze the complex structural architectures of microvascular networks

Microvascular networks play key roles in oxygen transport and nutrient delivery to meet the varied and dynamic metabolic needs of different tissues throughout the body, and their spatial architectures of interconnected blood vessel segments are highly complex. Moreover, functional adaptations of the microcirculation enabled by structural adaptations in microvascular network architecture are required for development, wound healing, and often invoked in disease conditions, including the top eight causes of death in the Unites States. Effective characterization of microvascular network architectures is not only limited by the available techniques to visualize microvessels but also reliant on the available quantitative metrics that accurately delineate between spatial patterns in altered networks. In this review, we survey models used for studying the microvasculature, methods to label and image microvessels, and the metrics and software packages used to quantify microvascular networks. These programs have provided researchers with invaluable tools, yet we estimate that they have collectively attained low adoption rates, possibly due to limitations with basic validation, segmentation performance, and nonstandard sets of quantification metrics. To address these existing constraints, we discuss opportunities to improve effectiveness, rigor, and reproducibility of microvascular network quantification to better serve the current and future needs of microvascular research.

Using high resolution X-ray computed tomography to create an image based model of a lymph node

Lymph nodes are an important part of the immune system. They filter the lymphatic fluid as it is transported from the tissues before being returned to the blood stream. The fluid flow through the nodes influences the behaviour of the immune cells that gather within the nodes and the structure of the node itself. Measuring the fluid flow in lymph nodes experimentally is challenging due to their small size and fragility. In this paper, we present high resolution X-ray computed tomography images of a murine lymph node. The impact of the resulting visualized structures on fluid transport are investigated using an image based model. The high contrast between different structures within the lymph node provided by phase contrast X-ray computed tomography reconstruction results in images that, when related to the permeability of the lymph node tissue, suggest an increased fluid velocity through the interstitial channels in the lymph node tissue. Fluid taking a direct path from the afferent to the efferent lymphatic vessel, through the centre of the node, moved faster than the fluid that flowed around the periphery of the lymph node. This is a possible mechanism for particles being moved into the cortex.

Investigation of microvascular morphological measures for skeletal muscle tissue oxygenation by image-based modelling in three dimensions

The supply of oxygen in sufficient quantity is vital for the correct functioning of all organs in the human body, especially for skeletal muscle during exercise. Traditionally, microvascular oxygen supply capability is assessed by the analysis of morphological measures on transverse cross-sections of muscle, e.g. capillary density or capillary-to-fibre ratio. In this work, we investigate the relationship between microvascular structure and muscle tissue oxygenation in mice. Phase contrast imaging was performed using synchrotron radiation computed tomography (SR CT) to visualize red blood cells (RBCs) within the microvasculature in mouse soleus muscle. Image-based mathematical modelling of the oxygen diffusion from the RBCs into the muscle tissue was subsequently performed, as well as a morphometric analysis of the microvasculature. The mean tissue oxygenation was then compared with the morphological measures of the microvasculature. RBC volume fraction and spacing (mean distance of any point in tissue to the closest RBC) emerged as the best predictors for muscle tissue oxygenation, followed by length density (summed RBC length over muscle volume). The two-dimensional measures of capillary density and capillary-to-fibre ratio ranked last. We, therefore, conclude that, in order to assess the states of health of muscle tissue, it is advisable to rely on three-dimensional morphological measures rather than on the traditional two-dimensional measures.

Simultaneous visualisation of calcified bone microstructure and intracortical vasculature using synchrotron X-ray phase contrast-enhanced tomography

3D imaging of the bone vasculature is of key importance in the understanding of skeletal disease. As blood vessels in bone are deeply encased in the calcified matrix, imaging techniques that are applicable to soft tissues are generally difficult or impossible to apply to the skeleton. While canals in cortical bone can readily be identified and characterised in X-ray computed tomographic data in 3D, the soft tissue comprising blood vessels that are putatively contained within the canal structures does not provide sufficient image contrast necessary for image segmentation. Here, we report an approach that allows for rapid, simultaneous visualisation of calcified bone tissue and the vasculature within the calcified bone matrix. Using synchrotron X-ray phase contrast-enhanced tomography we show exemplar data with intracortical capillaries uncovered at sub-micrometre level without the need for any staining or contrast agent. Using the tibiofibular junction of 15 week-old C57BL/6 mice post mortem, we show the bone cortical porosity simultaneously along with the soft tissue comprising the vasculature. Validation with histology confirms that we can resolve individual capillaries. This imaging approach could be easily applied to other skeletal sites and transgenic models, and could improve our understanding of the role the vasculature plays in bone disease.