Next Article in Journal
Exosomal microRNAs: Pleiotropic Impacts on Breast Cancer Metastasis and Their Clinical Perspectives
Next Article in Special Issue
Mechanisms of Venoarteriolar Reflex in Type 2 Diabetes with or without Peripheral Neuropathy
Previous Article in Journal
The Expression of Anti-Müllerian Hormone Type II Receptor (AMHRII) in Non-Gynecological Solid Tumors Offers Potential for Broad Therapeutic Intervention in Cancer
Previous Article in Special Issue
A Stress Syndrome Prototype Reflects Type 3 Diabetes and Ischemic Stroke Risk: The SABPA Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

3D Imaging and Quantitative Characterization of Mouse Capillary Coronary Network Architecture

Univ. Bordeaux, Inserm, UMR1034, Biology of Cardiovascular Diseases, F-33600 Pessac, France
*
Author to whom correspondence should be addressed.
Biology 2021, 10(4), 306; https://doi.org/10.3390/biology10040306
Submission received: 1 March 2021 / Revised: 1 April 2021 / Accepted: 4 April 2021 / Published: 7 April 2021
(This article belongs to the Special Issue Microcirculation in Health and Disease)

Abstract

:

Simple Summary

The cardiovascular system is composed of two physically and functionally linked parts, the heart and the vascular network, composed of the arteries, capillaries, and veins. In this system, the heart maintains the blood flow, including its own through the coronary vascular network, bringing oxygen and nutrients required for its activity. Microvascularization constitutes the distal part of this network and is composed of small vessels, i.e., arterioles and capillaries. Coronary microvascular diseases are associated with the development of cardiac dysfunction. Identification of the coronary microvascular network structure is hence of critical importance for the understanding of heart diseases. The aim of our study was to establish an accessible methodology for 3D imaging and 3D processing to quantitatively characterize the capillary coronary network architecture in mice, a species largely used as an animal model of cardiac failure. This study proposes a standardized methodology for 3D image processing that allows the quantification of the 3D capillary architecture. Based on an open source software, the whole process is easily accessible for biologists. It is hence useful for the study of cardiovascular diseases and the development of future therapies.

Abstract

Characterization of the cardiac capillary network structure is of critical importance to understand the normal coronary functional properties and coronary microvascular diseases. The aim of our study was to establish an accessible methodology for 3D imaging and 3D processing to quantitatively characterize the capillary coronary network architecture in mice. Experiments were done on C57BL/6J mice. 3D imaging was performed by light sheet microscopy and confocal microscopy on iDISCO+ optical cleared hearts after labelling of the capillary endothelium by lectin injection. 3D images were processed with the open source software ImageJ. Non-visual image segmentation was based of the frequency distribution of the voxel greyscale values, followed by skeletonization and distance mapping. Capillary networks in left and right ventricles and septum were characterized by the volume network density, the fractal dimension, the number of segments and nodes and their ratio, the total network length, and the average length, diameter, and tortuosity of the segments. Scale-dependent parameter values can be impacted by the resolution limit of the 3D imaging technique. The proposed standardized methodology for 3D image processing is easily accessible for a biologist in terms of investment and difficulty level, and allows the quantification of the 3D capillary architecture and its statistical comparison in different conditions.

1. Introduction

Oxygen and nutrient delivery to the heart is ensured by blood flow through the coronary vascular network, composed by arteries, capillaries, and veins. The structure of this network is a major determinant of heart metabolism and functional capacities. The microvascularization constitutes the distal part of this network and is composed of the arterioles and capillaries with a diameter lower than 100-150 µm. Exchanges between blood and cardiac tissue take place in the capillary network. Structural or functional alterations of the coronary microvascularization are referred as coronary microvascular diseases (CMDs) [1]. These CMDs are associated with the development of cardiac pathologies, including angina and heart failure [2].
Identification of the coronary microvascular network structure is hence of critical importance for the understanding of normal coronary functional properties and the understanding of several cardiac dysfunctions. Historically, most of the studies that intended to analyze human and animal coronary capillary networks have been done in 2D. In humans, the coronary microcirculation is hardly detectable by standard imaging techniques [1]. 2D angiography has been used to measure the large coronary vessel diameters, but not microvessels [3]. In mice, histological techniques have been used for 2D measurements of capillary diameter and density [4,5]. During the last few years, the use of 3D imaging techniques to study blood vessels has increased. For example, such techniques have been used to study and quantify the renal arterial and arteriolar system [6] and the brain vascular system [7,8]. However, to our knowledge, no studies have been published that have quantified the coronary capillary 3D architecture. The aim of our study was to establish an accessible methodology for 3D imaging and 3D processing to quantitatively characterize the capillary coronary network architecture in mice, a species largely used as an animal model of cardiac failure [9].
For 3D imaging, our strategy was to use 3D photonic imaging to visualize the coronary capillary system on intact mice hearts. A major problem is light scattering that limits light penetration in the tissue. Clearing methods that homogenize the refractive index have been developed in the past decade to overcome this difficulty [10]. Among the several techniques available for optical clearing, we chose the iDISCO+ method [11]. Its main advantages are easy access, short clearing time and high clearing capacity. Labelling of the capillary network was ensured by lectin coupled to a fluorophore, as lectin has been shown to be a specific marker of the capillary endothelium [12]. Cleared samples were imaged by light sheet microscopy (LSM) that allows large 3D visualization in intact samples [10]. Since confocal microscopy can also be used on cleared tissues, we also performed complementary analysis with this imaging technique.
For 3D image processing, we developed a standardized methodology to extract 3D quantitative data with the open source software ImageJ. Capillary network structure was characterized with nine parameters divided in two types: (1) global parameters related to the overall structure of the vascular system (vascular density, fractal dimension, number of segments, total length, number and percentage of node), and (2) architectural parameters related to the structure of each segment of the vascular system (length, diameter, and tortuosity). Comparisons between the left (LV) and right (RV) ventricles and the septum (S) were made, as no comparative information about mouse physiological coronary capillaries regarding the different parts of the heart was found in the literature. The only comparison found concerned the rat LV and feline LV and RV [13,14].

2. Materials and Methods

2.1. Animal Model

Experiments were performed on 8-week-old C57BL/6J mice (n = 7). Food and water were available ad libitum, with a 12 h dark/light cycle.

2.2. Arterial Pressure Measurements

Systolic (SAP), diastolic (DAP), and mean arterial pressure (MAP) were measured on vigil animals using the CODA® High Throughput System (Kent Scientific corporation, CT, USA). This non-invasive method is based on the analysis of the caudal blood flow by a sensor cuff after the complete occlusion of the blood flow and its progressive recovery controlled by an occlusion cuff placed at the basis of the tail. Each mouse was placed in an animal holder on a warming support in order to vasodilate the caudal vasculature. The tail external temperature was monitored and the recording began when it reached 32 °C. Measurements were realized on 15 consecutive cycles over 15 min, and the mean SAP, DAP, and MAP were calculated on these 15 cycles.

2.3. Mice Preparation

Each mouse received a retroorbital intravenous 100 µL injection of lectin coupled to a fluorophore (Lycopersicon esculentum Lectin DyLight649® 1 mg/mL, Vector Labs, Eurobio Scientific, Les Ulis, France) with an insulin syringe. Ten minutes after the injection, an intraperitoneal 100 µL injection of isosorbide dinitrate (Risordan® 10 mg/10mL, Medisol, Lyon, France) was made with 25G needle to dilate the vessels. The mouse was then euthanized by an intraperitoneal injection of 300 µL sodium pentobarbital (Exagon® à 400 mg/mL, Axience, Pantin, France) diluted in physiological saline solution. After death, a sternotomy was carried out to catheterize the left ventricle. A perfusion of physiological solution at 80 mm Hg pressure for 3 min was done to remove the blood from the vasculature. A second perfusion of 4% formalin (10% neutral buffered formalin, DiaPath, Martinengo, Italy) was done to fix the tissues. The heart was delicately removed and placed in paraformalin overnight at 4 °C. A negative control was done in similar conditions, but without lectin injection.

2.4. Optical Clearing

Optical clearing of the tissues was performed using the iDISCO+ method [11]. The technique consists of a methanol pre-treatment, followed by methanol and dichloromethane permeabilization and lipid removal, and RI matching in dibenzyl ether.
Each heart was pre-treated by successive immersions in 20%, 40%, 60%, and 80% methanol solutions for 1 h each at room temperature, then left overnight in a pure methanol solution (Methanol ≥ 99.5%, GPR RECTAPUR®, VWR Chemicals, Fontenay-sous-Bois, France) at room temperature. Permeabilization and lipid removal were then performed in a 2/3 dichloromethane (Dichloromethane, anhydro ≥ 99.8%, with 40–150 ppm amylene, Sigma-Aldrich, Lyon, France) and 1/3 methanol solution for 3 h at room temperature. The remaining methanol was removed by two 15 min baths in pure dichloromethane solutions at room temperature. Last, RI matching was done by leaving the sample in a dibenzyl ether solution (Benzyl ether 98%, Sigma-Aldrich) for a few hours at room temperature, then kept at 4 °C until imaging.

2.5. Shrinkage Measurement

The iDISCO+ technique has been shown to generate a shrinkage of the tissues [10]. On each sample, the shrinkage was estimated from two orthogonal 2D-images taken before and after clearing. Length (l), width (w), and thickness (t) were measured using ImageJ/Fiji software. These measurements were used to calculate an estimated heart volume, considered as an ellipsoid, thanks to the following equation:
C a r d i a c   V o l u m e = 4 3 × π × l 2 × w 2 × t 2
The shrinkage percentage was calculated for the length, the width, the thickness, and the ellipsoidal volume as the ratio of each parameter value after and before optical clearing on the parameter value before optical clearing.

2.6. Image Acquisition

2.6.1. Light Sheet Microscopy

Each heart was mounted in ethyl cinnamate solution in a specific device adapted to the light sheet microscope. The ultramicroscopy was done using the system from LaVision BioTec (Bielefeld, Germany) equipped with a 639 nm (70 mW) laser line, a sCMOS Andor camera, and a 0.5 NA 2× objective with a deeping lens, with 6.3 zoom. The tissue was illuminated laterally by three horizontal sheets. Exposition time was 200 ms. Emitted fluorescence was collected at 690 nm. At the end of the acquisition, the stack was automatically reconstructed in 16 bits format. Capillary network imaging was made by sampling three 1080 × 1280 × 300 µm3 parallelepipedic sections in the left ventricle (LV), the septum (S), and the right ventricle (RV) of the heart. The system spatial resolution was 1 µm (x, y) and 4 µm (z). Step size used was 2 µm. Voxel dimensions were 0.5 µm (x, y) and 2 µm (z).

2.6.2. Confocal Microscopy

3D imaging of cleared tissues is usually done either by light sheet microscopy or confocal microscopy [10]. To compare the potential variations due to the imaging technique, LV confocal microscopy 3D imaging was performed on one heart sample used for light sheet microscope imaging. A thin section of LV was cut and mounted in a specific device adapted to confocal microscopy for cleared samples. The tissue was illuminated by a laser at 639 nm and the emitted fluorescence at 690 nm was collected by a HCX PL APO CS 20.0 × 0.70 IMM UV objective. At the end of the acquisition, the stack was automatically reconstructed in 16-bits format. Capillary network imaging was made by sampling three 387 × 387 × 105 µm3 parallelepipedic sections. The system spatial resolution was 279 nm (x, y) and 1284 nm (z). The step size used was 630 nm. Voxel dimensions were 380 nm (x, y) and 630 nm (z).

2.7. Image Processing

Commands and functions for each step of the image processing are given in File S1.
The stacks were converted to 8 bits format (256 grey levels) and image processing was performed using the open-source ImageJ/Fiji software (ImageJ 2.1.0/1.53h/Java 1.8.0_66 (64 bit)) and several of its plugins.

2.7.1. Segmentations

The extraction of the two structures of interest, i.e., the whole cardiac tissue and the capillary network, was performed by a segmentation technique based on two greyscale-defined thresholds. For each sample, the threshold values were determined from the frequency distribution of the voxel greyscale values. This distribution showed the existence of three voxel populations corresponding, from 0 to 255 grey values, to (1) the image background, (2) the non-labelled cardiac tissue, and (3) the lectin-labelled capillary network, respectively. The first threshold, segmenting the cardiac tissue from the background, was set at the inflexion point of the distribution curve between the first and second population. Due to the overlap of the 2 voxel populations corresponding to the cardiac tissue and the lectin-labelled capillary network, no inflexion point was available for segmenting the capillary network, making difficult to identify an objective threshold value. However, the negative control showed that the non-labelled cardiac tissue is an homogenous population of voxels following a Gaussian distribution. Hence, on each sample, a Gaussian fit of the voxel frequency distribution was performed to determine the mean and standard deviation (SD) of the voxel population corresponding to the non-labelled cardiac tissue, and the threshold value was set at mean + 1.96 SD. Due to the mathematical properties of the normal distribution, such a threshold value ensured the exclusion of 97.5% of the voxels that do not correspond to the capillary network. Frequency distribution curves of the negative control and of one representative lectin-labelled sample are given in Figure 1a–d.

2.7.2. Filtering

After binarization of the segmented images, a 3D median filter set at 3 voxels for (x, y, z) was applied to suppress artefactual isolated voxels.

2.7.3. Skeletonization and Distance Mapping

After segmentation, a skeletonization of the segmented capillary network image was performed using the BoneJ plugin [15]. In the resulting skeleton, each segment of the network was transformed into a one voxel line corresponding to the center of the segment and identified by its own orthonormal coordinates (x, y, z). In parallel, distance mapping was performed on the segmented image, using the 3D Suites plugin [16]. Basically, distance mapping computes for each voxel of the segment axis the radius of the segment at that point. Skeleton and distance map images were fused before skeleton analysis. This provided information on the number of segments and the number of junctions between segments and, for each individual segment, the orthonormal coordinates of its extremities, its length, its mean radius, and its Euclidian distance, i.e., the linear distance between the two ends of the segment. Segments with a Euclidean distance equal to 0, are considered as artefactual loops as their starting and end points are the same, and segments with a radius inferior to 1 voxel, reaching the resolution limit, were removed from the data. The mean percentage of these artefactual loops was 4.62% of the total number of segments. Representative images of LV capillary network of image processing are shown in Figure 1e,f.

2.8. Data Analysis

The data obtained from skeletonization and distance mapping were used to determine 6 “global” parameters, i.e., considering the capillary network as a whole (vascular density, fractal dimension, number of segments, number and percentage of nodes, and total length), and 3 “topological” parameters, i.e., calculated from each segment of the network (length, diameter and tortuosity).

2.8.1. Cardiac Volumes and Vascular Density

The total cardiac volume (non-labelled cardiac tissue + lectin-labelled capillary network volumes) and the capillary network volume were obtained from their respective segmented binary images. For each object, the volume, in mm3, was calculated by multiplying the number of white voxels, corresponding to the volume object in voxels, by the volume of one voxel in mm3.
The volume vascular density, expressed as %, was calculated according to the following equation:
V a s c u l a r   d e n s i t y = C a p i l l a r y   v o l u m e T o t a l   v o l u m e

2.8.2. Fractal Dimension

The fractal dimension measures the ability of a self-similar structure to fill the three-dimensional space it occupies. The fractal dimension of the capillary network was calculated on the skeletonized image using the Minkowski–Bouligan method [17], also called Box Counting method, with BoneJ plugin.

2.8.3. Normalized Number of Segments

The number of segments obtained from the skeletonization was normalized to the total cardiac tissue volume (in mm3).

2.8.4. Normalized Total Capillary Length

The length of each segment of the skeleton was obtained in µm from the skeletonization, and the normalized total length of the network, expressed in m.mm-3, was calculated by summing up all segment lengths converted in meter, normalized to the total cardiac tissue volume in mm3.

2.8.5. Normalized Number and Percentage of Nodes

The normalized number of nodes was defined as the number of junctions formed by the capillary network normalized to the total cardiac tissue volume (mm3). The percentage of nodes was calculated by dividing the number of nodes by the number of segments. The node/segments ratio depends on the number of segments connected at each node, and can hence be used as an index of the connectivity of the network.

2.8.6. Segment Diameter

The segment mean radii were obtained in µm from the skeletonization in µm, and the segment mean diameters were calculated by multiplying the mean radii by 2.

2.8.7. Tortuosity

The tortuosity, used as an index representing the twisting of each capillary segment, was calculated using the following equation:
T o r t u o s i t y = S e g m e n t   l e n g h t E u c l i d i a n   d i s t a n c e

2.9. Statistical Analysis

Statistical analyses were done using GraphPad Prism® software (Prism 9.0.1 version), (San Diego, CA, USA).

2.9.1. Global Parameters

Vascular density, fractal dimension, number of segments, number and percentage of nodes, and total length are expressed as mean ± standard deviation (SD). Statistical comparisons between LV, S and RV were done with Kruskal–Wallis non-parametric tests followed, if relevant, by post-hoc Dunn’s test. Results were considered as statistically significant for p < 0.05.

2.9.2. Topological Parameters

For each topological parameter, i.e., segment length, diameter and tortuosity, frequency classes were defined and the relative frequency distribution was calculated. The frequency distributions were adjusted by non-linear regression. For each parameter, the best fitting curve of the frequency distribution was determined by statistical comparison of Gaussian versus one-phase exponential decay equations with F test using GraphPad Prism. Statistical comparisons between LV, S and RV were performed with F test. Results were considered as statistically significant for p < 0.05.

3. Results

3.1. Cardiovascular Parameters

Arterial pressures (SAP, DAP, MAP) and cardiac rhythm were measured on six mice. Mean ± SD are given in Table 1.

3.2. Shrinkage

The shrinkage percentage in length, width, thickness and ellipsoidal volume was calculated on seven hearts, and are given in Table 2. Shrinkage values were very similar in length, width, and thickness, and their statistical comparison by Kruskal–Wallis test showed no significant difference, indicating that the shrinkage was isotropic.

3.3. Global Parameters

3.3.1. Vascular Density

Mean vascular density was calculated for LV, S, and RV on 6 mice. Results are presented in Figure 2a and mean ± SD are given in Table 3. The vascular density was significantly higher in the LV compared to S, but not between LV and RV, nor between RV and S.

3.3.2. Fractal Dimension

The fractal dimension was calculated on the skeletonized images of the LV, the S, and the RV of six hearts. According to the Minkowski–Bouligan method, the box counting results were fitted by a simple linear regression (R² = 0.99), which allowed us to calculate the fractal dimension as the curve slope value. Results are presented in Figure 2b and mean ± SD are given in Table 3. Statistical comparison showed no significant difference between LV, S, and RV.

3.3.3. Normalized Number of Segments

The number of segments normalized to the volume of cardiac tissue was calculated for the LV, the S, and the RV of six hearts. Results are presented in Figure 2c and mean ± SD are given in Table 3. Statistical comparison showed no significant difference between LV, S, and RV.

3.3.4. Total Capillary Length

The total capillary length was calculated for the LV, the S, and the RV of six hearts. Results are presented in Figure 2d and mean ± SD are given in Table 3. Statistical comparison showed no significant difference between LV, S, and RV.

3.3.5. Number and Percentage of Nodes

The number and percentage of nodes were calculated for the LV, the S, and the RV of six hearts. Results are presented in Figure 2e,f, and mean ± SD are given in Table 3. For the normalized number of segments, no significant difference was found for the normalized number of nodes between LV, S and RV. However, the percentage on nodes was significantly higher in RV compared to S, but not between LV and RV, nor between LV and S.

3.4. Topological Parameters

3.4.1. Segment Length

Segment length distributions in LV, S and RV (n = 6) were fitted by the following one-phase exponential decay equation:
F = F 0 F × e K   ×   L + F
where L is the segment length (in µm), F is the relative frequency (in %), F 0 is F for L = 0, F   is F for L = ∞, and K is the rate constant (in µm−1). From the rate constant K the length constant λ was defined as 1/K, in µm, and used as an index of length distribution. The shorter λ is, the higher the proportion of sort segments is. Curves are presented in Figure 3a and λ values are given in Table 4. The median values are 12.34, 12.88, and 12.01 for LV, RV, and S, respectively. Statistical comparisons showed no significant differences for λ between RV, S, and LV.

3.4.2. Diameter

Diameter length distributions in LV, S, and RV (n = 6) were fitted by the following Gaussian equation:
F = A × e 0.5 × ( d µ / σ ) 2
where d is the segment diameter (in µm), F is the relative frequency (in %), A is the maximal frequency, µ is the diameter mean value (in µm), and σ is the standard deviation (in µm). Curves are given in Figure 3b and µ and σ are given in Table 4. Statistical comparisons showed that µ was significantly lower in S versus RV and LV, and lower is LV versus RV.

3.4.3. Tortuosity

Tortuosity distributions in LV, S and RV (n = 6) were fitted by the following one-phase exponential decay equation:
F = F 0 F × e K   ×   t + F
where t is the segment tortuosity (unitless), F is the relative frequency (in %), F 0 is F for t = 0, F   is F for t = ∞, and K is the rate constant (in µm1). From the rate constant K the tortuosity constant τ was defined as 1/K, and used as an index of tortuosity distribution. Curves are presented in Figure 3c and τ values are given in Table 4. τ was below 1, indicating that the capillaries were almost straight. The median values are 1.22, 1.33, and 1.24 for LV, RV, and S, respectively. Statistical comparisons showed non significant differences between RV, S, and LV.

3.5. Confocal Microscopy

On one heart, both global parameters and architectural parameters of the LV capillary network were calculated from images obtained by confocal microscopy, in order to be compared with the values obtained in the same LV heart from images obtained by light sheet microscopy. A 3D image of the capillary network and the architectural parameter curves obtained from confocal microscopy are given in Figure 4. Data obtained from both imaging methods are given in Table 5.

4. Discussion

Our study has established a standardized method of 3D image processing applicable to 3D photonic imaging of the coronary capillary network in mice. This method generates both global and segment-level data allowing 3D quantitative statistical analysis and comparison between samples. Our results provided quantitative characterization of LV, RV, and S capillary networks in normal mice.

4.1. 3D Imaging and Image Processing

The iDISCO+ optical clearing method used in this study is known, according to its developers, to induce an 11% volume shrinkage with minor deformation of the brain [11]. Our results showed an estimated heart volume shrinkage around 19%, without significant difference in each dimension, suggesting that, as in the brain, the iDISCO+ method induced a limited and isotropic heart reduction.
One key point for the quantitative analysis of structures such as the coronary network extracted from 3D images is unbiased image segmentation. Due to the poor grayscale resolution power of the eye, a non-visual method was applied for the determination of the threshold values, based on the statistical analysis of the voxel distribution. Because of the overlapping grayscale values of the different objects of interest, i.e., the cardiac tissue and the capillary network, objective identification of the voxel population corresponding to the vascular network was impossible. The segmentation method was hence based on the identification of the cardiac tissue. Indeed, the analysis of non-labelled samples showed that the cardiac tissue corresponded to a homogenous voxel population following a Gaussian distribution. Fluorescence of the cardiac tissue at 639nm has already been reported [18]. For a Gaussian distribution, the mean±1.96SD interval corresponds to 95% of the population. Hence, a threshold value set at mean+1.96SD, as the voxel population corresponding to non-labelled tissue has a normal distribution, ensured the exclusion of 97.5% of non-labelled (non-capillary) tissue. Though the voxel distribution may vary from sample to sample, and despite the overlap of voxel populations, this method ensured that, for each sample, the same proportion of non-capillary cardiac tissue is excluded from the capillary network, allowing unbiased statistical comparison.

4.2. Global and Archhitectural Parameters

Our LSM results showed a number of capillaries close to 350,000 capillaries/mm3 of cardiac tissue, with an average length around 15 µm, representing nearly 30% of the cardiac volume, and a total capillary length around 10m/mm3 of cardiac tissue. To our knowledge, there are no previously published 3D data on the coronary capillary volume and number density. The coronary capillary number density has been calculated in 2D with histological images, and was estimated at 2500 capillaries/mm² [5,19]. Considering an average capillary length of 15 µm, this corresponds to 170,000 capillaries/mm3, a value lower but in the same range of our own values. The total capillary length/mm3 of cardiac tissue has already been calculated on rats by 2D measurements and mathematical 3D reconstruction at 3.5 m/mm3 [20]. It has also been calculated on mouse cleared brain with values near to 1m/mm3 of health cerebral tissue [7]. More recently, total vessel length has been calculated on mouse heart by 3D imaging and found to be 3m/mm3 [18]. This value is consistent with our own, though a bit lower. The difference may be explained by differences in image resolution. Indeed, in these studies, voxel size was higher than in our study, and hence may have failed to identify small segments.
The percentage of nodes on the number of segments is a marker of the connectivity pattern of the network since it depends on the relative proportion of the different degrees of furcation (bifurcation, trifurcations, quadfurcation). Removal of artefactual loops may generate artefactual one-to-one junctions between two segments, which may overestimate the number of nodes and segments. However, these possible biases are negligible since the percentage of artefactual loops involved in any type of multifurcation is very low (around 4%).
The fractal dimension is a useful parameter to characterize the complexity of the 3D architecture of a network [21]. Our results showed a fractal dimension close to 2.4, with a relative standard deviation less than 5%, hence with less interindividual variation than the vascular density, whose relative standard deviation was around 35%. Some studies have calculated the fractal dimension of the arterial network in kidney, which was found to be close to 2 [6], and, to our knowledge, it has not yet been calculated for the coronary capillary network. As far as the comparison with the kidney arterial tree makes sense, it indicates that the capillary network is more dense and complex than the arterial one.
Considering the capillary diameters, independently of the species and organ, the value is usually found between 3 and 10 µm, and is estimated around 4 µm in mouse heart [4,14]. These estimations mainly result from 2D studies. Our 3D study aligns with the 2D results observed in the literature as the mean capillary diameter is set around 4 µm. Our results also acknowledge that our methodology based on a lectin injection allows the characterization of the capillary network, as more than 90% of the studied population is composed of vessels with a diameter set between 2 and 10 µm.
The tortuosity, calculated by different methods, has already been studied on large vessels and on microvessels as a pathological marker [22]. The ratio obtained in our study was characterized by a tortuosity constant inferior to 1, and median values were close to 1, the theoretical value of a straight line, indicating that the capillaries were not tortuous.
Finally, comparison between LV, RV, and S showed that the capillary network is broadly similar in these different parts of the heart, with a lower volume and node density in the S.

4.3. Light Sheet Microscopy Versus Confocal Microscopy

Data obtained from CM compared to LSM showed that for scale-free parameters, such as the fractal dimension, the percentage of nodes and the tortuosity constant, the difference between LSM and CM in the same heart were of similar amplitude than the SD of the parameters obtained in LSM. By contrast, the number of segments and nodes, and the total segment length, were much greater in CM, whereas the average length and diameter were lower. These results suggest that CM imaging was able to identify small segments that were not detected by LSM imaging. The difference might be due to the resolution limit of LSM, close to the capillary diameter.

5. Conclusions

In conclusion, our study proposes a standardized methodology for 3D image processing of a coronary capillary network that allows the quantification of the 3D capillary architecture and its statistical comparison in different conditions. Thus, it might be used to study pathological coronary capillary networks, and compare them with healthy ones. It is applicable to different imaging techniques, e.g., LSM and CM, and, being based on an open source software, the whole process is easily accessible for a biologist in terms of investment and difficulty level. It also provides a set of parameter values for normal C57BL/6J mice. This study should be completed and confirmed by further studies, taking into account that the scale-dependent parameter values can be impacted by the resolution limit of the 3D imaging technique.

Supplementary Materials

The following is available online at https://www.mdpi.com/article/10.3390/biology10040306/s1, File S1: Image processing commands and functions.

Author Contributions

Conceptualization, E.R. and N.N.; methodology, E.R. and N.N.; formal analysis, E.R. and N.N.; investigation, E.R. and N.N.; resources, E.R. and N.N.; data curation, E.R. and N.N.; writing—original draft preparation, N.N.; writing—review and editing, E.R. and N.N.; visualization, E.R. and N.N.; supervision, E.R.; project administration, E.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the declaration of Helsinki, and approved by the Ethics Committee of Bordeaux (Agreement n°14896 April 2019).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The ultramicroscopy was done in the Bordeaux Imaging Center a service unit of the CNRS-INSERM and Bordeaux University, member of the national infrastructure France BioImaging supported by the French National Research Agency (ANR-10-INBS-04). The help of Jérémie Teillon is acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vancheri, F.; Longo, G.; Vancheri, S.; Henein, M. Coronary Microvascular Dysfunction. JCM 2020, 9, 2880. [Google Scholar] [CrossRef] [PubMed]
  2. Taqueti, V.R.; Di Carli, M.F. Coronary Microvascular Disease Pathogenic Mechanisms and Therapeutic Options. J. Am. Coll. Cardiol. 2018, 72, 2625–2641. [Google Scholar] [CrossRef] [PubMed]
  3. de Feyter, P.J.; Serruys, P.W.; Davies, M.J.; Richardson, P.; Lubsen, J.; Oliver, M.F. Quantitative Coronary Angiography to Measure Progression and Regression of Coronary Atherosclerosis. Value, Limitations, and Implications for Clinical Trials. Circulation 1991, 84, 412–423. [Google Scholar] [CrossRef] [Green Version]
  4. Müller, B.; Lang, S.; Dominietto, M.; Rudin, M.; Schulz, G.; Deyhle, H.; Germann, M.; Pfeiffer, F.; David, C.; Weitkamp, T. High-Resolution Tomographic Imaging of Microvessels. Proc. SPIE 2008, 7078, 70780B. [Google Scholar] [CrossRef] [Green Version]
  5. Wu, Y.; Ip, J.E.; Huang, J.; Zhang, L.; Matsushita, K.; Liew, C.-C.; Pratt, R.E.; Dzau, V.J. Essential Role of ICAM-1/CD18 in Mediating EPC Recruitment, Angiogenesis, and Repair to the Infarcted Myocardium. Circ. Res. 2006, 99, 315–322. [Google Scholar] [CrossRef]
  6. Markovič, R.; Peltan, J.; Gosak, M.; Horvat, D.; Žalik, B.; Seguy, B.; Chauvel, R.; Malandain, G.; Couffinhal, T.; Duplàa, C.; et al. Planar Cell Polarity Genes Frizzled4 and Frizzled6 Exert Patterning Influence on Arterial Vessel Morphogenesis. PLoS ONE 2017, 12, e0171033. [Google Scholar] [CrossRef]
  7. Lugo-Hernandez, E.; Squire, A.; Hagemann, N.; Brenzel, A.; Sardari, M.; Schlechter, J.; Sanchez-Mendoza, E.H.; Gunzer, M.; Faissner, A.; Hermann, D.M. 3D Visualization and Quantification of Microvessels in the Whole Ischemic Mouse Brain Using Solvent-Based Clearing and Light Sheet Microscopy. J. Cereb. Blood Flow Metab. 2017, 37, 3355–3367. [Google Scholar] [CrossRef]
  8. Smith, A.F.; Doyeux, V.; Berg, M.; Peyrounette, M.; Haft-Javaherian, M.; Larue, A.-E.; Slater, J.H.; Lauwers, F.; Blinder, P.; Tsai, P.; et al. Brain Capillary Networks Across Species: A Few Simple Organizational Requirements Are Sufficient to Reproduce Both Structure and Function. Front. Physiol. 2019, 10, 233. [Google Scholar] [CrossRef] [Green Version]
  9. Valero-Muñoz, M.; Backman, W.; Sam, F. Murine Models of Heart Failure With Preserved Ejection Fraction. JACC Basic Transl. Sci. 2017, 2, 770–789. [Google Scholar] [CrossRef] [PubMed]
  10. Costa, E.C.; Silva, D.N.; Moreira, A.F.; Correia, I.J. Optical Clearing Methods: An Overview of the Techniques Used for the Imaging of 3D Spheroids. Biotechnol. Bioeng. 2019, 116, 2742–2763. [Google Scholar] [CrossRef]
  11. Renier, N.; Adams, E.L.; Kirst, C.; Wu, Z.; Azevedo, R.; Kohl, J.; Autry, A.E.; Kadiri, L.; Umadevi Venkataraju, K.; Zhou, Y.; et al. Mapping of Brain Activity by Automated Volume Analysis of Immediate Early Genes. Cell 2016, 165, 1789–1802. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Robertson, R.T.; Levine, S.T.; Haynes, S.M.; Gutierrez, P.; Baratta, J.L.; Tan, Z.; Longmuir, K.J. Use of Labeled Tomato Lectin for Imaging Vasculature Structures. Histochem. Cell Biol. 2015, 143, 225–234. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Breisch, E.A.; Houser, S.R.; Carey, R.A.; Spann, J.F.; Bove, A.A. Myocardial Blood Flow and Capillary Density in Chronic Pressure Overload of the Feline Left Ventricle. Cardiovasc. Res. 1980, 14, 469–475. [Google Scholar] [CrossRef] [PubMed]
  14. Potter, R.F.; Groom, A.C. Capillary Diameter and Geometry in Cardiac and Skeletal Muscle Studied by Means of Corrosion Casts. Microvasc. Res. 1983, 25, 68–84. [Google Scholar] [CrossRef]
  15. Doube, M.; Kłosowski, M.M.; Arganda-Carreras, I.; Cordelières, F.P.; Dougherty, R.P.; Jackson, J.S.; Schmid, B.; Hutchinson, J.R.; Shefelbine, S.J. BoneJ: Free and Extensible Bone Image Analysis in ImageJ. Bone 2010, 47, 1076–1079. [Google Scholar] [CrossRef] [Green Version]
  16. Ollion, J.; Cochennec, J.; Loll, F.; Escudé, C.; Boudier, T. TANGO: A Generic Tool for High-Throughput 3D Image Analysis for Studying Nuclear Organization. Bioinformatics 2013, 29, 1840–1841. [Google Scholar] [CrossRef] [PubMed]
  17. Luppe, M. Fractal Dimension Based on Minkowski-Bouligand Method Using Exponential Dilations. Electron. Lett. 2015, 51, 475–477. [Google Scholar] [CrossRef]
  18. Merz, S.F.; Korste, S.; Bornemann, L.; Michel, L.; Stock, P.; Squire, A.; Soun, C.; Engel, D.R.; Detzer, J.; Lörchner, H.; et al. Contemporaneous 3D Characterization of Acute and Chronic Myocardial I/R Injury and Response. Nat. Commun. 2019, 10, 2312. [Google Scholar] [CrossRef]
  19. Schmidt-Nielsen, K.; Pennycuik, P. Capillary Density in Mammals in Relation to Body Size and Oxygen Consumption. Am. J. Physiol. Leg. Content 1961, 200, 746–750. [Google Scholar] [CrossRef]
  20. Amann, K.; Wiest, G.; Zimmer, G.; Gretz, N.; Ritz, E.; Mall, G. Reduced Capillary Density in the Myocardium of Uremic Rats—A Stereological Study. Kidney Int. 1992, 42, 1079–1085. [Google Scholar] [CrossRef] [Green Version]
  21. Reishofer, G.; Koschutnig, K.; Enzinger, C.; Ebner, F.; Ahammer, H. Fractal Dimension and Vessel Complexity in Patients with Cerebral Arteriovenous Malformations. PLoS ONE 2012, 7, e41148. [Google Scholar] [CrossRef] [PubMed]
  22. Han, H.-C. Twisted Blood Vessels: Symptoms, Etiology and Biomechanical Mechanisms. J. Vasc. Res. 2012, 49, 185–197. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Image processing. (a,b) Representative frequency distribution of pixel population of non-labelled (negative control) (a) and lectin-labelled (b) left ventricle after 3D imaging by light sheet microscopy. 1: background. 2: non-labelled cardiac tissue. 3: lectin-labelled capillary network. (c), Gaussian fit of the subset of pixels corresponding to the cardiac tissue from curve (a). (d), Gaussian fit of the subset of pixels corresponding to the cardiac tissue from curve (b). (e), Representative segmented image of left ventricle capillary network. The image was obtained by overlaying the original capillaries image and the binarized image of the capillaries obtain after segmentation. (f), Skeletonization and image mapping of the capillary network of (e). Average radii of each segment are encoded in false colors from blue to red. The image was obtained by overlaying the skeleton and the distance map images of the capillary network.
Figure 1. Image processing. (a,b) Representative frequency distribution of pixel population of non-labelled (negative control) (a) and lectin-labelled (b) left ventricle after 3D imaging by light sheet microscopy. 1: background. 2: non-labelled cardiac tissue. 3: lectin-labelled capillary network. (c), Gaussian fit of the subset of pixels corresponding to the cardiac tissue from curve (a). (d), Gaussian fit of the subset of pixels corresponding to the cardiac tissue from curve (b). (e), Representative segmented image of left ventricle capillary network. The image was obtained by overlaying the original capillaries image and the binarized image of the capillaries obtain after segmentation. (f), Skeletonization and image mapping of the capillary network of (e). Average radii of each segment are encoded in false colors from blue to red. The image was obtained by overlaying the skeleton and the distance map images of the capillary network.
Biology 10 00306 g001
Figure 2. Global parameters. Data obtained from the left ventricle (blue), septum (red) and right ventricle (green) of six hearts. Columns are means and error bars are standard deviations. Black dots are individual values. Data were compared by Kruskal–Wallis test with post-hoc Dunn’s test. ns = non significant. * p < 0.05. (a). Volume capillary density, in %. (b). Fractal dimension. (c). Number of capillary segments per mm3 of cardiac tissue. (d). Total length of capillary segments (m) per mm3 of cardiac tissue. (e). Number of nodes per mm3 of cardiac tissue. (f). Percentage of nodes on the number of segments.
Figure 2. Global parameters. Data obtained from the left ventricle (blue), septum (red) and right ventricle (green) of six hearts. Columns are means and error bars are standard deviations. Black dots are individual values. Data were compared by Kruskal–Wallis test with post-hoc Dunn’s test. ns = non significant. * p < 0.05. (a). Volume capillary density, in %. (b). Fractal dimension. (c). Number of capillary segments per mm3 of cardiac tissue. (d). Total length of capillary segments (m) per mm3 of cardiac tissue. (e). Number of nodes per mm3 of cardiac tissue. (f). Percentage of nodes on the number of segments.
Biology 10 00306 g002
Figure 3. Architectural parameters. Data obtained from the left ventricle (blue), septum (red), and right ventricle (green) of six hearts. (a). Relative frequency distribution of segment length. Data were fitted by one-phase exponential decay. (b). Relative frequency distribution of segment diameter. Data were fitted by Gaussian equation. (c). Relative frequency distribution of segment tortuosity. Data were fitted by one-phase exponential decay. Curves between LV, S, and RV were compared by F tests.
Figure 3. Architectural parameters. Data obtained from the left ventricle (blue), septum (red), and right ventricle (green) of six hearts. (a). Relative frequency distribution of segment length. Data were fitted by one-phase exponential decay. (b). Relative frequency distribution of segment diameter. Data were fitted by Gaussian equation. (c). Relative frequency distribution of segment tortuosity. Data were fitted by one-phase exponential decay. Curves between LV, S, and RV were compared by F tests.
Biology 10 00306 g003
Figure 4. Confocal microscopy. (a). Segmented image of left ventricle capillary network. (b). Relative frequency distribution of segment length. Data were fitted by one-phase exponential decay. (c). Relative frequency distribution of segment diameter. Data were fitted by Gaussian equation. (d). Relative frequency distribution of segment tortuosity. Data were fitted by one-phase exponential decay.
Figure 4. Confocal microscopy. (a). Segmented image of left ventricle capillary network. (b). Relative frequency distribution of segment length. Data were fitted by one-phase exponential decay. (c). Relative frequency distribution of segment diameter. Data were fitted by Gaussian equation. (d). Relative frequency distribution of segment tortuosity. Data were fitted by one-phase exponential decay.
Biology 10 00306 g004
Table 1. Cardiovascular parameters. Systolic (ASP), diastolic (DAP), and mean (MAP) vascular pressure in mmHg and cardiac rhythm in beat per minute (BPM) were measured on six mice.
Table 1. Cardiovascular parameters. Systolic (ASP), diastolic (DAP), and mean (MAP) vascular pressure in mmHg and cardiac rhythm in beat per minute (BPM) were measured on six mice.
SAP (Mean ± SD)DAP (Mean ± SD)MAP (Mean ± SD)BPM (Mean ± SD)
123.1 ± 15 mm Hg89.7 ± 6.9 mm Hg100.5 ± 9.5 mm Hg303.6 ± 31
Table 2. Shrinkage. Length (L), width (W), thickness (T), and ellipsoid volume (V) reduction (in %) after optical clearing were measured on seven mouse hearts.
Table 2. Shrinkage. Length (L), width (W), thickness (T), and ellipsoid volume (V) reduction (in %) after optical clearing were measured on seven mouse hearts.
L (Mean ± SD)W (Mean ± SD)T (Mean ± SD)Volume (Mean ± SD)
6.70 ± 4.2%10.7 ± 4.2%2.91 ± 8.2%19.18 ± 7.7%
Table 3. Global parameters. Values were calculated from left ventricle (LV), septum (S), and right ventricle (RV) of six mouse hearts.
Table 3. Global parameters. Values were calculated from left ventricle (LV), septum (S), and right ventricle (RV) of six mouse hearts.
ParameterLVSRV
Vascular density34.4 ± 11%18.9 ± 4.7%27.8 ± 11%
Fractal dimension2.46 ± 0.052.32 ± 0.112.33 ± 0.07
Number of segments/mm3 of cardiac tissue615,784 ± 220,000399,922 ± 230,000387,457 ± 200,000
Total length/mm3 of cardiac tissue13.01 ± 6.1 m11.5 ± 9.2 m7.11 ± 2.7 m
Number of nodes/mm3 of cardiac tissue289,878 ± 96,000169,886 ± 92,000190,392 ± 100,000
Number of nodes/number of segments47.56 ± 3.0%42.9 ± 2.9%48.72 ± 2.9%
Table 4. Architectural parameters. Values were obtained by non-linear regression of frequency distribution of data from left ventricle (LV), septum (S), and right ventricle (RV) of six mouse hearts (see Figure 3). λ: length constant of one-phase exponential decay. μ and σ: mean and standard deviation of Gaussian distribution. τ: tortuosity constant of one-phase exponential decay.
Table 4. Architectural parameters. Values were obtained by non-linear regression of frequency distribution of data from left ventricle (LV), septum (S), and right ventricle (RV) of six mouse hearts (see Figure 3). λ: length constant of one-phase exponential decay. μ and σ: mean and standard deviation of Gaussian distribution. τ: tortuosity constant of one-phase exponential decay.
ParameterLVSRV
Length
λ (mean ± SEM)16.9 ± 0.6 µm15.6 ± 0.6 µm17.9 ± 4.4 µm
Diameter
μ (mean ± SEM)4.81 ± 0.24 µm3.78 ± 0.51 µm5.12 ± 0.25 µm
σ (mean ± SEM)2.52 ± 0.27 µm2.75 ± 0.44 µm2.17 ± 0.26 µm
Tortuosity
τ (mean ± SEM)0.32 ± 0.020.35 ± 0.020.35 ± 0.02
Table 5. Light sheet microscopy versus confocal microscopy. Global and architectural parameters obtained on light sheet microscopy (LSM, voxel size 0.5 × 0.5 × 2µm) and confocal microscopy (CM, voxel size 380 × 380 × 630 nm). The difference between LSM and CM (LSM-CM) values was calculated.
Table 5. Light sheet microscopy versus confocal microscopy. Global and architectural parameters obtained on light sheet microscopy (LSM, voxel size 0.5 × 0.5 × 2µm) and confocal microscopy (CM, voxel size 380 × 380 × 630 nm). The difference between LSM and CM (LSM-CM) values was calculated.
ParametersLSMCM(LSM-CM)
Vascular density43.3%31.2%12.1%
Fractal dimension2.492.530.04
Number of segments/mm3 of cardiac tissue794,0324,297,3303,503,298
Number of nodes/mm3 of cardiac tissue394,2961,992,2011,597,905
Number of nodes/number of segments49.7%46.4%3.3%
Total length/mm3 of cardiac tissue17.7 m27.1 m9.4 m
Length
λ (mean ± SEM)28.9 ± 9.8 µm6.89 ± 1.0 µm22.01 µm
Diameter
μ (mean ± SEM)4.61 ± 0.68 µm3.05 ± 0.29 µm1.56 µm
σ (mean ± SEM)2.25 ± 0.61 µm2.07 ± 0.26 µm0.18 µm
Tortuosity
τ (mean ± SEM)0.31 ± 0.070.27 ± 0.020.04
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nicolas, N.; Roux, E. 3D Imaging and Quantitative Characterization of Mouse Capillary Coronary Network Architecture. Biology 2021, 10, 306. https://doi.org/10.3390/biology10040306

AMA Style

Nicolas N, Roux E. 3D Imaging and Quantitative Characterization of Mouse Capillary Coronary Network Architecture. Biology. 2021; 10(4):306. https://doi.org/10.3390/biology10040306

Chicago/Turabian Style

Nicolas, Nabil, and Etienne Roux. 2021. "3D Imaging and Quantitative Characterization of Mouse Capillary Coronary Network Architecture" Biology 10, no. 4: 306. https://doi.org/10.3390/biology10040306

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop