# MiR-146a is over-expressed and controls IL-6 production in cystic fibrosis macrophages

Francesco R. Luly, Manuella Leveque, Valerio Licursi, Giuseppe Cimino, Corinne Martin-Chouly, Nathalie Théret, Rodolfo Negri, Luca Cavinato, Fiorentina Ascenzioni, Paola del Porto

► To cite this version:

Francesco R. Luly, Manuella Leveque, Valerio Licursi, Giuseppe Cimino, Corinne Martin-Chouly, et al.. MiR-146a is over-expressed and controls IL-6 production in cystic fibrosis macrophages. Scientific Reports, 2019, 9 (1), pp.16259. 10.1038/s41598-019-52770-w. hal-02391717

HAL Id: hal-02391717

https://univ-rennes.hal.science/hal-02391717

Submitted on 9 Jul 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC
REPORTS

nature research

OPEN

# MiR-146a is over-expressed and controls IL-6 production in cystic fibrosis macrophages

Francesco R. Luly<fnref n="1" />,⁶, Manuella Lévêque<fnref n="1" />,⁵,⁶, Valerio Licursi<fnref n="1" />, Giuseppe Cimino², Corinne Martin-Chouly³, Nathalie Théret³, Rodolfo Negri<fnref n="1" />,⁴, Luca Cavinato<fnref n="1" />, Fiorentina Ascenzioni<fnref n="1" /> & Paola Del Porto<fnref n="1" />*

Cystic fibrosis (CF) is an inherited disease that is characterised by susceptibility to bacterial infections and chronic lung inflammation. Recently, it was suggested that macrophages contribute to impaired host defence and excessive inflammatory responses in CF. Indeed, dysfunction attributed to CF macrophages includes decreased bacterial killing and exaggerated inflammatory responses. However, the mechanisms behind such defects have only been partially defined. MicroRNAs (miRNAs) have emerged as key regulators of several macrophage functions, including their activation, differentiation and polarisation. The goal of this study was to investigate whether miRNA dysregulation underlies the functional abnormalities of CF macrophages. MiRNA profiling of macrophages was performed, with 22 miRNAs identified as differentially expressed between CF and non-CF individuals. Among these, miR-146a was associated with significant enrichment of validated target genes involved in responses to microorganisms and inflammation. As miR-146a dysregulation has been reported in several human inflammatory diseases, we analysed the impact of increased miR-146a expression on inflammatory responses of CF macrophages. These data show that inhibition of miR-146a in lipopolysaccharide-stimulated CF macrophages results in increased interleukin-6 production, which suggests that miR-146a overexpression in CF is functional, to restrict inflammatory responses.

Cystic fibrosis (CF) is an autosomal genetic disease that is caused by mutations in the CF transmembrane conductance regulator (CFTR) gene<fnref n="1" />. CFTR is an epithelial channel that regulates anion transport and mucociliary clearance in the airways, and its absence results in mucus plugging and chronic bacterial infections, the persistence of which can lead to chronic inflammation and tissue destruction, and finally to respiratory failure²,³. In recent years, it has emerged that the lack of control of both infections and excessive inflammatory responses in CF patients is due, at least in part, to dysfunction of the innate immune system⁴. In particular, abnormalities in several functions of human macrophages in CF have been reported, including altered phagocytosis and decreased killing of lung pathogens, such as *Pseudomonas aeruginosa* and *Burkholderia cepacia*⁵⁻⁹. Moreover, CFTR deficiency in human macrophages has been associated with exaggerated inflammatory responses to bacterial stimuli and defective M2 polarisation¹⁰,¹¹.

MicroRNAs (miRNAs) are short non-coding RNA molecules that usually comprise 18 to 25 nucleotides that can inhibit expression of target genes by binding to complementary sequences in the 3'-untranslated regions of mRNAs. This targets these molecules to translational repression or degradation¹². miRNAs have emerged as key regulators of several aspects of macrophage biology, including their activation, differentiation and polarisation. Accordingly, M1 and M2 macrophages appear to be dominated by different miRNAs. MiR-155 is highly expressed in classical activated macrophages (M1), where it contributes to drive M1 polarisation by directly targeting interleukin (IL)-13 receptor α1, and it indirectly inhibits expression of M2-related genes¹³. In contrast, miR-511-3p and mir-223 have been shown to promote M2 macrophage polarisation. Mir-511-3p targets Rho-associated coiled-coil containing protein kinase 2 (Rock2), which is a serine threonine kinase that phosphorylates IRF4,

<footnote n="1">¹Department of Biology and Biotechnology, "C. Darwin" Sapienza University, Rome, Italy. ²Cystic Fibrosis Centre, Policlinico Umberto I, Rome, Italy. ³Inserm, EHESP, Irset (Institut de recherche en santé, environnement et travail) - UMR S 1085, F-35000 University of Rennes, Rennes, France. ⁴Institute of Molecular Biology and Pathology, Italian National Research Council, Rome, Italy. ⁵Present address: Service de Dermatologie CHU de Poitiers, University of Poitiers, Poitiers, France. ⁶These authors contributed equally: Francesco R. Luly and Manuella Lévêque. *email: paola.delporto@uniroma1.it</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** MiRNA profiling in CF and non-CF macrophages. (a) Heatmap showing the differentially expressed miRNAs. Each row indicates the expression level of the indicated single miRNA (z-score scaled), and each column is from one single subject. The row dendrogram represents the miRNA clustering using the hierarchical clustering method, with the Pearson correlations as the similarity measure. (b) Enriched categories of the miR-146a target genes. Bar plot showing the top significantly enriched categories of the GO biological processes of the miR-146a target genes. X-axis, significance of the categories; bars are ordered according to their −log₁₀(q-value); dashed line, −log₁₀(0.05) value.</figcaption></figure>

and thus promotes the expression of M2-related genes¹⁴. MiR-223 targets the transcription factor Pknox1, which suppresses the pro-inflammatory activation of macrophages¹⁵.

In addition, following macrophage activation, miRNAs contribute to prompt induction of inflammation and its rapid resolution through a complex network of interactions¹⁶. Accordingly, in a murine model of acute inflammation, it was shown that an NF-κB-miRNAs network controls the robustness of the inflammatory response and its time-dependent restriction. In particular, it has been proposed that the temporal asymmetry of miR-146a and miR-155 production and activity creates the correct balance between positive and negative regulators of NF-κB, which in turn regulates inflammatory responses. Indeed, it has been reported that miR-155 is rapidly and highly transcribed in response to an inflammatory stimulus, and that through repression of SHIP1 and SOCS1, miR-155 amplifies NF-κB activity. As an inflammatory response develops, miR-146a levels accumulate, which leads to repression of IRAK1 and TRAF6, and to subsequent attenuation of NF-κB activation. Thus, it appears that not only the miRNA levels, but also their temporal induction are important in the correct regulation of inflammatory responses.

To determine whether miRNA dysregulation indeed underlies the functional abnormalities of CF macrophages, we performed miRNA profiling in macrophages from CF and non-CF individuals. This led to the identification of a panel of differentially expressed miRNAs in CF macrophages, as compared to non-CF macrophages. Among these, miR-146a was associated with significant enrichment of validated target genes that are involved in responses to microorganisms and inflammation. In addition, as miR-146a dysregulation has been reported for several human inflammatory diseases, we analysed the impact of increased miR-146a expression on the inflammatory response of CF macrophages¹⁷⁻¹⁹. These data show that inhibition of miR-146a in lipopolysaccharide (LPS)-stimulated CF macrophages results in increased IL-6 production, which suggests that the function of miR-146a overexpression in CF is to restrict the inflammatory response.

## Results

**MiRNA profiling in CF macrophages and in-silico functional analysis.** MiRNA-specific RNA sequencing was performed in monocyte-derived macrophages from six patients with CF and from six age- and sex-matched non-CF individuals (controls). In this analysis, 785 miRNAs were identified and quantified in these macrophages from CF patients and controls. The most differentially expressed miRNAs were selected according to the following criteria: threshold of ≥1.5-fold or ≤−1.5-fold change, and false-discovery-rate-adjusted *P* (q) ≤ 0.05. From this analysis, 22 miRNAs that were differentially expressed in CF vs non-CF macrophages (10 up-regulated, 12 down-regulated) were selected (Fig. 1a).

To predict potential targets and to allow functional analysis of the differentially expressed miRNAs, *in-silico* analysis was performed, which included miRNA-target interactions, Gene Ontology (GO) enrichment, and network interactions. Thus, each of the 22 selected miRNAs was subjected to analysis of their miRNA targets, using a manually curated database of experimentally validated microRNA-target interactions: miRTarBase²⁰. Next, for the validated targets, enrichment analysis of the GO terms in biological processes and the Reactome pathways was performed²¹,²².

Through this approach, a significant number of validated targets was identified from miRTarBase for a subset of the 22 differentially expressed miRNAs, and significantly enriched categories (FDR < 0.1) were identified for five of the up-regulated miRNAs and six of the down-regulated miRNAs (Supplementary Figs S1 and S2). Collectively, each miRNA defined a unique set of categories that comprised many different terms and pathways. Among the differentially expressed miRNAs, miR-146a was the only one that showed enrichment in GO terms
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Up-regulation of miR-146a and down-regulation of TRAF6 mRNA in CF macrophages. (a) MiR-146a levels were measured using RT-qPCR in non-CF (n = 16) and CF (n = 11) macrophages. The data were normalised to expression of the RNU6B endogenous control. Each symbol represents a single individual.

**p** < 0.01 (Mann-Whitney tests) (b,c) TRAF6 and IRAK1 mRNAs were evaluated using RT-qPCR in non-CF (n = 8) and CF (n = 8) macrophages. The data were normalised to expression of the endogenous β-actin control.

Data are means ± SEM. *p* < 0.05 (unpaired t-tests).</figcaption></figure>

related to cellular responses to pathogens (Fig. 1b; Supplementary Fig. S1A). In accordance, the Reactome analysis confirmed the involvement of miR-146a in the Toll-like receptor (TLR) cascade, and in interferon and interleukin signalling (Supplementary Fig. S1B).

**MiR-146a is up-regulated in CF macrophages.** Among the differentially expressed miRNAs in the CF macrophages, miR-146a was the miRNA mostly implicated in the regulation of inflammation in innate immune cells, as determined by the GO terms and the Reactome enrichment analyses. Thus, considering the hyper-inflammatory phenotype of the CF macrophages, miR-146a was selected for further investigation²³,²⁴.

To validate the sequencing data, the analysis of miR-146a expression was extended to macrophages from 11 CF and 16 non-CF individuals, using RT-qPCR. These data are reported in Fig. 2a, and they show median miR-146a expression for CF macrophages as 3.11-fold that for non-CF macrophages (p = 0.0056), which confirms its significant up-regulation in these CF macrophages.

As negative regulators with short mRNA-binding sequences, miRNAs have the potential to regulate numerous target genes. To investigate whether miR-146a up-regulation had any effects on the expression of its targets in CF macrophages, the expression was examined for two of the most relevant validated miR-146a targets, *IRAK1* and *TRAF6*, which are crucial adaptors for TLR-mediated NF-κB signalling²⁵. Determination of TRAF6 and IRAK1 mRNA levels in these CF and non-CF macrophages using RT-qPCR revealed that the CF macrophages had significantly lower levels of TRAF6 mRNA (p = 0.0103) than the non-CF macrophages, but similar levels of IRAK1 mRNA (Fig. 2b,c).

Thus, these data demonstrated that increased miR-146a in CF macrophages is associated with mRNA decrease of at least one of its validated targets, TRAF6.

**MiR-146a inhibition in macrophages.** MiR-146a is known to promote resolution of inflammation in response to LPS¹⁶. To functionally assess the role of miR-146a in these CF macrophages, the response to LPS was examined for the macrophages following inhibition of endogenous miR-146a. To accomplish this, the macrophages were stimulated for 2 h with 100 ng/mL LPS, and then transfected with either miR-146a or a control inhibitor at 50 nM, as previously reported¹⁶. The knock-down efficiency was determined through evaluation of miR-146a expression in the non-CF and CF macrophages using RT-qPCR. As shown in Fig. 3, miR-146a was down-regulated to a similar degree in both the non-CF and CF macrophages. In particular, 16 h after transfection and compared to the control inhibitor, the miR-146a reduction was about 78% and 76% in these non-CF and CF macrophages, respectively (non-CF, *p* = 0.0007; CF, *p* = 0.0477; Fig. 3a,b). The expression of the unrelated miR-155 was not affected by the miR-146a inhibitor, which confirmed the specificity of this miR-146a inhibition (Fig. 3c). In addition, RT-qPCR and Western blotting of non-CF macrophages transfected with the miR-146a inhibitor demonstrated that miR-146a inhibition was associated with significant increases in *TRAF6* mRNA (*p* = 0.0395) and protein expression (*p* = 0.0273), compared to the control cells (Fig. 4a,b, Supplementary Fig. 3). At variance with this, analysis of IRAK1 protein expression did not show any significant differences between non-CF macrophages transfected with miR-146a or the control inhibitor (data not shown).

**MiR-146a reduces IL-6 production in CF macrophages.** Based on previous data in murine models that have shown that miR-146a modulates IL-6 and IL-10 production in response to LPS, the effects of miR-146a inhibition on the production of IL-6, IL-10 and IL-8 were compared in CF macrophages transfected with the miR-146a inhibitor¹⁶,²⁶. Here, macrophages isolated from seven CF patients and seven non-CF controls were stimulated with LPS and transfected with the miR-146a inhibitor, and the production of these cytokines was quantified in the supernatants using ELISA.
<figure><img src="image_4.png" /><figcaption>**Figure 5.** Increased production of IL-6 by miR-146a knock-down in CF macrophages in response to LPS. Cytokines produced by LPS-stimulated macrophages transfected with the miR-146a inhibitor (α-miR-146a-5p) or the control inhibitor (Ctr). IL-6 (top), IL-10 (middle) and IL-8 (bottom) were quantified in the supernatants of non-CF (n = 7) and CF (n = 7) macrophages 24 h after transfection, using ELISA. \*p < 0.05 (Wilcoxon tests).</figcaption></figure>

These data are reported in Fig. 5, and they demonstrated that inhibition of miR-146a in CF macrophages led to significant increases in production of IL-6, with respect to the same CF macrophages treated with the control inhibitor (p = 0.0156); however, this did not significantly affect the cytokine production in non-CF macrophages. In contrast, inhibition of miR-146a did not have any effects on the production of IL-10 and IL-8 by either the CF or the non-CF macrophages, which suggests that the anti-inflammatory effects of miR-146a on IL-6 production are specific (Fig. 5). Further support for the role of miR-146a in the modulation of IL-6 production in LPS-stimulated
<figure><img src="image_5.png" /><figcaption>**Figure 6.** Increased IL-6 mRNA levels in miR-146a knock-down CF macrophages in response to LPS. Relative quantification of IL-6 mRNA in LPS-stimulated non-CF (n = 5) or CF (n = 5) macrophages transfected with the miR-146a-5p inhibitor (α-miR-146a-5p) or the control inhibitor (Ctr). Data are means ± SEM; \*p < 0.05 (unpaired t-tests).</figcaption></figure>

CF macrophages was shown by the evidence that the inhibition of miR-146a leads to significant increases in IL-6 mRNA levels in the CF macrophages, as compared to the non-CF macrophages (Fig. 6, *p* = 0.045).

Taken together, these data show that in CF macrophages, miR-146a controls IL-6 production in response to bacterial stimuli.

Discussion

At present, there is general consensus on the involvement of immune cells in the progression of cystic fibrosis, and in particular of phagocytic cells<sup>27</sup>. This relatively new concept is supported by the demonstration that human monocytes and macrophages express CFTR and that phagocytosis and inflammatory responses to bacteria are altered in CF macrophages<sup>28-30</sup>. It is therefore possible that macrophage dysfunction contributes to disease progression and severity by impairing the ability of CF patients to clear lung infections and by sustaining inflammation. Although many studies have reported on the functional abnormalities of CF macrophages, the molecular mechanisms that underlie such defects have not been fully elucidated to date.

Evidence that has accumulated over the last decade has suggested that miRNAs are key regulators of multiple processes of macrophage biology, which includes their differentiation, infiltration and activation<sup>31</sup>. Thus, to determine whether dysregulation of miRNA expression occurs in CF macrophages, we performed global miRNA profiling of macrophages from non-CF and CF individuals, and used computational analysis to identify 22 differentially expressed miRNAs in the CF macrophages compared to the control macrophages. Similarly, in well-differentiated primary cultures of human CF airway epithelia cells and in CF bronchial brushing, altered expression of numerous microRNAs has been reported<sup>32,33</sup>. In these previous studies, miR-145, miR-223, miR-494 and miR-509-3p were shown to be over-expressed in CF cells compared to controls, and their expression was inversely correlated with CFTR levels, which suggested that these miRNAs can directly regulate CFTR expression. Furthermore, aberrant expression of miRNAs involved in the regulation of inflammatory pathways has been observed in CF epithelial cells. In CF bronchial brushings, down-regulation of miR-126 was correlated with significant up-regulation of Target of Myb1 (TOM1), a Tollip-binding protein that negatively regulates TLR2, TLR4 and IL-1RI signalling<sup>34</sup>. More recently, a parallel miRNome and transcriptome analyses of primary *in-vitro* differentiated epithelial cells from CF and non-CF individuals demonstrated that about half of the dysregulated miRNAs in CF samples was predicted to be associated with the NF-κB pathway. In particular, down-regulation of miR-199a-3p inversely correlated with increased expression of IKKβ and IL-8 in CF cells<sup>35</sup>. In the present study, our *in-silico* approaches identified common biological processes and molecular pathways that are associated with dysregulated miRNAs in CF macrophages, whereby we identified miR-146a as the miRNA mainly involved in inflammation.

MiR-146a dysregulation in innate immune cells has been implicated in several human inflammatory diseases, including rheumatoid arthritis, systemic lupus erythematosus and Sjogren's syndrome, which has suggested a role for miR-146a in CF macrophages<sup>17-19</sup>. Here, we have shown that expression of miR-146a was consistently and reproducibly increased in CF macrophages *versus* non-CF, and that it was associated with lower levels of its validated target TRAF6. Indeed, among the several miR-146a target genes in immune cells, which include *CXCR4*, *STAT1*, *IRAK1* and *TRAF6*, the last two of these have been the most frequently validated<sup>36,37</sup>. Dysregulation of miR-146a has been previously documented in CF. Indeed, Montanini *et al*. demonstrated that miR-146a was up-regulated in CF serum, and that its dysregulation was related to the onset of CF-related diabetes<sup>38</sup>. Moreover, miRNA profiling in peripheral blood mononuclear cells cultured with CF plasma identified miR-146a as one of the top-ranked up-regulated miRNAs<sup>39</sup>. A significant increase in miR-146a expression has also been reported in peripheral blood mononuclear cells and monocytes from rheumatoid arthritis patients, but no differences in the levels of TRAF6 and IRAK1 were shown between patients and controls, thus leading to the hypothesis that defective negative regulation of the target genes by miR-146a contributes to sustained pro-inflammatory cytokine production in patients<sup>18</sup>.
At variance with this, the evidence from the present study that inhibition of miR-146a in CF macrophages leads to significant increases in the production of IL-6 indicates that miR-146a is functional and contributes to the limiting of IL-6 production in these macrophages. It is interesting to note that miR-146a inhibition did not significantly affect the production of IL-6 in the control macrophages stimulated with LPS. These data are in agreement with previous findings that have shown that miR-146a inhibition does not significantly alter IL-6 production by the human monocytic THP1 cell line in response to LPS stimulation¹⁹.

There are several possibilities to explain why miR-146a specifically modulates IL-6 production in CF macrophages and does not appear to do the same in the cells from healthy donors. First, it might be a direct consequence of the increased expression of miR-146a in the CF macrophages, which through target-gene modulation can re-programme the cellular responses to endotoxins. Indeed, it has been shown that miR-146a ⁷/- bone-marrow-derived dendritic cells pre-treated with exosomes isolated from wild-type cells produce less IL-6 following LPS stimulation than cells that received exosomes from miR-146a ⁷/- mice²⁶. This suggests that even the modest increase in miR-146a as that induced by exosome transfer can change the immune responses of recipient cells.

Alternatively, it can be speculated that in CF macrophages, miR-146a up-regulation is functional to down-regulate the otherwise hyper-production of LPS-induced IL-6 that is due to an increase in TLR4 signalling in the CF macrophages²³. In CF macrophages, dysregulation of TLR signalling has been attributed to defective up-regulation of Caveolin-1, which allows trafficking of HO-1 to the caveolae in a p38 MAPK-dependent manner, and leads to down-regulation of TLR-4 pro-inflammatory signalling. Reduced Caveolin-1 expression was caused by high miR-199a-5p levels as a consequence of blunted phosphatidylinositol 3-kinase/ protein kinase B (AKT) pathway in CF macrophages⁴⁰.

Together with the evidence that miR-146a is up-regulated and contributes to the limiting of IL-6 production in CF macrophages, these data underline the importance of miRNAs in shaping inflammatory responses in CF macrophages. In particular, in response to LPS stimulation, IL-6 production by CF macrophages appears to be the result of simultaneous and opposing actions of pro-inflammatory and anti-inflammatory miRNAs, where their dysregulated expression leads to reduced or enhanced TLR signalling through modulation of receptor-associated molecules.

## Methods

**Study subjects.** Sixteen patients with CF (from the Regional Cystic Fibrosis Centre, Sapienza University, Rome Italy) that was confirmed by positive sweat tests and genotyping were enrolled in the study. Fourteen patients were F508del homozygous, and two were F508del/L558S and F508del/E409X. Blood samples from 18 age- and sex-matched healthy donors were used as controls.

MiRNA sequencing was performed for the macrophages from six CF individuals (three females, three males; five F508del/F508del, one F508del/L558S; mean age, 35.0 ± 14 years) and six sex-matched controls (mean age, 35.3 ± 11). All of the patients were clinically stable and blood samples were collected for isolation of CD14⁺ cells when they attended the clinic for routine evaluations. Informed written consent was obtained from all of the participants after approval of the study by the local Ethics Committee (Comitato Etico, Azienda Policlinico Umberto I, Rome, Italy; 1233/2016). All of the methods were carried out in accordance with the relevant guidelines and regulations.

**Isolation and differentiation of human monocytes.** Macrophages were generated as previously described⁴¹. In brief, CD14⁺ cells were purified from peripheral blood mononuclear cells by positive selection with an anti-CD14 monoclonal antibody coupled to magnetic beads (Miltenyi Biotec, Bergisch Gladbach, Germany). CD14⁺ cells were differentiated for 7 days in RPMI 1640 (Gibco-BRL, Invitrogren Corporation Carlsbad, CA, USA) supplemented with 20% foetal calf serum and 100 ng/mL recombinant macrophage-colony stimulating factor (PeproTech Inch, Rocky Hill, NY, USA).

**RNA extraction and real-time PCR.** The isolation of total RNA and the DNAse treatments were performed using miRNeasy kits and RNAse-free DNAse sets, respectively (Qiagen, Germany). For miRNA quantification, total RNA was reverse transcribed using miScript II RT kits (Qiagen, Germany), and real-time qPCR was performed with miScript SYBR Green PCR kits and miScript primer assays for miR-146a-5p (MS00003535; Qiagen, Germany) and miR-155 (MS00031486; Qiagen, Germany). For mRNAs, total RNAs were reverse transcribed with Reverse Transcription System kits (Promega, USA) and amplified using PowerUp SYBR Green Master Mix (Applied Biosystems, USA). The endogenous controls were RNU6B (MS00033740; Qiagen, Germany) and β-actin. The following primers were used: IRAK1-F (5'-TGAGGAACACGGTGTATG-3'), IRAK1-R (5'-GTTTGGGTGACGAAACCTGGA-3'), TRAF6-F (5'-TTGCCATGAAAAGATGCAGAGG-3'), TRAF6-R (5'-CGCTCAAACTATGAACAGCCTG-3'), IL-6-F (5'-AGGCACTGGCAGAAAACAAC-3'), IL-6-R (5'-TTTTCACCAGGCAAGTCTCC-3'), β-ACTIN-F (5'-GCCGGGACCTGACTGACTA-3'), β-ACTIN-R (5'-TGGTGATGACCTGGCCGT-3').

**MicroRNA-Seq.** MicroRNA-Seq library preparation and sequencing was performed by IGA Technology Services (Udine, Italy) using the Illumina TruSeq Small RNA Library Preparation kits (Illumina, San Diego, USA), according to manufacturer instructions. An average number of 15 million reads was obtained per sample. The reads quality was evaluated using FastQC (version 0.11.2; Babraham Institute, Cambridge, UK). After adapter cleaning, the resulting reads ranged from 94.7% to 97.7%. A very high percentage (98.4%-99.2%) of these reads was mapped to the human genome (hg38 build; Genome Reference Consortium Human GRCh38) using the *sRNAbench* command line tool with 'genome mapping mode' from the *sRNAtoolbox* suite of software⁴².
Briefly, the *sRNAbench* pipeline first aligns all of the reads to the reference genome using the *Bowtie* aligner, then the coordinates are compared to the miRBase annotation version 21<sup>43,44</sup>. A read is assigned to the reference RNA if its coordinates lie completely within the chromosome coordinates of the reference RNA. Subsequently, the differentially expressed microRNAs analysis was performed with the *sRNAde* pipeline from the *sRNAtoolbox* suite. The module generates an expression matrix and uses the R/Bioconductor package *edgeR* to infer differential expression<sup>45</sup>. By using *edgeR*, *sRNAbench* applies TMM normalisation for detection of differentially expressed microRNAs, which has been reported to be among the most stable of the methods available<sup>46</sup>. To understand the biological meaning of the differentially expressed microRNAs, we performed an enrichment analysis of experimentally validated microRNA-target interactions from the manually curated miRTarBase database<sup>20</sup>. Enrichment analysis was performed with hypergeometric tests from Bioconductor R package clusterProfiler using the GO Biological Process and the Reactome database categories<sup>21,22,47</sup>.

The miRNA-seq raw data files have been deposited in the NCBI Gene Expression Omnibus under GEO: GSE134702 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE134702).

**MiRNA inhibitor transfection.** Macrophages were cultured in RPMI 1640 with 20% foetal calf serum and 100 ng/mL macrophage-colony stimulating factor. On day seven, 3 × 10⁵ cells were stimulated with LPS (100 ng/ mL). After 2 h of stimulation, the cells were transfected with either miR-146a or the control inhibitor (Exiqon, Denmark), at a final concentration of 50 nM, using Lipofectamine RNAiMAX Transfection Reagent (13778; Invitrogen, USA), according to the manufacturer instructions. The supernatants were removed and stored at -80 °C for cytokine quantification, and the cells were collected for RNA extraction and quantitative PCR analysis.

**Immunoblotting.** Proteins were isolated using lysis buffer containing 50 mM HEPES, 150 mM NaCl, 1% Triton X100, 20 mM EDTA, and added protease inhibitor cocktail. Total protein levels were quantified using Bio-Rad protein assays, and equal amounts of protein were loaded and separated using 10% SDS-PAGE, followed by immunoblotting with the appropriate antibodies. These antibodies included: anti-TRAF6 (D21G3; Cell Signaling Technology, USA) diluted 1:1000, and anti-GAPDH (sc-25778; Santa Cruz, USA) diluted 1:400. The membranes were then washed, incubated with horseradish-peroxidase-linked secondary antibodies (GE Healthcare, UK), and visualised (Chemi Doc XRS system; Bio-Rad Laboratories Ltd, Hemel Hempstead, UK). Western blotting was quantified using the ImageJ software.

**ELISA.** The IL-6, IL-10 and IL-8 in cell supernatants were quantified using ELISA kits (R&D Systems, UK), following the manufacturer instructions. The sensitivity of these tests was 9.38 pg/mL, 31.3 pg/mL and 31.3 pg/ mL, respectively.

**Statistics.** All of the statistical analyses were carried out in the Graphpad Prism software using Student's t-tests, Mann-Whitney tests and Wilcoxon tests. *p* ≤ 0.05 was considered significant. In the Figures, statistically significant differences are indicated for *p* ≤ 0.05 (*), *p* ≤ 0.01 (**) and *p* ≤ 0.001 (**).

**Data availability**

The datasets generated and analysed during the current study are available from the corresponding author on
reasonable request.

Received: 5 August 2019; Accepted: 17 October 2019;

Published online: 07 November 2019

**References**

1. Elborn, J. S. Cystic fibrosis. *Lancet* **388**, 2519–2531 (2016).

2. Boucher, R. C. Airway surface dehydration in cystic fibrosis: pathogenesis and therapy. *Annu. Rev. Med.* **58**, 157–170 (2007).

3. Stoltz, D. A., Meyerholz, D. K. & Welsh, M. J. Origins of cystic fibrosis lung disease. *N. Engl. J. Med.* **372**, 1574–1575 (2015).

4. Nichols, D. P. & Chmiel, J. F. Inflammation and its genesis in cystic fibrosis. *Pediatr. Pulmonol.* **50** (Suppl 40), S39–56 (2015).

5. Del Porto, P. et al. Dysfunctional CFTR alters the bactericidal activity of human macrophages against Pseudomonas aeruginosa. *PLoS ONE* **6**, e19970 (2011).

6. Di Pietro, C. et al. Ezrin links CFTR to TLR4 signaling to orchestrate anti-bacterial immune response in macrophages. *Sci Rep* **7**, 10882 (2017).

7. Lévêque, M., Le Trionnaire, S., Del Porto, P. & Martin-Chouly, C. The impact of impaired macrophage functions in cystic fibrosis disease progression. *J. Cyst. Fibros.* **16**, 443–453 (2017).

8. Lévêque, M. et al. Phagocytosis depends on TRPV2-mediated calcium influx and requires TRPV2 in lipids rafts: alteration in macrophages from patients with cystic fibrosis. *Sci Rep* **8**, 4310 (2018).

9. Simonin-Le Jeune, K. et al. Impaired functions of macrophage from cystic fibrosis patients: CD11b, TLR-5 decrease and sCD14, inflammatory cytokines increase. *PLoS ONE* **8**, e75667 (2013).

10. Bruscia, E. M. & Bonfield, T. L. Cystic Fibrosis Lung Immunity: The Role of the Macrophage. *J Innate Immun* **8**, 550–563 (2016).

11. Tarique, A. A. et al. CFTR-dependent defect in alternatively-activated macrophages in cystic fibrosis. *J. Cyst. Fibros.* **16**, 475–482 (2017).

12. Iwakawa, H.-O. & Tomari, Y. The Functions of MicroRNAs: mRNA Decay and Translational Repression. *Trends Cell Biol.* **25**, 651–665 (2015).

13. Martinez-Nunez, R. T., Louafi, F. & Sanchez-Elsner, T. The interleukin 13 (IL-13) pathway in human macrophages is modulated by microRNA-155 via direct targeting of interleukin 13 receptor alpha1 (IL13Ralpha1). *J. Biol. Chem.* **286**, 1786–1794 (2011).

14. Squadrito, M. L. et al. miR-511-3p modulates genetic programs of tumor-associated macrophages. *Cell Rep* **1**, 141–154 (2012).

15. Zhuang, G. et al. A novel regulator of macrophage activation: miR-223 in obesity-associated adipose tissue inflammation. *Circulation* **125**, 2892–2903 (2012).

16. Mann, M. et al. An NF-κB-microRNA regulatory network tunes macrophage inflammatory responses. *Nat Commun* **8**, 851 (2017).

17. Tang, Y. et al. MicroRNA-146A contributes to abnormal activation of the type I interferon pathway in human lupus by targeting the key signaling proteins. *Arthritis Rheum.* **60**, 1065–1075 (2009).
18. Pauley, K. M. et al. Upregulated miR-146a expression in peripheral blood mononuclear cells from rheumatoid arthritis patients. *Arthritis Res. Ther.* **10**, R101 (2008).

19. Pauley, K. M. et al. Altered miR-146a expression in Sjögren's syndrome and its functional role in innate immunity. *Eur. J. Immunol.* **41**, 2029–2039 (2011).

20. Chou, C.-H. et al. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. *Nucleic Acids Res.* **46**, D296–D302 (2018).

21. Ashburner, M. et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. *Nat. Genet.* **25**, 25–29 (2000).

22. Croft, D. et al. The Reactome pathway knowledgebase. *Nucleic Acids Res.* **42**, D472–477 (2014).

23. Bruscia, E. M. et al. Abnormal trafficking and degradation of TLR4 underlie the elevated inflammatory response in cystic fibrosis. *J. Immunol.* **186**, 6990–6998 (2011).

24. Kopp, B. T. et al. Exaggerated inflammatory responses mediated by Burkholderia cenocepacia in human macrophages derived from Cystic fibrosis patients. *Biochem. Biophys. Res. Commun.* **424**, 221–227 (2012).

25. Taganov, K. D., Boldin, M. P., Chang, K.-J. & Baltimore, D. NF-kappaB-dependent induction of microRNA miR-146, an inhibitor targeted to signaling proteins of innate immune responses. *Proc. Natl. Acad. Sci. USA* **103**, 12481–12486 (2006).

26. Alexander, M. et al. Exosome-delivered microRNAs modulate the inflammatory response to endotoxin. *Nat Commun* **6**, 7321 (2015).

27. Munder, A. & Tümmler, B. Origins of cystic fibrosis lung disease. *N. Engl. J. Med.* **372**, 1574 (2015).

28. Barnaby, R. et al. Lumacaftor ( VX-809 ) restores the ability of CF macrophages to phagocytose and kill Pseudomonas aeruginosa. *Am. J. Physiol. Lung Cell Mol. Physiol.* **314**, L432–L438 (2018).

29. Van de Weert-van Leeuwen, P. B. et al. Optimal complement-mediated phagocytosis of Pseudomonas aeruginosa by monocytes is cystic fibrosis transmembrane conductance regulator-dependent. *Am. J. Respir. Cell Mol. Biol.* **49**, 463–470 (2013).

30. Zhang, S., Shrestha, C. L. & Kopp, B. T. Cystic fibrosis transmembrane conductance regulator (CFTR) modulators have differential effects on cystic fibrosis macrophage function. *Sci Rep* **8**, 17066 (2018).

31. Curtale, G., Rubino, M. & Locati, M. MicroRNAs as Molecular Switches in Macrophage Activation. *Front Immunol* **10**, 799 (2019).

32. Ramachandran, S. et al. Post-transcriptional regulation of cystic fibrosis transmembrane conductance regulator expression and function by microRNAs. *Am. J. Respir. Cell Mol. Biol.* **49**, 544–551 (2013).

33. Oglesby, I. K., Chotirmall, S. H., McElvaney, N. G. & Greene, C. M. Regulation of cystic fibrosis transmembrane conductance regulator by microRNA-145, -223, and -494 is altered in ΔF508 cystic fibrosis airway epithelium. *J. Immunol.* **190**, 3354–3362 (2013).

34. Oglesby, I. K. et al. miR-126 is downregulated in cystic fibrosis airway epithelial cells and regulates TOM1 expression. *J. Immunol.* **184**, 1702–1709 (2010).

35. Bardin, P. et al. Small RNA and transcriptome sequencing reveal the role of miR-199a-3p in inflammatory processes in cystic fibrosis airways. *J. Pathol.* **245**, 410–420 (2018).

36. Labbaye, C. et al. A three-step pathway comprising PLZF/miR-146a/CXCR4 controls megakaryopoiesis. *Nat. Cell Biol.* **10**, 788–801 (2008).

37. Wang, S. et al. MicroRNA-146a feedback suppresses T cell immune function by targeting Stat1 in patients with chronic hepatitis B. *J. Immunol.* **191**, 293–301 (2013).

38. Montanini, L. et al. miR-146a, miR-155, miR-370, and miR-708 Are CFTR-Dependent, Predicted FOXO1 Regulators and Change at Onset of CFRDs. *J. Clin. Endocrinol. Metab.* **101**, 4955–4963 (2016).

39. Zhang, X. et al. Cystic Fibrosis Plasma Blunts the Immune Response to Bacterial Infection. *Am. J. Respir. Cell Mol. Biol.* **61**, 301–311 (2019).

40. Zhang, P.-X. et al. Pharmacological modulation of the AKT/microRNA-199a-5p/CAV1 pathway ameliorates cystic fibrosis lung hyper-inflammation. *Nat Commun* **6**, 6221 (2015).

41. Cifani, N. et al. Reactive-oxygen-species-mediated P. aeruginosa killing is functional in human cystic fibrosis macrophages. *PLoS ONE* **8**, e71717 (2013).

42. Rueda, A. et al. sRNAtoolbox: an integrated collection of small RNA research tools. *Nucleic Acids Res.* **43**, W467–473 (2015).

43. Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. *Nat. Methods* **9**, 357–359 (2012).

44. Griffiths-Jones, S. The microRNA Registry. *Nucleic Acids Res.* **32**, D109–111 (2004).

45. Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. *Bioinformatics* **26**, 139–140 (2010).

46. Dillies, M.-A. et al. A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis. *Brief. Bioinformatics* **14**, 671–683 (2013).

47. Yu, G., Wang, L.-G., Han, Y. & He, Q.-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. *OMICS* **16**, 284–287 (2012).

**Acknowledgements**

This study was funded by grants from the non-profit organisation Vaincre la Mucoviscidose, Paris, France, and Sapienza University C26A14LFNN and C26A15SJ7N to PDP. ML received a Post-doctoral Fellowship from Vaincre la Mucoviscidose.

**Author contributions**

F.R.L., M.L. and L.C. performed and designed the experiments; V.L., R.N. and N.T. interpreted the data and wrote the manuscript; G.C. identified and provided clinical samples; P.D.P., F.A. and C.M.-C. designed experiments, interpreted data, obtained funding, and wrote the manuscript.

**Competing interests**

The authors declare no competing interests.

**Additional information**

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-019-52770-w.

Correspondence and requests for materials should be addressed to P.D.P.

Reprints and permissions information is available at www.nature.com/reprints.
**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre-
ative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons license and your intended use is not per-
mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the
copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019