# Club Cell Loss as a Feature of Bronchiolization in ILD

Paul Reynaud, Engi Ahmed, Isabelle Serre, Lucie Knabe, Sébastien Bommart, Carey Suehs, Isabelle Vachier, Jean Philippe Berthet, Micaela Romagnoli, Charlotte Vernisse, et al.

► To cite this version:

Paul Reynaud, Engi Ahmed, Isabelle Serre, Lucie Knabe, Sébastien Bommart, et al. Club Cell Loss as a Feature of Bronchiolization in ILD. Frontiers in Immunology, 2021, 12, pp.630096. 10.3389/fimmu.2021.630096. [https://hal.science/hal-03170575v2](https://hal.science/hal-03170575v2)

HAL Id: hal-03170575

https://hal.science/hal-03170575v2

Submitted on 18 May 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

Check for
updates

# Club Cell Loss as a Feature of Bronchiolization in ILD

Paul Reynaud<fnref n="1" />, Engi Ahmed<fnref n="1" />,², Isabelle Serre³, Lucie Knabe<fnref n="1" />, Sébastien Bommart⁴,
Carey Suehs<fnref n="1" />, Isabelle Vachier<fnref n="1" />, Jean Philippe Berthet⁵, Micaela Romagnoli⁶,
Charlotte Vernisse<fnref n="1" />,², Jean Pierre Mallet<fnref n="1" />, Anne Sophie Gamez<fnref n="1" /> and Arnaud Bourdin<fnref n="1" />,²*

## OPEN ACCESS

**Edited by:**
Antoine Froidure,
Cliniques Universitaires
Saint-Luc, Belgium

**Reviewed by:**
Arnaud Mailleux,
Institut National de la Santé et de la
Recherche Médicale
(INSERM), France
Charles Pilette,
Catholic University of
Louvain, Belgium

*Correspondence:
Arnaud Bourdin
a-bourdin@chu-montpellier.fr

**Specialty section:**
This article was submitted to
Mucosal Immunity,
a section of the journal
Frontiers in Immunology

**Received:** 16 November 2020

**Accepted:** 25 January 2021

**Published:** 26 February 2021

**Citation:**
Reynaud P, Ahmed E, Serre I, Knabe L, Bommart S, Suehs C, Vachier I, Berthet JP, Romagnoli M, Vernisse C, Mallet JP, Gamez AS and Bourdin A (2021) Club Cell Loss as a Feature of Bronchiolization in ILD. *Front. Immunol.* 12:630096.
doi: 10.3389/fimmu.2021.630096

**Background:** Distal airway metaplasia may precede honeycombing in progressive fibrosing interstitial lung disease (ILD). The SCGB1A1⁺ bronchiolar-specific club cell may play a role in this aberrant regenerative process.

**Objective:** To assess the presence of club cells in the small airways of patients suffering from ILD.

**Methods:** Small airways (internal diameter <2 mm) in lung samples [surgical lung biopsy (SLB) and/or transbronchial lung cryobiopsy (TBLC)] from 14 patients suffering from ILD and 10 controls were morphologically assessed and stained for SCGB1A1. SCGB1A1 was weighted by epithelial height as a marker of airway generation (SCGB1A1/EH). Correlations between clinical, functional, and high-resolution CT (HRCT) prognostic factors and histomorphometry were assessed.

**Results:** Small airways from samples with ILD patterns were significantly less dense in terms of SCGB1A1⁺ cells [0.064 (0.020–0.172)] as compared to controls’ sample’s small airways [0.393 (0.082–0.698), *p* < 0.0001]. Usual interstitial pneumonia (UIP) patterns most frequently contained small airways with limited or absent SCGB1A1 expression (SCGB1A1/EH <0.025): UIP (18/33; 55%) as compared with non-UIP patterns (4/31; 13%) or controls (0/29; 0%): *p* < 0.0001. In addition, correlations with HRCT indicated a significant negative relationship between SCGB1A1 and bronchiectasis as a feature of bronchiolization (Rho −0.63, *p* < 0.001) and a positive relationship with both forced vital capacity (FVC) and Hounsfield unit (HU)-distribution pattern in kurtosis (Rho 0.38 and 0.50, respectively, both *p* < 0.001) as markers of fibrotic changes.

**Conclusion:** Compared with controls, the small airways of patients with ILD more often lack SCGB1A1, especially so in UIP. Low densities of SCGB1A1-marked cells correlate with bronchiectasis and fibrotic changes. Further research investigating SCGB1A1 staining as a pathological feature of the bronchiolization process is merited.

**Keywords:** club (clara) cell, idiopathic pulmonary fibrosis, metaplasia, SCGB1A1, bronchiolization, interstitial lung disease

<footnote n="1">¹ Department of Respiratory Diseases, Univ Montpellier, CHU Montpellier, Montpellier, France, ² PhyMedExp, Univ Montpellier, CNRS, INSERM, CHU Montpellier, Montpellier, France, ³ Department of Pathology, Univ Montpellier, CHU Montpellier, Montpellier, France, ⁴ Department of Radiology, Univ Montpellier, CHU Montpellier, Montpellier, France, ⁵ Pulmonology Unit, Ospedale Ca' Foncello, AULSS2 Marca Trevigiana, Trevise, Italy, ⁶ Department of Cardiac, Thoracic and Vascular Surgery, Univ Nice, CHU Nice, Nice, France</footnote>
# INTRODUCTION

Interstitial lung disease (ILD) are devastating affections nearly constantly leading to respiratory failure and death. The early phases of the disease are hallmarked by progressive fibrosis in the peripheral zones of the lung, which is associated with lung function decline (1). High-resolution CT (HRCT) has become the gold standard for diagnosing such fibrotic changes and characterizing their progression (2). Notably, the presence of traction bronchiectasis (the enlargement of small airways by traction) is now considered a red flag for subpleural fibrotic changes. Interestingly, histomorphological assessments of these fibrotic areas consistently report bronchiolar metaplasia of unclear origin (3). The latter "bronchiolization" may represent aberrant regeneration of the lung in a profibrotic environment, where distal airway precursors never reach their potential as functional alveoli (4). The result is the presence of proximal features in the distal compartment currently detectable only *via* histopathology (5).

Club cells are thought to be highly expressed at the small airway level (6), where they are supposed to play a protective role against noxious inhaled particles. Among the latter, cigarette smoke and occupational dust have also been implicated in the pathogenesis of idiopathic pulmonary fibrosis (IPF) (7, 8). Interestingly, club cells were also attributed an accessory stem cell role at the small airway level, but this role has mostly been evidenced in animal models at the bronchoalveolar duct junction, where they are preferentially located (9). Due to recent advances in single-cell biology, it is now clear that these cells are heterogeneous and not all are playing the same role. Currently, they are considered as differentiation of basal cells, and the transition from basal to club cells has been ontogenically evidenced. Additionally, club cells are situated at the crossroads of the epithelial phenotype. Their natural trend is to progressively acquire features of goblet cells, but some can also move toward ciliogenesis. Whether or not appropriate transition pathways are failing during the bronchiolisation process remain unclear.

New opportunities for lung tissue sampling are provided by the arrival of transbronchial lung cryobiopsies (TBLC). Beyond clinical and technical considerations, it seems that TBLC, as

compared to traditional surgical lung biopsies (SLB), samples airway-centered, more proximal tissue, and this is likely to change not only the pathological patterns and eventually the final pathological diagnosis but also the representation of the different compartments. In particular, airways contained in TBLC samples are more likely to represent earlier generations than those taken in the more peripheral SLB samples (10). Intuitively, the small airways present in TBLC samples might then be less frequently affected by bronchiolization. Given the relationship between this process and the fibrotic changes affecting the lung in progressing ILD, these pathological changes may have some diagnostic and/or prognostic value identifiable by clinical correlations in a cross-sectional study.

The objective of this study was to assess the presence of club cells in the small airways of patients suffering from ILD by combining pathology, physiology, and HRCT data.

We hypothesized that "protective" club cells characterized by their specific protein SCGB1A1 might be deficient in fibrotic diseases and particularly at the site where bronchiolar metaplasia occurs.

# MATERIALS AND METHODS

CryoPID [NCT02763540] was a prospective trial evaluating the diagnostic concordance between TBLC and SLB, which were performed sequentially in the same patients between January 2016 and March 2018 at the Montpellier University Hospital (11). The study included patients with ILD requiring an SLB as per a multidisciplinary approach decision, given an undefined HRCT pattern. For the current histology assessment, 14 CryoPID patients were compared with 10 control subjects randomly selected from the bio-bank of the Montpellier University Hospital. Pathological patterns reported from SLB and TBLC together by a pathologist involved in the present study were used. Controls had normal lung function, were free of chronic respiratory disease, and matched, as well as a possibility for smoking, in gender and age with the ILD population. They had undergone a lobar resection for a colic metastasis in the last 3 years. Lung parenchyma at a large distance from the metastasis was analyzed.

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | After SCGB1A1 staining, small airway epithelia were isolated *via* thresholding, perimeters, thickness/height, and by the area quantified.</figcaption></figure>
Briefly, lung samples were fixed in 4% formalin at pH 7.4 and embedded in paraffin. Three-millimeter-thick sections were stained by H&E. Lung samples were anonymized and blindly assessed by an expert pathologist (TVC). Serial sections were used for SCGB1A1 immunostaining (dilution 1/10,000; Biovendor, Czech Republic). For some slides, a Bleu Alcian staining could also be performed.

Small airways (identified by an internal diameter of 2 mm or less) were assessed in all samples. Morphometry was assessed using the Cell P software (Olympus, Tokyo, Japan). Bronchiolization was referred to small airways with evidence of proximal features, including the presence of ciliated and/or goblet cells, as per the previously published standards (12, 13). Epithelial area, epithelial height (EH), and the percentage of SCGB1A1⁺ immune-stained area were determined for each small airway. The SCGB1A1⁺/EH ratio was computed to adjust for variations known to affect club cell density through airway generation (Figure 1). The small airway density for each sample was based on a semi-quantitative scale ranging from 0 (absence of small airways) to 5 (highest small airway density). The absence or near-absence of club cells was defined by a SCGB1A1/EH ratio <0.025 by taking into account the background noise.

Prognostic factors were recorded and included clinical (age, gender, smoking history), biological/laboratory [arterial oxygen pressure (PaO₂), bronchoalveolar lavage cellularity (BAL)], and functional [forced vital capacity (FVC), diffusing capacity of the lungs for carbon monoxide (DLCO)] criteria.

A visual tomodensitometric analysis scored fibrotic changes (honeycombing and reticulations), intensity, and the spread of ground-glass opacity (GGO) (14), as well as the extent of bronchiectasis [ranking the severity of bronchiectasis by a lobe on a scale ranging from 0 to 3 (15) by acquiring consensus among two experienced thoracic senior radiologists (SB and GD)].

Automatic quantification via the Thoracic VCAR software (GE HealthCare, USA) was also performed, classifying pixels from the whole lung (WL) volume into four groups as proposed by Shin et al. (16) (< -950 HU: "emphysema," -950 to -700 HU "functional parenchyma," -700 to -500 HU "interstitial disease," > -500 HU "extra parenchymal"). Three prognostic variables were then derived: total lung capacity (TLC), percentage of normal lung (%NL), and mean lung attenuation (MLA) (17). Pulmonary attenuation distribution histograms were obtained for assessing skewness and kurtosis, previously reported as having the prognostic value (18).

## Statistical Analysis

Parametric and non-parametric tests were used to describe and compare normally and non-normally distributed data as appropriate. Spearman's correlation coefficients were computed for testing variable relationships.

## RESULTS

### Population Description

Fourteen patients with ILD and 10 controls were studied (Table 1). As expected, patients with ILD had worse pulmonary function tests and higher HRCT semi-quantitative scores. Of

<table><caption>TABLE 1 | Population characteristics.</caption><thead><tr><th></th><th>ILD [Median (IQR25-75)]</th><th>Controls [Median (IQR25-75)]</th><th>p-value</th></tr></thead><tbody><tr><td>N</td><td>14</td><td>10</td><td></td></tr><tr><td>Age (years)</td><td>63 [43.5-65]</td><td>64 [57.5-67.5]</td><td>0.21</td></tr><tr><td>Gender (W/M)</td><td>10 vs. 04</td><td>3 vs. 7</td><td>0.0987</td></tr><tr><td>Weight (kg)</td><td>65.5 [60.65-81.75]</td><td>65 [62.25-83.75]</td><td>0.49</td></tr><tr><td>Height (cm)</td><td>160 [158.5-170]</td><td>166.5 [160-170]</td><td>0.44</td></tr><tr><td>Smoking history (p.y)</td><td>2.5 [0-28.75]</td><td>19 [1.25-30]</td><td>0.77</td></tr><tr><td>Smoking status (current-former-never)</td><td>1/5/8</td><td>1/5/4</td><td></td></tr><tr><td>Positive autoantibodies</td><td>0</td><td>0</td><td></td></tr><tr><td>PaO<sub>2</sub> (mmHg)</td><td>84.2 [75.75-88.28]</td><td>88.65 [81.29-93.5]</td><td>0.15</td></tr><tr><td>FVC (l)</td><td>2.3 [1.74-3.31]</td><td>3.43 [2.75-4.07]</td><td>0.027*</td></tr><tr><td>FVC (% pred values)</td><td>84.5 [68.5-2.77]</td><td>100 [92.25-107.5]</td><td>0.004*</td></tr><tr><td>FEV1 (l)</td><td>2.1 [1.59-2.77]</td><td>2.39 [2.02-3.39]</td><td>0.077</td></tr><tr><td>FEV1 (% pred values)</td><td>84 [72.75-90]</td><td>95.5 [86-110.8]</td><td>0.010*</td></tr><tr><td>TLC (l)</td><td>3.69 [3.55-4.49]</td><td>6.23 [6.16-6.59]</td><td>0.001*</td></tr><tr><td>TLC (% pred values)</td><td>77.5 [66.25-84.75]</td><td>103 [101-106]</td><td>0.00005*</td></tr><tr><td>RV (l)</td><td>1.57 [1.31-1.79]</td><td>2.64 [2.58-3.12]</td><td>0.0002*</td></tr><tr><td>RV (% pred values)</td><td>89 [65-102]</td><td>116 [100-121]</td><td>0.0049*</td></tr><tr><td>TLCO (% pred values)</td><td>56 [40-60]</td><td>76 [70-80]</td><td>0.0097*</td></tr><tr><td>KCO</td><td>79 [73-96]</td><td>85 [78-103]</td><td>0.58</td></tr><tr><td>Kurtosis</td><td>6.33 [4.32-7.48]</td><td>14.77 [13.78-17.16]</td><td>0.0017*</td></tr><tr><td>Skewness</td><td>2.29 [2.02-2.53]</td><td>3.59 [3.39-3.85]</td><td>0.00014*</td></tr><tr><td>Whole lung (WL)</td><td>3.44 [2.79-4.12]</td><td>5.61 [3.54-6.17]</td><td>0.1095</td></tr><tr><td>% Normal lung (NL)</td><td>79.98 [71.14-82.04]</td><td>88.72 [88.2-90.25]</td><td>0.00045</td></tr><tr><td>Mean lung attenuation (Hounsfield)</td><td>-764 [-796 to -720]</td><td>-838 [-859 to -819]</td><td>0.0012*</td></tr><tr><td>Fibrotic score</td><td>15 [10-22.5]</td><td>0 [0-0]</td><td>0.00006*</td></tr><tr><td>Ground glass opacity</td><td>5 [5-5]</td><td>0 [0-0]</td><td>0.059</td></tr><tr><td>Bronchiectasis score</td><td>6 [5.25-8]</td><td>0 [0-0]</td><td>0.00003*</td></tr></tbody></table>

%NL, percentage of normal lung; BMI, body mass index; CT, computed tomography; DLCO, diffusing capacity of the lungs for carbon monoxide; FEV1, forced expiratory volume in 1 s; FVC, forced vital capacity; ILD, interstitial lung disease; IQR, interquartile range; KCO, transfer coefficient of the lung for carbon monoxide; MLA, mean lung attenuation; N, number; PaO₂, arterial oxygen pressure; RV, residual volume; SD, standard deviation; TLC, total lung capacity.

note, as these patients were deemed eligible for SLB, most had a mild functional impairment, suggesting that many were at an early stage of the disease. As previously reported (11), discrepancies in the pathological diagnosis between SLB and TBLC existed, potentially leading to different final diagnoses during a second multidisciplinary assessment (Table 2). For subsequent classification as usual interstitial pneumonia (UIP) and non-UIP, only samples with a clear and confirmed (high level of confidence) UIP pattern were retained.

Among patients with ILD, the mean SLB sample size was 49.5 mm [41.5-60] compared to 7 mm (6-8) (*p* < 0.00001) for TBLC. A significant difference in the semi-quantitative scoring

<footnote>*p < 0.05.</footnote>
of small airway density between SLB and TBLC sample was found (Table 3). No difference was found between TBLC and SLB

<table><caption>**TABLE 2** | Blinded histology and multidisciplinary diagnostic results for each patient.</caption><thead><tr><th>Patient</th><th>Histology diagnosis based on TBLC</th><th>Histology diagnosis based on SLB</th><th>Diagnosis based on MDA 1 (before histology)</th><th>Diagnosis based on MDA 2 (after histology)</th></tr></thead><tbody><tr><td>1</td><td>UIP</td><td>UIP</td><td>Possible UIP</td><td>CHP</td></tr><tr><td>2</td><td>Non-diagnostic</td><td>PLCH</td><td>DIP</td><td>PLCH</td></tr><tr><td>3</td><td>UIP</td><td>CHP</td><td>NSIP</td><td>IPF</td></tr><tr><td>4</td><td>UIP</td><td>NSIP</td><td>No classification</td><td>NSIP</td></tr><tr><td>5</td><td>RB-ILD</td><td>RB-ILD</td><td>RB-ILD</td><td>RB-ILD</td></tr><tr><td>6</td><td>Non-diagnostic</td><td>UIP</td><td>NSIP</td><td>IPF</td></tr><tr><td>7</td><td>PLCH</td><td>UIP</td><td>Possible UIP</td><td>IPF</td></tr><tr><td>8</td><td>NSIP</td><td>NSIP</td><td>Fibrotic NSIP</td><td>NOS</td></tr><tr><td>9</td><td>UIP</td><td>UIP</td><td>NSIP</td><td>IPF</td></tr><tr><td>10</td><td>Non-diagnostic</td><td>ALI</td><td>NSIP</td><td>NOS</td></tr><tr><td>11</td><td>NSIP</td><td>LP</td><td>Sarcoidosis</td><td>CVID</td></tr><tr><td>12</td><td>UIP</td><td>NSIP</td><td>NSIP</td><td>NSIP</td></tr><tr><td>13</td><td>Non-diagnostic</td><td>UIP</td><td>CHP</td><td>IPF</td></tr><tr><td>14</td><td>UIP</td><td>UIP</td><td>Possible UIP</td><td>IPF</td></tr></tbody></table>

CHP, chronic hypersensitivity pneumonitis; CVID, common variable immune deficiency; DIP, desquematative interstitial pneumonia; HP, hypersensitivity pneumonitis; IPF, idiopathic pulmonary fibrosis; LP, lymphoid process; MDA1, a first multidisciplinary assessment; MDA2, the final multidisciplinary assessment; NOS, not otherwise specified; NSIP, non-specific interstitial pneumonia; OP, organizing pneumonia; PLCH, pulmonary Langerhans cell histiocytosis; RB-ILD, respiratory bronchiolitis-associated interstitial lung disease; SLB, surgical lung biopsy; TBLC, transbronchial lung cryobiopsy; UIP, usual interstitial pneumonia.

<table><caption>**TABLE 3** | Samples characteristics.</caption><thead><tr><th></th><th>SLB (controls)</th><th>SLB (ILD)</th><th>TBLC (ILD)</th><th>p-value</th></tr></thead><tbody><tr><td>N</td><td>10</td><td>14</td><td>14</td><td></td></tr><tr><td>Number of small airways assessed</td><td>32</td><td>35</td><td>31</td><td></td></tr><tr><td>Small airway density [semi quantitative score (0 = minimal; 5 = maximal)]</td><td>2.1 ± 0.1</td><td>2.4 ± 0.6</td><td>0.4 ± 0.3</td><td>0.04</td></tr><tr><td>Epithelial area (mean ± SD per bronchiole. µm²)</td><td>20,818 ± 9,040</td><td>26,738 ± 11,733</td><td>19,638 ± 11,216</td><td>ns</td></tr><tr><td>Epithelial height (mean ± SD per bronchiole. µm)</td><td>37.8 ± 7.8</td><td>39.5 ± 11.9</td><td>35.8 ± 11.4</td><td>ns</td></tr><tr><td>SCGB1A1⁺ stained area</td><td>19.7 (4.2–27.7)</td><td>3.6 (1.6–7.4)</td><td>0.9 (0.3–2.9)</td><td>&lt;0.0001</td></tr><tr><td>SCGB1A1⁺/EH</td><td>0.5 (0.1–0.7)</td><td>0.1 (0.0–0.2)</td><td>0.0 (0.0–0.1)</td><td>0.02</td></tr><tr><td>% airway without club cells<sup>§</sup></td><td>0 (0/32)</td><td>14 (5/35)</td><td>55 (17/31)</td><td>0.01</td></tr></tbody></table>

regarding club cell representation in the small airways. Fewer airways positively stained for SCGB1A1 in ILD, in particular, for those with a pathological UIP pattern, and in general, this staining was less intense (Figure 2A). This decrease persisted after adjusting on epithelial height as a surrogate marker of the generation of the airway assessed. Hence, SCGB1A1/EH was statistically lower in the small airways of samples from patients with ILD vs. controls (Figure 2B; the Kruskal-Wallis test, $p = 7.598e-07$).

Moreover, sparsely stained airways (defined as SCGB1A1/EH <0.025) were significantly more represented in the UIP pathological patterns [UIP 18/33 vs. non-UIP 4/31 ($p = 0.0027$)] (Table 2 and Figure 3).

The tested prognostic correlations are provided in Figure 4. Interestingly, statistically significant non-parametric correlations were found between SCGB1A1/EH and two surrogate markers of lung fibrosis. First, a correlation was found with the kurtosis value of the histogram distribution of lung attenuation assessed at HRCT (Rho = 0.50, $p < 0.001$). The less the lung is involved by fibrosis, the greater the club cell density found in the small airways, as the kurtosis value decreases with the extent of fibrosis in ILD (18). Second, SCGB1A1/EH also positively correlated with FVC (Rho = 0.38, $p < 0.001$), as expected.

The negative correlation found between SCGB1A1/EH and the bronchiectasis scores (Rho = -0.63, $p < 0.001$) is also represented in Figure 4. The latter indicates that fewer club cells were present in the small airways of patients with more intense bronchiectasis.

## DISCUSSION

With the aim of understanding the mechanisms of bronchiolization, we found that club cells are absent or nearly absent from these aberrant small airways in the most subpleural areas of fibrotic lungs. This was particularly true in ILD where a UIP pattern could be identified. Bronchiolized areas seem to fit with small airways saliently featuring proximal differentiation patterns very close to terminal airspaces. As expected, both TBLC and SLB could identify these patterns, even though the density of small airways in TBLC is largely lower. Interestingly, the correlation between these bronchiolization patterns (highlighted by the absence of club cells) and clinical prognostic markers (such as HRCT kurtosis and bronchiectasis), as well as physiology (FVC), was demonstrated.

Most of our results were driven by a particularly high number of unstained or sparsely stained airways in lung samples from patients with a UIP pattern. Most of these pauci-SCGB1A1 small airways displayed features of bronchiolar metaplasia. Interestingly, they were referred to as “bronchiolizing areas” in the routine pathology report.

Fukumoto et al. hypothesized that this phenomenon may indicate an impaired differentiation with aberrant proliferation (19). The latter hypothesis is supported by the morphologic appearance of these areas where proximal features are found where they should not be.

The underexpression of the club cell-specific protein SCGB1A1 may be directly linked to the bronchiolization process, as already supposed by Jensen-Taubman (20). A blockade in the

<footnote>§ Defined by SCGB1A1⁺/EH ratio < 0.025.</footnote>
<figure><img src="image_4.png" /></figure>

**FIGURE 2** | (A) Small airways are immune-stained for SCGB1A1 in controls, UIP, and non-UIP patterns. Upper left: (A) SLB from a control (no TBLC were available in this group). Middle row: (B) Very limited staining in some (left lower part) but not all small airways sampled by SLB; Inserts are: Serial cut with a Blue Alcian staining suggestive of the presence of mucin-costaining and a higher magnification suggestive of goblet cell morphology for SCGB1A1 positively stained epithelial cells. Ectasic small airway from a TBLC UIP pattern with bronchiolar metaplasia nearly free of staining despite abundant mucostasis. Right: SCGB1A1 stained small airways obtained by SLB (D) and by TBLC (E) in non-UIP patterns. (B) Box plots demonstrating group differences (Control vs. ILD for SLB and TBLC) for SCGB1A1 staining variables. TBLC, transbronchial lung cryobiopsies; SLB, surgical lung biopsy; UIP, usual interstitial pneumonia; ILD, interstitial lung disease.

<figure><img src="image_5.png" /></figure>

**FIGURE 3** | The proportion of samples with SCGB1A1/EH <0.025 was used to define the absence of positive staining [UIP 18/33 vs. non-UIP 4/31 (p = 0.0027)].
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | Scatterplots demonstrating the relationships between selected pulmonary function variables (bronchiectasis score (A), % predicted forced vital capacity (FVC, % predicted) (B), and kurtosis of lung attenuation scores (C)] and % epithelial area of small airways that stains SCGB1B1⁺ weighted by epithelial height (EA<sup>SCGB1B1+</sup>/EH). Smoothing is performed via a local loess estimator (in black). Spearman's correlation statistics (Rho and FDR-corrected p-value) are provided.</figcaption></figure>

notch pathway might be involved in this ontogenically impaired process, as the notch is critically involved in basal to club cell transitioning (21). Wnt/β-catenin pathways were also shown to be involved.

This mechanism may help accelerate the paradigm shift that is already in motion, where instead of envisioning bronchiolization and bronchiectasis as a consequence of fibrosis, they are now understood to be an active part of the fibrotic process [as again recently suggested (4)].

Senescence may play a key role in the fibrotic process characteristic of UIP (22) possibly *via* this bronchiolization. Aberrantly proliferating cells can lead to both replicative and premature senescence (23), possibly leading to these poorly differentiated airways. The secretion of senescence-associated secretory protein (SASP), mainly interleukin-8, has already been shown to be involved (24).

Our study has several limits. Few patients were included due to the particular procedures used in the CryoPID study. Our patient group was so heterogeneous that artificially gathering patients from a “non-UIP” group was required. Moreover, staining for MUC5B, a validated bronchiolization marker, would have been of great interest to increase confidence in the mechanism and further characterize the phenotype of the remaining club cells. As shown by Best et al. (25), a higher proportion of SCGB1A1⁺/MUC5B⁺ cells is observed in IPF, highly expressing genes related to mucins and chemoattractant cytokines for immune cells. Though we missed this opportunity, we nonetheless had attempted the Bleu Alcian staining in a few slides when some unused material was available; the latter suggested convergent findings and merits further investigation.

Of note, we decided to focus our analysis at the sample level (as opposed to the patient level) to take into account the regional heterogeneity of the fibrotic process where zones of normal appearance are classically neighboring the involved tissue (6). This was also justified by the aim of comparing

TBLC and SLB in this regard. Correlations were tested at the individual level by expressing SCGB1A1 as a mean and our findings remained consistent.

There are a few steps left to definitively demonstrate that bronchiolization precedes the fibrotic process. Recently, by matching micro-CT imaging with histology, Verleden et al. showed that small airway involvement is the most striking finding in the early IPF stages, suggesting that the disease starts at this level (4).

Our study also supports that the bronchiolization process is key in the course of the disease as we could find correlations with two different HRCT-derived parameters. Kurtosis refers to the tailedness of the distribution histogram of lung attenuation, and this value was suggested to be a surrogate of interstitial changes (25). A significant negative relationship between the expression of the club cell marker SCGB1A1 and the bronchiectasis score was also found; an intuitively-expected finding enhancing the concept that bronchiectasis found at HRCT is related to the bronchiolization process hallmarked by the absence of club cells. Convincingly, a significant correlation was also found with FVC, a physiological marker more than relevant in interstitial lung diseases.

In summary, we demonstrated that bronchiolization, a typical characteristic of ILD and of UIP in particular, is also hallmarked by the absence of club cells. This suggests that an impaired transition process affected epithelial cell fate.

Pauci-SCGB1A1 staining status merits further study as a potential marker of bronchiolization in ILD.

## DATA AVAILABILITY STATEMENT

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
## ETHICS STATEMENT

The studies involving human participants were reviewed and approved by CPP Ile de France. The patients/participants provided their written informed consent to participate in this study.

## AUTHOR CONTRIBUTIONS

AB elaborated the study, enrolled patients, performed biopsies, participated to the interpretation of the data, and the writing of the manuscript. PR participated to the

study, enrolled patients, performed biopsies, acquisition of data, immunostaining, morphometry, HRCT analysis, and the writing of the manuscript. EA participated to the study, enrolled patients, performed biopsies, acquisition of data, and the writing of the manuscript. IS participated to the study, assessed biopsies, and contributed to the immunostaining. LK, CV, CS, and IV participated to the study, interpretation of the data, and the writing of the manuscript. SB participated to HRCT analysis. MR, JB, JM, and AG participated to the study, enrolled patients, and performed biopsies. All authors contributed to the article and approved the submitted version.

## REFERENCES

1. Olson AL, Gifford AH, Inase N, Fernández Pérez ER, Suda T. The epidemiology of idiopathic pulmonary fibrosis and interstitial lung diseases at risk of a progressive-fibrosing phenotype. *Eur Respir Rev.* (2018) 27:180077. doi: 10.1183/16000617.0077-2018

2. Sverzellati N, Brillet P-Y. When Deep Blue first defeated Kasparov: is a machine stronger than a radiologist at predicting prognosis in idiopathic pulmonary fibrosis? *Euro Respir J.* (2017) 49:1602144. doi: 10.1183/13993003.02144-2016

3. Visscher DW, Myers JL. Histologic spectrum of idiopathic interstitial pneumonias. *Proc Am Thorac Soc.* (2006) 3:322-9. doi: 10.1513/pats.200602-019TK

4. Verleden SE, Tanabe N, McDonough JE, Vasilescu DM, Xu F, Wuys WA, et al. Small airways pathology in idiopathic pulmonary fibrosis: a retrospective cohort study. *Lancet Respir Med.* (2020) 8:573-84. doi: 10.1016/S2213-2600(19)30356-X

5. Plantier L, Crestani B, Wert SE, Dehoux M, Zweytick B, Guenther A, et al. Ectopic respiratory epithelial cell differentiation in bronchiolised distal airspaces in idiopathic pulmonary fibrosis. *Thorax.* (2011) 66:651-7. doi: 10.1136/thx.2010.151555

6. Boers JE, Ambergen AW, Thunnissen FB. Number and proliferation of clara cells in normal human airway epithelium. *Am J Respir Crit Care Med.* (1999) 159:1585-91. doi: 10.1164/ajrccm.159.5.9806044

7. Akram KM, Lomas NJ, Spiteri MA, Forsyth NR. Club cells inhibit alveolar epithelial wound repair via TRAIL-dependent apoptosis. *Euro Respir J.* (2013) 41:683-94. doi: 10.1183/09031936.00213411

8. Singh G, Katyal SL. Clara cells and clara cell 10 kD protein (CC10). *Am J Respir Cell Mol Biol.* (1997) 17:141-3. doi: 10.1165/ajrcmb.17.2.f138

9. Rawlins EL, Okubo T, Xue Y, Brass DM, Auten RL, Hasegawa H, et al. The role of Scgb1a1+ Clara cells in the long-term maintenance and repair of lung airway, but not alveolar, epithelium. *Cell Stem Cell.* (2009) 4:525-34. doi: 10.1016/j stems.2009.04.002

10. Colby TV, Tomassetti S, Cavazza A, Dubini A, Poletti V. Transbronchial cryobiopsy in diffuse lung disease: update for the pathologist. *Arch Pathol Lab Med.* (2017) 141:891-900. doi: 10.5858/arpa.2016-0233-RA

11. Romagnoli M, Colby TV, Berthet J.-P, Gamez AS, Mallet J.-P, Serre I, et al. Poor concordance between sequential transbronchial lung cryobiopsy and surgical lung biopsy in the diagnosis of diffuse interstitial lung diseases. *Am J Respir Crit Care Med.* (2019) 199:1249-56. doi: 10.1164/rccm.201810-1947OC

12. Nettesheim P, Szakal AK. Morphogenesis of alveolar bronchiolization. *Lab Invest.* (1972) 26:210-9.

13. Allen TC. Pathology of small airways disease. *Arch Pathol Lab Med.* (2010) 134:702-18. doi: 10.5858/134.5.702

14. Shin KM, Lee KS, Chung MP, Han J, Bae YA, Kim TS, et al. Prognostic determinants among clinical, thin-section CT, and histopathologic findings for fibrotic idiopathic interstitial pneumonias: tertiary hospital study. *Radiology.* (2008) 249:328-37. doi: 10.1148/radiol.2483071378

15. Edey AJ, Devaraj AA, Barker RP, Nicholson AG, Wells AU, Hansell DM. Fibrotic idiopathic interstitial pneumonias: HRCT findings that predict mortality. *Eur Radiol.* (2011) 21:1586-93. doi: 10.1007/s00330-011-2098-2

16. Shin KE, Chung MJ, Jung MP, Choe BK, Lee KS. Quantitative computed tomographic indexes in diffuse interstitial lung disease: correlation with physiologic tests and computed tomography visual scores. *J Comput Assist Tomogr.* (2011) 35:266-71. doi: 10.1097/RCT.0b013e31820ccf18

17. Ohkubo H, Kanemitsu Y, Uemura T, Takakuwa O, Takemura M, Maeno K, et al. Normal lung quantification in usual interstitial pneumonia pattern: the impact of threshold-based volumetric CT analysis for the staging of idiopathic pulmonary fibrosis. *PLoS ONE.* (2016) 11:e0152505. doi: 10.1371/journal.pone.0152505

18. Lynch DA. Quantitative CT of fibrotic interstitial lung disease. *Chest.* (2007) 131:643-4. doi: 10.1378/chest.06-2955

19. Fukumoto J, Soundararajan R, Leung J, Cox R, Mahendrasah S, Muthavarapu N, et al. The role of club cell phenocroversion and migration in idiopathic pulmonary fibrosis. *Aging.* (2016) 8:3091-105. doi: 10.18632/aging.1 01115

20. Jensen-Taubman SM, Steinberg SM, Linnoila RI. Bronchiolization of the alveoli in lung cancer: pathology, patterns of differentiation and oncogene expression. *Int J Cancer.* (1998) 75:489-96. doi: 10.1002/(SICI)1097-0215(19980209)75:4<489::AID-ICJ1>3.0.CO;2-P

21. Gras D, Chanez P, Vachier I, Petit A, Bourdin A. Bronchial epithelium as a target for innovative treatments in asthma. *Pharmacol Ther.* (2013) 140:290-305. doi: 10.1016/j.pharmthera.2013.07.008

22. Alder JK, Chen JJ.-L, Lancaster L, Danoff S, Su S.-c, Cogan JD, et al. Short telomeres are a risk factor for idiopathic pulmonary fibrosis. *Proc Natl Acad Sci USA.* (2008) 105:13051-6. doi: 10.1073/pnas.08042 80105

23. Adnot S, Amsellem V, Boyer L, Marcos E, Saker M, Houssaini A, et al. Telomere dysfunction and cell senescence in chronic lung diseases: therapeutic potential. *Pharmacol Ther.* (2015) 153:125-34. doi: 10.1016/j.pharmthera.2015.06.007

24. Chilosi M, Carloni A, Rossi A, Poletti V. Premature lung aging and cellular senescence in the pathogenesis of idiopathic pulmonary fibrosis and COPD/emphysema. *Transl Res.* (2013) 162:156-73. doi: 10.1016/j.trsl.2013.06.004

25. Best AC, Meng J, Lynch AM, Bozic CM, Miller D, Grunwald GK, et al. Idiopathic pulmonary fibrosis: physiologic tests, quantitative CT indexes, and CT visual scores as predictors of mortality. *Radiology.* (2008) 246:935-40. doi: 10.1148/radiol.2463062200

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Reynaud, Ahmed, Serre, Knabe, Bommart, Suehs, Vachier, Berthet, Romagnoli, Vernisse, Mallet, Gamez and Bourdin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.