Skip to main content


Sputum and blood transcriptomics characterisation of the inhaled PDE4 inhibitor CHF6001 on top of triple therapy in patients with chronic bronchitis



Although phosphodiesterase-4 (PDE4) inhibitors have been shown to reduce COPD exacerbation rate, their biological mechanism of action is not completely elucidated at the molecular level. We aimed to characterise the whole genome gene expression profile of the inhaled PDE4-inhibitor CHF6001 on top of triple therapy in sputum cells and whole blood of patients with COPD and chronic bronchitis.


Whole genome gene expression analysis was carried out by microarray in 54 patients before and after 32 days treatment with CHF6001 800 and 1600 μg and placebo twice daily (BID) in a randomised crossover study.


CHF6001 had a strong effect in sputum, with 1471 and 2598 significantly differentially-expressed probe-sets relative to placebo (p-adjusted for False Discovery Rate < 0.05) with 800 and 1600 μg BID, respectively. Functional enrichment analysis showed significant modulation of key inflammatory pathways involved in cytokine activity, pathogen-associated-pattern-recognition activity, oxidative stress and vitamin D with associated inhibition of downstream inflammatory effectors. A large number of pro-inflammatory genes coding for cytokines and matrix-metalloproteinases were significantly differentially expressed for both doses; the majority (> 87%) were downregulated, including macrophage inflammatory protein-1-alpha and 1-beta, interleukin-27-beta, interleukin-12-beta, interleukin-32, tumour necrosis factor-alpha-induced-protein-8, ligand-superfamily-member-15, and matrix-metalloproteinases-7,12 and 14. The effect in blood was not significant.


Inhaled PDE4 inhibition by CHF6001 on top of triple therapy in patients with COPD and chronic bronchitis significantly modulated key inflammatory targets and pathways in the lung but not in blood. Mechanistically these findings support a targeted effect in the lung while minimising unwanted systemic class-effects.

Trial registration, EudraCT, 2015–005550-35. Registered 15 July 2016.


The orally administered phosphodiesterase 4 (PDE4) inhibitor roflumilast has been shown to reduce the rate of exacerbations in patients with COPD who have a chronic bronchitis phenotype [1]. Clinical trials have demonstrated the greatest benefit is in patients who are at high risk of frequent or severe exacerbations despite inhaled corticosteroid (ICS) and long-acting β2 agonist therapy (LABA), with or without a long-acting muscarinic antagonist (LAMA) [2, 3]. However, systemic exposure after oral administration of roflumilast can cause side effects such as nausea, weight loss and gastrointestinal disturbance [4]. Furthermore, although in clinical practice a PDE4 inhibitor is often administered in addition to inhaled triple ICS/LABA/LAMA therapy [5], there is a lack of information regarding the anti-inflammatory effects of a PDE4 inhibitor administered in addition to maintenance triple therapy.

CHF6001 is a novel PDE4 inhibitor currently in clinical development that has been specifically formulated as an inhaled extrafine formulation (i.e., mass median aerodynamic diameter ≤ 2 μm) [6], and designed to have high protein binding and rapid elimination from the systemic circulation [7] thus maximising exposure in the lung and avoiding systemic adverse effects [8]. In a 24-week dose finding study, when administered in addition to LABA therapy, CHF6001 numerically reduced the exacerbation rate in patients with chronic bronchitis, with all doses being well tolerated [9].

We previously reported that CHF6001 significantly decreased the levels of various inflammatory mediators in the sputum of patients with COPD and a chronic bronchitis phenotype who are receiving inhaled maintenance triple therapy [10]. The current analyses aimed to validate and further characterise the biological effect of CHF6001 using high throughput gene expression analysis of the whole protein coding genome in sputum cells and blood in the patients recruited into the inflammatory mediator study [10].


Study objective and design

Samples were collected from a multicentre, three-way, placebo-controlled, double-blind crossover study, the results for which have been previously reported [10]. This study aimed to evaluate the effect of CHF6001 on inflammatory biomarkers in induced sputum and blood. The gene expression analyses described in this manuscript were a prespecified exploratory objective of the study.

After randomisation, eligible patients commenced three, 32-day treatment periods during which they received CHF6001 800 or 1600 μg twice daily (BID, total daily doses of 1600 or 3200 μg) or matching placebo, all via multi-dose dry-powder inhaler (NEXThaler). Treatment periods were separated by a 28–42 day washout. Induced sputum was collected pre-dose on Day 1, and 2 h post-dose on Days 20, 26 and 32 of each period, with blood samples collected pre-dose on Day 1 and 2 h post-dose on Day 32 (see the supplement for further detail, including methods for the sputum and blood collection, and processing for the ribonucleic acid [RNA] assessments, extraction and amplification, sample profiling, microarray data quality control and pre-processing). The trial was approved by the independent ethics committees at each institution, and was performed in accordance with the principles of the Declaration of Helsinki, and the International Conference on Harmonisation notes for guidance on Good Clinical Practice (ICH/CPMP/135/95). The trial is registered on (NCT03004417) and EudraCT (2015–005550-35#). There were no protocol amendments. The data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus and are accessible through GEO Series accession number GSE133513 (


Eligible patients were male or female, ≥40 years of age, current or ex-smokers with a smoking history ≥10 pack-years, a diagnosis of COPD, post-bronchodilator forced expiratory volume in 1 s (FEV1) ≥ 30 and < 70% predicted, ratio of FEV1 to forced vital capacity < 0.70, COPD Assessment Test score ≥ 10, a history of chronic bronchitis (defined as chronic cough and sputum production for more than 3 months per year for at least two consecutive years), and treated with inhaled triple ICS/LABA/LAMA therapy for at least 2 months prior to enrolment. All patients provided written informed consent prior to any study-related procedure.

The key exclusion criteria were a moderate or severe COPD exacerbation within 6 weeks prior to entry or between screening and randomisation, and the use of PDE4 inhibitors within 2 months prior to entry. Full inclusion and exclusion criteria can be found in the supplement.

Statistical analysis

For the sputum analyses, patients had to have a minimum of one sample on Day 20, 26 or 32 that passed quality control, with the latest sample used for the analyses. If no matching pre-dose sample was available, the pre-dose sample from the closest available period was used.

Any probe sets with statistically significantly different pre-dose expression between treatments were identified and excluded (see ANOVA analysis, supplementary methods). An ANCOVA model was fitted to identify probe sets responding differentially to treatment, with change from pre-dose to post-dose expression as dependent variable, and subject, period, treatment and pre-dose expression as independent variables. For all analyses, statistical significance was considered as a p-value adjusted for False Discovery Rate using the Benjamini-Hochberg method (pFDR) below 0.05 [11].

Functional enrichment analysis

Differentially expressed probe set lists were annotated with gene identifiers using the latest annotation provided by Affymetrix for the Plus 2.0 array. Differential gene lists smaller than 10 genes were not used for functional enrichment; lists larger than 1000 genes were filtered based on fold change (FC) using a cut-off of 1.3. For multiple probe sets associated with one gene, the probe set with the lowest pFDR value was used for analysis. Using Gene Ontology (GO) annotations and commercial software Ingenuity® Pathway Analysis (IPA®, QIAGEN Inc.,, functional enrichment analysis was conducted to identify and rank biological entities associated with gene sets of interest. Entities were ranked according to a statistically derived enrichment score and were adjusted for multiple testing (pFDR < 0.05). Simulation of perturbations for significant IPA® canonical pathways was carried out using Molecule Activity Predictor (MAP) to simulate the downstream consequences of up or downregulation of molecules mediated by the treatment.

Gene interaction network analysis

The genes of interest were input into the STRING database (Search Tool for Retrieval of Interacting Genes/proteins) version 11.0 ( [12]. When multiple probe sets were associated with the same gene of interest, the probe set with the lowest pFDR value was considered. For probe sets associated with more than one gene, all genes were considered for the network analysis.


Patients and samples

Fifty-nine patients gave consent for gene expression investigations and were randomised, 54 of whom completed the study and were eligible for analysis, having pre- and post-dose samples available for at least one active dose and placebo in blood and/or sputum (Table 1). Sputum weight and the number of sputum cells/g were similar for the two CHF6001 doses (Fig. S1). A small number of post-dose samples did not have matching pre-dose samples. An analysis of principal components on pre-dose samples showed no association between period (1, 2 or 3) and gene expression data (principal component [PC] 1, p = 0.71 in sputum, p = 0.48 in blood), or between treatment and gene expression data (PC1, p = 0.94 in sputum, p = 0.38 in blood) (Fig. S2). Missing pre-dose samples were therefore replaced with the pre-dose samples from the closest available period (see supplement). Following quality and RNA quantity assessment, 50 patients had pre- and post-dose blood samples available for analysis with both CHF6001 doses and placebo, 37 of whom had pre- and post-dose sputum samples available for analysis with CHF6001 800 μg BID and placebo, with 41 having samples available with CHF6001 1600 μg BID and placebo (see supplement).

Table 1 Baseline demographics and disease characteristics

Differential expression analysis in response to treatment

Filtering to remove any uninformative transcripts (see supplement) resulted in 45,163 and 44,355 reliably detected probe sets in blood and sputum, respectively, being available for analysis, corresponding to nearly the whole protein coding genome of 19,000 genes [13]. Notably, no probe sets had a pre-dose expression for either CHF6001 dose that significantly varied from the placebo pre-dose value.

The impact of treatment on gene expression in blood was minimal with either dose (Fig. 1a). In contrast, in sputum a large number of probe sets were differentially expressed between CHF6001 and placebo, the 1600 μg BID dose having the greatest effect (Fig. 1b). With CHF6001 800 μg BID, 1471 probe sets were significantly differentially expressed vs. placebo (pFDR< 0.05), 390 (27%) of which were downregulated and 1081 (73%) upregulated. CHF6001 1600 μg BID led to 3598 significantly differentially expressed probe sets vs. placebo, with 1226 (34%) downregulated and 2372 (66%) upregulated. Further, of the 1471 probe sets significantly differentially expressed with CHF6001 800 μg BID, 1229 (84%) were also significantly differentially expressed with CHF6001 1600 μg BID.

Fig. 1

Differentially expressed probe sets for CHF6001 800 μg and 1600 μg BID relative to placebo in (A) blood and (B) sputum. Volcano plot depicting all detected probe sets and coloured by fold change (FC) and adjusted p-values: orange, |FC| > 1.3; red, pFDR< 0.05; green, |FC| > 1.3 and pFDR< 0.05. Abbreviation: BID, twice daily

Among the probe sets differentially expressed with at least one CHF6001 dose, > 99% followed the same direction for both treatments, i.e., either both up-regulated or both down-regulated (Fig. 2), with 90.1% of the upregulated probe sets and 79.1% of the downregulated probe sets showing a dose-response relationship.

Fig. 2

Log2 fold change with CHF6001 1600 μg BID vs. log2 fold change with CHF6001 800 μg BID. Equality line in red, fitted loess curve in blue. Only differentially expressed probe sets (pFDR< 0.05 for at least one dose) are presented. Abbreviation: BID, twice daily

Functional enrichment analysis

In blood, functional enrichment analysis was not performed due to insufficient numbers of probe sets being available. In sputum, more than 1000 probe sets were significantly differentially expressed with at least one CHF6001 dose, so the analysis was restricted to probe sets with |FC| > 1.3. A total of 24 and 92 ingenuity canonical pathways were significantly modulated by CHF6001 800 and 1600 μg BID, respectively (Table S1), the most commonly associated terms being: immune response, cytokine signalling and inflammation, growth and development. Sixteen pathways associated with COPD pathophysiology were identified for CHF6001 800 μg BID, and 65 for CHF6001 1600 μg BID, covering cytokines, cyclic adenosine monophosphate (cAMP), T-helper (Th) 1 and 2 cells, dendritic cells, natural killer cells, oxidative stress, and vitamin D pathways (Table S1), with nine pathways in common between the two CHF6001 doses.

Simulation of perturbations for all significant IPA® canonical pathways was carried out for CHF6001 1600 μg BID using Molecule Activity Predictor to simulate the downstream consequences of treatment. This identified several key COPD pathways associated with consistent actual downregulation and/or predicted inhibition of downstream inflammatory mediators or mechanisms (Fig. S3) namely, C-C chemokine receptor type 5 (CCR5) signalling in macrophages, chemokine signalling, high mobility group-B1 (HMGB1) signalling, triggering receptor expressed on myeloid cells 1 (TREM1) signalling, interleukin (IL)22 signalling, ceramide signalling, protein kinase c theta (PKCθ) signalling in T lymphocytes, role of pattern recognition receptors in recognition of bacteria and viruses, toll-like receptors signalling, N-formylmethionyl-leucyl-phenylalanine (fMLP) signalling in neutrophils, production of nitric oxide and reactive oxygen species in macrophages and vitamin D receptor activation. In contrast, no key COPD pathways were identified that are associated with consistent actual upregulation and/or predicted activation of downstream inflammatory mediators.

Consistent with the canonical pathway functional enrichment analysis, in sputum the most common terms associated with the top 10 significant gene ontology (GO) biological processes and molecular mechanisms were immune response, programmed cell death and cytokine signalling for both doses (Table S2). Overall, 347 and 409 GO biological processes and 11 and 26 GO molecular functions were found to be significantly affected (pFDR< 0.05) by treatment with CHF6001 800 and 1600 μg BID, respectively.

Inflammatory gene interaction network analysis

Gene networks comprising all pro-inflammatory cytokines and matrix-metalloproteinases differentially expressed after treatment with CHF6001 relative to placebo (pFDR< 0.05 or p < 0.05 and |FC| > 1.3) are shown in Fig. 3a and b. CHF6001 800 μg BID differentially expressed 25 pro-inflammatory cytokines and matrix metalloproteinases (16 with pFDR< 0.05 and nine with p < 0.05 and |FC| > 1.3), 23 of which were downregulated (Table S3). CHF6001 1600 μg BID differentially expressed 33 pro-inflammatory cytokines and matrix metalloproteinases (25 pFDR< 0.05 and eight with p < 0.05 and |FC| > 1.3), 29 of which were downregulated (Table S4). Most of the differentially expressed genes for CHF6001 800 μg BID were also differentially expressed for CHF6001 1600 μg BID (19 out of 25 genes, 76%), with all common genes differentially regulated in the same direction. The gene with the most network interactions for both CHF6001 doses was tumour necrosis factor (TNF; 20 interactions with 800 μg BID and 28 with 1600 μg BID); notably this was the only inflammatory mediator with direct interaction with the PDE4 genes (Fig. S4).

Fig. 3

Network of pro-inflammatory cytokines and matrix metalloproteinases differentially expressed after treatment with CHF6001 A) 800 μg and B) 1600 μg BID relative to placebo. Each node represents all the proteins produced by a single, protein-coding gene locus, edges represent proteins that jointly contribute to a shared function, and the information inside the circle describes protein structure. Edges: a red line indicates the presence of fusion evidence; a green line, neighbourhood evidence; a blue line, co-occurrence evidence; a magenta line, experimental evidence; a yellow line, text mining evidence; a light blue line, database evidence; a black line, co-expression evidence; a purple line, protein homology evidence. A green halo around the nodes: downregulation, a red halo: upregulation. Unmarked nodes: pFDR< 0.05; * on the nodes: p < 0.05 and |FC| > 1.3

In addition to the network and pathway analyses, WikiPathways [14] was used to visualise the biological pathways affected by differentially regulated genes. In particular, the pathway “Cytokines and Inflammatory Response (Homo sapiens)” was used to match the genes that were differentially expressed with both doses (i.e. pFDR< 0.05 or p < 0.05 and |FC| > 1.3). Only downregulated genes matched the pathway (Fig. 4): One was downregulated with CHF6001 800 μg BID (IL-15), four with CHF6001 1600 μg BID (transforming growth factor β1, [TGFB1], colony stimulating factor 1 [CSF1], interferon γ and IL-1A), and three with both CHF6001 doses (TNF, platelet-derived growth factor alpha polypeptide [PDGFA], and IL-12B).

Fig. 4

Cytokines and inflammation response from WikiPathways. Downregulated genes are coloured in different shades of green

Notably, CHF6001 significantly reduced the expression of many inflammatory genes known to be involved in the pathophysiology of COPD (Fig. S5, Table S5, Fig. S6 and Table S6). In particular, both doses led to a significant (pFDR< 0.05) downregulation of genes coding for pro-inflammatory TNF superfamily members, interferon gamma receptor (IFNGR2) [15], purinergic receptor (P2RX7) [16], endothelin 1 (EDN1) [17, 18], complement system (C3) [19, 20] and the profibrotics SERPINE1 [21] and platelet-derived growth factor alpha and beta (PDGFA, PDGFB) [20]. Furthermore, both doses also led to significant upregulation of the gene coding for the anti-inflammatory interleukin 10 receptor (IL-10RB) [15] and of the suppressor of cytokine signalling 3 (SOCS3) [22], and significantly modulated genes toward a positive regulation of the vitamin D pathway [23,24,25,26,27] (upregulation of VDR, RXR, DHCR7) and toward a negative regulation of oxidative stress [23, 28, 29] (upregulation of SESN2, HP, CYGB and downregulation of ATG7, CD1B, NCF1), Th2 cytokine production [23, 30] (upregulation of SCGB1A1 and TNFRSF21), B cells activation/proliferation [23, 31] (upregulation of SAMSN1, INPP5D) and significant downregulation of the eosinophil major basic protein (PRG2) [32]. With both doses, genes associated with eosinophils and basophils were significantly downregulated (CLC, OLIG1, OLIG2 and PRSS33) [33]. Treatment with the highest dose also significantly downregulated genes coding for the pro-inflammatory and mucus producer mitogen-activated protein kinase 13 (MAPK13) [23, 34], the cytotoxicity natural killer cell triggering receptor NCR3 [23], significantly upregulated the metalloproteinase inhibitor TIMP1 [23, 35] and the neutrophil elastase inhibitor SERPINB1 [36], and significantly modulated genes toward positive regulation of ciliary function (IFT88, DNAI2, DNAI1) [23, 37, 38].


Here we show that CHF6001, on top of inhaled maintenance triple therapy, significantly modulated key pathophysiological inflammatory processes in sputum, specifically cytokine and matrix metalloprotease (MMP) activity, pathogen-associated pattern-recognition activity, oxidative stress and vitamin D (Fig. 5), demonstrating an add-on biological effect in the lung. In contrast, CHF6001 showed a non-significant pharmacodynamic effect in the blood. Mechanistically, these findings support the potential for CHF6001 to provide an additional beneficial effect in patients with chronic bronchitis who are still symptomatic despite regular use of ICS/LABA/LAMA, with a favourable systemic tolerability profile.

Fig. 5

Key pathophysiological pathways significantly modulated by CHF6001 in sputum with associated inhibition of downstream inflammatory effectors. Created with

The functional enrichment analysis in sputum cells highlighted an effect of CHF6001 on growth, development, inflammation, immune response and cytokine signalling processes. The majority of significant pathways modulated by CHF6001 are associated with the pathophysiology of COPD including cytokines, cAMP, Th1, Th2, dendritic cells, natural killer cells, oxidative stress, and vitamin D pathways. Furthermore, cytokines and MMPs play a major role in the regulation of immune response and inflammatory conditions [15, 39]. An analysis of the differentially regulated genes in our database clearly showed that the number of downregulated inflammatory cytokines and MMPs outweighed the number of those upregulated (92 and 88% of the differentially regulated genes were downregulated with CHF6001 800 and 1600 μg BID, respectively). In addition, the canonical pathway analysis showed significant modulation associated with actual downregulation and/or predicted inhibition of downstream inflammatory mediators for key inflammatory pathways involved in cytokine activity (specifically CCR5 [40], chemokine, high mobility group-B1 [41], TREM1 [42], IL-22 [15], ceramide [43] and PKCθ signalling [44]), pathogen-associated pattern-recognition activity (pattern recognition receptors for the recognition of bacteria [45, 46] and viruses and toll-like receptors signalling [47]), oxidative stress (fMLP signalling in neutrophils [48], and nitric oxide and reactive oxygen species production in macrophages [29]), and vitamin D signalling [25,26,27].

Importantly, analysis of canonical pathways was carried out using a mixed model of actual and/or predicted regulation to simulate the downstream consequences of molecules mediated by the treatment. This simulation showed a straightforward concordance with the actual downregulation of pro-inflammatory cytokines and MMPs (Fig. 3) which comprises most of the downstream effectors of each identified pathway (Fig. S3).

The trend toward decreased inflammation was apparent in the WikiPathway analyses, in which eight differentially expressed genes in our database matched the pro-inflammatory mediators included in the pathways, all of which were downregulated by both CHF6001 doses.

Among the pro-inflammatory cytokine genes significantly differentially expressed, IL-32, IL-12B, EBI3 (IL-27B protein), TNFAIP8 (tumour necrosis factor alpha-induced-protein-8), TNFSF15 (tumour necrosis factor ligand-superfamily-member-15), CCL3 (MIP-α) and CCL4 (MIP-1β) were downregulated with both CHF6001 doses. Importantly, IL-12, IL-32 and IL-27 induce a variety of inflammatory signalling, and are linked to IFN-γ which is the predominant cytokine produced by Th1 and Tc1 cells and plays a role in inflammation in COPD [15, 49, 50]. Notably the IFN-γ gene was downregulated by the lowest dose (p < 0.05, |FC| > 1.3), as were CXCL10 (p < 0.05, |FC| > 1.3) and SPP1 (osteopontin protein), which enhance production of IFN-γ and IL-12 [51, 52]. Furthermore, expression of CCL3 and CCL4 are elevated in the lungs of individuals with COPD, both of which code for cytokines that activate CCR5 receptors and contribute to the recruitment of macrophage and T cells into the airways [15]. Interestingly, CCR5 (significantly downregulated with the lower dose) and its ligand CCL5 (RANTES proteins; significantly downregulated with the higher dose) are increased in airways and sputum of patients with COPD during exacerbations [15]. Other key cytokine genes significantly downregulated by CHF6001 were CCL15 for the lower dose, which is a potent inflammatory cell chemotactic agent, and, for the high dose, IL-1A, TNF, IL-18, CSF1, TGFB1, LTA (lymphotoxin-alpha or tumour necrosis factor-beta protein), and LTB (lymphotoxin-beta or tumour necrosis factor C protein). Proteins coding for these genes are involved in chemotaxis, positive regulation of cytokine activity, chronic inflammation and oxidative stress [15, 53, 54]. CHF6001 also significantly downregulated a wide range of metalloprotease genes (MMP7, MMP12, MMP14, and, with the higher dose, MMP1) which cause morphological changes in the lungs and contribute significantly to the COPD state [25,26,27].

Since the unbiased analysis of all the differentially expressed inflammatory cytokines and MMPs in our dataset was used as a tool to probe the overall level of inflammation after treatment in sputum, it is important to recognise that the effect was not consistent for all the mediators. Indeed, the gene coding for the chemokine (C-X-C motif) ligand 6 (CXCL6), which is a chemoattractant for neutrophilic granulocytes, and the tumor necrosis factor ligand superfamily member 13B (TNFSF13B), which acts as a potent B cell activator, did not follow the same direction as the other pro-inflammatory mediators, being consistently upregulated with both doses of CHF6001.

Importantly, in this study we show a remarkable consistency between doses. Most of the probe sets significantly differentially expressed with the low dose were also differentially expressed with the high dose in the same direction, i.e. either both up-regulated or both down-regulated. In addition, we demonstrated a clear dose-response relationship both for the number of probe sets significantly modulated by CHF6001 and the differentially expressed fold change effect size.

Overall, the current results are in agreement with the analysis of biomarker data from these patients, showing that CHF6001 significantly decreased the levels of various target soluble inflammatory mediators in the sputum [10], namely leukotriene B4, C-X-C motif chemokine ligand 8, macrophage inflammatory protein 1β, matrix metalloproteinase 9, and tumour necrosis factor α (TNFα). The current analyses build on these findings, by providing a complete pharmacodynamic characterisation of the biological effect of CHF6001 in sputum cells and whole blood, covering nearly the whole genome coding for proteins. This approach provides an holistic picture of complex anti-inflammatory effects, which would be otherwise missed when measuring target biomarkers of inflammation at protein level.

To our knowledge this is the first ‘omics’ high-throughput placebo controlled biologic characterisation of an anti-inflammatory agent in COPD in the airways and systemic compartment. In addition, differentially regulated genes, pathways, biological processes and molecular functions were considered significant only upon correction for multiplicity (pFDR< 0.05). Few previous studies have investigated gene expression by microarray in blood after treatment of COPD patients with anti-inflammatory agents [33, 55, 56] and most of them did not correct for multiplicity [55, 56]. One study evaluated the effect of azithromycin by microarray both in blood and sputum of COPD patients, although the analyses were not corrected for placebo [57]. Taken together, these findings highlight the difficulties in generating significant results in complex transcriptomics studies, in which correction for multiplicity for thousands of concomitant assessments and the high variability makes obtaining statistically significant results prohibitive. To decrease variability, our trial was conducted in highly standardised conditions in a crossover design with an adequate washout to ensure there was no carry-over effect between the different periods. All clinical sites were highly skilled with blood and sputum collection and processing, leading to an average viability of sputum cells of 92.5%. In addition, the quality control analysis highlighted the high quality of RNA, with 100% of the blood samples and 99.3% of the sputum samples being acceptable. We acknowledge that the analyses have some limitations. A small number of matching pre-dose samples were replaced with the pre-dose samples from the closest available period, assuming there was no period effect on pre-dose expression. We recognise that this might increase uncertainty on the robustness of our results, but the validity of this assumption was supported by i) the lack of significant associations in the principal component analysis between pre-dose gene expression data and period/treatment; ii) the lack of significant differences between treatments in pre-dose expression for any of the analysed probe sets (see ANOVA analysis in supplement, results section); and iii) the high concordance of study results between different doses of CHF6001.


The current analyses demonstrate that in patients with chronic bronchitis receiving maintenance inhaled triple therapy, CHF6001 inhaled BID significantly modulates immune system processes in sputum but not in blood. Canonical pathway and network analyses highlight downregulation of key pro-inflammatory pathways and mediators involved in the pathophysiology of COPD in sputum. This is the profile that CHF6001 was designed to demonstrate: extrafine inhaled administration with high systemic clearance and high protein binding to minimise the pharmacodynamic effect in blood (and thus optimise systemic safety) and maximise the biological effect in the lung.

Availability of data and materials

Chiesi commits to sharing with qualified scientific and medical Researchers, conducting legitimate research, patient-level data, study-level data, the clinical protocol and the full clinical study report of Chiesi Farmaceutici S.p.A.-sponsored interventional clinical trials in patients for medicines and indications approved by the European Medicines Agency and/or the US Food and Drug Administration after 1st January 2015, following the approval of any received research proposal and the signature of a Data Sharing Agreement. Chiesi provides access to clinical trial information consistently with the principle of safeguarding commercially confidential information and patient privacy. To date, the current study is out of scope of the Chiesi policy on Clinical Data Sharing.

Other information on Chiesi’s data sharing commitment, access and research request’s approval process are available in the Clinical Trial Transparency section of


  1. 1.

    Rennard SI, Calverley PMA, Goehring UM, Bredenbröker D, Martinez FJ. Reduction of exacerbations by the PDE4 inhibitor roflumilast - the importance of defining different subsets of patients with COPD. Respir Res. 2011;12:18.

  2. 2.

    Martinez FJ, Rabe KF, Sethi S, Pizzichini E, McIvor A, Anzueto A, et al. Effect of roflumilast and inhaled corticosteroid/long-acting β2-agonist on chronic obstructive pulmonary disease exacerbations (RE2SPOND). A randomized clinical trial. Am. J. Respir. Crit. Care Med. 2016;194:559–67.

  3. 3.

    Martinez FJ, Calverley PMA, Goehring U-M, Brose M, Fabbri LM, Rabe KF. Effect of roflumilast on exacerbations in patients with severe chronic obstructive pulmonary disease uncontrolled by combination therapy (REACT): a multicentre randomised controlled trial. Lancet Elsevier. 2015;385:857–66.

  4. 4.

    Gupta S. Side-effects of roflumilast. Lancet. 2012;379:710–1.

  5. 5.

    Kardos P, Mokros I, Sauer R, Vogelmeier CF. Health status in patients with COPD treated with roflumilast: two large noninterventional real-life studies: DINO and DACOTA. Int. J. Chron. Obstruct. Pulmon. Dis. 2018;13:1455–68.

  6. 6.

    Mariotti F, Govoni M, Lucci G, Santoro D, Nandeuil MA. Safety, tolerability, and pharmacokinetics of single and repeat ascending doses of CHF6001, a novel inhaled phosphodiesterase-4 inhibitor: two randomized trials in healthy volunteers. Int. J. Chron. Obstruct. Pulmon. Dis. Dove Press. 2018;13:3399–410.

  7. 7.

    Cenacchi V, Battaglia R, Cinato F, Riccardi B, Spinabelli D, Brogin G, et al. In vitro and in vivo metabolism of CHF 6001, a selective phosphodiesterase (PDE4) inhibitor. Xenobiotica. 2015;45:693–710.

  8. 8.

    Villetti G, Carnini C, Battipaglia L, Preynat L, Bolzoni PT, Bassani F, et al. CHF6001 II: a novel phosphodiesterase 4 inhibitor, suitable for topical pulmonary administration - in vivo preclinical pharmacology profile defines a potent anti-inflammatory compound with a wide therapeutic window. J Pharmacol Exp Ther. 2015;352:568–78.

  9. 9.

    Singh D, Nandeuil MA, Pigeon-Francisco C, Emirova A, Santoro D, Biondaro S, et al. Efficacy and safety of CHF6001, a novel inhaled PDE4 inhibitor in COPD: the PIONEER dose finding study. Am J Respir Crit Care Med. 2019;199:A4529.

  10. 10.

    Singh D, Beeh KM, Colgan B, Kornmann O, Leaker B, Watz H, et al. Effect of the inhaled PDE4 inhibitor CHF6001 on biomarkers of inflammation in COPD. Respir Res. 2019;20:180.

  11. 11.

    Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57:289–300.

  12. 12.

    Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43:D447–52.

  13. 13.

    Ezkurdia I, Juan D, Rodriguez JM, Frankish A, Diekhans M, Harrow J, et al. Multiple evidence strands suggest that there may be as few as 19,000 human protein-coding genes. Hum Mol Genet. 2014;23:5866–78.

  14. 14.

    Slenter DN, Kutmon M, Hanspers K, Riutta A, Windsor J, Nunes N, et al. WikiPathways: a multifaceted pathway database bridging metabolomics to other omics research. Nucleic Acids Res. 2018;46:D661–7.

  15. 15.

    Barnes PJ. The cytokine network in asthma and chronic obstructive pulmonary disease. J Clin Invest. 2008;118:3546–56.

  16. 16.

    Arulkumaran N, Unwin RJ, Tam FW. A potential therapeutic role for P2X7 receptor (P2X7R) antagonists in the treatment of inflammatory diseases. Expert Opin Investig Drugs. 2011;20:897–915.

  17. 17.

    Fagan KA, McMurtry IF, Rodman DM. Role of endothelin-1 in lung disease. Respir Res. 2001;2:90–101.

  18. 18.

    Roland M, Bhowmik A, Sapsford RJ, Seemungal TA, Jeffries DJ, Warner TD, et al. Sputum and plasma endothelin-1 levels in exacerbations of chronic obstructive pulmonary disease. Thorax. 2001;56:30–5.

  19. 19.

    Westwood J-P, Mackay AJ, Donaldson G, Machin SJ, Wedzicha JA, Scully M. The role of complement activation in COPD exacerbation recovery. ERJ Open Research. 2016;2:00027–2016.

  20. 20.

    Barnes PJ, Drazen JM, Rennard SI, Thomson NC, editors. Asthma and COPD: Basic mechanisms and clinical management. 2nd ed. [Oxford, UK]: Academic Press; 2008.

  21. 21.

    Liu R-M. Oxidative stress, plasminogen activator inhibitor 1, and lung fibrosis. Antioxid Redox Signal. 2008;10:303–19.

  22. 22.

    Springer J, Scholz FR, Peiser C, Dinh QT, Fischer A, Quarcoo D, et al. Transcriptional down-regulation of suppressor of cytokine signaling (SOCS)-3 in chronic obstructive pulmonary disease. J Occup Med Toxicol. 2013;8:29.

  23. 23.

    UniProt Consortium. UniProt data [Internet]. [cited 2019 Mar 5]. Available from:

  24. 24.

    Prabhu AV, Luu W, Li D, Sharpe LJ, Brown AJ. DHCR7: a vital enzyme switch between cholesterol and vitamin D production. Prog. Lipid Res. 2016;64:138–51.

  25. 25.

    Sundar IK, Hwang J-W, Wu S, Sun J, Rahman I. Deletion of vitamin D receptor leads to premature emphysema/COPD by increased matrix metalloproteinases and lymphoid aggregates formation. Biochem Biophys Res Commun. 2011;406:127–33.

  26. 26.

    Ishii M, Yamaguchi Y, Nakamura T, Akishita M. The vitamin D receptors may function as antiinflammatory effects in patients with COPD. Chest. 2015;148:690.

  27. 27.

    García de Tena J, El Hachem Debek A, Hernández Gutiérrez C, Izquierdo Alonso JL. The role of vitamin D in chronic obstructive pulmonary disease, asthma and other respiratory diseases. Arch. Bronconeumol. 2014;50:179–84.

  28. 28.

    Gerritsen WBM, Asin J, Zanen P, Bosch JMM va. den, Haas FJLM. Markers of inflammation and oxidative stress in exacerbated chronic obstructive pulmonary disease patients. Respir. Med. 2005;99:84–90. doi:

  29. 29.

    Domej W, Oettl K, Renner W. Oxidative stress and free radicals in COPD--implications and relevance for treatment. Int. J. Chron. Obstruct. Pulmon. Dis. 2014;9:1207–24.

  30. 30.

    Barnes PJ. Inflammatory mechanisms in patients with chronic obstructive pulmonary disease. J Allergy Clin Immunol. 2016;138:16–27.

  31. 31.

    Nixon J, Newbold P, Mustelin T, Anderson GP, Kolbeck R. Monoclonal antibody therapy for the treatment of asthma and chronic obstructive pulmonary disease with eosinophilic inflammation. Pharmacol Ther. 2017;169:57–77.

  32. 32.

    George L, Brightling CE. Eosinophilic airway inflammation: role in asthma and chronic obstructive pulmonary disease. Ther Adv Chronic Dis. 2016;7:34–51.

  33. 33.

    Sridhar S, Liu H, Pham T-H, Damera G, Newbold P. Modulation of blood inflammatory markers by benralizumab in patients with eosinophilic airway diseases. Respir Res. 2019;20:14.

  34. 34.

    Alevy YG, Patel AC, Romero AG, Patel DA, Tucker J, Roswit WT, et al. IL-13-induced airway mucus production is attenuated by MAPK13 inhibition. J Clin Invest. 2012;122:4555–68.

  35. 35.

    Beeh KM, Beier J, Kornmann O, Buhl R. Sputum matrix metalloproteinase-9, tissue inhibitor of metalloprotinease-1, and their molar ratio in patients with chronic obstructive pulmonary disease, idiopathic pulmonary fibrosis and healthy subjects. Respir Med. 2003;97:634–9.

  36. 36.

    Groutas WC, Dou D, Alliston KR. Neutrophil elastase inhibitors. Expert Opin Ther Pat. 2011;21:339–54.

  37. 37.

    Hessel J, Heldrich J, Fuller J, Staudt MR, Radisch S, Hollmann C, et al. Intraflagellar transport gene expression associated with short cilia in smoking and COPD. PLoS One. 2014;9:e85453.

  38. 38.

    Yaghi A, Dolovich MB. Airway epithelial cell cilia and obstructive lung disease. Cells. 2016;5:40.

  39. 39.

    Srivastava PK, Dastidar SG, Ray A. Chronic obstructive pulmonary disease: role of matrix metalloproteases and future challenges of drug therapy. Expert Opin Investig Drugs. 2007;16:1069–78.

  40. 40.

    Costa C, Traves SL, Tudhope SJ, Fenwick PS, Belchamber KBR, Russell REK, et al. Enhanced monocyte migration to CXCR3 and CCR5 chemokines in COPD. Eur Respir J. 2016;47:1093–102.

  41. 41.

    Gangemi S, Casciaro M, Trapani G, Quartuccio S, Navarra M, Pioggia G, et al. Association between HMGB1 and COPD: a systematic review. Mediat Inflamm. 2015;2015:164913.

  42. 42.

    Radsak MP, Taube C, Haselmayer P, Tenzer S, Salih HR, Wiewrodt R, et al. Soluble triggering receptor expressed on myeloid cells 1 is released in patients with stable chronic obstructive pulmonary disease. Clin Dev Immunol. 2007;2007:52040.

  43. 43.

    Petrusca DN, Gu Y, Adamowicz JJ, Rush NI, Hubbard WC, Smith PA, et al. Sphingolipid-mediated inhibition of apoptotic cell clearance by alveolar macrophages. J Biol Chem. 2010;285:40322–32.

  44. 44.

    Chand S, Mehta N, Bahia MS, Dixit A, Silakari O. Protein kinase C-theta inhibitors: a novel therapy for inflammatory disorders. Curr Pharm Des. 2012;18:4725–46.

  45. 45.

    Pouwels SD, van Geffen WH, Jonker MR, Kerstjens HAM, Nawijn MC, Heijink IH. Increased neutrophil expression of pattern recognition receptors during COPD exacerbations. Respirology. 2017;22:401–4.

  46. 46.

    Kinose D, Ogawa E, Kudo M, Marumo S, Kiyokawa H, Hoshino Y, et al. Association of COPD exacerbation frequency with gene expression of pattern recognition receptors in inflammatory cells in induced sputum. Clin Respir J. 2016;10:11–21.

  47. 47.

    Pomerenke A, Lea SR, Herrick S, Lindsay MA, Singh D. Characterization of TLR-induced inflammatory responses in COPD and control lung tissue explants. Int J Chron Obstruct Pulmon Dis. 2016;11:2409–17.

  48. 48.

    Yoshikawa T, Dent G, Ward J, Angco G, Nong G, Nomura N, et al. Impaired neutrophil chemotaxis in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2007;175:473–9.

  49. 49.

    Calabrese F, Baraldo S, Bazzan E, Lunardi F, Rea F, Maestrelli P, et al. IL-32, a novel proinflammatory cytokine in chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2008;178:894–901.

  50. 50.

    Cao J, Zhang L, Li D, Xu F, Huang S, Xiang Y, et al. IL-27 is elevated in patients with COPD and patients with pulmonary TB and induces human bronchial epithelial cells to produce CXCL10. Chest. 2012;141:121–30.

  51. 51.

    Vasquez RE, Xin L, Soong L. Effects of CXCL10 on dendritic cell and CD4 T-cell functions during Leishmania amazonensis infection. Infect Immun. 2008;76:161–9.

  52. 52.

    Li X, O’Regan AW, Berman JS. IFN- γ induction of osteopontin expression in human monocytoid cells. J Interf Cytokine Res. 2004;23:259–65.

  53. 53.

    Vlahos R, Bozinovski S. Role of alveolar macrophages in chronic obstructive pulmonary disease. Front Immunol. 2014;5:435.

  54. 54.

    John-Schuster G, Conlon TM, Gunter S, Eickelberg O, Heikenwalder M, Yildirim AO. Lymphotoxin-β-receptor blockade prevents cigarette smoke-induced COPD. Am J Respir Crit Care Med. 2017;195:A4297.

  55. 55.

    Betts JC, Mayer RJ, Tal-singer R, Warnock L, Clayton C, Bates S, et al. Gene expression changes caused by the p38 MAPK inhibitor dilmapimod in COPD patients : analysis of blood and sputum samples from a randomized, placebo-controlled clinical trial. Pharmacol Res Perspect. 2015;3:e00094.

  56. 56.

    van den Berge M, Steiling K, Timens W, Hiemstra PS, Sterk PJ, Heijink IH, et al. Airway gene expression in COPD is dynamic with inhaled corticosteroid treatment and reflects biological pathways associated with disease activity. Thorax. 2014;69:14–23.

  57. 57.

    Baines KJ, Wright TK, Gibson PG, Powell H, Hansbro PM, Simpson JL. Azithromycin treatment modifies airway and blood gene expression networks in neutrophilic COPD. ERJ Open Research. 2018;4:00031–2018.

Download references


Results of this study have been presented at the European Respiratory Society 2019 annual conference.

The authors would like to thank the investigators and patients at the investigative sites for their support of this study. Dave Singh is supported by the National Institute for Health Research (NIHR) Manchester Biomedical Research Centre (BRC).


This study was funded by Chiesi Farmaceutici SpA. Chiesi Farmaceutici SpA, was responsible for the design and analysis of the study, oversaw its conduct and was responsible for study report preparation. Writing support was provided by David Young of Young Medical Communications and Consulting Ltd. This support was funded by Chiesi Farmaceutici SpA.

Author information

The study was conceived by MG, HW and DS. BC, BL, OK, KMB, HW and DS contributed to data acquisition. Data were analysed by MG, MB, SV, GLJ and EKA. All authors interpreted the data, and revised the manuscript for intellectual content and approved the submitted version.

Correspondence to Mirco Govoni.

Ethics declarations

Ethics approval and consent to participate

The study was approved by independent ethics committees for each institution (see the supplement for details). All patients provided written informed consent prior to study start.

Consent for publication

Not applicable.

Competing interests

MG, MB, SV, GL, AE, MAN and SP are employees of Chiesi, the sponsor of this study.

GLJ and EKA are employees of the Almac group; Almac received funding from Chiesi as part of a consultation fee for the work reported in this manuscript.

BC and BL have nothing to disclose.

OK’s institution received fees from Chiesi for conducting this study as a participating site. He reports personal fees as speaker or Advisory Board member from AstraZeneca, Boehringer Ingelheim, GlaxoSmithKline, Sanofi and Novartis.

KMB declares that no personal payments were received from any pharmaceutical entity in the past 5 yrs. He is a full time employee of insaf Respiratory Research Institute. The institution has received compensation for services on advisory boards or consulting for Ablynx, Almirall, AstraZeneca, Berlin Chemie, Boehringer, Chiesi, Cytos, Mundipharma, Novartis, Pohl Boskamp, Zentiva. The institution has received compensation for speaker activities in scientific meetings supported by Almirall, AstraZeneca, Berlin Chemie, Boehringer, Cytos, ERT, GSK, Novartis, Pfizer, Pohl Boskamp, Takeda. The institution has received compensation for design and performance of clinical trials from Almirall, Altana/Nycomed, AstraZeneca, Boehringer, Cytos, GSK, Infinity, Medapharma, MSD, Mundipharma, Novartis, Parexel, Pearl Therapeutics, Pfizer, Revotar, Teva, Sterna, and Zentiva.

HW reports personal fees from Chiesi during the conduct of the study. Outside the submitted work, Dr. Watz reports personal fees from Bayer, personal fees from GSK, personal fees from Boehringer Ingelheim, personal fees from Novartis, personal fees from AstraZeneca, personal fees from BerlinChemie, and personal fees from Roche.

DS received personal fees from Chiesi during the conduct of this study. Outside the submitted work, he reports grants and personal fees from AstraZeneca, Boehringer Ingelheim, Chiesi, GlaxoSmithKline, Glenmark, Menarini, Mundipharma, Novartis, Pfizer, Pulmatrix, Therevance, and Verona, and personal fees from Cipla, Genentech and Peptinnovate.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Govoni, M., Bassi, M., Vezzoli, S. et al. Sputum and blood transcriptomics characterisation of the inhaled PDE4 inhibitor CHF6001 on top of triple therapy in patients with chronic bronchitis. Respir Res 21, 72 (2020).

Download citation


  • Gene expression
  • Inflammation
  • Phosphodiesterase 4 inhibitors
  • Biomarkers
  • Chronic obstructive pulmonary disease