Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Quantification of Pulmonary Fibrosis in a Bleomycin Mouse Model Using Automated Histological Image Analysis

Abstract

Current literature on pulmonary fibrosis induced in animal models highlights the need of an accurate, reliable and reproducible histological quantitative analysis. One of the major limits of histological scoring concerns the fact that it is observer-dependent and consequently subject to variability, which may preclude comparative studies between different laboratories. To achieve a reliable and observer-independent quantification of lung fibrosis we developed an automated software histological image analysis performed from digital image of entire lung sections. This automated analysis was compared to standard evaluation methods with regard to its validation as an end-point measure of fibrosis. Lung fibrosis was induced in mice by intratracheal administration of bleomycin (BLM) at 0.25, 0.5, 0.75 and 1 mg/kg. A detailed characterization of BLM-induced fibrosis was performed 14 days after BLM administration using lung function testing, micro-computed tomography and Ashcroft scoring analysis. Quantification of fibrosis by automated analysis was assessed based on pulmonary tissue density measured from thousands of micro-tiles processed from digital images of entire lung sections. Prior to analysis, large bronchi and vessels were manually excluded from the original images. Measurement of fibrosis has been expressed by two indexes: the mean pulmonary tissue density and the high pulmonary tissue density frequency. We showed that tissue density indexes gave access to a very accurate and reliable quantification of morphological changes induced by BLM even for the lowest concentration used (0.25 mg/kg). A reconstructed 2D-image of the entire lung section at high resolution (3.6 μm/pixel) has been performed from tissue density values allowing the visualization of their distribution throughout fibrotic and non-fibrotic regions. A significant correlation (p<0.0001) was found between automated analysis and the above standard evaluation methods. This correlation establishes automated analysis as a novel end-point measure of BLM-induced lung fibrosis in mice, which will be very valuable for future preclinical drug explorations.

Introduction

Pulmonary fibrosis is a severe, often progressive pathologic condition in many respiratory diseases with idiopathic pulmonary fibrosis (IPF) being the most common with high morbidity and mortality [1]. Only two treatment options (nintedanib and pirfenidone) were recently approved both showing a partial reduction of the disease progression but no cure or even regression of the disease [2,3]. Hence, the preclinical evaluation of novel compounds in animal models of lung fibrosis represents further on a critical step in drug development. The animal model of BLM-induced lung fibrosis is widely used to characterize the potential inhibitory effect of newly developed drugs. Although this model has important limitations [4] and does only resemble aspects of pulmonary fibrosis in humans it helped also to characterize the currently available drugs for IPF [5,6,7]. Such characterization requires accurate quantification of morphological changes occurring in pulmonary tissue to determine the severity of pulmonary fibrosis and the efficacy of preclinical drugs. Several methods are currently used and considered as end-point measures in BLM-induced lung fibrosis in mice as micro-computed tomography (micro-CT), magnetic resonance imaging (MRI), gene expression analysis, biochemical analysis (hydroxyproline), blood biomarkers and lung function measurements [818]. Quantification of the severity of fibrosis is also accessed by means of histological methods based either on grading score [10, 1924] or by digital imaging analysis [10, 12, 14, 2527]. Histological scoring has been developed by assigning numerical scores in relationship with the amount of fibrosis in histological samples. Despite the fact that histological scoring is extensively used it does not allow, due to the limited grades of the numerical scale, to achieve a full quantitative analysis of fibrosis. It is noteworthy that the evaluation of scores is observer-dependent and consequently may be subjected to intra- and inter-variability [20]. Histological scoring is performed at a high magnification (x10 or x20) on a limited number of fields (10 to 15) which only accounts for a variable part of the total lung section and is therefore sensitive to bias considering the heterogeneous distribution of fibrosis. Digital imaging analysis is also used to quantify pulmonary morphological changes elicited by BLM-induced fibrosis. Most of the established digital imaging analyses evaluate the severity of BLM-induced fibrosis by quantifying collagen content stained by picrosirius red or to a lesser extent Masson trichrome or even more selectively by immunostaining [12, 14, 2531]. The expression of collagen is determined from the ratio of the stained area versus the total area of the lung section. The main limitation of histological scoring and digital image analyses is related to the fact that quantitative data refer to an average value of morphological changes of collagen content. Due to the high heterogeneity of pulmonary fibrosis this may result in an underestimation of the severity of fibrosis particularly if analyzing mild pathologies. This may hamper an accurate and reliable evaluation of the efficacy of pharmacologically active compounds.

To overcome this main limitation of grading scores and digital imaging methods we developed a suitable automated software image analysis (automated analysis) based on a new assessment of the pulmonary tissue density from entire histological lung sections. Such a new digital imaging method allowed us to determine the distribution of density values in affected and non-affected alveolar parenchyma regions and then to quantify directly and objectively for the first time the BLM-induced fibrotic alterations. With the aim to establish the automated analysis as a quantitative end-point measure of pulmonary fibrosis in the BLM mouse model, the automated analysis has been compared and correlated to Ashcroft scoring [19], in vivo micro-CT analysis and lung function measurements on the same animals. We have provided evidence, from a dose-response study with increasing doses of BLM, that automated analysis allows a robust evaluation of the severity of pulmonary fibrosis. It should be highlighted that automated analysis could be carried out using a single dose administration of BLM, at a lower concentration than those usually used in most bleomycin studies. The lower dose of BLM reduces the stress for the animals used and increases the chance for successful and reliable drug testing.

Materials and Methods

Animals

Eight to 12 week-old C57/BL6 male mice (n = 54) (Charles River, Sulzfeld, Germany), 22 to 30 g, were used throughout the study. Animal husbandry and experimentation was conducted in accordance with German national guidelines and legal regulations and the guidelines of the Association for Accreditation of Laboratory Animal Care. The permission for the animal experiments was issued by the Regierungspräsidium Tübingen, Germany.

BLM-induced pulmonary fibrosis

To facilitate BLM administration animals were anaesthetized for a short period of time using inhaled isoflurane (3–5%). A single dose of BLM (Calbiochem, Darmstadt, Germany) 0.25, 0.5, 0.75 and 1 mg/kg (1–4 U/mg) (n = 12 mice per group) in sterile isotonic saline (50 μL per animal) was intratracheally instilled by means of a 22 gauge plastic cannula (Vasofix, 0.5x25 mm, B. Braun Melsungen, Germany) coupled to a 1 mL syringe to each animal at the start of the study (day 0). The same volume of sterile saline was administered to control mice (n = 6 mice). Body weights were monitored throughout the whole experiment. Non of the animals dropped out of the study due to weight loss till the end of the study at day 14.

Lung function measurement

On day 14 mice were anesthetized with pentobarbital 48 mg/kg combined with xylazine 2.32 mg/kg injected i.p.. After tracheotomy and intubation, a tracheal cannula was connected to a FlexiVent system (SCIREQ, Montreal, PQ, Canada) for pulmonary function testing. To prevent spontaneous breathing the animals received pancuronium bromide 0.64 mg/kg i.v. (Inresa Arzneimittel GmbH, Freiburg, Germany). Ventilation was conducted with 150 breath/min, a tidal volume of 10 mL/kg and an end-expiratory pressure of 3 cm H2O. After recruiting total lung capacity (30 cm H2O) dynamic lung compliance (Cdyn) and forced vital capacity (FVC) was measured.

In vivo micro-CT analysis

Fourteen days after BLM administration in vivo micro-CT analysis of the whole lung was performed by means of a volumetric assessment of BLM-induced changes in x-ray absorption. Animals were anesthetized with isoflurane 1.5% and fixed in prone position. Micro-CT images were acquired on a Quantum FX micro-CT system (Perkin Elmer, Waltham, MA) with cardiac gating (without respiratory gating), using the following parameters: 90 kV; 160 μA; FOV, 60 x 60 x 60 mm; spatial resolution, 0.11 mm, resulting in a total acquisition time of 4–5 minutes. Resulting images were analyzed with MicroView 2.0 software (GE Healthcare, Amersham, UK), which allows a semi-automatic segmentation of left and right lungs. Hounsfield unit (HU) histograms were obtained for left and right lungs using bins of 10-HU width. In the absence of a well-established gold standard to quantify fibrotic changes in rodents with micro-CT and to avoid relying on an arbitrary threshold to identify fibrotic regions [14,15] the HU corresponding to the peak of the HU-histogram for the segmented pixels were used as a measure of the progression of fibrosis. Total lung volume was also computed.

Histopathology score

Following micro-CT analysis and lung function measurements the animals were euthanized at day 14 by an overdose of pentobarbital i.p., lungs were excised and lung wet weight was determined. After cannulation of the trachea the lungs were inflated with paraformaldehyde 4% (PFA) for 20 min at a pressure of 20 cm H2O. The filled lung was then sealed by a ligature and immersed in PFA for at least 24 h. Subsequently PFA fixed lungs were embedded in paraffin following standard procedures. Three serial slices of 3 μm thickness each 50 μm apart of both lung lobes were stained with Masson Trichrome. Fibrotic scoring was performed by a blinded pathologist according to the scale defined by Ashcroft [19]. The slices were analyzed with a microscope using an x10 objective.

Automated histological image analysis (automated analysis)

The total sets of lung slices used for the above standard fibrosis evaluation were scanned at x20 magnification using a NanoZoomer-SQ and digital images of entire lung sections were captured using the NDP.view 2 software (both from Hamamatsu Corporation, Hamamatsu, Japan). To focus the analysis on alveolar parenchyma the walls of large bronchi located in the vicinity of the lung lobes, small bronchi and vessels (diameter >200 μm) associated with alveolar parenchyma as well as their surrounding collagen fibers were manually deleted using Gimp 2.8 software (Free Software Foundation Inc.). Digital images were then reduced from x20 to x2.5 magnification with a pixel size of 3.632 μm allowing both a high-resolution visualization of morphological structure of pulmonary tissue and a very short processing time for the software analysis (<1 sec per entire lung section). The quantification of the BLM-induced fibrotic alterations has been assessed on the basis of pulmonary tissue density. For this purpose Biocellvia has developed a proprietary software program allowing the determination of pulmonary tissue density from thousands of micro-tiles (30–56 μm2) crisscrossing the selected pulmonary tissue of entire lung sections. The density of the pulmonary tissue was evaluated for all individual micro-tiles corresponding to the ratio of the area of the lung tissue inside the micro-tile and the total area of the micro-tile. To quantify the distribution of pulmonary tissue densities they were graded in 20 classes of increasing values in increments of 0.05 (Fig 1). The frequency of tissue density was calculated dividing the number of tissue density values in a given class by the total number of density values in all 20 classes. A distribution of the frequency of tissue density according to their classification was then determined for a comparison between saline-treated control lungs and BLM-treated lungs. To visualize the distribution of density values 2D-reconstructed images of lung sections were composed by assigning pseudocolors to tissue density values according to their classification (Fig 1). These 2D-reconstructions were compared with the Masson trichrome-stained slices.

thumbnail
Fig 1. Representative images of original and 2D-reconstructed lung sections for the quantification and the mapping of pulmonary tissue density by means of digital automatic analysis.

Masson trichrome-stained images (A, E) and their respective 2D-reconstructed images (C, G) correspond to a lung section from saline control (A, C) and BLM-treated (0.75 mg/kg) lungs (E, G). Panels B, D, F, H show details from their original stained (B, F) and 2D-reconstructed (D, H) images at higher magnification. Lung tissue density was determined from thousands of micro-tiles crisscrossing entire lung sections. For mapping thousands tissue density values throughout lung section density values were graded in 20 classes of increasing values (I) and pseudocolours were assigned from blue (low density values) to yellow (high density values) according to their classification. Note that high density values (yellow) were restricted in alveolar parenchyma of BLM-treated lung (G, H) and located in fibrotic lesions evidenced in the respective original stained image (E, F). The frequency of tissue density (I) was determined from the classification of the whole unitary density values obtained in each lung section (A, E). HDFm index corresponds to the sum of the frequencies of the highest tissue density (classes 12 to 20) expressed in fibrotic foci. The mean tissue density (Dm) (J) was evaluated for each lung section from thousands of micro-tiles. Scale bars: 1 mm (A, C, E, G), 100 μm (B, D, F, H).

https://doi.org/10.1371/journal.pone.0170561.g001

Two pulmonary tissue density indexes were defined with the aim to quantify and compare pulmonary fibrotic alterations in saline-treated control lungs and BLM-treated groups: (i) the mean tissue density (Dm) and (ii) the high tissue density frequency (HDFm). Dm corresponds to the mean density value of alveolar parenchyma determined from individual density values of the micro-tiles. Dm can be expressed per lung section, per animal or per group. HDFm was established to specifically quantify fibrotic alterations in alveolar parenchyma. HDFm corresponds to the sum of frequencies of high tissue densities limited to fibrotic alterations. The range of classes allocated to HDFm was established from the saline-treated control group wherein fibrotic alterations were missing and for which high tissue densities frequency per class were less than 1%. In the present study HDFm corresponded to the sum of frequencies ranging from class 12 to 20 (Fig 1).

Statistical analysis

Data are presented as mean ± standard error of mean. Statistical differences between groups were analyzed by one-way analysis of variance (ANOVA) with subsequent Dunnett’s multiple comparison test for all parametric data and Kruskal-Wallis test followed by Dunn’s multiple comparison test for nonparametric data (GraphPad Prism 6.0; GraphPad Software, Inc. La Jolla, CA), Spearman correlation coefficient has been calculated to determine the statistical significance of correlated data. A p-value <0.05 was considered statistically significant.

Results

Quantification of lung fibrosis by automatic histological analysis

Histological lung sections stained with Masson trichrome were scanned and high resolution digital images (3.632 μm/pixel) of entire lung sections were obtained. Following a manual selection of alveolar parenchyma automatic histological analysis was applied to assess the mean pulmonary tissue density (Dm) in the saline-treated control lungs and BLM-treated lungs. Dm was evaluated from thousands of micro-tiles (30–56 μm2) crisscrossing the entire pulmonary tissue of digital images of each lung section.

At day 14, BLM induced a dose-dependent statistical significant increase of Dm in all BLM-treated lungs compared to control lungs (Fig 2A). The observed increase of Dm were in agreement with pathological alterations of lung structures 14 days after administration of BLM in particular thickening of alveolar walls, formation of fibrous bands and fibrous foci (Figs 1 and 3). Such pathological alterations were extensively described in various mouse models of BLM-induced lung fibrosis [4, 8, 19, 20, 22]. It is noteworthy that no statistically significant differences in Dm were found between the BLM-treated lungs of the 0.25 mg/kg and 0.5 mg/kg group and between the 0.75 mg/kg and 1.0 mg/kg group.

thumbnail
Fig 2. Bleomycin (BLM) administration resulted in elevated lung tissue density and fibrotic Ashcroft scoring.

14 days after intratracheal BLM administration at doses of 0.25, 0.5, 0.75 and 1.0 mg/kg lung tissue alterations was assessed on the same whole lung slices stained with Masson trichrome by automatic histological analysis (A) and Ashcroft scoring (B) and in vivo in the same animals by means of micro-CT analysis (C). BLM-treated mice (n = 12/dosing group) were compared to saline control mice (n = 6). Data are shown as mean ± s.e.m. The mean tissue density (Dm) determined by automatic histological analysis (A) was assessed from thousands of micro-tiles covering the whole lung sections of each group of mice. *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001 was considered statistically significant in comparison to the non BLM-treated control group. ns: non-significant.

https://doi.org/10.1371/journal.pone.0170561.g002

thumbnail
Fig 3. Representative original images of Masson trichrome stained lung sections and their corresponding automated 2D-reconstructed images according to their respective lung tissue density values.

Original images (A, C, E, G, I) correspond to representative Masson trichrome stained lung sections from saline control (A) and BLM-treated lungs at 0.25 mg/kg (C), 0.50 mg/kg (E), 0.75 mg/kg (G) and 1 mg/kg (I). Reconstructed-2D lung sections of their respective original images were obtained by grading their tissue density values, obtained from thousands micro-tiles, in 20 classes of increasing values. Pseudocolours were then assigned to tissue densities according to their class from light blue (low density values) to yellow (high density values) (K). For each original image of saline control and BLM-treated lungs the value of their Ascroft score (mean ± s.e.m.) was indicated. Regarding their corresponding 2D-reconstructed images the value of their Dm (mean ± s.e.m.) and HDFm was indicated, respectively. Note that high density values visualized in 2D-reconstructed images were focalized in fibrotic regions evidenced in their corresponding Masson trichrome stained lung sections. Scale bars: 1 mm.

https://doi.org/10.1371/journal.pone.0170561.g003

With the aim to map the topographical distribution of fibrosis, 2D-reconstructed images of entire lung sections were composed on the basis of their tissue density values determined from thousands individual micro-tiles. For this purpose pulmonary tissue densities were graded in 20 classes of increasing values in increments of 0.05 and pseudocolours were assigned from blue (lower density values) to yellow (higher density values) according to their classification (Figs 1 and 3). The mapping of pulmonary tissue density in saline-control lung sections was characterized essentially by low tissue density values (Figs 1 and 3). In lung sections of BLM-treated groups in addition to low density values we observed the presence of high tissue density values which were specifically located in fibrotic regions as evidenced by comparing Masson trichrome images with their respective 2D-reconstructed images (Figs 1 and 3). These fibrotic regions were characterized by small to large fibrotic masses with progressive obliteration of alveolar septa. According to the scale defined by Ascroft et al [19] and Hübner et al [20] these fibrotic regions belonged to elevated scores ranging from 5 to 8.

With the aim to quantify specifically fibrotic lesions in alveolar parenchyma HDFm index was assessed in saline-treated control and BLM-treated lungs. HDFm corresponds to the sum of frequencies of high tissue densities expressed specifically in BLM-treated lungs (Figs 1 and 3). Under these conditions the range of HDFm classes was established from the distribution of the frequency of tissue density in the saline-treated control group. In the present study the threshold of the range was set up at class 12 for which high density frequency value was below 1% (Figs 1, 3 and 4). HDFm value corresponded to the sum of the frequencies of high density values ranged from class 12 to 20. The 2D-reconstructed images of saline control lungs revealed that density values ranged from class 12 to 20 were not related to fibrotic alterations but to a noise background formed by some remaining collagen fibers and/or erythrocytes which were not completely deleted from the manually selection of the alveolar parenchyma (Figs 1 and 3). At day 14, BLM induced a dose-dependent statistical significant increase (p<0.0001) of HDFm in all BLM-treated lungs compared to control lungs (Fig 4B). These quantitative data were in line with those assessed by the Dm index. HDFm turned out to be a greater discriminator of pulmonary fibrosis than Dm. This notion is supported if the percentage of the increase of Dm and HDFm between saline control and BLM-treated lungs is compared. For Dm the BLM-induced increase ranged between 18% to 41% (Fig 2A) but for HDFm between 425% and 1100% (Fig 4B). As expected BLM administration induced fibrotic changes in the alveolar parenchyma resulting in a marked decrease of the percentages of lower frequencies of tissue density and in parallel a marked increase of the higher frequencies. Such changes in the distribution of tissue density elicited by fibrotic alteration of the alveolar parenchyma were accurately quantified by Dm and HDFm indexes independent of the BLM concentration used (Fig 4).

thumbnail
Fig 4. Distribution of tissue density frequency and determination of high tissue density frequency (HDFm) in saline control and bleomycin (BLM)-treated lungs.

Lung tissue density was determined from several thousands of micro-tiles covering the whole lung sections of saline control (n = 6) and BLM-treated (n = 12/dosing group) mice. A: Tissue density values were graded in 20 classes of increasing values (mean ± s.e.m) and their frequency per class was expressed in percent (compared to the total number of density values). B: BLM administration induced a dose-dependent increase of HDFm which corresponded to the sum of the percentage of tissue densities from class 12 to 20. ****p<0.0001 was considered statistically significant in comparison to the non BLM-treated control group.

https://doi.org/10.1371/journal.pone.0170561.g004

Validation of automatic histological analysis in mice lung fibrosis

In order to define the potential value of Dm and HDFm indexes as outcome measures in preclinical studies of lung fibrosis we have explored their correlation with established standards namely Ashcroft scoring, micro-CT imaging and lung function measurements.

The severity of the BLM-induced alterations was examined by comparing the body weight and the lung wet weight of saline-treated control mice with BLM-treated (Fig 5A and 5B). At day 14, BLM at the concentration of 0.75 and 1 mg/kg caused a statistical significant reduction of the body weight gain during the 14 day experimentation period by 6 and 12%, respectively. No significant body weight reduction was observed at the lower BLM concentrations 0.25 and 0.50 mg/kg. Already the lowest dose of 0.25 mg/kg BLM induced a statistical significant increase (p<0.001) of the lung wet weight by 56% (Fig 5B). A dose-dependent increase of lung wet weight was observed for higher BLM concentrations 0.5, 0.75 and 1.0 mg/kg corresponding to 58, 75 and 71% increase, respectively. Such increase of lung wet weight in BLM-treated mice was in agreement with histopathological observations and Ashcroft scoring analysis we performed in Masson trichrome stained lung sections obtained from the same mice. Independent from the BLM dose used all lungs from BLM-treated animals showed peribronchial and subpleural fibrosis. BLM administration induced a dose-dependent statistically significant increase in the Ashcroft scores (Fig 2B).

thumbnail
Fig 5. Bleomycin (BLM) administration resulted in a reduced body weight gain and increased lung wet weight and in parallel in impaired lung functions.

14 days after intratracheal BLM administration at doses of 0.25, 0.5, 0.75 and 1.0 mg/kg the mean ± s.e.m. body weight (A), lung wet weight (B), dynamic lung compliance (Cdyn) (C) and force vital capacity (FVC) (D) were determined in saline control (n = 6) and BLM-treated (n = 12/dosing group) mice. *p<0.05, **p<0.01, ***p<0.001, ****p<0.0001 was considered statistically significant in comparison to the non bleomycin-treated control group.

https://doi.org/10.1371/journal.pone.0170561.g005

BLM administration induced a dose-dependent increase in lung tissue density assessed by micro-CT analysis at day 14 (Figs 2C and 6). At a dose of 0.25, 0.5, 0.75 and 1.0 mg/kg BLM-induced a statistically significant increase of lung tissue density by 18.3% (p<0.05), 133% (p<0.05), 267% (p<0.01) and 366% (p<0.0001), respectively.

thumbnail
Fig 6. Representative micro-CT images of lung sections from saline control and bleomycin-treated mice.

Transverse (top row) and corresponding coronal (bottom row) micro-CT images acquired at 14 days after saline (control) (A, B) and bleomycin administration at the concentration of 0.25 mg/kg (C, D), 0.50 mg/kg (E, F), 0.75 mg/kg (G, H) and 1mg/kg (I, J). Total lung volume is given as mean ± s.e.m. (n = 12/dosing group). **p<0.01, ***p<0.001 was considered statistically significant in comparison to the non bleomycin-treated control group.

https://doi.org/10.1371/journal.pone.0170561.g006

At day 14, BLM administration induced a dose-dependent statistically significant decrease of dynamic lung compliance (Fig 5C) and forced vital capacity (Fig 5D). Even at the lowest BLM dose of 0.25 mg/kg significant differences were detected compared to the control group.

The correlation between the Dm and HDFm indexes with Ashcroft scoring, micro-CT analysis and lung function was evaluated using linear regression analysis (Fig 7). Dm and HDFm mean values demonstrated good correlation with Ascroft score mean values (Dm: r = 0.94, p<0.0001, HDFm: r = 0.96, p<0.0001), mean tissue density values assessed by micro-CT (Dm: r = 0.83, p<0.0001, HDFm: r = 0.84, p<0.0001), dynamic lung compliance (Dm: r = - 0.67, p<0.0001, HDFm: r = - 0.71, p<0.0001) and force vital capacity (Dm: r = - 0.64, p<0.0001, HDFm: r = - 0.68, p<0.0001).

thumbnail
Fig 7. Correlation of tissue density (Dm) and high frequency density (HDFm) indexes with Ashcroft score, micro-CT and lung function measurements in the bleomycin mice model.

The agreement of the Dm and HDFm measurements with those obtained with Ashcroft scoring (A, B), micro-CT analysis (C, D), dynamic lung compliance (Cdyn) (E, F) and force vital capacity (FVC) (G, H) was evaluated using linear regression analysis. Data correspond to mean values per animal of saline control (n = 6) and BLM-treated (0.25, 0.50, 0.75, 1 mg/kg) (n = 12/dosing group) mice. r = Spearman correlation coefficient.

https://doi.org/10.1371/journal.pone.0170561.g007

Discussion

The present study has shown that the automated histological image analysis (automated analysis) allowed a robust histological evaluation of pulmonary fibrosis in the BLM mouse model. Automated analysis was performed by assessing lung tissue density determined from thousands of microtiles crisscrossing sections from the entire lung. With the aim to quantify BLM-induced lung fibrosis two tissue density indexes have been established from individual tissue density values: the mean tissue density (Dm) and the high tissue density frequency (HDFm). In addition, pseudocoloured 2D-reconstructed images of lung slices were composed reflecting the tissue density values. These 2D-reconstructed images enabled visualization and mapping of the distribution of tissue density values in fibrotic and non-fibrotic regions of the lung slices and allowed to select target regions (e.g. fibrotic foci). The potential value of both density indexes as an outcome measure of the severity of lung fibrotic alterations has been validated by their correlation with current standard histological scoring (Ashcroft scores), in vivo imaging technique (micro-CT) and lung function testing.

Because the main bronchi, vascular walls and their surrounding collagen fibers were removed before the images were analyzed the morphological changes assessed by Dm and HDFm indexes were related to fibrotic alterations restricted to the alveolar parenchyma. The Dm index provided access to the mean value of morphological changes occurring in the entire lung sections whereas HDFm index to the mean value of fibrotic alterations. Contrary to the HDFm the Dm index reflects the distribution of fibrotic and non-fibrotic regions in lung sections and may underestimate the percentage of small fibrotic alterations induced either by low BLM concentrations or pharmacological interventions. Although Dm index was less discriminating that HDFm index it allowed to quantify accurately small morphological changes induced by a single dose of BLM (0.25 mg/kg) a concentration lower than the usual dose in rodents [8].

In the present study we sought to validate Dm and HDFm indexes as a direct read out of morphological changes induced in the BLM mouse model comparable to that currently performed by means of established histological scoring methods [19, 20]. The automated analysis and Ashcroft scoring methods were carried out on the same histological sections stained with Masson trichrome to include the changes in collagen in the overall quantification of the morphological changes in the lung. Collagen plays an important role in the fibrotic remodeling of the alveolar parenchyma and accordingly its quantification is frequently assessed by hydroxyproline assays or quantification using picrosirius red or Masson trichrome staining. However, in the present study hydroxyproline method would not allow to pick up subtle fibrotic alterations induced in the alveolar parenchyma [25]. Considering the available histological digital imaging methods they allow an overall quantification of collagen content in the alveolar parenchyma, bronchi and vessels [12, 14] but would not allow to quantify the distribution of collagen content restricted to fibrous and non-fibrous regions of the alveolar parenchyma in view of its correlation with density indexes.

A substantial portion of the quantitative histological analyses of fibrotic lung alterations is still studied by means of scoring methods [19, 20, 22, 32]. The comparison between the automated analysis and the Ashcroft scoring on the same lung sections has led to a highly statistical significant correlation which contributes to validate Dm and HDFm indexes as a direct read out of morphological changes in lung tissue. However, Ashcroft scoring, as well as histological scoring, has major limitations that were not exhibited by automated analysis. Histological scoring is not observer-independent and then subject to intra- and inter-variability [20]. Histological scoring is in general not performed on the entire lung section but on a limited number of microscopy fields, which did not account for the total lung section area and therefore is sensitive to variability considering the heterogeneous distribution of fibrosis. These limitations may also restrict comparative studies between different laboratories.

A highly significant correlation has been established between the automated analysis and the in vivo micro-CT analysis which, in addition with Ashcroft scoring, contributes to validate automated analysis as an end-point measure for quantifying the severity of BLM-induced pulmonary fibrotic alterations in mice. It is noteworthy that in vivo micro-CT analysis is currently one of the key investigative methods for studies exploring pulmonary morphological changes longitudinally [4,8,20,32,33]. Magnetic resonance imaging (MRI) is also a valuable technique in the longitudinal follow-up of fibrosis progression in BLM mouse models [10,12,1416] which might be used in future studies as a complementary approach to the present ex vivo automated analysis. In a recent study it was demonstrated that data obtained with micro-CT analysis are in agreement with those obtained by MRI analysis [16]. Assessment of the pulmonary tissue density by automated analysis and micro-CT requires very different software programs and imaging procedures to achieve the objectives. Accordingly, marked differences have been observed notably for resolution and image acquisition time. In vivo micro-CT resolutions are usually in the range of 25–90 μm/pixel considering multiple constrains due to respiratory and cardiac motion and scanning/reconstruction time [11,3437]. Hence, in vivo micro-CT resolutions are much lower compared to the high resolution (3.5 μm/pixel) used in the automated histological analysis. An improved micro-CT analysis using higher sensitivity involving high resolution scanning [13,34,38,39] might also be used to enhance accurate quantification of low level of fibrotic alterations. High resolutions (1–2 μm/pixel) have also been achieved with micro-CT analysis on ex vivo fixed lungs [34]. However, such resolutions can only be applied on restricted 3D fields of some millimeters. In contrast, the automated analysis is applied to the entire lung section. It should be noted that micro-CT high resolutions are associated with sometimes prohibitively long scan/reconstruction times of several hours per sample precluding high throughput [13,34]. In contrast the automated analysis acquisition/analysis time is much shorter in the range of 1 hour for 100 entire lung sections. Ex vivo micro-CT may also be used for a 3D evaluation of severity of lung fibrosis providing a more sensitive and quantitative measure as well as a marked reduction of scanning time [13].

A potential limitation of the automated analysis approach is related to the exclusion of the large bronchi and vessels manually before starting the automated image processing. Such preparation is simple and requires less than 3 hours for 100 lung sections. However, the manual exclusion of large bronchi and vessels might be subject to bias in their delimitation that might impact on the assessment of Dm index. Nevertheless, considering the elevated surface of lung tissue analyzed per group such impact can be considered insignificant. We also underlined that manual exclusion did not interfere on the evaluation of HDFm index since it was determined exclusively from fibrotic alterations located in the alveolar parenchyma. The manual exclusion of large bronchi and vessels from the automated analysis limits the use of the automated analysis to pharmacological principles not affecting bronchi and large vessels.

In conclusion, the presented automated analysis of images of lung sections from mice treated with bleomycin represents a fully accurate and reliable quantification of bleomycin-induced fibrotic lung alterations. The sensitivity to quantify precisely the fibrotic lung alterations at low bleomycin concentrations is a major improvement. It should be highlighted that using low BLM concentration may be very attractive for future studies since it allows to test the efficacy of compounds in animals in much better shape, eliminating the morbidity, and contributing therefore to a higher reliability of the preclinical drug evaluation. The high correlation between the newly developed automated analysis and standard evaluation methods like Ashcroft scoring, micro-CT and lung function testing establishes it as a new, reliable, accurate and fast alternative to determine fibrotic lung alterations. In line with micro-CT, automated image analysis may be regarded as a new additional or complementary robust measurement for future preclinical drug testing.

Author Contributions

  1. Conceptualization: JCG YJ LW.
  2. Data curation: SK BS DS LW JCG YJ.
  3. Formal analysis: JCG YJ LW.
  4. Investigation: SK BS DS LW JCG YJ.
  5. Methodology: SK BS DS JCG.
  6. Project administration: JCG YJ LW.
  7. Resources: SK BS DS LW JCG YJ.
  8. Software: JCG YJ.
  9. Supervision: JCG YJ LW.
  10. Validation: JCG YJ LW.
  11. Visualization: JCG YJ LW.
  12. Writing – original draft: YJ LW.
  13. Writing – review & editing: YJ LW.

References

  1. 1. Raghu G, Rochwerg B, Zhang Y, Garcia CA, Azuma A, Behr J, et al. (2015) An Official ATS/ERS/JRS/ALAT Clinical Practice Guideline: Treatment of Idiopathic Pulmonary Fibrosis. An Update of the 2011 Clinical Practice Guideline. Am J Respir Crit Care Med 192: 3–19.
  2. 2. Richeldi L, du Bois RM, Raghu G, Azuma A, Brown KK, Costabel U, et al. (2014) Efficacy and safety of nintedanib in idiopathic pulmonary fibrosis. N Engl J Med 370: 2071–2082. pmid:24836310
  3. 3. King TE, Bradford WZ, Castro-Bernardini S, Fagan EA, Glaspole I, Glassberg MK, et al. (2014) A phase 3 trial of pirfenidone in patients with idiopathic pulmonary fibrosis. N Engl J Med 370: 2083–2092. pmid:24836312
  4. 4. Moeller A, Ask K, Warburton D, Gauldie J, Kolb M (2008) The bleomycin animal model: a useful tool to investigate treatment options for idiopathic pulmonary fibrosis? Int J Biochem Cell Biol 40: 362–382. pmid:17936056
  5. 5. Schaefer CJ, Ruhrmund DW, Pan L, Seivert SD, Kossen K (2011) Antifibrotic activities of pirfenidone in animal models. Eur Respir Rev 20: 85–97. pmid:21632796
  6. 6. Wollin L, Maillet I, Quesniaux V, Holweg A, Ryffel B (2014) Antifibrotic and anti-inflammatory activity of the tyrosine kinase inhibitor nintedanib in experimental models of lung fibrosis. J Pharmacol Exp Ther 349: 209–220. pmid:24556663
  7. 7. Borzone G, Moreno R, Urrea R, Meneses M, Oyarzun M, Lisboa C (2001) BLM-induced chronic lung damage does not resemble human idiopathic pulmonary fibrosis. Am J Respir Crit Care Med 163: 1648–1653.
  8. 8. Scotton CJ, Chamber RC (2010) BLM revisited: towards a more representative model of IPF?. Am J Physiol Lung Cell Mol Physiol 299: L439–L441.
  9. 9. Manali ED, Moschos C, Triantafillidou C, Kotanidou A, Psallidas I, Karabela SP, et al. (2011) Static and dynamic mechanics of the murine lung after intratracheal bleomycin. BMC Pulmonary Medicine 11: 33–47. pmid:21627835
  10. 10. Babin AL, Cannet C, Gérard C, Saint-Mezard P, Page CP, Sparrer H, et al. (2012) Bleomycin-Induced Lung Injury in Mice Investigated by MRI: Model Assessment for Target Analysis. Magn Reson Med 67: 499–509. pmid:21656559
  11. 11. De Langhe E, Vande Velde G, Hostens J, Himmelreich U, Emery B, Luyten FP, et al. (2012) Quantification of lung fibrosis and emphysema in mice using automated micro-computed tomography. Plos One 7: e43123. pmid:22912805
  12. 12. Egger C, Cannet C, Gérard C, Jarman E, Jarai G, Feige A, et al. (2013) Administration of bleomycin via the oropharyngeal aspiration route leads to sustained lung fibrosis in mice and rats as quantified by UTE-MRI and histology. Plos One 8: e63432. pmid:23667616
  13. 13. Scotton CJ, Hayes B, Alexander R, Datta A, Forty EJ, Mercer PF, et al. (2013) Ex vivo micro-computed tomography analysis of BLM-induced lung fibrosis for preclinical drug evaluation. Eur Respir 42: 1633–1645
  14. 14. Egger C, Gérard C, Vidotto N, Accart N, Cannet C, Dunbar A, et al. (2014) Lung volume quantified by MRI reflects extracellular-matrix deposition and altered pulmonary function in bleomycin models of fibrosis: effects of SOM230. Am J Physiol Lung Cell Mol Physiol 306: L1024–L1077.
  15. 15. Gammon ST, Foje N, Brewer EM, Owers E, Downs CA, Budde MD, et al. (2014) Preclinical anatomical, molecular, and functional imaging of the lung with multiple modalities. Am J Physiol Lung Cell Mol Physiol 306: L897–L914. pmid:24658139
  16. 16. Vande Velde G, De Langhe E, Poelmans J, Dresselaers T, Lories RJ, Himmelreich U (2014) Cross-validation and comparison of different magnetic resonance imaging protocols with micro-computed tomography and histology in the bleomycin-induced mouse model. Invest Radiol 49: 691–698. pmid:24872004
  17. 17. Fernandez IE, Amarie OV, Mutze K, Königshoff M, Yildirim AO, Eickelberg O (2016) Systematic phenotying and correlation of biomarkers with lung function and histology in lung fibrosis. Am J Physiol Lung Cell Mol Physiol 310: L919–L927. pmid:26993522
  18. 18. Vande Velde G, Poelmans J, De Langhe E, Hillen A, Vanoirbeek J, Himmelreich U, et al. (2016) Longitudinal micro-CT provides biomarkers of lung disease that can be used to assess the effect of therapy in preclinical mouse models, and reveal compensatory changes in lung volume. Disease Models Mechanisms 9: 91–98. pmid:26563390
  19. 19. Ashcroft T, Simpson JM, Timbrell V (1988) Simple method of estimating severity of pulmonary fibrosis on a numerical scale. J Clin Pathol 41: 467–470. pmid:3366935
  20. 20. Hübner R, Gitter W, El Mokhtari NE, Mathiak M, Both M, Bolte H, et al. (2008) Standardized quantification of pulmonary fibrosis in histological samples. BioTechniques 299: L439–L441.
  21. 21. Izumo T, Kondo M, Nagai A (2009) Effects of a leukotriene B4 receptor antagonist on bleomycin-induced pulmonary fibrosis. Eur Respir J 34: 1444–1451. pmid:19460787
  22. 22. Phillips EJ, Peng R, Burns L, Harris P, Garrido R, Tyagi G, et al. (2012) BLM- induced lung fibrosis increases work of breathing in the mouse. Pul Pharmacol Ther 25: 281–285.
  23. 23. Weng T, Poth JS, Karmouty-Quintana H, Garcia-Morales LJ, Melicoff E, Luo F, et al. (2014) Hypoxia-induced deoxycytidine kinase contributes to epithelial proliferation in pulmonary fibrosis. Am J Respir Crit Care Med 190: 1402–1412. pmid:25358054
  24. 24. Bei Y, Hua-Huy T, Nicco C, Duong-Quy S, Le-Dong NN, Tiev KP, et al. (2016) RhoA/Rho-kinase activation promotes lung fibrosis in an animal model of systemic sclerosis. Exp Lung Res 42: 44–55. pmid:26873329
  25. 25. Rydell-Törmanen K, Andréasson K, Hesselstrand R, Risteli J, Heinegard D, Saxne T, et al. (2012) Extracellular matrix alterations and acute inflammation; developing in parallel during early induction of pulmonary fibrosis. Lab Invest 92: 917–925. pmid:22469699
  26. 26. Cicko S, Grimm M, Ayata K, Beckert J, Meyer A, Hossfeld M, et al. (2015) Uridine supplementation exerts anti-inflammatory and anti-fibrotic effects in an animal model of pulmonary fibrosis. Respir Res 16: 105–115. pmid:26369416
  27. 27. Chow LN, Schreiner P, Ng BY, Lo B, Hughes MR, Wilder Scott R, et al. (2016) Impact of a CXCL12/CXCR4 antagonist in bleomycin (BLM) induced pulmonary fibrosis and carbon tetrachloride (CCL4) induced hepatic fibrosis in mice. Plos one 11: e0151765. pmid:26998906
  28. 28. Pickering JG, Boughner DR (1991) Quantitative assessment of the age of fibrotic lesions using polarized light microscopy and digital imaging analysis. Am J Pathol 138: 1225–1231. pmid:2024708
  29. 29. Decologne N, Wettstein G, Kolb M, Margetts P, Garrido C, Camus P, et al. (2010) Bleomycin induces pleural and subpleural fibrosis in the presence of carbon particles. Eur Respir J 35: 176–185. pmid:19574324
  30. 30. Jarman ER, Khambata VS, Cope C, Jones P, Roger J, Ye LY, et al. (2014) An inhibitor of NADPH oxydase-4 attenuates established pulmonary fibrosis in a rodent disease model. Am J Physiol Lung Cell Mol biol 50: 158–169.
  31. 31. Stumm CL, Wettlaufer SH, Jancar S, Peters-Golden M, et al. (2011) Airway remodeling in murine asthma correlates with a defect in PGE2 synthesis by lung fibroblasts. Am J Physiol Lung Cell Mol Physiol 301: L636–L644. pmid:21873451
  32. 32. Ask K, Labiris R, Farkas L, Moeller A, Froese A, Farncombe T, et al. (2008) Comparison between conventional and clinical assessment of experimental lung fibrosis. J Translational Medicine 6: 1–10.
  33. 33. Artaechevarria X, Blanco D, Perez-Martin D, de Biurrun G, Montuenga LM, de Torres JP, et al. (2010) Longitudinal study of a mouse model of chronic pulmonary inflammation using breath hold gated micro-CT. Eur Radiol 20: 2600–2608. pmid:20574632
  34. 34. Vasilescu DM, Knudsen L, Ochs M, Weibel ER, Hoffman EA (2012) Optimized murine lung preparation for detailed structural evaluation via micro-computed tomography. J Appl Physiol 112: 159–166. pmid:21817110
  35. 35. Rodt T, von Falk C, Dettmer S, Halter R, Maus R, Ask K, et al. (2010) Micro-computed tomography of pulmonary fibrosis in mice induced by adenoviral gene transfert of biologically active transforming growth factor-β1. Respir Res 11: 181–188. pmid:21176193
  36. 36. Cavanaugh D, Travis EL, Price RE, Gladish G, White RA, Wang M, et al. (2006) Quantification of BLM-induced lung damage in vivo with micro-computed tomography. Acamedic Radiology 13: 1505–1512.
  37. 37. Ford NL, Martin EL, Lewis JF, Veldhuizen RA, Holdsworth DW, Drangova M (2009) Quantifying lung morphology with respiratory-gated micro-CT in a murine model of emphysema. Phys Med Biol 54: 2121–2130. pmid:19287083
  38. 38. Thiesse J, Namati E, Sieren JC, Smith AR, Reinhardt JM, Hoffman EA, et al. (2010) Lung structure phenotype variation in inbred mouse strains revealed through in vivo micro-CT imaging. J Appl Physiol 109: 1960–1968. pmid:20671036
  39. 39. Hsia CCW, Hyde DM, Ochs M, Weibel ER (2010) An official research policy statement of the american thoracic society/american respiratory society: standards for quantitative assessment of lung structure. Am J Respir Crit Care Med 181: 394–418. pmid:20130146