key: cord-0066429-kmvane1b
authors: Bolourani, Siavash; Sari, Ezgi; Brenner, Max; Wang, Ping
title: Extracellular CIRP Induces an Inflammatory Phenotype in Pulmonary Fibroblasts via TLR4
date: 2021-07-23
journal: Front Immunol
DOI: 10.3389/fimmu.2021.721970
sha: e7a0029dc2ec09b18f7c42cddfc687bdef0b3ace
doc_id: 66429
cord_uid: kmvane1b

Extracellular cold-inducible RNA-binding protein (eCIRP), a new damage-associated molecular pattern (DAMP), has been recently shown to play a critical role in promoting the development of bleomycin-induced pulmonary fibrosis. Although fibroblast activation is a critical component of the fibrotic process, the direct effects of eCIRP on fibroblasts have never been examined. We studied eCIRP’s role in the induction of inflammatory phenotype in pulmonary fibroblasts and its connection to bleomycin-induced pulmonary fibrosis in mice. We found that eCIRP causes the induction of proinflammatory cytokines and differentially expression-related pathways in a TLR4-dependent manner in pulmonary fibroblasts. Our analysis further showed that the accessory pathways MD2 and Myd88 are involved in the induction of inflammatory phenotype. In order to study the connection of the enrichment of these pathways in priming the microenvironment for pulmonary fibrosis, we investigated the gene expression profile of lung tissues from mice subjected to bleomycin-induced pulmonary fibrosis collected at various time points. We found that at day 14, which corresponds to the inflammatory-to-fibrotic transition phase after bleomycin injection, TLR4, MD2, and Myd88 were induced, and the transcriptome was differentially enriched for genes in those pathways. Furthermore, we also found that inflammatory cytokines gene expressions were induced, and the cellular responses to these inflammatory cytokines were differentially enriched on day 14. Overall, our results show that eCIRP induces inflammatory phenotype in pulmonary fibroblasts in a TLR4 dependent manner. This study sheds light on the mechanism by which eCIRP induced inflammatory fibroblasts, contributing to pulmonary fibrosis.

Extracellular cold-inducible RNA-binding protein (eCIRP), a novel damage-associated molecular pattern (DAMP) discovered by our lab, plays a key role in sepsis, hemorrhagic shock, ischemiareperfusion injuries (1) (2) (3) (4) (5) . eCIRP induces these effects via its binding to and activation of the tolllike receptor 4 (TLR4) and myeloid differentiation factor 2 (MD-2) receptor complex (6) . The role of eCIRP in the activation of pulmonary cell populations such as macrophages (1) , alveolar epithelial cells (7, 8) , and endothelial cells (9) have been studied. However, the effect of eCIRP on pulmonary fibroblasts has yet to be studied.

Recently, we have discovered that eCIRP plays a role in the development of pulmonary fibrosis (publication forthcoming). Furthermore, a growing body of evidence suggests that blocking other DAMPs that also stimulate TLR4 can ameliorate pulmonary fibrosis (10) (11) (12) . These factors have been shown to induce an inflammatory phenotype in fibroblasts through a TLR4 dependent process (13, 14) and that halting TLR4 dependent pathway alleviates the process of inflammatory fibroblast transformation (15) (16) (17) . Fibroblasts are the key cells in the development of fibrosis in the lung. These mesenchymal cells are not terminally differentiated, and the principal regulator and potent inducer of fibroblast-to myofibroblast differentiation are transforming growth factor (TGF)-b (18) (19) (20) . However, for persistent activation of fibroblasts to occur, TGF-b signaling needs to be enhanced by induction of a primed cellular microenvironment (21) (22) (23) (24) . A variety of signals and receptors have been shown to aid fibroblasts proliferation and cellular differentiation (25) (26) (27) (28) (29) . However, the role of eCIRP on pulmonary fibroblasts has yet to be elucidated.

Inflammatory fibroblasts have been studied in various pathological processes such as neoplastic differentiation, autoimmune diseases, and fibrosis (30) (31) (32) . Multiple factors have been implicated in the induction of inflammatory fibroblasts, such as mechanical stress, immunoglobulins, cytokines, and DAMPs (33) (34) (35) (36) (37) (38) (39) (40) . A growing body of evidence suggests that many factors that induce inflammatory fibroblast transformation operate in a TLR4 dependent manner (13) (14) (15) (16) (17) . In this work, we examined the transcriptional response of pulmonary fibroblasts stimulated with eCIRP in vitro, focusing on proinflammatory, TLR4-induced pathways. We then evaluated how these pathways change over the course of bleomycin-induced pulmonary fibrosis.

Wild type (WT) C57BL/6 mice and TLR4 -/mice (Jackson Labs, on a C57BL/6 genetic background) were used in this study. All experiments involving live animals were carried out in accordance with the National Institutes of Health guidelines for the use of experimental animals (41) and were reviewed and approved by the Institutional Animal Care and Use Committee (IACUC) at the Feinstein Institutes for Medical Research. The lungs of 4-8-week-old WT and TLR4 -/mice were explanted, and each group of pulmonary fibroblasts was isolated from a pool of 3-4 lungs, as previously described (42) . The isolated fibroblasts were then cultured in DMEM supplemented with 10% heat-inactivated fetal bovine serum (FBS, MP Biomedicals), 2% penicillin-streptomycin, and L-glutamine (Gibco, ThermoFisher Scientific).

At passages 3-6, the cultured medium was replaced with reduced serum media Opti-MEM. After incubation overnight in Opti-MEM, the pulmonary fibroblasts were treated with phosphatebuffered saline (PBS), 1 mg/mL recombinant mouse (rm) CIRP, 2 ng/mL rmTGF-b1 (R&D Systems), or rmCIRP plus rmTGF-b1. The rmCIRP was produced in our lab, as previously described (1) . The cell lysates collected at 24 hours were used for high throughput mRNA sequencing and real-time reverse transcription-polymerase chain reaction (RT-PCR), the cell lysates collected at 48 hours were used for Western blotting (WB). Supernatants were collected at 48 hours and used for enzyme-linked immunosorbent assay (ELISA) quantification.

Isolation of the mRNA, RNA-Seq, and Analysis RNA was isolated using Ambion's mirVana kit (catalog #AM1561, Thermo Fisher) following the manufacturer's recommendations and quantified on a Nanodrop 2000 (catalog # ND-2000, Thermo Fisher). All samples had 260/280 and 260/230 ratios between 1.8 and 2.1. RNA quality was further assessed with an RNA Nano Chip using a Bioanalyzer 2100 (Agilent). All samples had a RIN number of 8 or higher. RNA was stored at -80 until Library preparation. Libraries were prepared using Illumina TruSeq Stranded Total RNA Library Prep Kit (catalog # RS-122-2203, Illumina Inc). The mRNA expression was obtained by RNA-Seq libraries run on NextSeq 550 following the manufacturer's instructions. FastQC was used to ensure there was no adapter contamination and reads had acceptable quality scores (43) . Sequenced segment alignments were performed using STAR2 (44) to the mouse GENCODE reference (45) . The mRNA sequencing raw counts were obtained using ht-seq counts (46) .

After high-throughput mRNA sequencing raw counts were obtained, the data were analyzed and normalized using the DESeq2 standard workflow (47) . M (log ratio) and A (mean average) (MA) plots were obtained using the function plotMA, to visually represent the differences between PBS and rmCIRPtreated cells) by transforming the data onto logarithmic of base 2 (log 2 ) fold changes vs. mean normalized count. The mRNA normalized counts were obtained using plotCounts function, which normalizes counts by sequencing depth. Principal component analysis was performed after variance stabilizing transformation (using vst function) to deal with the sampling variability of low counts and included correction for size factor (48) . The pathway and time series analysis was done using the Goexpress workflow (49) . The gene ontologies (GO) were first scored and ranked using the GO_analyse function, which consists of a random forest analysis to evaluate the ability of each gene to cluster samples according to the treatment groups in each analysis: pulmonary fibroblast mRNA sequencing data, PBS, eCIRP, TGF-b1, or eCIRP+TGF-b1; GSE132869 mRNA sequencing data, PBS, or bleomycin. The heatmap for each GO was created using the heatmap_GO function and Z-scores were scaled on each gene included in GO. Time-series data were analyzed using the expression_plot function, which plots the expression profile of a specific gene smoothed across time points (days 7, 14, 21, 28, and 42) while representing the mean and confidence interval of groups of lung tissue samples (bleomycin and PBS injected groups).

For Western blotting, the cultured WT and TLR4 -/primary pulmonary fibroblasts cells were scraped from 6 well culture plates and homogenized in a radioimmunoprecipitation assay buffer (RIPA) buffer containing phenylmethylsulfonyl fluoride (PMSF), Na-orthovanadate, and protease inhibitors. After centrifugation of cell lysates at 4°C and 12,000 rpm for 15 min, the total protein concentration was measured. The 10 mg protein of each sample was run in NuPAGE ™ 4-12% Bis-Tris Gel (Invitrogen, Thermo Fisher Scientific), and gels were transferred to a nitrocellulose membrane. Membranes were blocked with 0.1% casein 1h at the room temperature, after that, they were incubated with primary antibodies such as CIRP (catalog # 10209-2-AP, Protein Tech), b-actin (clone AC-15, catalog # A5441, MilliporeSigma), tumor necrosis factor (TNFa, catalog # 3707S, Cell Signalling Technology) and interleukin 1b (IL-1b, catalog #16806-1-AP, Protein Tech), MD2 (catalog # B100-56655, Novus Biologicals) and glyceraldehyde 3-phosphate dehydrogenase (GAPDH, catalog # 60004-1-lg, Proteintech) at 4°C overnight. After washing the membranes with Tris-Buffered Saline (TBS) containing 0.1% Tween-20 (TBST), and TBS; they were incubated with secondary antibodies for 1 h at room temperature. The signal of protein bands was detected and measured by using Odyssey CLx Imaging Machine and Image Studio Ver 5.2 Software. The bands were normalized with GAPDH. For ELISA, the levels of IL-6 of cell culture supernatant were measured using the mouse IL-16 DuoSet ELISA kit (DY206, R&D Systems), and it was performed according to the company's protocol.

Total RNA was isolated using Trizol reagent (Invitrogen) and reverse transcribed to cDNA using M-MLV (Moloney murine leukemia virus) reverse transcriptase (Applied Biosystems, Foster City, CA). Polymerase chain reactions (PCR) were carried out in a final volume of 10.5 ml, which included a 5 ml SYBR Green PCR master mix (Applied Biosystems) and 31 nM of both reverse and forward primers. GAPDH mRNA was used to normalize the amplification data, and fold changes were calculated in comparison with WT Control mice using the 2 -(DDCt) method. The sequence of forward and reverse primers used were as follows:

GAPDH: 5'-CATCACTGCCACCCAGAAGACTG-3' (forward) and 5'-ATGCCAGTGAGCTTCCCGTTCAG-3'(reverse); TNF-a: 5'-AGACCCTCACACTCAGATCATCTTC-3' (forward) and 5'-TTGCTACGACGTGGGCTACA-3' (reverse); IL-1b: 5'-CAGGATGAGGACATGAGCACC-3' (forward) and 5'-CTCTGCAGACTCAAACTCCAC-3'(reverse); IL-6: 5'-CCGGAGAGGAGACTTCACAG-3' (forward) and 5'-GGAAATTGGGGTAGGAAGGA-3'(reverse); MD2 (Ly96): 5'-CGCTGCTTTCTCCCATATTGA-3' (forward) and 5'-CCTCAGTCTTATGCAGGGTTCA-3'(reverse).

The expression datasets were was obtained directly from the gene expression omnibus website (https://www.ncbi.nlm.nih.gov/geo) using the getGEO function of R (http://www.r-project.org/) 4.0.3 library called GEOquery (50). After obtaining the dataset for GSE132869, the columns from lung tissues were selected. This GEO contains RNA sequenced from lung tissues of mice injected with bleomycin (10 mg/g/day, 100 ml of 1 mg/ml solution per injection of a~20 g mouse) or PBS five times per week for two weeks. The tissues were collected on 7, 14, 21, 28, and 42 days after the first injection (51) . While the exact timing of peak inflammation and fibrosis in the lung, in this model, is not clear, based on previous studies, inflammation occurs between days 7 and 14, and peak fibrosis likely happens between days 21 and 28 (52) (53) (54) . Day 14, which we call the inflammatory-to-fibrotic transition phase of exposure to bleomycin, is a time when both processes are active.

Data from mRNA sequencing normalized count were analyzed using the results function of DESeq2 workflow, which estimates the log 2 fold change of maximum likelihood estimates (MLE) between treated and untreated samples per unit of change on the specific genes, and p-values provided are Benjamini-Hochberg (BH)-adjusted p-values with false discovery cut-off set to 0.1 (55) . The p-values are set to 0 if they are too small to be represented (the smallest floating-point value to be represented in R is 2.225074e-308). Data for RT-PCR, WB, and ELISA are shown in the bar graphs as mean and standard deviation (SD), with individual values depicted by colored points. The comparisons were performed using two-way analysis of variance (ANOVA) followed by Tukey's multiple comparison test with the factors of two-way ANOVA being TLR4 gene expression status (WT vs. TLR4 -/-) and fibroblast treatment (PBS vs. CIRP). The groups were checked for Shapiro-Wilk test for normality prior to the application of two-way ANOVA.

To examine the effects of eCIRP on the expression of proinflammatory cytokines, we examined the normalized mRNA count (RNA-Seq), RT-PCR expression, and supernatant protein levels of TNF-a, IL-1b, and IL-6 in pulmonary fibroblasts isolated from the lungs of WT and TLR4 -/mice and cultured for 24h with rmCIRP or PBS. Compared with PBS, rmCIRP significantly increased the number of the mRNA normalized counts of TNF-a by 4.9 log 2 fold (adjusted p-value≈0, Figure 1A ), IL-1b by 9.9 log 2 fold (adjusted p-value≈0, Figure 1D ), and IL-6 by 2.1 log 2 fold (adjusted p-value = 2.31e-33, Figure 1G ). However, there were no significant increases in the mRNA normalized counts of TNF-a, IL-1b, and IL-6 in TLR4 -/cells stimulated with rmCIRP as compared to PBS (Figures 1A, D, G ; for more details, please refer to Supplementary Table 1) . We confirmed the RNA-Seq results using RT-PCR amplification of mRNA obtained from a separate set of cells. Compared with PBS, rmCIRP significantly induced the WT fibroblast RT-PCR expression of TNF-a (23.6 fold, p<0.0001, Figure 1B ), IL-1b (3368.7 fold, p= 0.0002, Figure 1E ), and IL-6 (14.0 fold, p<0.0001, Figure 1H ) expression by rmCIRP. The mRNA levels of TNF-a, IL-1b, and IL-6 measured by RT-PCR in TLR4 -/pulmonary fibroblasts stimulated with rmCIRP, however, were significantly lower than those of WT stimulated with rmCIRP ( Figures 1B, E, H) . We also measured the effects of stimulation with rmCIRP on the pulmonary fibroblast protein expression of the above-mentioned cytokines (TNF-a and IL-1b by WB, IL-6 by ELISA). Compared with PBS, rmCIRP significantly induced the WT fibroblast protein expression of TNF-a (1.7 fold, p=0.0021, Figure 1C ) and IL-6 (53.4 fold, p<0.0001, Figure 1I ) expression by rmCIRP. The protein levels of TNF-a and IL-6 in TLR4 -/pulmonary fibroblasts stimulated with rmCIRP, however, were significantly lower than those of WT stimulated with rmCIRP ( Figures 1C, I) . Stimulation with rmCIRP increased the protein expression of IL-1b in WT pulmonary fibroblasts by 2.5 fold, which showed a trend towards statistical significance (p=0.0884, Figure 1F ). However, the protein expression of IL-1b in rmCIRPstimulated TLR4 -/pulmonary fibroblasts was 75% lower than that of rmCIRP-stimulated WT cells (p=0.0191, Figure 1F ). Taken together, these results indicate that eCIRP induces proinflammatory cytokines in pulmonary fibroblasts in a TLR4 dependent manner.

eCIRP Differentially Enriches TLR4 and its Dependent Pathways in Pulmonary Fibroblasts eCIRP binds to and activates TLR4 (1). To determine the importance of TLR4 in eCIRP's induction of proinflammatory pulmonary fibroblasts, we examined the MA plots of genomic data from high throughput sequencing of mRNA data in WT and TLR4 -/pulmonary fibroblasts. Our analysis showed that the eCIRP-induced translational profile was restricted to WT pulmonary fibroblasts and did not occur in TLR4 -/pulmonary fibroblasts (Figures 2A, B) . In order to determine the effects of rmCIRP on the TLR4 and its related pathways, we analyzed the enrichment for transcripts associated with the TLR4 (GO:0034142) and Myd88 (GO:0002755) signaling pathways in WT pulmonary fibroblasts treated with PBS, rmCIRP, rmTGF-b1, or the combination of the rmCIRP and rmTGF-b1. Our analysis showed that stimulation with rmCIRP was associated with differential enrichment for both GO:0034142 and GO:0002755 (p-values; GO:0034142, <0.0001, Figure 2C ; GO:0002755, 0.0007, Figure 2D increased the MD2 mRNA normalized count in WT cells by 0.3 log 2 fold (adjusted p=0.004, Figure 3A ). This induction, however, was not observed in the mRNA normalized count of TLR4 -/cells ( Figure 3A) . Stimulation with rmCIRP also increased the mRNA expression of MD2 in WT pulmonary fibroblasts by 1.9 fold (p=0.047), but a similar increase was not seen in TLR4 -/pulmonary fibroblasts ( Figure 3B) . Although there was an increase in protein expression of MD-2 in WT pulmonary fibroblasts treated with eCIRP, the increase was not statistically significant. Conversely, the MD2 protein expression of TLR4 -/pulmonary fibroblasts stimulated with rmCIRP was significantly lower than that of WT pulmonary fibroblasts stimulated with rmCIRP (p=0.3741, Figure 3C ). These data suggest that eCIRP induces MD2 expression in pulmonary fibroblasts in a TLR4 dependent manner.

In the previous section, we showed that TGF-b1 did not significantly affect the transcriptional profile of TLR4 and Myd88 signaling pathways in WT pulmonary fibroblasts ( Figures 2C, D) .

To determine whether eCIRP requires TGF-b1 to induce the proinflammatory phenotype, we conducted a principal component (PC) analysis of mRNA profile sequenced from WT pulmonary fibroblasts treated with PBS, rmCIRP, rmTGF-b1, and the combination of the rmCIRP and rmTGF-b1. Our analysis showed that two PCs explain 95% of the variation in mRNA genomic data in WT pulmonary fibroblasts: PC1 explained 84% of the variance and PC2, 11%. Furthermore, the variance explained by PC1 aligned well with stimulation by rmCIRP, whereas PC2 aligned well with stimulation by TGF-b1 ( Figure 4A ). We then examined selected genes to see the direction of contribution and observed that TNF-a, IL-1a, IL-1b, IL-6, and CxCl2 contributed to the variance towards the eCIRP-aligned PC1. Furthermore, MD2 (Ly96) and Myd88 also contributed to the variance along PC1, albeit to a lesser degree ( Figure 4A) . To see the full list of genes examined, please see Supplementary Figure 1 . Taken together, these results indicate that eCIRP induces a proinflammatory phenotype of pulmonary fibroblasts independently of TGF-b1.

Stimulation with rmCIRP also caused a transcriptional reprogramming that was differentially enriched for genes associated with the TNF-a mediated signaling pathway (GO:0033209, p-value = 0.008, Figure 5A ) and positive regulation of IL-6 production (GO:00032755, p-value <0.0001, Figure 5B ) in WT pulmonary fibroblasts stimulated with rmCIRP or the combination of the rmCIRP and rmTGF-b1. However, differential enrichment for these two pathways was not present in the transcriptome of pulmonary fibroblasts cultured with rmTGF-b1. To see high-resolution figures to examine genes involved in GOs analyzed, please see Supplementary Figures 2, 3 . As such, eCIRP differentially enriches the proinflammatory transcriptional profile in pulmonary fibroblasts that includes TNF-a signaling and IL-6 production.

To compare the transcriptomic changes after stimulation with eCIRP with those in the lungs of mice exposed to bleomycin, we examined a publicly available Gene Expression Omnibus (GEO) profile (GSE132869) containing high throughput RNA sequences from lung tissue of WT mice collected on days 7, 14, 21, 28, or 42 from the start of bleomycin injection. Similar to our observations in eCIRP-stimulated lung fibroblasts, the lung transcriptome on day 14 after bleomycin showed significant upregulation of TLR4, MD2, and Myd88 mRNA counts ( Figures 6A, C, E) . The lung transcriptome on day 14 was also differentially enriched for genes associated with the TLR signaling pathway (GO:0002224, p-value = 0.005, Figure 6B ), TRIF dependent TLR signaling pathway (GO:0035666, p-value = 0.013, Figure 6D ), and Myd88 dependent TLR signaling pathway (GO:0002755, p-value = 0.04, Figure 6F ). To see a higher resolution of Figure 6B , please refer to Supplementary Figure 4 . As a whole, the results indicate that the TLR4 and related genes and pathways activated by eCIRP in pulmonary fibroblasts are also transcriptionally active on day 14 after bleomycin injection.

To determine the timeline at which inflammatory cytokines and their cellular response pathways are induced, we examined the mRNA of cytokines TNF-a, IL-6, and the cellular response to these cytokines. Similar to our observations in eCIRP-stimulated lung fibroblasts, the lung transcriptome on day 14 after bleomycin showed significant upregulation of TNF-a and IL-6 mRNA counts (Figures 7A, C) . The lung transcriptome on day 14 was also differentially enriched for genes associated with the cellular response to TNF-a (GO:0071356, p-value = 0.009, Figure 7B ) and cellular response to IL-6 (GO:0071354, p-value= 0.021, Figure 7D ). To see a higher resolution of Figure 7B , please refer to Supplementary Figure 5 . Taken together, the results indicate that the TLR4 and related genes and pathways activated by eCIRP in pulmonary fibroblasts are also transcriptionally active on day 14 after bleomycin injection. 

Fibroblasts are key cells involved in the fibrotic process (56) . In order to illuminate the possible contribution of eCIRP to the fibrotic process, we investigated how eCIRP influences the transcriptional profile of the pulmonary fibroblasts. We evaluated the transcriptional response of lung fibroblasts stimulated with eCIRP with a particular focus on proinflammatory and TLR4-induced pathways. We also evaluated how these pathways change over the course of bleomycin-induced pulmonary fibrosis. What we observed was that eCIRP amplifies proinflammatory cytokines in a TLR4 dependent manner and differentially enriches the inflammatory pathways. Furthermore, we found that this activation pattern recapitulates that of the inflammatory-tofibrotic transition phase in the bleomycin-induced pulmonary fibrosis model. We have previously shown that CIRP binds with high affinity to the TLR4-MD2 complex, as well as to TLR4 and MD2 individually (1). Furthermore, TLR4 expression is known to be elevated in lung biopsies of patients with pulmonary fibrosis caused by chronic inflammatory diseases (21, 57, 58) . While TGF-b1 is required for profibrotic activation of the fibroblasts, this activation and induction of organ fibrosis require a primed cellular microenvironment, which can be induced by TLR4 activation (22, 59, 60) . The role of inflammatory fibroblasts in the pathogenesis of fibrosis has been studied (36, 61, 62) . Inflammatory fibroblasts persistently activate the fibroblasts in an autocrine/paracrine manner in the lungs, resulting in pulmonary fibrosis (63) (64) (65) . The proinflammatory role of TLR4-MD2 complex and Myd88 signaling pathway in the lungs has also been evaluated (66) (67) (68) (69) and associated with the fibrotic process leading to pulmonary fibrosis (60, (70) (71) (72) (73) (74) . In our study, we found that eCIRP induces proinflammatory cytokines and their pathways in pulmonary fibroblasts in a TLR4 dependent manner. Furthermore, we showed that the TLR coreceptor MD2 and the downstream Myd88 pathway are enriched by eCIRP in a TLR4 dependent manner. These results shed light on the mechanism by which eCIRP induces proinflammatory pulmonary fibroblasts during the inflammatory phase in the pathogenesis of bleomycin-induced pulmonary fibrosis.

It is well known that TGF-b1 is the master regulator, and its signaling is required for the promotion and maintenance of the fibrotic process in the lungs (75) (76) (77) . The recent works by Bhattacharyya et al. showed that TLR4 signaling synergizes with TGF-b in the activation of fibroblasts in the setting of pulmonary fibrosis (21, 78, 79) . After illustrating that eCIRP proinflammatory programming of fibroblasts is done in a TLR4 dependent, we set out to see the effect of TGF-b in this process. Our principal component analysis showed that the principal component that explains the variation in inflammatory cytokines and chemokines aligns with the treatment of eCIRP and not TGF-b1 (Figure 4 and Supplementary Figure 6 ). This confirms that the proinflammatory programming of fibroblasts is dependent on eCIRP, not TGF-b1.

In order to establish how eCIRP's induction of proinflammatory phenotype may contribute to the pathogenesis of bleomycin- induced pulmonary fibrosis, we examined the pathways induced by eCIRP in fibroblast in the lung tissues of mice. We showed that TLR4, MD2, Myd88, and related pathways are significantly induced in lungs collected on day 14 of bleomycin-induced pulmonary fibrosis. Proinflammatory cytokines and the related were also significantly overrepresented in day 14 of pulmonary fibrosis lung tissues, showing that the changes in the lung transcriptome during the inflammatory-to-fibrotic transition phase after bleomycin injection recapitulate those pulmonary fibroblasts stimulated with eCIRP. This is in line with previous studies showing an increase in the lung infiltration by total white blood cells, neutrophils, and macrophages peaks at day 14, the expressions of profibrotic markers peak at day 21-28 (53, 54, 80, 81) . We have previously shown that eCIRP promotes acute lung injury via activation of macrophages, neutrophils, pneumocytes, and lung vascular endothelial cells in the contexts of sepsis, hemorrhagic shock, and intestinal ischemia/ reperfusion injury (1, 5, 6, (82) (83) (84) . However, little was known about eCIRP's role in chronic inflammatory diseases of the lung. Our results show that eCIRP activates a proinflammatory phenotype in pulmonary fibroblasts, illuminating a potential mechanism by which eCIRP in the lung environment in the setting of chronic inflammation plays a critical priming role leading to lung fibrosis.

Our results should be interpreted taking into consideration some particularities and limitations. Our RNASeq analysis focused on TLR4 and proinflammatory genes, not on other high-ranking but unrelated genes and pathways. Our study was underpowered to detect changes of small magnitude and high variance. We did not perform a time-course analysis in our in vitro experiments, which might further show how fibroblast transcription can fluctuate and change over time, as other sequencing analyses have demonstrated (85) (86) (87) . In our in vivo time series analysis, while we were able to measure total CIRP levels in pulmonary tissues and show it is elevated prior to day 14 ( Supplementary Figure 7) , we were not able to measure the levels of eCIRP in the lung interstitial space, which would further clarify the link between eCIRP and pulmonary fibrosis. This model only considered the subcutaneous bleomycin-induced pulmonary fibrosis, which models global inflammatory diseases that cause interstitial lung fibrosis (52, 60, 88, 89) . The study of other models can shed light on the role of DAMPs such as eCIRP on other pulmonary diseases causing lung fibrosis.

In summary, we showed that eCIRP induces an inflammatory phenotype in pulmonary fibroblasts, which recapitulates that of lung tissue during the inflammatory phase of bleomycin-induced pulmonary fibrosis. Our results not only shed light on the mechanism by which eCIRP promotes inflammatory pulmonary fibroblasts but also show a possible stage at which this phenotype contributes to the most commonly used preclinical model of pulmonary fibrosis. Our study suggests that eCIRP is critical for the proinflammatory fibroblasts reprogramming and may represent a druggable target to stop the progression of pulmonary fibrosis in chronic inflammatory diseases affecting the lung.

The data has been uploaded to NCBI -accession number GSE178255 (https://www.ncbi.nlm.nih.gov/gds/?term= GSE178255[Accession]).

The animal study was reviewed and approved by Feinstein Institute's IACUC.

SB, MB, and PW designed the experiments. SB and ES performed all experiments and analyzed the data. MB and PW critically reviewed the manuscript. MB and PW conceived the idea and supervised the project. All authors contributed to the article and approved the submitted version.

This study was supported by the National Institutes of Health (NIH) grants R01HL076179 and R35GM118337 (PW).

Cold-Inducible RNA-binding Protein (CIRP) Triggers Inflammatory Responses in Hemorrhagic Shock and Sepsis

A Cold-Inducible RNA-binding Protein (CIRP)-Derived Peptide Attenuates Inflammation and Organ Injury in Septic Mice

Deficiency of Cold-Inducible Ribonucleic Acid-Binding Protein Reduces Renal Injury After Ischemia-Reperfusion

Extracellular Cold Inducible RNA-binding Protein Mediates Binge Alcohol-Induced Brain Hypoactivity and Impaired Cognition in Mice

Deficiency in Cold-Inducible RNA-binding Protein Attenuates Acute Respiratory Distress Syndrome Induced by Intestinal Ischemia-Reperfusion

Extracellular CIRP (eCIRP) and Inflammation

Extracellular CIRP Induces Inflammation in Alveolar Type II Cells Via TREM-1

Extracellular microRNA 130b-3p Inhibits eCIRP-induced Inflammation

Cold-Inducible RNAbinding Protein Causes Endothelial Dysfunction Via Activation of Nlrp3 Inflammasome

The Role of High Mobility Group Box1 in Pulmonary Fibrosis

The Small Heat-Shock Protein ab-Crystallin Is Essential for the Nuclear Localization of Smad4: Impact on Pulmonary Fibrosis

Clinical Significance of Serum S100 Calcium-Binding Protein A4 in Idiopathic Pulmonary Fibrosis

Tetraand Penta-Acylated Lipid A Structures of Porphyromonas Gingivalis LPS Differentially Activate TLR4-mediated Nf-kb Signal Transduction Cascade and Immuno-Inflammatory Response in Human Gingival Fibroblasts

Salivary Pellets Induce a Pro-Inflammatory Response Involving the TLR4-NF-kB Pathway in Gingival Fibroblasts

Sirt1 Alleviates LPS Induced Inflammation of Periodontal Ligament Fibroblasts Via Downregulation of TLR4

Tormentic Acid Inhibits LPS-induced Inflammatory Response in Human Gingival Fibroblasts Via Inhibition of TLR4-mediated Nf-kb and MAPK Signalling Pathway

Antifibroblast Antibodies in Systemic Sclerosis Induce Fibroblasts to Produce Profibrotic Chemokines, With Partial Exploitation of Toll-Like Receptor 4

Tgf-b: The Master Regulator of Fibrosis

Effects of Pirfenidone on Transforming Growth Factor-Beta Gene Expression at the Transcriptional Level in Bleomycin Hamster Model of Lung Fibrosis

Regulation of TGF-b Family Signaling by Inhibitory Smads

Toll-Like Receptor 4 Signaling Augments Transforming Growth Factor-b Responses: A Novel Mechanism for Maintaining and Amplifying Fibrosis in Scleroderma

TLR4 Enhances TGF-beta Signaling and Hepatic Fibrosis

Signaling Pathway Cooperation in TGF-b-Induced Epithelial-Mesenchymal Transition

Tgf-b Pathway Activation by Idiopathic Pulmonary Fibrosis (IPF) Fibroblast Derived Soluble Factors is Mediated by IL-6 Trans-Signaling

Tnf-a-Induced NF-kb Activation Promotes Myofibroblast Differentiation of LR-MSCs and Exacerbates Bleomycin-Induced Pulmonary Fibrosis

Novel High-Throughput Myofibroblast Assays Identify Agonists With Therapeutic Potential in Pulmonary Fibrosis That Act Via EP2 and EP4 Receptors

Recent Developments in the Pathobiology of Lung Myofibroblasts

Expression Analysis of microRNAs and mRNAs in Myofibroblast Differentiation of Lung Resident Mesenchymal Stem Cells

Inflammatory Fibroblasts in Cancer

Functionally Distinct Disease-Associated Fibroblast Subsets in Rheumatoid Arthritis

Fibroblast-Extracellular Matrix Interactions in Tissue Fibrosis

Porphyromonas Gingivalis Lipopolysaccharide Induces a Pro-Inflammatory Human Gingival Fibroblast Phenotype

Ifn-b Plays Both Pro-and Anti-Inflammatory Roles in the Rat Cardiac Fibroblast Through Differential STAT Protein Activation

Pulsed Electromagnetic Fields Decrease Proinflammatory Cytokine Secretion (IL-1b and TNF-a) on Human Fibroblast-Like Cell Culture

Autoantibodies to Fibroblasts Induce a Proadhesive and Proinflammatory Fibroblast Phenotype in Patients With Systemic Sclerosis

The B1-agonist [des-Arg10]-kallidin Activates Transcription Factor NF-kappaB and Induces Homologous Upregulation of the Bradykinin B1-receptor in Cultured Human Lung Fibroblasts

Proinflammatory Cytokine Profiles in Pulp Fibroblasts Stimulated With Lipopolysaccharide and Methyl Mercaptan

Regulation of Cytokine Production in Human Gingival Fibroblasts Following Treatment With Nicotine and Lipopolysaccharide

IL-1 Induces Proinflammatory Leukocyte Infiltration and Regulates Fibroblast Phenotype in the Infarcted Myocardium

The Care and Feeding of an IACUC: The Organization and Management of an Institutional Animal Care and Use Committee

Establishing Primary Adult Fibroblast Cultures From Rodents

Babraham Bioinformatics -Fastqc A Quality Control Tool for High Throughput Sequence Data

Optimizing RNA-Seq Mapping With STAR

GENCODE Reference Annotation for the Human and Mouse Genomes

Htseq-a Python Framework to Work With High-Throughput Sequencing Data

Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data With Deseq2

A Variance-Stabilizing Transformation for Gene-Expression Microarray Data

Goexpress: An R/Bioconductor Package for the Identification and Visualisation of Robust Gene Ontology Signatures Through Supervised Learning of Gene Expression Data

Geoquery: A Bridge Between the Gene Expression Omnibus (GEO) and Bioconductor

Systemic Bleomycin Model of Scleroderma

Modeling Pulmonary Fibrosis Through Bleomycin Delivered by Osmotic Minipump: A New Histomorphometric Method of Evaluation

Bleomycin Induces Drug Efflux in Lungs. A Pitfall for Pharmacological Studies of Pulmonary Fibrosis

Mouse Models of Bleomycin-Induced Pulmonary Fibrosis

Discovering the False Discovery Rate

Fibroblasts in Fibrosis: Novel Roles and Mediators

An Essential Role for Fibronectin Extra Type III Domain A in Pulmonary Fibrosis

Toll-Like Receptor 4 Activation Attenuates Profibrotic Response in Control Lung Fibroblasts But Not in Fibroblasts From Patients With IPF

Transcriptional Repression of the Transforming Growth Factor b (Tgf-b) Pseudoreceptor BMP and Activin Membrane-Bound Inhibitor (BAMBI) by Nuclear Factor kb (Nf-kb) p50 Enhances TGF-b Signaling in Hepatic Stellate Cells

TLR4-Dependent Fibroblast Activation Drives Persistent Organ Fibrosis in Skin and Lung

Aicar Treatment Reduces Interstitial Fibrosis in Aging Mice: Suppression of the Inflammatory Fibroblast

Exuberant Fibroblast Activity Compromises Lung Function Via ADAMTS4

Tnf-a-Stimulated Fibroblasts Secrete Lumican to Promote Fibrocyte Differentiation

Fibroblast Paracrine TNF-a Signaling Elevates Integrin A5 Expression in Idiopathic Pulmonary Fibrosis (IPF)

HMGB1 Accelerates Alveolar Epithelial Repair Via an IL-1b-and avb6 Integrin-Dependent Activation of TGF-b1

The TLR4-MyD88 Signaling Axis Regulates Lung Monocyte Differentiation Pathways in Response to

TLR5 Participates in the TLR4 Receptor Complex and Promotes MyD88-Dependent Signaling in Environmental Lung Injury

TLR4 Signalling in Pulmonary Stromal Cells Is Critical for Inflammation and Immunity in the Airways

Discovery of New MD2 Inhibitor From Chalcone Derivatives With Anti-Inflammatory Effects in LPSinduced

Targeting the TLR4-MD2 Axis in Systemic Sclerosis

Ligustilide Modulates Oxidative Stress, Apoptosis, and Immunity to Avoid Pathological Damages in Bleomycin Induced Pulmonary Fibrosis Rats Via Inactivating TLR4/MyD88/NF-KB P65

Modified Citrus Pectin Ameliorates Myocardial Fibrosis and Inflammation Via Suppressing Galectin-3 and TLR4/Myd88/NF-kb Signaling Pathway

Bixin Attenuates Carbon Tetrachloride Induced Oxidative Stress, Inflammation and Fibrosis in Kidney by Regulating the Nrf2/TLR4/MyD88 and PPAR-g/Tgf-b1/Smad3 Pathway

Pathogenesis of Systemic Sclerosis: Recent Insights of Molecular and Cellular Mechanisms and Therapeutic Opportunities

Tgf-b Activation and Lung Fibrosis

The Impact of TGF-b on Lung Fibrosis: From Targeting to Biomarkers

Developmental Pathways in the Pathogenesis of Lung Fibrosis

FibronectinEDA Promotes Chronic Cutaneous Fibrosis Through Toll-like Receptor Signaling

Tenascin-C Drives Persistence of Organ Fibrosis

Time Course of Bleomycin-Induced Lung Fibrosis

Preparedness of Personal Protective Equipment and Implementation of New CPR Strategies for Patients With Out-of-Hospital Cardiac Arrest in the COVID-19 Era

Cold-Inducible RNA-binding Protein-Derived Peptide C23 Attenuates Inflammation and Tissue Injury in a Murine Model of Intestinal Ischemia-Reperfusion

CIRP Increases ICAM-1 Phenotype of Neutrophils Exhibiting Elevated iNOS and NETs in Sepsis

Cold-Inducible RNA-binding Protein Mediates Neuroinflammation in Cerebral Ischemia

Time-Series Single-Cell RNA-Seq Data Reveal Auxin Fluctuation During Endocycle

RNA-Seq: Technical Variability and Sampling

RNA-Seq Differential Expression Studies: More Sequence or More Replication?

Bleomycin Delivery by Osmotic Minipump: Similarity to Human Scleroderma Interstitial Lung Disease

Repetitive Intradermal Bleomycin Injections Evoke T-helper Cell 2 Cytokine-Driven Pulmonary

We would like to thank Monowar Aziz, PhD for valuable discussions and assistance.

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.721970/ full#supplementary-material Conflict of Interest: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.Publisher's Note: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.Copyright © 2021 Bolourani, Sari, Brenner and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.