# Development of a multivariate prediction model for early-onset bronchiolitis obliterans syndrome and restrictive allograft syndrome in lung transplantation

Angela Koutsokera, Pierre J. Royer, Jean P. Antonietti, Andreas Fritz, Christian Benden, John D. Aubert, Adrien Tissot, Karine Botturi, Antoine Roux, Martine L. Reynaud-Gaubert, et al.

► To cite this version:

Angela Koutsokera, Pierre J. Royer, Jean P. Antonietti, Andreas Fritz, Christian Benden, et al. Development of a multivariate prediction model for early-onset bronchiolitis obliterans syndrome and restrictive allograft syndrome in lung transplantation. Frontiers in Medicine, 2017, 4, 10.3389/fmed.2017.00109. [hal-01608943](https://hal.science/hal-01608943)

HAL Id: hal-01608943

https://hal.science/hal-01608943

Submitted on 3 Oct 2017

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>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Development of a Multivariate Prediction Model for Early-Onset Bronchiolitis Obliterans Syndrome and Restrictive Allograft Syndrome in Lung Transplantation

**OPEN ACCESS**

**Edited by:**

Kian Fan Chung,
Imperial College London,
United Kingdom

**Reviewed by:**

Mark Wewers,
The Ohio State University Columbus,
United States
Gunnar N. Hillerdal,
Karolinska University Hospital,
Sweden

**\*Correspondence:**

Angela Koutsokera
angela.koutsokera@chuv.ch

¹The collaboration of the French-Belgian Cohort of Lung Transplantation (COLT) and the Swiss Transplant Cohort Study (STCS).

**Specialty section:**

This article was submitted to Pulmonary Medicine, a section of the journal *Frontiers in Medicine*

**Received:** 22 March 2017

**Accepted:** 30 June 2017

**Published:** 17 July 2017

**Citation:**

Koutsokera A, Royer PJ, Antonietti JP, Fritz A, Benden C, Aubert JD, Tissot A, Botturi K, Roux A, Reynaud-Gaubert ML, Kessler R, Dromer C, Mussot S, Mal H, Mornex J-F, Guillemain R, Knoop C, Dahan M, Soccal PM, Claustre J, Sage E, Gomez C. Magnan A, Pison C and Nicod LP (2017) Development of a Multivariate Prediction Model for Early-Onset Bronchiolitis Obliterans Syndrome and Restrictive Allograft Syndrome in Lung Transplantation. *Front. Med.* 4:109. doi: 10.3389/fmed.2017.00109

Angela Koutsokera¹*, Pierre J. Royer², Jean P. Antonietti¹, Andreas Fritz³, Christian Benden⁴, John D. Aubert¹, Adrien Tissot², Karine Botturi², Antoine Roux⁵, Martine L. Reynaud-Gaubert⁶, Romain Kessler⁷, Claire Dromer⁸, Sacha Mussot⁹, Hervé Mal¹⁰, Jean-François Mornex¹¹, Romain Guillemain¹², Christiane Knoop¹³, Marcel Dahan¹⁴, Paola M. Soccal¹⁵, Johanna Claustre¹⁶, Edouard Sage¹⁷, Carine Gomez⁶, Antoine Magnan², Christophe Pison¹⁶, Laurent P. Nicod¹ and The SysCLAD Consortium¹⁸†

¹Division of Pulmonary Medicine, Centre Hospitalier Universitaire Vaudois (CHUV), University of Lausanne, Lausanne, Switzerland, ²Institut du thorax, INSERM UMR 1087/CNRS UMR 6291, CHU de Nantes, Université de Nantes, Nantes, France, ³Biomax Informatics AG, Planegg, Germany, ⁴Division of Pulmonary Medicine, University Hospital Zurich, Zurich, Switzerland, ⁵Pneumology, Adult CF Center and Lung transplantation Department, Foch Hospital, Université Versailles Saint-Quentin-en-Yvelines, UPRES EA220, Suresnes, France, ⁶Pulmonary Medicine, CF Center and Lung Transplantation Department, Centre Hospitalier Universitaire Nord, CNRS UMR 6236 Aix-Marseille Université, Marseille, France, ⁷Lung Transplant Center, Hôpitaux universitaires de Strasbourg, Strasbourg, France, ⁸Service des Maladies respiratoires, Hôpital Haut Lévèque, Pessac, France, ⁹Service de Chirurgie Thoracique, Vasculaire et Transplantation Cardiopulmonaire, Hôpital Marie Lannelongue, Le Plessis Robinson, France, ¹⁰Service de Pneumologie et Transplantation pulmonaire, Hôpital Bichat, Université Denis Diderot, INSERM UMR1152, Paris, France, ¹¹Université de Lyon, INRA UMR 754, Hospices civils de Lyon, Lyon, France, ¹²Assistance Publique Hôpitaux de Paris, Paris, France, ¹³Department of Chest Medicine, Erasme University Hospital, Brussels, Belgium, ¹⁴CHU Larrey, Toulouse, France, ¹⁵Division of Pulmonary Medicine, Geneva University Hospitals, Geneva, Switzerland, ¹⁶Clinique Universitaire de Pneumologie, Pôle Thorax et Vaisseaux, CHU Grenoble, INSERM 1055, Université Grenoble Alpes, Grenoble, France, ¹⁷Thoracic Surgery Department, Foch Hospital, Université Versailles Saint-Quentin-en-Yvelines, UPRES EA220, Suresnes, France, ¹⁸Full List of the Authors Is Given in the Section SysCLAD Consortium Members

**Background:** Chronic lung allograft dysfunction and its main phenotypes, bronchiolitis obliterans syndrome (BOS) and restrictive allograft syndrome (RAS), are major causes of mortality after lung transplantation (LT). RAS and early-onset BOS, developing within 3 years after LT, are associated with particularly inferior clinical outcomes. Prediction models for early-onset BOS and RAS have not been previously described.

**Methods:** LT recipients of the French and Swiss transplant cohorts were eligible for inclusion in the SysCLAD cohort if they were alive with at least 2 years of follow-up but less than 3 years, or if they died or were retransplanted at any time less than 3 years. These patients were assessed for early-onset BOS, RAS, or stable allograft function by an adjudication committee. Baseline characteristics, data on surgery, immunosuppression, and year-1 follow-up were collected. Prediction models for BOS and RAS were developed using multivariate logistic regression and multivariate multinomial analysis.

**Results:** Among patients fulfilling the eligibility criteria, we identified 149 stable, 51 BOS, and 30 RAS subjects. The best prediction model for early-onset BOS and RAS included the underlying diagnosis, induction treatment, immunosuppression, and year-1 class
II donor-specific antibodies (DSAs). Within this model, class II DSAs were associated with BOS and RAS, whereas pre-LT diagnoses of interstitial lung disease and chronic obstructive pulmonary disease were associated with RAS.

**Conclusion:** Although these findings need further validation, results indicate that specific baseline and year-1 parameters may serve as predictors of BOS or RAS by 3 years post-LT. Their identification may allow intervention or guide risk stratification, aiming for an individualized patient management approach.

**Keywords:** chronic lung allograft dysfunction, bronchiolitis obliterans syndrome, restrictive allograft syndrome, chronic rejection, predictive model

## INTRODUCTION

Chronic lung allograft dysfunction (CLAD) is the principal cause of poor long-term survival in lung transplantation (LT) and, although no international consensus definition has been developed to date, CLAD refers to the persistent decline of the forced expiratory volume in one second (FEV₁) that cannot be attributed to a specific cause other than chronic graft rejection (1, 2). Bronchiolitis obliterans syndrome (BOS) and restrictive allograft syndrome (RAS) are considered to be the two main phenotypes of CLAD (2, 3). BOS is characterized by persistent airflow obstruction in the absence of a restrictive ventilation defect and imaging studies that may be unremarkable or show air trapping (4). Early-onset BOS, developing within 2 or 3 years after LT, is associated with particularly unfavorable outcomes, resulting in high morbidity and mortality soon after LT (5, 6). For RAS, various diagnostic criteria have been used in different studies, but overall, RAS is characterized by a restrictive ventilation defect and radiological signs of fibrosis or infiltrates (7-9). Prognosis is worse for RAS than BOS, and for RAS, time of diagnosis after transplant does not seem to influence survival (10). Diagnosis of BOS or RAS requires exclusion of alternative diagnoses, and this may be challenging (1).

Although the pathogenic mechanisms and the risk factors implicated in BOS and RAS are not fully elucidated (4, 11), recent literature provides increasing evidence and novel insights. Concerning BOS, most studies emphasize the role of multiple allo-immune and non-immune mechanisms (4, 11-13), but information on the risk factors of RAS or a comparative analysis of RAS and BOS is limited (10, 14). Identification of clinical risk factors associated with specific CLAD phenotypes is of particular clinical importance as it may assist patient risk stratification, optimize follow-up, or allow early intervention for potentially modifiable factors.

Multivariate prediction models are being increasingly developed to help health-care providers estimate the probability of a future disease (15). A few studies described clinical predictive models for BOS (16-19), but to our knowledge, no study has used this approach for BOS and RAS. The SysCLAD (systems prediction of CLAD) study is a collaborative project merging data from two large European LT cohorts. It aims to develop prediction models for early-onset BOS and RAS by implementing a multilevel approach and incorporating multiple clinical and laboratory parameters (20). The present study concerns the clinical arm of SysCLAD and describes the development of a multivariate predictive model for early-onset BOS and RAS. Our methodology and findings are reported in accordance with the TRIPOD recommendations (15). Some of the results of this study have been previously reported as conference abstracts (21, 22).

## MATERIALS AND METHODS

### Patient Population

This multicenter cohort study used data from two prospective European LT cohorts, the French-Belgian Cohort of LT (COLT) and the Swiss Transplant Cohort Study (STCS). COLT includes 11 French centers (Bordeaux, Grenoble, Lyon, Marseille, Nantes, Georges Pompidou Hospital, Hospital Bichat in Paris, Centre Chirurgical Marie Lennelongue in Le Plessis-Robinson, Strasbourg, Toulouse, Foch Hospital in Suresnes) and the Erasme center in Brussels, Belgium. It was established in October 2009. STCS is a multicenter cohort collecting data from all Swiss transplant programs (23, 24). The two Swiss LT centers, Zurich and Lausanne/Geneva, participate in STCS since its establishment in April 2008.

COLT and STCS patients were eligible for the SysCLAD cohort if they were alive with at least 2 years of follow-up but less than 3 years, or if they died or were retransplanted at any time less than 3 years. Only patients with a first LT were included in the study. Data collection concerned the periods between October 2009 (for COLT) or April 2008 (for STCS) and February 2014 (for both cohorts). Data harmonization, completion, and quality control were performed during the following 1.5 years. The corresponding national and local ethics committees approved the study protocol, and all participants provided written informed consent.

<footnote>**Abbreviations:** AIC, akaike information criterion; AUC, area under the curve; BMI, body mass index; BOS, bronchiolitis obliterans syndrome; CLAD, chronic lung allograft dysfunction; DSAs, donor-specific antibodies; FEV₁, forced expiratory volume in one second; FVC, forced vital capacity; HLA, human leukocyte antigen; ISHLT, International Society for Heart and Lung Transplantation; ILD/IPF, interstitial lung disease, idiopathic pulmonary fibrosis; LRA, logistic regression analysis; LT, lung transplantation; PGD, primary graft dysfunction; OR, odds ratio; RAS, restrictive allograft syndrome; rATG, rabbit antithymocyte globulin; ROC, receiver operating characteristic; t-AR, treated acute rejections; t-CMV, treated CMV infections; t-infections, treated infections; TLC, total lung capacity.</footnote>
# Assessment by the Adjudication Committee

Outcomes were established by an adjudication committee of LT specialists (Antoine Magnan, Christophe Pison, Antoine Roux, Martine L. Reynaud-Gaubert, Laurent P. Nicod, John D. Aubert, and Christian Benden) who evaluated pulmonary function tests (PFTs), imaging studies, and confounding factors. This committee formed a group opinion and identified LT recipients who remained stable or developed a definite BOS or RAS. Actions to blind assessment included data anonymization and an initial evaluation of the PFTs before assessment of imaging studies and confounding factors.

Bronchiolitis obliterans syndrome was defined as the persistent drop of FEV₁ in the absence of a restrictive defect and in the absence of confounding factors (4, 8). RAS was defined as the persistent decline of FEV₁ in the presence of a restrictive defect [i.e., decline of total lung capacity (TLC) and/or forced vital capacity (FVC)] and in the absence of confounding factors (8, 25). Imaging studies compatible with RAS were used as an additional diagnostic criterion. Patients were censored either at death/retransplantation or they were scored as uncensored at their last available assessment between 2 and 3 years of follow-up. More details on the elements assessed are presented below.

## Pulmonary Function Tests

In addition to specific time-point values recorded in the databases, PFTs were retrieved from LT centers. Spirometry was performed at every patient visit. The frequency of TLC measurements varied among centers: Zurich conducted TLC measurements at each patient visit, whereas in other centers, this was done at predefined assessment visits and when clinically indicated.

All hospitals calibrated their spirometers and body plethysmographers routinely according to the ATS/ERS guidelines. Home spirometry measurements were not used for diagnostic purposes in the present study. All pediatric patients included were able to perform PFTs reliably.

## Imaging Studies

Chest X-rays and chest CT scans of patients with declining PFTs were assessed to identify confounding factors for this decline or to support the PFT-based diagnosis of BOS or RAS. Unremarkable imaging studies or signs of air trapping (4) were considered supportive of BOS, whereas fibrosis, pleural thickening, or infiltrates in the absence of a confirmed infection were considered supportive of RAS (8, 25).

## Confounding Factors

Before establishing the diagnosis of BOS or RAS, we excluded allograft or extra-allograft factors that could cause the decline of PFTs (1, 4). Anastomotic, parenchymal, or thoracic wall abnormalities were identified using imaging studies and bronchoscopy reports. Whenever necessary, LT centers were contacted to obtain additional information.

## Establishment of the Diagnosis

A flow chart diagram of the assessed population is shown in **Figure 1**. During level-1 analysis, FEV₁ baseline was calculated

as the mean of the two best FEV₁ measured 3 weeks apart. At the last assessment time-point, FEV₁% of baseline was calculated to identify two groups of patients: (1) LT recipients without a persistent FEV₁ decline (stable patients) and (2) LT recipients with a persistent decline of FEV₁ < 80% of baseline. In the absence of confounding factors the diagnosis of CLAD was established and PFTs within 6 months were assessed to confirm the diagnosis.

During level-2 analysis, among CLAD cases, a decline of TLC to <90% of baseline was used to indicate RAS (8). When the number of available TLC values was insufficient for a confident diagnosis, a decline of FVC to <80% of baseline was used in combination with compatible imaging studies. For RAS cases, the FEV₁/FVC ratio was calculated to identify a purely restrictive (≥0.70) or mixed ventilation defect (ratio <0.70). The final study population was re-assessed for an improvement of FEV1 of ≥10% after 3 months of macrolide therapy.

Patients were characterized as "Not stable/No definite CLAD" if: (1) they were not stable but did not have CLAD, i.e., deceased or retransplanted patients for causes other than CLAD (including patients who died within 3 months after LT) and (2) they had one or more confounding factors contributing to the decrease of PFTs. Cases with confounding factors were assessed in detail by the adjudication committee and whenever possible patient evolution was followed over time. These cases were characterized as "no definite CLAD" when the predefined diagnostic criteria of BOS or RAS were not met despite this detailed assessment. After this evaluation, three groups (stable, BOS, and RAS) were created and used as outcomes for the multivariate models.

## Collection of Predictors

Clinical variables were chosen for their relevance regarding CLAD (1, 12), and subsequent data harmonization targeted to overcome discrepancies between the cohorts and among different centers. With a focus on those objectives, the following variables were collected and served as predictors in the multivariate models: (a) recipients' baseline characteristics, (b) donors' baseline characteristics, (c) data on induction treatment and surgery, and (d) data on follow-up: stage 3 primary graft dysfunction (PGD), maintenance immunosuppression, number of treated acute cellular rejection episodes during year 1 (Y1 t-AR), number of treated infections during year 1 (Y1 t-infections), number of treated CMV infections during year 1 (Y1 t-CMV), and donor-specific antibodies (DSAs). Details for specific variables are provided below.

Human leukocyte antigen (HLA) mismatches concerned the sum of A1, A2, B1, B2, DR1, and DR2 mismatches between recipients and donors (maximum possible sum: 6). Baseline DQ mismatches were not included in this sum as they were not available for all patients since the creation of the two cohorts.

Regarding PGD, for each patient, the worse stage developing within the first 72 h post-operatively (26, 27) is captured in the databases. Charts were reviewed to confirm accurate characterization according to established criteria (PaO₂/FiO₂ ratio < 200, radiological infiltrates and/or ECMO). Although we were able to confirm severe PGD, milder PGD stages were difficult to establish unequivocally in all cases and for that reason they were not included in the analysis.
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Flow chart diagram of the evaluated population.</figcaption></figure>
Y1 t-AR refers to the total number of treated events of grade A rejection or lymphocytic bronchiolitis (biopsy proven or clinically suspected). For harmonization reasons, we used this parameter instead of the cumulative A score [the latter depended on the frequency of conducted biopsies and the occurrence of Ax (inconclusive) results which were variable among centers].

Y1 t-infections refer to the total number of treated microbial, viral, and fungal infections (pulmonary and extra-pulmonary). Data collection especially for viral infections was not consistent between the databases and among centers, and for this reason, we did not use this variable. However, data on treated infections (including treated CMV infections) were recorded consistently and in detail.

Donor-specific antibody measurements were conducted at each study center. In order to alleviate methodological discrepancies, DSAs were considered positive when the result was validated as such by the corresponding immunology laboratory and no unique cut-off for mean fluorescence intensity (MFI) was used (Table S1 in Supplementary Material).

For some parameters, clinical practices, definitions, or data collection were very heterogeneous or inconsistent among centers, not allowing harmonization under one common definition. These parameters did not pass the quality control to be included in the analysis and were considered missing variables. Diagnosis of antibody-mediated rejection was not possible to establish with certainty for all patients since the creation of the cohorts, not only due to the evolving definition criteria but also due to the variability of clinical practices among centers and over time. The lack of an international standardization during data collection [International Society for Heart and Lung Transplantation (ISHLT) guideline published in 2016 (28)] led to subsequent discrepancies in the interpretation and management of positive results. Investigations and definitions of gastroesophageal reflux varied significantly among centers, depending on different local clinical practices and era of transplant. For that reason, it was not possible to harmonize this variable under one common definition.

## Data Quality Control and Harmonization

Data quality control strategies were implemented by three teams: (a) COLT coordinators, COLT datacenter (Informatique Données Base Centralisées—IDBC, St Luce/Loire), (b) STCS coordinators, STCS datacenter (Basel, Switzerland), and (c) SysCLAD coordinators, SysCLAD datacenter (Biomax, Germany).

Within each cohort, data were collected prospectively by the clinical research assistants according to the cohort-specific data dictionary and a predefined acquisition methodology. The corresponding datacenters conducted regular, automatic data quality validation assessments including but not restricted to: completeness, availability of the required data, data type check, data ranges based on realistic expected values, and rule-based inconsistency checks (e.g., diagnosis of “stable” is in conflict with an FEV1 < 80% of baseline). Database coordinators addressed datacenter originating queries to verify or complete data.

To identify inconsistencies in the definitions and in data completion between COLT and STCS, a harmonization working group (Christophe Pison, Laurent P. Nicod, Pierre J. Royer, Angela Koutsokera, and Andreas Fritz) defined which data were sufficiently harmonized, which needed refinement or completion

before achieving harmonization, and which were not possible to harmonize (latter considered as missing variables). The conclusions and suggestions of this working team were presented at a regular basis to the SysCLAD consortium members who validated the final decisions. Details on the collected clinical parameters, quality control, and missing variables are reported in Tables S1–S3 in Supplementary Material.

## Statistical Analysis

Continuous variables were reported as median (IQR) and categorical as *n* (%). Mann–Whitney *U*-test was used for two-group comparisons, Kruskal–Wallis for three groups, and chi-square for categorical variables.

Univariate and multivariate logistic regression analyses (LRAs) were used to develop a prediction model for early-onset CLAD. Univariate and multivariate multinomial analyses were used to develop a prediction model for early-onset BOS or RAS. Independent explanatory variables were identified by backward-forward and forward-backward elimination techniques (29). The entry and removal criteria from the equation were a probability of likelihood ratio <0.05 and >0.10, respectively. The “center of transplantation” was coded as an indicator variable and, during multivariate analysis, systematically forced into the model as one of the explanatory variables to provide adjustment for center effect. For model performance and to test the overall fit of the models, we used the Hosmer–Lemeshow goodness-of-fit and the $R^2$ (Cox–Snell and McFadden) tests. Results were reported as odds ratios (ORs, 95% CI).

Variables of the best performing multivariate LRA model were used to create a receiver operating characteristic (ROC) curve. The equation for the prediction of CLAD corresponding to this ROC curve was created for the studied population. Probabilities of stability, BOS, or RAS were calculated for all modalities of each significant independent variable in the multivariate multinomial analysis.

Cox proportional hazards analysis was used to identify risk factors for mortality in the studied population. Variables with *p*-value of 0.1 or less in the univariate analysis were included in the multivariate model. A backward/forward selection procedure with Akaike information criterion was used to identify risk factors for mortality, and survival curves were obtained from Kaplan–Meier estimates.

Statistical significance was defined as *p* < 0.05. All statistical analyses were performed using R 3.2.0 and SPSS 22 software versions.

## RESULTS

### Study Population and Subgroups of Patients

As shown in **Figure 1**, from 1,263 LT included in COLT and STCS, 422 adult and pediatric recipients fulfilled the SysCLAD eligibility criteria. Following adjudication, a definite diagnosis was established for 230 subjects, transplanted in 12 different centers, and these patients were analyzed further. Patients’ characteristics are presented in **Table 1**.

The following groups were analyzed: (1) stable recipients (*n* = 149), CLAD (*n* = 81), and (2) CLAD patients were further
<table><caption>TABLE 1 | Characteristics of the studied population and comparisons between stable recipients and recipients diagnosed with BOS or RAS by 3 years post-LT.</caption><thead><tr><th></th><th></th><th>All patients (n = 230)</th><th>Stable recipients (n = 149)</th><th>BOS (n = 51)</th><th>RAS (n = 30)</th><th>p-Value</th></tr></thead><tbody><tr><td colspan="7"><strong>Recipients' characteristics</strong></td></tr><tr><td>Age (years)</td><td></td><td>46.5 [29, 58]<sup>a</sup></td><td>43 [30, 56]</td><td>47 [27, 57.5]</td><td>57 [38, 62]</td><td>0.0517</td></tr><tr><td>Gender</td><td>Male</td><td>118 (51.3)</td><td>78 (52)</td><td>23 (45.1)</td><td>17 (56.7)</td><td>0.55</td></tr><tr><td>BMI (kg/m<sup>2</sup>)</td><td></td><td>19.9 [17.9, 24.2]</td><td>19.6 [17.8, 23.7]</td><td>20.1 [17.9, 23.1]</td><td>23.6 [18.6, 26.8]</td><td><strong>0.0341</strong></td></tr><tr><td rowspan="4">Blood group</td><td>A</td><td>97 (42.4)</td><td>66 (44.3)</td><td>20 (40)</td><td>11 (36.7)</td><td rowspan="4">0.63</td></tr><tr><td>AB</td><td>15 (6.6)</td><td>11 (7.4)</td><td>3 (6)</td><td>1 (3.3)</td></tr><tr><td>B</td><td>29 (12.7)</td><td>20 (13.4)</td><td>7 (14)</td><td>2 (6.7)</td></tr><tr><td>O</td><td>88 (38.4)</td><td>52 (34.9)</td><td>20 (40)</td><td>16 (53.3)</td></tr><tr><td rowspan="4">Underlying diagnosis</td><td>COPD</td><td>61 (26.5)</td><td>37 (24.8)</td><td>14 (27.5)</td><td>10 (33.3)</td><td rowspan="4"><strong>&lt;0.001</strong></td></tr><tr><td>CF</td><td>88 (38.3)</td><td>67 (45)</td><td>15 (29.4)</td><td>6 (20)</td></tr><tr><td>ILD/IPF</td><td>43 (18.7)</td><td>20 (13.4)</td><td>10 (19.6)</td><td>13 (43.3)</td></tr><tr><td>Others</td><td>38 (16.5)</td><td>25 (16.8)</td><td>12 (23.5)</td><td>1 (3.3)</td></tr><tr><td>Smoking history</td><td>Yes</td><td>108 (47.8)</td><td>63 (43.2)</td><td>26 (52)</td><td>19 (63.3)</td><td>0.10</td></tr><tr><td>HLA mismatches</td><td></td><td>5 [4,5]</td><td>5 [4, 5]</td><td>5 [4, 5]</td><td>5 [4, 5.75]</td><td>0.81</td></tr><tr><td colspan="7"><strong>Donors' characteristics</strong></td></tr><tr><td>Age (years)</td><td></td><td>45 [32, 55]<sup>b</sup></td><td>43 [31, 56]</td><td>47 [37.5, 54]</td><td>49 [37, 53]</td><td>0.61</td></tr><tr><td>Gender</td><td>Male</td><td>124 (53.9)</td><td>83 (55.7)</td><td>23 (45.1)</td><td>18 (60)</td><td>0.33</td></tr><tr><td rowspan="4">Blood group</td><td>A</td><td>95 (45.3)</td><td>64 (43)</td><td>21 (41.2)</td><td>10 (33.3)</td><td rowspan="4">0.89</td></tr><tr><td>AB</td><td>10 (4.3)</td><td>6 (4)</td><td>3 (5.9)</td><td>1 (3.3)</td></tr><tr><td>B</td><td>29 (12.6)</td><td>19 (12.8)</td><td>7 (13.7)</td><td>3 (10)</td></tr><tr><td>O</td><td>96 (41.7)</td><td>60 (40.3)</td><td>20 (39.2)</td><td>16 (53.3)</td></tr><tr><td>Smoking history</td><td>Yes</td><td>95 (42.8)</td><td>59 (40.7)</td><td>24 (50)</td><td>12 (41.4)</td><td>0.521</td></tr><tr><td colspan="7"><strong>Intervention</strong></td></tr><tr><td rowspan="4">Type of intervention</td><td>Double lung</td><td>192 (83.5)</td><td>127 (85.2)</td><td>37 (72.5)</td><td>28 (93.3)</td><td rowspan="4">0.20</td></tr><tr><td>Single lung</td><td>5 (2.2)</td><td>16 (10.7)</td><td>11 (21.6)</td><td>2 (6.7)</td></tr><tr><td>Lobar</td><td>4 (1.7)</td><td>2 (1.3)</td><td>2 (3.9)</td><td>0 (0)</td></tr><tr><td>Heart-lung</td><td>29 (12.6)</td><td>4 (2.7)</td><td>1 (2)</td><td>0 (0)</td></tr><tr><td>Max cold ischemia time (min)</td><td></td><td>320 [275, 380]</td><td>330 [281, 380]</td><td>320 [247, 397]</td><td>304 [284, 371]</td><td>0.62</td></tr><tr><td rowspan="3">Induction treatment</td><td>Basiliximab</td><td>108 (47)</td><td>80 (53.7)</td><td>14 (27.5)</td><td>14 (46.7)</td><td rowspan="3"><strong>0.0016</strong></td></tr><tr><td>None</td><td>65 (28.3)</td><td>40 (26.8)</td><td>14 (27.5)</td><td>11 (36.7)</td></tr><tr><td>rATG</td><td>57 (24.8)</td><td>29 (19.5)</td><td>23 (45.1)</td><td>5 (16.7)</td></tr><tr><td colspan="7"><strong>Follow-up</strong></td></tr><tr><td>PGD stage 3</td><td>Severe</td><td>7 (3.1)</td><td>5 (3.4)</td><td>1 (2)</td><td>1 (3.3)</td><td>0.87</td></tr><tr><td rowspan="2">Immunosuppression</td><td>Cyclosporin</td><td>91 (39.6)</td><td>54 (36.2)</td><td>19 (37.3)</td><td>18 (60)</td><td rowspan="2"><strong>0.0488</strong></td></tr><tr><td>Tacrolimus</td><td>139 (60.4)</td><td>95 (63.8)</td><td>32 (62.7)</td><td>12 (40)</td></tr><tr><td>Y1 t-AR</td><td></td><td>0 [0, 1]</td><td>0 [0, 1]</td><td>0 [0, 1]</td><td>0 [0, 1]</td><td>0.65</td></tr><tr><td>Y1 t-infections</td><td></td><td>1 [0, 2]</td><td>1 [0, 2]</td><td>1 [0, 2]</td><td>1 [0, 2]</td><td>0.99</td></tr><tr><td>Y1 t-CMV</td><td></td><td>0 [0, 0]</td><td>0 [0, 0]</td><td>0 [0, 1]</td><td>0 [0, 0]</td><td>0.37</td></tr><tr><td>DSA before LT</td><td>Yes</td><td>42 (20.3)</td><td>20 (14.9)</td><td>15 (33.3)</td><td>7 (25.0)</td><td><strong>0.024</strong></td></tr><tr><td>Y1 DSAs (I or II)</td><td>Yes</td><td>48 (21.3)</td><td>21 (14.3)</td><td>17 (35.4)</td><td>10 (33.3)</td><td><strong>0.002</strong></td></tr><tr><td>Y1 DSAs I</td><td>Yes</td><td>24 (10.6)</td><td>11 (7.5)</td><td>10 (20.4)</td><td>3 (10)</td><td><strong>0.039</strong></td></tr><tr><td>Y1 DSAs II</td><td>Yes</td><td>42 (18.7)</td><td>17 (11.6)</td><td>16 (33.3)</td><td>9 (30)</td><td><strong>0.001</strong></td></tr></tbody></table>

Statistically significant results are highlighted in bold.

Continuous variables are presented as median [IQR] and categorical variables as n (%).

<sup>a</sup>Recipient's age range = 14–68 years.

<sup>b</sup>Donor's age range = 10–78 years.

AR, acute cellular rejection episodes; BMI, body mass index; CF, cystic fibrosis; CMV, cytomegalovirus; COPD, chronic obstructive pulmonary disease; DSAs, donor-specific antibodies; HLA, human leukocyte antigen; ILD/IPF, interstitial lung disease/idiopathic pulmonary fibrosis; LT, lung transplantation; rATG, rabbit antithymocyte globulin; t, treated; PGD, primary graft dysfunction; Y1, year 1; BOS, bronchiolitis obliterans syndrome; RAS, restrictive allograft syndrome.

subcategorized as BOS (51) or RAS (n = 30). Among the BOS patients, 17.6% were stage 1, 29.4% stage 2, and 53% stage 3. Concerning RAS, 37% (n = 11) had a FEV₁/FVC ratio <0.70 indicating a mixed pattern (example case in Supplementary Material).

The final study population was assessed for an FEV₁ improvement of at least 10% after 3 months of macrolide therapy (either azithromycin or clarithromycin). FEV₁ of all CLAD patients remained below the 80% of baseline threshold, but six (7.4%) had a reversibility of at least 10% after 3 months of macrolides. Within the stable population, 74 (49.7%) received macrolides for at least

3 months (either for a non-sustained decline of FEV₁ after an acute event, or as part of a treatment for non-tuberculous mycobacteria). Twenty-two patients (14.8% of the stable patients) had an improvement of at least 10% of FEV₁ at 3 months of treatment. In all stable patients, FEV1 remained above the 80% of baseline threshold. In macrolide-responsive cases, the improvement of FEV₁ could not be attributed to macrolides alone, since during the same period, most patients had concomitant treatments for a viral or bacterial infection or an acute cellular rejection.

Among the CLAD patients, 35 died (for 19, BOS or RAS was the main cause of death) and 8 were retransplanted (5 due to BOS and
3 due to RAS). The “Not stable/No definite CLAD” group included 147 patients who died or were retransplanted for causes other than CLAD and 45 patients with at least 1 confounding factor contributing to the decrease of FEV₁ and/or TLC (detailed in **Figure 1**).

## Model Development and Model Regression Diagnostics

The parameters tested by backward-forward and forward-backward elimination techniques to identify independent explanatory variables were the following: recipient's age, donor's age, difference of recipient and donor's age, underlying diagnosis, recipient's smoking history before LT, recipient's body mass index (BMI), stage 3 PGD, max cold ischemia time, number of HLA mismatches, induction treatment, maintenance immunosuppression (cyclosporine vs. tacrolimus), Y1 t-AR, Y1 t-infections, Y1 t-CMV infections, DSAs before LT, Y1 DSAs, Y1 DSAs class I, and Y1 DSAs class II.

For the analyzed parameters, only a small number of data was missing (Table S2 in Supplementary Material). Single imputation was used to handle missing data. For nominal variables, missing values were replaced by the variable mode (i.e., most frequent value), and for numerical variables, missing values were replaced by the median value. The only exception was “DSAs before LT” which had the most missing data (23 of 230). In this case, missing data were coded as unknown.

For induction treatment, we used “basiliximab” as the baseline and for “underlying diagnosis” we chose “CF.” For the former, a

preliminary exploratory univariate analysis using “none” as the baseline did not provide statistically significant differences [basiliximab vs. none OR (95% CI) 0.882 (0.763, 1.020), *p* = 0.09 and rabbit antithymocyte globulin (rATG) vs. none 1.113 (0.941, 1.315), *p* = 0.214]. We considered that the statistically significant difference observed using basiliximab as the baseline was of clinical interest. In addition, basiliximab was the most frequently used induction agent within this study (Table S3 in Supplementary Material). The choice of CF as the baseline for the “underlying diagnosis” was done for the same reasons (statistical significant results in the preliminary analysis, most frequent diagnosis in the study population). We considered that this comparison would be of particular clinical interest.

With regard to model regression diagnostics: (a) the dependent variable was either binary (stable, BOS) or ordinal (stable, BOS, and RAS), (b) the dependent variables were coded appropriately, (c) the model fitted correctly after using a stepwise method, the Hosmer-Lemeshow goodness-of-fit and the R² test, (d) the data originated from independent samples (no paired samples) and we used the stepwise approach to avoid multicollinearity, (e) linearity of independent variables and log odds was observed in the graphics diagnostics, and (f) we used a sample size of 230 cases to test 18 variables (i.e., more than 10 samples per independent variable).

## Prediction Model for Early-Onset CLAD

**Table 2** shows the results of the univariate and multivariate LRA. The variables of the best performing multivariate LRA

<table><caption>**TABLE 2** | Risk factors for the development of CLAD (multivariate analysis adjusted for center effect) by 3 years post-LT.</caption><thead><tr><th rowspan="2">Variable</th><th rowspan="2"></th><th colspan="2">Univariate analysis</th><th colspan="2">Multivariate analysis</th></tr><tr><th>OR (95% CI)</th><th>p-Value</th><th>OR (95% CI)</th><th>p-Value</th></tr></thead><tbody><tr><td>Recipient age</td><td></td><td>1.002 (0.998, 1.006)</td><td>0.2288</td><td>0.971 (0.938, 1.006)</td><td>0.102</td></tr><tr><td>Donor age</td><td></td><td>1.002 (0.998, 1.006)</td><td>0.3172</td><td></td><td></td></tr><tr><td>Difference of R/D age</td><td></td><td>1.000 (0.997, 1.004)</td><td>0.9041</td><td></td><td></td></tr><tr><td>Recipient smoking</td><td>Yes</td><td>1.128 (0.996, 1.277)</td><td>0.0598</td><td></td><td></td></tr><tr><td>Recipient BMI</td><td></td><td>1.012 (0.998, 1.025)</td><td>0.0924</td><td></td><td></td></tr><tr><td rowspan="4">Underlying diagnosis</td><td>CF</td><td>Baseline</td><td></td><td>Baseline</td><td></td></tr><tr><td>COPD</td><td>1.167 (1.002, 1.361)</td><td><strong>0.049</strong></td><td>5.158 (1.444, 18.426)</td><td><strong>0.012</strong></td></tr><tr><td>ILD/IPF</td><td>1.345 (1.133, 1.596)</td><td><strong>&lt;0.001</strong></td><td>9.429 (2.291, 38.807)</td><td><strong>0.002</strong></td></tr><tr><td>Others</td><td>1.109 (0.928, 1.326)</td><td>0.2573</td><td>2.435 (0.783, 7.569)</td><td>0.124</td></tr><tr><td>Sum of HLA mismatches</td><td></td><td>1.011 (0.948, 1.077)</td><td>0.7396</td><td></td><td></td></tr><tr><td>Max cold ischemia time</td><td></td><td>1.000 (0.999, 1.001)</td><td>0.6245</td><td></td><td></td></tr><tr><td rowspan="2">Induction treatment</td><td>Basiliximab</td><td>Baseline</td><td></td><td>Baseline</td><td></td></tr><tr><td>None</td><td>1.134 (0.981, 1.310)</td><td>0.0914</td><td>1.382 (0.451, 4.236)</td><td>0.572</td></tr><tr><td>PGD stage 3</td><td>Yes</td><td>1.261 (1.084, 1.467)</td><td><strong>0.0029</strong></td><td>3.519 (0.946, 13.085)</td><td>0.060</td></tr><tr><td rowspan="2">Immunosuppression</td><td>Cyclosporin</td><td>Baseline</td><td></td><td></td><td></td></tr><tr><td>Tacrolimus</td><td>0.914 (0.805, 1.037)</td><td>0.1635</td><td></td><td></td></tr><tr><td>Y1 t-AR</td><td></td><td>1.003 (0.938, 1.072)</td><td>0.9392</td><td></td><td></td></tr><tr><td>Y1 t-infections</td><td></td><td>1.003 (0.970, 1.036)</td><td>0.8784</td><td></td><td></td></tr><tr><td>Y1 t-CMV</td><td></td><td>1.019 (0.899, 1.155)</td><td>0.764</td><td></td><td></td></tr><tr><td>DSA before LT</td><td>Yes</td><td>1.240 (1.056, 1.455)</td><td><strong>0.0092</strong></td><td></td><td></td></tr><tr><td>Y1 DSAs (I or II)</td><td>Yes</td><td>1.316 (1.135, 1.526)</td><td><strong>0.0004</strong></td><td></td><td></td></tr><tr><td>Y1 DSAs I</td><td>Yes</td><td>1.240 (1.014, 1.515)</td><td><strong>0.0370</strong></td><td></td><td></td></tr><tr><td>Y1 DSAs II</td><td>Yes</td><td>1.357 (1.162, 1.585)</td><td><strong>0.0001</strong></td><td>4.221 (1.784, 9.991)</td><td><strong>0.001</strong></td></tr></tbody></table>

Statistically significant results are highlighted in bold.

Results are presented as OR (95% confidence intervals). For binary variables, the “No” group has been considered the baseline group.

AR, acute cellular rejection episodes; BMI, body mass index; CF, cystic fibrosis; CMV, cytomegalovirus; COPD, chronic obstructive pulmonary disease; DSAs, donor-specific antibodies; HLA, human leukocyte antigen; ILD/IPF, interstitial lung disease/idiopathic pulmonary fibrosis; LT, lung transplantation; t, treated; PGD, primary graft dysfunction; rATG, rabbit antithymocyte globulin; Y1, year 1; CLAD, chronic lung allograft dysfunction; OR, odds ratio.
model for early-onset CLAD were used to create a multivariate
prediction model for estimating probability of developing early-
onset CLAD. The variables included recipient age, underlying
diagnosis, induction treatment and presence of Y1 class II DSAs,
and the prediction model generated a ROC curve (Figure 2A)
with an area under the curve (AUC) of 0.766. The equation cor-
responding to this ROC curve, its limitations, and some clinical
examples for the prediction of early-onset CLAD in the studied
population are presented in Table S4 in Supplementary Material.
The ROC curve of the same model but without using Y1 class II
DSA provided an AUC = 0.730. The difference between the AUC
of the complete model and the model without Y1 class II DSA did
not reach statistical significance (p = 0.074, Figure 2B).

Of note, increasing recipient age had a weak protective effect
over CLAD development (OR 0.971, *p* = 0.102 in the multivariate
LRA, factor −0.03 in the ROC curve equation). A graphic repre-
sentation of early-onset CLAD diagnosis according to recipient
age showed a tendency for a U-shaped distribution. However,
when the recipients' age was tested for a quadratic (i.e., non-
linear) effect, *p*-value was not statistically significant (Figure S1
in Supplementary Material).

Prediction Models for Early-Onset BOS
or RAS

**Tables 3 and 4** present the results of the univariate and multivariate multinomial analysis. The multivariate multinomial predictive model included the underlying diagnosis, maintenance immunosuppression, induction treatment, and Y1 class II DSAs. Within the model, Y1 class II DSAs were associated with BOS (OR 3.83) and RAS (OR 6.97). Compared to CF, interstitial lung disease/idiopathic pulmonary fibrosis (ILD/IPF) (OR 5.47) and chronic obstructive pulmonary disease (COPD) (OR 3.86) were associated with a higher risk for RAS. Induction treatment and

maintenance immunosuppression were included in the best
prediction model, although they were not statistically significant.
Probabilities for BOS and RAS for each significant independent
variable are shown in **Figure 3**.

**Survival Analysis**

Risk factors for mortality within 3 years after LT are shown in Table S5 in Supplementary Material, and survival curves obtained from Kaplan-Meier estimates are displayed in Figure S2 in Supplementary Material. In the multivariate analysis, the recipient previous smoking status, a diagnosis of IPF/ILD as compared to a diagnosis of CF, the sum of HLA mismatches, and the presence of DSA during year 1 (either class I or II) were associated with a higher risk of mortality.

DISCUSSION

This study analyzed clinical data of two large European LT
cohorts aiming to develop multivariate prediction models for
early-onset BOS or RAS. The created models provide prob-
abilities for a binary outcome (stable vs. CLAD), not taking into
account the probability of being in the “not stable/not CLAD”
group. Our main findings are: first, the multivariate prediction
model for CLAD included recipient age, underlying diagnosis,
type of induction treatment, and Y1 class II DSAs. A model,
using baseline variables only (i.e., exclusion of Y1 class II DSAs),
had a predictive capacity similar to the complete model. Second,
the multivariate prediction models for BOS and RAS included
underlying diagnosis, type of induction treatment, maintenance
immunosuppression, and Y1 class II DSAs, but were not identi-
cal. Year-1 class II DSAs were associated with both BOS and
RAS, whereas pre-LT diagnoses of ILD/IPF and COPD were
associated with RAS.

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

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

**FIGURE 2** | (A) Receiver operating characteristic (ROC) analysis of the best performing model for the prediction of early-onset chronic lung allograft dysfunction after adjusting for center effect [area under the curve (AUC) = 0.766, SE = 0.0325, 95% CI 0.703–0.830]. (B) ROC curve of the model, with and without the parameter of “Y1 class II DSA” (AUC 0.766 vs. 0.730, respectively, *p* = 0.074).
<table><caption>**TABLE 3** | Risk factors for BOS and RAS by 3 years post-LT as compared to stable recipients (univariate multinomial analysis).</caption><thead><tr><th rowspan="2">Variable</th><th rowspan="2"></th><th colspan="2">BOS</th><th colspan="2">RAS</th></tr><tr><th>OR (95% CI)</th><th>p-Value</th><th>OR (95% CI)</th><th>p-Value</th></tr></thead><tbody><tr><td>Recipient age</td><td></td><td>1.000 (0.980, 1.021)</td><td>0.974</td><td>1.031 (1.003, 1.060)</td><td><strong>0.030</strong></td></tr><tr><td>Donor age</td><td></td><td>1.005 (0.985, 1.025)</td><td>0.623</td><td>1.015 (0.990, 1.040)</td><td>0.235</td></tr><tr><td>Difference of R/D age</td><td></td><td>0.996 (0.979, 1.014)</td><td>0.683</td><td>1.008 (0.987, 1.030)</td><td>0.433</td></tr><tr><td>Recipient smoking history</td><td>Yes</td><td>1.420 (0.750, 2.687)</td><td>0.282</td><td>2.358 (1.048, 5.303)</td><td><strong>0.038</strong></td></tr><tr><td>Recipient BMI</td><td></td><td>1.014 (0.944, 1.089)</td><td>0.704</td><td>1.107 (1.024, 1.197)</td><td><strong>0.011</strong></td></tr><tr><td rowspan="4">Underlying diagnosis</td><td>CF</td><td>Baseline</td><td></td><td>Baseline</td><td></td></tr><tr><td>COPD</td><td>1.690 (0.736, 3.882)</td><td>0.216</td><td>3.018 (1.016, 8.966)</td><td><strong>0.047</strong></td></tr><tr><td>ILD/IPF</td><td>2.233 (0.870, 5.736)</td><td>0.095</td><td>7.258 (2.444, 21.559)</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Other</td><td>2.144 (0.833, 5.207)</td><td>0.092</td><td>0.447 (0.051, 3.897)</td><td>0.466</td></tr><tr><td>Sum of HLA mismatches</td><td></td><td>0.973 (0.706, 1.342)</td><td>0.870</td><td>1.180 (0.771, 1.805)</td><td>0.446</td></tr><tr><td>Max cold ischemia time</td><td></td><td>1.000 (0.996, 1.003)</td><td>0.904</td><td>0.998 (0.994, 1.003)</td><td>0.472</td></tr><tr><td rowspan="3">Induction treatment</td><td>Basiliximab</td><td>Baseline</td><td></td><td>Baseline</td><td></td></tr><tr><td>None</td><td>2.000 (0.870, 4.598)</td><td>0.103</td><td>1.571 (0.654, 3.774)</td><td>0.312</td></tr><tr><td>rATG</td><td>4.532 (2.060, 9.972)</td><td><strong>&lt;0.001</strong></td><td>0.985 (0.326, 2.977)</td><td>0.979</td></tr><tr><td>PGD stage 3</td><td>Yes</td><td>0.576 (0.066, 5.050)</td><td>0.618</td><td>0.993 (0.112, 8.819)</td><td>0.995</td></tr><tr><td rowspan="2">Immunosuppression</td><td>Cyclosporin</td><td>Baseline</td><td></td><td>Baseline</td><td></td></tr><tr><td>Tacrolimus</td><td>0.957 (0.495, 1.850)</td><td>0.897</td><td>0.379 (0.170, 0.846)</td><td><strong>0.018</strong></td></tr><tr><td>Y1 t-AR</td><td></td><td>1.110 (0.802, 1.536)</td><td>0.530</td><td>0.826 (0.505, 1.352)</td><td>0.448</td></tr><tr><td>Y1 t-infections</td><td></td><td>0.962 (0.803, 1.153)</td><td>0.678</td><td>1.087 (0.904, 1.308)</td><td>0.375</td></tr><tr><td>Y1 t-CMV</td><td></td><td>1.303 (0.721, 2.355)</td><td>0.380</td><td>0.763 (0.304, 1.915)</td><td>0.565</td></tr><tr><td>DSA before LT</td><td>Yes</td><td>2.850 (1.305, 6.223)</td><td><strong>0.009</strong></td><td>1.900 (0.714, 5.055)</td><td>0.199</td></tr><tr><td>Y1 DSAs (I or II)</td><td>Yes</td><td>3.048 (1.450, 6.406)</td><td><strong>0.003</strong></td><td>3.048 (1.254, 7.409)</td><td><strong>0.014</strong></td></tr><tr><td>Y1 DSAs I</td><td>Yes</td><td>3.060 (1.214, 7.714)</td><td><strong>0.018</strong></td><td>1.394 (0.364, 5.332)</td><td>0.628</td></tr><tr><td>Y1 DSAs II</td><td>Yes</td><td>3.550 (1.631, 7.726)</td><td><strong>0.001</strong></td><td>3.328 (1.313, 8.434)</td><td><strong>0.011</strong></td></tr></tbody></table>

Statistically significant results are highlighted in bold.

Results are presented as OR (95% confidence intervals). For binary variables, the “No” group has been considered the baseline group.

AR, acute cellular rejection episodes; BMI, body mass index; CF, cystic fibrosis; CMV, cytomegalovirus; COPD, chronic obstructive pulmonary disease; D, donor; DSAs, donor-specific antibodies; HLA, human leukocyte antigen; ILD/IPF, interstitial lung disease/idiopathic pulmonary fibrosis; LT, lung transplantation; R, recipient; t, treated; PGD, primary graft dysfunction; rATG, rabbit antithymocyte globulin; Y1, year 1; BOS, bronchiolitis obliterans syndrome; RAS, restrictive allograft syndrome; OR, odds ratio.

<table><caption>**TABLE 4** | Risk factors for BOS and RAS by 3 years post-LT as compared to stable recipients (multivariate multinomial analysis).</caption><thead><tr><th rowspan="2">Variable</th><th rowspan="2"></th><th colspan="2">BOS</th><th colspan="2">RAS</th></tr><tr><th>OR (95% CI)</th><th>p-Value</th><th>OR (95% CI)</th><th>p-Value</th></tr></thead><tbody><tr><td rowspan="4">Underlying diagnosis</td><td>CF</td><td>Baseline</td><td></td><td>Baseline</td><td></td></tr><tr><td>COPD</td><td>1.606 (0.559, 4.610)</td><td>0.379</td><td>3.857 (1.041, 14.289)</td><td><strong>0.043</strong></td></tr><tr><td>ILD/IPF</td><td>2.436 (0.738, 8.036)</td><td>0.144</td><td>5.467 (1.482, 20.170)</td><td><strong>0.011</strong></td></tr><tr><td>Other</td><td>2.589 (0.848, 7.905)</td><td>0.095</td><td>0.230 (0.022, 2.398)</td><td>0.219</td></tr><tr><td rowspan="2">Immunosuppression</td><td>Cyclosporin</td><td>Baseline</td><td></td><td>Baseline</td><td></td></tr><tr><td>Tacrolimus</td><td>3.179 (0.704, 14.357)</td><td>0.133</td><td>0.670 (0.080, 5.590)</td><td>0.711</td></tr><tr><td rowspan="3">Induction treatment</td><td>Basiliximab</td><td>Baseline</td><td></td><td>Baseline</td><td></td></tr><tr><td>None</td><td>0.541 (0.117, 2.505)</td><td>0.432</td><td>4.528 (0.888, 23.076)</td><td>0.069</td></tr><tr><td>rATG</td><td>3.101 (0.681, 14.123)</td><td>0.144</td><td>2.393 (0.307, 18.674)</td><td>0.405</td></tr><tr><td>Y1 DSAs II</td><td>Yes</td><td>3.827 (1.459, 10.040)</td><td><strong>0.006</strong></td><td>6.965 (1.839, 26.376)</td><td><strong>0.004</strong></td></tr></tbody></table>

Statistically significant results are highlighted in bold.

Results are presented as OR (95% confidence intervals). For binary variables, the “No” group has been considered the baseline group.

CF, cystic fibrosis; COPD, chronic obstructive pulmonary disease; DSAs, donor-specific antibodies; ILD/IPF, interstitial lung disease/idiopathic pulmonary fibrosis; rATG, rabbit antithymocyte globulin; t, treated; Y1, year 1; BOS, bronchiolitis obliterans syndrome; RAS, restrictive allograft syndrome; OR, odds ratio.

Establishment of a definite diagnosis of BOS or RAS is often challenging. Despite the detailed evaluation by the adjudication committee, 10% (n = 45) of the assessed patients had confounding factors not allowing a definite diagnosis. For the differential diagnosis of BOS and RAS, TLC is considered the gold standard; however, in case of insufficient TLC values, FVC may be used (1, 8, 25). Imaging studies are not included in the current diagnostic algorithms, but may support or contradict the PFT-based diagnosis. This may be particularly helpful when FVC decrease

is associated with air trapping (pseudo-restriction) (30). The isolated evaluation of the FEV₁/FVC ratio in these cases may be misleading and not allow identification of mixed ventilation defects. In our study, 37% (n = 11) of RAS patients had an FEV₁/FVC ratio <0.70 suggesting a mixed PFTs pattern.

Prognostic prediction models for LT outcomes may assist patient risk stratification and improve follow-up strategies. Some previously described predictive models identified donor-specific risk factors for BOS (age ≥60 years, high PaO₂, smoking,
<figure><img src="image_6.png" /><figcaption>**FIGURE 3** | Probabilities of stability, bronchiolitis obliterans syndrome (BOS), and restrictive allograft syndrome (RAS) calculated for all modalities of each significant independent variable for the studied population (n = 230).</figcaption></figure>

pulmonary infection, and HLA mismatch) (16) or recipient-related acute events occurring during follow-up (acute cellular rejection, infections, and fungal pathogens) (18, 19) as risk factors for BOS. So far, no study described a prediction model specifically for early-onset BOS and RAS. In our study, although factors of the BOS and RAS models were the same, significant differences were observed within the models. Y1 class II DSAs were associated with RAS and BOS, but the OR were 6.97 and 3.83, respectively. Compared to CF, ILD/IPF was COPD were associated with RAS. In the univariate analysis, use of rATG (vs. basiliximab) was a statistically significant risk factor for BOS, whereas tacrolimus was protective for RAS. Although neither induction nor maintenance immunosuppression were statistically significant in the prediction model, they both participated in it.

To date, no definite association between BOS development and underlying diagnoses has been established, but a shorter

time to BOS has been reported for emphysema patients as compared to CF (31). Interestingly, IPF was associated with worse pulmonary function after BOS onset (32) and patients undergoing retransplantation for RAS were less likely to have CF and more likely to have IPF (33). In our study, as compared to CF, ILD/IPF (OR 5.467) and COPD (OR 3.857) were associated with RAS. To our knowledge, the later finding regarding COPD has not been previously described. Further studies are needed to confirm the underlying mechanisms of this association, but it can be hypothesized that pathological remodeling processes, as, for example, TGF-β signaling pathways or circulating fibroblasts, associated with ILD/IPF (34, 35) or COPD (36–39), may persist or be activated in some patients after LT and thus contribute to the higher risk of RAS.

Concerning recipients' age and CLAD, different cutoff points have been used with discordant results (40–43). Data from the
ISHLT Registry showed that 5-year incidence of BOS was 35.9% in pediatric and 41% adult patients (44, 45). In our study, the pediatric population was too small (3.5%) to be analyzed separately. A graphic representation of early-onset CLAD diagnosis according to recipient age showed a tendency for a U-shaped distribution, but this was not statistically significant. In the final model, increasing recipient age had a weak protective effect for CLAD, suggesting that age-related factors, such as worse adherence in children and adolescents or functional changes in the immune system during growth, may influence alloreactivity (46, 47). However, recipients' age was not a component of the best predictive model for BOS or RAS. Regarding donors' age, it was not associated with BOS/RAS or mortality; and, although literature is inconclusive (31, 48, 49), the observed lack of association in our cohort study is reassuring in regard to current donor selection criteria.

Evidence is not conclusive concerning induction treatment and the risk for CLAD (42, 43, 50–53), but recent data indicate that induction with basiliximab (or alemtuzumab) may be protective against BOS (54). In our study, use of basiliximab, as compared to rATG was also protective against BOS. It may be assumed that center effect or patient selection bias (e.g., high-risk patients receiving rATG) could account for the observed results. However, this factor remained significant even after controlling for center effect and independently of other perioperative risk factors (pre-LT DSAs, HLA mismatches, or cold ischemia time). Based on our results, causality cannot be established, but differences in the mechanisms of action of these molecules may explain our findings. The characteristics of post-depletion T cells and the susceptibility of individual T-cell subsets may vary between these agents (55–59). These results need to be further investigated as they concern a potentially modifiable factor.

Growing evidence suggests that DSAs are a risk factor for BOS and BOS-related mortality (31, 60, 61). DSAs that develop early after LT and persistent DSAs have been associated with worse outcomes (31). In our study, detection of type II DSAs at least once during year-1 was a risk factor of CLAD, BOS, and RAS consistently. As previously mentioned, important discrepancies existed in DSA measurements among different centers, reflecting a common problem of multicenter studies. In order to alleviate these discrepancies, we used the interpretation of the corresponding specialized laboratory rather than MFI cutoff points. Ongoing data collection in our cohort, specifically focusing on this parameter, is expected to provide additional information on this factor.

In line with evidence from transplantation of other solid organs, recent studies indicate that total HLA mismatches are associated with an increased risk for CLAD (54, 62). Hayes et al. described an increased risk for BOS, with HLA-A mismatches being associated with a greater hazard (54), whereas Walton et al. demonstrated the importance of eplet mismatches as a risk factor for RAS (62). In our study, HLA mismatches were an independent risk factor of mortality in the multivariate model but not for BOS or RAS. The observed discrepancy in these results may be due to differences in the definition of HLA mismatches and notably the non-inclusion of baseline DQ mismatches in our study since they were not available for all patients since the creation of the two cohorts.

Concerning other parameters, such as BMI (54), infections (63), or stage 3 PGD (64) and their association with BOS or RAS, data are limited. Clear associations may be hampered by the use of preemptive therapy or the decreased number of severe PGD survivors studied longitudinally. Acute cellular rejections (42, 65) and lymphocytic bronchiolitis (14, 66, 67) have been associated with BOS, but different definitions have been used in the literature. Concerning maintenance immunosuppression, two prospective studies showed a higher risk of BOS for cyclosporine vs. tacrolimus (68, 69). In our study, none of these clinical parameters had statistically significant associations with CLAD, with the exception of BMI and tacrolimus. BMI was associated with RAS in the univariate analysis but not in the multivariate analysis. For tacrolimus, a protective effect over RAS was shown in the univariate analysis, and this variable participated in the multivariate prediction model. Finally, for harmonization purposes, Y1-treated acute cellular rejections and lymphocytic bronchiolitis were studied together and Y1 t-AR was not an independent risk factor for BOS or RAS.

Our study has limitations. First, the choice of a 2-year minimum follow-up resulted in the exclusion of a large patient population from both cohorts. Although this approach may be associated with a patient selection bias, it allowed an unequivocal diagnosis in a number of patients sufficient for the construction of the models, while avoiding a bias associated with diagnostic uncertainty. Second, a specific time-point analysis was chosen over a time-dependent design to facilitate the subsequent multilevel clinical and laboratory data integration of the SysCLAD project. Follow-up was limited at 3 years because the main focus was early-onset CLAD. Although for RAS outcome seems to be inferior independently of the time-point of onset after LT, early BOS development has been associated with worse outcomes compared to late-onset BOS. Collectively, early-onset CLAD is associated with increased morbidity and mortality making the identification of associated risk factors even more clinically relevant. Finally, although there was not a significant number of missing values during model construction, very heterogeneous variables were not possible to harmonize (missing variables). These discrepancies reflect daily clinical practice and need to be taken into account when designing and analyzing multicenter databases. Adjudication of study patients and rigorous data harmonization are two of the main strengths of this project which differentiate it from registry studies handling more heterogeneous and difficult to control data (70). These results will need validation in a separate group of patients for whom data collection is ongoing within both cohorts.

In conclusion, the initial clinical data analysis of the SysCLAD cohort identified clinical factors as potential predictors of early-onset CLAD within 3 years post-LT. Among these factors, the underlying diagnosis, induction treatment, and the presence of Y1 class II DSAs were consistently associated with the development of CLAD and its main phenotypes, BOS or RAS. Validation of these results in subsequent patient populations and integration in the ongoing genetic, biological, and microbiological analyses of the SysCLAD project may assist patient risk stratification allowing a better understanding of the mechanisms propagating CLAD.
# SysCLAD CONSORTIUM

## Cohort of Lung Transplantation—COLT

**Bordeaux:** J. Jougon, J-F. Velly, H. Rozé, E. Blanchard, C. Dromer; **Bruxelles:** M. Antoine, M. Cappello, M. Ruiz, Y. Sokolow, F. Vanden Eynden, G. Van Nooten, L. Barvais, J. Berré, S. Brimioulle, D. De Backer, J. Créteur, E. Engelman, I. Huybrechts, B. Ickx, T. J. C. Preiser, T. Tuna, L. Van Obberghe, N. Vancutsem, J-L. Vincent, P. De Vuyst, I. Etienne, F. Féry, F. Jacobs, C. Knoop, J. L. Vachiéry, P. Van den Borne, I. Wellemans, G. Amand, L. Collignon, M. Giroux; **Grenoble:** E. Arnaud-Crozat, V. Bach, P-Y. Brichon, P. Chaffanjon, O. Chavanon, A. de Lambert, J-P. Fleury, S. Guigard, R. Hacini, K. Hireche, A. Pirvu, P. Porcu, P. Albaladejo, C. Allègre, A. Bataillard, D. Bedague, E. Briot, M. Casez-Brasseur, D. Colas, G. Dessertaine, M. Durand, G. Francony, A. Hebrard, M. R. Marino, B. Oummahan, D. Protar, D. Rehm, S. Robin, M. Rossi-Blancher, C. Augier, P. Bedouch, A. Boignard, H. Bouvaist, E. Brambilla, A. Briault, B. Camara, J. Claustre, S. Chanoine, M. Dubuc, S. Quétant, J. Maurizi, P. Pavèse, C. Pison, C. Saint-Raymond, N. Wion, C. Chérion; **Lyon:** R. Grima, O. Jegaden, J-M. Maury, F. Tronc, C. Flamens, S. Paulus, J-F. Mornex, F. Philit, A. Senechal, -C. Glérant, S. Turquier, D. Gamondes, L. Chalabresse, F. Thivolet-Bejui, C. Barnel, C. Dubois, A. Tiberghien; **Paris, Hôpital Européen Georges Pompidou:** F. Le Pimpec-Barthes, A. Bel, P. Mordant, P. Achouh, V. Boussaud, R. Guillemain, D. Méléard, M. O. Bricourt, B. Cholley, V. Pezella; **Marseille:** G. Brioude, X. B. D'Journo, C. Doddoli, P. Thomas, D. Trousse, S. Dizier, M. Leone, L. Papazian, F. Bregeon, A. Basire, B. Coltey, N. Dufeu, H. Dutau, S. Garcia, J. Y. Gaubert, C. Gomez, S. Laroumagne, A. Nieves, L. C. Picard, M. Reynaud-Gaubert, V. Secq, G. Mouton; **Nantes:** O. Baron, C. Brossaud, E. Durand, M. Durand, P. Lacoste, C. Perigaud, J. C. Roussel, I. Danner, A. Haloun, A. Magnan, A. Tissot, T. Lepoitvre, M. Treilhaud, K. Botturi-Cavaillès, S. Brouard, R. Danger, J. Loy, M. Morisset, M. Pain, S. Pares, D. Reboulleau, P.-J. Royer; **Le Plessis Robinson, Hôpital Marie Lannelongue:** Ph. Dartefelle, D. Fabre, E. Fadel, O. Mercier, S. Mussot, F. Stephan, P. Viard, J. Cerrina, P. Dorfmuller, S. Feuillet, M. Ghigna, Ph. Hervén F. Le Roy Ladurie, J. Le Pavec, V. Thomas de Montpreville, L. Lamrani; **Paris, Hôpital Bichat:** Y. Castier, P. Mordant, P. Cerceau, P. Augustin, S. Jean-Baptiste, S. Boudinet, P. Montravers, O. Brugière, G. Dauriat, G. Jébrak, H. Mal, A. Marceau, A-C. Métivier, G. Thabut, E. Lhuillier, C. Dupin, V. Bunel; **Strasbourg:** P. Falcoz, G. Massard, N. Santelmo, G. Ajob, O. Collange, O. Helms, J. Hentz, A. Roche, B. Bakouboula, T. Degot, A. Dory, S. Hirschi, S. Ohlmann-Caillard, L. Kessler, R. Kessler, A. Schuller, B. Renaud-Picard, K. Bennedif, S. Vargas; **Suresnes:** P. Bonnette, A. Chapelier, P. Puyo, E. Sage, J. Bresson, V. Caille, C. Cerf, J. Devaquet, V. Dumans-Nizard, M. L. Felten, M. Fischler, A. G. Si Larbi, M. Leguen, L. Ley, N. Liu, G. Trebbia, S. De Miranda, B. Douvry, F. Gonin, D. Grenet, A. M. Hamid, H. Neveu, F. Parquin, C. Picard, A. Roux, M. Stern, F. Bouillioud, P. Cahen, M. Colombat, C. Dautricourt, M. Delahousse, B. D'Urso, J. Gravisse, A. Guth, S. Hillaire, P. Honderlick, M. Lequintrec, E. Longchampt, F. Mellot, A. Scherrer, L. Temagoult, L. Tricot, M. Vasse, C. Veyrie, L. Zemoura; **Toulouse:** J. Berjaud, L. Brouchet, M. Dahan, F. O. Mathe,

H. Benahoua, M. DaCosta, I. Serres, V. Merlet-Dupuy, M. Grigoli, A. Didier, M. Murris, L. Crognier, O. Fourcade.

## Swiss Lung Transplant Centers

**Lausanne-Geneva:** T. Krueger, H. B. Ris, M. Gonzalez, J.-D. Aubert, L. P. Nicod, B. J. Marsland, C. Berutto, T. Rochat, P. Soccal, Ph. Jolliet, A. Koutsokera, C. Marcucci, O. Manuel, E. Bernasconi, M. Chollet, F. Gronchi, C. Courbon; **Zurich:** S. Hillinger, I. Inci, P. Kestenholz, W. Weder; R. Schuepbach, M. Zalunardo, C. Benden, U. Buergi, L. C. Huber, B. Isenring, M. M. Schuurmans, A. Gaspert, D. Holzmann, N. Müller, T. Rechsteiner, C. Schmid, B. Vrugt.

We would like to thank the following collaborators of the Swiss Lung Transplantation Centers for their contribution in data collection and/or project coordination (alphabetical order): E. Catana, C. Cowaloosur-Noirat, M. F. Derkenne, JL Dreifuss, P. Grendelmeier, J. Hartwig, N. Lourenco, M. Magno, H. Muller-McKenna, E. Perret, and K. Zangger.

## Swiss Transplant Cohort Study—STCS

The members of the Swiss Transplant Cohort Study are: Rita Achermann, Patrizia Amico, John-David Aubert, Philippe Baumann, Guido Beldi, Christian Benden, Christoph Berger, Isabelle Binet, Pierre-Yves Bochud, Elsa Boely, Heiner Bucher, Leo Bühler, Thierry Carell, Emmanuelle Catana, Yves Chalandon, Sabina de Geest, Olivier de Rougemont, Michael Dickenmann, Michel Duchosal, Laure Elkrief, Thomas Fehr, Sylvie Ferrari-Lacraz, Christian Garzoni, Paola Gasche Soccal, Christophe Gaudet, Emiliano Giostra, Délà Golshayan, Karine Hadaya, Jörg Halter, Dominik Heim, Christoph Hess, Sven Hillinger, Hans H. Hirsch, Günther Hofbauer, Uyen Huynh-Do, Franz Immer, Richard Klaghofer, Michael Koller (Head of the data center), Bettina Laesser, Roger Lehmann, Christian Lovis, Oriol Manuel, Hans-Peter Marti, Pierre Yves Martin, Luca Martinolli, Pascal Meylan, (Head, Biological samples management group), Paul Mohacsi, Philippe Morel, Ulrike Mueller, Nicolas J. Mueller (Chairman Scientific Committee), Helen Mueller-McKenna (Head of local data management), Antonia Müller, Thomas Müller, Beat Müllhaupt, David Nadal, Manuel Pascual (Executive office), Jakob Passweg, Juliane Rick, Eddy Roosnek, Anne Rosselet, Silvia Rothlin, Frank Ruschitzka, Urs Schanz, Stefan Schaub, Aurelia Schnyder, Christian Seiler, Susanne Stampf, Jürg Steiger (Head, Executive Office), Guido Stirnimann, Christian Toso, Christian Van Delden (Executive office), Jean-Pierre Venetz, Jean Villard, Madeleine Wick (STCS coordinator), Markus Wilhelm, and Patrick Yerly.

## SME and Platforms

**Biomax (Munich, Germany):** A. Fritz, D. Maier; **Finovatis (Lyon, France):** K. Desplanche, D. Koubi; **GATC (Germany):** F. Ernst, T. Paprotka, M. Schmitt, B. Wahl; **Novasdcovery (Lyon, France):** J.-P. Boissel, G. Olivera-Botello; **Prométhée Proteomics Platform (Grenoble, France):** C. Trocmé, B. Toussaint, S. Bourgoin-Voillard, M. Séve; **Inserm U823, Université Joseph Fourier (Grenoble, France):** M. Benmerad, V. Siroux, R. Slama; **European Institute for Systems Biology & Medicine (Lyon, France):** C. Auffray, D. Charron, J. Pellet, C. Pison.
## ETHICS STATEMENT

This study was carried out in accordance with the recommendations of national and local ethics committees of the participating cohorts with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the national and local ethics committees of the participating hospitals.

## AUTHOR CONTRIBUTIONS

AK, PR, AF, AT, and KB performed data collection and harmonization. CB, JDA, AR, MR-G, AM, CP, and LN participated in the adjudication committee. JPA and AK performed the statistical analysis. AK drafted the manuscript. CB, JDA, AR, MR-G, RK, CD, SM, HM, J-FM, RG, CK, MD, PS, JC, ES, CG, AM, CP, and LN included the patients and revised the manuscript for important intellectual content.

## FUNDING

The SysCLAD study is an EU-funded project, HEALTH-F5-2012 (grant agreement #305457) under the Seventh Framework Programme (FP7). The authors are indebted to the “Programme Hospitalier de Recherche Clinique 2008,” to “Vaincre la

## REFERENCES

1. Verleden GM, Raghu G, Meyer KC, Glanville AR, Corris P. A new classification system for chronic lung allograft dysfunction. *J Heart Lung Transplant* (2014) 33(2):127-33. doi:10.1016/j.healun.2013.10.022

2. Sato M. Chronic lung allograft dysfunction after lung transplantation: the moving target. *Gen Thorac Cardiovasc Surg* (2013) 61(2):67-78. doi:10.1007/s11748-012-0167-3

3. Verleden GM, Vos R, Vanaudenaerde B, Dupont L, Yserbyt J, Van Raemdonck D, et al. Current views on chronic rejection after lung transplantation. *Transpl Int* (2015) 28(10):1131-9. doi:10.1111/tri.12579

4. Meyer KC, Raghu G, Verleden GM, Corris PA, Aurora P, Wilson KC, et al. An international ISHLT/ATS/ERS clinical practice guideline: diagnosis and management of bronchiolitis obliterans syndrome. *Eur Respir J* (2014) 44(6):1479-503. doi:10.1183/09031936.00107514

5. Jackson CH, Sharples LD, McNeil K, Stewart S, Wallwork J. Acute and chronic onset of bronchiolitis obliterans syndrome (BOS): are they different entities? *J Heart Lung Transplant* (2002) 21(6):658-66. doi:10.1016/S1053-2498(02)00381-9

6. Brugiere O, Pessione F, Thabut G, Mal H, Jebrak G, Leseche G, et al. Bronchiolitis obliterans syndrome after single-lung transplantation: impact of time to onset on functional pattern and survival. *Chest* (2002) 121(6):1883-9. doi:10.1378/chest.121.6.1883

7. Verleden SE, Vandermeulen E, Ruttens D, Vos R, Vaneylen A, Dupont LJ, et al. Neutrophilic reversible allograft dysfunction (NRAD) and restrictive allograft syndrome (RAS). *Semin Respir Crit Care Med* (2013) 34(3):352-60. doi:10.1055/s-0033-1348463

8. Sato M, Waddell TK, Wagnetz U, Roberts HC, Hwang DM, Haroon A, et al. Restrictive allograft syndrome (RAS): a novel form of chronic lung allograft dysfunction. *J Heart Lung Transplant* (2011) 30(7):735-42. doi:10.1016/j.healun.2011.01.712

9. Suhling H, Dettmer S, Greer M, Fuehner T, Avsar M, Haverich A, et al. Phenotyping chronic lung allograft dysfunction using body plethysmography and computed tomography. *Am J Transplant* (2016) 16(11):3163-70. doi:10.1111/ajt.13876

10. Verleden SE, Ruttens D, Vandermeulen E, Bellon H, Dubbeldam A, De Wever W, et al. Predictors of survival in restrictive chronic lung allograft dysfunction

Mucoviscidose" and Association Grégory Lemarchal for supporting this project from its beginning, when 11 French lung transplantation centers gave rise to the Cohort in Lung Transplantation (COLT), "Programme Transplantation 2008," PRTP-13, http://ClinicalTrials.gov Identifier: NCT00980967. This work was realized in the context of the IHU-Cesti project thanks to the French government financial support managed by the National Research Agency *via* the "Investment into the Future" program ANR-10-IBHU-005. The IHU-Cesti project is also supported by Nantes Métropole and Région Pays de la Loire. This study has been conducted in the framework of the Swiss Transplant Cohort Study, supported by the Swiss National Science Foundation and the Swiss University Hospitals (G15) and transplant centers. The authors thank the Swiss National Research Foundation for supporting the STCS and in particular the lung transplant section (No 3347CO-108795) and the Juchum foundation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at http://journal.frontiersin.org/article/10.3389/fmed.2017.00109/full#supplementary-material.

after lung transplantation. *J Heart Lung Transplant* (2016) 35(9):1078-84. doi:10.1016/j.healun.2016.03.022

11. Royer PJ, Olivera-Botello G, Koutsokera A, Aubert JD, Bernasconi E, Tissot A, et al. Chronic lung allograft dysfunction: a systematic review of mechanisms. *Transplantation* (2016) 100(9):1803-14. doi:10.1097/TP.0000000000001215

12. Sharples LD, McNeil K, Stewart S, Wallwork J. Risk factors for bronchiolitis obliterans: a systematic review of recent publications. *J Heart Lung Transplant* (2002) 21(2):271-81. doi:10.1016/S1053-2498(01)00360-6

13. Meyer K. Bronchiolitis obliterans syndrome and chronic lung allograft dysfunction: evolving concepts and nomenclature. In: Meyer KC, Glanville AR, editors. *Bronchiolitis Obliterans Syndrome in Lung Transplantation*. Rounds, SI: Springer Science and Business Media (2013). p. 1-19.

14. Verleden SE, Ruttens D, Vandermeulen E, Vaneylen A, Dupont LJ, Van Raemdonck DE, et al. Bronchiolitis obliterans syndrome and restrictive allograft syndrome: do risk factors differ? *Transplantation* (2013) 95(9): 1167-72. doi:10.1097/TP.0b013e318286e076

15. Collins GS, Reitsma JB, Altman DG, Moons KG. Transparent Reporting of a multivariable prediction model for Individual Prognosis or Diagnosis (TRIPOD): the TRIPOD statement. *Ann Intern Med* (2015) 162(1):55-63. doi:10.7326/M14-0697

16. Hennessy SA, Hranjec T, Swenson BR, Kozower BD, Jones DR, Ailawadi G, et al. Donor factors are associated with bronchiolitis obliterans syndrome after lung transplantation. *Ann Thorac Surg* (2010) 89(5):1555-62. doi:10.1016/j.athoracsur.2010.01.060

17. Jaksch P, Taghavi S, Klepetko W, Salama M. Pretransplant serum human chitinase-like glycoprotein YKL-40 concentrations independently predict bronchiolitis obliterans development in lung transplant recipients. *J Thorac Cardiovasc Surg* (2014) 148(1):273-81. doi:10.1016/j.jtcvs.2014.02.059

18. Valentine VG, Gupta MR, Walker JE Jr, Seoane L, Bonvillain RW, Lombard GA, et al. Effect of etiology and timing of respiratory tract infections on development of bronchiolitis obliterans syndrome. *J Heart Lung Transplant* (2009) 28(2):163-9. doi:10.1016/j.healun.2008.11.907

19. Weigt SS, Elashoff RM, Huang C, Ardehali A, Gregson AL, Kubak B, et al. *Aspergillus* colonization of the lung allograft is a risk factor for bronchiolitis obliterans syndrome. *Am J Transplant* (2009) 9(8):1903-11. doi:10.1111/j.1600-6143.2009.02635.x
20. Pison C, Magnan A, Botturi K, Seve M, Brouard S, Marsland BJ, et al. Prediction of chronic lung allograft dysfunction: a systems medicine challenge. *Eur Respir J* (2014) 43(3):689-93. doi:10.1183/09031936.00161313

21. Koutsokera A, Royer PJ, Fritz A, Benden C, Tissot A, Aubert JD, et al. Risk factors for chronic lung allograft dysfunction (CLAD) in the SysCLAD cohort (conference abstract). *Eur Respir J* (2015) 46:A1800. doi:10.1183/13993003. congress-2015.PA1800

22. Pison C, Koutsokera A, Tissot A, Botturi K, Benden C, Aubert JD, et al. Cohort characteristics and chronic lung allograft dysfunction (CLAD) adjudication in a FP-7 funded project, systems prediction of CLAD (SysCLAD) (conference abstract). *Eur Respir J* (2014) 44(S58):1429.

23. Berger C, Bochud PY, Boggian K, Cusini A, Egli A, Garzoni C, et al. The Swiss Transplant Cohort Study: lessons from the first 6 years. *Curr Infect Dis Rep* (2015) 17(6):486. doi:10.1007/s11908-015-0486-5

24. Koller MT, van Delden C, Muller NJ, Baumann P, Lovis C, Marti HP, et al. Design and methodology of the Swiss Transplant Cohort Study (STCS): a comprehensive prospective nationwide long-term follow-up cohort. *Eur J Epidemiol* (2013) 28(4):347-55. doi:10.1007/s10654-012-9754-y

25. Verleden SE, Ruttens D, Vandermeulen E, Bellon H, Van Raemdonck DE, Dupont LJ, et al. Restrictive chronic lung allograft dysfunction: where are we now? *J Heart Lung Transplant* (2015) 34(5):625-30. doi:10.1016/j. healun.2014.11.007

26. Christie JD, Carby M, Bag R, Corris P, Hertz M, Weill D, et al. Report of the ISHLT working group on primary lung graft dysfunction part II: definition. A consensus statement of the International Society for Heart and Lung Transplantation. *J Heart Lung Transplant* (2005) 24(10):1454-9. doi:10.1016/j. healun.2004.11.049

27. Christie JD, Bavaria JE, Palevsky HI, Litzky L, Blumenthal NP, Kaiser LR, et al. Primary graft failure following lung transplantation. *Chest* (1998) 114(1):51-60. doi:10.1378/chest.114.1.51

28. Levine DJ, Glanville AR, Aboeyoun C, Belperio J, Benden C, Berry GJ, et al. Antibody-mediated rejection of the lung: a consensus report of the International Society for Heart and Lung Transplantation. *J Heart Lung Transplant* (2016) 35(4):397-406. doi:10.1016/j.healun.2016.01.1223

29. Hosmer DL, Lemeshow S, Sturdivant RX. *Applied Logistic Regression*. New Jersey: John Wiley and Sons (2000). p. 116-28.

30. Al-Ashkar F, Mehra R, Mazzone PJ. Interpreting pulmonary function tests: recognize the pattern, and the diagnosis will follow. *Cleve Clin J Med* (2003) 70(10): 866, 868, 871-3. doi:10.3949/ccjm.70.10.866

31. Safavi S, Robinson DR, Soresi S, Carby M, Smith JD. De novo donor HLA-specific antibodies predict development of bronchiolitis obliterans syndrome after lung transplantation. *J Heart Lung Transplant* (2014) 33(12):1273-81. doi:10.1016/j.healun.2014.07.012

32. Lama VN, Murray S, Lonigro RJ, Toews GB, Chang A, Lau C, et al. Course of FEV1(1) after onset of bronchiolitis obliterans syndrome in lung transplant recipients. *Am J Respir Crit Care Med* (2007) 175(11):1192-8. doi:10.1164/rccm.200609-1344OC

33. Verleden SE, Todd JL, Sato M, Palmer SM, Martinu T, Pavlisko EN, et al. Impact of CLAD phenotype on survival after lung retransplantation: a multicenter study. *Am J Transplant* (2015) 15(8):2223-30. doi:10.1111/ajt.13281.

34. LaPar DJ, Burdick MD, Emaminia A, Harris DA, Strieter BA, Liu L, et al. Circulating fibrocytes correlate with bronchiolitis obliterans syndrome development after lung transplantation: a novel clinical biomarker. *Ann Thorac Surg* (2011) 92(2):470-7; discussion 7. doi:10.1016/j.athoracsur.2011.04.065

35. Strieter RM, Keeley EC, Hughes MA, Burdick MD, Mehrad B. The role of circulating mesenchymal progenitor cells (fibrocytes) in the pathogenesis of pulmonary fibrosis. *J Leukoc Biol* (2009) 86(5):1111-8. doi:10.1189/jlb.0309132

36. Dupin I, Allard B, Ozier A, Maurat E, Ousova O, Delbrel E, et al. Blood fibrocytes are recruited during acute exacerbations of chronic obstructive pulmonary disease through a CXCR4-dependent pathway. *J Allergy Clin Immunol* (2016) 137(4): 1036-42.e1-7. doi:10.1016/j.jaci.2015.08.043

37. Wong SL, Sukkar MB. SPARC in the lung: an overview of its role in lung cancer and pulmonary fibrosis and its potential role in chronic airways disease. *Br J Pharmacol* (2017) 174(1):3-14. doi:10.1111/bph.13653

38. Mahmood MQ, Reid D, Ward C, Muller HK, Knight DA, Sohal SS, et al. Transforming growth factor (TGF) beta1 and Smad signalling pathways: a likely key to EMT-associated COPD pathogenesis. *Respirology* (2017) 22(1):133-140. doi:10.1111/resp.12882

39. Kurundkar AR, Kurundkar D, Rangarajan S, Locy ML, Zhou Y, Liu RM, et al. The matricellular protein CCN1 enhances TGF-beta1/SMAD3-dependent profibrotic signaling in fibroblasts and contributes to fibrogenic responses to lung injury. *FASEB J* (2016) 30(6):2135-50. doi:10.1096/fj.201500173

40. Weiss ES, Allen JG, Modi MN, Merlo CA, Conte JV, Shah AS. Lung transplantation in older patients with cystic fibrosis: analysis of UNOS data. *J Heart Lung Transplant* (2009) 28(2):135-40. doi:10.1016/j.healun.2008.11.903

41. Mahidhara R, Bastani S, Ross DJ, Saggar R, Lynch J III, Schnickel GT, et al. Lung transplantation in older patients? *J Thorac Cardiovasc Surg* (2008) 135(2):412-20. doi:10.1016/j.jtcvs.2007.09.030

42. Hachem RR, Khalifah AP, Chakinala MM, Yusen RD, Aloush AA, Mohanakumar T, et al. The significance of a single episode of minimal acute rejection after lung transplantation. *Transplantation* (2005) 80(10):1406-13. doi:10.1097/01.tp.0000181161.60638.fa

43. Khalifah AP, Hachem RR, Chakinala MM, Yusen RD, Aloush A, Patterson GA, et al. Minimal acute rejection after lung transplantation: a risk for bronchiolitis obliterans syndrome. *Am J Transplant* (2005) 5(8):2022-30. doi:10.1111/j.1600-6143.2005.00953.x

44. Yusen RD, Edwards LB, Kucheryavaya AY, Benden C, Dipchand AI, Goldfarb SB, et al. The registry of the international society for heart and lung transplantation: thirty-second official adult lung and heart-lung transplantation report – 2015; focus theme: early graft failure. *J Heart Lung Transplant* (2015) 34(10):1264-77. doi:10.1016/j.healun.2015.08.014

45. Goldfarb SB, Benden C, Edwards LB, Kucheryavaya AY, Dipchand AI, Levvey BJ, et al. The registry of the international society for heart and lung transplantation: eighteenth official pediatric lung and heart-lung transplantation report – 2015; focus theme: early graft failure. *J Heart Lung Transplant* (2015) 34(10):1255-63. doi:10.1016/j.healun.2015.08.005

46. Schmid FA, Benden C. Special considerations for the use of lung transplantation in pediatrics. *Expert Rev Respir Med* (2016) 10(6):655-62. doi:10.1586/17476348.2016.1168298

47. Denecke C, Bedi DS, Ge X, Kim IK, Jurisch A, Weiland A, et al. Prolonged graft survival in older recipient mice is determined by impaired effector T-cell but intact regulatory T-cell responses. *PLoS One* (2010) 5(2):e9232. doi:10.1371/journal.pone.0009232

48. De Perrot M, Waddell TK, Shargall Y, Pierre AF, Fadel E, Uy K, et al. Impact of donors aged 60 years or more on outcome after lung transplantation: results of an 11-year single-center experience. *J Thorac Cardiovasc Surg* (2007) 133(2):525-31. doi:10.1016/j.jtcvs.2006.09.054

49. Pizanis N, Heckmann J, Tsagakis K, Tossios P, Massoudy P, Wendt D, et al. Lung transplantation using donors 55 years and older: is it safe or just a way out of organ shortage? *Eur J Cardiothorac Surg* (2010) 38(2):192-7. doi:10.1016/j.ejcts.2010.01.054

50. Penninga L, Moller CH, Penninga EI, Iversen M, Gluud C, Steinbruchel DA. Antibody induction therapy for lung transplant recipients. *Cochrane Database Syst Rev* (2013) 11:CD008927. doi:10.1002/14618588.CD008927.pub2

51. Borro JM, De la Torre M, Miguelez C, Fernandez R, Gonzalez D, Lemos C. Comparative study of basiliximab treatment in lung transplantation. *Transplant Proc* (2005) 37(9):3996-8. doi:10.1016/j.transproceed.2005.09.192

52. de la Torre M, Pena E, Calvin M, Miguelez C, Gonzalez D, Fernandez R, et al. Basiliximab in lung transplantation: preliminary experience. *Transplant Proc* (2005) 37(3):1534-6. doi:10.1016/j.transproceed.2005.02.026

53. Hachem RR, Chakinala MM, Yusen RD, Lynch JP, Aloush AA, Patterson GA, et al. A comparison of basiliximab and anti-thymocyte globulin as induction agents after lung transplantation. *J Heart Lung Transplant* (2005) 24(9):1320-6. doi:10.1016/j.healun.2004.09.002

54. Hayes D Jr, Auletta JJ, Whitson BA, Black SM, Kirkby S, Tobias JD, et al. Human leukocyte antigen mismatching and survival after lung transplantation in adult and pediatric patients with cystic fibrosis. *J Thorac Cardiovasc Surg* (2016) 151(2):549-57.e1. doi:10.1016/j.jtcvs.2015.08.022

55. Pearl JP, Parris J, Hale DA, Hoffmann SC, Bernstein WB, McCoy KL, et al. Immunocompetent T-cells with a memory-like phenotype are the dominant cell type following antibody-mediated T-cell depletion. *Am J Transplant* (2005) 5(3):465-74. doi:10.1111/j.1600-6143.2005.00759.x

56. Gurkan S, Luan Y, Dhillon N, Allam SR, Montague T, Bromberg JS, et al. Immune reconstitution following rabbit antithymocyte globulin. *Am J Transplant* (2010) 10(9):2132-41. doi:10.1111/j.1600-6143.2010.03210.x

57. LaCorcia G, Swistak M, Lawendowski C, Duan S, Weeden T, Nahill S, et al. Cyclonal rabbit antithymocyte globulin exhibits consistent
immunosuppressive capabilities beyond cell depletion. *Transplantation* (2009) 87(7):966-74. doi:10.1097/TP.0b013e31819c84b8

58. Leitner J, Grabmeier-Pfistershammer K, Majdic O, Zlabinger G, Steinberger P. Interaction of antithymocyte globulins with dendritic cell antigens. *Am J Transplant* (2011) 11(1):138-45. doi:10.1111/j.1600-6143.2010.03322.x

59. Lopez M, Clarkson MR, Albin M, Sayegh MH, Najafian N. A novel mechanism of action for anti-thymocyte globulin: induction of CD4+CD25+Foxp3+ regulatory T cells. *J Am Soc Nephrol* (2006) 17(10):2844-53. doi:10.1681/ASN.2006050422

60. Morrell MR, Pilewski JM, Gries CJ, Pipeling MR, Crespo MM, Ensor CR, et al. De novo donor-specific HLA antibodies are associated with early and high-grade bronchiolitis obliterans syndrome and death after lung transplantation. *J Heart Lung Transplant* (2014) 33(12):1288-94. doi:10.1016/j.healun.2014.07.018

61. Kauke T, Kneidinger N, Martin B, Dick A, Schneider C, Schramm R, et al. Bronchiolitis obliterans syndrome due to donor-specific HLA-antibodies. *Tissue Antigens* (2015) 86(3):178-85. doi:10.1111/tan.12626

62. Walton DC, Hiho SJ, Cantwell LS, Diviney MB, Wright ST, Snell GI, et al. HLA matching at the eplet level protects against chronic lung allograft dysfunction. *Am J Transplant* (2016) 16(9):2695-703. doi:10.1111/ajt.13798

63. Avery R. The role of infections in BOS. In: Meyer KC, Glanville AR, editors. Bronchiolitis obliterans syndrome in lung transplantation. *Respir Med* (2013) 8:199-218.

64. Kreisel D, Krupnick AS, Puri V, Guthrie TJ, Trulock EP, Meyers BF, et al. Short- and long-term outcomes of 1000 adult lung transplant recipients at a single center. *J Thorac Cardiovasc Surg* (2011) 141(1):215-22. doi:10.1016/j.jtcvs.2010.09.009

65. Hadjiliadis D, Davis RD, Palmer SM. Is transplant operation important in determining posttransplant risk of bronchiolitis obliterans syndrome in lung transplant recipients? *Chest* (2002) 122(4):1168-75. doi:10.1378/chest.122.4.1168

66. Burton CM, Iversen M, Scheike T, Carlsen J, Andersen CB. Is lymphocytic bronchiolitis a marker of acute rejection? An analysis of 2,697

transbronchial biopsies after lung transplantation. *J Heart Lung Transplant* (2008) 27(10):1128-34. doi:10.1016/j.healun.2008.06.014

67. Glanville AR, Aboyoun CL, Havryk A, Plit M, Rainer S, Malouf MA. Severity of lymphocytic bronchiolitis predicts long-term outcome after lung transplantation. *Am J Respir Crit Care Med* (2008) 177(9):1033-40. doi:10.1164/rccm.200706-951OC

68. Keenan RJ, Konishi H, Kawai A, Paradis IL, Nunley DR, Iacono AT, et al. Clinical trial of tacrolimus versus cyclosporine in lung transplantation. *Ann Thorac Surg* (1995) 60(3):580-4; discussion 4-5. doi:10.1016/0003-4975(95)00407-C

69. Treede H, Glanville AR, Klepetko W, Aboyoun C, Vettorazzi E, Lama R, et al. Tacrolimus and cyclosporine have differential effects on the risk of development of bronchiolitis obliterans syndrome: results of a prospective, randomized international trial in lung transplantation. *J Heart Lung Transplant* (2012) 31(8):797-804. doi:10.1016/j.healun.2012.03.008

70. Urschel S. Apples, oranges, and statistical magic: limitations of registry studies and need for collaborative studies. *J Heart Lung Transplant* (2015) 34(9):1136-8. doi:10.1016/j.healun.2015.05.013

**Conflict of Interest Statement:** CP received fees for symposia and financial support to attend medical meetings from Astellas, Novartis, Sanofi-Genzyme, and Therakos. AR has served as a consultant for Novartis France concerning CMV in solid organ transplantation. J-FM acknowledges consultancy, research grant and support for meeting attendance from LFB Biomedicament, CSL Behring, Actelion, Pfizer, Pioproject, and Bayer. The authors have declared that no competing interests exist.

Copyright © 2017 Koutsokera, Royer, Antonietti, Fritz, Benden, Aubert, Tissot, Botturi, Roux, Reynaud-Gaubert, Kessler, Dromer, Mussot, Mal, Mornex, Guillemain, Knoop, Dahan, Soccal, Claustre, Sage, Gomez, Magnan, Pison and Nicod. 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) or licensor 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.