ADAR1 Isoforms Regulate Let-7d Processing in Idiopathic Pulmonary Fibrosis

Double-stranded RNA adenosine deaminase 1 (ADAR1) is significantly down-regulated in fibroblasts derived from Idiopathic Pulmonary Fibrosis (IPF) patients, and its overexpression restored levels of miRNA-21, PELI1, and SPRY2. There are two ADAR1 isoforms in humans, ADAR1-p110 and ADAR1-p150, generated by an alternative promoter. Let-7d is considered an essential microRNA in Pulmonary Fibrosis (PF). In silico analysis revealed COL3A1 and SMAD2, proteins involved in the development of IPF, as Let-7d targets. We analyzed the role of ADAR1-p110 and ADAR1-p150 isoforms in the regulation of Let-7d maturation and the effect of this regulation on the expression of COL3A1 and SMAD2 in IPF fibroblast. We demonstrated that differential expression and subcellular distribution of ADAR1 isoforms in fibroblasts contribute to the up-regulation of pri-miR-Let-7d and down-regulation of mature Let-7d. Induction of overexpression of ADAR1 reestablishes the expression of pri-miR-Let-7d and Let-7d in lung fibroblasts. The reduction of mature Let-7d upregulates the expression of COL3A1 and SMAD2. Thus, ADAR1 isoforms and Let-7d could have a synergistic role in IPF, which is a promising explanation in the mechanisms of fibrosis development, and the regulation of both molecules could be used as a therapeutic approach in IPF.


Introduction
Idiopathic Pulmonary Fibrosis (IPF) is the most frequent and aggressive interstitial lung disease. It is considered a consequence of several risk factors such as smoking habits, genetic determinants, environmental factors, gastroesophageal reflux, viral infections, and aging; however, its origin and further progression are not fully understood [1][2][3][4][5]. In situ, fibroblasts differentiate into myofibroblasts and form fibroblastic foci, which are responsible for secreting excessive amounts of extracellular matrix (ECM) components such as fibrillar 2 of 13 collagens [6][7][8]. The excessive accumulation of these proteins entails the destruction of the pulmonary architecture [1][2][3][4][5].
The molecular mechanisms of IPF have not been fully described. Some evidence suggests the possible participation of some microRNAs [9][10][11][12][13]. Several editing microRNAs by Adenosine Deaminase acting on RNA (ADAR) are affected in their maturation process mediated by Drosha and Dicer, or they are redirecting to new mRNAs target because their seed region nucleotides are modified [14][15][16][17][18][19][20][21][22][23]. This edition represents a specific tissue regulatory mechanism of microRNA expression [22]. ADAR1 has two isoforms generated by the alternative promoter, ADAR1-p110 and ADAR1-p150. ADAR1-p110 is the amino-terminally truncated isoform located in nuclear regions and constitutively expressed. ADAR1-p150 is the full-length isoform, which is interferon-inducible and is preferentially located in cytoplasmic regions [14,21]. ADAR1 isoforms have 2-3 double-stranded RNA binding domains (dsRNA) in the amino-terminal region, followed by a deaminase domain in the carboxyl-terminal region and defined domains of subcellular localization that allow them to shift between the nucleolus, nucleus, and cytoplasm [19]. Sub-cellular location differences suggest that ADAR1 isoforms display differential RNA regulation [3]. ADAR1 is significantly down-regulated in IPF fibroblasts, and the over-expression of ADAR1 and ADAR2 reestablishes the expression levels of miRNA-21, PELI1, and SPRY2 in fibroblasts from IPF patients [14]. Different works showed a Let-7d down-regulation in lung fibroblasts of IPF patients [9,10,12]. SMAD2 and 3, positive regulator of the TGF-β1 pathway, represents one of the main Let-7d targets, while Let-7d low expression promotes an increase in mesenchymal markers such as HMGA2, vimentin, N-cadherin, and ACTA2 [10,13]. In leukemia stem cells, ADAR1 can edit pri-miRNA-Let-7d transcript forming a secondary double-stranded structure that induces a negative regulation of Let-7d and self-renewal [23].
Dysregulation of the microRNA maturation processing caused by ADAR1 editing, and its effects on pulmonary fibroblasts, are not well known. We studied the role of the ADAR1-p110 and ADAR1-p150 isoforms in the regulation of Let-7d maturation and the effect of these modifications on the expression of collagen-3, a putative Let-7d target related to IPF.

Expression of ADAR1 Isoforms in Controls and IPF Fibroblasts
ADAR1-p110 isoform in IPF fibroblasts (n = 10) was downregulated (** p = 0.008, Figure 1A), whereas ADAR1-p150 isoform was overexpressed (*** p = 0.00054; Figure 1B), compared with control fibroblasts (n = 10). Additionally, we verified that there were no differences in the basal levels of ADAR1 expression between both control and fibrotic fibroblasts ( Figure S1). We further confirmed these observations at the protein level by Western blot analysis, and we observed that ADAR1-p110 isoform levels were less than 66% (* p = 0.05; Figure 1C,D) in total protein extracts in IPF fibroblasts compared to controls. The ADAR1-p150 isoform did not show changes ( Figure 1C,E). Different distribution of fluorescence intensity of ADAR1 was detected in the nuclei of the IPF fibroblasts (n = 6) compared to control cells (n = 6) ( Figure 2A). ADAR1 was in the nucleolar region mainly in control fibroblasts, whereas in IPF cells, ADAR1 was exclusively detected in the nucleoplasm Figure 2A,B). These preliminary observations suggested an altered subnuclear location of ADAR1 in IPF. It was confirmed by the analysis of fluorescence intensity. The ratio of ADAR1 nucleolar to nuclear expression is significantly higher in control fibroblasts than in IPF fibroblasts ( Figure 2C).

ADAR1 Is Hyperactive in IPF Fibroblasts
RESS-qPCR analysis [24] was performed to determine RNA editing mediated by ADAR1 in lung fibroblasts. We analyzed six controls and six IPF fibroblasts; IPF fibroblasts possess increased editing activity in specific adenosine of APOBEC3D compared with control fibroblasts (* p = 0.0193; Figure 3). The editing activity on APOBEC3D in control and fibrotic fibroblasts transfected with ADAR1 isoforms showed that the catalytic activity is increased after transfection with active plasmids. In contrast, the mutant plasmids do not modify the activity. In control fibroblasts (black bars), ADAR1-p110 wt shows a significant increase in editing activity (** p = 0.003737). A remarkable increase in editing activity is observed with the ADAR1-p150 isoform (** p = 0.000139), while the most significant increase the editing activity in IPF fibroblasts are observed after overexpression of the ADAR1-p110 (** p = 0.007030). An increase in ADAR1-p150 wt activity was also significant (** p = 0.009158). Editing-deficient ADAR-p110 mt and ADAR-p150 mt constructs did not affect the editing activity. ADAR editing activity on healthy control and fibrotic fibroblasts considerably decreases upon Pentostatin treatment, which is a specific inhibitor of the editing activity of Adenosine Deaminases (ADA) ( Figure 4). Pentostatin, as a specific ADA inhibitor, decreases the catalytic activity in control fibroblasts (** p = 0.0059) and IPF fibroblasts (** p = 0.0014). Over-expression of ADAR1-p150 wt (** p = 0.0058) and ADAR1-p110 wt (** p = 0.0018) increased the catalytic activity of ADAR in control fibroblasts. On the contrary, the overexpression of ADAR1-p150 (p = 0.8) and ADAR1-p110 mutant (p = 0.72) did not induce significant changes in the catalytic activity in the control fibroblasts. However, ADAR1-p150 wt over-expression in IPF fibroblasts increases its activity (** p = 0.0042), similar to the overexpression of ADAR1-p110 (** p = 0.0067). Over-expression of ADAR1 mutant plasmids does not modify the activity. ADAR1-p150 mutant (p = 0.2) and ADAR1-p110 mutant (p = 0.02).

Basal Expression of Let7d and pri-miR-Let7d in Normal and IPF Fibroblasts and Correlation with ADAR1p110 and p150
The basal expression of pri-miR-Let7d is higher in IPF fibroblasts than in normal fibroblasts (p = 0.05, Figure 5A). On the other hand, Let7d expression is lower in IPF fibroblasts than in Normal fibroblasts (p = 0.0028, Figure 5D). The correlation of ADAR1 isoforms does not show significant differences between ADAR1-p110 ( Figure 5B) and ADAR1-p150 ( Figure 5C). Let7d displayed a negative correlation with ADAR1-p150 (p < 0.0001) ( Figure 5E) and a positive correlation with ADAR1-p110 (p = 0.0014) ( Figure 5F).

Figure 3.
The catalytic activity of ADAR1 on APOBEC3D in control fibroblasts compared with fibrotic fibroblasts. ADAR1 catalytic activity in basal conditions in IPF fibroblasts concerning control fibroblasts (** p = 0.0018). Editing activity increase is shown as the ratio of the edition divided by the ratio of wild-type (n = 6). The catalytic activity is increased after transfection with the catalytically active plasmids, whereas the mutant plasmids do not modify the activity. In control fibroblasts (black bars),  . The catalytic activity of ADAR1 on APOBEC3D in control fibroblasts compared with fibrotic fibroblasts. ADAR1 catalytic activity in basal conditions in IPF fibroblasts concerning control fibroblasts (** p = 0.0018). Editing activity increase is shown as the ratio of the edition divided by the ratio of wild-type (n = 6). The catalytic activity is increased after transfection with the catalytically active plasmids, whereas the mutant plasmids do not modify the activity. In control fibroblasts (black bars), ADAR1 p110 wt shows a significant increase in editing activity (** p = 0.003737). Remarkable increase in editing activity is observed with the ADAR1-p150 isoform (** p = 0.000139), while the greatest increase the editing activity in IPF fibroblasts are observed after over-expression of the ADAR1-p110 (** p = 0.007030). An increase in ADAR1-p150 wt activity was also significant (** p = 0.009158).

Basal Expression of Let7d and pri-miR-Let7d in Normal and IPF Fibroblasts and Correlation
with ADAR1p110 and p150.

Expression of COL3A1 and SMAD2 Targets of Let-7d in IPF
The expression of COL3A1 and SMAD2 mRNA targets of Let-7d potentially involved with the progression of IPF was analyzed in lung fibroblasts of control and IPF samples. The Let-7d target, COL3A1 obtained by in silico analysis, was upregulated (** p = 0.006, Figure  7A) in IPF when Let-7d expression was decreased. Over-expression of ADAR1p110 wt isoform showed a negative effect on the expression levels of COL3A1 (p < 0.05) and SMAD2 (p < 0.05) ( Figure 7B,C) in both control and IPF fibroblasts, however, ADARp150 wt showed an increased effect in the expression of COL3A1 in normal and IPF fibroblasts (p < 0.0001) ( Figure 7B). However, expression levels of SMAD2 did not show changes by ADARp150 ( Figure 7C). Figure 6. The expression of pri-miR-Let-7d is modified by the over-expression of ADAR1-p110 wt and ADAR1-p150 wt. In (A), the over-expression of ADAR1-p110 wt in control fibroblasts significantly decreases the expression of pri-miR-Let-7d (* p = 0.054479, black bars). In IPF fibroblasts, both isoforms promote the decrease in pri-miR-Let-7d. ADAR1-p110 wt (** p = 0.0002108, gray bars) and ADAR1-p150 wt (** p = 0.00038554) in a higher proportion compared to the control fibroblasts. The expression of Let-7d in (B) was increased in the fibroblasts transfected with ADAR1. In control fibroblasts transfected with ADAR1, the expression of Let-7d is increased significantly ADAR1-p110 wt (** p = 0.001477, black bars) and ADAR1-p150 wt (** p = 0.001946, black bars). In IPF fibroblasts, the over-expression of catalytically active isoforms promotes a significant increase in Let-7d. ADAR1-p110 wt (*** p = 0.000151, gray bars) and ADAR1-p150 wt (*** p = 0.000855, gray bars).

Expression of COL3A1 and SMAD2 Targets of Let-7d in IPF
The expression of COL3A1 and SMAD2 mRNA targets of Let-7d potentially involved with the progression of IPF was analyzed in lung fibroblasts of control and IPF samples. The Let-7d target, COL3A1 obtained by in silico analysis, was upregulated (** p = 0.006, Figure 7A) in IPF when Let-7d expression was decreased. Over-expression of ADAR1p110 wt isoform showed a negative effect on the expression levels of COL3A1 (p < 0.05) and SMAD2 (p < 0.05) ( Figure 7B,C) in both control and IPF fibroblasts, however, ADARp150 wt showed an increased effect in the expression of COL3A1 in normal and IPF fibroblasts (p < 0.0001) ( Figure 7B). However, expression levels of SMAD2 did not show changes by ADARp150 ( Figure 7C).  Figure 6. The expression of pri-miR-Let-7d is modified by the over-expression of ADAR1-p110 wt and ADAR1-p150 wt. In (A), the over-expression of ADAR1-p110 wt in control fibroblasts significantly decreases the expression of pri-miR-Let-7d (* p = 0.054479, black bars). In IPF fibroblasts, both isoforms promote the decrease in pri-miR-Let-7d. ADAR1-p110 wt (** p = 0.0002108, gray bars) and ADAR1-p150 wt (** p = 0.00038554) in a higher proportion compared to the control fibroblasts. The expression of Let-7d in (B) was increased in the fibroblasts transfected with ADAR1. In control fibroblasts transfected with ADAR1, the expression of Let-7d is increased significantly ADAR1-p110 wt (** p = 0.001477, black bars) and ADAR1-p150 wt (** p = 0.001946, black bars). In IPF fibroblasts, the over-expression of catalytically active isoforms promotes a significant increase in Let-7d. ADAR1-p110 wt (*** p = 0.000151, gray bars) and ADAR1-p150 wt (*** p = 0.000855, gray bars).

Expression of COL3A1 and SMAD2 Targets of Let-7d in IPF
The expression of COL3A1 and SMAD2 mRNA targets of Let-7d potentially involved with the progression of IPF was analyzed in lung fibroblasts of control and IPF samples. The Let-7d target, COL3A1 obtained by in silico analysis, was upregulated (** p = 0.006, Figure  7A) in IPF when Let-7d expression was decreased. Over-expression of ADAR1p110 wt isoform showed a negative effect on the expression levels of COL3A1 (p < 0.05) and SMAD2 (p < 0.05) ( Figure 7B,C) in both control and IPF fibroblasts, however, ADARp150 wt showed an increased effect in the expression of COL3A1 in normal and IPF fibroblasts (p < 0.0001) ( Figure 7B). However, expression levels of SMAD2 did not show changes by ADARp150 ( Figure 7C).

Discussion
The ADAR1-p110 isoform is down-regulated both transcript and protein level, while the ADAR1-p150 isoform is overexpressed in IPF fibroblasts, even though at protein ADAR1-p150 level, there are no differences. This discrepancy may be due in part to the fact that the messenger RNA encoding ADAR1p150 has IRES-type elements in its sequence, which in pathological or stress conditions could respond to the presence of IRES transactivating factors, such as the polypyrimidine track binding protein (PTBP1). This results in the use of different translation starts, producing different isoforms of ADAR1 including ADARp110 [25]. Additionally, the distribution of these proteins is different in control fibroblasts; ADAR1-p110 is found mainly in the nucleolus, while in fibrotic cells, it is found in the nucleoplasm (Figure 2A-C). Previous reports showed that nucleolar ADAR2 location decreases its editing activity, and we hypothesize that the same regulation would be present in ADAR1 since both ADARs have dsRNA-binding motifs (DRBMs) and binding to rRNA [26]. Then, we propose that the ADAR1-p110 isoform is not correctly performing the editing in controls. In the case of fibrotic fibroblasts, this isoform is in the nucleoplasm, a region in which it has been demonstrated that the catalytic activity of this protein is executed [20], and even though the fibrotic fibroblasts have a diminished expression of the ADAR1-p110 isoform, their editing activity is increased. A report has mentioned that the ADAR editing activity is independent of their protein levels and expression levels [21], which supports our hypothesis that the ADAR1-editing activity depends mainly on its subcellular location. The expression levels of the ADAR1 isoforms and their location indicate that the regulatory mechanisms of each isoform are specific and complex. Differential locations for each ADAR1 isoform make us think that ADAR1-p110 is responsible for editing pri-miR-Let-7d and that cytoplasmic ADAR1-p150 for editing Let-7d transcript.
When the basal expression of an ADAR1 isoform is low, as in the case of ADAR1-p110, and is exogenously overexpressed, its edit activity levels do increase. Hence, we could suggest the existence of a negative feedback mechanism, which is essential to maintain specific expression levels in ADAR1 proteins to achieve an optimal function of the cells. However, we observed that the edition activity significantly increased in both control and fibrotic primary cultures when transfected with the wild-type isoforms, suggesting the existence of other regulatory mechanisms involved in the catalytic activity of ADAR in IPF fibroblasts. Analysis in silico showed that pri-miR-Let-7d has more than 70% adenosines edited by ADAR1; it is considered hyper-edition. Hyper-edited primary transcript is accumulated and rapidly degraded by Tudor-SN and is a less efficient substrate for Drosha and Dicer [26][27][28]. Thus, ADAR1 may inhibit the processing of the pri-miR-Let-7d through the loss of its double-stranded secondary structure. Consequently, Dicer and Drosha do not recognize it, redirecting it to another target or its degradation [29][30][31][32]. This degradation or deficiency in the processing of the pri-miR-Let-7d induced by the hyper-edition of the same transcript could explain the decrease in the Let-7d mature transcript in IPF since a down-regulation was observed along with ADAR1 expression.
The COL3A1 and SMAD2 mRNAs, putative targets of Let-7d, were quantified in IPF and controls. High levels of COL3A1 and SMAD2 were corroborated in IPF. A decrease in Let-7d expression induces an increase in COL3A1 and SMAD2 expression, which are involved in the development of pulmonary fibrosis.
When ADAR1 isoforms are over-expressed, the expression levels of COL3A1 and SMAD2 mRNA are modified, ADARp-110 induces downregulation of COL3A1 and SMAD2, and interestingly ADARp-150 induced an increase in COL3A1 on both normal and IPF fibroblasts. In this sense, we propose that the ADAR1-p150 isoform is editing pri-miR-Let-7d, avoiding its maturation to Let-7d. However, this has no effect on SMAD2.
To demonstrate the direct relationship between ADAR1 activity and the expression of Let-7d, we performed the experiment treating fibrotic fibroblasts with pentostatin, an inhibitor of ADAR activity ( Figure S2). We were able to observe that Let-7d levels increased and pri-miR-Let-7d levels decreased, which demonstrates that ADAR1 activity directly influences the inhibition of maturation of pri-miR-Let-7d. These observations correlate with what is shown in Figure 4, where we observe a decrease in ADAR activity, which is reestablished with the transfection of the ADAR1 isoforms.
ADAR1-p110 and ADAR1-p150 regulate the processing of the pri-miR-Let-7d and mature Let-7d in such a way that the sub-expression of the latter deregulated, in turn, the expression of COL3A1 and SMAD2, proteins involved in the development of idiopathic pulmonary fibrosis. In summary, the regulation of biosynthesis and the processing of Let-7d using the ADAR1 isoforms should be the subject of new studies and, in the future, could be considered important therapeutic targets for idiopathic pulmonary fibrosis.
Our study has several limitations. At first, we did know how many ADAR1 editing sites have pri-miR-Let-7d and if these sites could be promoting its maturation or its degradation, which has been reported for pri-miR-142, whose degradation is ADAR1mediated-hyper-edition [31]. On the other hand, we did not explore the role of Let-7d in the conformation of the nucleolus in normal fibroblasts compared with fibrotic fibroblasts because it has been reported that Let-7d takes part in the nucleolar organization in mouse lung fibroblasts (MLg and MFML4), mouse lung epithelial cells (MLE-12), mouse mammary gland epithelial cells (NMuMG), and human lung epithelial carcinoma cells (A549); participate in the epigenetic regulation of several genes [32]. Future work will be necessary to investigate the role of ADAR-2 as a possible regulator of Let-7d and determine the role of ADAR-3, an isoform without catalytic activity but with a role of inhibitor of ADAR1 and 2 because it can form inactive dimers in combination with both proteins.

Cell Culture
The study population consisted of primary fibroblasts isolated from control (n = 10) or IPF (n = 10) lung tissue biopsies collected in the frame of the Mexican and the European IPF. All cells were cultured with Ham's/F-12 medium (Thermo-Fischer, San Francisco, CA, USA) supplemented with 10% Fetal Bovine Serum (Thermo-Fischer, San Francisco, CA, USA). Patient characteristics are shown in Table S1.

RNA Isolation, Reverse Transcription, Quantitative PCR
The total RNA of fibroblast and mice lung tissue was extracted with Trizol (Thermo-Fischer, San Francisco, CA, USA) and quantified using the Nanodrop system (Thermo-Scientific, San Francisco, CA, USA). The cDNA synthesis was performed using total RNA and the High-Capacity cDNA Reverse Transcriptase kit (Applied Biosystems, San Francisco, CA, USA). Quantitative real-time PCR reactions were performed using SYBR ® Green (Applied Biosystems, San Francisco, CA, USA) and TaqMan gene expression Real-Time PCR system (Applied Biosystems, San Francisco, CA, USA) on the Step-One-plus Real-time PCR system (Applied Biosystems, San Francisco, CA, USA) according to the manufacturer's instructions as previously described [32]. The housekeeping gene HPRT was used to normalize gene expression. All primer pairs and TaqMan assay used are described in Table S2.

Western Blot
Western blots of total protein were performed as previously described. Proteins were isolated using RIPA lysis buffer (Santa-Cruz Biotechnology, Dallas, TX, USA). Proteins were quantified using the Bradford method (Pierce, CA, USA) in a Beckman DU640 spectrophotometer [32].
All samples were incubated at 95 • C for 5 m with Laemmli buffer, followed by SDS-PAGE. Samples were transferred to nitrocellulose membranes, which were blocked for 1 h at room temperature with 5% fat-free milk in PBS-T 0.05%. Membranes were incubated overnight with a primary antibody at 4 • C (ADAR1, sc-271854; Lamin B, sc-17810, Santa-Cruz Biotechnology), washed with 0.05% PBS-T, and incubated for 1 h at room temperature with an HRP-linked secondary antibody, which was detected with the Super Signal West Dura Extended Duration Substrate detection solutions (Thermo-Scientific, San Francisco, CA, USA). The signal was detected and analyzed with the Chemi-Doc technology system (Bio-Rad, Hercules, CA, USA).

ADAR Catalytic Activity Assay
Editing of RNA in control fibroblasts compared with IPF fibroblasts was measured using the RESS-qPCR [30] as follows. As an edition template, APOBEC3D, the ADAR1mediated editing target, was selected. cDNA of 6 control lines and 6 fibrotic lines (2 µL cDNA per sample) was amplified by qPCR (Table S1), by serial dilutions, and the efficiency of the primers for APOBEC3D wt and APOBEC3D edit was tested, and in a final volume of 10 µL (5 µL #1176202K SYBR GreenER Super Mix Life Technologies, San Francisco, CA, USA. 0.2 µL of primers, calibrated to 10 µL with nuclease-free water). Edition ratios were calculated using the values of the efficiency of the primers in the following way: Relative ratio of RNA edition = E −ct Edit E −ct wt (1) Control and IPF primary cultures were transfected with overexpression constructs for ADAR1-p110 wt, ADAR1-p110 mutant, ADAR1-p150 wt, and ADAR1-p150 mutant, as well as the empty vector. Cells were treated with 0.5 µM of pentostatin (#sc-204177 Santa-Cruz Biotechnology, Santa Cruz, CA, USA) 24 h post-transfection to inhibit the catalytic activity of Adenosine Deaminase (ADA) in order to increase the specificity of this assay. Cells were collected 48 h post-transfection and used for the catalytic activity test.
The formation of inosine (from the hydrolysis of adenosine) was measured in aliquots of the proteins extracted from the cells transfected with specific plasmids that were used for each isoform wild-type and mutant of ADAR1, the colorimetric assay Adenosine Deaminase Activity Assay Kit (#K321-100 Biosvision, Milpitas, CA, USA) was used. The protein extract of the cells was prepared by adding 150µL of cold ADA 1× assay buffer with protease inhibitor (#14865751 Thermo-Scientific, San Francisco, CA, USA). Dilutions of inosine standard were used as a standard curve. As a negative control, 50 µL of 1× ADA Assay Buffer was used. The reading plate was incubated at 37 • C for 5 m, after which it was read at a wavelength of 293 nm using the Infinite 200 plate reader (Tecan, Männedorf, Switzerland) at time points 0, 10, 20, and 40 m. All data were analyzed using the following formula: Editing activity of ADAR per sample = B ∆T · µg de proteína · DF = nmol/min/µg = mU/µg (2) Inosine standard curves were set at 4 time points for the analysis.

Selection of Target mRNA
In silico analysis was performed to investigate Let-7d targets, using miRWalk, miRanda, and TargetScan algorithms. Obtained readouts were compared to select the top 250 mRNA with differential expression in IPF, and Collagen III alpha 1 (COL3A1) and SMAD2 were selected.

Pentostatin Treatment in IPF Fibroblasts
IPF fibroblasts were treated with 0.5 µM of pentostatin (#sc-204177 Santa-Cruz Bio technology, Dallas, TX, USA). At 24 h after, Let-7d and pri-mir-Let-7d expressions were determined by real-time PCR as described before.

Statistical Methods
All data follow a normal distribution, the Student s test was used to compare RNA expression, and values are expressed as mean ± SD. All statistical analyses were performed using GraphPad Prism (San Diego, CA, USA).