Design-Based Stereology of the Lung in the Hyperoxic Preterm Rabbit Model of Bronchopulmonary Dysplasia

Mühlfeld, Christian; Schulte, Henri; Jansing, Johanna Christine; Casiraghi, Costanza; Ricci, Francesca; Catozzi, Chiara; Ochs, Matthias; Salomone, Fabrizio; Brandenberger, Christina

doi:https://doi.org/10.1155/2021/4293279

Oxidative Medicine and Cellular Longevity

On this page

Abstract Introduction Methods Results Discussion Conclusions Data Availability Disclosure Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Pathophysiology and Therapeutic Strategies for Perinatal Organ Injury Caused by Disturbed Oxygenation

View this Special Issue

Research Article | Open Access

Volume 2021 | Article ID 4293279 | https://doi.org/10.1155/2021/4293279

Design-Based Stereology of the Lung in the Hyperoxic Preterm Rabbit Model of Bronchopulmonary Dysplasia

Christian Mühlfeld,^1,2Henri Schulte,¹Johanna Christine Jansing,¹Costanza Casiraghi,³Francesca Ricci,³Chiara Catozzi,³Matthias Ochs,^1,4,5Fabrizio Salomone,³and Christina Brandenberger^1,2

Academic Editor: Ana Cipak Gasparovic

Received16 Jun 2021

Accepted03 Sept 2021

Published06 Oct 2021

Abstract

Bronchopulmonary dysplasia (BPD) is a complex condition frequently occurring in preterm newborns, and different animal models are currently used to mimic the pathophysiology of BPD. The comparability of animal models depends on the availability of quantitative data obtained by minimally biased methods. Therefore, the aim of this study was to provide the first design-based stereological analysis of the lungs in the hyperoxia-based model of BPD in the preterm rabbit. Rabbit pups were obtained on gestation day 28 (three days before term) by cesarean section and exposed to normoxic (21% O₂, ) or hyperoxic (95% O₂, ) conditions. After seven days of exposure, lung function testing was performed, and lungs were taken for stereological analysis. In addition, the ratio between pulmonary arterial acceleration and ejection time (PAAT/PAET) was measured. Inspiratory capacity and static compliance were reduced whereas tissue elastance and resistance were increased in hyperoxic animals compared with normoxic controls. Hyperoxic animals showed signs of pulmonary hypertension indicated by the decreased PAAT/PAET ratio. In hyperoxic animals, the number of alveoli and the alveolar surface area were reduced by one-third or by approximately 50% of control values, respectively. However, neither the mean linear intercept length nor the mean alveolar volume was significantly different between both groups. Hyperoxic pups had thickened alveolar septa and intra-alveolar accumulation of edema fluid and inflammatory cells. Nonparenchymal blood vessels had thickened walls, enlarged perivascular space, and smaller lumen in hyperoxic rabbits in comparison with normoxic ones. In conclusion, the findings are in line with the pathological features of human BPD. The stereological data may serve as a reference to compare this model with BPD models in other species or future therapeutic interventions.

1. Introduction

Bronchopulmonary dysplasia (BPD) is a chronic human lung disease associated with preterm birth. In the early days, BPD was described to be the result of a combination of respiratory distress syndrome and the effects of therapeutic interventions, including ventilation and oxygen supply [1]. With the possibility of treating preterm babies with exogenous surfactant, this form of BPD (“old BPD”) has decreased; however, the progress in neonatal care has helped to keep alive very immature neonates, born between gestation weeks 24 and 28 [2]. At this stage, lung development is incomplete and needs to progress to gain full functional capacity, while at the same time, the immature lung is required to work as a gas exchanger. The conditions for further lung development are aggravated by perinatal inflammation/infection, mechanical ventilation, and oxygen supply [3]. As a result, lung development gets disrupted which is characterized by a decrease in alveolarization and (micro-)vascular maturation as well as interstitial changes, including deposition of fibrotic tissue [4, 5]. In order to distinguish this form from the previously known condition, it has been termed “new BPD.” The pathophysiology of the new BPD, in particular, the mechanisms by which alveolarization is halted, is still not fully understood. In addition, due to the chronic nature of BPD, children with this disease are affected lifelong by the consequences of disrupted lung development, and therapeutic approaches are needed to promote extrauterine lung development in extremely preterm newborns [2, 6–8].

For these purposes, preclinical experiments addressing the pathophysiology of arrested lung development and the effects of new therapeutic interventions are essential [9, 10]. Currently, several animal models are used in BPD research to mimic parts of the pathophysiological sequelae of human preterm birth. The most frequently used is the term mouse [11] or rat [12] model of BPD in which developmental arrest is achieved by exposure of the pups to a hyperoxic environment. At birth, the lungs of mice and rats are physiologically in the saccular stage of lung development which corresponds to the 24^th to 30^th gestational week in humans [11, 13]. Hyperoxia disrupts the alveolarization process which takes place after birth in rodents [14]. However, lung biology between humans and mice differs greatly, and the translation of therapeutics tested in rodent models into the clinical situation has led to disappointing results [9]. In contrast, preterm models of BPD are used in chronically ventilated lambs and baboons [15]. These models are valuable as they are much closer to the situation of human preterm newborns; however, they are not only complex, technically demanding, and expensive but also ethically challenging. In particular, the use of the baboon model is restricted to questions that cannot be addressed otherwise, e.g., last testing of a therapeutic intervention in an animal model before entering a clinical trial [16]. As a bridge between the hyperoxia-based and preterm animal models, the rabbit model of preterm birth and subsequent exposure to hyperoxic conditions has been introduced. In this model, the pups are born three days before term by a cesarean section and then exposed to 85%-95% ambient oxygen conditions [17, 18].

Morphometric methods are extremely valuable to characterize the lungs of different animal models in quantitative terms [19, 20]. At best, the resulting data can be used for statistical testing of robust parameters allowing comparisons among different animal models and therapeutic interventions. The gold standard of lung morphometry is design-based stereology, a set of methods that is aimed at obtaining, in theory, unbiased and, in practice, minimally biased quantitative data about the volume, surface area, length, and number of certain structures [21]. While useful stereological estimations have been extensively performed in the mouse model of BPD [22–24], a formal stereological characterization of the rabbit preterm/hyperoxia model is missing so far. Although morphometric analyses in the rabbit model have been carried out [25, 26], comparison of the pathophysiological events taking place in different models requires the use of a comparable morphometric methodology, i.e., design-based stereology. Therefore, the present study used a formal stereological approach according to current American Thoracic Society/European Respiratory Society standards [21] to analyse pulmonary changes in the rabbit model of BPD and to provide reference data for future studies using this animal model. The stereological data were further complemented by lung function testing to allow structure-function analysis.

2. Material and Methods

2.1. BPD Rabbit Model

Time-mated New Zealand white rabbits were obtained from Charles River Laboratories (France). All experimental procedures were approved by the local animal ethics committee and complied with the standard European regulations on animal research (no. 744/2017). For cesarean section at the 28^th day of gestation (early saccular stage of lung development, term = 31^st days of gestation), the does were sedated with intramuscular (i.m.) medetomidine 0.25 mg/kg (Domitor®, Orion Pharma, Finland). Ten minutes later, they received i.m. 20 mg/kg of ketamine (Nimatek®, Dechra, Italy) and 5 mg/kg of xylazine (Rompun®, Bayer, Germany). Subsequently, the does were euthanized with an overdose of pentothal sodium (100 mg/kg, MSD Animal Health, USA). After death occurred, the abdomen was immediately opened, and the uterus exposed to extract all pups through hysterectomy. At delivery, pups were dried, stimulated, and divided into two groups similarly as described by Salaets et al. [27]: the normoxia group where pups were housed in 21% oxygen for 7 days and the hyperoxia group where pups were housed in 95% oxygen for 7 days. All pups were placed in custom-made incubators (Okolab, Italy) at 32°C and 50-60% relative humidity. Oxygen level, humidity, and temperature were continuously monitored for the whole experiment. Pups were fed twice daily via an orogastric tube with increasing quantities of a milk replacer (Day One®, Protein 30%, Fat 50%; Fox Valley, USA), supplemented with vitamins and probiotics (Bio-Lapis®; Probiotics International Ltd, UK), whereas immunoglobulins were added only during the first 2 days of life (Col-o-Cat®, SanoBest, Netherlands). Furthermore, intramuscular vitamin K was injected on postnatal day 2 (0.002 mg/kg, Konakion Paediatric; Roche).

2.2. Pulmonary Artery Microultrasound Doppler Analysis

After 7 days of exposure with normoxic or hyperoxic air, the pups were anesthetized with isoflurane (2.5% in pure oxygen) and were placed in the supine position on the Vevo Imaging station (VisualSonics, Canada). Two-dimensional images of the pulmonary artery were obtained from a parasternal short axis view at the level of the aortic valve using the MX-550D transducer (40 MHz, broadband frequency 22 MHz-55 MHz, Vevo 3100). Measurements were performed offline (Vevo LAB software package V3.0). The following variables were measured: pulmonary artery acceleration time (PAT), ejection time (PET), and PAAT/PAET ratio. An average of three cardiac cycles was used for the analyses.

2.3. Necropsy and Lung Function Measurements

At the end of the experimental period, pups were anesthetized with ketamine (35 mg/kg) and xylazine (6 mg/kg) i.p., and lung function testing was performed using a forced oscillation technique with the flexiVent system (flexiVent; SCIREQ, Montreal, Canada) [28]. The analysis was done with the flexiVent module 2 that is suited for rodents and other small animals weighing not more than 80-85 g and that has been used previously in the preterm rabbit model by other research groups [25, 27, 28]. The following parameters were assessed: airway resistance (Rn), tissue damping (resistance, G), and tissue elasticity (elastance, H) using Primewave-8 forced oscillation and the inspiratory capacity (IC) and static compliance (Cst) using the pressure-volume perturbation. At the end of the lung function measurements, pups were euthanized with an overdose of pentothal sodium (100 mg/kg, MSD Animal Health, USA) i.p., and whole lung samples were harvested for stereological analysis.

2.4. Lung Tissue Preparation

Lungs of preterm born rabbits were fixed by instillation at a pressure of 25 cmH₂O with a fixative containing 1.5% glutaraldehyde and 1.5% paraformaldehyde in 0.15 M HEPES buffer. Our inflation fixation protocol was tested and optimized in advance for inflation fixation in preterm rabbits. After at least 24 h submersion in fixative, lungs were further processed for stereological analysis. Prior to embedding in Technovit resin, lung volume measurements were assessed by volume displacement (based on Archimedes’ principle), and lung tissue was sampled by applying systematic uniform random sampling [29]. This was done by cutting the lung tissue into 2 mm thick slices, and every other slice, with a random start, was embedded in Technovit resin as described previously [29, 30]. Resin-embedded lung tissue was cut either into 1.5 μm thick sections and stained with toluidine blue or into two sequential sections with a distance of 4.5 μm thickness and stained with orcein for disector analysis and alveolar counting [31]. Histological slides with tissue sections were then scanned and digitalized with a microscopic slide scanner (AxioScan.Z1, Zeiss, Germany) at a 20x objective lens magnification.

2.5. Stereological Analysis

Stereological analysis was done with the newCAST™ software (Visiopharm®, Denmark) on digitalized slides. For the estimation of parenchymal, nonparenchymal, and ductal airspace volumes, automated systematic uniform random subsampling (SURS) was performed at a 5x objective lens magnification and with a sampling fraction of 5-7% to collect the images for further analysis. A point grid consisting of 36 points with a 4-fold subsampling fraction was superimposed over the collected images and used to estimate the respective volume densities as described previously in detail [30]. Respective total volumes were assessed by multiplying the densities with the lung volume.

For the estimation of alveolar and ductal airspace, septal volume, septal thickness, and septal surface area, an automated SURS was performed at a 20x objective lens magnification and with a sampling fraction of 0.5-1%. A test line grid with 12 lines and a length per point of 9.37 μm was superimposed over the images, and points and intersections were counted. Respective volume and surface densities were estimated as described previously [30], and total volume and surface areas were obtained by multiplying the densities with the parenchymal lung volume. The presence of alveolar edema was also assessed as an additional parameter within the alveolar airspace.

For the estimation of the alveolar number, the consecutive tissue sections were aligned in the Visiopharm software. An automated SURS was performed at a 20x objective lens magnification with a sampling fraction of 0.9-1.8%, resulting in approximately 100 to 140 image pairs per lung. A counting frame with an area of 35,570.40 μm² was superimposed over the disector image pairs, and the number of alveoli was estimated as described previously [31, 32]. Mean linear intercept lengths (Lm) were calculated by dividing four times the total parenchymal airspace volume by the septal surface area as described previously [33].

The vascular volume, vascular wall thickness, and perivascular area were estimated using a sampling fraction of 4-9% with approximately 170-200 images at a 10x objective lens magnification. A length grid with lines and a length per point of 18.74 μm was used to calculate the vascular volume and surface area in vessels with a diameter of more than 25 μm. A point grid with points was used to estimate the mean volume and thickness of the perivascular area and vascular wall in vessels with diameter of more than 25 μm. Calculations were done as described previously [34].

2.6. Statistical Analysis

A number of 9-10 rabbits per experimental group were included in the study. Lungs that did not properly inflate during organ harvesting due to technical issues were excluded from morphometric analysis, resulting in stereological analysis of 8 lungs per experimental group. Statistical comparison between the two groups was done with Student’s -test using the SigmaPlot® software (SYSTAT® Software Inc.). If normality testing failed, the Mann–Whitney rank sum test was used. Results were considered significant if .

3. Results

3.1. Lung Function Measurements and Pulmonary Artery Analysis with Microultrasound

The pups were weighed, and pulmonary artery microultrasound Doppler analysis and lung function measurements were performed prior to sacrifice. The average body weight did not differ between the two experimental groups after 7 days of exposure to hyperoxic or normoxic air and was 50 g (±4.9 g) for the normoxia-exposed pups and 49 g (±9 g) for the hyperoxia-exposed pups. However, significant differences were found in inspiratory capacity (IC), static compliance (Cst), tissue resistance (G), and tissue elastance (H) between the two experimental groups (Figures 1(a)–1(d)). While IC and Cst were significantly reduced in hyperoxia- compared with normoxia-exposed pups, tissue elastance and resistance were significantly elevated in the hyperoxic group. No differences between the groups were detected in airway resistance (Rn, Figure 1(e)), suggesting that hyperoxia treatment affected the pulmonary parenchyma but not the airways. Microultrasound Doppler analysis of the pulmonary artery furthermore revealed a decreased ratio of pulmonary artery acceleration time to pulmonary artery ejection time (PAAT/PAET) with hyperoxia treatment, suggesting development of pulmonary hypertension (Figure 1(f)). The data hence provides evidence of impaired lung function and pulmonary hypertension due to hyperoxia treatment without having an impact on body weight within the experimental duration.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 1

Effects of hyperoxia on lung function measurement, pulmonary arterial acceleration, and ejection time. Lung function measurements were assessed with a rodent flexiVent, measuring inspiratory capacity (a), static compliance (b), tissue resistance (c), tissue elastance (d), and airway resistance (Rn). Pulmonary artery acceleration per ejection time ratio (PAAT/PAET) was assessed with Doppler microultrasound analysis. Lines indicate statistical differences between groups (), assessed with -test or Mann–Whitney rank sum test if normality distribution failed.

3.2. Structural Changes in the Lung Parenchyma

The lung volume, assessed with volume displacement after instillation fixation, was also significantly reduced with exposure to hyperoxia (Figure 2(a)) — similarly as the IC and the Cst. This was primarily due to changes in the parenchymal compartment of the lung that declined with hyperoxia treatment (Figure 2(b)). The total volume of the nonparenchymal lung tissue (conducting airways, vasculature, and peribronchiolar and perivascular tissues), however, was not affected by hyperoxia exposure (Figure 2(c)). Representative images providing an overview of the alterations in pulmonary structure and histopathology are shown in Figures 2(d) (normoxia) and 2(e) (hyperoxia). Stereological analysis of the parenchyma further revealed that the alveolar number as well as the septal surface area significantly declined in hyperoxia-treated rabbits (Figures 3(a) and 3(b)). Similarly, the ductal and alveolar airspace volumes were significantly reduced with hyperoxia treatment (Figures 3(c) and 3(d)). However, there was no significant difference in the mean alveolar volume and mean linear intercept lengths (Lm) between the two treatment groups (Figures 3(e) and 3(f)). The total septal volume (Figure 3(g)) did not decline with hyperoxia exposure, which, however, was mostly due to an increased septal thickening in the hyperoxia-treated lungs (Figure 3(h)). Figure 4 shows the parenchymal lung structure in the normoxia- (Figures 4(a) and 4(c)) and hyperoxia- (Figures 4(b) and 4(d)) exposed pups at higher magnification and provides evidence of thickened alveolar septa and the presence of intra-alveolar inflammation with inflammatory cells after hyperoxia treatment.

(a)

(b)

(c)

(d)

(e)

Figure 2

Effect of hyperoxia treatment on lung volume and parenchymal volume. Total lung volume (a) was assessed with volume displacement. Changes in the parenchymal (b) and nonparenchymal (c) lung compartments were further quantified using stereology. Representative microimages of lung samples embedded in glycol methacrylate and stained with toluidine blue are shown in images (d) (normoxia) and (e) (hyperoxia). Lines indicate statistical differences between groups (), assessed with -test or Mann–Whitney rank sum test if normality distribution failed.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 3

Stereological data of alveolar and septal characteristics in the lung parenchyma. Stereology was used to quantify alveolar number (a), septal surface area (b), ductal airspace (c) and alveolar airspace (d), and mean alveolar number (e) in normoxic and hyperoxic rabbit lungs. The mean linear intercept length (f) was calculated for better comparison with previous studies. Stereological quantification of septal volume (g) and mean septal thickness (h) was also performed. Lines indicate statistical differences between groups (), assessed with -test or Mann–Whitney rank sum test if normality distribution failed.

(a)

(b)

(c)

(d)

3.3. Structural Changes in the Vasculature and Perivascular Area

Structural changes in the vasculature and the perivascular area included a decline in the vascular luminal volume (Figure 5(a)) and the mean vascular diameter (Figure 5(b)) in hyperoxia-treated lungs. The endothelial surface area of vessels with a diameter larger than 25 μm (Figure 5(c)) was not affected by hyperoxia exposure. However, the mean thickness of the vessel wall, the perivascular area, and the perivascular volume (Figures 5(d)–5(f)) significantly increased with hyperoxia exposure in preterm born rabbits. The apparent differences in perivascular and vascular wall thickness are shown in representative images of pulmonary arteries of comparable size of normoxia- (Figure 5(g)) and hyperoxia- (Figure 5(h)) treated rabbits.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 5

Alterations in vascular and perivascular structures. Vascular luminal volume (a), mean diameter (b), and surface area (c) as well as the mean thickness of the vascular wall (d) and the perivascular area (e) and thickness (f) were estimated using stereology. All data were collected from vessels with a diameter larger than 25 μm. Lines indicate statistical differences between groups (), assessed with -test or Mann–Whitney rank sum test if normality distribution failed. Representative images of vessels (here, pulmonary arteries) of normoxic (g) and hyperoxic (h) lungs show increased thickening of vascular wall (arrow) and perivascular area (arrowheads). b: bronchiole; v: vessel.

A summary of all acquired stereological data, including estimation of densities and total values, is provided in Table 1.

4. Discussion

The present study provides stereological reference data on an important animal model of BPD. In comparison with the normoxic controls, the lungs of preterm rabbit pups exposed to 95% hyperoxia were characterized by reduced alveolarization (smaller total number and surface area of alveoli, smaller total volume of parenchyma) and septal edema and/or fibroproliferation (thicker alveolar septa). In addition, the luminal volume of nonparenchymal blood vessels was reduced whereas the volume of their walls and perivascular connective tissue space was enlarged. Functionally, the lungs showed significant signs of parenchymal damage (increase in tissue resistance and elastance) in the absence of strong effects on larger airways (similar airway resistance). However, signs of pulmonary hypertension were present (decreased PAAT/PAET ratio).

Design-based stereology is the current gold standard of lung morphometry [21] because it does not depend on assumptions and provides robust data that allow for statistical testing. Design-based stereological methods have been extensively used in mouse models of BPD [22–24, 35–37] but not in the preterm rabbit model so far.

Most of the data reported in the current study are in accordance with previous morphometric studies in this animal model [25, 26] which reported similar functional data and morphometric parameters suggesting disrupted alveolarization. Although the previous studies were in line with the pathological concepts of BPD and allow comparisons between experimental groups of the same study, they are not suited to allow cross-study and cross-species comparisons. For example, Jiménez and colleagues reported a significantly higher mean linear intercept length (Lm) in the hyperoxic pups, a parameter that is mostly used as a measure of alveolar size [25]. In the present study, neither the Lm nor the mean volume of the alveoli was altered by hyperoxia. All measurements of alveolar size (Lm or mean alveolar volume) critically depend on the elastic properties of the lungs in relation to the applied fixation pressures [33]; thus, they strongly depend on the experimental setting. In addition, they are not suited to allow cross-species comparisons because of different alveolar size and elastic properties in differently sized animals. As the elastic properties may also differ among experimental groups, it is also possible that alveolar size measures can differ in the absence of alveolar number of surface area alterations.

Alveolar surface area is an important parameter to evaluate the functionally relevant diffusion area for gas exchange [38]. Its estimation, however, depends on the applied resolution and magnification, illustrated by the “Coast of Britain” phenomenon [39]. For example, surface area estimations of human lung alveoli range between 50-70 m² and 130-140 m² depending on whether the estimation based on light microscopy/microcomputed X-ray tomography [40, 41] or electron microscopy [38]. Tissue deformation/shrinkage during embedding can also significantly and differentially (between experimental groups) affect the results of surface area estimations [42, 43].

The most straightforward analysis of alveolarization is to estimate the total alveolar number. While originally based on geometric assumptions [44], an efficient and unbiased way of estimating alveolar number based on the disector principle [45] and the Euler-Poincaré characteristic [46] was established several years ago [32, 47]. The advantage over other morphometric methods is the possibility to compare BPD models or the influence of an experimental intervention in a certain model between laboratories and across animal species. In a rat hyperoxia model of BPD, two weeks of exposure to 60% O₂ resulted in a significant decrease in alveolar number whereas both total alveolar surface area and mean linear intercept length were not significantly different from those of the control group [48]. In a systematic analysis of various levels of hyperoxia (21%, 40%, 60%, and 85%) or differing periods of exposure to 85% oxygen, Nardiello et al. [14] provided stereological data for the mouse model of BPD. After exposure to 85% hyperoxia from postnatal days P4 to P7, the number of alveoli was reduced to 68% of controls, whereas alveolar surface area and Lm were not significantly different. Longer hyperoxia exposure (from P1 to P7), however, decreased alveolar number to 43% of controls which was accompanied by decreased alveolar surface area and increased Lm.

In the present study, the alveolar number was reduced to 68% of control lungs with no changes in Lm or mean alveolar volume, which corresponds well with the data for P4 to P7 hyperoxia exposure in the mouse study by Nardiello et al. [14]. In contrast to their study, however, the total surface area was significantly smaller in our hyperoxia group than in the control group. The results of the hyperoxic group in the present study — significantly reduced alveolar number and surface area but no increase in Lm or mean alveolar volume — are probably the result of a smaller degree of lung inflation during fixation. Less pronounced inflation leaves the alveolar number unaffected but leads to reduced surface area and airway volume estimations. The most probable reason for this was hyperoxia-induced inflammation and edema formation. Both factors are known to inactivate pulmonary surfactant and, thereby, reduce lung distensibility [49, 50]. This is well reflected by the higher intra-alveolar volumes of inflammatory cells and edematous fluid in combination with reduced static compliance and elevated resistance and elastance. Hyperoxia-induced lung injury in adult baboons led to similar morphological alterations and was relieved by exogenous surfactant treatment [51].

Besides the parenchymal changes, the present study also provides evidence for an involvement of larger sized arteries and veins in the hyperoxia model of BPD. Arteries and veins larger than 25 μm in diameter were grouped into a single compartment because previous studies using 3D reconstructions of the vasculature have shown that it may be impossible to distinguish morphologically between small arteries and veins in single 2D sections [52, 53]. Nevertheless, the data clearly show that the lumen of the blood vessels was narrower and the thickness of the wall as well as the perivascular space were increased in hyperoxic rabbits. These data complement the observations by Jiménez et al. [25] who reported a smaller peripheral vessel density and a thicker media layer of peripheral arteries. Our current data indicate that besides the peripheral arteries, remodeling of more centrally located blood vessels also takes place in the hyperoxia model and may additionally contribute to the pulmonary hypertension observed in this model [25, 54].

5. Conclusions

The present study provides robust stereological reference data for the preterm rabbit hyperoxia model of BPD. The data confirm and extend previous morphometric studies in this animal model and fit well with the pathological features of human BPD. In contrast to their normoxic counterparts, the lungs of preterm rabbits exposed to hyperoxia were characterized by reduced alveolarization, septal edema and/or fibroproliferation, thickened vessel walls, functional impairment of the parenchyma, and pulmonary hypertension. Together, these findings may serve as a reference to compare this model with BPD models in other species or future therapeutic interventions.

Data Availability

Data are available on request.

Disclosure

The company did not influence the study design, the interpretation of the results, and the revision of the manuscript.

Conflicts of Interest

Chiara Catozzi, Francesca Ricci, Costanza Casiraghi, and Fabrizio Salomone are Chiesi Farmaceutici employees. The rest of the authors declare that they do not have any conflict of interest.

Acknowledgments

We thank Susanne Kuhlmann for embedding and preparing of the histological section. The study was funded by Chiesi Farmaceutici S.p.A, Via Palermo 26/a, 43122, Parma, Italy. The work was further supported by the Deutsche Forschungsgemeinschaft (DFG MU 3118/8-1).

References

W. Northway, R. Rosan, and D. Porter, “Pulmonary disease following respirator therapy of hyaline-membrane disease,” The New England Journal of Medicine, vol. 276, no. 7, pp. 357–368, 1967.
View at: Publisher Site | Google Scholar
S. J. Cassady, A. Lasso-Pirot, and J. Deepak, “Phenotypes of bronchopulmonary dysplasia in adults,” Chest, vol. 158, no. 5, pp. 2074–2081, 2020.
View at: Publisher Site | Google Scholar
B. Thébaud, K. N. Goss, M. Laughon et al., “Bronchopulmonary dysplasia,” Nature Reviews Disease Primers, vol. 5, no. 1, 2019.
View at: Publisher Site | Google Scholar
J. J. Coalson, “Pathology of new bronchopulmonary dysplasia,” Seminars in Neonatology, vol. 8, no. 1, pp. 73–81, 2003.
View at: Publisher Site | Google Scholar
S. Niedermaier and A. Hilgendorff, “Bronchopulmonary dysplasia - an overview about pathophysiologic concepts,” Molecular and cellular pediatrics, vol. 2, no. 1, pp. 1–7, 2015.
View at: Publisher Site | Google Scholar
M. L. Choukroun, H. Feghali, S. Vautrat et al., “Pulmonary outcome and its correlates in school-aged children born with a gestational age ≤ 32 weeks,” Respiratory Medicine, vol. 107, no. 12, pp. 1966–1976, 2013.
View at: Publisher Site | Google Scholar
M. O'Reilly, F. Sozo, and R. Harding, “Impact of preterm birth and bronchopulmonary dysplasia on the developing lung: long-term consequences for respiratory health,” Clinical and Experimental Pharmacology & Physiology, vol. 40, no. 11, pp. 765–773, 2013.
View at: Publisher Site | Google Scholar
M. Vom Hove, F. Prenzel, H. H. Uhlig, and E. Robel-Tillig, “Pulmonary outcome in former preterm, very low birth weight children with bronchopulmonary dysplasia: a case-control follow-up at school age,” The Journal of Pediatrics, vol. 164, no. 1, pp. 40–45.e4, 2014.
View at: Publisher Site | Google Scholar
A. H. Jobe, “Animal models, learning lessons to prevent and treat neonatal chronic lung disease,” Frontiers in Medicine, vol. 2, pp. 1–13, 2015.
View at: Publisher Site | Google Scholar
R. E. Morty, “Using experimental models to identify pathogenic pathways and putative disease management targets in bronchopulmonary dysplasia,” Neonatology, vol. 117, no. 2, pp. 233–239, 2020.
View at: Publisher Site | Google Scholar
J. Berger and V. Bhandari, “Animal models of bronchopulmonary dysplasia. The term mouse models,” American Journal of Physiology-Lung Cellular and Molecular Physiology, vol. 307, no. 12, pp. L936–L947, 2014.
View at: Publisher Site | Google Scholar
M. O'Reilly and B. Thébaud, “Animal models of bronchopulmonary dysplasia. The term rat models,” American Journal of Physiology-Lung Cellular and Molecular Physiology, vol. 307, no. 12, pp. L948–L958, 2014.
View at: Publisher Site | Google Scholar
J. C. Schittny, “Development of the lung,” Cell and Tissue Research, vol. 367, no. 3, pp. 427–444, 2017.
View at: Publisher Site | Google Scholar
C. Nardiello, I. Mižíková, and R. E. Morty, “Looking ahead: where to next for animal models of bronchopulmonary dysplasia?” Cell and Tissue Research, vol. 367, no. 3, pp. 457–468, 2017.
View at: Publisher Site | Google Scholar
K. H. Albertine, “Progress in understanding the pathogenesis of BPD using the baboon and sheep models,” Seminars in Perinatology, vol. 37, no. 2, pp. 60–68, 2013.
View at: Publisher Site | Google Scholar
B. A. Yoder and J. J. Coalson, “Animal models of bronchopulmonary dysplasia. The preterm baboon models,” American Journal of Physiology-Lung Cellular and Molecular Physiology, vol. 307, no. 12, pp. L970–L977, 2014.
View at: Publisher Site | Google Scholar
C. T. D’Angio and R. M. Ryan, “Animal models of bronchopulmonary dysplasia. The preterm and term rabbit models,” American Journal of Physiology-Lung Cellular and Molecular Physiology, vol. 307, no. 12, pp. L959–L969, 2014.
View at: Publisher Site | Google Scholar
T. Salaets, A. Gie, B. Tack, J. Deprest, and J. Toelen, “Modelling bronchopulmonary dysplasia in animals: arguments for the preterm rabbit model,” Current Pharmaceutical Design, vol. 23, no. 38, pp. 5887–5901, 2018.
View at: Publisher Site | Google Scholar
C. Mühlfeld and M. Ochs, “Quantitative microscopy of the lung: a problem-based approach. Part 2: stereological parameters and study designs in various diseases of the respiratory tract,” American Journal of Physiology. Lung Cellular and Molecular Physiology, vol. 305, no. 3, pp. L205–L221, 2013.
View at: Publisher Site | Google Scholar
M. Ochs and C. Mühlfeld, “Quantitative microscopy of the lung: a problem-based approach. Part 1: basic principles of lung stereology,” American Journal of Physiology. Lung Cellular and Molecular Physiology, vol. 305, no. 1, pp. L15–L22, 2013.
View at: Publisher Site | Google Scholar
C. C. W. Hsia, D. M. Hyde, M. Ochs, and E. R. Weibel, “An official research policy statement of the American Thoracic Society/European Respiratory Society: standards for quantitative assessment of lung structure,” American Journal of Respiratory and Critical Care Medicine, vol. 181, no. 4, pp. 394–418, 2010.
View at: Publisher Site | Google Scholar
S. V. Appuhn, S. Siebert, D. Myti et al., “Capillary changes precede disordered alveolarization in a mouse model of bronchopulmonary dysplasia,” American Journal of Respiratory Cell and Molecular Biology, vol. 65, no. 1, pp. 81–91, 2021.
View at: Publisher Site | Google Scholar
D. Pérez-Bravo, D. Myti, I. Mižíková et al., “A comparison of airway pressures for inflation fixation of developing mouse lungs for stereological analyses,” Histochemistry and Cell Biology, vol. 155, no. 2, pp. 203–214, 2021.
View at: Publisher Site | Google Scholar
A. Pozarska, J. A. Rodríguez-Castillo, D. E. Surate Solaligue et al., “Stereological monitoring of mouse lung alveolarization from the early postnatal period to adulthood,” American Journal of Physiology-Lung Cellular and Molecular Physiology, vol. 312, no. 6, pp. L882–L895, 2017.
View at: Publisher Site | Google Scholar
J. Jiménez, J. Richter, T. Nagatomo et al., “Progressive vascular functional and structural damage in a bronchopulmonary dysplasia model in preterm rabbits exposed to hyperoxia,” International Journal of Molecular Sciences, vol. 17, no. 10, p. 1776, 2016.
View at: Publisher Site | Google Scholar
T. Salaets, M. Aertgeerts, A. Gie et al., “Preterm birth impairs postnatal lung development in the neonatal rabbit model,” Respiratory Research, vol. 21, no. 1, pp. 1–13, 2020.
View at: Publisher Site | Google Scholar
T. Salaets, A. Gie, J. Jimenez et al., “Local pulmonary drug delivery in the preterm rabbit: feasibility and efficacy of daily intratracheal injections,” American Journal of Physiology-Lung Cellular and Molecular Physiology, vol. 316, no. 4, pp. L589–L597, 2019.
View at: Publisher Site | Google Scholar
J. Richter, J. Toelen, J. Vanoirbeek et al., “Functional assessment of hyperoxia-induced lung injury after preterm birth in the rabbit,” American Journal of Physiology-Lung Cellular and Molecular Physiology, vol. 306, no. 3, pp. L277–L283, 2014.
View at: Publisher Site | Google Scholar
J. P. Schneider and M. Ochs, “Stereology of the lung,” in Methods in Cell Biology, P. M. Conn, Ed., pp. 257–294, Academic Press Inc., Oxford, 2013.
View at: Google Scholar
K. M. Kling, E. Lopez-Rodriguez, C. Pfarrer, C. Mühlfeld, and C. Brandenberger, “Aging exacerbates acute lung injury induced changes of the air-blood barrier, lung function and inflammation in the mouse,” American Journal of Physiology. Lung Cellular and Molecular Physiology, vol. 312, no. 1, pp. L1–12, 2017.
View at: Publisher Site | Google Scholar
H. Schulte, M. Mühlfeld, and C. Brandenberger, “Age-related structural and functional changes in the mouse lung,” Frontiers in Physiology, vol. 10, p. 1466, 2019.
View at: Publisher Site | Google Scholar
M. Ochs, J. R. Nyengaard, A. Jung et al., “The number of alveoli in the human lung,” American Journal of Respiratory and Critical Care Medicine, vol. 169, no. 1, pp. 120–124, 2004.
View at: Publisher Site | Google Scholar
L. Knudsen, E. R. Weibel, H. J. G. Gundersen, F. V. Weinstein, and M. Ochs, “Assessment of air space size characteristics by intercept (chord) measurement: an accurate and efficient stereological approach,” Journal of Applied Physiology, vol. 108, no. 2, pp. 412–421, 2010.
View at: Publisher Site | Google Scholar
C. Brandenberger, M. Ochs, and C. Mühlfeld, “Assessing particle and fiber toxicology in the respiratory system: the stereology toolbox,” Particle and Fibre Toxicology, vol. 12, no. 35, pp. 1–15, 2015.
View at: Publisher Site | Google Scholar
S. Feddersen, C. Nardiello, B. Selvakumar et al., “Impact of litter size on survival, growth and lung alveolarization of newborn mouse pups,” Annals of Anatomy, vol. 232, p. 151579, 2020.
View at: Publisher Site | Google Scholar
E. Lignelli, F. Palumbo, S. Bayindir et al., “The H₂S-generating enzyme 3-mercaptopyruvate sulfurtransferase regulates pulmonary vascular smooth muscle cell migration and proliferation but does not impact normal or aberrant lung development,” Nitric Oxide, vol. 107, pp. 31–45, 2021.
View at: Publisher Site | Google Scholar
J. Tiono, D. E. Surate Solaligue, I. Mižíková et al., “Mouse genetic background impacts susceptibility to hyperoxia-driven perturbations to lung maturation,” Pediatric Pulmonology, vol. 54, no. 7, pp. 1060–1077, 2019.
View at: Publisher Site | Google Scholar
P. Gehr, M. Bachofen, and E. R. Weibel, “The normal human lung: ultrastructure and morphometric estimation of diffusion capacity,” Respiration Physiology, vol. 32, no. 2, pp. 121–140, 1978.
View at: Publisher Site | Google Scholar
B. Mandelbrot, “How long is the coast of Britain? Statistical self-similarity and fractional dimension,” Science, vol. 156, no. 3775, pp. 636–638, 1967.
View at: Publisher Site | Google Scholar
D. M. Vasilescu, A. B. Phillion, D. Kinose et al., “Comprehensive stereological assessment of the human lung using multiresolution computed tomography,” Journal of Applied Physiology, vol. 128, no. 6, pp. 1604–1616, 2020.
View at: Publisher Site | Google Scholar
B. M. Wiebe, C. M. Burton, N. Milman, M. Iversen, and C. B. Andersen, “Morphometric examination of native lungs in human lung allograft recipients,” APMIS, vol. 114, no. 11, pp. 795–804, 2006.
View at: Publisher Site | Google Scholar
K. A. Dorph-Petersen, J. R. Nyengaard, and H. J. G. Gundersen, “Tissue shrinkage and unbiased stereological estimation of particle number and size,” Journal of Microscopy, vol. 204, no. 3, pp. 232–246, 2001.
View at: Publisher Site | Google Scholar
J. P. Schneider and M. Ochs, “Alterations of mouse lung tissue dimensions during processing for morphometry: a comparison of methods,” American Journal of Physiology-Lung Cellular and Molecular Physiology, vol. 306, no. 4, pp. L341–L350, 2014.
View at: Publisher Site | Google Scholar
E. R. Weibel and D. M. Gomez, “A principle for counting tissue structures on random sections,” Journal of Applied Physiology, vol. 17, no. 2, pp. 343–348, 1962.
View at: Publisher Site | Google Scholar
D. Sterio, “The unbiased estimation of number and sizes of arbitrary particles using the disector,” Journal of Microscopy, vol. 134, no. 2, pp. 127–136, 1984.
View at: Publisher Site | Google Scholar
H. J. G. Gundersen, R. W. Boyce, J. R. Nyengaard, and A. Odgaard, “The connEulor: unbiased estimation of connectivity using physical disectors under projection,” Bone, vol. 14, no. 3, pp. 217–222, 1993.
View at: Publisher Site | Google Scholar
D. M. Hyde, N. K. Tyler, L. F. Putney, P. Singh, and H. J. G. Gundersen, “Total number and mean size of alveoli in mammalian lung estimated using fractionator sampling and unbiased estimates of the Euler characteristic of alveolar openings,” The Anatomical Record. Part A, Discoveries in Molecular, Cellular, and Evolutionary Biology, vol. 277A, no. 1, pp. 216–226, 2004.
View at: Publisher Site | Google Scholar
A. Porzionato, D. Guidolin, V. Macchi et al., “Fractal analysis of alveolarization in hyperoxia-induced rat models of bronchopulmonary dysplasia,” American Journal of Physiology-Lung Cellular and Molecular Physiology, vol. 310, no. 7, pp. L680–L688, 2016.
View at: Publisher Site | Google Scholar
C. Mühlfeld, L. Becker, C. Bussinger et al., “Exogenous surfactant in ischemia/reperfusion: effects on endogenous surfactant pools,” The Journal of Heart and Lung Transplantation, vol. 29, no. 3, pp. 327–334, 2010.
View at: Publisher Site | Google Scholar
M. Ochs, I. Nenadic, A. Fehrenbach et al., “Ultrastructural alterations in intraalveolar surfactant subtypes after experimental ischemia and reperfusion,” American Journal of Respiratory and Critical Care Medicine, vol. 160, no. 2, pp. 718–724, 1999.
View at: Publisher Site | Google Scholar
C. A. Piantadosi, P. J. Fracica, F. G. Duhaylongsod et al., “Artificial surfactant attenuates hyperoxic lung injury in primates. II. Morphometric analysis,” Journal of Applied Physiology, vol. 78, no. 5, pp. 1823–1831, 1995.
View at: Publisher Site | Google Scholar
R. Grothausmann, L. Knudsen, M. Ochs, and C. Mühlfeld, “Digital 3D reconstructions using histological serial sections of lung tissue including the alveolar capillary network,” American Journal of Physiology-Lung Cellular and Molecular Physiology, vol. 312, no. 2, pp. L243–L257, 2017.
View at: Publisher Site | Google Scholar
R. Grothausmann, J. Labode, P. Hernandez-Cerdan et al., “Combination of μCT and light microscopy for generation-specific stereological analysis of pulmonary arterial branches: a proof-of-concept study,” Histochemistry and Cell Biology, vol. 155, no. 2, pp. 227–239, 2021.
View at: Publisher Site | Google Scholar
J. Jiménez, F. Lesage, J. Richter et al., “Upregulation of vascular endothelial growth factor in amniotic fluid stem cells enhances their potential to attenuate lung injury in a preterm rabbit model of bronchopulmonary dysplasia,” Neonatology, vol. 113, no. 3, pp. 275–285, 2018.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Christian Mühlfeld et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

440

Downloads

601

Citations