Regulation of microRNA-221, -222, -21 and -27 in articular cartilage subjected to abnormal compressive forces.

KEY POINTS: microRNAs (miRs) are small non-coding molecules that regulate post-transcriptional target gene expression. miRs are involved in regulating cellular activities in response to mechanical loading in all physiological systems, although it is largely unknown whether this response differs with increasing magnitudes of load. miR-221, miR-222, miR-21-5p and miR-27a-5p were significantly increased in ex vivo cartilage explants subjected to increasing load magnitude and in in vivo joint cartilage exposed to abnormal loading. TIMP3 and CPEB3 are putative miR targets in chondrocytes Identification of mechanically regulated miRs that have potential to impact on tissue homeostasis provides a mechanism by which load-induced tissue behaviour is regulated, in both health and pathology, in all physiological systems. ABSTRACT: MicroRNAs (miRs) are small non-coding molecules that regulate post-transcriptional target gene expression and are involved in mechano-regulation of cellular activities in all physiological systems. It is unknown whether such epigenetic mechanisms are regulated in response to increasing magnitudes of load. The present study investigated mechano-regulation of miRs in articular cartilage subjected to 'physiological' and 'non-physiological' compressive loads in vitro as a model system and validated findings in an in vivo model of abnormal joint loading. Bovine full-depth articular cartilage explants were loaded to 2.5 MPa (physiological) or 7 MPa (non-physiological) (1 Hz, 15 min) and mechanically-regulated miRs identified using next generation sequencing and verified using a quantitative PCR. Downstream targets were verified using miR-specific mimics or inhibitors in conjunction with 3'-UTR luciferase activity assays. A subset of miRs were mechanically-regulated in ex vivo cartilage explants and in vivo joint cartilage. miR-221, miR-222, miR-21-5p and miR-27a-5p were increased and miR-483 levels decreased with increasing load magnitude. Tissue inhibitor of metalloproteinase 3 (TIMP3) and cytoplasmic polyadenylation element binding protein 3 (CPEB3) were identified as putative downstream targets. Our data confirm miR-221 and -222 mechano-regulation and demonstrates novel mechano-regulation of miR-21-5p and miR-27a-5p in ex vivo and in vivo cartilage loading models. TIMP3 and CPEB3 are putative miR targets in chondrocytes. Identification of specific miRs that are regulated by increasing load magnitude, as well as their potential to impact on tissue homeostasis, has direct relevance to other mechano-sensitive physiological systems and provides a mechanism by which load-induced tissue behaviour is regulated, in both health and pathology.

OATargets: a knowledge base of genes associated with osteoarthritis joint damage in animals.

OBJECTIVES: To collate the genes experimentally modulated in animal models of osteoarthritis (OA) and compare these data with OA transcriptomics data to identify potential therapeutic targets. METHODS: PubMed searches were conducted to identify publications describing gene modulations in animal models. Analysed gene expression data were retrieved from the SkeletalVis database of analysed skeletal microarray and RNA-Seq expression data. A network diffusion approach was used to predict new genes associated with OA joint damage. RESULTS: A total of 459 genes were identified as having been modulated in animal models of OA, with ageing and post-traumatic (surgical) models the most prominent. Ninety-eight of the 143 genes (69%) genetically modulated more than once had a consistent effect on OA joint damage severity. Several discrepancies between different studies were identified, providing lessons on interpretation of these data. We used the data collected along with OA gene expression data to expand existing annotations and prioritise the most promising therapeutic targets, which we validated using the latest reported associations. We constructed an online database OATargets to allow researchers to explore the collated data and integrate it with existing OA and skeletal gene expression data. CONCLUSIONS: We present a comprehensive survey and online resource for understanding gene regulation of animal model OA pathogenesis.

Histone ChIP-Seq identifies differential enhancer usage during chondrogenesis as critical for defining cell-type specificity.

Epigenetic mechanisms are known to regulate gene expression during chondrogenesis. In this study, we have characterized the epigenome during the in vitro differentiation of human mesenchymal stem cells (hMSCs) into chondrocytes. Chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) was used to assess a range of N-terminal posttranscriptional modifications (marks) to histone H3 lysines (H3K4me3, H3K4me1, H3K27ac, H3K27me3, and H3K36me3) in both hMSCs and differentiated chondrocytes. Chromatin states were characterized using histone ChIP-seq and cis-regulatory elements were identified in chondrocytes. Chondrocyte enhancers were associated with chondrogenesis-related gene ontology (GO) terms. In silico analysis and integration of DNA methylation data with chondrogenesis chromatin states revealed that enhancers marked by histone marks H3K4me1 and H3K27ac were de-methylated during in vitro chondrogenesis. Similarity analysis between hMSC and chondrocyte chromatin states defined in this study with epigenomes of cell-types defined by the Roadmap Epigenomics project revealed that enhancers are more distinct between cell-types compared to other chromatin states. Motif analysis revealed that the transcription factor SOX9 is enriched in chondrocyte enhancers. Luciferase reporter assays confirmed that chondrocyte enhancers characterized in this study exhibited enhancer activity which may be modulated by DNA methylation and SOX9 overexpression. Altogether, these integrated data illustrate the cross-talk between different epigenetic mechanisms during chondrocyte differentiation.

The function of microRNAs in cartilage and osteoarthritis.

MicroRNAs are small double-stranded RNAs, which negatively regulate gene expression and have been shown to have key roles in both chondrocyte development and cartilage homeostasis with age. Deletion of all microRNAs in chondrocytes leads to skeletal growth defects in mice, whilst deletion of specific microRNAs, e.g. miR-140, leads to premature articular cartilage degradation and increased susceptibility to posttraumatic osteoarthritis. Studies comparing microRNA expression in normal human articular cartilage compared to osteoarthritic cartilage show differential expression, but varying sample groups make interpretation difficult. MicroRNAs have been proposed as circulating biomarkers of osteoarthritis, but again, this differs amongst patient cohorts. Many micro-RNAs have been shown to have roles in chondrocyte phenotype via signalling pathways, apoptosis, autophagy and senescence. Modulating microRNAs in the joint has been shown to reduce osteoarthritis in animal models and translating this to man as a novel therapeutic strategy will be key.

Epigenetic mechanisms of osteoarthritis risk in human skeletal development.

The epigenome, including the methylation of cytosine bases at CG dinucleotides, is intrinsically linked to transcriptional regulation. The tight regulation of gene expression during skeletal development is essential, with ~1/500 individuals born with skeletal abnormalities. Furthermore, increasing evidence is emerging to link age-associated complex genetic musculoskeletal diseases, including osteoarthritis (OA), to developmental factors including joint shape. Multiple studies have shown a functional role for DNA methylation in the genetic mechanisms of OA risk using articular cartilage samples taken from aged patients. Despite this, our knowledge of temporal changes to the methylome during human cartilage development has been limited. We quantified DNA methylation at ~700,000 individual CpGs across the epigenome of developing human articular cartilage in 72 samples ranging from 7-21 post-conception weeks, a time period that includes cavitation of the developing knee joint. We identified significant changes in 8% of all CpGs, and >9400 developmental differentially methylated regions (dDMRs). The largest hypermethylated dDMRs mapped to transcriptional regulators of early skeletal patterning including MEIS1 and IRX1. Conversely, the largest hypomethylated dDMRs mapped to genes encoding extracellular matrix proteins including SPON2 and TNXB and were enriched in chondrocyte enhancers. Significant correlations were identified between the expression of these genes and methylation within the hypomethylated dDMRs. We further identified 811 CpGs at which significant dimorphism was present between the male and female samples, with the majority (68%) being hypermethylated in female samples. Following imputation, we captured the genotype of these samples at >5 million variants and performed epigenome-wide methylation quantitative trait locus (mQTL) analysis. Colocalization analysis identified 26 loci at which genetic variants exhibited shared impacts upon methylation and OA genetic risk. This included loci which have been previously reported to harbour OA-mQTLs (including GDF5 and ALDH1A2), yet the majority (73%) were novel (including those mapping to CHST3, FGF1 and TEAD1). To our knowledge, this is the first extensive study of DNA methylation across human articular cartilage development. We identify considerable methylomic plasticity within the development of knee cartilage and report active epigenomic mediators of OA risk operating in prenatal joint tissues.

The expression and function of microRNAs in chondrogenesis and osteoarthritis.

OBJECTIVE: To use an in vitro model of chondrogenesis to identify microRNAs (miRNAs) with a functional role in cartilage homeostasis. METHODS: The expression of miRNAs was measured in the ATDC5 cell model of chondrogenesis using microarray and was verified using quantitative reverse transcription-polymerase chain reaction. MicroRNA expression was localized by in situ hybridization. Predicted miRNA target genes were validated using 3'-untranslated region-Luc reporter plasmids containing either wild-type sequences or mutants of the miRNA target sequence. Signaling through the Smad pathway was measured using a (CAGA)(12) -Luc reporter. RESULTS: The expression of several miRNAs was regulated during chondrogenesis. These included 39 miRNAs that are coexpressed with miRNA-140 (miR-140), which is known to be involved in cartilage homeostasis and osteoarthritis (OA). Of these miRNAs, miR-455 resides within an intron of COL27A1 that encodes a cartilage collagen. When human OA cartilage was compared with cartilage obtained from patients with femoral neck fractures, the expression of both miR-140-5p and miR-455-3p was increased in OA cartilage. In situ hybridization showed miR-455-3p expression in the developing limbs of chicks and mice and in human OA cartilage. The expression of miR-455-3p was regulated by transforming growth factor ß (TGFß) ligands, and miRNA regulated TGFß signaling. ACVR2B, SMAD2, and CHRDL1 were direct targets of miR-455-3p and may mediate its functional impact on TGFß signaling. CONCLUSION: MicroRNA-455 is expressed during chondrogenesis and in adult articular cartilage, where it can regulate TGFß signaling, suppressing the Smad2/3 pathway. Diminished signaling through this pathway during the aging process and in OA chondrocytes is known to contribute to cartilage destruction. We propose that the increased expression of miR-455 in OA exacerbates this process and contributes to disease pathology.

Increased serum miR-193a-5p during non-alcoholic fatty liver disease progression: Diagnostic and mechanistic relevance.

BACKGROUND & AIMS: Serum microRNA (miRNA) levels are known to change in non-alcoholic fatty liver disease (NAFLD) and may serve as useful biomarkers. This study aimed to profile miRNAs comprehensively at all NAFLD stages. METHODS: We profiled 2,083 serum miRNAs in a discovery cohort (183 cases with NAFLD representing the complete NAFLD spectrum and 10 population controls). miRNA libraries generated by HTG EdgeSeq were sequenced by Illumina NextSeq. Selected serum miRNAs were profiled in 372 additional cases with NAFLD and 15 population controls by quantitative reverse transcriptase PCR. RESULTS: Levels of 275 miRNAs differed between cases and population controls. Fewer differences were seen within individual NAFLD stages, but miR-193a-5p consistently showed increased levels in all comparisons. Relative to NAFL/non-alcoholic steatohepatitis (NASH) with mild fibrosis (stage 0/1), 3 miRNAs (miR-193a-5p, miR-378d, and miR378d) were increased in cases with NASH and clinically significant fibrosis (stages 2-4), 7 (miR193a-5p, miR-378d, miR-378e, miR-320b, miR-320c, miR-320d, and miR-320e) increased in cases with NAFLD activity score (NAS) 5-8 compared with lower NAS, and 3 (miR-193a-5p, miR-378d, and miR-378e) increased but 1 (miR-19b-3p) decreased in steatosis, activity, and fibrosis (SAF) activity score 2-4 compared with lower SAF activity. The significant findings for miR-193a-5p were replicated in the additional cohort with NAFLD. Studies in Hep G2 cells showed that following palmitic acid treatment, miR-193a-5p expression decreased significantly. Gene targets for miR-193a-5p were investigated in liver RNAseq data for a case subgroup (n = 80); liver GPX8 levels correlated positively with serum miR-193a-5p. CONCLUSIONS: Serum miR-193a-5p levels correlate strongly with NAFLD activity grade and fibrosis stage. MiR-193a-5p may have a role in the hepatic response to oxidative stress and is a potential clinically tractable circulating biomarker for progressive NAFLD. LAY SUMMARY: MicroRNAs (miRNAs) are small pieces of nucleic acid that may turn expression of genes on or off. These molecules can be detected in the blood circulation, and their levels in blood may change in liver disease including non-alcoholic fatty liver disease (NAFLD). To see if we could detect specific miRNA associated with advanced stages of NAFLD, we carried out miRNA sequencing in a group of 183 patients with NAFLD of varying severity together with 10 population controls. We found that a number of miRNAs showed changes, mainly increases, in serum levels but that 1 particular miRNA miR-193a-5p consistently increased. We confirmed this increase in a second group of cases with NAFLD. Measuring this miRNA in a blood sample may be a useful way to determine whether a patient has advanced NAFLD without an invasive liver biopsy.

The role of microRNA-3085 in chondrocyte function.

MicroRNAs have been shown to play a role in cartilage development, homeostasis and breakdown during osteoarthritis. We previously identified miR-3085 in humans as a chondrocyte-selective microRNA, however it could not be detected by Northern blot. The aim of the current study was to prove that miR-3085 is a microRNA and to investigate the function of miR-3085 in signaling pathways relevant to cartilage homeostasis and osteoarthritis. Here, we confirm that miR-3085 is a microRNA and not another class of small RNA using (1) a pre-miR hairpin maturation assay, (2) expression levels in a Dicer null cell line, and (3) Ago2 pulldown. MicroRNA-3085-3p is expressed more highly in micromass than monolayer cultured chondrocytes. Transfection of miR-3085-3p into chondrocytes decreases expression of COL2A1 and ACAN, both of which are validated as direct targets of miR-3085-3p. Interleukin-1 induces the expression of miR-3085-3p, at least in part via NFκB. In a feed-forward mechanism, miR-3085-3p then potentiates NFκB signaling. However, at early time points after transfection, its action appears to be inhibitory. MyD88 has been shown to be a direct target of miR-3085-3p and may be responsible for the early inhibition of NFκB signaling. However, at later time points, MyD88 knockdown remains inhibitory and so other functions of miR-3085-3p are clearly dominant. TGFß1 also induces the expression of miR-3085-3p, but in this instance, it exerts a feedback inhibition on signaling with SMAD3 and SMAD4 shown to be direct targets. This in vitro analysis shows that miR-3085-3p functions in chondrocytes to induce IL-1-signaling, reduce TGFß1 signaling, and inhibit expression of matrix genes. These data suggest that miR-3085-3p has a role in chondrocyte function and could contribute to the process of osteoarthritis.

The microRNA-29 family in cartilage homeostasis and osteoarthritis.

UNLABELLED: MicroRNAs have been shown to function in cartilage development and homeostasis, as well as in progression of osteoarthritis. The objective of the current study was to identify microRNAs involved in the onset or early progression of osteoarthritis and characterise their function in chondrocytes. MicroRNA expression in mouse knee joints post-DMM surgery was measured over 7 days. Expression of miR-29b-3p was increased at day 1 and regulated in the opposite direction to its potential targets. In a mouse model of cartilage injury and in end-stage human OA cartilage, the miR-29 family was also regulated. SOX9 repressed expression of miR-29a-3p and miR-29b-3p via the 29a/b1 promoter. TGFß1 decreased expression of miR-29a, b, and c (3p) in primary chondrocytes, whilst IL-1ß increased (but LPS decreased) their expression. The miR-29 family negatively regulated Smad, NFκB, and canonical WNT signalling pathways. Expression profiles revealed regulation of new WNT-related genes. Amongst these, FZD3, FZD5, DVL3, FRAT2, and CK2A2 were validated as direct targets of the miR-29 family. These data identify the miR-29 family as microRNAs acting across development and progression of OA. They are regulated by factors which are important in OA and impact on relevant signalling pathways. KEY MESSAGES: Expression of the miR-29 family is regulated in cartilage during osteoarthritis. SOX9 represses expression of the miR-29 family in chondrocytes. The miR-29 family is regulated by TGF-ß1 and IL-1 in chondrocytes. The miR-29 family negatively regulates Smad, NFκB, and canonical Wnt signalling. Several Wnt-related genes are direct targets of the miR-29 family.