Bronchial Thermoplasty in severe asthma: a real-world study on efficacy and gene profiling

doi:10.21203/rs.3.rs-1311690/v1

Download PDF

Research

Bronchial Thermoplasty in severe asthma: a real-world study on efficacy and gene profiling

https://doi.org/10.21203/rs.3.rs-1311690/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background: Bronchial Thermoplasty (BT) is an effective treatment in severe asthma. How to select patients who more likely benefit from BT is an unmet clinical need. Moreover, mechanisms of BT efficacy are still largely unknown. We sought to determine BT efficacy and to identify potential mechanisms of response.

Methods: This retrospective cohort study evaluated clinical outcomes in 27 patients with severe asthma: 13 with T2-high and 14 with T2-low endotype. Expression levels of 20 genes were compared by real-time PCR in bronchial biopsies performed at the third BT session versus baseline. Clinical response was measured based on Asthma Control Questionnaire (ACQ) score < 1.5, asthma exacerbations < 2, oral corticosteroids reduction of at least 50% at 12 months post-BT. Patients were classified as responders when they had at least 2 of 3 outcome measures.

Results: 81% of patients were defined as responders. BT induced a reduction in alpha smooth muscle actin (ACTA2) and an increase in CD68, fibroblast activation protein-alpha (FAP), alpha-1 and alpha-2 type I collagen (COL1A1, COL1A2) gene expression in the majority of patients. A higher reduction in ubiquitin carboxy-terminal-hydrolase L1 (PGP9.5) mRNA correlated with a better response based on Asthma Quality of Life Questionnaire (AQLQ). Lower changes in CD68 and FAP mRNAs correlated with a better response based on ACQ. Lower levels of occludin (OCLN), CD68, connective tissue growth factor (CTGF), higher levels of secretory leukocyte protease inhibitor (SLPI) and lower changes in CD68 and CTGF mRNAs were observed in patients who had less than 2 exacerbations post-BT. Lower levels of COL1A2 at baseline were observed in patients who had ACQ < 1.5 at 12 months post-BT.

Conclusions: BT is effective irrespective of the asthma endotypes and seems associated with airway remodelling. Quantification of OCLN, CD68, CTGF, SLPI, COL1A2 mRNAs could be useful to identify patients with better results.

Trial registration: The study protocol was approved by the Local Ethics Committee (Azienda USL-IRCCS of Reggio Emilia - Comitato Etico Area Vasta Nord of Emilia Romagna; protocol number: 2019/0014076) and all the patients provided written informed consent before participating in the study.

Bronchial thermoplasty

severe asthma

gene expression

real-time PCR

smooth muscle

bronchial biopsies

The heterogeneous and inflammatory background of severe asthma is today richly documented [1]. High-throughput analyses contributed to better define the current and well known multiple endotypes and phenotypes of severe refractory asthma (SRA) [2]. However, it is still debated how to choose the best therapeutic options for the right asthma phenotype. Immunologic profile includes two severe asthma endotypes: T helper type 2 cell high (T2-high), associated to serum biomarkers as eosinophils and immunoglobulin E (IgE), and T helper type 2 cell low (T2-low) characterized by neutrophilic or paucigranulocytic inflammatory pattern and airway damage associated with mucus gland hyperplasia and hyper secretion, airway hyperreactivity, remodelling, and corticosteroid insensitivity [3].

Bronchial Thermoplasty (BT) is considered an effective treatment in the T2-low endotype [4], however it could be a therapeutic option for all patients whose symptoms are not controlled with standard of care and despite monoclonal antibody treatments [5]. BT is an endoscopic procedure based on the local delivery of radio frequency at 65°C to the airways with a device called Alair™ Catheter (Boston Scientific, Natick, MA, USA) [6]. Several clinical trials have demonstrated BT effectiveness especially on clinical outcomes, particularly in the reduction of exacerbations [7–9]. Long-term effects of BT have been documented by observational studies [10].

However, the discussion concerning effectiveness of BT and its mechanism of action remains open [11]. These include changes in function and structure of airway smooth muscles, nerve fibers, inflammatory cells/mediators and components of the extracellular matrix [12–17]. Furthermore, a key aspect yet to be clarified is the identification of potential predictive response factors [18], to identify those patients who more likely benefit from BT. An operator-dependent variable that can influence response may be the number of activations delivered during BT sessions, which appears to play an important role in determining favorable clinical outcomes [19].

The aims of this study were to determine efficacy of BT, to identify genes modulated by BT and associations between gene expression levels in bronchial biopsy specimens and clinical outcomes. We investigated the expression of 20 genes, selected following a hypothesis-driven approach starting from and expanding literature data [6, 13–17, 20–24]. In particular, we analyzed expression of genes reflecting airway remodelling: alpha smooth muscle actin (ACTA2) marker of airway smooth muscle cells, cadherin 1 (CDH1) marker of epithelial cells, connective tissue growth factor (CTGF), fibroblast activation protein alpha (FAP), alpha-1 and alpha-2 chain type I collagen (COL1A1, COL1A2), galectin-3 (LGAL3) markers of stroma. We investigated expression of genes mirroring airway inflammation: CD45, CD68, interleukin (IL)-4, IL-5, IL-6, IL-13, IL-17. We studied expression of innervation markers: ubiquitin carboxy-terminal hydrolase L1 (UCHL1/PGP9.5) and occludin (OCLN) as well as genes reported to be associated with asthma: periostin (POSTN), galectin-3 (LGAL3), secretory leukocyte protease inhibitor (SLPI), IL-4, IL-5, IL-13, IL-17. We investigated genes coding high temperature receptors (transient receptor potential cation channel subfamily V member 1 and member 2: TRPV1 and TRPV2) [22, 23] because BT is based on heat delivery. Finally we analyzed the expression of SLPI which has emerged as the most promising potential marker of response to BT in the work by Ano S et al [24].

Study population

A total of 27 adult patients who underwent BT at Pulmonology Unit of the Azienda USL-IRCCS of Reggio Emilia were enrolled. Patients were affected by severe refractory asthma according to ERS/ATS classification [25].

The study protocol was approved by the Local Ethics Committee (Azienda USL-IRCCS of Reggio Emilia - Comitato Etico Area Vasta Nord of Emilia Romagna; protocol number: 2019/0014076) and all the patients provided written informed consent before participating in the study.

Clinical data, blood biomarkers, asthma medications, patient reported outcomes in terms of asthma quality of life questionnaire (AQLQ) and asthma control questionnaire (ACQ) were collected. Data were assessed at baseline (before BT, T0), at 2 months during the last BT session (T2) and 12 months after the last BT session (T12).

On the basis of the cut-off for anti-IL-5 therapy (blood eosinophils count ≥ 300 cells/µl), patients were classified in the T2-high and T2-low endotypes. The first was defined as eosinophilic (blood eosinophils count ≥ 300 cells/µl), eosinophilic and allergic (blood eosinophils count ≥ 300 cells/µl, IgE ≥ 76 kU/l, perennial inhalant allergens), allergic (IgE ≥ 76 kU/l and perennial inhalant allergens). The T2-low endotype was defined by blood eosinophil count < 300 cells/ µl.

Measures to determine BT response

Outcome measures included: 1) ACQ score < 1.5 at T12; 2) < 2 asthma exacerbations at T12 versus T2; 3) oral corticosteroids (OCSs) reduction of at least 50% at T12 versus baseline. Patients were classified as responders when they had at least 2 of the 3 outcome measures. Patients who did not receive OCSs at baseline were classified as responders when they had one of the other outcome measures.

BT procedure

BT was performed with the device called Alair™ Catheter (Boston Scientific, Natick, MA, USA). All patients underwent 3 sessions of BT which consisted in the treatment of right lower lobe (first session), left lower lobe (second session) and right and left upper lobes (third session). Sessions were performed every four weeks. No heat activations were released in the middle lobe [13].

Airway sample collection

Bronchial biopsy specimens were performed in the left lower lobe at baseline (as T0) and in the right lower lobe after two months during the third BT session (T2) which was administered in the upper lobes. Bronchial biopsies were collected before heat delivery. At least two samples were embedded in Tissue Tek II OCT (Miles Scientific, Naperville, IL, USA), frozen within 15 minutes at -20°C, and finally stored at –80°C.

RNA extraction and cDNA preparation

Biopsies were recovered from thawed OCT and RNA was extracted with the RNA/DNA/Protein Purification Plus Kit (Norgen Biotek). Tissue homogenization was made effective by the use of pestles. NanoDrop instrument (Thermo Scientific) was used for RNA quantification, then 62.5 ng of RNA were reverse transcribed with the PrimeScrip RT Reagent Kit with gDNA Eraser (Takara) according to the manufacturer's instructions. cDNA was diluted 1:2 with DNase and RNase-free water and stored at -20° C.

Real-time PCR

Real-time PCR was performed with the SYBR Premix Ex Taq II (Tli RNaseH Plus) containing the ROX Reference Dye (Takara) and QuantiTect primer assays (Qiagen). Catalog numbers of the primers are listed in Supplementary Table 1. cDNA (1 µL) was amplified in 25 µL PCR mix. Amplification products were checked by dissociation curve analysis. Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as housekeeper gene for data normalization. Real-time PCR results are shown as normalized expression: 2-ΔCt (ΔCt = Ct mRNA target – Ct GAPDH) so that they might be comparable among different laboratories. Amplification was carried out with the 7300 Applied Biosystem instrument (Thermofisher).

Statistical analysis

Statistical analysis was performed with GraphPad Prism 6 software using non-parametric tests. For comparison between two groups, Wilcoxon test was applied for paired data, while Mann-Whitney U test was applied for unpaired data. For comparison between three groups, Kruskal-Wallis test with Dunns correction was used. Spearman test was used to determine correlations between two variables. Categorical values were compared with Fisher exact test. P values < 0.05 were considered statistically significant.

Demographic results at baseline

A total of 27 patients with severe refractory asthma according to ERS/ATS 2014 guidelines and GINA 2018 classification underwent BT from 2012 to 2018. Demographic, clinical features and comorbidities of the cohort of patients are shown in Table 1.

Table 1

Demographic and clinical characteristics
Demographic and Clinical Characteristics	T0
No. of patients	27
Age y, mean ± SD	52 ± 10
Female, n (%)	14 (52)
Age onset y, mean ± SD	36 ± 12
Ex smokers, n (%)	5 (18)
Atopics, n (%)	6 (22)
Switch from Omalizumab, n (%)	5 (18)
Phenotypes
T2-high, n (%)	13 (48)
Eosinophilic, n (%)	7 (26)
Eosinophilic and Allergic, n (%)	3 (11)
Only Allergic, n (%)	3 (11)
T2-low, n (%)	14 (52)
Neutrophilic, n (%)	14 (52)
Comorbidities
Gastroesophageal reflux disease (GERD), n (%)	14 (52)
Chronic rhinosinusitis with nasal polyps (CRSwNP), n (%)	9 (33)
BT Activations mean (SD)	193 ± 39
GERD gastro-esophageal reflux disease, CRSwNP chronic rhinosinusitis with nasal polyposis

BT clinical outcomes

According to outcome measures, 81% of the patients were defined as responders to BT treatment (Table 2). In particular, 77% of the patients with the T2-high endotype and 86% with the T2-low endotype were responders.

Table 2 Clinical criteria for definition of patients with a better response to BT

Clinical and airway functional parameters before and 12 months after BT are described in Table 3. Significant improvements in patients reported outcomes (ACQ and AQLQ scores), rate of exacerbations and airway functional data were observed (Table 3). The overall consumption of OCSs scaled down from 20 mg to 5 mg daily after 1 year (p<0.0001). At baseline, 89% of patients received OCSs and among these, 40% stopped OCSs at T12 and 55% reduced OCS daily dose by 50% (p< 0.0058).

Table 3

Clinical characteristics of the patients before and after BT
Parameter	T0	T12	p-value
No. of patients	27	27
Total serum IgE KU/l, median [IQR]	65 [30-160]	87 [27-272]	0.2552
Blood eosinophil count cells/µl, median [IQR]	140 [20-410]	340 [120-640]	0.1533
Lung function, pre-bronchodilator
FEV₁%,median [IQR]	71 [63-94]	67 [59-87]	0.0179
FEV₁ l, median [IQR]	2.43 [1.86-3.22]	1.97 [1.5-2.68]	0.0103
FVC %, median [IQR]	99 [89-114]	86 [77-102]	0.0011
FVC l, median [IQR]	3.62 [2.90-4.62]	3.19 [2.57-3.89]	0.0051
Asthma concomitant medication
Oral corticosteroid-prednisone or the equivalent mg, median [IQR]	20 [5–25]	5 [0-12.5]	< 0.0001
Patients taking OCS, n (%)	24 (89)	15 (55)	0.0058
Patients taking LAMA, n (%)	12 (44)	21 (77)	0.0244
Patients taking LTRA, n (%)	10 (37)	6 (22)	1
Patients taking theophylline, n (%)	4 (15)	0 (0)	0.1110
*Assigned adapted GINA step**
Step 4, n (%)	2 (7)	3 (11)	1
Step 5, n (%)	25 (93)	24 (89)	1
Asthma control
Exacerbations/y, median [IQR]	6 [3–10]	2 [1–3]	< 0.0001
ED visits/hospitalizations/y, median [IQR]	0 [0-1]	0 [0-0]	0.0022
Patient reported outcomes
AQLQ score, median [IQR]	2.6 [2.1-3.0]	5.5 [4.4-6.0]	< 0.0001
ACQ score, median [IQR]	3.8 [3.2-4.6]	1.6 [1.2-2.0]	< 0.0001
*GINA 2019 asthma guidelines
Data were compared with Mann Withney U test (continuous variables) and Fisher exact test (categorical variables).
P < 0.05 were considered statistically significant. IQR interquartile range

Characterization of bronchial gene expression at baseline

IL-4, IL-13 and IL-17 mRNAs were not detected in bronchial biopsies at baseline. IL-5 was detected in 9/27 samples while TRPV2 was detected in 21/27 samples with low expression (median Ct of IL-5: 34.2, IQR: 33.4-34.4; median Ct of TRPV2: 32.9, IQR: 31.3-33.9). The other investigated transcripts were expressed in all the biopsies, showing variable expression levels among different patients. Coefficient of variations (CVs) ranged from 42.1–168.0% (Supplementary Table 2). The genes that showed more variability (CV>100%) were FAP and IL-6, while those which showed less variability (CV< 50%) were CD68, LGAL3, CDH1.

Since patients were receiving OCSs, we determined if gene expression levels were associated with the dose of OCSs. Patients were divided in three groups based on OCS dose: ≤ 5 mg; 5–25 mg; ≥ 25 mg. OCSs did not affect the expression of the investigated genes (Supplementary Figure S1).

No differences were observed in gene expression levels when patients were stratified in endotypes (data not shown).

Effects of BT treatment on gene expression

To determine the effects of BT treatment, gene expression levels were compared in bronchial biopsies performed during the last BT session (T2) versus baseline. Patients showed to be heterogeneous regarding changes in gene expression induced by BT treatment (Figure 1). Overall, ACTA2 showed a significant reduction and CD68, FAP, COL1A1 and COL1A2 showed a significant increase at T2 versus T0 (P value < 0.05, Figure 1). No significant changes were observed in the other markers analyzed. Supplementary Table 3 reports the fraction of patients who had changes greater or equal to two-fold. Specifically, 52% of patients showed at least a two-fold decrease in ACTA2 while 70%, 89%, 63% of patients showed at least a two-fold increase in FAP, COL1A1, COL1A2 respectively. Changes in CD68 gene expression were mainly lower than two-fold, although statistically significant.

Correlations between changes in gene expression induced by BT treatment and differences in patient reported outcomes

To understand which alterations induced by BT treatment could be linked with patients reported outcomes, changes of gene expression in bronchial biopsies were correlated with changes in the AQLQ and ACQ scores at T2 versus T0. A better response to treatment in terms of ΔACQ was associated with lower changes in CD68 (P = 0.0154; r = 0.5651; 95% CI = 0.2160 to 0.7861) and FAP (P = 0.0126; r = 0.4736; 95% CI = 0.1021 to 0.7290) (Figure 2A). Instead, a higher reduction of PGP9.5 gene expression correlated with a better response to treatment in terms of ΔAQLQ (P = 0.0033; r = -0.5451; 95% CI = -0.7712 to -0.1969) (Figure 2B). No correlations were observed between changes in the expression of the other genes and changes in AQLQ or ACQ (Supplementary Figure S2 and S3).

Association between gene expression and clinical parameters 12 months post-BT

In our cohort only 5 patients were classified as non-responders to BT, thus it was not feasible to stratify patients in responders (n = 22) versus non-responders (n = 5) to find associations between clinical and molecular data. However, we determined if there were associations between gene expression levels at T0 and T2 and the clinical parameters used to evaluate BT efficacy: number of exacerbations in the 12 months post-BT and ACQ at T12. It was not feasible to consider OCS reduction because only two patients did not reduce OCSs by at least 50%.

Based on the number of exacerbations, statistically significant lower levels of OCLN in bronchial biopsies at T0; lower levels of OCLN, CD68, CTGF and higher levels of SLPI in bronchial biopsies at T2 were observed in patients who had less than 2 exacerbations (p < 0.05, Figure 3). Supplementary Figure S4 and S5 shows the expression of all the investigated genes at T0 and T2 grouping patients according to the numbers of exacerbations.

Based on ACQ at T12, statistically significant lower levels of COL1A2 in bronchial biopsies at T0 were observed in patients with ACQ < 1.5 (considered as controlled asthma) (Supplementary Figure S6). Instead no differences were found analyzing gene expression at T2 (data not shown).

Finally, classification of patients according to the fold changes in gene expression at T2 versus T0 revealed that lower changes in CD68 and CTGF mRNAs were observed in patients who had less than 2 exacerbations post-BT (Figure 4 and Supplementary Figure S7 to evaluate all the investigated genes). Instead no differences in the fold changes were found stratifying patients on ACQ values at T12 (data not shown).

This study is divided in two different parts: clinical effects of BT and evaluation of gene expression levels associated to clinical outcomes. Our results confirmed the efficacy of BT in the treatment of severe asthma, as we observed 81% of responders, in line with recent results published by Langton’s group [26].

Notably, our research was conducted in real life and it follows that patients who underwent BT had more severe asthma than patients included in earlier trials in term of rate of exacerbations (6 vs 0.7), AQLQ score (2.5 vs 4.7) and OCS maintenance dose (89% vs 41%) [8, 9].

Regarding clinical outcomes, a significant improvement in asthma symptoms in terms of ACQ scores, AQLQ scores and rate of exacerbations was observed. These results confirm the clinical benefits of BT in line with previous researches [20, 21, 26]. AQLQ improvement at 12 months was statistically significant (p< 0.0001) which is different from what was presented in AIR2 Trial, where AQLQ was statistically borderline and questioned [12, 25, 27–29]. Moreover, we observed a significant reduction in maintenance doses of OCS.

In the scenario of multiple therapeutic approaches for severe asthma treatments, BT still remains the most controversial. One of the biggest objections is the absence of a specific biomarker able to define target population for BT [12]. Currently, BT is used primary in non-eosinophilic phenotype, according to guidelines as there is no usable biological therapy [25]. In our site, BT started in 2012, when the only available biological treatment was omalizumab. Based on clinical criteria for responder definition described in Table 2, we observed 86% of responders in the T2-low cohort and 77% of responders in the T2-high cohort. This finding confirms that BT is a good therapeutic option in T2-low asthma patients. Moreover, it may be effective also in T2-high asthma patients (18% of patients were switched from omalizumab). Thus, our data emerge in contrast with the conception on BT being a treatment only in case of T2-low endotypes or when biological drugs fails [25]. In the real practice, BT represents a therapeutic option for patients with severe asthma despite of asthma endotypes.

The study expands the knowledge regarding the mechanisms of action of BT using gene expression profiling in bronchial biopsy samples. Genes to be investigated were selected through a hypothesis-driven approach based on and expanding data available in the literature [6, 13–17, 20–24].

BT induced a reduction of ACTA2 and an increase of CD68, FAP, COL1A1, COL1A2 mRNAs in the majority of patients suggesting tissue remodelling with a likely decrease in smooth muscle cells and an increase in monocytes/macrophages/dendritic cells and connective tissue. However, these changes seemed to occur irrespective of clinical outcomes, with the exception of the increase in CD68 mRNA levels which resulted associated with higher numbers of exacerbations post-BT.

ACTA2 gene (also known as α-SMA) encodes for smooth muscle actin. Its expression is characteristic of smooth muscle cells and myofibroblasts. There are some reports showing that one of the effects of BT treatment is the reduction of airway α-SMA [6, 13–15, 20, 21, 30, 31]. Herein data confirm this finding at mRNA level in bronchial biopsies from a wide cohort of asthmatic patients after BT treatment.

CD68 is mainly expressed by monocytes, macrophages (both M1 and M2 macrophages) and dendritic cells. The increase of CD68 positive cells in bronchial tissues from asthmatic patients compared with healthy controls is known [32]. Moreover, CD68 expression has been previously reported as increased after BT treatment at protein level in bronchial biopsies from 12 asthmatic patients by our group [6], herein confirmed at mRNA level by real-time PCR in a larger cohort of patients.

FAP expression is low in most adult tissues under physiological conditions but it is usually high in reactive stromal fibroblasts. FAP is thought to be involved in the control of fibroblast growth and epithelial-mesenchymal transitions during development, tissue repair, and carcinogenesis. In our knowledge, this is the first report that documented the up-regulation of FAP following BT treatment.

COL1A1 and COL1A2 mRNAs encode pro-alpha I and pro-alpha II chains of type I collagen. Chakir J. and collaborators have reported a reduction in type I collagen thickness below the basement membrane by immunohistochemistry in 9 bronchial biopsies collected from patients at third BT session in comparison to baseline [15]. This alteration persisted after 27 months of follow up [16]. In contrast, we found a huge increase of mRNAs for type I collagen in 89% (COL1A1) and 63% (COL1A2) of bronchial biopsies during BT treatment. Our results at mRNA levels agree with the observations of Pretolani and collaborators showing significantly larger and more intense area stained for collagen in bronchial biopsies at the last BT session versus baseline [20]. The explanation for the discordant results could be found in the fact that Chakir J et al [15]. and Salem J.H. et al [16]. focused the attention only on the collagen localized under the reticular layer of the basement membrane, while we analyzed the total content of collagen, as well as Pretolani et al did [20].

Asthma is characterized by airway hyperresponsiveness and it is known that there is a direct correlation with airway remodelling [33]. At cellular level, airway remodelling consists in smooth muscle hyperplasia characterized by a thicker α-SMA layer, epithelial damage and collagen deposition in the subepithelial basement membrane [21, 33–35]. The emerging idea is that the deposition of connective tissue around the airways can determine tissue stiffness, dampening the contractility of α-SMA layer [35, 36]. Based on the present data we speculate that BT treatment could reduce airway hyperresponsiveness decreasing α-SMA layer and increasing fibroblasts with collagen production.

PGP9.5 is a protein expressed in the innervation of mucosal glands, smooth muscle and mucosal blood vessels of the lungs [37]. Previous studies, which used PGP9.5 as a marker for the identification and count of nerve C fibers by immunohistochemistry, have reported a reduction of neural innervation after BT treatment [6, 20, 38]. Herein, changes in PGP9.5 mRNA levels during BT treatment were heterogeneous among patients. However, a higher reduction in PGP9.5 mRNA levels correlated with better patients reported outcome in terms of ΔAQLQ, suggesting that a reduction of innervation might be linked with BT efficacy.

To identify predictors of BT efficacy we associated normalized gene expression levels at baseline and T2 as well as changes in gene expression at T2 versus baseline with clinical parameters 12 months after BT. At baseline, lower mRNA levels of OCLN and COL1A2 resulted associated with parameters of better outcome: lower numbers of exacerbations and lower ACQ values post-BT. Moreover, it appeared valuable to quantify gene expression during BT, at the third session of BT. Indeed, lower mRNA levels of OCLN, CD68, CTGF and higher mRNA levels of SLPI at T2, lower changes in CD68 and CTGF mRNA levels at T2 versus T0 resulted associated with fewer exacerbations post-BT. We speculate that the quantification of OCLN, CD68, CTGF and SLPI mRNA could be useful to identify patients at higher risk of exacerbations following BT treatment.

OCLN is a key component of tight junctions, it is involved in the creation of a barrier in airway epithelial cells reducing the transport from apical to basolateral surface. Reduction in OCLN levels has been associated with impaired barrier function in asthma patients [39, 40]. OCLN is also expressed in tight junctions by myelinating Schwann cells protecting peripheral nerves and promoting neurons communication [41]. This is the first study that investigated OCLN in the context of BT treatment and further in-depth analyses will be necessary to elucidate the relationship between tight junction functionality and development of exacerbations.

CTGF gene expression has been linked to airway smooth muscle cells [42] and fibroblast to myofibroblast transition [43]. Blocking CTGF has been suggested as a rationale option in patients with asthma [44].

SLPI gene expression has emerged as the most promising marker at baseline to predict response to BT in the work by Ano S et al [24]. SLPI can inhibit leukocyte elastase, cathepsin G, trypsin and mast cell chymase. Based on Figure 3 shown in that manuscript [24], it seems even possible to define a cut-off in SLPI expression in bronchial specimens which could distinguish responders from non-responders to BT [24]. We were not able to confirm the usefulness of the quantification of SLPI gene expression at baseline to predict response to BT. However, the quantification of SLPI gene expression levels during BT resulted the best marker associated with the numbers of exacerbations post-BT, thus confirming a potential value for SLPI.

We did not detect any expression of IL-4, IL-13 and IL-17A mRNAs and a limited expression of IL-5 mRNA in bronchial samples. These cytokines, markers of Th2 and Th17 lymphocytes, have been reported as expressed in samples from asthmatic patients (bronchial biopsies, epithelial brushings, sputum, bronchoalveolar lavage) [2] but usually by means of different techniques: immunohistochemistry and in situ hybridization. One study has investigated the expression of IL-4, IL-5 and IL-13 by PCR: no expression has been detected unless tissue fragments were stimulated ex vivo [45]. Recently, Dr Pretolani and collaborators have reported a decrease in IL-13-positive cells per mm² by immunohistochemistry in bronchial biopsies after BT [31]. Differences in the results can be explained hypothesizing a scarce correlation between mRNA and protein expression and/or a different ability to detect these markers in tissue samples between PCR and immunohistochemistry.

The limit of the present study is the hypothesis-driven selection of the investigated genes. However, the results are in part consistent with those of previous studies, in part open up new possibilities. A high-throughput approach would allow to discover additional mechanisms of BT treatment. Analysis of the transcriptome has been recently performed on samples from really small cohorts of patients who underwent BT treatment: on cytology brushes from tracheal walls of 5 patients [46], and on bronchial biopsies from 8 patients [24]. Furthermore, Sun Q. et al. used a high-throughput approach on 9 primary airway epithelial cell lines obtained from patients before and after BT treatment and cultured in vitro [47]. In our knowledge this study analyzed the largest cohort of patients who underwent BT, coupling molecular and clinical data at baseline, at the last BT session and after 1 year of follow up. Moreover, the use of real-time PCR allowed us to analyze in each sample a higher number of markers compared with previous studies [6, 13–17, 20, 21, 30, 31].

In conclusion, our study revealed that BT was effective in 81% of patients irrespective of asthma endotypes and was associated with a reduction in ACTA2 and an increase in FAP, COL1A1, COL1A2 and CD68 gene expression. Quantification of COL1A2, OCLN, CD68, CTGF, SLPI at baseline and at the third BT session might be useful to identify patients at higher risk for exacerbations post-BT. These results are exploratory, further studies are needed to test whether the investigation of the expression of any of these genes could be useful to support clinical practice.

SRA

Severe refractory asthma

IgE

Immunoglobulin E

Bronchial Thermoplasty

AQLQ

Asthma quality of life questionnaire

ACQ, Asthma control questionnaire

TRPV1

transient receptor potential cation channel subfamily V member 1

TRPV2

transient receptor potential cation channel subfamily V member 2

UCHL1

Ubiquitin carboxy-terminal hydrolase L1 (also known as PGP9.5)

ACTA2

alpha smooth muscle actin

CTGF

Connective tissue growth factor

FAP

Fibroblast activation protein alpha

IL-4

Interleukin-4

IL-5

Interleukin-5

IL-6

Interleukin-6

IL-8

IL-13:Interleukin-13

IL-17

Interleukin-17

COL1A1

Alpha-1 type I collagen

COL1A2

Alpha-2 type I collagen

POSTN

Periostin

LGAL3

Galectin-3

OCLN

Occludin

CDH1

Cadherin 1

GAPDH

Glyceraldehyde 3-phosphate dehydrogenase

SLPI, Secretory leukocyte protease inhibitor

OCSs

Oral corticosteroids.

Acknowledgements

Not applicable

Authors' contributions

Nicola Facciolongo had full access to all of data in the study and takes responsibility for the integrity of data and the accuracy of data analysis, including any adverse effects. Nicola Facciolongo and Stefania Croci conceived the study. Martina Bonacini and Stefania Croci performed gene expression profiling. Martina Bonacini, Carla Galeone, Patrizia Ruggiero, Francesco Menzella, Giulia Ghidoni, Roberto Piro, Chiara Scelfo, Chiara Catellani, Alessandro Zerbini, Stefania Croci contributed substantially to the study design, data analysis, interpretation and writing of the manuscript.

Funding

There is no funding to report.

Availability of data and materials

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Declarations

Ethics approval and consent to participate

Written informed consent was obtained from all subjects who participated in the study. The study protocol was approved by the Local Ethics Committee (Azienda USL-IRCCS of Reggio Emilia - Comitato Etico Area Vasta Nord of Emilia Romagna; protocol number: 2019/0014076).

Consent for publication

All authors are aware of and agree to the content of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Authors details

¹Pneumology Unit, Azienda Unità Sanitaria Locale-IRCCS di Reggio Emilia, Reggio Emilia, Italy.

²Unit of Clinical Immunology, Allergy and Advanced Biotechnologies, Azienda Unità Sanitaria Locale-IRCCS di Reggio Emilia, Reggio Emilia, Italy.

Corren J. Asthma phenotypes and endotypes: an evolving paradigm for classification. Discov Med. 2013;15(83):243–9.
Galeone C, Scelfo C, Bertolini F, Caminati M, Ruggiero P, Facciolongo N, et al. Precision Medicine in Targeted Therapies for Severe Asthma: Is There Any Place for "Omics" Technology? Biomed Res Int. 2018;2018:4617565.
Thomson NC. Novel approaches to the management of noneosinophilic asthma. Ther Adv Respir Dis. 2016;10(3):211–34.
https://ginasthma.org/wp-content/uploads/2021/05/GINA-Main-Report-2021-V2-WMS.pdf Accessed date: 28 december 2021.
Niven RM, Simmonds MR, Cangelosi MJ, Tilden DP, Cottrell S, Shargill NS. Indirect comparison of bronchial thermoplasty versus omalizumab for uncontrolled severe asthma. J Asthma. 2018;55(4):443–51.
Facciolongo N, Di Stefano A, Pietrini V, Galeone C, Bellanova F, Menzella F, et al. Nerve ablation after bronchial thermoplasty and sustained improvement in severe asthma. BMC Pulm Med. 2018;18(1):29.
Cox G, Thomson NC, Rubin AS, Niven RM, Corris PA, Siersted HC, et al. Asthma control during the year after bronchial thermoplasty. N Engl J Med. 2007;356(13):1327–37.
Pavord ID, Cox G, Thomson NC, Rubin AS, Corris PA, Niven RM, et al RISA Trial Study Group. Safety and efficacy of bronchial thermoplasty in symptomatic, severe asthma. Am J Respir Crit Care Med. 2007;176(12):1185–91.
Castro M, Rubin AS, Laviolette M, Fiterman J, De Andrade Lima M, Shah PL, et al. Effectiveness and safety of bronchial thermoplasty in the treatment of severe asthma: a multicenter, randomized, double-blind, sham-controlled clinical trial. Am J Respir Crit Care Med. 2010;181(2):116–24.
Chupp G, Laviolette M, Cohn L, McEvoy C, Bansal S, Shifren A, et al. Other members of the PAS2 Study Group. Long-term outcomes of bronchial thermoplasty in subjects with severe asthma: a comparison of 3-year follow-up results from two prospective multicentre studies. Eur Respir J. 2017;50(2):1700017.
Russell RJ, Brightling CE. Bronchial thermoplasty: what we know, what we don't know, and what we need to know. Eur Respir J. 2021 Dec;31(1):2102018. 59(.
Thomson NC. Bronchial thermoplasty as a treatment for severe asthma: controversies, progress and uncertainties. Exp Rev Respir Med. 2018;12(4):269–82.
Pretolani M, Dombret MC, Thabut G, Knap D, Hamidi F, Debray MP, et al. Reduction of airway smooth muscle mass by bronchial thermoplasty in patients with severe asthma. Am J Respir Crit Care Med. 2014;190(12):1452–4.
Denner DR, Doeing DC, Hogarth DK, Dugan K, Naureckas ET, White SR. Airway Inflammation after Bronchial Thermoplasty for Severe Asthma. Ann Am Thorac Soc. 2015;12(9):1302–9.
Chakir J, Haj-Salem I, Gras D, Joubert P, Beaudoin ÈL, Biardel S, et al. Effects of Bronchial Thermoplasty on Airway Smooth Muscle and Collagen Deposition in Asthma. Ann Am Thorac Soc. 2015;12(11):1612–8.
Salem IH, Boulet LP, Biardel S, ampron N, Martel S, Laviolette M, et al. Long-Term Effects of Bronchial Thermoplasty on Airway Smooth Muscle and Reticular Basement Membrane Thickness in Severe Asthma. Ann Am Thorac Soc. 2016;13(8):1426–8.
Haj Salem I, Gras D, Joubert P, Boulet LP, Lampron N, Martel S, Godbout K, Chanez P, Laviolette M, Chakir J. Persistent Reduction of Mucin Production after Bronchial Thermoplasty in Severe Asthma. Am J Respir Crit Care Med. 2019;199(4):536–8.
Hall CS, Castro M. Predicting response to bronchial thermoplasty in patients with severe uncontrolled asthma: an elusive goal. Respirology. 2019;24(1):11–12. http://dx.doi.org/10.1111/resp.13424 PMID: 30324712.
Langton D, Sha J, Ing A, Fielding D, Thien F, Plummer V. Bronchial thermoplasty: activations predict response. Respir Res. 2017;18(1):134. http://dx.doi.org/10.1186/s12931-017-0617-7 PMID: 28676053.
Pretolani M, Bergqvist A, Thabut G, Dombret MC, Knapp D, Hamidi F, et al. Effectiveness of bronchial thermoplasty in patients with severe refractory asthma: Clinical and histopathologic correlations. J Allergy Clin Immunol. 2017;139(4):1176–85.
d'Hooghe JNS, Ten Hacken NHT, Weersink EJM, Sterk PJ, Annema JT, Bonta PI. Emerging understanding of the mechanism of action of Bronchial Thermoplasty in asthma. Pharmacol Ther. 2018;181:101–7.
Geppetti P, Materazzi S, Nicoletti P. The transient receptor potential vanilloid 1: Role in airway inflammation and disease. Eur J Pharm. 2006;533:207–14.
Yanlin Jia L-Y, Lee. Role of TRPV receptors in respiratory diseases. Biochem Biophys Acta. 2007;1772:915–27.
Ano S, Kikuchi N, Matsuyama M, Nakajima M, Kondo Y, Masuda M, Osawa H, Ishii Y, Hizawa N. Transcriptome genetic differences between responders and non-responders before bronchial thermoplasty. J Asthma. 2021 Jul;8:1–11.
Holguin F, Cardet JC, Chung KF, Diver S, Ferreira DS, Fitzpatrick A, et al. Management of severe asthma: a European Respiratory Society/American Thoracic Society guideline. Eur Respir J. 2020;55(1):1900588.
Langton D, Wang W, Sha J, Ing A, Fielding D, Hersch N, et al. Predicting the Response to Bronchial Thermoplasty. J Allergy Clin Immunol Pract. 2020;8(4):1253–60.e2.
Bel EH. Bronchial thermoplasty: has the promise been met? Am J Respir Crit Care Med. 2010;181(2):101–2.
Bel EH, Zwinderman AH. Outcome reporting in asthma research. Am J Respir Crit Care Med. 2011;183(1):132.
Iyer VN, Lim KG. Bronchial thermoplasty: reappraising the evidence (or lack thereof). Chest. 2014;146(1):17–21.
Ichikawa T, Panariti A, Audusseau S, Mogas AK, Olivenstein R, Chakir J, et al. Effect of bronchial thermoplasty on structural changes and inflammatory mediators in the airways of subjects with severe asthma. Respir Med. 2019;150:165–72.
Ladjemi MZ, Di Candia L, Heddebaut N, Techoueyres C, Airaud E, Soussan D, et al. Clinical and histopathologic predictors of therapeutic response to bronchial thermoplasty in severe refractory asthma. J Allergy Clin Immunol. 2021:S0091-6749(21)00005-1..
Eguíluz-Gracia I, Malmstrom K, Dheyauldeen SA, Lohi J, Sajantila A, Aaløkken R, et al. Monocytes accumulate in the airways of children with fatal asthma. Clin Exp Allergy. 2018;48(12):1631–9.
Boulet LP, Turmel J, Côté A. Asthma and exercise-induced respiratory symptoms in the athlete: new insights. Curr Opin Pulm Med. 2017;23(1):71–7.
Niimi A, Matsumoto H, Takemura M, Ueda T, Chin K, Mishima M. Relationship of airway wall thickness to airway sensitivity and airway reactivity in asthma. Am J Respir Crit Care Med. 2003;168(8):983–8.
Paré PD. Airway hyperresponsiveness in asthma: geometry is not everything! Am J Respir Crit Care Med. 2003;168(8):913–4.
Draijer C, Robbe P, Boorsma CE, Hylkema MN, Melgert BN. Dual role of YM1+ M2 macrophages in allergic lung inflammation. Sci Rep. 2018;8(1):5105.
West PW, Canning BJ, Merlo-Pich E, Woodcock AA, Smith JA. Morphologic Characterization of Nerves in Whole-Mount Airway Biopsies. Am J Respir Crit Care Med. 2015;192(1):30–9.
Luo YL, Cheng YQ, Zhou ZQ, Fan MY, Chen DF, Chen Y, Chen XB, Zhong CH, Tang CL, Li SY, Su ZQ. A clinical and canine experimental study in small-airway response to bronchial thermoplasty: Role of the neuronal effect. Allergol Int. 2022;71(1):66–72.
Wawrzyniak P, Wawrzyniak M, Wanke K, Sokolowska M, Bendelja K, Rückert B, et al. Regulation of bronchial epithelial barrier integrity by type 2 cytokines and histone deacetylases in asthmatic patients. J Allergy Clin Immunol. 2017;139(1):93–103.
Steelant B, Farré R, Wawrzyniak P, Belmans J, Dekimpe E, Vanheel H, et al. Impaired barrier function in patients with house dust mite-induced allergic rhinitis is accompanied by decreased occludin and zonula occludens-1 expression. J Allergy Clin Immunol. 2016;137(4):1043–53.e5.
Reinhold AK, Rittner HL. Barrier function in the peripheral and central nervous system-a review. Pflugers Arch. 2017;469(1):123–34.
Wang J, Faiz A, Ge Q, Vermeulen CJ, Van der Velden J, Snibson KJ, et al. Unique mechanisms of connective tissue growth factor regulation in airway smooth muscle in asthma: Relationship with airway remodelling. J Cell Mol Med. 2018;22(5):2826–37.
Wójcik-Pszczoła K, Jakieła B, Plutecka H, Koczurkiewicz P, Madeja Z, Michalik M, et al. Connective tissue growth factor regulates transition of primary bronchial fibroblasts to myofibroblasts in asthmatic subjects. Cytokine. 2018;102:187–90.
Gao W, Cai L, Xu X, Fan J, Xue X, Yan X, et al. Anti-CTGF single-chain variable fragment dimers inhibit human airway smooth muscle (ASM) cell proliferation by down-regulating p-Akt and p-mTOR levels. PLoS One. 2014;9(12):e113980.
Jaffar Z, Roberts K, Pandit A, Linsley P, Djukanovic R, Holgate S. B7 costimulation is required for IL-5 and IL-13 secretion by bronchial biopsy tissue of atopic asthmatic subjects in response to allergen stimulation. Am J Respir Cell Mol Biol. 1999;20(1):153–62.
Liao SY, Linderholm AL, Yoneda KY, Kenyon NJ, Harper RW. Airway transcriptomic profiling after bronchial thermoplasty. ERJ Open Res. 2019;5(1).
Sun Q, Fang L, Roth M, Tang X, Papakonstantinou E, Zhai W, et al. Bronchial thermoplasty decreases airway remodelling by blocking epithelium-derived heat shock protein-60 secretion and protein arginine methyltransferase-1 in fibroblasts. Eur Respir J. 2019;54(6).

SupplementaryTable.docx
SupplementaryFigureS1.tif
Supplementary figure S1 Gene expression levels at baseline according to OCS dose.Expressions of the investigated genes in bronchial biopsies at T0 from patients classified in three groups based on the OCS dose: ≤ 5 mg (n = 8); 5 – 25 mg (n = 7); ≥ 25 mg (n = 12). Gene expressions were calculated by the 2^{− ΔCt} method using the GAPDH as housekeeper gene. Horizontal lines show the median ± interquartile range (IQR). Data were analysed by Kruskal-Wallis test. No statistically significant differences were found.
SupplementaryFigureS2.tif
Supplementary Figure S2 Correlations between fold changes in gene expression and ΔAQLQ during BT treatmentDot plot visualization of correlations between fold changes in gene expression and differences in AQLQ scores between T2 and T0. Fold changes in gene expression were determined by real-time PCR relative to gene expression at T0 through the 2^-ΔΔCt method. Data were analyzed by Spearman's correlation test (n = 27). r = correlation coefficient; 95% CI = 95% confidence interval.
SupplementaryFigureS3.tif
Supplementary Figure S3 Correlations between fold changes in gene expression and ΔACQ during BT treatmentDot plot visualization of correlations between fold changes in gene expression and differences in ACQ scores between T2 and T0. Fold changes in gene expression were determined by real-time PCR relative to gene expression at T0 through the 2^-ΔΔCt method. Data were analyzed by Spearman's correlation test (n = 27). r = correlation coefficient; 95% CI = 95% confidence interval.
SupplementaryFigureS4.tif
Supplementary figure S4 Gene expressions at baseline in patients stratified on exacerbations post-BTExpression of the investigated genes in bronchial biopsies at T0 from patients classified according to the number of exacerbations during the 12 months after BT (n = 27). Gene expressions were calculated by the 2^{− ΔCt} method using the GAPDH as housekeeper gene. Horizontal lines show the median ± interquartile range (IQR). Data were analyzed by Mann-Whitney U test.
SupplementaryFigureS5.tif
Supplementary figure S5 Gene expressions at T2 in patients stratified on exacerbations post-BTExpression of the investigated genes in bronchial biopsies at T2 from patients classified according to the number of exacerbations during the 12 months after BT (n = 27). Gene expressions were calculated by the 2^{− ΔCt} method using the GAPDH as housekeeper gene. Horizontal lines show the median ± interquartile range (IQR). Data were analyzed by Mann-Whitney U test.
SupplementaryFigureS6.tif
Supplementary Figure S6 Gene expressions at baseline in patients stratified on the ACQ questionnaire scores post-BTExpression of the investigated genes in bronchial biopsies at T0 from patients classified on the ACQ questionnaire scores 12 months after BT (n = 27). ACQ < 1.5 controlled asthma; ACQ ≥ 1.5 uncontrolled asthma. Gene expressions were calculated by the 2^{− ΔCt} method using the GAPDH as housekeeper gene. Horizontal lines show the median ± interquartile range (IQR). Data were analyzed by Mann-Whitney U test.
SupplementaryFigureS7.tif
Supplementary Figure S7 Fold changes in gene expressions in patients stratified on exacerbations post-BTFold changes in gene expression (T2 versus T0) in bronchial biopsies grouping patients according to the numbers of exacerbations experienced during the 12 months of follow up post-BT (n = 27). Horizontal lines show the median ± interquartile range (IQR). Data were analyzed by Mann-Whitney U test.

Download PDF

Reviews received at journal
08 Mar, 2022
Reviewers invited by journal
08 Mar, 2022
Editor assigned by journal
02 Feb, 2022
Submission checks completed at journal
01 Feb, 2022
Editor invited by journal
01 Feb, 2022
First submitted to journal
30 Jan, 2022

You are reading this latest preprint version

Bronchial Thermoplasty in severe asthma: a real-world study on efficacy and gene profiling

Status:

Version 1

Abstract

Figures

Background

Methods

Study population

Measures to determine BT response

BT procedure

Airway sample collection

RNA extraction and cDNA preparation

Real-time PCR

Statistical analysis

Results

Demographic results at baseline

BT clinical outcomes

Characterization of bronchial gene expression at baseline

Effects of BT treatment on gene expression

Association between gene expression and clinical parameters 12 months post-BT

Discussion

Conclusions

Abbreviations

Declarations

References

Supplementary Files

Status:

Version 1