Skip to main content

Main menu

  • Home
  • Current issue
  • ERJ Early View
  • Past issues
  • ERS Guidelines
  • Authors/reviewers
    • Instructions for authors
    • Submit a manuscript
    • Open access
    • Peer reviewer login
    • WoS Reviewer Recognition Service
  • Alerts
  • Subscriptions
  • ERS Publications
    • European Respiratory Journal
    • ERJ Open Research
    • European Respiratory Review
    • Breathe
    • ERS Books
    • ERS publications home

User menu

  • Log in
  • Subscribe
  • Contact Us
  • My Cart

Search

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

Login

European Respiratory Society

Advanced Search

  • Home
  • Current issue
  • ERJ Early View
  • Past issues
  • ERS Guidelines
  • Authors/reviewers
    • Instructions for authors
    • Submit a manuscript
    • Open access
    • Peer reviewer login
    • WoS Reviewer Recognition Service
  • Alerts
  • Subscriptions

Identification of asthma-associated microRNAs in bronchial biopsies

Mirjam P. Roffel, Ilse M. Boudewijn, Jos L.L. van Nijnatten, Alen Faiz, Corneel J. Vermeulen, Antoon J. van Oosterhout, Karen Affleck, Wim Timens, Ken R. Bracke, Tania Maes, Irene H. Heijink, Corry-Anke Brandsma, Maarten van den Berge
European Respiratory Journal 2022 59: 2101294; DOI: 10.1183/13993003.01294-2021
Mirjam P. Roffel
1Dept of Pathology and Medical Biology, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
2Groningen Research Institute for Asthma and COPD, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
3Dept of Respiratory Medicine, Laboratory for Translational Research in Obstructive Pulmonary Diseases, Ghent University, University Hospital Ghent, Ghent, Belgium
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ilse M. Boudewijn
2Groningen Research Institute for Asthma and COPD, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
4Dept of Pulmonology, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jos L.L. van Nijnatten
1Dept of Pathology and Medical Biology, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
2Groningen Research Institute for Asthma and COPD, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
5Faculty of Science, Respiratory Bioinformatics and Molecular Biology (RBMB), University of Technology Sydney, Sydney, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jos L.L. van Nijnatten
Alen Faiz
1Dept of Pathology and Medical Biology, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
2Groningen Research Institute for Asthma and COPD, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
5Faculty of Science, Respiratory Bioinformatics and Molecular Biology (RBMB), University of Technology Sydney, Sydney, Australia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Corneel J. Vermeulen
2Groningen Research Institute for Asthma and COPD, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
4Dept of Pulmonology, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Corneel J. Vermeulen
Antoon J. van Oosterhout
6Allergic Inflammation Discovery Performance Unit, GlaxoSmithKline, Stevenage, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Karen Affleck
7Adaptive Immunity Research Unit, GlaxoSmithKline, Stevenage, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Wim Timens
1Dept of Pathology and Medical Biology, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
2Groningen Research Institute for Asthma and COPD, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Wim Timens
Ken R. Bracke
3Dept of Respiratory Medicine, Laboratory for Translational Research in Obstructive Pulmonary Diseases, Ghent University, University Hospital Ghent, Ghent, Belgium
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tania Maes
3Dept of Respiratory Medicine, Laboratory for Translational Research in Obstructive Pulmonary Diseases, Ghent University, University Hospital Ghent, Ghent, Belgium
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Irene H. Heijink
1Dept of Pathology and Medical Biology, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
2Groningen Research Institute for Asthma and COPD, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
4Dept of Pulmonology, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Corry-Anke Brandsma
1Dept of Pathology and Medical Biology, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
2Groningen Research Institute for Asthma and COPD, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
8Both C-A. Brandsma and M. van den Berge are senior authors and contributed equally to this work
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Maarten van den Berge
2Groningen Research Institute for Asthma and COPD, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
4Dept of Pulmonology, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
8Both C-A. Brandsma and M. van den Berge are senior authors and contributed equally to this work
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: m.van.den.berge@umcg.nl
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Abstract

Background Changes in microRNA (miRNA) expression can contribute to the pathogenesis of many diseases, including asthma. We aimed to identify miRNAs that are differentially expressed between asthma patients and healthy controls, and explore their association with clinical and inflammatory parameters of asthma.

Methods Differentially expressed miRNAs were determined by small RNA sequencing on bronchial biopsies of 79 asthma patients and 82 healthy controls using linear regression models. Differentially expressed miRNAs were associated with clinical and inflammatory asthma features. Potential miRNA–mRNA interactions were analysed using mRNA data available from the same bronchial biopsies, and enrichment of pathways was identified with Enrichr and g:Profiler.

Results In total, 78 differentially expressed miRNAs were identified in bronchial biopsies of asthma patients compared with controls, of which 60 remained differentially expressed after controlling for smoking and inhaled corticosteroid treatment. We identified several asthma-associated miRNAs, including miR-125b-5p and miR-223-3p, based on a significant association with multiple clinical and inflammatory asthma features and their negative correlation with genes associated with the presence of asthma. The most enriched biological pathway(s) affected by miR-125b-5p and miR-223-3p were inflammatory response and cilium assembly/organisation. Of interest, we identified that lower expression of miR-26a-5p was linked to more severe eosinophilic inflammation as measured in blood, sputum as well as bronchial biopsies.

Conclusion Collectively, we identified miR-125b-5p, miR-223-3p and miR-26a-5p as potential regulators that could contribute to the pathogenesis of asthma.

Abstract

Using small RNA sequencing on bronchial biopsies of asthma patients and healthy controls, this study identified microRNAs associated with clinical and inflammatory features of asthma that are potential regulators of asthma-associated gene expression https://bit.ly/3rZXMDf

Introduction

Asthma is a chronic inflammatory airway disease characterised by airway obstruction and airway hyperresponsiveness resulting in shortness of breath, cough, wheezing and chest tightness. Around 300 million people worldwide are affected by asthma [1]. Although asthma symptoms can be controlled by inhaled corticosteroids (ICS) and long-acting β-agonists, there is no cure available and optimal asthma control cannot be obtained in all patients with currently available treatments [1]. Therefore, there is a need to better understand the pathogenesis of asthma to identify new therapeutic targets and improve the management of asthma patients.

Genetics, epigenetics and microRNAs (miRNAs) are involved in the pathogenesis of asthma [2]. miRNAs are small noncoding RNA molecules that target specific mRNAs leading to their degradation or translational repression. Several miRNAs can target the same mRNA and each miRNA can regulate hundreds of mRNAs [3]. Therefore, miRNAs have the potential to regulate many cellular functions, including inflammation, cell differentiation and cell death [4]. miRNAs can provide novel targets for new therapeutic treatment strategies, because miRNAs are small, highly conserved and have a specific sequence. As an example, blocking the function of miR-155 with an miR-155 antagomir has been shown to overcome cancer cell drug resistance in a mouse model [5].

Several groups have investigated miRNA profiles in bronchial biopsies, brushings, blood and sputum of asthma patients compared with healthy controls [6–11]. However, generally small cohorts were used and most studies were based on microarray technology, which is less sensitive and specific compared with sequencing. In this study, we performed small RNA sequencing with matched RNA sequencing in the same bronchial biopsies of 79 well-characterised asthma patients and 82 matched healthy individuals. We aimed to identify miRNAs involved in the pathogenesis of asthma. We focused on differentially expressed miRNAs and miRNA–mRNA correlations that are associated with clinical and inflammatory asthma features.

Material and methods

Subject characteristics

Asthma patients (n=79) had a doctor's diagnosis of asthma and presence of bronchial hyperresponsiveness to AMP or histamine. Healthy subjects (n=82) had normal pulmonary function, normal spirometry with forced expiratory volume in 1 s (FEV1) >80% predicted and FEV1/forced vital capacity (FVC) greater than the lower limit of normal, absence of reversibility (FEV1 % pred to salbutamol <10%), no bronchial hyperresponsiveness (provocative concentration causing a 20% fall in FEV1 (PC20) methacholine >16 mg·mL−1), and no respiratory symptoms [12]. Patients were included from three studies [12, 13], all carried out in the University Medical Center Groningen (Groningen, The Netherlands) between 2001 and 2012, and approved by the local Medical Ethics Committee (METc 2009/007, 2001/074 and 2004/271). All subjects provided written informed consent.

Clinical measurements

Sputum induction, blood sampling, lung function and bronchial biopsies were performed in all subjects as previously described [12].

RNA isolation

Total RNA was extracted from frozen bronchial biopsies using the AllPrep DNA/RNA/miRNA Universal Kit (Qiagen, Venlo, The Netherlands), according to the manufacturer's instructions. Quality of RNA was assessed using the NanoDrop 1000 and Labchip GX platform (PerkinElmer, Waltham, MA, USA). Total RNA was extracted at the same time, and asthma patients and controls were evenly distributed over batches of 12 samples based on disease/control status, sex, age and smoking status.

miRNA sequencing

Small RNA sequencing library preparation was performed using the NEXTFLEX Small RNA-Seq Kit V3 (Bioo-Scientific, Austin, TX, USA). Clean-up beads were used to eliminate the small RNA fraction and quality control of the RNA library was checked using LabChip GX (Perkin Elmer, Waltham, MA, USA). Sequencing was performed using the HiSeq 2500 system (Illumina, San Diego, CA, USA). Quality control of raw RNA sequencing data was performed using FastQC (version 0.11.5) and quality control of raw reads of adapter sequences was done using TrimGalore (version 0.3.7). NEXTFLEX Small RNA-Seq Kit V3 processing was done using custom scripts documented in the GitHub repository. Alignment and quantification were performed using miRDeep2 (version 2.0.0.8) with Bowtie (version 0.12.7). During library preparation and sequencing, the samples were evenly distributed over the batches based on disease/control status, sex, age and smoking status.

mRNA sequencing

Sample preparation and mRNA sequencing were performed as previously described [14]. Library preparation for mRNA and miRNA sequencing was performed at the same time for asthma patients and healthy controls, and samples were evenly distributed over the batches.

Statistical analyses of small RNA sequencing data

The (small) RNA sequencing data analyses were performed with R (version 3.5.3). miRNAs with an average read count <100 in asthma patients or healthy controls were filtered out. Differential expression of microRNAs was analysed with DESeq2 (version 1.24.0) adjusting for sex, age, smoking and library batch, and corrected for multiple testing using the Benjamini–Hochberg false discovery rate (FDR); FDR adjusted p-values <0.05 were considered as statistically significant. A heatmap was generated and ordered from highest to lowest negative fold change miRNA in asthma patients compared with healthy controls. To exclude potential effects of ICS and smoking, a sensitivity analysis was performed in a subset of asthma patients (n=30) compared with healthy controls (n=42) excluding participants that smoked or used ICS.

Linear regression was performed to associate differentially expressed miRNAs with relevant clinical and inflammatory features of asthma: FEV1 % pred, PC20 AMP, eosinophil and neutrophil counts in blood and tissue, and percentages of eosinophils and neutrophils in sputum in all asthma patients adjusting for age, sex, smoking status and ICS treatment.

To assess for miRNA–mRNA correlations, expression levels were log2 transformed into counts per million (logCPM). Pearson's correlations for miRNA–mRNA were calculated for miRNAs with one or more associations with asthma features and genomewide gene expression data available in matched subjects (FDR adjusted p-values <0.05 were considered statistically significant). Negatively correlated genes were checked for their association with asthma using the same statistical model in matched biopsies; these were defined as asthma genes. For those miRNAs, we checked for experimentally validated target genes in miRTarBase (version 7.0) [15].

Gene Ontology and pathway analyses

For each miRNA, the negatively or the positively and negatively correlated asthma-associated genes were analysed with Enrichr [16, 17] and g:Profiler [18] to identify enriched biological processes (Gene Ontology (GO) Biological Process 2018 in Enrichr and GO:BP releases/2020-12-15 in g:Profiler) and Reactome pathways (Reactome_2016 in Enrichr and Reactome_2020-12-15 in g:Profiler).

Results

Patients characteristics

The clinical characteristics of 79 asthma patients and 82 healthy controls are shown in table 1. There were no significant differences between asthma patients and healthy controls with respect to sex and pack-years of smoking. Asthma patients had significantly higher age, were less often current smokers, had a lower FEV1 % pred and FEV1/FVC, and had more severe bronchial hyperresponsiveness to AMP.

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

Characteristics of asthma patients and healthy controls

Differentially expressed miRNAs in asthma patients compared with healthy controls

Of the 1860 miRNAs that were detected in bronchial biopsies, 136 miRNAs remained after filtering for low abundance. Of these 136, 78 miRNAs were differentially expressed in asthma patients compared with healthy controls: 48 having higher expression and 30 having lower expression in asthma (FDR adjusted p-value <0.05) (figure 1a). miR-451a had the highest fold change, whereas miR-125b-5p had the lowest negative fold change in asthma. A volcano plot illustrates all differentially expressed miRNAs, with the top five most significant highest and lowest miRNAs indicated (FDR adjusted p-value <0.05) (figure 1b). To assess whether smoking and ICS treatment influenced the 78 differentially expressed miRNAs, we conducted a sensitivity analysis on a subgroup of individuals that did not smoke and did not use ICS. Of the 78 miRNAs that were differentially expressed between asthma patients and healthy controls in the complete cohort, 60 miRNAs remained significant in the same direction in the sensitivity analysis, showing that our findings are robust (FDR adjusted p-value <0.05) (figure 1c). Supplementary tables S1 and S2 present the full lists of significantly differentially expressed miRNAs in asthma patients compared with healthy controls, including and excluding smokers and ICS-treated patients.

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

MicroRNAs (miRNAs) that were differentially expressed in bronchial biopsies of asthma patients compared with bronchial biopsies of healthy controls. a) Heatmap showing 78 differentially expressed miRNAs in asthma patients (n=79) compared with healthy controls (n=82). b) Volcano plot showing differentially expressed miRNAs (false discovery rate (FDR) adjusted p<0.05) between asthma patients and healthy controls including smokers and inhaled corticosteroid (ICS) treatment. c) Volcano plot showing the differentially expressed miRNAs (FDR adjusted p<0.05) between asthma patients (n=30) and healthy controls (n=42) excluding smokers and ICS treatment. The top five upregulated (red) and downregulated (blue) miRNAs are indicated. The dotted line indicates the significance level.

Expression of 25 differentially expressed miRNAs is associated with asthma features

To determine whether the 60 differentially expressed miRNAs in bronchial biopsies are relevant in asthma, we investigated the association of differentially expressed miRNAs with FEV1 % pred, PC20 AMP, eosinophil and neutrophil counts in blood and tissue, and percentages of eosinophils and neutrophils in sputum of all asthma patients. We observed that 25 out of the 60 differentially expressed miRNAs were associated with one or more clinical and inflammatory feature(s). Figure 2 illustrates the association between miRNAs and features of asthma with T-values of the linear regression (p<0.05). Higher miR-21-5p expression in biopsies had the most significant association with more neutrophils in blood, and it was also associated with more neutrophils and eosinophils in tissue. Furthermore, miR-26a-5p had three significant associations with asthma features: lower miR-26a-5p expression was associated with higher eosinophil levels in blood, sputum and tissue. Moreover, higher expression levels of miR-223-3p were linked to higher eosinophils and neutrophils in tissue. Additionally, we observed that lower miR-125b-5p levels, the lowest negative fold change miRNA in asthma, were associated with lower FEV1 % pred and higher neutrophil levels in blood. For miR-451a, the miRNA with the highest fold change in asthma, we did not observe any associations with asthma features in our study.

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

Linear regression between differentially expressed microRNAs (miRNAs) in bronchial biopsies and clinical parameters. Positive associations are displayed in red and negative associations are displayed in blue (n=79, all asthma patients). FEV1: forced expiratory volume in 1 s; PC20: provocative concentration causing a 20% fall in FEV1. *: p<0.05; **: p<0.01; ***: p<0.001.

Identification of key asthma-associated miRNA–mRNA targets

To investigate potential interactions between miRNA and mRNA expression, we correlated expression levels of the 25 miRNAs that were associated with one or more clinical and inflammatory asthma features with mRNA sequencing data available from the same biopsies from all asthma patients and healthy controls. Hereby, we focused on the negatively correlated genes. The expression of 23 of the 25 miRNAs was negatively correlated with mRNA expression of one or more genes (table 2). Thereafter, we investigated which of the negatively correlated genes were also associated with the presence of asthma. For miR-199b-5p, miR-223-3p, miR-199b-3p, miR-142-5p, miR-181b-5p, miR-195-5p and miR-125b-5p, we observed that ∼10% of the negatively correlated genes were also associated with asthma. A full list of negatively correlated genes per miRNA is provided in supplementary table S3. As the expression of these seven miRNAs (miR-199b-5p, miR-223-3p, miR-199b-3p, miR-142-5p, miR-181b-5p, miR-195-5p and miR-125-5p) was associated with one or more clinical and inflammatory asthma features and ∼10% of their negatively correlated genes were associated with asthma genes, we selected those miRNAs as our candidate asthma miRNAs. For the five lower expressed miRNAs in asthma in our study, i.e. miR-199b-5p, miR-199b-3p, miR-181b-5p, miR-195-5p and miR-125b-5p, we observed that the genes HBA1 and RGS18 had the most significant negative correlation (figure 3a and b). HBA1 and RGS18 expression levels were higher in bronchial biopsies of asthma patients compared with healthy controls. Furthermore, for the two higher expressed asthma miRNAs in our study, i.e. miR-223-3p and miR-142-5p, we observed that GREM2 and RPS3AP5 were the most significant negatively correlated genes (figure 3c and d). GREM2 expression levels were lower in bronchial biopsies of asthma patients compared with healthy controls, while no expression differences were found for RPS3AP5. In addition, among the negatively correlated asthma-associated genes we identified three validated target genes for miR-199b-5p (AGTRAP, RNF11 and SNTB1), one validated target gene for miR-125b-5p (S100A8) and 13 validated target genes for miR-223-3p (ARL8B, ARTN, CAPRIN1, CDS1, CHMP2B, CHUK, MSMO1, NSUN3, PDZD8, SECISBP2L, TWF1, WASL and ZBTB18).

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

Differentially expressed microRNAs (miRNAs) in asthma patients compared with healthy controls that have a correlation with clinical and inflammatory asthma features: negatively correlated genes per miRNA

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

Correlation between asthma-associated genes and microRNA (miRNA) expression. Pearson's correlation coefficient (r) between a) HBA1, b) RGS18, c) GREM2 and d) RPS3AP5 expression correlated with miRNA expression in all subjects (n=79 asthma patients and n=82 healthy controls). The x- and y-axes represent normalised read counts.

Enrichment of biological processes and pathways in negatively correlated asthma genes

To identify the potential mechanisms underlying the involvement of our seven key candidate asthma miRNAs in disease pathogenesis, pathway enrichment analyses were performed on the negatively correlated asthma-associated genes, whereby we included all asthma patients and healthy controls. GO revealed in total 187 biological processes and 212 Reactome pathways that were significantly enriched among the negatively correlated mRNA transcripts (FDR adjusted p-value <0.05). The most significant biological process for miR-199b-5p and miR-223-3p was “cilium assembly”, for miR-181b-5p, miR-199b-3p and miR-195-5p “neutrophil degranulation”, and for miR-125b-5p “inflammatory response”. For miR-142-5p, no significant enrichment was found using Enrichr; however, using g:Profiler, four biological processes were obtained. The most significant biological process for miR-142-5p was “protein transport along microtubules”. The top five biological processes and Reactome pathways for each of the six asthma miRNAs are displayed in tables 3 and 4, respectively; a list of all biological processes and Reactome pathways is provided in supplementary tables S4 and S5. Additionally, we also performed pathway enrichment analyses on positively and negatively correlated asthma-associated genes. Using Enrichr, in total 56 biological processes and 51 pathways were significantly enriched among the positively and negatively correlated genes. The most significant biological process for miR-199b-5p and miR-223-3p remained “cilium assembly”, for miR-199b-3p “potassium ion transport”, for miR-142-5p “neutrophil degranulation”, and for miR-181b-5p, miR-195-5p and miR-125b-5p no significantly enriched biological processes were found. A full list of all biological processes and Reactome pathways identified by Enrichr and g:Profiler is provided in supplementary tables S6 and S7.

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

Top five Gene Ontology (GO) biological processes of negatively correlated asthma genes per microRNA

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

Top five Reactome pathways of negatively correlated asthma genes per microRNA

Discussion

In this study, we identified 78 differentially expressed miRNAs in bronchial biopsies of asthma patients compared with healthy controls. Associations between miRNA expression and clinical and inflammatory asthma features identified seven candidate “asthma miRNAs”, i.e. miR-199b-5p, miR-223-3p, miR-199b-3p, miR-142-5p, miR-181b-5p, miR-195-5p and miR-125b-5p, that were significantly associated with one or more clinical and inflammatory asthma features as well as asthma-associated gene expression. The biological pathways affected by these miRNAs included cilium assembly/organisation, neutrophil activation/degranulation and inflammatory response. Of interest, we identified that lower miR-26a-5p expression was strongly related to more eosinophilic inflammation.

The most significant lower expressed miRNA in asthma was miR-125b-5p, which was associated with more severe airflow obstruction and increased blood neutrophil counts. Enrichment of biological pathways revealed that miR-125b-5p is linked to inflammatory responses. Previously, Liu et al. [19] showed that miR-125b-5p is also lower expressed in sputum of childhood asthma patients versus controls, especially in those with more severe eosinophilic inflammation. Furthermore, in current smokers with or without chronic obstructive pulmonary disease, lower levels of miR-125b-5p were also linked with more airflow obstruction [20]. In our study, miR-125b-5p still remained lower expressed in asthma in the sensitivity analysis excluding smokers and patients that used ICS, indicating that in our cohort smoking had no major effect on the expression of miR-125b-5p. Additionally, functional studies demonstrated that intranasal administration of miR-125b can attenuate asthma features by reducing interleukin (IL)-4 and IL-13 levels, goblet cell differentiation, and mucus production in a murine asthma model [19].

Since asthma patients had higher miR-223-3p levels, and those levels were associated with higher tissue eosinophils and neutrophils and negatively correlated with asthma-associated genes, we suggest that miR-223-3p could play a role in the pathogenesis of asthma. In line with this, several studies showed higher levels of miR-223-3p in bronchial airway epithelial cells and sputum of asthma patients, especially in severe (neutrophilic) asthma patients [10, 21]. Furthermore, several studies demonstrated that overexpression of miR-223-3p reduced pro-inflammatory responses, while depletion of miR-223-3p enhanced inflammatory responses (reviewed in [22]). In our study, higher levels of miR-223-3p in biopsies were strongly associated with lower GREM2 expression. GREM2, which encodes for gremlin-2 and is lower expressed in asthma patients, is an antagonist of bone morphogenetic proteins [23]. Gremlin-2 is involved in abnormal tissue damage and repair responses, including lung fibrosis [24]. Additionally, the most significantly enriched biological process in our study for miR-223-3p is cilium assembly/organisation. Cilium assembly/organisation is important for the formation of cilia on bronchial epithelia cells, leading to mucociliary function and clearance of mucus, and preventing infections and inflammation [25]. We also observed that miR-223-3p is enriched in several neutrophilic processes, which is in line with previous studies [10, 21]. However, in our study these processes do not belong to the top 50 most significantly enriched biological pathways, which could be due to the differences in sources (biopsies versus sputum) and patient cohorts (mild to moderate versus severe asthma patients). Overall, this indicates that miR-223-3p may contribute to aberrant airway inflammation and repair responses by regulating (among others) GREM2 expression, suggesting a mechanistic role of miR-223-3p in asthma pathogenesis.

In addition, we identified miR-199b-5p, miR-181b-5p, miR-199b-3p and miR-195-5p as key candidate asthma miRNAs, since these miRNAs are associated with clinical and inflammatory asthma features and can play a potential role as regulators of asthma-associated gene expression. In this study, lower levels of miR-199b-5p, miR-181b-5p, miR-199b-3p and miR-195-5p were associated with higher levels of HBA1 in bronchial biopsies of asthma patients. Erythrocyte precursors are the major producer of HBA1, which encodes for the α subunit of haemoglobin [26]. Additionally, lung epithelial cells can produce the haemoglobin α subunit [27]. As yet, little is known about the role of HBA1 in clinical and inflammatory features of asthma. Furthermore, enrichment analyses with the negatively correlated asthma-associated genes of these miRNAs (miR-199b-5p, miR-181b-5p, miR-199b-3p and miR-195-5p) revealed that miR-199b-5p was involved in cilium assembly/organisation, and miR-181b-5p, miR-199b-3p and miR-195-5p were involved in neutrophil activation/degranulation. Regarding neutrophil activation/degranulation, neutrophilic inflammation is linked to severe, steroid-resistant asthma patients and to smoking asthmatic subjects [28, 29]. Using both negatively and positively correlated asthma-associated genes for enrichment analyses, it remains that miR-199b-5p was involved in cilium assembly/organisation. However, for miR-181b-5p and miR-195-5p no significant biological processes were found, and miR-199b-3p was revealed to be involved in potassium ion transport. Future functional studies should help to elucidate the role of these miRNAs in cilium organisation/assembly, neutrophil activation/degranulation and potassium ion transport.

Another miRNA of interest was miR-26a-5p. The expression of miR-26a-5p was lower in asthma patients and associated with more severe eosinophilic inflammation as measured in blood, sputum as well as bronchial biopsies. Other studies have reported that miR-26a was lower expressed in bronchial epithelial brushings, serum and exosomes of bronchoalveolar lavage fluid derived from asthma patients compared with controls [8, 30, 31]. Furthermore, a previous study linked lower miR-26a expression in serum of asthma patients to lower levels of FEV1 % pred [31]; however, we did not observe this association for miR-26a-5p in bronchial biopsies. An in vivo mice study showed that miR-26a/b regulates allergic inflammation by reducing the levels of cyclooxygenase-2 [32]. While the expression of cyclooxygenase-2 is increased in airways of asthma patients [33], we did not observe differences in cyclooxygenase-2 expression in bronchial biopsies of asthma patients compared with healthy controls and also found no correlation between miR-26a-5p and cyclooxygenase-2. Altogether, these data indicate that there might be a potential role for miR-26a-5p in the pathogenesis of asthma, especially eosinophilic asthma.

Of interest, miR-451a was the miRNA with the highest fold change in bronchial biopsies of asthma patients, while no correlation was found with clinical features of asthma. This may suggest that the association of this miRNA with asthma is not driven by any of the well-known asthma features assessed in our study and may point towards a different underlying mechanism. Previous studies have shown decreased miR-451a expression in peripheral blood lymphocytes from asthmatic children [34] and that in vivo overexpression of miR-451a inhibits airway remodelling by targeting cadherin-11 [35], suggesting a potential role in (controlling) airway remodelling.

One of the main strengths of this study is the size of the cohort and the careful characterisation of the subjects, which made it possible to make the association of miRNAs with clinical and inflammatory asthma features. Moreover, we had matched miRNA and mRNA data available from the same bronchial biopsies. Therefore, we could directly correlate global miRNA and gene expression levels to identify miRNAs involved in the regulation of gene expression changes in asthma. One limitation of this study is the lack of suitable replication datasets to validate the differentially expressed miRNAs in asthma and their association with gene expression. There are only a few studies that have determined miRNA expression in asthma patients; however, to the best of our knowledge, no miRNA–mRNA data are available of bronchial biopsies.

In conclusion, we have profiled miRNA expression in bronchial biopsies from asthma patients and controls, and identified several candidate “asthma miRNAs”, including miR-223-3p and miR-125b-5p, that were associated with multiple clinical and inflammatory asthma features and their negatively correlated genes are linked with the presence of asthma. Furthermore, miR-26a-5p is linked with inflammatory asthma features, especially eosinophilic inflammation. Therefore, these miRNAs are important candidates for future studies to unravel their mechanistic role in the pathogenesis of 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 material ERJ-01294-2021.Supplement

Shareable PDF

Supplementary Material

This one-page PDF can be shared freely online.

Shareable PDF ERJ-01294-2021.Shareable

Footnotes

  • Conflict of interest: M.P. Roffel has nothing to disclose.

  • Conflict of interest: I.M. Boudewijn has nothing to disclose.

  • Conflict of interest: J.L.L. van Nijnatten has nothing to disclose.

  • Conflict of interest: A. Faiz has nothing to disclose.

  • Conflict of interest: C.J. Vermeulen has nothing to disclose.

  • Conflict of interest: A.J. van Oosterhout owns GlaxoSmithKline shares.

  • Conflict of interest: K. Affleck is an employee of GlaxoSmithKline.

  • Conflict of interest: W. Timens reports personal fees for consultancy from Roche Diagnostics/Ventana, Bristol-Myers-Squibb and Diaceutics, personal fees for consultancy and lectures from MSD, outside the submitted work.

  • Conflict of interest: K.R. Bracke has nothing to disclose.

  • Conflict of interest: T. Maes reports grants from Ghent University (BOF19-GOA-008, BOF U4 funding; 05U40516) and Fund for Scientific Research in Flanders, during the conduct of the study; personal fees for advisory board work from GlaxoSmithKline, grants from Chiesi, outside the submitted work; and is shareholder of Oryzon Genomics and Mendelion Lifesciences SL.

  • Conflict of interest: I.H. Heijink has nothing to disclose.

  • Conflict of interest: C-A. Brandsma has nothing to disclose.

  • Conflict of interest: M. van den Berge reports research grants (Paid to institution) from GlaxoSmithKline, Sanofi, Novartis, Genentech and Roche, outside the submitted work.

  • Support statement: This work was supported through a scientific research collaborative agreement with GlaxoSmithKline and co-financed by the Ministry of Economic Affairs and Climate Policy by means of the PPP allowance made available by the Top Sector Life Sciences to stimulate public–private partnerships. M.P. Roffel was funded by U4 project (Ghent University Grant BOF 05U40516). Funding information for this article has been deposited with the Crossref Funder Registry.

  • Received May 5, 2021.
  • Accepted July 30, 2021.
  • Copyright ©The authors 2022. For reproduction rights and permissions contact permissions{at}ersnet.org
https://www.ersjournals.com/user-licence

References

  1. ↵
    1. Lambrecht BN,
    2. Hammad H
    . The airway epithelium in asthma. Nat Med 2012; 18: 684–692. doi:10.1038/nm.2737
    OpenUrlCrossRefPubMed
  2. ↵
    1. Moheimani F,
    2. Hsu ACY,
    3. Reid AT, et al.
    The genetic and epigenetic landscapes of the epithelium in asthma. Respir Res 2016; 17: 119. doi:10.1186/s12931-016-0434-4
    OpenUrlCrossRef
  3. ↵
    1. Selbach M,
    2. Schwanhäusser B,
    3. Thierfelder N, et al.
    Widespread changes in protein synthesis induced by microRNAs. Nature 2008; 455: 58–63. doi:10.1038/nature07228
    OpenUrlCrossRefPubMedWeb of Science
  4. ↵
    1. O'Connell RM,
    2. Rao DS,
    3. Chaudhuri AA, et al.
    Physiological and pathological roles for microRNAs in the immune system. Nat Rev Immunol 2010; 10: 111–122. doi:10.1038/nri2708
    OpenUrlCrossRefPubMedWeb of Science
  5. ↵
    1. Van Roosbroeck K,
    2. Fanini F,
    3. Setoyama T, et al.
    Combining anti-miR-155 with chemotherapy for the treatment of lung cancers. Clin Cancer Res 2017; 23: 2891–2904. doi:10.1158/1078-0432.CCR-16-1025
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Williams AE,
    2. Larner-Svensson H,
    3. Perry MM, et al.
    MicroRNA expression profiling in mild asthmatic human airways and effect of corticosteroid therapy. PLoS One 2009; 4: e5889. doi:10.1371/journal.pone.0005889
    OpenUrlCrossRefPubMed
    1. Jardim MJ,
    2. Dailey L,
    3. Silbajoris R, et al.
    Distinct microRNA expression in human airway cells of asthmatic donors identifies a novel asthma-associated gene. Am J Respir Cell Mol Biol 2012; 47: 536–542. doi:10.1165/rcmb.2011-0160OC
    OpenUrlCrossRefPubMed
  7. ↵
    1. Solberg OD,
    2. Ostrin EJ,
    3. Love MI, et al.
    Airway epithelial miRNA expression is altered in asthma. Am J Respir Crit Care Med 2012; 186: 965–974. doi:10.1164/rccm.201201-0027OC
    OpenUrlCrossRefPubMed
    1. Panganiban RP,
    2. Wang Y,
    3. Howrylak J, et al.
    Circulating microRNAs as biomarkers in patients with allergic rhinitis and asthma. J Allergy Clin Immunol 2016; 137: 1423–1432. doi:10.1016/j.jaci.2016.01.029
    OpenUrlCrossRefPubMed
  8. ↵
    1. Maes T,
    2. Cobos FA,
    3. Schleich F, et al.
    Asthma inflammatory phenotypes show differential microRNA expression in sputum. J Allergy Clin Immunol 2016; 137: 1433–1446. doi:10.1016/j.jaci.2016.02.018
    OpenUrlPubMed
  9. ↵
    1. Weidner J,
    2. Bartel S,
    3. Kılıç A, et al.
    Spotlight on microRNAs in allergy and asthma. Allergy 2021; 76: 1661–1678. doi:10.1111/all.14646
    OpenUrlPubMed
  10. ↵
    1. Broekema M,
    2. ten Hacken NHT,
    3. Volbeda F, et al.
    Airway epithelial changes in smokers but not in ex-smokers with asthma. Am J Respir Crit Care Med 2009; 180: 1170–1178. doi:10.1164/rccm.200906-0828OC
    OpenUrlCrossRefPubMedWeb of Science
  11. ↵
    1. Hoonhorst SJM,
    2. ten Hacken NHT,
    3. Lo Tam Loi AT, et al.
    Lower corticosteroid skin blanching response is associated with severe COPD. PLoS One 2014; 9: e91788. doi:10.1371/journal.pone.0091788
    OpenUrl
  12. ↵
    1. Vermeulen CJ,
    2. Xu C-J,
    3. Vonk JM, et al.
    Differential DNA methylation in bronchial biopsies between persistent asthma and asthma in remission. Eur Respir J 2019; 55: 1901280. doi:10.1183/13993003.01280-2019
    OpenUrl
  13. ↵
    1. Chou C-H,
    2. Shrestha S,
    3. Yang C-D, et al.
    miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res 2018; 46: D296–D302. doi:10.1093/nar/gkx1067
    OpenUrlCrossRefPubMed
  14. ↵
    1. Chen EY,
    2. Tan CM,
    3. Kou Y, et al.
    Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics 2013; 14: 128. doi:10.1186/1471-2105-14-128
    OpenUrlCrossRefPubMed
  15. ↵
    1. Kuleshov MV,
    2. Jones MR,
    3. Rouillard AD, et al.
    Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res 2016; 44: W90–W97. doi:10.1093/nar/gkw377
    OpenUrlCrossRefPubMed
  16. ↵
    1. Raudvere U,
    2. Kolberg L,
    3. Kuzmin I, et al.
    g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res 2019; 47: W191–W198. doi:10.1093/nar/gkz369
    OpenUrlCrossRefPubMed
  17. ↵
    1. Liu Z,
    2. Chen X,
    3. Wu Q, et al.
    miR-125b inhibits goblet cell differentiation in allergic airway inflammation by targeting SPDEF. Eur J Pharmacol 2016; 782: 14–20. doi:10.1016/j.ejphar.2016.04.044
    OpenUrlCrossRefPubMed
  18. ↵
    1. Van Pottelberge GR,
    2. Mestdagh P,
    3. Bracke KR, et al.
    MicroRNA expression in induced sputum of smokers and patients with chronic obstructive pulmonary disease. Am J Respir Crit Care Med 2011; 183: 898–906. doi:10.1164/rccm.201002-0304OC
    OpenUrlCrossRefPubMedWeb of Science
  19. ↵
    1. Gomez JL,
    2. Chen A,
    3. Diaz MP, et al.
    A network of sputum microRNAs is associated with neutrophilic airway inflammation in asthma. Am J Respir Crit Care Med 2020; 202: 51–64. doi:10.1164/rccm.201912-2360OC
    OpenUrl
  20. ↵
    1. Roffel MP,
    2. Bracke KR,
    3. Heijink IH, et al.
    miR-223: a key regulator in the innate immune response in asthma and COPD. Front Med 2020; 7: 196. doi:10.3389/fmed.2020.00196
    OpenUrl
  21. ↵
    1. Costello CM,
    2. Cahill E,
    3. Martin F, et al.
    Role of gremlin in the lung: development and disease. Am J Respir Cell Mol Biol 2010; 42: 517–523. doi:10.1165/rcmb.2009-0101TR
    OpenUrlCrossRefPubMedWeb of Science
  22. ↵
    1. Myllärniemi M,
    2. Lindholm P,
    3. Ryynänen MJ, et al.
    Gremlin-mediated decrease in bone morphogenetic protein signaling promotes pulmonary fibrosis. Am J Respir Crit Care Med 2008; 177: 321–329. doi:10.1164/rccm.200706-945OC
    OpenUrlCrossRefPubMedWeb of Science
  23. ↵
    1. Thomas B,
    2. Rutman A,
    3. Hirst RA, et al.
    Ciliary dysfunction and ultrastructural abnormalities are features of severe asthma. J Allergy Clin Immunol 2010; 126: 722–729. doi:10.1016/j.jaci.2010.05.046
    OpenUrlCrossRefPubMedWeb of Science
  24. ↵
    1. Yu X,
    2. Kong Y,
    3. Dore LC, et al.
    An erythroid chaperone that facilitates folding of α-globin subunits for hemoglobin synthesis. J Clin Invest 2007; 117: 1856–1865. doi:10.1172/JCI31664
    OpenUrlCrossRefPubMedWeb of Science
  25. ↵
    1. Newton DA,
    2. Rao KMK,
    3. Dluhy RA, et al.
    Hemoglobin is expressed by alveolar epithelial cells. J Biol Chem 2006; 281: 5668–5676. doi:10.1074/jbc.M509314200
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Hansbro PM,
    2. Kim RY,
    3. Starkey MR, et al.
    Mechanisms and treatments for severe, steroid-resistant allergic airway disease and asthma. Immunol Rev 2017; 278: 41–62. doi:10.1111/imr.12543
    OpenUrlCrossRefPubMed
  27. ↵
    1. Chalmers GW,
    2. MacLeod KJ,
    3. Thomson L, et al.
    Smoking and airway inflammation in patients with mild asthma. Chest 2001; 120: 1917–1922. doi:10.1378/chest.120.6.1917
    OpenUrlCrossRefPubMedWeb of Science
  28. ↵
    1. Levänen B,
    2. Bhakta NR,
    3. Torregrosa Paredes P, et al.
    Altered microRNA profiles in bronchoalveolar lavage fluid exosomes in asthmatic patients. J Allergy Clin Immunol 2013; 131: 894–903. doi:10.1016/j.jaci.2012.11.039
    OpenUrlCrossRefPubMed
  29. ↵
    1. Panganiban RPL,
    2. Pinkerton MH,
    3. Maru SY, et al.
    Differential microRNA epression in asthma and the role of miR-1248 in regulation of IL-5. Am J Clin Exp Immunol 2012; 1: 154–165.
    OpenUrlPubMed
  30. ↵
    1. Kwon Y,
    2. Kim Y,
    3. Eom S, et al.
    MicroRNA-26a/-26b-COX-2-MIP-2 loop regulates allergic inflammation and allergic inflammation-promoted enhanced tumorigenic and metastatic potential of cancer cells. J Biol Chem 2015; 290: 14245–14266. doi:10.1074/jbc.M115.645580
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Dileepan M,
    2. Rastle-Simpson S,
    3. Greenberg Y, et al.
    Effect of dual SEH/COX-2 inhibition on allergen-induced airway inflammation. Front Pharmacol 2019; 10: 1118. doi:10.3389/fphar.2019.01118
    OpenUrl
  32. ↵
    1. Wang T,
    2. Zhou Q,
    3. Shang Y
    . Downregulation of miRNA-451a promotes the differentiation of CD4+ T cells towards Th2 cells by upregulating ETS1 in childhood asthma. J Innate Immun 2021; 13: 38–48. doi:10.1159/000509714
    OpenUrl
  33. ↵
    1. Wang T,
    2. Zhou Q,
    3. Shang Y
    . MiRNA-451a inhibits airway remodeling by targeting Cadherin 11 in an allergic asthma model of neonatal mice. Int Immunopharmacol 2020; 83: 106440. doi:10.1016/j.intimp.2020.106440
    OpenUrlCrossRef
PreviousNext
Back to top
View this article with LENS
Vol 59 Issue 3 Table of Contents
European Respiratory Journal: 59 (3)
  • Table of Contents
  • Index by author
Email

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

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

Enter multiple addresses on separate lines or separate them with commas.
Identification of asthma-associated microRNAs in bronchial biopsies
(Your Name) has sent you a message from European Respiratory Society
(Your Name) thought you would like to see the European Respiratory Society web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
Citation Tools
Identification of asthma-associated microRNAs in bronchial biopsies
Mirjam P. Roffel, Ilse M. Boudewijn, Jos L.L. van Nijnatten, Alen Faiz, Corneel J. Vermeulen, Antoon J. van Oosterhout, Karen Affleck, Wim Timens, Ken R. Bracke, Tania Maes, Irene H. Heijink, Corry-Anke Brandsma, Maarten van den Berge
European Respiratory Journal Mar 2022, 59 (3) 2101294; DOI: 10.1183/13993003.01294-2021

Citation Manager Formats

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

Share
Identification of asthma-associated microRNAs in bronchial biopsies
Mirjam P. Roffel, Ilse M. Boudewijn, Jos L.L. van Nijnatten, Alen Faiz, Corneel J. Vermeulen, Antoon J. van Oosterhout, Karen Affleck, Wim Timens, Ken R. Bracke, Tania Maes, Irene H. Heijink, Corry-Anke Brandsma, Maarten van den Berge
European Respiratory Journal Mar 2022, 59 (3) 2101294; DOI: 10.1183/13993003.01294-2021
del.icio.us logo Digg logo Reddit logo Technorati logo Twitter logo CiteULike logo Connotea logo Facebook logo Google logo Mendeley logo
Full Text (PDF)

Jump To

  • Article
    • Abstract
    • Abstract
    • Introduction
    • Material and methods
    • Results
    • Discussion
    • Supplementary material
    • Shareable PDF
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF
  • Tweet Widget
  • Facebook Like
  • Google Plus One

More in this TOC Section

Original Research Articles

  • ELX/TEZ/IVA in people with CF and at least one F508del allele
  • Hypoxic burden to guide CPAP treatment in patients with OSA
  • Asthma and incident CHD: an observational and Mendelian Srandomisation study
Show more Original Research Articles

Asthma

  • MCID for impulse oscillometry in adults with asthma
  • Tool to detect small airways dysfunction in asthma clinical practice
  • Ecleralimab blocks responses to allergen in mild asthma
Show more Asthma

Related Articles

Navigate

  • Home
  • Current issue
  • Archive

About the ERJ

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

The European Respiratory Society

  • Society home
  • myERS
  • Privacy policy
  • Accessibility

ERS publications

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

Help

  • Feedback

For authors

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

For readers

  • Alerts
  • Subjects
  • Podcasts
  • RSS

Subscriptions

  • Accessing the ERS publications

Contact us

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

ISSN

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

Copyright © 2023 by the European Respiratory Society