Experimental Animals
Female inbred C57Bl/6j (8-weeks old, 20g of body weight) mice were purchased from Envigo, Italy (San Pietro al Natisone, Udine, Italy). Animals were housed five per cage at Chiesi Farmaceutici animal facility. Prior to use, animals were acclimatized for at least 7 days to the local vivarium conditions (room temperature: 20–24°C; relative humidity: 40–70%; 12-h light–dark cycle), having free access to regular rodent chow and softened tap water. Sterile sunflower seeds and hydrogel were supplied as diet integration to prevent excessive body weight loss. All animal experiments described below were approved by the intramural animal-welfare committee for animal experimentation of Chiesi Farmaceutici under protocol number: 809/2020-PR and comply with the European Directive 2010/63 UE, Italian D.Lgs 26/2014 and with the Animal Research: Reporting of In Vivo Experiments (ARRIVE) guidelines [10]. A Visual Analogue scale (0-10) for pain assessment was assessed daily by a designated veterinarian or trained technicians. VAS ≥7 and/or body weight loss ≥ 20% were considered as humane endpoints, as well as signs of dyspnoea or apathy evaluated by a designated veterinarian.
Oropharyngeal administration
The inhalation solution was composed of bleomycin (BLM) hydrochloride (Baxter Oncology GmbH, Germany) and Indocyanine Green (ICG, Sigma-Aldrich, St. Louis, MO) dissolved in 0.9% NaCl, at 1mg/mL and 5mg/mL, respectively. BLM and ICG solutions (BLM+ICG) were mixed (volume ratio 1:1) and administered to 36 mice. ICG solution was diluted with saline (2.5 mg/mL) and administered to 15 mice enrolled as a control group (Sal+ICG).
Administration of either solution (day 0) was performed after mice were slightly anesthetized with 2.5% isoflurane (IsoFlo, Zoetis Inc., New Jersey, USA) delivered in a box. The animals were placed on an intubation platform, positioning their incisor on the wire, the tongue was then pulled out with forceps using a small laryngoscope and 50 µL of solution was administered via oropharyngeal aspiration (OA) using a micropipette [11].
The operator kept the laryngoscope in place for 5-10 seconds to ascertain the correct and complete aspiration of the instilled liquid. Mice were held upright for 10-15 seconds before being placed back into the cage.
Three independent experiments were carried out. The whole experimental procedure is summarized in supplementary figure 1.
Fluorescence imaging
Animals were slightly anesthetized with 2% isoflurane, shaved and imaged using IVIS Lumina II (PerkinElmer Inc., Waltham, MA) [12].
ICG (λEx=790nm, λEm=840nm) fluorescent signal was quantified in calibrated (radiant efficiency) units ([photons/sec/cm2/str]/[µW/cm2]) in Sal+ICG and BLM+ICG mice using the software Living Image® version 4.3.1. (PerkinElmer Inc.).
Mice were imaged in prone position at 0, 7, 14 and 21 days after OA treatment; an average of total fluorescence signal emitted from the chest region was calculated for each mouse.
Micro-Computed Tomography:
Acquisition protocol and image post-processing analysis. Micro-Computed tomography (Micro-CT) was performed at day 7, 14 and 21 using a Quantum GX Micro-CT (PerkinElmer, Inc. Waltham, MA). Mice were slightly anesthetized with 2% isoflurane and images were acquired with the following parameters: 90KV, 88µA, total scan time of 4 minutes (over a total angle of 360°). The ‘high speed’ acquisition protocol was used, and a respiratory gated technique was applied. The entire set of projection radiographs was reconstructed using a filtered back-projection algorithm with a Ram-Lak filter and resulted in two stacks of 512 slices with a nominal resolution of 50 µm.
The reconstructed datasets were analysed with Analyze software (Analyze 12.0; Copyright 1986-2017, Biomedical Imaging Resource, Mayo Clinic, Rochester, MN). Following the analysis protocol, CT scans were filtered and then converted from grey levels to CT numbers (Hounsfield Units- HU).
For the quantitative assessment of the aeration degrees, the total lung volume was extracted from the reconstructed image [11, 13] through manual segmentation. the lung aeration compartments were then determined applying ‘HU preclinical ranges’ [13] dividing the whole parenchyma in normo-aerated ([–860, -435] HU), hypo-aerated ([-435, -121] HU), non-aerated ([-121, +121] HU) and hyper-inflated areas ([-1040, -860] HU), expressed as percentage of total lung volume.
The hypo- and non-aerated lung tissues areas refer to those with a low gas/tissue ratio, which was previously developed to quantify lung fibrosis progression and evaluate antifibrotic drug efficacy [13].
Micro CT Visual categorisation.
Micro-CT scans were visually evaluated by three independent radiologists: four patterns were defined to categorize the lung parenchyma features (see supplementary figure 2). The final decisions reached by consensus.
Healthy lung was defined by homogenous density of lung parenchyma with regular-shaped blood vessels and airways. The reconstructed datasets revealed an entirely detectable lung volume, with normally aerated areas covering around 80% of total lung volume.
Typical BLM fibrosis (TBF) was defined by a patchy pattern with high attenuation areas randomly distributed within the lung parenchyma without visible low attenuation areas [11, 13, 16].
Fibrotic lung associated with moderate airspace enlargement (FAE) was defined by a patchy pattern with heterogeneous density from hazy fibrotic pattern with moderate presence of low attenuation areas, which are difficult to isolate from the surrounding tissue.
Fibrotic lung with severe airspace enlargement (FAE+) was defined by a patchy pattern with heterogeneous density into an overt mosaic pattern characterized by multiple and markedly detectable low-density structures, with well-defined margins appearing in a bubble-like shape in lung parenchyma.
Each Micro-CT scan was attributed to one of these categories at days 7, 14 and 21 and the number of animals assigned to each category was calculated at every time point.
Histological assessment of fibrosis
At 7, 14 and 21 days, after the in-vivo imaging, subsets of 5 (Sal+ICG) and 12 (BLM+ICG) mice were euthanized with an overdose of anaesthetic followed by bleeding from the abdominal aorta. The lungs were removed and inflated with a cannula through the trachea by gentle infusion of 0.6 mL of 10% neutral-buffered formalin and fixed for 24h. For histological assessment, the samples were dehydrated in a graded ethanol series, clarified in xylene, and paraffin embedded. Sections of 5 μm thickness were cut with a rotary microtome (Slee Cut 6062; Slee Medical, Mainz, Germany) and then stained with hematoxylin and eosin (H&E) and Masson’s trichrome (TM), according to the manufacturer’s specifications (Histo-Line Laboratories). The whole slide images (WSI) were acquired by the NanoZoomer S-60 Digital slide scanner (Hamamatsu, Japan) for analysis. Three sections for each lung sample were stained with TM and scored on a scale of 0 to 8 by two independent investigators who were blinded to the treatments. Fibrotic modifications were assessed morphologically and semi-quantitatively graded according to the scale defined by Ashcroft et al. [14] and modified by Hübner et al [15]. The final score was expressed as a mean of individual scores observed across all microscopic fields. To quantify the distribution of pulmonary fibrosis, the Ashcroft scores were graded in three classes of increasing values ranging from 0 to 3 (mild), 4 (moderate), and ≥5 (severe).
Histo-morphometric assessment of airspace enlargements
Alveolar airspace enlargement and alveolar destruction, in accordance with emphysema definition [17], were histomorphometrically evaluated by the parameters of Mean Linear Intercept (MLI) and alveolar airspace area (AAA) [18-23]. The alveolar air spaces area was dimensionally categorized for the analysis as normal (0-100 µm), medium (101-300 µm), large (≥ 300 µm), labeled in white, grey and black, respectively. The distribution of AAA during the time course was expressed as percentage of air content normalized to the total lung parenchyma. See supplementary methods for more details.
Statistics
Statistical analyses were performed using Prism 8 software (GraphPad Software Inc., San Diego, CA, United States). All data are presented as mean ± SEM. Two-way analysis of variance (ANOVA) was performed, followed by Sidak’s multiple comparison post hoc test. Normal distribution was assessed through Shapiro-Wilk test accompanied by visual inspection of QQ-plots. Sample size was calculated with A-priori Power Analysis (G*Power Version 3.1.2) considering Ashcroft Score as endpoint. For all the applied tests, a p-value <0.05 (*) was considered as statistically significant.