Abstract
Bronchial thermoplasty (BT) is to date the only therapy that provides a lasting reduction in airway wall remodelling. However, the mechanism of action of BT is not well understood. This study aimed to characterise the changes of remodelling regulating signalling pathways by BT in asthma.
Bronchoalveolar lavage fluid (BALF) was obtained from eight patients with severe asthma before and after BT. Primary bronchial epithelial cells were isolated from 23 patients before (n=66) and after (n=62) BT. Epithelial cell culture supernatant (Epi.S) was collected and applied to primary fibroblasts.
Epithelial cells obtained from asthma patients after BT proliferated significantly faster compared with epithelial cells obtained before BT. In airway fibroblasts, BALF or Epi.S obtained before BT increased CCAAT enhancer-binding protein-β (C/EBPβ) expression, thereby downregulating microRNA-19a. This upregulated extracellular signal-regulated kinase-1/2 (ERK1/2) expression, protein arginine methyltransferase-1 (PRMT1) expression, cell proliferation and mitochondrial mass. BALF or Epi.S obtained after BT reduced the expression of C/EBPβ, ERK1/2, peroxisome proliferator-activated receptor-γ coactivator-1α (PGC1α), PRMT1 and mitochondrial mass in airway fibroblasts. Proteome and transcriptome analyses indicated that epithelial cell-derived heat shock protein-60 (HSP60) is the main mediator of BT effects on fibroblasts. Further analysis suggested that HSP60 regulated PRMT1 expression, which was responsible for the increased mitochondrial mass and α-smooth muscle actin expression by asthmatic fibroblasts. These effects were ablated after BT. These results imply that BT reduces fibroblast remodelling through modifying the function of epithelial cells, especially by reducing HSP60 secretion and subsequent signalling pathways that regulate PRMT1 expression.
We therefore hypothesise that BT decreases airway remodelling by blocking epithelium-derived HSP60 secretion and PRMT1 in fibroblasts.
Abstract
Bronchial thermoplasty reduces asthma-associated airway wall remodelling by epithelial cell-controlled downregulation of PRMT1 expression in fibroblasts http://bit.ly/2TLhCkC
Introduction
Asthma is characterised by irreversible structural changes of the airway wall [1, 2], resulting in hyperresponsiveness and stiffness of the tissue [3, 4]. Remodelling is not reduced by asthma therapies, except bronchial thermoplasty (BT) [5, 6]. In clinical trials, BT improved quality of life, and reduced exacerbation rates and hospitalisation, which correlated with reduced airway remodelling [6, 7].
The function of the airways is regulated by the interaction of epithelial and subepithelial mesenchymal cells. In asthma, this communication is disrupted, leading to irreversible tissue structural changes [8, 9]. Airway wall remodelling is resistant to all available asthma drugs [10, 11], whereas BT was found to modify the function of subepithelial mesenchymal cells and the composition of the extracellular matrix, and to reduce inflammation [5, 6, 12, 13]. BT was also found to reduce the mass of airway smooth muscle [5, 6, 14], and decrease neuroendocrine-activated epithelial cells and inflammation [11–13]. However, it is unclear if all these effects resulted from “burning off smooth muscle”.
Thus, BT might achieve its beneficial effects by modulating the methylation of DNA, histone-3 or noncoding RNA, which have been linked to asthma-associated remodelling [15, 16]. A major enzyme controlling methylation is protein arginine methyltransferase-1 (PRMT1), which was constitutively upregulated in airway smooth muscle cells of asthma patients [17]. PRMT1 contributes to the pathogenesis of lung cancer, pulmonary fibrosis, pulmonary hypertension, chronic obstructive pulmonary disease (COPD) and asthma [18]. In an animal model of allergic asthma, PRMT1 was cell type specifically regulated in airway epithelial and mesenchymal cells during acute and chronic inflammation [19–22]. Platelet-derived growth factor-BB stimulated PRMT1 expression and remodelling by downregulating microRNA-19a (miR-19a), thereby upregulating extracellular signal-regulated kinase-1/2 (ERK1/2) mitogen-activated protein kinase (MAPK) expression and activity [17, 23].
This study investigated if BT affects the regulatory mechanism of PRMT1 on mitochondria activity and subsequent remodelling in airway fibroblasts.
Methods
Patient cohorts
In total, 66 endobronchial tissue biopsies were obtained from 51 patients with severe asthma (according to Global Initiative for Asthma guidelines) by the lung clinic (University Hospital Basel, Basel, Switzerland) before BT and 62 biopsies were obtained from 23 patients at 4, 8 and 12 weeks after BT. BT was performed according to a standardised protocol. Table 1 presents the clinical characteristics of the patients.
Clinical characteristics of the patient cohorts before and after bronchial thermoplasty (BT)
All patients provided written informed consent. The study was approved by the local Institutional Ethical Committee of Northwest and Central Switzerland (EKNZ 2016-01057).
Study design
Endobronchial biopsies and bronchoalveolar lavage fluid (BALF) were collected before and after BT. Tissue biopsies were used to isolate primary epithelial cells and fibroblasts (see the Cell isolation and characterisation section). BALF was aliquoted and used for proteomic analysis or was stored (−20°C) for experiments with fibroblasts. A graphic summary of the study design is provided in figure 1.
Study design. BT: bronchial thermoplasty; BALF: bronchoalveolar lavage fluid; ECM: extracellular matrix; siRNA: small interfering RNA; NGS: next-generation sequencing; KEGG: Kyoto Encyclopedia of Genes and Genomes. Endobronchial biopsies and BALF were obtained from 23 patients with severe asthma before and after undergoing BT. Primary epithelial cells and fibroblasts were isolated by selective medium from each biopsy. Experimental BT (10 s, 65°C submerged in a water bath) was performed with epithelial cells, and cell culture medium was collected after 24 h and used to treat fibroblasts. In parallel, fibroblasts were exposed to BALF (diluted 1:1 in fibroblast medium) over different time intervals.
BALF and bronchoscopy
BALF was collected and bronchoscopy was performed according to standard procedures before and after BT.
BALF proteomic analysis and Kyoto Encyclopedia of Genes and Genomes pathway analysis
Proteomic analyses were performed in BALF samples of seven asthma patients before and after BT and compared in pairs. Mass spectrometric analysis was performed in 200 µL BALF by the SWATH method on a Q-Exactive HF platform (Thermo Fisher Scientific, Waltham, MA, USA) and analysed by Spectronaut software (https://biognosys.com) against the Proteomics Core Facility internal HEK cell line spectral library (Biocentre, University of Basel, Basel, Switzerland). The retention time was calibrated by internal indexed retention time standard peptides (iRT Kit; Biognosys). All data are provided as supplementary material.
Statistical analysis of BALF proteomics was performed using OmicShare tools (www.omicshare.com/tools) after removing results with low abundance values. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was conducted to identify BT-modified signalling pathways by comparing protein profiles before and after BT.
Cell isolation and characterisation
Primary bronchial epithelial cells were isolated from biopsies by cell type selective medium CnT-PR-A (CellnTec, Bern, Switzerland), and were characterised by positive staining for E-cadherin and cytokeratin-13 and negative staining for fibronectin (supplementary figure S1). The culture epithelial cell supernatant (Epi.S) was collected from confluent epithelial cells.
Subepithelial bronchial fibroblasts were generated from nonasthma (n=8) and asthma (n=10) patients in selective medium (CnT-PR-F; CellnTec), and were characterised by positive staining for fibronectin and negative staining for E-cadherin and fibrilar α-smooth muscle actin (α-SMA) (supplementary figure S1) [20]. All experiments were performed between passage 2 and 5.
Next-generation sequencing
Next-generation sequencing (NGS) analysis was performed in primary airway epithelial cell lines (n=9) obtained before and after experimental BT. Confluent cell layers (25 cm flasks) were exposed for 10 s to either 37°C (control) or 65°C (experimental BT) in a water bath. The flasks were then placed in a humidified incubator for 24 h at 37°C, 5% CO2. Total RNA was isolated and assessed for differential mRNA expression of the entire human genome by NGS. RNA libraries were prepared using a TruSeq Stranded mRNA Sample Kit (Illumina, San Diego, CA, USA). mRNA was enriched from total RNA (500 ng) by an oligo-dT bead system and fragmented enzymatically. First- and second-strand syntheses were performed, and double-stranded cDNA was purified (AMPure XP; Thermo Fisher Scientific). The cDNA was 3′-adenylated and Illumina sequencing adaptors ligated onto the fragments. mRNA libraries were pre-amplified by PCR and purified (AMPure XP), and the size distribution was validated and quality inspected on a BioAnalyzer 4200 TapeStation (Agilent Technologies, Santa Clara, CA, USA). High-quality libraries were pooled at equimolar concentrations (BioAnalyzer Smear Analysis; Agilent Technologies) and quantified by quantitative PCR. Optimised library pools were used to generate clusters on the surface of a flow cell before sequencing on a NextSeq 500 instrument (75 cycles; Illumina). The NGS analysis report is provided as supplementary material.
Fibroblast stimulation and small interfering RNA treatment
Fibroblasts grown in six-well plates were serum deprived (0.1% serum) overnight before being treated with Epi.S or BALF (mixed 1:1 with basic medium) for 24 and 48 h.
PRMT1 activity was inhibited by 2 h pre-incubation with 10 µM of the protein arginine methyltransferase inhibitor AMI-1 (Santa Cruz Biotechnology, Dallas, TX, USA).
Quiescent fibroblasts were treated with heat shock protein-60 (HSP60) (1 μg·mL−1; Sino Biological, Beijing, China) or transfected by HiPerFect (TaKaRa, Kusatsu, Japan) for 48 h with CCAAT enhancer-binding protein-β (C/EBPβ) small interfering RNAs (siRNAs) or control sequences (all 50 nM) [19–23].
miR-19a real-time quantitative PCR
miR-19a was detected in 500 ng of total RNA as described previously [17, 23] using the primer sequence 5′-TGT GCA AAT CTA TGC AAA ACT GA-3′. miR-19a expression was normalised to U6 small nuclear RNA and analysed by the 2−ΔΔCt (relative quantification) method [17].
Western blotting
Fibroblasts were lysed in RIPA buffer and protein concentration was quantified by bicinchoninic acid (Thermo Fisher Scientific). Equal amounts of denatured proteins (20 µg) were size fractionated by electrophoresis (8–16% SDS–PAGE; Thermo Fisher Scientific) and transferred onto polyvinylidene fluoride membranes. Proteins were detected by antibodies (PRMT1, peroxisome proliferator-activated receptor-γ coactivator-1α (PGC1α) or HSP60 (Abcam, Cambridge, UK), or phosphorylated ERK1/2 MAPK, total ERK1/2 MAPK, C/EBPβ or glyceraldehyde 3-phosphate dehydrogenase (Cell Signaling Technology, Danvers, MA, USA)). Protein bands were visualised after binding of species-specific secondary horseradish peroxidase-conjugated antibodies (Abcam) by chemiluminescence (Thermo Fisher Scientific) [19–23].
Immunofluorescence
Fibroblasts, seeded on 10 mm sterile glass slides and serum deprived overnight, were treated with Epi.S or BALF for 0–48 h, fixed in 4% paraformaldehyde (in PBS, 2×5 min) and immunofluorescence staining was performed as described previously [17, 23]. Nuclei were stained by 4′,6-diamidino-2-phenylindole (DAPI).
MitoTracker staining
MitoTracker (Thermo Fisher Scientific) staining was performed in quiescent fibroblasts (80% confluence) as described previously [23]. MitoTracker imaging was captured by an Olympus BX61 fluorescence microscope (Olympus Optical, Tokyo, Japan). Nuclei were stained by DAPI. Image analysis was performed by ImageJ version 2.0 (https://imagej.nih.gov/ij).
Fibroblast activity and proliferation
Fibroblast activity was determined after incubation (4 h, 37°C) with 500 µg·mL−1 MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide), before being lysed with dimethyl sulfoxide and absorbance measured at 560 nm (microplate reader; BioTek, Winooski, VT, USA).
Manual cell counts (Neubauer chamber slide) were performed after 48 h and in the presence or absence of chemical inhibitors for signalling proteins or PRMT1 (AMI-1) [17, 23].
Epithelial cell numbers were determined by calculating the number of cells within a circle of 0.5 cm diameter every second day over 10 days.
ELISA of collagen type I and fibronectin deposition
Deposition of collagen type I and fibronectin was measured by ELISA [23]. Confluent fibroblasts in 96-well plates were starved for 24 h and stimulated with Epi.S. Collagen type I and fibronectin were detected after 24 h (both antibodies 1:200; Santa Cruz Biotechnology).
Immunohistochemistry
Bronchial biopsies were embedded in paraffin, cut into 5 µm slices and fixed onto glass slides. After de-paraffinisation, rehydration and antigen unmasking by standard procedures, slides were blocked in 5% bovine serum albumin (BSA), stained for PRMT1 and phosphorylated ERK1/2 MAPK (4°C, overnight), followed by antibody detection [19–23]. Images were examined by an Olympus BX61 white balance microscope. The percentage of positive cells was determined by double-blind visual quantification.
Confocal microscopy
Fibroblasts seeded onto coverslips were incubated according to experiments, followed by fixation and cell permeabilisation [17, 23]. Coverslips were blocked by 5% BSA in PBS, and then incubated with anti-cytochrome c antibody (BD Biosciences, San Jose, CA, USA) and Alexa Fluor 488-labelled secondary antibody (Thermo Fisher Scientific). F-actin was stained with TRIC-phalloidin and nuclei were stained with DAPI. Images were acquired by a Nikon Confocal A1 microscope (×40 1.3 NA FI oil objective of 0.225 μm Z-stack step; Nikon, Tokyo, Japan). Images are presented as Z-stack projections by imaging software FIJI (https://imagej.net/Fiji). Single-cell cytochrome c expression was calculated by FIJI with the “3D objects counter” plugin.
Statistical analysis
The null hypothesis was that no treatment had any effect on signalling or PRMT1 expression. All data are expressed as mean±sem. The statistical analysis was performed by ANOVA followed by the t-test for paired sample data. The effect of inhibitors on PRMT1 and remodelling were analysed by one-way ANOVA. p-values <0.05 were regarded as significant.
Results
BT altered the function of bronchial epithelial cells
The generation of primary bronchial epithelial cells from tissues obtained after BT was significantly faster compared with cells grown from tissue obtained before BT. The increase of cells was determined by direct cell counts over 10 days (figure 2a and b).
Bronchial thermoplasty (BT) altered the function of bronchial epithelial cells. BBT: before BT; ABT: after BT. a) Proliferation of epithelial cells isolated ABT was increased compared with epithelial cells isolated from the same patients BBT (n=18). b) The motility of epithelial cells increased ABT as monitored by counting epithelial cell colonies at a distance >1 cm from the biopsy (n=12). Data are presented as mean±sem. *: p<0.05; ***: p<0.001, comparing BBT versus ABT. c) Representative microscopic image of epithelial cell growth out from a tissue biopsy obtained BBT. d) Single epithelial cell at a distance >1 cm from tissue on the first day after starting the cell cultures of biopsies obtained ABT. e) Same single epithelial cell as shown in (d) 2 days later, forming a colony.
In addition to growing in direct contact to the biopsy tissue, epithelial cells growing at a distance >1 cm were defined as a “colony”. Before BT, 92% of all epithelial cells grew directly from the biopsy and expanded by forming a continuous cell monolayer. Thus, proliferation occurred mainly on the front of the epithelial cell layer (figure 2c). After BT, epithelial cells were significantly more often observed to grow at a distance from the biopsy, as soon as 1 day after placing the tissue in a flask (figure 2d). The position of such single epithelial cells was marked on the culture dish and cell numbers were counted daily. Within 2 days, single cells multiplied locally forming a “colony” by definition (figure 2e). Within 1 week, the colonies formed a homogenous monolayer of epithelial cells.
BT reduced ERK1/2 MAPK activation, PRMT1 expression and mitochondrial mass in lung tissue and primary fibroblasts
In tissue sections of asthma patients obtained before BT, ERK1/2 MAPK was phosphorylated in most epithelial cells and in the subepithelial mesenchymal cells. The ERK1/2 MAPK phosphorylation was significantly reduced in tissues obtained after BT (figure 3a). PRMT1 was mainly located in the nuclei (figure 3a; dark brown/black) and only a few cells showed cytosolic expression (figure 3a; brown), and both were reduced after BT, even 2 years after BT (supplementary figure S2). In isolated fibroblasts, BALF obtained before BT significantly decreased the expression of miR-19a, while BALF obtained after BT did not have this effect (figure 3b). As shown by representative Western blots and image analysis, BALF obtained before BT significantly increased ERK1/2 MAPK phosphorylation, and C/EBPβ and PRMT1 expression, while BALF obtained after BT did not (figure 3c).
Bronchial thermoplasty (BT) reduced extracellular signal-regulated kinase-1/2 (ERK1/2) activation, protein arginine methyltransferase-1 (PRMT1) expression and mitochondrial mass in lung tissue and primary fibroblasts. BBT: before BT; ABT: after BT; miR-19a: microRNA-19a; CM: control medium; C/EBPβ: CCAAT enhancer-binding protein-β; p-: phosphorylated; t-: total; GAPDH: glyceraldehyde 3-phosphate dehydrogenase; BALF: bronchoalveolar lavage fluid. a) Representative immunohistochemistry of p-ERK1/2 mitogen-activated protein kinase (MAPK) and PRMT1 in airway tissue sections obtained BBT and ABT. Black arrows: epithelium; yellow arrows: mesenchymal cells. Bar charts summarise image analyses results. Data are presented as median and 95% confidence interval of the positive percentile (n=10). Scale bars: 10 µm. *: p<0.05; **: p<0.01. b) miR-19a expression in control fibroblasts stimulated with BALF from asthma patients obtained BBT or ABT. Data are presented as mean±sem (n=9). c) Representative Western blot of the stimulatory effect of BALF from asthma patients BBT and ABT on p-ERK1/2 MAPK, t-ERK1/2 MAPK and PRMT1 in fibroblasts. Bar charts summarise image analyses results. The horizontal lines indicate the level for CM. Data are presented as mean±sem of five asthmatic and five nonasthma BALF samples, each tested in three fibroblast lines. *: p<0.05; **: p<0.01; ***: p<0.001, comparing healthy versus other fibroblasts; #: p<0.05; ##: p<0.01, comparing BBT versus ABT. d) Representative confocal microscopy images of the BALF effect on cytochrome c and F-actin (phalloidin). Nuclei were stained by 4′,6-diamidino-2-phenylindole (DAPI). Data are presented as mean±sem (n=11) for cytochrome c staining. Scale bars: 10 µm. *: p<0.05; ***: p<0.001.
Significantly increased mitochondrial mass in fibroblasts as monitored by cytochrome c expression was observed in BALF obtained before BT (figure 3d; middle row) compared with BALF obtained after BT (figure 3d; bottom row). No difference was found comparing the mitochondrial mass in nonasthma patients (n=5) with asthma patients (n=6). Thus, the data were combined for image analysis, which is presented in figure 3d as a bar chart.
Similar to BALF obtained before BT, Epi.S obtained from epithelial cells before BT stimulated PRMT1 expression by fibroblasts significantly when compared with Epi.S obtained from epithelial cells after BT (figure 4). Epi.S obtained after BT had a reduced stimulating effect on PRMT1 expression in fibroblasts. The effect on PRMT1 expression was time dependent and maximal after 48 h, but did not show a disease-specific effect (figure 4).
The effect of epithelial cell culture supernatant (Epi.S) obtained before and after bronchial thermoplasty (BBT and ABT, respectively) on protein arginine methyltransferase-1 (PRMT1) expression in fibroblasts: control fibroblast medium (CM), nonasthma epithelial cell culture medium (EM) or Epi.S obtained BBT and ABT was applied to primary fibroblasts from a) healthy or b) asthmatic donors. GAPDH: glyceraldehyde 3-phosphate dehydrogenase. Representative Western blots and quantitative image analysis of PRMT1 expression in fibroblasts of nonasthma and asthma patients at 24 and 48 h. Data are presented as mean±sem (n=4 in each group). **: p<0.01; ***: p<0.001, comparing CM versus others; ##: p<0.01, comparing BBT versus ABT.
BT reduced the stimulatory effect of Epi.S on fibroblast proliferation and mitochondrial mass
Epi.S of epithelial cells obtained before BT induced ERK1/2 MAPK phosphorylation, and PRMT1 and PGC1α expression. Western blot data was summarised as a bar chart (figure 5a) and confirmed by immunochemistry (figure 5b and c). The expression levels of PRMT1 and PGC1α correlated with MitoTracker signals (figure 5d).
Bronchial thermoplasty (BT) reduced epithelial cell-induced increase of mitochondrial mass in fibroblasts. p-: phosphorylated; ERK1/2: extracellular signal-regulated kinase-1/2; t-: total; PRMT1: protein arginine methyltransferase-1; PGC1α: peroxisome proliferator-activated receptor-γ coactivator-1α; GAPDH: glyceraldehyde 3-phosphate dehydrogenase; CM: control medium; BBT: before BT; ABT: after BT; Epi.S: epithelial cell culture supernatant; DAPI: 4′,6-diamidino-2-phenylindole. a) ABT Epi.S did not induce the expression and phosphorylation of ERK1/2 mitogen-activated protein kinase (p-ERK1/2 or t-ERK1/2), PRMT1 and PGC1α compared with BBT Epi.S. Five Western blots of each group were analysed by image quantification. Data are presented as mean±sem. *: p<0.05; **: p<0.01; ***: p<0.001. b, c) Representative immunofluorescence of b) nuclear PRMT1 accumulation and c) PGC1α in healthy fibroblasts in the presence of BBT Epi.S and ABT Epi.S. d) Representative immunofluorescence of mitochondrial mass in healthy fibroblasts, determined by MitoTracker 48 h after exposure to BBT Epi.S and ABT Epi.S. Enlarged images of single cells are depicted in the right panel. Scale bars: 10 µm.
Epi.S of cells obtained before BT significantly increased fibroblast proliferation (figure 6a) and their overall cell activity (figure 6b) when compared with Epi.S of cells obtained after BT. There was no difference comparing the response of fibroblasts from nonasthma with asthma cell donors and therefore the data were combined. The deposition of collagen type I and fibronectin was not different comparing Epi.S obtained before and after BT. Under both conditions, Epi.S significantly increased the deposition of collagen type I and fibronectin (figure 6c and d).
Bronchial thermoplasty (BT) reduced fibroblast proliferation and cell activation by bronchial epithelial cells. CM: control medium; BBT: before BT; ABT: after BT; Epi.S: epithelial cell culture supernatant; AU: arbitrary units. a) Direct cell counts after 48 h of incubating fibroblasts with either BBT Epi.S or ABT Epi.S (n=6). b) Fibroblast activity measured by MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide) over 48 h (n=4). c, d) BT had no effect on the deposition of c) collagen type I (COL1A1) (n=4) or d) fibronectin (n=4). Data are presented as mean±sem. **: p<0.01; ***: p<0.001, comparing control versus other groups; #: p<0.05, comparing BBT versus ABT.
BT reduced mitochondrial mass in fibroblasts by reducing C/EBPβ and PRMT1 expression
Epi.S obtained before BT significantly increased the expression of the mitochondria activators C/EBPβ, PGC1α and PRMT1 by fibroblasts of healthy tissue donors (figure 7a; for quantitative analysis, see supplementary figure S5a). The stimulatory effect of Epi.S obtained before BT was significantly reduced in fibroblasts pre-treated with C/EBPβ siRNA, as shown by Western blot (figure 7a) and immunofluorescence (figure 7b). Downregulation of C/EBPβ also reduced the stimulatory effect of Epi.S obtained before BT on mitochondrial mass (MitoTracker) in control fibroblasts (figure 7b).
Bronchial thermoplasty (BT) reduced the stimulatory effect of epithelial cell culture supernatant (Epi.S) obtained before BT (BBT) on mitochondrial mass in fibroblasts through CCAAT enhancer-binding protein-β (C/EBPβ) and protein arginine methyltransferase-1 (PRMT1). PGC1α: peroxisome proliferator-activated receptor-γ coactivator-1α; GAPDH: glyceraldehyde 3-phosphate dehydrogenase; si: small interfering; Con: control. a) In healthy fibroblasts (n=3) BBT Epi.S increased the expression of C/EBPβ, PGC1α and PRMT1, which was downregulated by siRNA to C/EBPβ (siC/EBPβ) within 48 h. The quantitative analysis of all Western blots is provided in supplementary figure S5a. b) Representative immunofluorescence of healthy fibroblasts in the presence of siC/EBPβ with and without BBT Epi.S, showing the expression of C/EBPβ, PRMT1 and MitoTracker. Scale bars: 10 µm. c) Downregulation of C/EBPβ significantly reduced the proliferative effect of BBT Epi.S on fibroblasts of asthma patients and controls (n=3). Data are presented as mean±sem. *: p<0.05; **: p<0.01, comparing control versus other groups.
Furthermore, downregulation of C/EBPβ significantly reduced the proliferative effect of Epi.S obtained before BT on fibroblasts of asthma patients and controls (figure 7c).
Epi.S obtained before BT significantly upregulated C/EBPβ, PRMT1, PGC1α, cytochrome c, HSP60 and α-SMA expression by fibroblasts of both asthma patients and controls (figure 8a; for quantitative analysis, see supplementary figure S5b). In fibroblasts treated with Epi.S obtained before BT, AMI-1 significantly reduced the expression of HSP60, cytochrome c and α-SMA (figure 8a), as well as C/EBPβ expression and MitoTracker signals (figure 8b). Similarly, downregulation of PRMT1 by siRNA reduced the expression of PGC1α and C/EBPβ in control fibroblasts (figure 8c).
Protein arginine methyltransferase-1 (PRMT1) inhibition reduced mitochondrial mass and regulators. C/EBPβ: CCAAT enhancer-binding protein-β; PGC1α: peroxisome proliferator-activated receptor-γ coactivator-1α; HSP60: heat shock protein-60; α-SMA: smooth muscle actin; GAPDH: glyceraldehyde 3-phosphate dehydrogenase; Epi.S: epithelial cell culture supernatant; BT: bronchial thermoplasty; BBT: before BT; ABT: after BT; si: small interfering; Con: control; CM: control medium. a) BBT Epi.S-induced expression of mitochondria activators was reduced by AMI-1. The quantitative analysis of all Western blots (n=4 in each group) is provided in supplementary figure S5b. b) Representative immunofluorescence of BBT Epi.S-induced expression of C/EBPβ and MitoTracker in the presence and absence of AMI-1. Scale bars: 10 µm. c) Representative Western blot of siRNA for PRMT1 (siPRMT1) on the expression of PGC1α and C/EBPβ. The quantitative analysis of all Western blots (n=3) is shown in the bar chart below the blot. d) The effect of AMI-1 on BBT Epi.S-induced fibroblast proliferation (n=3). Data are presented as mean±sem. *: p<0.05; **: p<0.01, comparing control versus other groups.
Blocking PRMT1 activity by AMI-1 also significantly reduced the pro-proliferative effect of Epi.S obtained before BT in fibroblasts from asthma patients and controls (figure 8d).
Proteomics analysis of BALF and NGS analysis in primary epithelial cells indicated that HSP60 expression is reduced after BT
14 BALF samples were analysed by high-resolution tandem mass spectrometry that included extensive peptide fractionation and phospho-peptide enrichment. In total, 1137 proteins were identified in the BALF (supplementary figures S3 and S4). Proteins with low abundance values (mean <1) were excluded before data analysis by OmicShare tools. Compared with BALF obtained before BT, 49 differentially expressed proteins (DEPs) were identified in BALF obtained after BT. Among these, 28 DEPs were upregulated and 21 DEPs were downregulated (table 2). Most significantly, HSP60 (HSPD1) was downregulated in BALF obtained after BT compared with BALF obtained before BT; in contrast, mucin (MUC5A) and complement C5 were upregulated (figure 9a).
Mass spectrometric analysis of the identified differentially expressed proteins in bronchoalveolar lavage fluid from patients before and after bronchial thermoplasty (BBT and ABT)
Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of differentially expressed proteins from bronchoalveolar lavage fluid
Proteomics analysis of bronchoalveolar lavage fluid (BALF) using the OmicShare tools indicated that heat shock protein-60 (HSP60) controls protein arginine methyltransferase-1 (PRMT1) expression. BT: bronchial thermoplasty; BBT: before BT; ABT: after BT; p-: phosphorylated; ERK1/2: extracellular signal-regulated kinase-1/2; t-: total; C/EBPβ: CCAAT enhancer-binding protein-β; GAPDH: glyceraldehyde 3-phosphate dehydrogenase; PGC1α: peroxisome proliferator-activated receptor-γ coactivator-1α; miR-19a: microRNA-19a; si: small interfering. a) Heat map of differentially expressed proteins contained in BBT BALF from seven patients compared with ABT BALF. Each row shows the expression of a single protein and each column represents a single BALF sample. The colour scale indicates high, medium and low protein expression. b) Kyoto Encyclopedia of Genes and Genomes analysis of 49 differentially expressed proteins (table 3). c) Cytoscape map of representative signalling pathway networks and related genes. d) Kinetics of HSP60-stimulated p-ERK1/2, t-ERK1/2 and C/EBPβ expression in control fibroblasts demonstrated by Western blot. The quantitative analysis of all Western blots is provided in supplementary figure S5c. e) Primary fibroblasts of control patients (n=3) were treated with 1 μg·mL−1 recombinant human HSP60 for 48 h. The expression of C/EBPβ, PRMT1 and PGC1α was detected by Western blot. The quantitative analysis of all Western blots is provided in supplementary figure S5c. f) Graphic summary of the study showing that BT interrupts pro-remodelling signalling through PRMT1.
To determine the effect of BT on protein/gene regulatory pathways, KEGG analysis was performed. KEGG analysis suggested that most BT-regulated DEPs were linked to signal transduction (EEF2, YWHAE, FN1, CAT, RUVBL1 and ATP1A2), infectious diseases (HSPD1, VCP, SERPINB1, FN1, HSPA9, C5 and YWHAE) and the immune system (S100A9, MUC5AC and C5) (figure 9b and c, and table 3).
Transcriptome analysis of primary epithelial cells confirmed the modification of HSP60 expression (table 4). Therefore, the role of HSP60 in the regulation of PRMT1 expression was further assessed in lung fibroblasts. Stimulation of fibroblasts with human recombinant HSP60 (1 μg·mL−1) significantly increased the C/EBPβ and ERK1/2 MAPK signalling cascade within 2 h, which stimulated the expression of PRMT1 and PGC1α at 48 h (figure 9d and e; for quantitative Western blot analysis, see supplementary figure S5c).
Transcriptome analysis
Figure 9f summarises the interaction of epithelial cells and fibroblasts, and its modification by BT.
Discussion
BT is the only therapy shown to reduce remodelling in patients with severe asthma; however, the mechanism is unknown [6, 7, 12]. This study provides the first evidence that BT produces lasting effects on bronchial epithelial cell function and thereby changes some remodelling properties of airway wall fibroblasts. BALF and Epi.S of cells isolated from patients with severe asthma before BT reduced the expression of miR-19a in fibroblasts, thereby stimulating the expression of ERK1/2 MAPK, PGC1α and PRMT1, and thus increasing mitochondrial mass and fibroblast proliferation. Importantly, these effects were lost after BT.
In an animal asthma model, PRMT1 was upregulated first in epithelial cells and later in mesenchymal cells, by cell type specific cytokines [19, 20, 22]. In tissue sections of asthma patients obtained before BT, phosphorylation of ERK1/2 MAPK and PRMT1 expression were high in epithelial cells and lower in subepithelial mesenchymal cells. ERK1/2 MAPK phosphorylation and PRMT1 expression were lost in tissue sections obtained after BT. It should be noted that the tissue biopsies were obtained 4–12 weeks after BT and thus the effect was lasting. Our findings suggest a central role for PRMT1 in the regulation of asthma-associated airway fibroblast remodelling.
Bronchial epithelial cells isolated after BT proliferated significantly faster than epithelial cells isolated from the same patients before BT, indicating improved epithelium repair capacity. These findings may explain improved re-epithelisation and structure adjustment, which were reported in clinical studies after BT [12, 24]. The response of fibroblasts to BALF or Epi.S suggested that the beneficial effect of BT depends on the modified function of epithelial cells.
In human asthmatic mesenchymal cells, the constitutive expression of PRMT1 was caused by decreased levels of miR-19a and upregulated ERK1/2 MAPK, subsequently activating signal transducer and activator of transcription-1 (STAT1) that increased both PRMT1 expression and mitochondrial mass [17, 21, 23]. Similarly, BALF and Epi.S obtained before BT downregulated miR-19a, leading to increased ERK1/2 MAPK in fibroblasts. This study also confirmed that C/EBPβ is a negative regulator of miR-19a, as reported in other cell types [24, 25]. Furthermore, our observations suggest a feedback mechanism (C/EBPβ→miR-19a→ERK1/2 MAPK→PRMT1→C/EBPβ) leading to fibroblast remodelling in asthma. Most importantly, such a mechanism can be lastingly broken by BT.
Mitochondria function and biogenesis were required to drive airway wall remodelling in asthma and COPD [26–28]. Here, the mitochondria regulating proteins PGC1α, cytochrome c and HSP60 were increased in fibroblasts by Epi.S obtained before BT and were downregulated after BT.
In vitro, BT reduced smooth muscle mass and α-SMA in myofibroblasts [13]. Salem et al. [5] described that BT, over a period of 27 months, reduced airway smooth muscle cell mass, α-SMA expression and basement membrane thickness in 10 patients with severe asthma. In this study, Epi.S obtained before BT increased α-SMA expression and mitochondrial mass in fibroblasts, and these effects were lost in Epi.S obtained after BT. This observation may help to explain the lasting reduction of subepithelial mesenchymal cell numbers reported after BT [5, 11, 12]. The fact that in this study epithelial cells isolated after BT lost their stimulatory properties suggests that BT induced a lasting “re-programming”.
Proteome analysis of BALF and transcriptome analyses in isolated epithelial cells suggested that HSP60 mediated most BT effects. HSP60 contributed to inflammation by activating immune cells [29] and inducing cell dysfunction in obese patients [30, 31]. In this study, exogenous HSP60 increased PRMT1 expression, the corresponding signalling cascade and mitochondrial mass in fibroblasts. Our results indicated that BALF and Epi.S of asthma patients obtained before BT caused mitochondrial dysfunction through upregulation of PGC1α→PRMT1, which is consistent with earlier reports [32, 33].
The strength of this study is to provide a molecular biological mechanism that explains the lasting beneficial effect of BT on some parameters of airway wall remodelling.
The limitations of this study include: 1) further studies assessing the consequences of BT-induced changes on protein expression on target cell function, such as smooth muscle cells, have to be performed; 2) the observed BT-induced changes in the proteome and transcriptome should be confirmed in tissue sections from a larger number of asthma patients; 3) mass spectrometry could not be performed for all conditions, due to the presence of bovine albumin in the cell culture medium which obscured the reading (therefore we were not able to confirm that the secreted factors found in BALF originated from epithelial cells); and 4) we could not investigate all BT-induced changes due to limited availability of tissue sections after BT.
In conclusion, we present a novel self-sustaining, remodelling driving mechanism in fibroblasts of asthma patients. Importantly, this mechanism could be lastingly broken by BT. Furthermore, the findings suggest that HSP60, miR-19a or PRMT1 may present new therapeutic targets to control airway wall remodelling in asthma.
Supplementary material
Supplementary Material
Please note: supplementary material is not edited by the Editorial Office, and is uploaded as it has been supplied by the author.
Supplementary figure S1. Representative immuno-cytochemical characterisation of primary human bronchial epithelial cells, fibroblasts and smooth muscle cells by phenotype using phase-contrast microscopy, fibronectin, α-smooth muscle actin, and cytokeratin-14 (only expressed by epithelial cells), the white scale bar: 10 μm. ERJ-00300-2019.Figure_S1
Supplementary figure S2. Immuno-histochemical staining for PRMT1 in bronchial tissue sections obtained a) before BT and b) two years after BT in the same patient. Similar results were obtained in a second patient. ERJ-00300-2019.Figure_S2
Supplementary figure S3. 14 BALF samples were selected for proteomic analysis. The cohort included BALF samples obtained from seven asthma patients BBT and ABT, and shows the median and the 95% confidence intervals. Samples were analysed by high-resolution, accurate mass, tandem mass spectrometry (MS) that included extensive peptide fractionation and phosphopeptide enrichment, the protein concentration and peptide concentration are qualified for MS. ERJ-00300-2019.Figure_S3
Supplementary figure S4: The heat map of all proteins in BALF sample from 13 patients before or after BT. ERJ-00300-2019.Figure_S4
Supplementary figure S5: Quantitative analysis of protein expressions for Western-blots in figures 7a, as well as for figure 8a, 8d and 8e. a) Protein expression was quantified for all Western-blots that contributed to figure 7a for: C/EBP-β, PRMT1 and PGC-1α. b) Protein expression was quantified for all Western-blots that contributed to figure 8a for: C/EBP-β, PRMT1, PGC-1α, Cyt C, HSP60 and α-SMA. c) Protein expression was quantified for all Western-blots that contributed to figures 8d and 8e for: p-ERK/total ERK, C/EBP-β, PRMT1 and PGC-1α. Bars present mean±SEM of corresponding measurements (“n” as stated above). t-test: *: p<0.05, **: p<0.01, ***: p<0.001 comparing control and other groups, or indicated groups. ERJ-00300-2019.Figure_S5
Shareable PDF
Supplementary Material
This one-page PDF can be shared freely online.
Shareable PDF ERJ-00300-2019.Shareable
Acknowledgements
We are grateful to Xinggang Wang (Reproductive Medical Centre, Jiaotong University, Shanghai, China) and Hui Wang (Stem Cell Research, University Hospital Basel, Basel, Switzerland) for their expertise assistance, helpful discussions and productive critiques. We also thank C.T. S'ng (PersCellMed S'ng, Sissach, Switzerland) for his help in preparing this manuscript.
Footnotes
This article has been revised according to the correction published in the July 2021 issue of the European Respiratory Journal.
This article has supplementary material available from erj.ersjournals.com
Author contributions: Q. Sun, L. Fang, D. Stolz, S. Lu and M. Roth conceived and designed the study. Q. Sun, L. Fang, X. Tang, M. Roth and E. Papakonstantinou performed experiments, and analysed the data. W. Zhai and Q. Sun analysed the data of BALF proteins by KEGG pathway analysis. R. Louis, V. Heinen, F.N. Schleich, D. Stolz and M. Tamm provided the lung tissues, patients' clinical data, and contributed to manuscript preparation and writing. S. Savic contributed lung tissue section for histology. Q. Sun, S. Lu, M. Tamm, D. Stolz and M. Roth had scientific discussion of the study and wrote the paper. All authors read and approved the manuscript.
Conflict of interest: Q. Sun has nothing to disclose.
Conflict of interest: L. Fang has nothing to disclose.
Conflict of interest: M. Roth has nothing to disclose.
Conflict of interest: X. Tang has nothing to disclose.
Conflict of interest: E. Papakonstantinou has nothing to disclose.
Conflict of interest: W. Zhai has nothing to disclose.
Conflict of interest: R. Louis reports grants and personal fees for advisory board work from GSK, AstraZeneca and Novartis, and grants from Chiesi, outside the submitted work.
Conflict of interest: V. Heinen has nothing to disclose.
Conflict of interest: F.N. Schleich has nothing to disclose.
Conflict of interest: S. Lu has nothing to disclose.
Conflict of interest: S. Savic has nothing to disclose.
Conflict of interest: M. Tamm has nothing to disclose.
Conflict of interest: D. Stolz has nothing to disclose.
Support statement: This work was supported by Swiss National Foundation grant 31003A_176248/1, Swiss Lung Association grant 2017-02_Roth to M. Roth, a Swiss Government Excellence Scholarship to Q. Sun 2014.0405, and National Natural Science Foundation of China grants 31871314 and 31601033 to S. Lu. Funding information for this article has been deposited with the Crossref Funder Registry.
- Received May 16, 2018.
- Accepted August 13, 2019.
- Copyright ©ERS 2019