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.
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).
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).
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).⇓
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