Skip to main content

Main menu

  • Home
  • Current issue
  • ERJ Early View
  • Past issues
  • Authors/reviewers
    • Instructions for authors
    • Submit a manuscript
    • Open access
    • COVID-19 submission information
    • Peer reviewer login
  • Alerts
  • Podcasts
  • Subscriptions
  • ERS Publications
    • European Respiratory Journal
    • ERJ Open Research
    • European Respiratory Review
    • Breathe
    • ERS Books
    • ERS publications home

User menu

  • Log in
  • Subscribe
  • Contact Us
  • My Cart

Search

  • Advanced search
  • ERS Publications
    • European Respiratory Journal
    • ERJ Open Research
    • European Respiratory Review
    • Breathe
    • ERS Books
    • ERS publications home

Login

European Respiratory Society

Advanced Search

  • Home
  • Current issue
  • ERJ Early View
  • Past issues
  • Authors/reviewers
    • Instructions for authors
    • Submit a manuscript
    • Open access
    • COVID-19 submission information
    • Peer reviewer login
  • Alerts
  • Podcasts
  • Subscriptions

Transcriptomic gene signatures associated with persistent airflow limitation in patients with severe asthma

Pieter-Paul Hekking, Matthew J. Loza, Stelios Pavlidis, Bertrand De Meulder, Diane Lefaudeux, Frederic Baribaud, Charles Auffray, Ariane H. Wagener, Paul Brinkman, René Lutter, Aruna T. Bansal, Ana R. Sousa, Stewart A. Bates, Ioannis Pandis, Louise J. Fleming, Dominick E. Shaw, Stephen J. Fowler, Yike Guo, Andrea Meiser, Kai Sun, Julie Corfield, Peter Howarth, Elisabeth H. Bel, Ian M. Adcock, Kian Fan Chung, Ratko Djukanovic, Peter J. Sterk the U-BIOPRED Study Group
European Respiratory Journal 2017 50: 1602298; DOI: 10.1183/13993003.02298-2016
Pieter-Paul Hekking
1Dept of Respiratory Medicine, Academic Medical Center, Amsterdam, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: p.w.hekking@amc.uva.nl
Matthew J. Loza
2Janssen R&D, Johnson & Johnson, Springhouse, PA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stelios Pavlidis
3Data Science Institute, Imperial College London, London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bertrand De Meulder
4European Institute for Systems Biology and Medicine, CIRI UMR5308, CNRS-ENS-UCBL-INSERM, Université de Lyon, Lyon, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Diane Lefaudeux
4European Institute for Systems Biology and Medicine, CIRI UMR5308, CNRS-ENS-UCBL-INSERM, Université de Lyon, Lyon, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Frederic Baribaud
2Janssen R&D, Johnson & Johnson, Springhouse, PA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Charles Auffray
4European Institute for Systems Biology and Medicine, CIRI UMR5308, CNRS-ENS-UCBL-INSERM, Université de Lyon, Lyon, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ariane H. Wagener
1Dept of Respiratory Medicine, Academic Medical Center, Amsterdam, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paul Brinkman
1Dept of Respiratory Medicine, Academic Medical Center, Amsterdam, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
René Lutter
1Dept of Respiratory Medicine, Academic Medical Center, Amsterdam, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Aruna T. Bansal
5Acclarogen, Cambridge, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ana R. Sousa
6Discovery Medicine, GlaxoSmithKline, Brentford, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stewart A. Bates
6Discovery Medicine, GlaxoSmithKline, Brentford, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ioannis Pandis
3Data Science Institute, Imperial College London, London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Louise J. Fleming
3Data Science Institute, Imperial College London, London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Dominick E. Shaw
7Respiratory Research Unit, University of Nottingham, Nottingham, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stephen J. Fowler
8Centre for Respiratory Medicine and Allergy, Institute of Inflammation and Repair, University of Manchester and University Hospital of South Manchester, Manchester Academic Health Sciences Centre, Manchester, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yike Guo
3Data Science Institute, Imperial College London, London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Andrea Meiser
3Data Science Institute, Imperial College London, London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kai Sun
3Data Science Institute, Imperial College London, London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Julie Corfield
9Areteva, Nottingham, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Peter Howarth
10NIHR Southampton Centre for Biomedical Research, University of Southampton, Southampton, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Elisabeth H. Bel
1Dept of Respiratory Medicine, Academic Medical Center, Amsterdam, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ian M. Adcock
11National Heart and Lung Institute, Imperial College London, London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Kian Fan Chung
11National Heart and Lung Institute, Imperial College London, London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ratko Djukanovic
10NIHR Southampton Centre for Biomedical Research, University of Southampton, Southampton, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Peter J. Sterk
1Dept of Respiratory Medicine, Academic Medical Center, Amsterdam, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Abstract

A proportion of severe asthma patients suffers from persistent airflow limitation (PAL), often associated with more symptoms and exacerbations. Little is known about the underlying mechanisms. Here, our aim was to discover unexplored potential mechanisms using Gene Set Variation Analysis (GSVA), a sensitive technique that can detect underlying pathways in heterogeneous samples.

Severe asthma patients from the U-BIOPRED cohort with PAL (post-bronchodilator forced expiratory volume in 1 s/forced vital capacity ratio below the lower limit of normal) were compared with those without PAL. Gene expression was assessed on the total RNA of sputum cells, nasal brushings, and endobronchial brushings and biopsies. GSVA was applied to identify differentially enriched predefined gene signatures based on all available gene expression publications and data on airways disease.

Differentially enriched gene signatures were identified in nasal brushings (n=1), sputum (n=9), bronchial brushings (n=1) and bronchial biopsies (n=4) that were associated with response to inhaled steroids, eosinophils, interleukin-13, interferon-α, specific CD4+ T-cells and airway remodelling.

PAL in severe asthma has distinguishable underlying gene networks that are associated with treatment, inflammatory pathways and airway remodelling. These findings point towards targets for the therapy of PAL in severe asthma.

Abstract

Persistent airflow limitation in severe asthma is associated with a mechanism involving IL-13 and remodelling http://ow.ly/JYcC30daSRf

Introduction

A small proportion of asthma patients suffers from uncontrolled disease despite treatment with high-dose inhaled or oral steroids, which is defined as severe asthma [1, 2]. These patients are in need of effective and steroid-sparing treatments. However, severe asthma is a heterogeneous disease, which complicates the development of these treatments [3, 4]. Emerging targeted treatments appear to be successful in selected subgroups of patients with specific characteristics, such as eosinophilic airway inflammation and allergy, emphasising the need for precision medicine in severe asthma [5–7].

Agusti et al. [8] recently described a novel strategy for the treatment of chronic airways disease, focusing on treatable traits based on clinical, physiological and/or biological characteristics rather than solely on a diagnosis of asthma or chronic obstructive pulmonary disease (COPD). Persistent airflow limitation (PAL) was postulated as a treatable trait of airways disease, being characterised by reduced expiratory flow even after the administration of a short-acting bronchodilator. Notably, PAL is not only a determinant for the diagnosis of COPD, but regardless of smoking history, has also long been implicated in patients with severe asthma [9].

PAL is clinically relevant in severe asthma as it may have prognostic implications. Even though PAL and exacerbations are different phenotypic characteristics, there is evidence of an association of PAL with more frequent exacerbations and more hospitalisations in subgroups of severe asthma patients [10, 11]. Clinical clustering studies showed that PAL is an important constituent for severe asthma phenotypes that are characterised by eosinophilic airway inflammation, poor quality of life and smoking history [11–14]. To understand this clinical trait, several studies have been conducted to identify risk factors for PAL in smoking and nonsmoking severe asthmatic subjects. Eosinophilic airway inflammation, longer duration of disease, adult onset of disease and a history of smoking appear to be the most important factors [9, 15, 16]. However, management strategies for the prevention or treatment of PAL in severe asthma are lacking, which is unsurprising given the lack of knowledge regarding its mechanism.

Airway remodelling, associated with chronic inflammation, may contribute to PAL in asthma. Thickening of the airway smooth muscle layer, changes in extracellular matrix such as subepithelial fibrosis and infiltration of inflammatory cells in the bronchial airway wall have been linked to increased lung function decline [15, 17, 18]. Even though studies have been conducted to identify underlying pathways of airway remodelling, little is known about the underlying mechanism of PAL in severe asthma [19, 20]. Notably, there are differences in airway remodelling between patients with severe and patients with nonsevere asthma, demonstrating the need to examine the underlying pathways of PAL in patients with severe asthma themselves [21].

Gene Set Variation Analysis (GSVA) is a sensitive statistical technique that can be used for identifying underlying pathways in heterogeneous samples [22]. Instead of more conventional analysis of single-gene expression, GSVA allows between-groups comparisons of the expression of gene sets that are known to be associated with biological networks of relevance. This identifies differences in underlying pathways in severe asthma patients with PAL compared with those without PAL.

We hypothesised that patients with severe asthma and PAL exhibit different gene expression profiles compared with those without PAL. The aim of the present study was to identify transcriptomic pathways of PAL in severe asthma by applying GSVA in sputum cells, nasal brushings, and endobronchial brushings and biopsies. This study was conducted in the severe asthma adult cohort of the U-BIOPRED study, which according to a recent cluster analysis based on clinical and physiological parameters includes subsets of patients with and without PAL [23].

Methods

Detailed methods are provided in the supplementary material.

Design and subjects

This was a cross-sectional, observational study using data from the severe asthma adult cohort of U-BIOPRED. Details on the design of U-BIOPRED have been published previously [24]. Severe asthma was defined according to the Innovative Medicine Initiative consensus statement [25]. Patients were included who had a confirmed diagnosis of asthma that was uncontrolled despite high-intensity asthma treatment (≥1000 µg fluticasone equivalent plus a second controller) or could only be controlled with treatment by systemic oral corticosteroids or omalizumab, including current smokers, ex-smokers and never-smokers. The U-BIOPRED study was registered at ClinicalTrials.gov with identifier NCT01976767 and was approved by the Medical Ethics Boards of all participating centres. All patients provided written informed consent.

Persistent airflow limitation

PAL was defined as post-bronchodilator forced expiratory volume in 1 s (FEV1)/forced vital capacity (FVC) ratio below the lower limit of normal (LLN). The LLN was calculated according to formulae provided by Quanjer et al. [26] (for male subjects: −0.18*AGE+75.41; for female subjects: −0.19*AGE+78.4), consistent with the LLN defined by the Global Lung Function Initiative (GLI) [27]. In order to examine consistency with other criteria for PAL we additionally performed an analysis based on PAL defined by post-bronchodilator FEV1 <LLN.

Measurements of mRNA samples

Gene expression was assessed on the total RNA in unselected samples of severe asthma patients in the following four airway compartments: 1) nasal brushings (n=37), 2) induced sputum (n=79), 3) endobronchial brushings (n=62) and 4) endobronchial biopsies (n=50). There was no complete overlap between the samples, mainly because of limited sputum induction success rate (figure 1). For details on sample collection, see the supplementary material.

FIGURE 1
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 1

Venn diagram indicating the overlap between the samples: nasal brushings (NBr) n=37, sputum n=79, endobronchial brushings (EBr) n=62 and endobronchial biopsies (EBx) n=50. Overlap of NBr/EBx n=0; overlap of sputum/EBx/NBr n=3.

Microarray analysis was performed with the HT HG-U133+ PM microarray platform (Affymetrix, Santa Clara, CA, USA). For details, see the supplementary material.

Statistical analysis

For details on the statistical analysis on clinical variables, see the supplementary material.

The GSVA technique allows sensitive identification of differences in expression of predefined sets of genes between heterogeneous groups and can be used to explore underlying pathways [22]. Sets of genes used were based on all available gene expression publications and data on airways disease, including human and murine models, from both in vivo and in vitro studies, to ensure sensitive pathway detection. These contained studies on the gene expression associated with 1) airways disease treatment, 2) immunological pathways and 3) induced lung injury or inflammation. The latter included gene signatures that were identified at several time-points after the admission of either polyinosinic–polycytidylic acid (poly(I:C)) as a model for exacerbations [28] or bleomycin, which is used as a model for mimicking the course of fibrotic processes in the lung [29]. For this discovery study, 105 predefined gene sets were entered into the statistical models (supplementary table E8).

Enrichment scores and false discovery

GSVA was used to calculate the enrichment score (ES) for each patient and for each of the gene signatures based on the gene expression of genes in the sets. ES values range from −1 to 1 [22]. Subsequently, mean ES values were calculated for patients with and without PAL, and Generalised Linear Models, including correction for smoking status, corticosteroid usage and duration of asthma, were applied to statistically compare ES values between the groups. To minimise the false discovery rate, only gene signatures that had p<0.05 and a difference of ES (dES) between the groups of ≥0.2 were considered significantly different, following the Microarray Consortium for Quality Control (MACQC) recommendations regarding the need for applying group-difference thresholds in order to stringently limit false discovery [30].

Results

Of 421 severe asthma patients in the U-BIOPRED cohort, 224 (58%) had PAL. These patients were significantly older and had longer asthma duration than severe asthma patients without PAL (table 1). In addition, they had significantly more eosinophilia in blood and sputum. Patients with and without PAL did not significantly differ in smoking behaviour, number of pack-years and exacerbation rate. Some distinct patient characteristics were found when patients were grouped according to site of mRNA collection. Significantly fewer females and higher blood neutrophil counts were observed in patients with PAL in the subsets of nasal and endobronchial brushings. In addition, the number of pack-years was higher in patients with PAL in the subsets of endobronchial brushings and biopsies (table 2 and see supplementary tables E1–E4 for p-values).

View this table:
  • View inline
  • View popup
TABLE 1

Characteristics of patients with persistent airflow limitation (PAL)

View this table:
  • View inline
  • View popup
TABLE 2

Characteristics of patients with persistent airflow limitation (PAL) according to site of mRNA collection

Identified gene signatures

In total, GSVA identified 14 unique signatures that were significantly differentially enriched in patients with PAL compared with those without PAL: n=1 in the nasal brushings, n=9 in sputum samples, n=1 in endobronchial brushings and n=4 in endobronchial biopsies (table 3 and figures 2–5).

View this table:
  • View inline
  • View popup
TABLE 3

Differentially enriched gene signatures in severe asthma patients with persistent airflow limitation (PAL) compared with those without PAL

FIGURE 2
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 2

Significant gene signatures identified in nasal brushing samples. GSVA: Gene Set Variation Analysis; ES: enrichment score; PAL: persistent airflow limitation. Immunological gene signature: interleukin-13 (type 2 high):  up.

Treatment gene signatures

Two gene signatures that were expressed to a lesser extent in asthmatic and COPD patients after treatment with fluticasone were significantly differentially enriched. The gene signature lowly expressed in fluticasone-treated asthma patients (supplementary table E7, signature 12) was significantly more enriched in the bronchial brushings of patients with PAL (dES=0.22; p=0.0012). Similarly, the lowly expressed gene signature after fluticasone treatment in COPD patients, containing a different set of genes (supplementary table E7, signature 34), was also enriched to a greater extent in patients with PAL (dES=0.22; p=0.045) (table 3, and figures 4 and 5).⇓

FIGURE 3
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 3

Significant gene signatures identified in sputum samples. GSVA: Gene Set Variation Analysis; ES: enrichment score; PAL: persistent airflow limitation. Immunological gene signatures: a) eosinophils: up and b) interferon-α: up. Induced lung injury gene signatures: c) polyinosinic–polycytidylic acid (poly(I:C)), 2 h: up, d) poly(I:C), 6 h: up, e) poly(I:C), 48 h: up, f) poly(I:C), 72 h: up, g) poly(I:C), 96 h: up, h) bleomycin, day 1: up and i) bleomycin, day 28: down.

FIGURE 4
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 4

Significant gene signatures identified in endobronchial brushing samples. GSVA: Gene Set Variation Analysis; ES: enrichment score; PAL: persistent airflow limitation. Treatment gene signature: fluticasone treatment in asthma: down.

FIGURE 5
  • Download figure
  • Open in new tab
  • Download powerpoint
FIGURE 5

Significant gene signatures identified in endobronchial biopsy samples. GSVA: Gene Set Variation Analysis; ES: enrichment score; PAL: persistent airflow limitation. Treatment gene signature: a) fluticasone treatment in chronic obstructive pulmonary disease: down. Immunological gene signature: b) CD4+ T-cells in rheumatoid arthritis: up. Induced lung injury gene signatures: c) bleomycin, day 1: up and d) bleomycin, day 28: down.

Immunological gene signatures

Four gene signatures that are associated with immunological pathways were identified to be significantly differentially enriched between the groups. 1) A gene signature that is upregulated in the presence of eosinophils (supplementary table E7, signature 55) was more enriched in the sputum of patients with PAL compared with those without PAL (dES=0.21; p=0.04) (figure 3). 2) A gene signature that is identified to be associated with the presence of interleukin (IL)-13 (supplementary table E7, signature 101) was more enriched in the nasal brushings of PAL patients (dES=0.27; p=0.045) (figure 2). 3) The ES for the interferon (IFN)-α gene signature (supplementary table E7, signature 36) was less enriched in the sputum of patients with PAL (dES=0.32; p<0.001) (figure 3). 4) The ES for CD4+ T-cells of rheumatoid arthritis patients (supplementary table E7, signature 27) was more enriched in patients with PAL (dES=0.21; p=0.029) (figure 3).

Induced lung injury gene signatures

In total, eight gene signatures that are associated with induced airway inflammation or lung injury were identified to be significantly differentially enriched. Five of them were discovered in a study in which gene expression was assessed at different time-points after the admission of poly(I:C) (supplementary table E7, signatures 57, 58, 62, 64 and 65), mimicking viral replication. Patients with PAL had significantly lower enrichment of these signatures in sputum, ranging from dES= −0.25 to −0.20 (figure 3). In addition, three gene signatures associated with bleomycin-induced lung injury (supplementary table E7, signatures 67, 68 and 77) were significantly differentially enriched, one associated with higher expression of genes during early development of injury (supplementary table E7, signature 67) demonstrating a significantly lower ES in sputum (dES= −0.24; p=0.012) and a higher ES in endobronchial biopsy (dES=0.22; p=0.045) in PAL patients. In addition, the gene signatures identified to be expressed higher at later stages of bleomycin-induced injury (supplementary table E7, signatures 67 and 77) had a lower ES in PAL patients (dES= −0.20; p=0.006 and dES= −0.41; p=0.026) (table 3 and figure 5).

Consistency when using different definitions of PAL

The data when using post-bronchodilator FEV1 for defining PAL are given in the supplementary material. The data showed that, in general, signatures identified between FEV1/FVC <LLN and FEV1 <LLN were similar, including those associated with treatment with fluticasone, eosinophilic inflammation, and involvement of T-helper type 2 (Th2) cells and IFN-α. In addition, induced lung injury and inflammation gene signatures were also identified in both analyses (supplementary table E6). However, specific involvement of interleukin (IL)-13 and CD4+ T-cells of rheumatoid arthritis was not consistent, and no significant gene signatures were identified in the biopsies when applying FEV1 <LLN.

Discussion

This study shows that severe asthma patients with PAL have differentially enriched gene sets compared with patients without PAL, suggesting that distinguishable underlying mechanisms are involved. The identified gene signatures may indicate that these patients respond or have responded differently to treatment with inhaled corticosteroids and that particular immunological pathways are involved. These include the already recognised eosinophilic airway inflammation, but also specific CD4+ T-cells and involvement of IL-13, possibly via Th2 cells. In addition, a differentially expressed gene network associated with remodelling was identified in the endobronchial airway wall, while not being identified in other compartments, underlining local involvement of remodelling. These results show that PAL is biologically distinguishable in severe asthma and point towards possible targets for the prevention and treatment of PAL.

This study is the first to assess gene signatures underlying PAL for the discovery of its mechanisms in severe asthma. Clinical characterisation of the present patients with PAL showed no association with exacerbation rate. Multiple studies in severe asthma have shown mixed results. Some found associations between more frequent exacerbations and PAL, especially in certain phenotypes of severe asthma [10, 11], whereas others observed no association [13] or even an inverse association indicative of more frequent exacerbations in patients without PAL [9]. This suggests that certain subgroups of severe asthma patients may suffer from both clinical traits. However, in a broad severe asthma population, such as the present study, PAL and exacerbations may need to be seen as separate clinical traits.

GSVA not only confirmed expected mechanisms, but also identified novel pathways in PAL. First, the results indicate that the airways of severe asthma patients with PAL respond differently to the administration of inhaled corticosteroids. Even though a relationship between severe asthma and resistance to steroids has been recognised [31], a specific association with PAL has not previously been reported. Second, not surprisingly, a gene set associated with eosinophilic airway inflammation was identified in the sputum of patients with PAL. Eosinophilic airway inflammation has been reported in patients with PAL by multiple studies in severe asthma and was also present in this study population [9, 14, 15]. Third, IL-13 was identified to be involved in PAL in severe asthma, which was in concordance with a study by Kaminska et al. [32] who identified IL-13 in bronchial biopsies of severe asthma patients with PAL.

Interestingly, in sputum samples, we identified a gene set associated with IFN-α that was less enriched in patients with PAL. Impaired or decreased IFN-α has been associated with severe asthma compared with controls [33, 34] and, in a murine asthma model, where IFN-α inhibits eosinophilic airway inflammation and hyperresponsiveness [35]. However, a direct association with PAL in (severe) asthma has not been described previously. In addition, patients with PAL appear to have lower enrichment of gene sets expressed after exposure to poly(I:C), mimicking viral replication upon infection. The recognised link between poly(I:C) and IFNs strengthens these results [36]. Impaired lung function has been found after exposure to poly(I:C) in mice [37]; however, it is unknown to what extent the ability to respond to viral infection influences lung function. Moreover, the present data indicate a role for CD4+ T-cells, which have been associated previously with reduced lung function in patients with asthma [18]. Finally, remodelling has been recognised to be associated with PAL in (severe) asthma [17, 18]; in concert with this, our data identified gene sets associated with bleomycin-induced lung injury and fibrotic processes in the endobronchial airway wall to be more enriched in severe asthma patients with PAL.

We believe the strengths of this study include the following. It was conducted in the severe asthma cohort of the pan-European U-BIOPRED study, which included a large number of well-selected and defined patients. The study included severe asthma patients with PAL and those without PAL with equally severe disease as an appropriate control group for analysing the distinguishing gene signatures of PAL among patients with severe asthma. The multicentre design allowed standardised collection of large numbers of samples from patients all over Europe. In addition, collection of these samples was standardised and performed according to standard operating procedures, while mRNA expression in these samples was analysed centrally to gain optimised data quality. Furthermore, using GSVA as a sensitive statistical technique enabled the identification of underlying pathways involved in PAL even in heterogeneous samples, such as in severe asthma. Finally, selection bias was minimised by correction for smoking status, asthma duration and usage of oral corticosteroids as such. However, correction for corticosteroid dose was not indicated as this was similar between the groups. Similarly, no correction for cell counts was applied as this could partly obscure the disease mechanism of PAL in severe asthma.

This study also has limitations. First, this was a cross-sectional study, whereas lung function may vary over time. Stringent selection criteria for PAL were applied to minimise possible misdiagnosis of patients, but we cannot exclude longitudinal variations in PAL. The definition of PAL was derived from the LLN on two measures of PAL based on the Quanjer et al. [26] values. Applying the GLI reference equations for the LLN to our dataset had very little influence on the distribution of patients with or without PAL, since in total only three patients included in the GSVA analyses would switch groups. Second, as shown in figure 1, there is incomplete overlap of patients from whom various samples had been collected, which was unintended but inevitable due to patient consent, feasibility and quality control (e.g. for bronchoscopies and induced sputum). This hampers comparison of the identified gene sets between the different compartments. The sample size for nasal brushings, for example, was modest. Nevertheless, concordance of identified gene signatures between the different compartments (e.g. on gene sets for fluticasone treatment and bleomycin-induced injury) and concordance in identified gene signatures by applying two definitions for PAL lend validity to the identified pathways in the different samples. Third, GSVA requires predefined gene sets to be put into the analysis. These gene signatures are acquired from the literature and are based on other studies reported thus far. Inevitably, those studies had their own flaws as well, which were not taken into account in this study. In addition, some of the gene signatures are derived from mouse model studies. These data from mouse models were purposely included, because only by broad, nonselective input of gene signatures into the analysis can novel mechanisms involved in PAL be discovered. However, it needs to be emphasised that murine models will certainly deviate from human disease and it remains to be established whether the murine signatures can be used to identify phenotypic differences in severe asthma. Furthermore, by applying predefined gene sets, we may have gene interactions that are not yet annotated and thus not identified in this study. Finally, this study did not include a separate validation set or internal validation. The sample sizes were not large enough to split into a training and validation set. However, we have applied the most stringent measures for the prevention of false discovery of gene signatures by only selecting gene signatures with a low p-value and a minimum dES, following the standardised MACQC recommendations for the importance of applying group-difference thresholds for gaining the best reproducible results [30]. In addition, the results were largely consistent when defining PAL by using FEV1 or FEV1/FVC, which contributes to the validity of the data.

Based on the present results, a few well-defined mechanisms are revealed in relation to PAL in severe asthma. First, involvement of IL-13 does not seem to be surprising. Not only because it has been directly associated with PAL previously [32], but also because it has been linked to the other identified pathways in this study. IL-13 has been associated with resistance to steroids in asthma [31, 38] and is recognised as a promotor for local eosinophil recruitment in the lungs [39]. In addition, IL-13 has been identified as a direct promotor for airway remodelling [40]. Thus, IL-13 may be an important hub for the several pathways involved in PAL in severe asthma.

Second, IFN-α deficiency, which may impair antiviral responses [41], has been associated with persistent asthma as such [34]. In addition, a mouse model study indicates that IFN-α inhibits eosinophilic inflammation [33]. The results of our study imply that only the subgroup of severe asthma patients suffering from PAL exhibits IFN-α deficiency or impairment, which may well be associated with lower expression of gene sets associated with poly(I:C). One could speculate that IFN-α deficiency could result in enhanced viral loads during an infection or exacerbation and in augmented inflammation, which in turn could result in airway wall remodelling and eventually PAL. Finally, different gene signatures appeared to be enriched in the four compartments. Gene expression related to airway remodelling was found in the more profound layers of the airway wall (biopsies), but not in the superficial layers and lumen. The opposite is true for gene signatures related to inflammatory involvement, such as eosinophilia and CD4+ T-cells. This seems to be highly relevant, perhaps not unexpectedly implying differential pathways at different locations in the respiratory tract.

PAL has been recognised as a treatable trait in chronic airways disease [8]. In particular, severe asthma patients with PAL have high disease burden and therefore are in need of more individualised management, so-called “precision medicine”. The present results suggest that patients with PAL may be less responsive to inhaled steroids. Taking into account that these patients have more symptoms and higher rates of exacerbations and hospitalisations [10], they seem to be the patients in most need for targeted medicines, such as newly developed biologicals. For instance, the involvement of IL-13 pathways in these patients may make them primary candidates for treatment with anti-IL-13.

In conclusion, severe asthma patients with PAL express different gene mechanisms compared with those without, being suggestive of specific disease mechanisms, including IL-13-associated disease pathways such as eosinophilic airway inflammation and remodelling with involvement of impaired IFN-α function. Therefore, these may be the patients to be selected for studies on emerging new treatments.

Supplementary material

Supplementary Material

Please note: supplementary material is not edited by the Editorial Office, and is uploaded as it has been supplied by the author.

Supplementary material ERJ-02298-2016_Supplement

Supplementary table E7 ERJ-02298-2016_Table_E7

Supplementary table E8 ERJ-02298-2016_Table_E8

Disclosures

Supplementary Material

I.M. Adcock ERJ-02298-2016_Adcock

C. Auffray ERJ-02298-2016_Auffray

F. Baribaud ERJ-02298-2016_Baribaud

S.A. Bates ERJ-02298-2016_Bates

E.H. Bel ERJ-02298-2016_Bel

P. Brinkman ERJ-02298-2016_Brinkman

K.F. Chung ERJ-02298-2016_Chung

B. De Meulder ERJ-02298-2016_De_Meulder

R. Djukanovic ERJ-02298-2016_Djukanovic

L.J. Fleming ERJ-02298-2016_Fleming

P-P. Hekking ERJ-02298-2016_Hekking

P. Howarth ERJ-02298-2016_Howarth

D. Lefaudeux ERJ-02298-2016_Lefaudeux

M.J. Loza ERJ-02298-2016_Loza

D.E. Shaw ERJ-02298-2016_Shaw

A.R. Sousa ERJ-02298-2016_Sousa

P.J. Sterk ERJ-02298-2016_Sterk

Footnotes

  • This article has supplementary material available from erj.ersjournals.com

  • Support statement: U-BIOPRED is funded by IMI. Funding information for this article has been deposited with the Crossref Funder Registry.

  • Conflict of interest: Disclosures can be found alongside this article at erj.ersjournals.com

  • Received November 23, 2016.
  • Accepted June 26, 2017.
  • Copyright ©ERS 2017

References

  1. ↵
    1. Chung KF,
    2. Wenzel SE,
    3. Brozek JL
    , et al. International ERS/ATS guidelines on definition, evaluation and treatment of severe asthma. Eur Respir J 2014; 43: 343–373.
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Hekking PP,
    2. Wener RR,
    3. Amelink M
    , et al. The prevalence of severe refractory asthma. J Allergy Clin Immunol 2015; 135: 896–902.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Hekking PP,
    2. Bel EH
    . Developing and emerging clinical asthma phenotypes. J Allergy Clin Immunol Pract 2014; 2: 671–680.
    OpenUrlPubMed
  4. ↵
    1. Moore WC,
    2. Hastie AT,
    3. Li X
    , et al. Sputum neutrophil counts are associated with more severe asthma phenotypes using cluster analysis. J Allergy Clin Immunol 2014; 133: 1557–1563.
    OpenUrlCrossRefPubMedWeb of Science
  5. ↵
    1. Bel EH,
    2. Wenzel SE,
    3. Thompson PJ
    , et al. Oral glucocorticoid-sparing effect of mepolizumab in eosinophilic asthma. N Engl J Med 2014; 371: 1189–1197.
    OpenUrlCrossRefPubMedWeb of Science
    1. Busse W,
    2. Corren J,
    3. Lanier BQ
    , et al. Omalizumab, anti-IgE recombinant humanized monoclonal antibody, for the treatment of severe allergic asthma. J Allergy Clin Immunol 2001; 108: 184–190.
    OpenUrlCrossRefPubMedWeb of Science
  6. ↵
    1. Chung KF
    . New treatments for severe treatment-resistant asthma: targeting the right patient. Lancet Respir Med 2013; 1: 639–652.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Agusti A,
    2. Bel E,
    3. Thomas M
    , et al. Treatable traits: toward precision medicine of chronic airway diseases. Eur Respir J 2016; 47: 410–419.
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. ten Brinke A,
    2. Zwinderman AH,
    3. Sterk PJ
    , et al. Factors associated with persistent airflow limitation in severe asthma. Am J Respir Crit Care Med 2001; 164: 744–748.
    OpenUrlCrossRefPubMedWeb of Science
  9. ↵
    1. Ciebiada M,
    2. Domagala M,
    3. Gorska-Ciebiada M
    , et al. Risk factors associated with irreversible airway obstruction in nonsmoking adult patients with severe asthma. Allergy Asthma Proc 2014; 35: 72–79.
    OpenUrl
  10. ↵
    1. Moore WC,
    2. Meyers DA,
    3. Wenzel SE
    , et al. Identification of asthma phenotypes using cluster analysis in the Severe Asthma Research Program. Am J Respir Crit Care Med 2010; 181: 315–323.
    OpenUrlCrossRefPubMedWeb of Science
    1. Amelink M,
    2. de Nijs SB,
    3. de Groot JC
    , et al. Three phenotypes of adult-onset asthma. Allergy 2013; 68: 674–680.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Schatz M,
    2. Hsu JW,
    3. Zeiger RS
    , et al. Phenotypes determined by cluster analysis in severe or difficult-to-treat asthma. J Allergy Clin Immunol 2014; 133: 1549–1556.
    OpenUrlCrossRef
  12. ↵
    1. Konno S,
    2. Taniguchi N,
    3. Makita H
    , et al. Distinct phenotypes of cigarette smokers identified by cluster analysis of patients with severe asthma. Ann Am Thorac Soc 2015; 12: 1771–1780.
    OpenUrl
  13. ↵
    1. Bumbacea D,
    2. Campbell D,
    3. Nguyen L
    , et al. Parameters associated with persistent airflow obstruction in chronic severe asthma. Eur Respir J 2004; 24: 122–128.
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Lee JH,
    2. Haselkorn T,
    3. Borish L
    , et al. Risk factors associated with persistent airflow limitation in severe or difficult-to-treat asthma: insights from the TENOR study. Chest 2007; 132: 1882–1889.
    OpenUrlCrossRefPubMedWeb of Science
  15. ↵
    1. Yamauchi K,
    2. Inoue H
    . Airway remodeling in asthma and irreversible airflow limitation – ECM deposition in airway and possible therapy for remodeling. Allergol Int 2007; 56: 321–329.
    OpenUrlCrossRefPubMed
  16. ↵
    1. den Otter I,
    2. Willems LN,
    3. van Schadewijk A
    , et al. Lung function decline in asthma patients with elevated bronchial CD8, CD4 and CD3 cells. Eur Respir J 2016; 48: 393–402.
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Kwak HJ,
    2. Park DW,
    3. Seo JY
    , et al. The Wnt/beta-catenin signaling pathway regulates the development of airway remodeling in patients with asthma. Exp Mol Med 2015; 47: e198.
    OpenUrlPubMed
  18. ↵
    1. Perry MM,
    2. Durham AL,
    3. Austin PJ
    , et al. BET bromodomains regulate transforming growth factor-beta-induced proliferation and cytokine release in asthmatic airway smooth muscle. J Biol Chem 2015; 290: 9111–9121.
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Pepe C,
    2. Foley S,
    3. Shannon J
    , et al. Differences in airway remodeling between subjects with severe and moderate asthma. J Allergy Clin Immunol 2005; 116: 544–549.
    OpenUrlCrossRefPubMedWeb of Science
  20. ↵
    1. Hanzelmann S,
    2. Castelo R,
    3. Guinney J
    . GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics 2013; 14: 7.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Lefaudeux D,
    2. De Meulder B,
    3. Loza MJ
    , et al. U-BIOPRED clinical adult asthma clusters linked to a subset of sputum -omics. J Allergy Clin Immunol 2017; 139: 1797–1807
    OpenUrl
  22. ↵
    1. Shaw DE,
    2. Sousa AR,
    3. Fowler SJ
    , et al. Clinical and inflammatory characteristics of the European U-BIOPRED adult severe asthma cohort. Eur Respir J 2015; 46: 1308–1321.
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Bel EH,
    2. Sousa A,
    3. Fleming L
    , et al. Diagnosis and definition of severe refractory asthma: an international consensus statement from the Innovative Medicine Initiative (IMI). Thorax 2011; 66: 910–917.
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Quanjer PH,
    2. Tammeling GJ,
    3. Cotes JE
    , et al. Lung volumes and forced ventilatory flows. Eur Respir J 1993; 6: Suppl. 16, 5–40.
    OpenUrlFREE Full Text
  25. ↵
    1. Quanjer PH,
    2. Stanojevic S,
    3. Cole TJ
    , et al. Multi-ethnic reference values for spirometry for the 3–95-yr age range: the global lung function 2012 equations. Eur Respir J 2012; 40: 1324–1343.
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Harris P,
    2. Sridhar S,
    3. Peng R
    , et al. Double-stranded RNA induces molecular and inflammatory signatures that are directly relevant to COPD. Mucosal Immunol 2013; 6: 474–484.
    OpenUrl
  27. ↵
    1. Peng R,
    2. Sridhar S,
    3. Tyagi G
    , et al. Bleomycin induces molecular changes directly relevant to idiopathic pulmonary fibrosis: a model for “active” disease. PLoS One 2013; 8: e59348.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Shi L,
    2. Reid LH,
    3. Jones WD
    , et al. The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements. Nat Biotechnol 2006; 24: 1151–1161.
    OpenUrlCrossRefPubMedWeb of Science
  29. ↵
    1. Barnes PJ
    . Corticosteroid resistance in patients with asthma and chronic obstructive pulmonary disease. J Allergy Clin Immunol 2013; 131: 636–645.
    OpenUrlCrossRefPubMedWeb of Science
  30. ↵
    1. Kaminska M,
    2. Foley S,
    3. Maghni K
    , et al. Airway remodeling in subjects with severe asthma with or without chronic persistent airflow obstruction. J Allergy Clin Immunol 2009; 124: 45–51.
    OpenUrlCrossRefPubMedWeb of Science
  31. ↵
    1. Wright AK,
    2. Mistry V,
    3. Richardson M
    , et al. Toll-like receptor 9 dependent interferon-alpha release is impaired in severe asthma but is not associated with exacerbation frequency. Immunobiology 2015; 220: 859–864.
    OpenUrl
  32. ↵
    1. Rupani H,
    2. Martinez-Nunez RT,
    3. Dennison P
    , et al. Toll-like receptor 7 is reduced in severe asthma and linked to an altered microRNA profile. Am J Respir Crit Care Med 2016; 194: 26–37.
    OpenUrlCrossRef
  33. ↵
    1. Kikkawa Y,
    2. Sugiyama K,
    3. Obara K
    , et al. Interferon-alpha inhibits airway eosinophilia and hyperresponsiveness in an animal asthma model [corrected]. Asia Pac Allergy 2012; 2: 256–263.
    OpenUrlPubMed
  34. ↵
    1. Sionov E,
    2. Mayer-Barber KD,
    3. Chang YC
    , et al. Type I IFN induction via poly-ICLC protects mice against cryptococcosis. PLoS Pathog 2015; 11: e1005040.
    OpenUrlCrossRefPubMed
  35. ↵
    1. Stowell NC,
    2. Seideman J,
    3. Raymond HA
    , et al. Long-term activation of TLR3 by poly(I:C) induces inflammation and impairs lung function in mice. Respir Res 2009; 10: 43.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Saha SK,
    2. Berry MA,
    3. Parker D
    , et al. Increased sputum and bronchial biopsy IL-13 expression in severe asthma. J Allergy Clin Immunol 2008; 121: 685–691.
    OpenUrlCrossRefPubMedWeb of Science
  37. ↵
    1. Pope SM,
    2. Brandt EB,
    3. Mishra A
    , et al. IL-13 induces eosinophil recruitment into the lung by an IL-5- and eotaxin-dependent mechanism. J Allergy Clin Immunol 2001; 108: 594–601.
    OpenUrlCrossRefPubMedWeb of Science
  38. ↵
    1. Zhu Z,
    2. Homer RJ,
    3. Wang Z
    , et al. Pulmonary expression of interleukin-13 causes inflammation, mucus hypersecretion, subepithelial fibrosis, physiologic abnormalities, and eotaxin production. J Clin Invest 1999; 103: 779–788.
    OpenUrlCrossRefPubMedWeb of Science
  39. ↵
    1. Sykes A,
    2. Edwards MR,
    3. Macintyre J
    , et al. Rhinovirus 16-induced IFN-alpha and IFN-beta are deficient in bronchoalveolar lavage cells in asthmatic patients. J Allergy Clin Immunol 2012; 129: 1506–1514.
    OpenUrlCrossRefPubMedWeb of Science
View Abstract
PreviousNext
Back to top
View this article with LENS
Vol 50 Issue 3 Table of Contents
European Respiratory Journal: 50 (3)
  • Table of Contents
  • Index by author
Email

Thank you for your interest in spreading the word on European Respiratory Society .

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Transcriptomic gene signatures associated with persistent airflow limitation in patients with severe asthma
(Your Name) has sent you a message from European Respiratory Society
(Your Name) thought you would like to see the European Respiratory Society web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
Citation Tools
Transcriptomic gene signatures associated with persistent airflow limitation in patients with severe asthma
Pieter-Paul Hekking, Matthew J. Loza, Stelios Pavlidis, Bertrand De Meulder, Diane Lefaudeux, Frederic Baribaud, Charles Auffray, Ariane H. Wagener, Paul Brinkman, René Lutter, Aruna T. Bansal, Ana R. Sousa, Stewart A. Bates, Ioannis Pandis, Louise J. Fleming, Dominick E. Shaw, Stephen J. Fowler, Yike Guo, Andrea Meiser, Kai Sun, Julie Corfield, Peter Howarth, Elisabeth H. Bel, Ian M. Adcock, Kian Fan Chung, Ratko Djukanovic, Peter J. Sterk
European Respiratory Journal Sep 2017, 50 (3) 1602298; DOI: 10.1183/13993003.02298-2016

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero

Share
Transcriptomic gene signatures associated with persistent airflow limitation in patients with severe asthma
Pieter-Paul Hekking, Matthew J. Loza, Stelios Pavlidis, Bertrand De Meulder, Diane Lefaudeux, Frederic Baribaud, Charles Auffray, Ariane H. Wagener, Paul Brinkman, René Lutter, Aruna T. Bansal, Ana R. Sousa, Stewart A. Bates, Ioannis Pandis, Louise J. Fleming, Dominick E. Shaw, Stephen J. Fowler, Yike Guo, Andrea Meiser, Kai Sun, Julie Corfield, Peter Howarth, Elisabeth H. Bel, Ian M. Adcock, Kian Fan Chung, Ratko Djukanovic, Peter J. Sterk
European Respiratory Journal Sep 2017, 50 (3) 1602298; DOI: 10.1183/13993003.02298-2016
del.icio.us logo Digg logo Reddit logo Technorati logo Twitter logo CiteULike logo Connotea logo Facebook logo Google logo Mendeley logo
Full Text (PDF)

Jump To

  • Article
    • Abstract
    • Abstract
    • Introduction
    • Methods
    • Results
    • Discussion
    • Supplementary material
    • Disclosures
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF

Subjects

  • Asthma and allergy
  • Tweet Widget
  • Facebook Like
  • Google Plus One

More in this TOC Section

Original articles

  • Ambulatory management of secondary spontaneous pneumothorax
  • Systematic assessment of respiratory health in illness susceptible athletes
  • Identifying early PAH biomarkers in systemic sclerosis
Show more Original articles

Asthma

  • Zinc finger protein 33B and atopy-related markers in asthma
  • Green and blue spaces and lung function in the Generation XXI cohort
  • Adrenal function recovery after durable OCS sparing with benralizumab
Show more Asthma

Related Articles

Navigate

  • Home
  • Current issue
  • Archive

About the ERJ

  • Journal information
  • Editorial board
  • Reviewers
  • Press
  • Permissions and reprints
  • Advertising

The European Respiratory Society

  • Society home
  • myERS
  • Privacy policy
  • Accessibility

ERS publications

  • European Respiratory Journal
  • ERJ Open Research
  • European Respiratory Review
  • Breathe
  • ERS books online
  • ERS Bookshop

Help

  • Feedback

For authors

  • Instructions for authors
  • Publication ethics and malpractice
  • Submit a manuscript

For readers

  • Alerts
  • Subjects
  • Podcasts
  • RSS

Subscriptions

  • Accessing the ERS publications

Contact us

European Respiratory Society
442 Glossop Road
Sheffield S10 2PX
United Kingdom
Tel: +44 114 2672860
Email: journals@ersnet.org

ISSN

Print ISSN:  0903-1936
Online ISSN: 1399-3003

Copyright © 2023 by the European Respiratory Society