Computer-assisted measurement of vessel shape from 3T magnetic resonance angiography of mouse brain

Blood vessel morphology (vessel radius, branching pattern, and tortuosity) is altered by a multitude of diseases. Although murine models of human pathology are important to the investigation of many diseases, there are few publications that address quantitative measurements of murine vascular morphology. This report outlines methods of imaging mice in vivo using magnetic resonance angiograms obtained on a clinical 3T unit, of defining mouse vasculature from these images, and of quantifying measures of vessel shape. We provide examples of both healthy and diseased vasculature and illustrate how the approach can be used to assess pathology both visually and quantitatively. The method is amenable to the assessment of many diseases in both human beings and mice.     

Microfocal X-ray CT imaging and pulmonary arterial distensibility in excised rat lungs

The objective of this study was to develop an X-ray computed tomographic method for measuring pulmonary arterial dimensions and locations within the intact rat lung. Lungs were removed from rats and their pulmonary arterial trees were filled with perfluorooctyl bromide to enhance X-ray absorbance. The lungs were rotated within the cone of the X-ray beam projected from a microfocal X-ray source onto an image intensifier, and 360 images were obtained at 1 degrees increments. The three-dimensional image volumes were reconstructed with isotropic resolution using a cone beam reconstruction algorithm. The vessel diameters were obtained by fitting a functional form to the image of the vessel circular cross section. The functional form was chosen to take into account the point spread function of the image acquisition and reconstruction system. The diameter measurements obtained over a range of vascular pressures were used to characterize the distensibility of the rat pulmonary arteries. The distensibility coefficient alpha [defined by D(P) = D(0)(1 + alphaP), where D(P) is the diameter at intravascular pressure (P)] was approximately 2.8% mmHg and independent of vessel diameter in the diameter range (about 100 to 2,000 mm) studied.     

Numerical evaluation reveals the effect of branching morphology on vessel transport properties during angiogenesis

Blood flow governs transport of oxygen and nutrients into tissues. Hypoxic tissues secrete VEGFs to promote angiogenesis during development and in tissue homeostasis. In contrast, tumors enhance pathologic angiogenesis during growth and metastasis, suggesting suppression of tumor angiogenesis could limit tumor growth. In line with these observations, various factors have been identified to control vessel formation in the last decades. However, their impacts on the vascular transport properties of oxygen remain elusive. Here, we take a computational approach to examine the effects of vascular branching on blood flow in the growing vasculature. First of all, we reconstruct a 3D vascular model from the 2D confocal images of the growing vasculature at postnatal day 5 (P5) mouse retina, then simulate blood flow in the vasculatures, which are obtained from the gene targeting mouse models causing hypo- or hyper-branching vascular formation. Interestingly, hyper-branching morphology attenuates effective blood flow at the angiogenic front, likely promoting tissue hypoxia. In contrast, vascular hypo-branching enhances blood supply at the angiogenic front of the growing vasculature. Oxygen supply by newly formed blood vessels improves local hypoxia and decreases VEGF expression at the angiogenic front during angiogenesis. Consistent with the simulation results indicating improved blood flow in the hypo-branching vasculature, VEGF expression around the angiogenic front is reduced in those mouse retinas. Conversely, VEGF expression is enhanced in the angiogenic front of hyper-branching vasculature. Our results indicate the importance of detailed flow analysis in evaluating the vascular transport properties of branching morphology of the blood vessels.     

Nanoparticle transport and delivery in a heterogeneous pulmonary vasculature

Quantitative understanding of nanoparticles delivery in a complex vascular networks is very challenging because it involves interplay of transport, hydrodynamic force, and multivalent interactions across different scales. Heterogeneous pulmonary network includes up to 16 generations of vessels in its arterial tree. Modeling the complete pulmonary vascular system in 3D is computationally unrealistic. To save computational cost, a model reconstructed from MRI scanned images is cut into an arbitrary pathway consisting of the upper 4-generations. The remaining generations are represented by an artificially rebuilt pathway. Physiological data such as branch information and connectivity matrix are used for geometry reconstruction. A lumped model is used to model the flow resistance of the branches that are cut off from the truncated pathway. Moreover, since the nanoparticle binding process is stochastic in nature, a binding probability function is used to simplify the carrier attachment and detachment processes. The stitched realistic and artificial geometries coupled with the lumped model at the unresolved outlets are used to resolve the flow field within the truncated arterial tree. Then, the biodistribution of 200 nm, 700 nm and 2 µm particles at different vessel generations is studied. At the end, 0.2–0.5% nanocarrier deposition is predicted during one time passage of drug carriers through pulmonary vascular tree. Our truncated approach enabled us to efficiently model hemodynamics and accordingly particle distribution in a complex 3D vasculature providing a simple, yet efficient predictive tool to study drug delivery at organ level.     

In Vivo Quantitative Microcomputed Tomographic Analysis of Vasculature and Organs in a Normal and Diseased Mouse Model

Non-bone in vivo micro-CT imaging has many potential applications for preclinical evaluation. Specifically, the in vivo quantification of changes in the vascular network and organ morphology in small animals, associated with the emergence and progression of diseases like bone fracture, inflammation and cancer, would be critical to the development and evaluation of new therapies for the same. However, there are few published papers describing the in vivo vascular imaging in small animals, due to technical challenges, such as low image quality and low vessel contrast in surrounding tissues. These studies have primarily focused on lung, cardiovascular and brain imaging. In vivo vascular imaging of mouse hind limbs has not been reported. We have developed an in vivo CT imaging technique to visualize and quantify vasculature and organ structure in disease models, with the goal of improved quality images. With 1–2 minutes scanning by a high speed in vivo micro-CT scanner (Quantum CT), and injection of a highly efficient contrast agent (Exitron nano 12000), vasculature and organ structure were semi-automatically segmented and quantified via image analysis software (Analyze). Vessels of the head and hind limbs, and organs like the heart, liver, kidneys and spleen were visualized and segmented from density maps. In a mouse model of bone metastasis, neoangiogenesis was observed, and associated changes to vessel morphology were computed, along with associated enlargement of the spleen. The in vivo CT image quality, voxel size down to 20 μm, is sufficient to visualize and quantify mouse vascular morphology. With this technique, in vivo vascular monitoring becomes feasible for the preclinical evaluation of small animal disease models.     

Cerebral Neovascularization and Remodeling Patterns in Two Different Models of Type 2 Diabetes

We previously reported intense pial cerebral collateralization and arteriogenesis in a mild and lean model of type 2 diabetes (T2D), Goto-Kakizaki (GK) rats. Increased cerebral neovascularization differed regionally and was associated with poor vessel wall maturity. Building upon these findings, the goals of this study were to determine whether a) glycemic control prevents this erratic cerebral neovascularization in the GK model, and b) this pathological neovascularization pattern occurs in Lepr^db/db model, which is the most commonly used model of T2D for studies involving cerebral complications of diabetes. Vascular volume, surface area and structural parameters including microvessel/macrovessel ratio, non-FITC (fluorescein) perfusing vessel abundance, vessel tortuosity, and branch density were measured by 3D reconstruction of FITC stained vasculature in GK rats or Lepr^db/db mice. GK rats exhibited an increase in all of these parameters, which were prevented by glycemic control with metformin. In Lepr^db/db mice, microvascular density was increased but there was no change in nonFITC-perfusing vessels. Increased PA branch density was associated with reduced branch diameter. These results suggest that T2D leads to cerebral neovascularization and remodeling but some structural characteristics of newly formed vessels differ between these models of T2D. The prevention of dysfunctional cerebral neovascularization by early glucose control suggests that hyperglycemia is a mediator of this response.     

Automatic extraction of nuclei centroids of mouse embryonic cells from fluorescence microscopy images

Accurate identification of cell nuclei and their tracking using three dimensional (3D) microscopic images is a demanding task in many biological studies. Manual identification of nuclei centroids from images is an error-prone task, sometimes impossible to accomplish due to low contrast and the presence of noise. Nonetheless, only a few methods are available for 3D bioimaging applications, which sharply contrast with 2D analysis, where many methods already exist. In addition, most methods essentially adopt segmentation for which a reliable solution is still unknown, especially for 3D bio-images having juxtaposed cells. In this work, we propose a new method that can directly extract nuclei centroids from fluorescence microscopy images. This method involves three steps: (i) Pre-processing, (ii) Local enhancement, and (iii) Centroid extraction. The first step includes two variations: first variation (Variant-1) uses the whole 3D pre-processed image, whereas the second one (Variant-2) modifies the preprocessed image to the candidate regions or the candidate hybrid image for further processing. At the second step, a multiscale cube filtering is employed in order to locally enhance the pre-processed image. Centroid extraction in the third step consists of three stages. In Stage-1, we compute a local characteristic ratio at every voxel and extract local maxima regions as candidate centroids using a ratio threshold. Stage-2 processing removes spurious centroids from Stage-1 results by analyzing shapes of intensity profiles from the enhanced image. An iterative procedure based on the nearest neighborhood principle is then proposed to combine if there are fragmented nuclei. Both qualitative and quantitative analyses on a set of 100 images of 3D mouse embryo are performed. Investigations reveal a promising achievement of the technique presented in terms of average sensitivity and precision (i.e., 88.04% and 91.30% for Variant-1; 86.19% and 95.00% for Variant-2), when compared with an existing method (86.06% and 90.11%), originally developed for analyzing C. elegans images.     

Neuroanatomical phenotyping of the mouse brain with three-dimensional autofluorescence imaging

The structural organization of the brain is important for normal brain function and is critical to understand in order to evaluate changes that occur during disease processes. Three-dimensional (3D) imaging of the mouse brain is necessary to appreciate the spatial context of structures within the brain. In addition, the small scale of many brain structures necessitates resolution at the ～10 μm scale. 3D optical imaging techniques, such as optical projection tomography (OPT), have the ability to image intact large specimens (1 cm(3)) with ～5 μm resolution. In this work we assessed the potential of autofluorescence optical imaging methods, and specifically OPT, for phenotyping the mouse brain. We found that both specimen size and fixation methods affected the quality of the OPT image. Based on these findings we developed a specimen preparation method to improve the images. Using this method we assessed the potential of optical imaging for phenotyping. Phenotypic differences between wild-type male and female mice were quantified using computer-automated methods. We found that optical imaging of the endogenous autofluorescence in the mouse brain allows for 3D characterization of neuroanatomy and detailed analysis of brain phenotypes. This will be a powerful tool for understanding mouse models of disease and development and is a technology that fits easily within the workflow of biology and neuroscience labs.     

Quantitation of 3D ureteric branching morphogenesis in cultured embryonic mouse kidney

The growth and branching of the epithelial ureteric tree is critical for development of the permanent kidney (metanephros). Current methods of analysis of ureteric branching are mostly qualitative. We have developed a method for measuring the length of individual branches, and thereby the total length of the ureteric tree in 3 dimensions (3D). The method involves confocal microscopy of whole-mount immunostained metanephroi and computer-based image segmentation, skeletonisation and measurement. The algorithm performs semi-automatic segmentation of a set of confocal images and skeletonisation of the resulting binary object. Length measurements and number of branch points are automatically obtained. The final representation can be reconstructed providing a fully rotating 3D perspective of the skeletonised tree. After 36 h culture of E12 mouse metanephroi, the total length of the ureteric tree was 6103 +/- 291 microm (mean +/- SD), a four-fold increase compared with metanephroi cultured for just 6 h (1522 +/- 149 microm). Ureteric duct length increased at a rate of 153 microm/h over the first 30 h period and was maximal between 18 and 24 h at 325 microm/h. The distribution of branch lengths at the six time points studied was similar, suggesting tight control of ureteric lengthening and branching. This method will be of use in analysing ureteric growth in kidneys cultured in the presence of specific molecules suspected of regulating ureteric growth. The method can also be used to analyse in vivo kidneys and to quantify branching morphogenesis in other developing organs.     

Three-dimensional spatial reconstruction of coronary arteries based on fusion of intravascular optical coherence tomography and coronary angiography

We present a three-dimensional (3D) spatial reconstruction of coronary arteries based on fusion of intravascular optical coherence tomography (IVOCT) and digital subtraction angiography (DSA). Centerline of vessel in DSA images is exacted by multi-scale filtering, adaptive segmentation, morphology thinning and Dijkstra's shortest path algorithm. We apply the cross-correction between lumen shapes of IVOCT and DSA images and match their stenosis positions to realize co-registration. By matching the location and tangent direction of the vessel centerline of DSA images and segmented lumen coordinates of IVOCT along pullback path, 3D spatial models of vessel lumen are reconstructed. Using 1121 distinct positions selected from eight vessels, the correlation coefficient between 3D IVOCT model and DSA image in measuring lumen radius is 0.94% and 97.7% of the positions fall within the limit of agreement by Bland–Altman analysis, which means that the 3D spatial reconstruction IVOCT models and DSA images have high matching level.     

Flow and pressure distributions in vascular networks consisting of distensible vessels

We examine the influence of vessel distensibility on the fraction of the total network flow passing through each vessel of a model vascular network. An exact computational methodology is developed yielding an analytic proof. For a class of structurally heterogeneous asymmetric vascular networks, if all the individual vessels share a common distensibility relation when the total network flow is changed, this methodology proves that each vessel will continue to receive the same fraction of the total network flow. This constant flow partitioning occurs despite a redistribution of pressures, which may result in a decrease in the diameter of one and an increase in the diameter of the other of two vessels having a common diameter at a common pressure. This theoretical observation, taken along with published experimental observations on pulmonary vessel distensibilities, suggests that vessel diameter-independent distensibility in the pulmonary vasculature may be an evolutionary adaptation for preserving the spatial distribution of pulmonary blood flow in the face of large variations in cardiac output.     

An integrated geometric modelling framework for patient-specific computational haemodynamic study on wide-ranged vascular network

Patient-specific haemodynamic computations have been used as an effective tool in researches on cardiovascular disease associated with haemodynamics such as atherosclerosis and aneurysm. Recent development of computer resource has enabled 3D haemodynamic computations in wide-spread arterial network but there are still difficulties in modelling vascular geometry because of noise and limited resolution in medical images. In this paper, an integrated framework to model an arterial network tree for patient-specific computational haemodynamic study is developed. With this framework, 3D vascular geometry reconstruction of an arterial network and quantification of its geometric feature are aimed. The combination of 3D haemodynamic computation and vascular morphology quantification helps better understand the relationship between vascular morphology and haemodynamic force behind 'geometric risk factor' for cardiovascular diseases. The proposed method is applied to an intracranial arterial network to demonstrate its accuracy and effectiveness. The results are compared with the marching-cubes (MC) method. The comparison shows that the present modelling method can reconstruct a wide-ranged vascular network anatomically more accurate than the MC method, particularly in peripheral circulation where the image resolution is low in comparison to the vessel diameter, because of the recognition of an arterial network connectivity based on its centreline.     

An X-ray microtomography-based method for detailed analysis of the three-dimensional morphology of fungal pellets

Filamentous fungi are widely used in the production of biotechnological compounds. Since their morphology is strongly linked to productivity, it is a key parameter in industrial biotechnology. However, identifying the morphological properties of filamentous fungi is challenging. Owing to a lack of appropriate methods, the detailed three-dimensional morphology of filamentous pellets remains unexplored. In the present study, we used state-of-the-art X-ray microtomography (µCT) to develop a new method for detailed characterization of fungal pellets. µCT measurements were performed using freeze-dried pellets obtained from submerged cultivations. Three-dimensional images were generated and analyzed to locate and quantify hyphal material, tips, and branches. As a result, morphological properties including hyphal length, tip number, branch number, hyphal growth unit, porosity, and hyphal average diameter were ascertained. To validate the potential of the new method, two fungal pellets were studied—one from Aspergillus niger and the other from Penicillium chrysogenum. We show here that µCT analysis is a promising tool to study the three-dimensional structure of pellet-forming filamentous microorganisms in utmost detail. The knowledge gained can be used to understand and thus optimize pellet structures by means of appropriate process or genetic control in biotechnological applications.     

Distribution of vessel size,vessel density and xylem conducting efficiency within a crown of silver birch (<Emphasis Type="Italic">Betula pendula</Emphasis>)

Spatial patterns in vessel diameter, vessel density and xylem conducting efficiency within a crown were examined in closed-canopy trees of silver birch (Betula pendula). The variation in anatomical and hydraulic characteristics of branches was considered from three perspectives: vertically within a crown (lower, middle and upper crown), radially along main branches (proximal, middle and distal part), and with respect to branch orders (first-, second- and third-order branches). Hydraulically weighted mean diameter of vessels (D _h) and theoretical specific conductivity of the xylem (k _t) exhibited no vertical trend within the tree crown, whereas leaf-specific conductivity of the xylem (LSC_t) decreased acropetally. Variation in LSC_t was governed by sapwood area to leaf area ratio (Huber value) rather than by changes in xylem anatomy. The acropetal increase in soil-to-leaf conductance (G _T) within the birch canopy is attributable to longer path length within the lower-crown branches and higher hydraulic resistance of the shade leaves. D _h, k _t and LSC_t decreased, while vessel density (VD) and relative area of vessel lumina (VA) increased distally along main branches. A strong negative relationship between vessel diameter and VD implies a trade-off between hydraulic efficiency and mechanical stability of xylem. D _h and VD combined explained 85.4% of the total variation of k _t in the regression model applied to the whole data set. Xylem in fast-growing branches (primary branches) had greater area of vessel lumina per unit cross-sectional area of sapwood, resulting in a positive relationship between branch radial growth rate and k _t. D _h, k _t and LSC_t decreased, whereas VD increased with increasing branch order. This pattern promotes the hydraulic dominance of primary branches over the secondary branches and their dominance over tertiary branches. In this way crown architecture contributes to preferential water flow along the main axes, potentially providing better water supply for the branch apical bud and foliage located in the outer, better-insolated part of the crown.     

Retinal Vessel Width Measurement at Branchings Using an Improved Electric Field Theory-Based Graph Approach

The retinal vessel width relationship at vessel branch points in fundus images is an important biomarker of retinal and systemic disease. We propose a fully automatic method to measure the vessel widths at branch points in fundus images. The method is a graph-based method, in which a graph construction method based on electric field theory is applied which specifically deals with complex branching patterns. The vessel centerline image is used as the initial segmentation of the graph. Branching points are detected on the vessel centerline image using a set of detection kernels. Crossing points are distinguished from branch points and excluded. The electric field based graph method is applied to construct the graph. This method is inspired by the non-intersecting force lines in an electric field. At last, the method is further improved to give a consistent vessel width measurement for the whole vessel tree. The algorithm was validated on 100 artery branchings and 100 vein branchings selected from 50 fundus images by comparing with vessel width measurements from two human experts.     

A Modified Verhoeff-Van Gieson Elastin Histochemical Stain to Enable Pulmonary Arterial Hypertension Model Characterization

Optimal histochemical staining is critical to ensure excellent quality stained sections to enable light microscopic and histomorphometric image analysis. Verhoeff-van Gieson is the most widely used histochemical stain for the visualization of vascular elastic fibers. However, it is notoriously difficult to differentiate fine elastic fibers of small vasculature to enable histomorphometric image analysis, especially in organs such as the lung. A tissue fixation procedure of 10% neutral buffered formalin with subsequent fixation in 70% ethanol further compounds the problem of small vessel staining and identification. Therefore, a modified Verhoeff’s elastin stain was developed as a reliable method to optimally highlight the internal and external elastic laminae of small arteries (50-100 µm external diameter) and intra-acinar vessels (10-50 µm external diameter) in 3 µm thick lung tissue sections from models of pulmonary arterial hypertension. This modified Verhoeff’s elastin stain demonstrated well-defined staining of fine elastic fibers of pulmonary blood vessels enabling subsequent histomorphometric image analysis of vessel wall thickness in small arteries and intra-acinar vessels. In conclusion, modification of the standard Verhoeff-van Gieson histochemical stain is needed to visualize small caliber vessels’ elastic fibers especially in tissues fixed in 10% neutral buffered formalin followed by additional fixation in 70% ethanol.Key words: Histochemical stain, histomorphology, lung, Verhoeff-van Gieson, elastin     

Suppressive effects of vascular endothelial growth factor-B on tumor growth in a mouse model of pancreatic neuroendocrine tumorigenesis

Background

The family of vascular endothelial growth factors (VEGF) contains key regulators of blood and lymph vessel development, including VEGF-A, -B, -C, -D, and placental growth factor. The role of VEGF-B during physiological or pathological angiogenesis has not yet been conclusively delineated. Herein, we investigate the function of VEGF-B by the generation of mouse models of cancer with transgenic expression of VEGF-B or homozygous deletion of Vegfb.

Methodology/Principal Findings

Ectopic expression of VEGF-B in the insulin-producing β-cells of the pancreas did not alter the abundance or architecture of the islets of Langerhans. The vasculature from transgenic mice exhibited a dilated morphology, but was of similar density as that of wildtype mice. Unexpectedly, we found that transgenic expression of VEGF-B in the RIP1-Tag2 mouse model of pancreatic neuroendocrine tumorigenesis retarded tumor growth. Conversely, RIP1-Tag2 mice deficient for Vegfb presented with larger tumors. No differences in vascular density, perfusion or immune cell infiltration upon altered Vegfb gene dosage were noted. However, VEGF-B acted to increase blood vessel diameter both in normal pancreatic islets and in RIP1-Tag2 tumors.

Conclusions/Significance

Taken together, our results illustrate the differences in biological function between members of the VEGF family, and highlight the necessity of in-depth functional studies of VEGF-B to fully understand the effects of VEGFR-1 inhibitors currently used in the clinic.     

Rapid Reconstruction of 3D Neuronal Morphology from Light Microscopy Images with Augmented Rayburst Sampling

Digital reconstruction of three-dimensional (3D) neuronal morphology from light microscopy images provides a powerful technique for analysis of neural circuits. It is time-consuming to manually perform this process. Thus, efficient computer-assisted approaches are preferable. In this paper, we present an innovative method for the tracing and reconstruction of 3D neuronal morphology from light microscopy images. The method uses a prediction and refinement strategy that is based on exploration of local neuron structural features. We extended the rayburst sampling algorithm to a marching fashion, which starts from a single or a few seed points and marches recursively forward along neurite branches to trace and reconstruct the whole tree-like structure. A local radius-related but size-independent hemispherical sampling was used to predict the neurite centerline and detect branches. Iterative rayburst sampling was performed in the orthogonal plane, to refine the centerline location and to estimate the local radius. We implemented the method in a cooperative 3D interactive visualization-assisted system named flNeuronTool. The source code in C++ and the binaries are freely available at http://sourceforge.net/projects/flneurontool/. We validated and evaluated the proposed method using synthetic data and real datasets from the Digital Reconstruction of Axonal and Dendritic Morphology (DIADEM) challenge. Then, flNeuronTool was applied to mouse brain images acquired with the Micro-Optical Sectioning Tomography (MOST) system, to reconstruct single neurons and local neural circuits. The results showed that the system achieves a reasonable balance between fast speed and acceptable accuracy, which is promising for interactive applications in neuronal image analysis.     

Accurate Image Analysis of the Retina Using Hessian Matrix and Binarisation of Thresholded Entropy with Application of Texture Mapping

In this paper, we demonstrate a comprehensive method for segmenting the retinal vasculature in camera images of the fundus. This is of interest in the area of diagnostics for eye diseases that affect the blood vessels in the eye. In a departure from other state-of-the-art methods, vessels are first pre-grouped together with graph partitioning, using a spectral clustering technique based on morphological features. Local curvature is estimated over the whole image using eigenvalues of Hessian matrix in order to enhance the vessels, which appear as ridges in images of the retina. The result is combined with a binarized image, obtained using a threshold that maximizes entropy, to extract the retinal vessels from the background. Speckle type noise is reduced by applying a connectivity constraint on the extracted curvature based enhanced image. This constraint is varied over the image according to each region''s predominant blood vessel size. The resultant image exhibits the central light reflex of retinal arteries and veins, which prevents the segmentation of whole vessels. To address this, the earlier entropy-based binarization technique is repeated on the original image, but crucially, with a different threshold to incorporate the central reflex vessels. The final segmentation is achieved by combining the segmented vessels with and without central light reflex. We carry out our approach on DRIVE and REVIEW, two publicly available collections of retinal images for research purposes. The obtained results are compared with state-of-the-art methods in the literature using metrics such as sensitivity (true positive rate), selectivity (false positive rate) and accuracy rates for the DRIVE images and measured vessel widths for the REVIEW images. Our approach out-performs the methods in the literature.     

Branches of the left pulmonary artery vascularizing the left upper pulmonary lobe

The studies were carried out on 100 left lungs taken from dead human bodies of both sexes whose age varied from 16 to 80 years. The pulmonary artery and the bronchus were injected with a 65% solution of duracryl and then digested in sulfuric acid. The specimens obtained were examined to determine the number and dimensions of the branches of the left pulmonary artery penetrating into the upper lobe of the left lung as well as the places at which they branch off from this artery. It was found that in most cases 4 branches ramified from the left pulmonary artery. Their length was 30 mm at the most, and their diameter, 12 mm. In about 50% of the cases the branches which penetrated into the lobe were the apicoanterior trunk, the lingular branch and 1 or 2 subsegmental branches, in about 25% of the cases almost all segmental branches penetrated into the lobe separately. In about 20% of the cases the apicoposterior trunk and independent segmental or subsegmental branches were present. Only in about 5% of the cases did the branches under consideration include the apicoposteroanterior trunk and the remaining segmental and subsegmental branches.