Next Article in Journal
DNA Methylation Influences miRNA Expression in Gonadotroph Pituitary Tumors
Previous Article in Journal
Targeting the C-Terminal Domain Small Phosphatase 1
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Throughput Identification of MiR-145 Targets in Human Articular Chondrocytes

1
Kennedy Institute of Rheumatology, University of Oxford, Oxford OX3 7FY, UK
2
Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford OX3 7BN, UK
*
Authors to whom correspondence should be addressed.
Present address: Department of Metabolism, Digestion and Reproduction, Imperial College London, London W12 0NN, UK.
Life 2020, 10(5), 58; https://doi.org/10.3390/life10050058
Submission received: 30 March 2020 / Revised: 30 April 2020 / Accepted: 5 May 2020 / Published: 11 May 2020
(This article belongs to the Section Genetics and Genomics)

Abstract

:
MicroRNAs (miRNAs) play key roles in cartilage development and homeostasis and are dysregulated in osteoarthritis. MiR-145 modulation induces profound changes in the human articular chondrocyte (HAC) phenotype, partially through direct repression of SOX9. Since miRNAs can simultaneously silence multiple targets, we aimed to identify the whole targetome of miR-145 in HACs, critical if miR-145 is to be considered a target for cartilage repair. We performed RIP-seq (RNA-immunoprecipitation and high-throughput sequencing) of miRISC (miRNA-induced silencing complex) in HACs overexpressing miR-145 to identify miR-145 direct targets and used cWords to assess enrichment of miR-145 seed matches in the identified targets. Further validations were performed by RT-qPCR, Western immunoblot, and luciferase assays. MiR-145 affects the expression of over 350 genes and directly targets more than 50 mRNAs through the 3′UTR or, more commonly, the coding region. MiR-145 targets DUSP6, involved in cartilage organization and development, at the translational level. DUSP6 depletion leads to MMP13 upregulation, suggesting a contribution towards the effect of miR-145 on MMP13 expression. In conclusion, miR-145 directly targets several genes involved in the expression of the extracellular matrix and inflammation in primary chondrocytes. Thus, we propose miR-145 as an important regulator of chondrocyte function and a new target for cartilage repair.

1. Introduction

The articular chondrocyte is the only cell type present in the articular cartilage, solely responsible for synthetizing the components of the extracellular matrix (ECM) that allows an almost friction-free, painless movement of the joint. When the articular chondrocyte function is altered, the structure of the ECM can be compromised leading to the development of osteoarthritis (OA) [1]. Local inflammation within the joint may have a detrimental effect in chondrocyte function and further contribute to the development of OA [2]. Currently, there is no cure for OA, which affects over 250 million people worldwide, (WHO http://bjdonline.org/), and most common treatments can only aim to pain reduction.
A normal healthy articular chondrocyte maintains the equilibrium between anabolic molecules, such as collagen type II (COL2A1)—a main component of the ECM—and catabolic enzymes such as the collagen-degrading metalloproteinase MMP13 [1,3]. Understanding the mechanisms underlying the human articular chondrocyte (HAC) function is indispensable to develop therapies aiming to restore the capacity of the HAC to produce a balanced ECM. Failure to develop drugs to stop the development of the disease is partially due to the lack of a comprehensive knowledge of the factors leading to cartilage degeneration and to the poor intrinsic repair capacity of articular chondrocytes [4]. Moreover, joint damage often acts to alter chondrocyte homeostasis of the remaining viable chondrocytes causing further degeneration and/or hypertrophy and matrix calcification over time [5]. Recent research has identified mesenchymal stem cells with chondrogenic potential as a promising source for tissue engineering and replacement, which may suffer less hypertrophy/dedifferentiation following transplantation [6].
MicroRNAs (miRNAs) are small RNAs that silence gene expression post-transcriptionally [7]. Over two thousand miRNAs are expressed in humans [8] that play regulatory roles in most biological processes, development, and diseases [9,10,11]. Animal miRNAs recognize the target mRNAs via partial sequence complementarity and recruit the miRNA-induced silencing complex (miRISC). The MiRISC action results in decreased protein production as a consequence of reduced mRNA translation and/or increased mRNA degradation [7].
MiRNAs have been widely shown to not only regulate cartilage development but also to modulate the normal function of the adult articular chondrocyte [12,13,14]. We have previously identified miR-145 to profoundly impact the expression of essential components of the ECM, such as ACAN, COL2A1, or MMP13, at least partially by directly targeting the master regulator, transcription factor, SOX9 [15]. Interestingly, miR-145 levels are increased in OA chondrocytes [16] raising the possibility of inhibition of miR-145 as a promising approach for cartilage repair. Nevertheless, a single miRNA typically targets dozens to hundreds of mRNAs simultaneously [7], which could result in unpredicted phenotypes upon manipulation. Therefore, if miR-145 is going to be considered as a target for cartilage repair, it is fundamental to determine its whole targetome in HACs. Current available computational methods for miRNA target determination rely in the extent of the complementarity between the seed region of the miRNA (six nucleotides positions 2 ± 7) and the 3′UTR of the mRNAs, which results in hundreds of predictions, most of which do not represent bona-fide interactions [17]. Additionally, the efficiency of these tools tends to be poor for the detection of nonconserved binding sites, binding sites with poor pairing or for those located outside the 3′UTR of the mRNAs [17]. Moreover, computational approaches do not take into account the possibility of tissue-specific interactions [17,18].
Therefore, we aimed to determine miR-145 targets in HACs using an unbiased, high-throughput experimental approach capable of detecting actual miRNA–mRNA interactions [18]. Accordingly, in the present study we performed RNA/miRISC-immunoprecipitation experiments [19,20,21] followed by high-throughput sequencing (RIP-seq) in primary HACs isolated from healthy human cartilage and previously transfected with a miR-145 mimic. Using this method, we detected 59 mRNAs enriched in miRISC-immunoprecipitation upon miR-145 overexpression, indicating direct miRNA–mRNA interactions. Reporter assays confirmed that miR-145 silences the expression of several of these genes through sequences present in their 3′ untranslated region (3′UTR) or coding sequence (CDS). One such identified mRNA, DUSP6, had previously been reported to regulate cartilage function and inflammation of different tissues. The predicted miR-145 binding site in DUSP6 was validated by reporter assays and, as expected, DUSP6 protein levels were modulated in primary HACs upon miR-145 overexpression/inhibition. Both miR-145- and siRNA-mediated downregulation of DUSP6 resulted in increased expression of the catabolic gene MMP13, indicating that DUSP6 is at least partially responsible for the effect that miR-145 exerts in the MMP13 gene expression.

2. Materials and Methods

2.1. Cell Culture and Transfection

Following local ethics committee guidelines and after informed consent, healthy human articular cartilage was extracted and HACs were isolated from the femoral condyle and tibia plateau following amputation due to sarcomas not involving the join space. Cartilage samples were collected on the day of surgery, diced into 1–2 mm pieces and digested with 1.5 mg/mL of collagenase type II (Worthington) in DMEM—10% FCS during 18 h at 37 °C. Isolated human articular chondrocytes (HACs) were then strained, pelleted, washed, and seeded at a density of 1.2 × 104 cells/cm2 as previously described [15]. Primary cells (passage 0—P0) were harvested after 5–7 days of culture or subcultured for a further 5–7 days and passaged one, two, or three times (P1–P3). HACs, as well as HeLa cells, were kept in DMEM supplemented with 10% FCS in a humidified 5% CO2 atmosphere at 37 °C under normal (20%) or hypoxic (1%) oxygen tension, when indicated.

2.2. Plasmids

The sequence containing three consecutive perfect matches for hsa-miR-145-5p (AGGGATTCCTGGGAAAACTGGACAGCGTAGGGATTCCTGGGAAAACTGGACAGCGTAGGGATTCCTGGGAAAACTGGAC) was excised from the previously produced pSG5-Luc-3xmiR-145 15 and subcloned between XhoI–XbaI restriction sites in a pMirGlo dual-luciferase miRNA target expression vector (Promega) to generate pMiRGlo 3xmiR-145.
MOK 3′UTR, MOK CDS, CYR61 3′UTR, HSBP1 3′UTR, RTKN 3′UTR, CTBP1 3′UTR, CTBP1 CDS, FKBP2 CDS, HIST1H4E full-length, GLT8D1 3′UTR, RAB7A 3′UTR, RAB7A CDS, PLOD1 3′UTR, SEC11A 3′UTR, FN3K CDS, UXS1 3′UTR, and DUSP6 3′UTR were amplified by PCR using specific oligonucleotides (sequences available upon request) from HACS cDNA, inserted into the pGEM-T vector (Promega) and subcloned into pMiRGlo.
Three-point mutations in the predicted miR-145 seed binding site in DUSP6 were introduced in pGEM-T-DUSP6 3′UTR using a Phusion® site-directed mutagenesis kit (NEB) and the mutagenic primers: Forward: 5’-TGTGAGCATGGGTACCCATT-3’ and reverse 5’-CACACACACTTCGTCTTTTATACAAA-3’. Mutated DUSP6 3’UTR was subcloned into pMiRGlo. All the constructs were verified by sequencing.

2.3. Transfection of HACs

HACs were seeded at 1 × 104 cells/cm2 and, 12–24 h later, transfected using Lipofectamine 2000 (Thermo Fisher Scientific) and 5 nM control or miR-145 mimics (AM17110, AM17100-PM11480, respectively, Thermo Fisher Scientific), 10–50 nM control or hsa-miR-145 miRCURY LNA Power inhibitors (Cat No. 339136, YI00199006, YI04102423, respectively Qiagen, Netherlands) or 5 nM control or DUSP6 siRNA (On-Target Plus, J-003964-06 (DUSP6 siRNA #1), J-003964-07 (DUSP6 siRNA #2), D-001810-01 (Control), Horizon Discovery) following the manufacturer’s instructions and as previously described [22]. After transfection, the medium was replaced by a pre-equilibrated DMEM (in 20% or 1% oxygen tension) containing 10% FCS and antibiotics. Cells were kept at the appropriate oxygen tension for a further 24–48 h.
HeLa cells were seeded at 60% confluence and, 12–24 h later, transfected with 100 ng of plasmidic DNA together with 5 nM control or miR-145 mimics using Lipofectamine 2000 (Thermo Fisher Scientific). After 24 h, cells were lysed and luciferase activity determined with the Dual-Glo luciferase assay system (Promega, UK), following the manufacturer’s instructions.

2.4. Ribonucleoprotein Immunoprecipitation, mRNA–miRNA Isolation and miRNA Reverse Tanscription (RT), and Real-Time PCR (qPCR)

Ribonucleoprotein immunoprecipitation was performed as described previously [21]. Briefly, 30 × 106 P3 chondrocytes from three different patients (19 yo. female; 8 yo. male; 28 yo. male) were transfected with miR-145 or control mimics and, 24 h later, lysed in an equal volume (300–400 μL) of an ice cold polysome lysis buffer (5 mM MgCl2, 100 mM KCl, 10 mM Hepes, pH 7.0, and 0.5% Nonidet P-40) with freshly added 1 mM DTT, 100 units/mL RNase OUT (Invitrogen), and complete mini EDTA-free protease inhibitors mixture (Sigma-Aldrich) during 5 min, snap frozen, and kept at −80 °C for at least 18 h. Samples were then centrifuged at 14,000 g at 4 °C for 10 min and protein concentration was determined by the Bradford assay 2.
The supernatant containing 6 mg of total protein per condition was mixed with 2.7 mL of an ice cold NT2 buffer (50 mM Tris, pH 7.4, 150 mM NaCl, 1 mM MgCl2, 0.05% Nonidet P-40) containing freshly added 200 units/mL RNase OUT (Thermo Fisher Scientific), 0.5% vanadyl ribonucleoside (Thermo Fisher Scientific), 1 mM DTT, 15 mM EDTA, and 90 µl Sepharose G beads (Abcam, Cambridge, UK) precoated for at least 4 h with 6 μg of mouse antihuman Ago2 (clone2E12-1C9, Abnova). Incubation was carried out overnight at 4 °C on a rocking platform. On the following day, beads were washed six times with a 3 mL ice cold NT2 buffer and 5% of the beads were separated for Western blotting. The rest of the beads were incubated in a 1% SDS-NT2 buffer during 30 min at 50 °C in order to elute the bound RNA that was isolated using TRIzol and glycogen (Invitrogen, UK) as a carrier in the precipitation step. The mouse IgG1 isotype control (Abcam) was used as a negative control for the immunoprecipitation procedure.
In parallel, total RNA was extracted using TRIzol from a 10% volume of the samples submitted to immunoprecipitation, prior to incubation with the Ago2 antibody.
For miRNA detection, 5% of the recovered RNA was reverse transcribed using the TaqMan miRNA RT kit (Thermo Fisher Scientific) and miR-145, miR-140, and RNU24-specific primers (Thermo Fisher Scientific). QPCR was performed with a TaqMan PCR master mix (Thermo Fisher Scientific) and the appropriate TaqMan probes (Thermo Fisher Scientific). The remaining RNA was treated with DNAse I (Thermo Fisher Scientific) and sent for high-throughput sequencing.

2.5. mRNA-Seq Library Construction, Sequencing, and Analysis

RNA quantity and integrity were assessed using the Quant-IT RiboGreen Kit (Invitrogen) and Agilent Tapestation 2200 R6K. mRNA enrichment was achieved from 400–900 ng of total RNA using a Magnetic mRNA Isolation Kit (NEB, UK). Generation of double stranded cDNA and library construction were performed using the NEBNext® mRNA Sample Prep Reagent Set 1 (NEB). The ligation of adapters was performed using adapters prepared at the WTCHG according to the Illumina design (Multiplexing Sample Preparation Oligonucleotide Kit, NEB, UK). Each library was subsequently size selected with two AMpure Bead bindings. The following custom primers were used for PCR enrichment: Multiplex PCR primer 1.0: 5’-AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT-3’. Index primer: 5’-CAAGCAGAAGACGGCATACGAGAT[INDEX]CAGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT-3’. Indices were according to the eight bases tags developed by WTCHG [23]. The following modifications to the described workflow were applied to RNA recovered after ribonucleoprotein immunoprecipitation. Firstly, no poly(A) selection or ribodepletion was applied. Secondly, the number of PCR cycles was increased from 12 to 15. Amplified libraries were analyzed for size distribution using the Agilent Tapestation 2200 D1K. Libraries were quantified by quantitative RT-PCR using the Agilent qPCR Library Quantification Kit and a MX3005P instrument (Agilent, Germany) and relative volumes were pooled accordingly. Finally, a second RT-qPCR was performed to measure the relative concentration of the pool compared to a previously sequenced mRNA library in order to determine the volume to use for sequencing. Sequencing was performed as a 50 bp paired end read on a HiSeq2000 according to Illumina specifications.
FASTQ files were generated for each sample and initial data quality checks of the raw sequence data were performed. Reads were then mapped to the human genome (GRCh37) using TopHat [24]. HT-seq [25] was used to summarize the mapped reads into a gene count table, the final version of which contained one column per sample. Ensembl (v65) annotations were used to define gene models. All subsequent data analysis was performed in R using relevant Bioconductor packages [26].
Total RNA profiles were consistent and generated ~20 million reads mapping to Ensembl genes per sample. IP samples generated 2–3 million reads per sample and one potential outlier sample was observed on clustering and PCA plots (Figure S1). The R package ‘DESeq’ [27] was used to normalize the samples for library depth and perform a paired analysis to find differentially expressed genes between the samples overexpressing miR-145 and controls. Total RNA and IP samples were treated as separate datasets. The two resulting gene lists were compared to identify genes that were significantly upregulated in the IP samples while being significantly downregulated or unchanged (not significant, FDR > 0.05) in the total RNA samples as potential targets of miR-145. For the analysis performed with IP samples from only two patients, we investigated the effect of the outlier and the final IP list was generated excluding both samples from the patient (donor 3) with the outlier sample (miR-145 donor 3). Genes were filtered by expression levels and those with log2CPM > 2 from total RNA and IP samples were independently ranked by fold change of expression and subjected to the motif discovery method cWords in two separate analysis per list, one for enrichment in the 3′UTRs and the other in the CDS [28].

2.6. mRNA Reverse Transcription (RT) and Real-Time PCR (qPCR)

Total RNA was extracted from P0–P3 HACs using TRIzol according to the manufacturer’s instructions. An amount of 200–500 ng of RNA was reverse transcribed using the High Capacity cDNA RT kit (Thermo Fisher Scientific). QPCR followed, using a SYBR Green PCR master mix (Thermo Fisher Scientific) and specific primers (sequences available upon request). The age and gender of the HACs donors used for each experiment are indicated in the figure legends.

2.7. Western Blotting

P0–P3 HACs were cultured as monolayers and transfected 36–48 h before lysis. When appropriate, HACs were stimulated with a 25 nM human FGF-2 (Prepotech, UK). The age and gender of the HACs donors used for each experiment are indicated in the figure legends. HACs were lysed in a radio-immunoprecipitation assay buffer (RIPA—150 mM NaCl, 10 mM Tris, pH 7.2, 5 mM EDTA, 1% Triton X-100, 0.1% SDS, 1% deoxycholic acid) and 10–30 μg total protein extract subjected to SDS-polyacrilamide gel electrophoresis and transferred to PVDF membranes (Merck-Millipore, Germany). Mouse monoclonal anti-Ago2 (clone2E12-1C9, Abnova, 1:1000 dilution), rabbit monoclonal anti-DUSP6 (EPR129Y, Abcam, 1:1000 dilution), rabbit monoclonal anti-phospho-ERK (Cat No. 05-797R, Sigma-Aldrich, 1:1000 dilution), rabbit polyclonal anti-ERK (Cat No. sc-292838, SantaCruz Biotechnologies, 1:500 dilution), mouse monoclonal anti-Tubulin (clone B-5-1-2, Sigma-Aldrich, 1:5000 dilution), HRP-conjugated anti-mouse IgG secondary (Cat No. NA931, GE Healthcare, 1:5000 dilution), and HRP-conjugated anti-rabbit IgG secondary (Cat No. NA934, GE Healthcare, 1:5000 dilution), antibodies were used. Proteins were visualized by ECL fluorography (GE Healthcare, Germany).

2.8. Statistical Analysis

GraphPad Prism 6.0 was used for statistical analysis. Statistical significance was evaluated by the two tailed paired student’s t-test or two-way ANOVA and Fisher’s LSD test for multiple comparisons. All data are shown as mean ± standard error of the mean (SEM). p-values < 0.05 were considered statistically significant.

3. Results

As previously mentioned, miRNAs direct miRISC to the target mRNAs which, as a consequence, are less efficiently translated and/or degraded. Thus, regardless of the subsequent mechanism of repression, an mRNA targeted by miRNAs will be bound to the miRISC. Consequently, overexpression of a specific miRNA is expected to increase the binding of its mRNA targets to miRISC. MiRISC can be immunoprecipitated using antibodies against miRISC proteins such as Argonaute 2 (AGO2), a key component of the complex.
In order to determine miR-145 direct targets in human chondrocytes we performed mRNA/miRISC-immunoprecipitation followed by high-throughput sequencing of the RNA (RIP-seq) with an anti-Ago2 antibody and a lysate from P3 HACs previously transfected with a miR-145 mimic or a nontargeting control. In parallel, to assess the overall impact of miR-145 overexpression in gene expression in HACs, we sequenced total mRNA (T-RNA) of HACs transfected with miR-145 or control mimics. Upon miR-145 overexpression, miR-145 target mRNAs are expected to be enriched in Ago2 precipitates while their levels in the total RNA samples should remain unchanged or be reduced (depending on whether miR-145 suppresses target mRNA translation or mRNA stability, respectively).
The specificity of Ago2-immunoprecipitation was verified by Western blotting (Figure 1A) and hundreds of miRNAs were detected in Ago2 immunoprecipitates in comparison with an IgG control (data not shown). The amount of miR-145 associated with Ago2 was 40-fold higher in cells overexpressing miR-145 (Figure 1B). Immunoprecipitated RNA (IP-RNA), as well as total mRNA (T-RNA), was obtained from three independent experiments with chondrocytes isolated from three different donors and prepared for high-throughput sequencing.
The expression of 352 genes was significantly changed (FDR < 0.05) following miR-145 overexpression (T-RNA, Table S1). Out of these genes, only 80 were downregulated (DownT-RNA, Table 1) suggesting that most of the effects of miR-145 on gene expression are indirect. Ninety-six genes were enriched in IP-RNA from samples overexpressing miR-145. Out of these genes, 59 (Table 2) were simultaneously reduced or unchanged in the corresponding T-RNA samples. Since miRNAs normally act by decreasing target-mRNA stability or translation efficiency we considered those 59 mRNAs as miR-145 direct targets (Di-IP, Table 2).
Although miRNAs can also mediate repression through “seedless” sites [29] and/or through the CDS [30], miR-145 direct targets are generally expected to contain miR-145 seed matches in the 3′UTRs. TargetScan (TargetScan (v6.2) [31] http://www.targetscan.org/) detected miR-145 target sites in 17 out of our 59 candidates (Table 2) whereas RNAhybrid, a bioinformatic tool that gives little weight to the seed region and allows the analysis of any given mRNA sequence, including CDSs or 3′UTRs, identified miR-145 binding sites within the 3′UTR, the CDS, and both the CDS and 3′UTR of 8, 20, and 12 candidate genes, respectively (Table 2). We could not detect any predicted binding sites in 12 out of the 59 candidate genes (Table 2). Additionally, we used cWords [28] to analyze words enriched in 3′UTRs and CDSs of IP-RNAs (Figure 2). As expected, the word with the strongest enrichment in CDSs of upregulated genes corresponds to the 6mer-seed site of miR-145 (ACTGGA—Figure 2B) and several shifted variants of the seed site were also enriched [32]. Surprisingly, words corresponding to miR-145 seeds were not significantly enriched in 3′UTRs (Figure 2C,D). Altogether, our data suggest that genes upregulated in IP-RNA (Di-IP) are miR-145 direct targets and that miR-145-target interactions occur predominantly through their CDS.
Further supporting that most of the changes observed in gene expression following miR-145 modulation are indirect, cWords failed to detect an enrichment in miR-145 seed-matching sequences in genes upregulated or downregulated at the mRNA level (data not shown). Moreover, out of the 80 genes in which expression was downregulated by the action of miR-145 (DownT-RNA), only one was included in the Di-IP group. Due to the complexity of the immunoprecipitation protocol and the low amount of mRNA recovered, we could not rule out the possibility that the miR-145–mRNA weak interactions were lost. Indeed, clustering and PCA plots (Figure S1) suggested one potential outlier (IP-RNA, sample corresponding to miR-145 overexpression, donor 3). Thus, we expanded the analysis excluding the outlier (donor 3) and generated an additional list of target candidates (Table S2), with genes that were significantly enriched in IP-RNA from only two patients (1 and 2) and simultaneously reduced or unchanged in the T-RNA samples (donors 1–3). This additional list contained several genes, which expression was repressed by miR-145 action at the RNA level (DownT-RNA, Table 1, genes marked with *), and the gene Sox9, which we have previously demonstrated to be a direct miR-145 target, regulated exclusively at the protein level in HACs [15]. Consistently, words corresponding to the miR-145 seed were clearly enriched in both CDSs and 3′UTRs of these genes (Figure S2).
As mentioned above, our original RIP-seq results (with three HACs donors) suggest that all the direct target genes identified (Di-IP) but one are repressed by miR-145 at the level of translation, since no significant reduction in their mRNAs was observed upon miR-145 overexpression. Thus, to further validate our findings we set to confirm an effect of miR-145 at the protein level of 12 randomly selected genes from this list. Accordingly, we generated luciferase reporters containing the 3′UTR (CYR61, HSBP1, RTKN, GLT8D1, and PLOD1, with predicted sites only in the 3′UTR, Table 2), the CDS (FKBP2, CTBP1, and FN3K, with predicted sites only in the CDS, Table 2) or both (MOK, RAB7A, and SEC11A with predicted sites in both 3′UTR and CDS and HIST1H4E, lacking predicted sites) of these mRNAs downstream of the firefly luciferase ORF.
Supporting silencing at the protein level, miR-145 overexpression significantly reduced firefly luciferase activity (Figure 3) in the presence of CYR61, HSBP1, GLT8D1, RTKN and PLOD1 3′UTRs, FN3K, MOK and RAB7A CDS, HIST1H4E full-length sequence and SEC11A both 3′UTR and CDS. Only two of the sequences tested (CTBP1 and FKBP2 CDS) failed to mediate silencing. Interestingly, no target sites were predicted for HIST1H4E mRNA, but miR-145 overexpression significantly reduced luciferase expression in the presence of full-length HIST1H4E cDNA, suggesting that HIST1H4E is targeted by miR-145 and demonstrating that our approach can detect targets that would not be predicted by most bioinformatics tools.
It has been previously suggested that miRNA action on translational repression precedes mRNA degradation [33]. Therefore, we hypothesize that miR-145 action in some of our newly identified targets may eventually result in mRNA degradation. Thus, we set to quantify mRNA expression of 16 randomly selected genes from the Di-IP list (Table 2) in HACs transfected with miR-145 mimics for 48 h (instead of 24 h). HACs become gradually dedifferentiated with passage and, therefore, we decided to use freshly isolated (fully differentiated, P0) HACs for these experiments. Accordingly, we overexpressed miR-145 or a control mimic in P0 HACs isolated from seven different patients and analyzed the expression of the selected 16 miR-145 direct targets by qPCR (Figure 4A). Even though their mRNA levels remained unchanged 24 h after transfection (Table 1, Table S3), four genes were significantly (p < 0.5) downregulated and three showed a tendency to be downregulated (p < 0.1) 48 h after transfection. As expected, CC2DB1 mRNA, which was the only direct miR-145 target which mRNA was significantly downregulated at the mRNA level 24 h after transfection, was also significantly reduced at the 48 h mark. These data suggest that the miR-145-target interaction in HACs results first in translational repression and it is later on followed by mRNA destabilization.
We also validated the effect of miR-145 in three of the genes present in DownT-RNA but not in Di-IP, showing a significant effect for FSCN1 and UXS1, but not PPP3CA (Figure 4B). We additionally corroborated the miR-145 effect in COL2A1, which we have previously shown [15] to be indirectly regulated by miR-145 (Figure 4B).
In order to identify the biological processes and pathways in which miR-145 is involved, we submitted the full list of genes altered upon miR-145 overexpression (Table S1, T-RNA FDR < 0.05) together with those included in the Di-IP list (Table 2), to the Database for Annotation, Visualization and Integrated Discovery (DAVID v 6.7). Significantly enriched pathways included immune response and proteolytic processes, both characteristic of cartilage disease (Table 3) (p < 0.1).
One of the genes involved in immunity, DUSP6, encodes a dual-specificity phosphatase, MAP kinase phosphatase 3 (MKP3) that inhibits MAPK activity by dephosphorylating threonine and serine residues. MKP3 regulates the elongation of cartilage elements of digit 2 during development through repression of FGF signaling [34]. Interestingly, loss of DUSP6 in mouse embryos leads to dominant, incompletely penetrant, and variable phenotypes including hearing loss, coronal craniosynostosis, and skeletal dwarfism. The latter phenotype has also been observed upon mutations in Sox9 [35]. In addition, chondrocytes in the mutant’s growth plate-proliferating zone were disorganized [36]. Thus, we hypothesized that miR-145-mediated DUSP6 regulation affects adult chondrocyte function. As expected, miR-145 overexpression in HACs led to a strong reduction of DUSP6 protein (Figure 5A,B) in both regular (20%) and low (1%, hypoxia) oxygen tension, which is the normal physiologic condition for HACs and promotes their differentiated phenotype. Conversely, prevention of the miR-145 function with specific inhibitors resulted in increased DUSP6 expression (Figure 5C). As previously shown, miR-145 overexpression causes only a small [35], not significant, decrease in DUSP6 mRNA (Figure 4A), suggesting that its effect occurs mainly at the translational level. A miR-145 binding site was predicted in the position 1018–1025 of DUSP6 3′UTR. In order to validate the functionality of this binding site, we generated reporters containing the wild type DUSP6 3′UTR or a mutated version, with point mutations in three nucleotides of the target site, downstream of the luciferase ORF. Transfection of miR-145 mimics significantly reduced luciferase activity of the wild type construct but failed to silence luciferase activity in the presence of the mutations (Figure 5D). Our data confirm that miR-145 represses DUSP6 expression through the predicted binding site in its 3′UTR.
FGF2 (fibroblast growth factor 2) plays an anti-anabolic and/or catabolic role in HACS [37,38] at least partially via the activation of the Ras-Raf-MEK1/2-ERK1/2 pathways [39]. FGF2 is present in OA synovial fluid, where it activates Runx2 and promotes MMP13 expression [3]. Levels of FGF2 and MMP13 are increased in human OA cartilage and, in adult human articular cartilage, FGF2 stimulates MMP13 through the activation of multiple MAPKKs [40]. Additionally, mice in which increased Erk1/2 phosphorilation occurred upon heparin endosulfatases (Sulf1 and Sulf2) depletion, showed increased MMP13 levels and reduced Col2a1 and ACAN [41] expression. Interestingly, a raise in Runx2 and MMP13 and a reduction in Col2a1 and ACAN was observed upon miR-145 overexpression 15.
In mouse embryos, transcription of DUSP6 is dependent on FGF signaling and targeted inactivation of DUPS6 increased the levels of pERK, the main target described for DUSP6 [36]. Thus, we hypothesized that DUSP6 contributes to the effect of miR-145 in the expression of the catabolic/anabolic molecules mentioned above. Whereas DUSP6 depletion with specific siRNAs did not alter Col2a1 or ACAN levels (data not shown), expression of MMP13 was strongly upregulated (Figure 6A), indicating that miR-145-dependent MMP13 regulation may be partially mediated by DUSP6. Nevertheless, miR-145 overexpression or siRNA-mediated DUSP6 depletion failed to impair ERK phosphorylation in response to FGF2 (Figure 6B,C). Altogether, these data indicate that miR-145 function can be modified without significantly interfering with FGF2-dependant ERK activation and the effect of DUSP6 in MMP13 expression occurs through a still unidentified ERK1/2-independent pathway.

4. Discussion

We have applied RIP-seq in combination with total mRNA-seq following miR-145 overexpression to successfully identify > 50 new direct miR-145 targets in primary human articular chondrocytes. MiRNA-target interaction results in translational repression accompanied or not of mRNA degradation [42,43,44]. Therefore, experimental approaches that rely on transcriptome profiling alone fail to detect targets regulated exclusively at the level of translation or to differentiate between direct and indirect targets. Surprisingly, most miR-145-mRNA interactions in HACs did not result in mRNA degradation 24 h after overexpression. However, a reduction in several of these transcripts was observed at a later time point and luciferase assays further confirmed that miR-145 overexpression results in lower protein (luciferase) levels for the majority of the cases (10 out of 12). Therefore, our results suggest that miR-145 functions by inhibiting translation of most of its direct targets and this is later followed by mRNA decay. The main mechanism by which miRNAs repress gene expression remains highly controversial [42,45,46]. Interestingly, recent studies suggest that the mechanism of action might be cell-type specific [47] and dependent on the target site type. Thus, Liu et al. found that extensive complementarity to the seed region of the miRNA through the 3′UTR often results in mRNA decay instead of translational repression [43]. Nevertheless, it remains unclear which features favor mRNA translation over decay [43]. We demonstrate here that, in HACs, miR-145 binding occurs more often through CDSs. Recent studies have shown that interaction through the CDS occurs frequently [22,30,48]. Whether the region of interaction miR-145/mRNA influences the mechanism of silencing in HACs remains to be clarified. While most identified miR-145 targets contained predicted miR-145 binding sites for the miRNA, 12 genes (20%) were not predicted as miR-145 targets by TargetScan or RNAhybrid. One of these genes, HIST1H4E, was within the list of targets selected randomly for further validation. Overexpression of miR-145 reduced luciferase activity in the presence of the full HIST1H4E mRNA, suggesting that miR-145 acts on this sequence to repress translation, possibly via a noncanonical binding site. MiRNA action through noncanonical and seed-independent binding sites is well documented [49,50].
Even though our RIP-seq and luciferase reporter results suggest that our list of direct targets genes (Table 2) are silenced by direct binding of miR-145 to their mRNAs, similar results could potentially be obtained if miR-145 overexpression resulted in altered expression of other miRNAs. Luciferase assays with reporters containing mutations in the predicted miR-145 binding sites would further confirm a miR-145-target mRNA direct interaction. While we have performed these experiments for DUSP6, further mutation assays are outside the scope of this work. Moreover, whereas we have previously observed SOX9-dependent changes in miR-140 and miR-365 [15] expression 48 h after miR-145 overexpression, the RIP-seq and reporter experiments presented here have been performed within 24 h of miR-145 overexpression, to limit the occurrence of indirect effects.
Importantly, computational tools for miRNA target predictions typically observe hundreds to thousands of targets (905 transcripts contain conserved miR-145 binding sites according to TargetScan), including a high number of false-positive predictions [18]. Our approach allowed us to identify miR-145 targets in a cell-specific manner, reducing the number of false positives.
MiR-145 overexpression resulted in upregulation of a high number of genes (272—Table S1) and downregulation of > 80 genes (Table 1 and Table S1), suggesting indirect effects. This is not entirely surprising, since, as we have shown in the past, an important target of miR-145 is SOX9 [15], a transcription factor that regulates a wide range of genes in HACs. Further supporting this idea, only one of the downregulated genes was enriched in Ago-IPs although several additional genes from this list were enriched in Ago-IPs upon exclusion from the analysis of one outlier sample. Thus, our immunoprecipitation protocol may have failed to recover weak mRNA–miRISC interactions. mRNA crosslinking may increase the yield of recovered interactions, although it could lead to potential artifacts due to reduced cell lysis efficiency, introduce sequence biases, increase background, or be incompletely reversible [51,52]. Increasing the number of RIP-seq experiments would have definitely resulted in a more comprehensive list of actual targets, but this was limited both by the availability of healthy primary chondrocytes and the high cost of the HT-seq, rapidly decreasing as the technology is optimized. Importantly, mRNA-RISC is a dynamic interaction, which can be hugely affected by the cellular context: Our experiments have been performed in primary human chondrocytes, instead in the more accessible cell lines or mouse models that are commonly used in miRNA studies in cartilage, increasing the reliability and applicability of our results.
Although it remains to be demonstrated, several of the newly identified miR-145 targets in HACs may play important roles in cartilage. For example, POLDIP2 may modulate the production of MMP1 through NOX4 [53,54]. Moreover, PLOD1 could potentially be involved in collagen metabolism [55] and ERF, an ETS2 repressor factor, may indirectly modulate the expression of ECM components [56,57]. Interestingly, CYR61 has already been shown to play a role in chondrogenesis [58], although its function in adult cartilage remains to be clarified.
MiR-145 is abundantly expressed in cartilage, but is also highly expressed in various mouse tissues such as aorta, where miR-145 affects its structure [59]. Moreover, changes in miR-145 levels have been associated with various human diseases, including stroke [60], arteriosclerosis [61], or cancer [62,63]. Although miRNA-target interactions can be cell-specific, our list of targets in HACs could be of interest to researchers in other fields and some of the interactions may be relevant for other tissues or diseases. Supporting this idea, RTKN has been already proven to be targeted by miR-145 in breast cancer [62], FSNC1 in colon cancer and melanoma [64,65], and GOLM1 in prostate cancer [66].
A growing number of studies have associated changes in the expression of dozens of miRNAs to the development of OA [67]. Well studied examples include protective miR-140 [68] or detrimental miR-138 [69]. Moreover, miRNAs have also been proposed as biomarkers of the disease [67] that could facilitate early OA detection. This is important because symptoms are often shown once the cartilage is too degraded for intrinsic repair. Nevertheless, the specific function and gene targets of most of these miRNAs, and thus whether OA-related changes contribute or are consequential to the disease, remains unknown. The experimental approach used here can be applied to the study of any other chondrocyte miRNA and will serve to understand whether a given miRNA has therapeutic potential for OA.
OA pathology is characterized by progressive cellular and molecular changes in joint tissues. While it has always been considered a degenerative joint disease influenced by mechanical stresses and ageing causing cartilage breakdown, discoveries in the recent years point towards a critical contribution of the immune system and inflammatory mechanisms [2,70]. Chondrocytes in the joint can express and respond to cytokines and chemokines such as TNF-a, IL-1 or IL-6, and both IL-1β and Interferon-γ suppress COL2A1 transcription. Moreover, inflammatory signals can contribute to macrophage activation within the osteoarthritic joint, which can themselves produce cytokines such as IL-1 or TNFα [2]. Interestingly, Yang et al. [16] have recently found that miR-145 is upregulated in OA cartilage and chondrocytes, whereas HACs treated with IL-1β contained higher miR-145 levels. Moreover, inhibition of miR-145 partially reversed the IL-1β-mediated reduction in COL2A1 and ACAN as well as the stimulation of MMP-13. On the contrary, salidroside (SAL), a compound with pharmacological effects suggested to accelerate fracture healing [71], has been recently shown to alleviate LPS-mediated inflammatory injury in a chondrogenic cell line by, at least partially, stimulating miR-145 expression [72]. Our gene ontology analysis of the genes affected by overexpression of miR-145 show enrichment in pathways involved in immunity and inflammatory and interferon responses. Not surprisingly, a role for miR-145 in modulating interferon response has already been suggested in very different cellular contexts such as bladder cancer [73] or fish INF-γ-induced immune response [74]. Nevertheless, the mechanisms underlying the role of miR-145 in these pathways remain unknown. Whereas miR-145 induces INF-β by targeting the suppressor of cytokine signaling 7 (SOCS7), we did not detect an effect of miR-145 in SOCS7 expression or association with miRISC. Thus, miR-145 most probably modulates the inflammatory response in HACs through a SOCS7-independent mechanism. From our list of miR-145 targets in HACs, DUSP6 stood up as an interesting candidate to control inflammatory responses in HACs. DUSP6 has been previously shown to mediate the response of macrophages to LPS [75] and plays dual roles in the regulation of vascular inflammation [76]. Moreover, DUSP6 acts on FGF signaling to control cartilage digit development [34] and growth-plate chondrocytes organization [36] during mice development. Our work demonstrates that miR-145 silences DUSP6 expression at the protein level in primary adult HACs. In adult HACS, FGF2-ERK activation promotes the expression of MMP13 and other catabolic molecules [37,38], which are often increased in OA. Here, we show that the effect of miR-145 in MMP13 expression, previously described by Yang [16] and ourselves [15] might be at least partially mediated by DUSP6. Nevertheless, in our hands, the levels of pERK upon FGF2 stimulation were similar in the absence of DUSP6 and/or presence of miR-145, suggesting that the effect of DUSP6 on MMP13 is independent of FGF2-mediated ERK phosphorylation. Further studies, especially those mimicking cartilage diseases, are required to confirm these findings and to shed more light into the mechanism of action of this phosphatase in HACs.
In summary, we demonstrate that miR-145 directly targets dozens of genes in the human articular chondrocyte to strongly impact its function. We have previously proposed that miR-145 inhibition in OA cartilage results in increased synthesis of key ECM components, such as COL2A1, and silencing of degradative enzymes through its target SOX9 [22]. Our new findings support the argument that OA cartilage treated with miR-145 inhibitors would additionally benefit from a reduction in the inflammatory responses that contribute to cartilage degeneration. MiRNAs are hot candidates as drugs/drug targets and promising advances on the delivery methods and stability of both miRNA mimics and inhibitors (often referred to as antagoMiRs) have been made during the past decade [77]. Articular cartilage is poorly vascularized, alymphatic and aneural, and cartilage cells are embedded in a thick ECM, impairing effective pharmacologic intervention. Thus, autologous chondrocyte implantation has been suggested as a promising approach for cartilage repair while MSCs have been proposed as an exciting source of chondrocytes for autologous or allogenic transplantation. Suppression of miR-145 has been shown to promote chondrogenic differentiation of MSCs [78] in vitro. Our results suggest that inhibition of miR-145 could also prevent hypertrophy of the transplanted tissue. The effect of modulating the expression of miR-145 target genes for chondrocyte function will need to be carefully assessed prior to the development of therapeutic approaches targeting miR-145. Experiments performed in vivo will surely address the suitability of miR-145-based therapies for the treatment of OA in the near future.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-1729/10/5/58/s1, Figure S1: Principal component analysis (PCA), Figure S2: Enrichment of miR-145 seed sites in the CDSs (A, B) or 3’UTRs (C, D) of the precipitated mRNAs (outlier excluded, only 2 experiments). Table S1: Gene expression (Fold change miR-145 overexpression versus control) FDR < 0.05, Table S2: mRNAs significantly up-regulated in Ago2-IP from cells over-expressing miR-145 (Only un-changed or down-regulated in totals; outlier excluded, only 2 experiments), Table S3: Gene expression (Fold change miR-145 overexpression versus control) All detected Genes.

Author Contributions

Conceptualization, A.M.-S. and C.L.M.; methodology, A.M.-S.; software, A.M.-S., E.S. and H.L.; validation, A.M.-S. and S.L.; formal analysis, A.M.-S. and E.S.; resources, C.L.M. and H.L.; data curation, A.M.-S.; writing—original draft preparation, A.M.S.; writing—review and editing, A.M.-S. and C.L.M.; funding acquisition, C.L.M. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

We thank the High-Throughput Genomics Group, especially Paolo Piazza, at the Wellcome Trust Centre for Human Genetics (funded by the Wellcome Trust grant reference 090532/Z/09/Z and MRC Hub grant G0900747 91070) for the advice on the preparation of samples for HT-seq and the generation of the sequencing data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Billinghurst, R.C.; Dahlberg, L.; Ionescu, M.; Reiner, A.; Bourne, R.; Rorabeck, C.; Mitchell, P.; Hambor, J.; Diekmann, O.; Tschesche, H.; et al. Enhanced cleavage of type II collagen by collagenases in osteoarthritic articular cartilage. J. Clin. Investig. 1997, 99, 1534–1545. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Woodell-May, J.E.; Sommerfeld, S.D. Role of Inflammation, and the Immune System in the Progression of Osteoarthritis. J. Orthop. Res. 2020, 38, 253–257. [Google Scholar] [CrossRef]
  3. Wang, X.; Manner, P.A.; Horner, A.; Shum, L.; Tuan, R.S.; Nuckolls, G.H. Regulation of MMP-13 expression by RUNX2 and FGF2 in osteoarthritic cartilage. Osteoarthr. Cartil. 2004, 12, 963–973. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Caldwell, K.L.; Wang, J. Cell-based articular cartilage repair: The link between development and regeneration. Osteoarthr. Cartil. 2015, 23, 351–362. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Dreier, R. Hypertrophic differentiation of chondrocytes in osteoarthritis: The developmental aspect of degenerative joint disorders. Arthritis Res. Ther. 2010, 12, 216. [Google Scholar] [CrossRef] [Green Version]
  6. Colombini, A.; Perucca Orfei, C.; Kouroupis, D.; Ragni, E.; De Luca, P.; Vigano, M.; Correa, D.; de Girolamo, L. Mesenchymal stem cells in the treatment of articular cartilage degeneration: New biological insights for an old-timer cell. Cytotherapy 2019, 21, 1179–1197. [Google Scholar] [CrossRef]
  7. Bartel, D.P. microRNAs: Target recognition and regulatory functions. Cell 2009, 136, 215–233. [Google Scholar] [CrossRef] [Green Version]
  8. Kozomara, A.; Griffiths-Jones, S. miRBase: Annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014, 42, 68–73. [Google Scholar] [CrossRef] [Green Version]
  9. Sun, K.; Lai, E.C. Adult-specific functions of animal microRNAs. Nat. Rev. Genet. 2013, 14, 535–548. [Google Scholar] [CrossRef]
  10. Emde, A.; Hornstein, E. miRNAs at the interface of cellular stress and disease. EMBO J. 2014, 33, 1428–1437. [Google Scholar] [CrossRef] [Green Version]
  11. Kloosterman, W.P.; Plasterk, R.H. The diverse functions of microRNAs in animal development and disease. Dev. Cell 2006, 11, 441–450. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Miyaki, S.; Sato, T.; Inoue, A.; Otsuki, S.; Ito, Y.; Yokoyama, S.; Kato, Y.; Takemoto, F.; Nakasa, T.; Yamashita, S.; et al. microRNA-140 plays dual roles in both cartilage development and homeostasis. Genes Dev. 2010, 24, 1173–1185. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Swingler, T.E.; Wheeler, G.; Carmont, V.; Elliott, H.R.; Barter, M.J.; Abu-Elmagd, M.; Donell, S.T.; Boot-Handford, R.P.; Hajihosseini, M.K.; Münsterberg, A.; et al. The expression and function of microRNAs in chondrogenesis and osteoarthritis. Arthritis Rheum. 2012, 64, 1909–1919. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Shang, J.; Liu, H.; Zhou, Y. Roles of microRNAs in prenatal chondrogenesis, postnatal chondrogenesis and cartilage-related diseases. J. Cell Mol. Med. 2013, 17, 1515–1524. [Google Scholar] [CrossRef]
  15. Martinez-Sanchez, A.; Dudek, K.A.; Murphy, C.L. Regulation of human chondrocyte function through direct inhibition of cartilage master regulator SOX9 by microRNA-145 (miRNA-145). J. Biol. Chem. 2012, 287, 916–924. [Google Scholar] [CrossRef] [Green Version]
  16. Yang, B.; Kang, X.; Xing, Y.; Dou, C.; Kang, F.; Li, J.; Quan, Y.; Dong, S. Effect of microRNA-145 on IL-1β-induced cartilage degradation in human chondrocytes. FEBS Lett. 2014, 588, 2344–2352. [Google Scholar] [CrossRef] [Green Version]
  17. Hamzeiy, H.; Allmer, J.; Yousef, M. Computational methods for microRNA target prediction. Methods Mol. Biol. 2014, 1107, 207–221. [Google Scholar]
  18. Martinez-Sanchez, A.; Murphy, C.L. MicroRNA Target Identification-Experimental Approaches. Biology 2013, 2, 189–205. [Google Scholar] [CrossRef] [Green Version]
  19. Baroni, T.E.; Chittur, S.V.; George, A.D.; Tenenbaum, S.A. Advances in RIP-chip analysis: RNA-binding protein immunoprecipitation-microarray profiling. Methods Mol. Biol. 2008, 419, 93–108. [Google Scholar]
  20. Tan, L.P.; Seinen, E.; Duns, G.; de Jong, D.; Sibon, O.C.; Poppema, S.; Kroesen, B.J.; Kok, K.; van den Berg, A. A high throughput experimental approach to identify miRNA targets in human cells. Nucleic Acids Res. 2009, 37, e137. [Google Scholar] [CrossRef]
  21. Keene, J.D.; Komisarow, J.M.; Friedersdorf, M.B. RIP-Chip: The isolation and identification of mRNAs, microRNAs, and protein components of ribonucleoprotein complexes from cell extracts. Nat. Protoc. 2006, 1, 302–307. [Google Scholar] [CrossRef] [PubMed]
  22. Martinez-Sanchez, A.; Murphy, C.L. miR-1247 functions by targeting cartilage transcription factor SOX9. J. Biol. Chem. 2013, 288, 30802–30814. [Google Scholar] [CrossRef] [Green Version]
  23. Lamble, S.; Batty, E.; Attar, M.; Buck, D.; Bowden, R.; Lunter, G.; Crook, D.; El-Fahmawi, B.; Piazza, P. Improved workflows for high throughput library preparation using the transposome-based Nextera system. BMC Biotechnol. 2013, 13, 104. [Google Scholar] [CrossRef] [Green Version]
  24. Trapnell, C.; Roberts, A.; Goff, L.; Pertea, G.; Kim, D.; Kelley, D.R.; Pimentel, H.; Salzberg, S.L.; Rinn, J.L.; Pachter, L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 2012, 7, 562–578. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Anders, S.; Pyl, P.T.; Huber, W. HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics 2014, 31, 166–169. [Google Scholar] [CrossRef] [PubMed]
  26. Anders, S.; McCarthy, D.J.; Chen, Y.; Okoniewski, M.; Smyth, G.K.; Huber, W.; Robinson, M.D. Count-based differential expression analysis of RNA sequencing data using R and Bioconductor. Nat. Protoc. 2013, 8, 1765–1786. [Google Scholar] [CrossRef] [PubMed]
  27. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef] [Green Version]
  28. Rasmussen, S.H.; Jacobsen, A.; Krogh, A. cWords-systematic microRNA regulatory motif discovery from mRNA expression data. Silence 2013, 4, 2. [Google Scholar] [CrossRef] [Green Version]
  29. Lal, A.; Navarro, F.; Maher, C.A.; Maliszewski, L.E.; Yan, N.; O’Day, E.; Chowdhury, D.; Dykxhoorn, D.M.; Tsai, P.; Hofmann, O.; et al. miR-24 Inhibits cell proliferation by targeting E2F2, MYC, and other cell-cycle genes via binding to “seedless” 3’UTR microRNA recognition elements. Mol. Cell 2009, 35, 610–625. [Google Scholar] [CrossRef] [Green Version]
  30. Hafner, M.; Landthaler, M.; Burger, L.; Khorshid, M.; Hausser, J.; Berninger, P.; Rothballer, A.; Ascano, M., Jr.; Jungkamp, A.C.; Munschauer, M.; et al. Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP. Cell 2010, 141, 129–141. [Google Scholar] [CrossRef] [Green Version]
  31. Friedman, R.C.; Farh, K.K.; Burge, C.B.; Bartel, D.P. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19, 92–105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Lewis, B.P.; Burge, C.B.; Bartel, D.P. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 2005, 120, 15–20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Djuranovic, S.; Nahvi, A.; Green, R. miRNA-mediated gene silencing by translational repression followed by mRNA deadenylation and decay. Science 2012, 336, 237–240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Uejima, A.; Amano, T.; Nomura, N.; Noro, M.; Yasue, T.; Shiroishi, T.; Ohta, K.; Yokoyama, H.; Tamura, K. Anterior shift in gene expression precedes anteriormost digit formation in amniote limbs. Dev. Growth Differ. 2010, 52, 223–234. [Google Scholar] [CrossRef] [PubMed]
  35. Foster, J.W.; Dominguez-Steglich, M.A.; Guioli, S.; Kwok, C.; Weller, P.A.; Stevanović, M.; Weissenbach, J.; Mansour, S.; Young, I.D.; Goodfellow, P.N. Campomelic dysplasia and autosomal sex reversal caused by mutations in an SRY-related gene. Nature 1994, 372, 525–530. [Google Scholar] [CrossRef]
  36. Li, C.; Scott, D.A.; Hatch, E.; Tian, X.; Mansour, S.L. Dusp6 (Mkp3) is a negative feedback regulator of FGF-stimulated ERK signaling during mouse development. Development 2007, 134, 167–176. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Loeser, R.F.; Chubinskaya, S.; Pacione, C.; Im, H.J. Basic fibroblast growth factor inhibits the anabolic activity of insulin-like growth factor 1 and osteogenic protein 1 in adult human articular chondrocytes. Arthritis Rheum. 2005, 52, 3910–3917. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Kim, B.G.; Kim, H.J.; Park, H.J.; Kim, Y.J.; Yoon, W.J.; Lee, S.J.; Ryoo, H.M.; Cho, J.Y. Runx2 phosphorylation induced by fibroblast growth factor 2/protein kinase C pathways. Proteomics 2006, 6, 1166–1174. [Google Scholar] [CrossRef]
  39. Yan, D.; Chen, D.; Im, H.J. Fibroblast growth factor-2 promotes catabolism via FGFR1-Ras-Raf-MEK1/2-ERK1/2 axis that coordinates with the PKCδ pathway in human articular chondrocytes. J. Cell Biochem. 2012, 113, 2856–2865. [Google Scholar] [CrossRef] [Green Version]
  40. Im, H.J.; Muddasani, P.; Natarajan, V.; Schmid, T.M.; Block, J.A.; Davis, F.; van Wijnen, A.J.; Loeser, R.F. Basic fibroblast growth factor stimulates matrix metalloproteinase-13 via the molecular cross-talk between the mitogen-activated protein kinases and protein kinase Cdelta pathways in human adult articular chondrocytes. J. Biol. Chem. 2007, 282, 11110–11121. [Google Scholar] [CrossRef] [Green Version]
  41. Otsuki, S.; Hanson, S.R.; Miyaki, S.; Grogan, S.P.; Kinoshita, M.; Asahara, H.; Wong, C.H.; Lotz, M.K. Extracellular sulfatases support cartilage homeostasis by regulating BMP and FGF signaling pathways. Proc. Natl. Acad. Sci. USA 2010, 107, 10202–10207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Guo, H.; Ingolia, N.T.; Weissman, J.S.; Bartel, D.P. Mammalian microRNAs predominantly act to decrease target mRNA levels. Nature 2010, 466, 835–840. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Liu, Q.; Halvey, P.J.; Shyr, Y.; Slebos, R.J.; Liebler, D.C.; Zhang, B. Integrative omics analysis reveals the importance and scope of translational repression in microRNA-mediated regulation. Mol. Cell Proteomics 2013, 12, 1900–1911. [Google Scholar] [CrossRef] [Green Version]
  44. Fabian, M.R.; Sonenberg, N.; Filipowicz, W. Regulation of mRNA translation and stability by microRNAs. Annu. Rev. Biochem. 2010, 79, 351–379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Fabian, M.R.; Sonenberg, N. The mechanics of miRNA-mediated gene silencing: A look under the hood of miRISC. Nat. Struct Mol. Biol. 2012, 19, 586–593. [Google Scholar] [CrossRef] [PubMed]
  46. Zeng, Y.; Wagner, E.J.; Cullen, B.R. Both natural and designed micro RNAs can inhibit the expression of cognate mRNAs when expressed in human cells. Mol. Cell 2002, 9, 1327–1333. [Google Scholar] [CrossRef]
  47. Kundu, P.; Fabian, M.R.; Sonenberg, N.; Bhattacharyya, S.N.; Filipowicz, W. HuR protein attenuates miRNA-mediated repression by promoting miRISC dissociation from the target RNA. Nucleic Acids Res. 2012, 40, 5088–5100. [Google Scholar] [CrossRef]
  48. Easow, G.; Teleman, A.A.; Cohen, S.M. Isolation of microRNA targets by miRNP immunopurification. RNA 2007, 13, 1198–1204. [Google Scholar] [CrossRef] [Green Version]
  49. Chi, S.W.; Hannon, G.J.; Darnell, R.B. An alternative mode of microRNA target recognition. Nat. Struct Mol. Biol. 2012, 19, 321–327. [Google Scholar] [CrossRef]
  50. Seok, H.; Ham, J.; Jang, E.S.; Chi, S.W. MicroRNA Target Recognition: Insights from Transcriptome-Wide Non-Canonical Interactions. Mol. Cells 2016, 39, 375–381. [Google Scholar]
  51. Penalva, L.O.; Tenenbaum, S.A.; Keene, J.D. Gene expression analysis of messenger RNP complexes. Methods Mol. Biol. 2004, 257, 125–134. [Google Scholar] [PubMed]
  52. Penalva, L.O.; Burdick, M.D.; Lin, S.M.; Sutterluety, H.; Keene, J.D. RNA-binding proteins to assess gene expression states of co-cultivated cells in response to tumor cells. Mol. Cancer 2004, 3, 24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Lyle, A.N.; Deshpande, N.N.; Taniyama, Y.; Seidel-Rogol, B.; Pounkova, L.; Du, P.; Papaharalambus, C.; Lassègue, B.; Griendling, K.K. Poldip2, a novel regulator of Nox4 and cytoskeletal integrity in vascular smooth muscle cells. Circ. Res. 2009, 105, 249–259. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Grange, L.; Nguyen, M.V.; Lardy, B.; Derouazi, M.; Campion, Y.; Trocme, C.; Paclet, M.H.; Gaudin, P.; Morel, F. NAD(P)H oxidase activity of Nox4 in chondrocytes is both inducible and involved in collagenase expression. Antioxid. Redox Signal. 2006, 8, 1485–1496. [Google Scholar] [CrossRef] [PubMed]
  55. Eyre, D.; Shao, P.; Weis, M.A.; Steinmann, B. The kyphoscoliotic type of Ehlers-Danlos syndrome (type VI): Differential effects on the hydroxylation of lysine in collagens I and II revealed by analysis of cross-linked telopeptides from urine. Mol. Genet. Metab. 2002, 76, 211–216. [Google Scholar] [CrossRef]
  56. Trojanowska, M. Ets factors and regulation of the extracellular matrix. Oncogene 2000, 19, 6464–6471. [Google Scholar] [CrossRef] [Green Version]
  57. Sgouras, D.N.; Athanasiou, M.A.; Beal, G.J.; Fisher, R.J.; Blair, D.G.; Mavrothalassitis, G.J. ERF: An ETS domain protein with strong transcriptional repressor activity, can suppress ets-associated tumorigenesis and is regulated by phosphorylation during cell cycle and mitogenic stimulation. EMBO J. 1995, 14, 4781–4793. [Google Scholar] [CrossRef]
  58. Wong, M.; Kireeva, M.L.; Kolesnikova, T.V.; Lau, L.F. Cyr61, product of a growth factor-inducible immediate-early gene, regulates chondrogenesis in mouse limb bud mesenchymal cells. Dev. Biol. 1997, 192, 492–508. [Google Scholar] [CrossRef] [Green Version]
  59. Elia, L.; Quintavalle, M.; Zhang, J.; Contu, R.; Cossu, L.; Latronico, M.V.; Peterson, K.L.; Indolfi, C.; Catalucci, D.; Chen, J.; et al. The knockout of miR-143 and -145 alters smooth muscle cell maintenance and vascular homeostasis in mice: Correlates with human disease. Cell Death Differ. 2009, 16, 1590–1598. [Google Scholar] [CrossRef]
  60. Koutsis, G.; Siasos, G.; Spengos, K. The emerging role of microRNA in stroke. Curr. Top. Med. Chem. 2013, 13, 1573–1588. [Google Scholar] [CrossRef]
  61. Wei, Y.; Nazari-Jahantigh, M.; Neth, P.; Weber, C.; Schober, A. MicroRNA-126, -145, and -155: A therapeutic triad in atherosclerosis? Arterioscler Thromb Vasc. Biol. 2013, 33, 449–454. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Wang, S.; Bian, C.; Yang, Z.; Bo, Y.; Li, J.; Zeng, L.; Zhou, H.; Zhao, R.C. miR-145 inhibits breast cancer cell growth through RTKN. Int. J. Oncol. 2009, 34, 1461–1466. [Google Scholar] [PubMed]
  63. Wang, F.; Xia, J.; Wang, N.; Zong, H. miR-145 inhibits proliferation and invasion of esophageal squamous cell carcinoma in part by targeting c-Myc. Onkologie 2013, 36, 754–758. [Google Scholar] [CrossRef] [PubMed]
  64. Gregersen, L.H.; Jacobsen, A.B.; Frankel, L.B.; Wen, J.; Krogh, A.; Lund, A.H. MicroRNA-145 targets YES and STAT1 in colon cancer cells. PLoS ONE 2010, 5, e8836. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Dynoodt, P.; Speeckaert, R.; De Wever, O.; Chevolet, I.; Brochez, L.; Lambert, J.; Van Gele, M. miR-145 overexpression suppresses the migration and invasion of metastatic melanoma cells. Int. J. Oncol. 2013, 42, 1443–1451. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Kojima, S.; Enokida, H.; Yoshino, H.; Itesako, T.; Chiyomaru, T.; Kinoshita, T.; Fuse, M.; Nishikawa, R.; Goto, Y.; Naya, Y.; et al. The tumor-suppressive microRNA-143/145 cluster inhibits cell migration and invasion by targeting GOLM1 in prostate cancer. J. Hum. Genet. 2014, 59, 78–87. [Google Scholar] [CrossRef]
  67. Endisha, H.; Rockel, J.; Jurisica, I.; Kapoor, M. The complex landscape of microRNAs in articular cartilage: Biology, pathology, and therapeutic targets. JCI Insight 2018, 3. [Google Scholar] [CrossRef] [Green Version]
  68. Si, H.B.; Yang, T.M.; Li, L.; Tian, M.; Zhou, L.; Li, D.P.; Huang, Q.; Kang, P.D.; Yang, J.; Zhou, Z.K.; et al. miR-140 Attenuates the Progression of Early-Stage Osteoarthritis by Retarding Chondrocyte Senescence. Mol. Ther. Nucleic Acids 2020, 19, 15–30. [Google Scholar] [CrossRef]
  69. Seidl, C.I.; Martinez-Sanchez, A.; Murphy, C.L. Derepression of MicroRNA-138 Contributes to Loss of the Human Articular Chondrocyte Phenotype. Arthritis Rheumatol. 2016, 68, 398–409. [Google Scholar] [CrossRef] [Green Version]
  70. Kapoor, M.; Martel-Pelletier, J.; Lajeunesse, D.; Pelletier, J.P.; Fahmi, H. Role of proinflammatory cytokines in the pathophysiology of osteoarthritis. Nat. Rev. Rheumatol. 2011, 7, 33–42. [Google Scholar] [CrossRef]
  71. Guo, X.Q.; Qi, L.; Yang, J.; Wang, Y.; Wang, C.; Li, Z.M.; Li, L.; Qu, Y.; Wang, D.; Han, Z.M. Salidroside accelerates fracture healing through cell-autonomous and non-autonomous effects on osteoblasts. Cell Tissue Res. 2017, 367, 197–211. [Google Scholar] [CrossRef] [PubMed]
  72. Liu, M.; Zhang, J.; Liu, W.; Wang, W. Salidroside protects ATDC5 cells against lipopolysaccharide-induced injury through up-regulation of microRNA-145 in osteoarthritis. Int. Immunopharmacol 2019, 67, 441–448. [Google Scholar] [CrossRef] [PubMed]
  73. Noguchi, S.; Yamada, N.; Kumazaki, M.; Yasui, Y.; Iwasaki, J.; Naito, S.; Akao, Y. socs7, a target gene of microRNA-145, regulates interferon-β induction through STAT3 nuclear translocation in bladder cancer cells. Cell Death Dis. 2013, 4, e482. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Yi, S.; Lu, D.; Peng, W.; Wang, T.; Zhang, Y.; Lin, H. Differential expression profiling of spleen microRNAs in response to two distinct type II interferons in Tetraodon nigroviridis. PLoS ONE 2014, 9, e96336. [Google Scholar] [CrossRef] [Green Version]
  75. Zhang, F.; Tang, B.; Zhang, Z.; Xu, D.; Ma, G. DUSP6 Inhibitor (E/Z)-BCI Hydrochloride Attenuates Lipopolysaccharide-Induced Inflammatory Responses in Murine Macrophage Cells via Activating the Nrf2 Signaling Axis and Inhibiting the NF-κB Pathway. Inflammation 2019, 42, 672–681. [Google Scholar] [CrossRef]
  76. Bennett, A.M. DUSPs, twists and turns in the Journey to Vascular Inflammation. FEBS J. 2018, 285, 1589–1592. [Google Scholar] [CrossRef] [Green Version]
  77. Bonneau, E.; Neveu, B.; Kostantin, E.; Tsongalis, G.J.; De Guire, V. How close are miRNAs from clinical practice? A perspective on the diagnostic and therapeutic market. EJIFCC 2019, 30, 114–127. [Google Scholar]
  78. Yang, B.; Guo, H.; Zhang, Y.; Chen, L.; Ying, D.; Dong, S. MicroRNA-145 regulates chondrogenic differentiation of mesenchymal stem cells by targeting Sox9. PLoS ONE 2011, 6, e21679. [Google Scholar] [CrossRef]
Figure 1. Immunoblotting analysis of human articular chondrocyte (HAC) lysates following transfection with miR-145 mimics. HACs were transfected with control (C) or miR-145 (145) mimics and immunoprecipitation was performed with an anti-Ago2 antibody (Ago2) or an IgG control, as indicated. I (Input): 2.5% of the total lysate submitted to immunoprecipitation (IP); FT (Flow-through): 5% of the extract recovered after IP. (A) Samples were separated by SDS-PAGE and proven with the anti-Ago2 antibody to determine the presence of Ago proteins. (B) MiR-145 bound to Ago2 following immunoprecipitation, in comparison with IgG, was estimated by RT-qPCR. MiR-145 overexpression resulted in ~40 fold increased in the amount of miR-145 bound to Ago2.
Figure 1. Immunoblotting analysis of human articular chondrocyte (HAC) lysates following transfection with miR-145 mimics. HACs were transfected with control (C) or miR-145 (145) mimics and immunoprecipitation was performed with an anti-Ago2 antibody (Ago2) or an IgG control, as indicated. I (Input): 2.5% of the total lysate submitted to immunoprecipitation (IP); FT (Flow-through): 5% of the extract recovered after IP. (A) Samples were separated by SDS-PAGE and proven with the anti-Ago2 antibody to determine the presence of Ago proteins. (B) MiR-145 bound to Ago2 following immunoprecipitation, in comparison with IgG, was estimated by RT-qPCR. MiR-145 overexpression resulted in ~40 fold increased in the amount of miR-145 bound to Ago2.
Life 10 00058 g001
Figure 2. Enrichment of miR-145 seed sites in the CDSs (A,B) or 3′UTRs (C,D) of the precipitated mRNAs. (A,C) Scatter plot showing the maxima of enrichment profile. Each dot represents a word, the Y-position reflects the maximum score of an enriched word and X-position shows the gene-index where the Z-score is maximum. Top-ranked words with a left-shift associate with gene-expression change (lower gene index represents a stronger upregulation in Ago2 upon miR-145 transfection). Triangles annotate known seed sites of human miRNAs, red triangles show miR-145 seed sites. (B,D) Word-enrichment profile: The line plot shows enrichment through the gene rank for the top ten enriched words. Each line represents the running sum over all scores that quantify a degree of enrichment according to gene upregulation (from most upregulated to most downregulated).
Figure 2. Enrichment of miR-145 seed sites in the CDSs (A,B) or 3′UTRs (C,D) of the precipitated mRNAs. (A,C) Scatter plot showing the maxima of enrichment profile. Each dot represents a word, the Y-position reflects the maximum score of an enriched word and X-position shows the gene-index where the Z-score is maximum. Top-ranked words with a left-shift associate with gene-expression change (lower gene index represents a stronger upregulation in Ago2 upon miR-145 transfection). Triangles annotate known seed sites of human miRNAs, red triangles show miR-145 seed sites. (B,D) Word-enrichment profile: The line plot shows enrichment through the gene rank for the top ten enriched words. Each line represents the running sum over all scores that quantify a degree of enrichment according to gene upregulation (from most upregulated to most downregulated).
Life 10 00058 g002
Figure 3. Luciferase reporter assays to validate a miR-145 direct effect on the expression of selected genes. Hela cells were transfected with luciferase reporters containing full-length CDSs or 3′UTRs, or both, of 12 randomly selected mRNAs from Di-IP, as indicated, downstream the firefly luciferase ORF. Control (white bars) or miR-145 (black bars) mimics were cotransfected in the cells. Values are normalized to the levels of renilla luciferase, independently expressed by the same vector and are shown as relative to that obtained for each construct cotransfected with the control miRNA mimic (± SEM). As positive control, a construct containing three perfectly complementary binding sites for miR-145 (3 x miR-145) was included in the experiment. * p < 0.05; **** p < 0.0001; ns: Not significant.
Figure 3. Luciferase reporter assays to validate a miR-145 direct effect on the expression of selected genes. Hela cells were transfected with luciferase reporters containing full-length CDSs or 3′UTRs, or both, of 12 randomly selected mRNAs from Di-IP, as indicated, downstream the firefly luciferase ORF. Control (white bars) or miR-145 (black bars) mimics were cotransfected in the cells. Values are normalized to the levels of renilla luciferase, independently expressed by the same vector and are shown as relative to that obtained for each construct cotransfected with the control miRNA mimic (± SEM). As positive control, a construct containing three perfectly complementary binding sites for miR-145 (3 x miR-145) was included in the experiment. * p < 0.05; **** p < 0.0001; ns: Not significant.
Life 10 00058 g003
Figure 4. Analysis of the effect of miR-145 overexpression on mRNA levels of selected genes in freshly isolated HACs. (A) Sixteen randomly selected mRNAs from Di-IP and CD2DB1 from both Di-IP and DownT-RNA. (B) Three genes only present in DownT-RNA and COL2A1. Control (C) or miR-145 mimics were transfected in freshly isolated HACs and RNA was extracted 48 h after transfection. Values are presented as relative to that obtained in cells transfected with control mimics for each patient and normalized to RPLP0. Data represents average ± SEM from seven different experiments, each performed with different donor cells (19 yo. female; 8 yo. male; 28 yo. male; 45 yo. male; 16 yo. male; 13 yo. female; 31 yo. female). * p < 0.05; ** p < 0.01; *** p < 0.001; ns: Not significant.
Figure 4. Analysis of the effect of miR-145 overexpression on mRNA levels of selected genes in freshly isolated HACs. (A) Sixteen randomly selected mRNAs from Di-IP and CD2DB1 from both Di-IP and DownT-RNA. (B) Three genes only present in DownT-RNA and COL2A1. Control (C) or miR-145 mimics were transfected in freshly isolated HACs and RNA was extracted 48 h after transfection. Values are presented as relative to that obtained in cells transfected with control mimics for each patient and normalized to RPLP0. Data represents average ± SEM from seven different experiments, each performed with different donor cells (19 yo. female; 8 yo. male; 28 yo. male; 45 yo. male; 16 yo. male; 13 yo. female; 31 yo. female). * p < 0.05; ** p < 0.01; *** p < 0.001; ns: Not significant.
Life 10 00058 g004
Figure 5. MiR-145 directly represses DUSP6. Western blot showing reduced DUSP6 levels following miR-145 overexpression in P2 (A) and freshly isolated (B) HACs. (C) Western blot showing increased DUSP6 levels following miR-145 inhibition in freshly isolated HACs. DUSP6 protein bands were quantified by densitometry and are normalized against tubulin and expressed relative to the control (± SEM). * p < 0.05, ** p < 0.01 (A, right hand side panel, B and C, lower panels). In all experiments HACs were transfected with a relevant control (C) or miR-145 precursor or inhibitor, and subsequently cultured in 20% or, where indicated, 1% O2 tension for 44 h. (D) Hela cells were transfected with luciferase reporters containing three perfectly complementary binding sites for miR-145 (3 x miR-145), the DUSP6 3′UTR containing a putative miR-145 binding site (seed matching nucleotides 1018 to 1025; DUSP6 3′UTR), or a mutated seed-matching site (DUSP6 3′UTR miR-145 MUT) downstream the firefly luciferase open reading frame (ORF). Control (‘C’) or miR-145 mimics were cotransfected in the cells. Values were normalized to the levels of renilla luciferase, independently expressed by the same vector and are shown as relative to that obtained for each construct cotransfected with the control miRNA mimic (± SEM). **** p < 0.0001.
Figure 5. MiR-145 directly represses DUSP6. Western blot showing reduced DUSP6 levels following miR-145 overexpression in P2 (A) and freshly isolated (B) HACs. (C) Western blot showing increased DUSP6 levels following miR-145 inhibition in freshly isolated HACs. DUSP6 protein bands were quantified by densitometry and are normalized against tubulin and expressed relative to the control (± SEM). * p < 0.05, ** p < 0.01 (A, right hand side panel, B and C, lower panels). In all experiments HACs were transfected with a relevant control (C) or miR-145 precursor or inhibitor, and subsequently cultured in 20% or, where indicated, 1% O2 tension for 44 h. (D) Hela cells were transfected with luciferase reporters containing three perfectly complementary binding sites for miR-145 (3 x miR-145), the DUSP6 3′UTR containing a putative miR-145 binding site (seed matching nucleotides 1018 to 1025; DUSP6 3′UTR), or a mutated seed-matching site (DUSP6 3′UTR miR-145 MUT) downstream the firefly luciferase open reading frame (ORF). Control (‘C’) or miR-145 mimics were cotransfected in the cells. Values were normalized to the levels of renilla luciferase, independently expressed by the same vector and are shown as relative to that obtained for each construct cotransfected with the control miRNA mimic (± SEM). **** p < 0.0001.
Life 10 00058 g005
Figure 6. DUSP6 depletion results in MMP13 upregulation. (A) RT-qPCR analysis of DUSP6 and MMP13 levels after transfection of P2 HACs with 5 nM DUSP6 siRNA. Values are presented as relative to that obtained in cells transfected with control siRNA for each patient and normalized to RPLP0. Data represents average ± SEM from four different experiments, each performed with different donor cells (11 yo. male; 47 yo. male; 65 yo. male; 51 yo. male). * p < 0.05; ** p < 0.01; **** p < 0.0001. (B,C) Western blot showing ERK phosphorylation (pERK) after treatment with FGF2 (25 nM). Total ERK (ERK) is shown as a loading control. P1 HACs were transfected with control (C) or miR-145 (145) mimics (B), or control (C) or DUSP6 siRNAs (#1, #2, notice than only #2 siRNA results in DUSP6 depletion, therefore #1 acts as an additional control) (C). After 24 h, HACs were placed in a serum-free medium overnight and stimulated with 25 nM FGF2 for the indicated times. These experiments were performed three times, with three different donors; one representative experiment is shown.
Figure 6. DUSP6 depletion results in MMP13 upregulation. (A) RT-qPCR analysis of DUSP6 and MMP13 levels after transfection of P2 HACs with 5 nM DUSP6 siRNA. Values are presented as relative to that obtained in cells transfected with control siRNA for each patient and normalized to RPLP0. Data represents average ± SEM from four different experiments, each performed with different donor cells (11 yo. male; 47 yo. male; 65 yo. male; 51 yo. male). * p < 0.05; ** p < 0.01; **** p < 0.0001. (B,C) Western blot showing ERK phosphorylation (pERK) after treatment with FGF2 (25 nM). Total ERK (ERK) is shown as a loading control. P1 HACs were transfected with control (C) or miR-145 (145) mimics (B), or control (C) or DUSP6 siRNAs (#1, #2, notice than only #2 siRNA results in DUSP6 depletion, therefore #1 acts as an additional control) (C). After 24 h, HACs were placed in a serum-free medium overnight and stimulated with 25 nM FGF2 for the indicated times. These experiments were performed three times, with three different donors; one representative experiment is shown.
Life 10 00058 g006
Table 1. MicroRNAs (mRNAs) significantly downregulated (False Discovery Rate (FDR) < 0.05) upon miR-145 overexpression (DownT-RNA: Total downregulated mRNA). Control or miR-145 mimics were transfected in P3 HACs from three different patients. After transfection by 24 h, cells were lysed and RNA was extracted and submitted to high-throughput (HT)-seq. **: Genes significantly enriched in Di-IP; *: Genes significantly enriched in Ago2-IP following miR-145 overexpression when the outlier patient was excluded from the analysis (Table S2).
Table 1. MicroRNAs (mRNAs) significantly downregulated (False Discovery Rate (FDR) < 0.05) upon miR-145 overexpression (DownT-RNA: Total downregulated mRNA). Control or miR-145 mimics were transfected in P3 HACs from three different patients. After transfection by 24 h, cells were lysed and RNA was extracted and submitted to high-throughput (HT)-seq. **: Genes significantly enriched in Di-IP; *: Genes significantly enriched in Ago2-IP following miR-145 overexpression when the outlier patient was excluded from the analysis (Table S2).
GeneLogFCFDR
RPL3L−2.311490.018536
RP11-365D23.4−2.161090.02382
FSCN1 *−1.447511.52 × 10−10
HIST1H4L−1.399570.013354
RP11-498C9.15−1.399140.007266
STOX2−1.32540.034993
RP11-649A18.7−1.025320.024884
PKP1−1.000640.032246
ABRACL−0.989044.36 × 10−7
CXCL6−0.954450.012477
RP11-166P13.3−0.937630.024801
CTC-448F2.1−0.880350.000899
MALAT1−0.855540.00094
EID2B−0.847620.02382
DTD1−0.810240.000533
ADPGK−0.783870.002089
CLCN2−0.76620.017386
AC145291.1−0.752860.010585
RP11-620J15.3−0.750570.000152
FAM108C1 *−0.728290.011192
CXCL1−0.727830.010556
TMOD2−0.724330.011264
SERPINE1−0.722020.00137
TMEM9B−0.721660.000757
USP46−0.718850.004563
GCSH−0.709660.001641
NT5DC2−0.708820.001507
GMFB *−0.686070.004577
AL132709.8−0.686060.008869
HOMER2−0.683870.001653
DANCR−0.682070.015835
CTC-459F4.2−0.680750.034195
ZMYND10−0.66220.039097
MTRNR2L10−0.655940.028528
CCDC43 *−0.651810.004142
RPL13AP25−0.651050.039418
AKR1B1−0.650880.001828
SNHG9−0.641410.01704
BAALC−0.641030.011259
UXS1−0.636560.002733
AC005035.1−0.630710.025323
PIGF *−0.61810.045224
VASN−0.617220.000629
RP1-278E11.3−0.608750.004343
AKR1B10−0.608480.049377
ATP13A2−0.603510.002099
SNX24−0.59850.001368
ZNRD1−0.594940.014698
AL132709.5−0.59150.001592
GPRC5C−0.584010.048364
C15orf59−0.582560.014698
GCLM−0.557110.008631
CC2D1B **−0.556880.012839
PCBP4−0.54590.037638
CTB-63M22.1−0.536440.008388
PPP3CA−0.535750.007346
PMEPA1−0.53230.028545
TRIM65−0.529340.04805
FBF1−0.528850.039097
G0S2−0.525820.045454
FAM45A−0.523880.016243
RPL22L1−0.515430.02124
MAP3K11−0.514310.017156
RPS3AP5−0.493640.009962
MT1E−0.488050.02457
PSAT1−0.487160.01541
SMO−0.485850.049462
SWI5−0.482250.028528
H2AFX−0.474180.01541
PAFAH1B2−0.466440.037335
DTYMK−0.465830.042746
RP11-543P15.1−0.46240.027114
IMP3 *−0.4620.04281
RPS28−0.461820.013354
RP11-20O24.4−0.449760.034761
RP1-292B18.1−0.443620.041162
RPS2P55−0.437230.037335
RP11-864N7.2−0.436950.040095
RPL9−0.434890.044893
C7orf73 *−0.424720.042746
Table 2. Significantly enriched mRNAs (FDR < 0.05) in Ago2 immunoprecipitates following miR-145 overexpression (Di-IP: Direct targets (candidates)). Control or miR-145 mimics were transfected in P3 HACs from three different patients. After transfection by 24 h, cells were lysed and submitted to immunoprecipitation with anti-Ago2 antibodies. Immunoprecipitated RNA was extracted and submitted to HT-seq. The list only includes those genes, which expression was downregulated (*) or unchanged as a consequence of miR-145 overexpression. The last three columns indicate the presence of predicted binding sites (+) accordingly to TargetScan (3′ untranslated region (3′UTRs)) and RNAhybrid (3′UTRs or coding sequences (CDSs), as indicated).
Table 2. Significantly enriched mRNAs (FDR < 0.05) in Ago2 immunoprecipitates following miR-145 overexpression (Di-IP: Direct targets (candidates)). Control or miR-145 mimics were transfected in P3 HACs from three different patients. After transfection by 24 h, cells were lysed and submitted to immunoprecipitation with anti-Ago2 antibodies. Immunoprecipitated RNA was extracted and submitted to HT-seq. The list only includes those genes, which expression was downregulated (*) or unchanged as a consequence of miR-145 overexpression. The last three columns indicate the presence of predicted binding sites (+) accordingly to TargetScan (3′ untranslated region (3′UTRs)) and RNAhybrid (3′UTRs or coding sequences (CDSs), as indicated).
GeneLogFCFDRTargetScanRNAhybrid (3’UTR)RNAhybrid (CDS)
COA6 (C1orf31)4.7830814.8 × 10−6 +
EXOG3.418199.44 × 10−5
HIST1H4E3.4046454.49 × 10−19
RP11-539I5.12.9249150.0196
HSBP12.4352443.97 × 10−9+
HIRIP32.3034260.005801 +
DCAF132.2887772.13 × 10−6
BANF12.2503073.08 × 10−12 +
RAB7A2.2237535.22 × 10−5 +
TMEM2232.201130.007146 +
FN3K2.1779530.001289 +
SEC11A2.0400269.75 × 10−10 ++
MOK2.0057430.000539 ++
TALDO11.9890631.19 × 10−9 +
ULBP31.9634320.032833 +
ARHGEF10L1.932252.92 × 10−8 ++
CYR611.8689921.77 × 10−6++
TMOD31.8187920.001511++
GORAB1.8035280.030719++
C6orf1061.7726831.55 × 10−7 +
SUPT16H1.7643850.004327
DECR11.7438220.029152
CC2D1B*1.7435933.69 × 10−6+++
FARSA1.7378665.11 × 10−7 +
TPI1P11.7235065.57 × 10−7
GLT8D11.6988690.000477+
MRPL511.6590221.95 × 10−5 +
ATP6V0B1.6532340.000656++
FAM96A1.6326770.023186
FKBP21.604290.001637 +
SCAND11.5935820.014036
IFFO11.5517480.001511 ++
CDC61.5433870.012934 +
PSMD111.5320490.007204+++
MCM51.5310260.013155 +
MBD21.4955520.013169+ +
MRPL411.4730890.029152 +
TPI11.4515950.005256 ++
DUSP61.4242080.019572++
TIGD51.4192090.007146 ++
RPL391.4111360.025767
OGFR1.3992020.002343 +
RTKN1.3642580.04595+++
CTBP11.3513670.00019 +
ARL21.3328630.002363 +
TADA31.2394230.032833 +
PTGR11.2368090.016147
KLHDC101.1928910.048934+++
GOLM11.1416420.006373+ +
METRNL1.1003530.015776 +
PLOD11.0984260.047772+ +
ADM1.088390.012934
SIX11.0668970.042474++
RCN31.0460730.038849 +
POLDIP21.0417760.015349 +
COX6B11.0326690.030567 +
KIAA09301.0101460.030719+++
ERF1.0065070.047544+++
ATP5A10.9841980.029152 +
Table 3. Enriched Gene Ontology (GO) Terms. Genes altered at the mRNA level upon miR-145 overexpression (FDR < 0.05, Table S1) and direct miR-145 mRNA targets (Table 2) were submitted to the Database for Annotation, Visualization and Integrated Discovery (DAVID v 6.7) for enrichment analysis of biological processes (classification category Panther Biological Processes (Panther_BP_All)).
Table 3. Enriched Gene Ontology (GO) Terms. Genes altered at the mRNA level upon miR-145 overexpression (FDR < 0.05, Table S1) and direct miR-145 mRNA targets (Table 2) were submitted to the Database for Annotation, Visualization and Integrated Discovery (DAVID v 6.7) for enrichment analysis of biological processes (classification category Panther Biological Processes (Panther_BP_All)).
Termp-Value
Interferon-mediated immunity1.10 × 10−13
Immunity and defense3.30 × 10−5
Proteolysis8.80 × 10−3
Glycolysis1.20 × 10−2
T-cell mediated immunity1.30 × 10−2
Macrophage-mediated immunity2.30 × 10−2
MHCI-mediated immunity4.90 × 10−2
Protein metabolism and modification6.50 × 10−2
Protein biosynthesis9.00 × 10−2
AngiogenesisCOA6 (C1orf31)9.60 × 10−2

Share and Cite

MDPI and ACS Style

Martinez-Sanchez, A.; Lazzarano, S.; Sharma, E.; Lockstone, H.; Murphy, C.L. High-Throughput Identification of MiR-145 Targets in Human Articular Chondrocytes. Life 2020, 10, 58. https://doi.org/10.3390/life10050058

AMA Style

Martinez-Sanchez A, Lazzarano S, Sharma E, Lockstone H, Murphy CL. High-Throughput Identification of MiR-145 Targets in Human Articular Chondrocytes. Life. 2020; 10(5):58. https://doi.org/10.3390/life10050058

Chicago/Turabian Style

Martinez-Sanchez, Aida, Stefano Lazzarano, Eshita Sharma, Helen Lockstone, and Christopher L. Murphy. 2020. "High-Throughput Identification of MiR-145 Targets in Human Articular Chondrocytes" Life 10, no. 5: 58. https://doi.org/10.3390/life10050058

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop