Methylation quantitative trait locus analysis of osteoarthritis links epigenetics with genetic risk

Evolution and advancements in genomics and epigenomics in OA research: How far we have come

OBJECTIVE

Osteoarthritis (OA) is the most prevalent musculoskeletal disease affecting articulating joint tissues, resulting in local and systemic changes that contribute to increased pain and reduced function. Diverse technological advancements have culminated in the advent of high throughput "omic" technologies, enabling identification of comprehensive changes in molecular mediators associated with the disease. Amongst these technologies, genomics and epigenomics - including methylomics and miRNomics, have emerged as important tools to aid our biological understanding of disease.

DESIGN

In this narrative review, we selected articles discussing advancements and applications of these technologies to OA biology and pathology. We discuss how genomics, deoxyribonucleic acid (DNA) methylomics, and miRNomics have uncovered disease-related molecular markers in the local and systemic tissues or fluids of OA patients.

RESULTS

Genomics investigations into the genetic links of OA, including using genome-wide association studies, have evolved to identify 100+ genetic susceptibility markers of OA. Epigenomic investigations of gene methylation status have identified the importance of methylation to OA-related catabolic gene expression. Furthermore, miRNomic studies have identified key microRNA signatures in various tissues and fluids related to OA disease.

CONCLUSIONS

Sharing of standardized, well-annotated omic datasets in curated repositories will be key to enhancing statistical power to detect smaller and targetable changes in the biological signatures underlying OA pathogenesis. Additionally, continued technological developments and analysis methods, including using computational molecular and regulatory networks, are likely to facilitate improved detection of disease-relevant targets, in-turn, supporting precision medicine approaches and new treatment strategies for OA.

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.

Epigenomic profiling of the infrapatellar fat pad in osteoarthritis

Osteoarthritis is a prevalent, complex disease of the joints, and affects multiple intra-articular tissues. Here, we have examined genome-wide DNA methylation profiles of primary infrapatellar fat pad and matched blood samples from 70 osteoarthritis patients undergoing total knee replacement surgery. Comparing the DNA methylation profiles between these tissues reveal widespread epigenetic differences. We produce the first genome-wide methylation quantitative trait locus (mQTL) map of fat pad, and make the resource available to the wider community. Using two-sample Mendelian randomization and colocalization analyses, we resolve osteoarthritis GWAS signals and provide insights into the molecular mechanisms underpinning disease aetiopathology. Our findings provide the first view of the epigenetic landscape of infrapatellar fat pad primary tissue in osteoarthritis.

Osteoarthritis genetic risk acting on the galactosyltransferase gene COLGALT2 has opposing functional effects in articulating joint tissues

BACKGROUND

Investigation of cartilage and chondrocytes has revealed that the osteoarthritis risk marked by the independent DNA variants rs11583641 and rs1046934 mediate  their effects by decreasing the methylation status of CpG dinucleotides in enhancers and increasing the expression of shared target gene COLGALT2. We set out to investigate if these functional effects operate in a non-cartilaginous joint tissue.

METHODS

Nucleic acids were extracted from the synovium of osteoarthritis patients. Samples were genotyped, and DNA methylation was quantified by pyrosequencing at CpGs within the COLGALT2 enhancers. CpGs were tested for enhancer effects using a synovial cell line and a reporter gene assay. DNA methylation was altered using epigenetic editing, with the impact on gene expression determined using quantitative polymerase chain reaction. In silico analysis complemented laboratory experiments.

RESULTS

The rs1046934 genotype did not associate with DNA methylation or COLGALT2 expression in the synovium, whereas the rs11583641 genotype did. Surprisingly, the effects for rs11583641 were opposite to those previously observed in cartilage. Epigenetic editing in synovial cells revealed that enhancer methylation is causally linked to COLGALT2 expression.

CONCLUSIONS

This is the first direct demonstration for osteoarthritis genetic risk of a functional link between DNA methylation and gene expression operating in opposite directions between articular joint tissues. It highlights pleiotropy in the action of osteoarthritis risk and provides a cautionary note in the application of future genetically based osteoarthritis therapies: an intervention that decreases the detrimental effect of a risk allele in one joint tissue may inadvertently increase its detrimental effect in another joint tissue.

Intergenic Interactions of SBNO1, NFAT5 and GLT8D1 Determine the Susceptibility to Knee Osteoarthritis among Europeans of Russia

This study was conducted to examine the associations between genome-wide association studies (GWAS)-important single nucleotide polymorphisms (SNPs) and knee osteoarthritis (KOA) among Europeans of Russia. The present replicative study ("patient-control" design has been used) was carried out on 1000 DNA samples from KOA (n = 500) and KOA-free (n = 500) participants. Ten GWAS-important for KOA SNPs of eight candidate genes (LYPLAL1, GNL3, GLT8D1, SBNO1, WWP2, NFAT5, TGFA, GDF5) were studied. To assess the link between SNPs and KOA susceptibility, logistic regression (to establish independent SNP effects) and MB-MDR (to identify SNP-SNP interactions) were used. As a result of this genetic analysis, the associations of individual SNPs with KOA have not been proven. Eight loci out of ten tested SNPs interacted with each other (within twelve genetic models) and determined susceptibility to KOA. The greatest contribution to the disease development were made by three polymorphisms/genes such as rs6976 (C>T) GLT8D1, rs56116847 (G>A) SBNO1, rs6499244 (T>A) NFAT5 (each was included in 2/3 [8 out 12] KOA-responsible genetic interaction models). A two-locus epistatic interaction of rs56116847 (G >A) SBNO1 × rs6499244 (T>A) NFAT5 determined the maximum percentage (0.86%) of KOA entropy. KOA-associated SNPs are regulatory polymorphisms that affect the expression/splicing level, epigenetic modification of 72 genes in KOA-pathogenetically significant organs such as skeletal muscles, tibial arteries/nerves, thyroid, adipose tissue, etc. These putative KOA-effector genes are mainly involved in the organization/activity of the exoribonuclease complex and antigen processing/presentation pathways. In conclusion, KOA susceptibility among Europeans of Russia is mediated by intergenic interactions (but not the main effects) of GWAS-important SNPs.

The epigenetic players and the chromatin marks involved in the articular cartilage during osteoarthritis

Epigenetics defines the modifications of the genome that do not involve a change in the nucleotide sequence of DNA. These modifications constitute a mechanism of gene regulation poorly explored in the context of cartilage physiology. They are now intensively studied by the scientific community working on articular cartilage and its related pathology such as osteoarthritis. Indeed, epigenetic regulations can control the expression of crucial gene in the chondrocytes, the only resident cells of cartilage. Some epigenetic changes are considered as a possible cause of the abnormal gene expression and the subsequent alteration of the chondrocyte phenotype (hypertrophy, proliferation, senescence…) as observed in osteoarthritic cartilage. Osteoarthritis is a joint pathology, which results in impaired extracellular matrix homeostasis and leads ultimately to the progressive destruction of cartilage. To date, there is no pharmacological treatment and the exact causes have yet to be defined. Given that the epigenetic modifying enzymes can be controlled by pharmacological inhibitors, it is thus crucial to describe the epigenetic marks that enable the normal expression of extracellular matrix encoding genes, and those associated with the abnormal gene expression such as degradative enzyme or inflammatory cytokines encoding genes. In this review, only the DNA methylation and histone modifications will be detailed with regard to normal and osteoarthritic cartilage. Although frequently referred as epigenetic mechanisms, the regulatory mechanisms involving microRNAs will not be discussed. Altogether, this review will show how this nascent field influences our understanding of the pathogenesis of OA in terms of diagnosis and how controlling the epigenetic marks can help defining epigenetic therapies.

Epigenetic Regulation in Knee Osteoarthritis

Osteoarthritis (OA) is a complicated disease with both hereditary and environmental causes. Despite an increase in reports of possible OA risk loci, it has become clear that genetics is not the sole cause of osteoarthritis. Epigenetics, which can be triggered by environmental influences and result in transcriptional alterations, may have a role in OA pathogenesis. The majority of recent research on the epigenetics of OA has been focused on DNA methylation, histone modification, and non-coding RNAs. However, this study will explore epigenetic regulation in OA at the present stage. How genetics, environmental variables, and epigenetics interact will be researched, shedding light for future studies. Their possible interaction and control processes open up new avenues for the development of innovative osteoarthritis treatment and diagnostic techniques.

An epigenome-wide view of osteoarthritis in primary tissues

Osteoarthritis is a complex degenerative joint disease. Here, we investigate matched genotype and methylation profiles of primary chondrocytes from macroscopically intact (low-grade) and degraded (high-grade) osteoarthritis cartilage and from synoviocytes collected from 98 osteoarthritis-affected individuals undergoing knee replacement surgery. We perform an epigenome-wide association study of knee cartilage degeneration and report robustly replicating methylation markers, which reveal an etiologic mechanism linked to the migration of epithelial cells. Using machine learning, we derive methylation models of cartilage degeneration, which we validate with 82% accuracy in independent data. We report a genome-wide methylation quantitative trait locus (mQTL) map of articular cartilage and synovium and identify 18 disease-grade-specific mQTLs in osteoarthritis cartilage. We resolve osteoarthritis GWAS loci through causal inference and colocalization analyses and decipher the epigenetic mechanisms that mediate the effect of genotype on disease risk. Together, our findings provide enhanced insights into epigenetic mechanisms underlying osteoarthritis in primary tissues.

Multi-omics molecular biomarkers and database of osteoarthritis

Osteoarthritis (OA) is the most common form of arthritis in the adult population and is a leading cause of disability. OA-related genetic loci may play an important role in clinical diagnosis and disease progression. With the rapid development of diverse technologies and omics methods, many OA-related public data sets have been accumulated. Here, we retrieved a diverse set of omics experimental results from 159 publications, including genome-wide association study, differentially expressed genes and differential methylation regions, and 2405 classified OA-related gene markers. Meanwhile, based on recent single-cell RNA-seq data from different joints, 5459 cell-type gene markers of joints were collected. The information has been integrated into an online database named OAomics and molecular biomarkers (OAOB). The database (http://ibi.zju.edu.cn/oaobdb/) provides a web server for OA marker genes, omics features and so on. To our knowledge, this is the first database of molecular biomarkers for OA.

Genetics of osteoarthritis

Osteoarthritis genetics has been transformed in the past decade through the application of large-scale genome-wide association scans. So far, over 100 polymorphic DNA variants have been associated with this common and complex disease. These genetic risk variants account for over 20% of osteoarthritis heritability and the vast majority map to non-protein coding regions of the genome where they are presumed to act by regulating the expression of target genes. Statistical fine mapping, in silico analyses of genomics data, and laboratory-based functional studies have enabled the identification of some of these targets, which encode proteins with diverse roles, including extracellular signaling molecules, intracellular enzymes, transcription factors, and cytoskeletal proteins. A large number of the risk variants correlate with epigenetic factors, in particular cartilage DNA methylation changes in cis, implying that epigenetics may be a conduit through which genetic effects on gene expression are mediated. Some of the variants also appear to have been selected as humans adapted to bipedalism, suggesting that a proportion of osteoarthritis genetic susceptibility results from antagonistic pleiotropy, with risk variants having a positive role in joint formation but a negative role in the long-term health of the joint. Although data from an osteoarthritis genetic study has not yet directly led to a novel treatment, some of the osteoarthritis associated genes code for proteins that have available therapeutics. Genetic investigations are therefore revealing fascinating fundamental insights into osteoarthritis and can expose options for translational intervention.

Epigenetic Regulation of Chondrocytes and Subchondral Bone in Osteoarthritis

The aim of this review is to provide an updated review of the epigenetic factors involved in the onset and development of osteoarthritis (OA). OA is a prevalent degenerative joint disease characterized by chronic inflammation, ectopic bone formation within the joint, and physical and proteolytic cartilage degradation which result in chronic pain and loss of mobility. At present, no disease-modifying therapeutics exist for the prevention or treatment of the disease. Research has identified several OA risk factors including mechanical stressors, physical activity, obesity, traumatic joint injury, genetic predisposition, and age. Recently, there has been increased interest in identifying epigenetic factors involved in the pathogenesis of OA. In this review, we detail several of these epigenetic modifications with known functions in the onset and progression of the disease. We also review current therapeutics targeting aberrant epigenetic regulation as potential options for preventive or therapeutic treatment.

An eQTL variant of ALDH1A2 is associated with Kashin-Beck disease in Chinese population

INTRODUCTION

The aims of the study were to investigate the relationship between aldehyde dehydrogenase 1 family member A2 (ALDH1A2) and Kashin-Beck disease (KBD), explore the effects of the rs3204689 polymorphism and methylation status on the expression levels of ALDH1A2, and further clarify the pathogenesis of KBD.

MATERIALS AND METHODS

The genotype of ALDH1A2 rs3204689 was detected by PCR-RFLP in 103 KBD patients and 109 healthy controls in the whole blood. The mRNA level of ALDH1A2 was measured by qRT-PCR, and the protein expression was detected using IHC staining and Western blotting. The MSP-PCR was used to identify the ALDH1A2 methylation level.

RESULTS

There were significant differences in G/G, G/C, and C/C frequencies of ALDH1A2 rs3204689 between the KBD and control groups (χ² = 7.113, P = 0.029); the minor allele G of ALDH1A2 was associated with the risk of KBD (χ² = 5.984, P = 0.014). The mRNA and protein levels of ALDH1A2 were increased in the whole blood and cartilage of KBD patients compared with the controls (P = 0.049, P < 0.0001, P = 0.019). Meanwhile, a statistically significant difference was observed between G/G, G/C and C/C genotype on mRNA expression (P = 0.039). The methylation level of the ALDH1A2 gene promoter region showed no significant difference between the KBD and control groups (χ² = 0.317, P = 0.573).

CONCLUSION

Our case-control study indicates that the common variant rs3204689 near ALDH1A2 is associated with KBD in Chinese population. The risk allele G of rs3204689 is statistically linked to the high expression of ALDH1A2, which is up-regulated in the cartilage and whole blood of KBD patients. Our findings suggest a potential role of ALDH1A2 in the pathogenesis of KBD.

Associations between adipokines gene polymorphisms and knee osteoarthritis: a meta-analysis

Background

Adipokines gene polymorphisms are speculated to be associated with the risk of knee osteoarthritis (OA), but evidence remains conflicting. This study therefore aimed to examine whether associations exist between adipokines gene polymorphisms and knee OA by considering the evidence collected from eligible studies through a meta-analysis.

Methods

A systematic search was performed on PubMed, Embase, Web of Science, China National Knowledge Infrastructure (CNKI), and Wanfang up to March 31, 2020. Meta-analysis was carried out by focusing on the associations between adipokines gene polymorphisms and knee OA with the allele model, dominant model, and recessive model.

Results

The present meta-analysis included 5 eligible studies for ADIPOQ rs1501299 with 1,021 cases and 1,097 controls, 3 eligible studies for ADIPOQ rs2241766 with 549 cases and 544 controls, 3 eligible studies for LEPR rs1137101 with 808 cases and 856 controls, 2 eligible studies for VISFATIN rs4730153 with 339 cases and 680 controls and 2 eligible studies for VISFATIN rs16872158 with 339 cases and 680 controls. Significant association was observed between LEPR rs1137101 and knee OA in the overall population (recessive: OR = 0.40, 95% CI 0.21–0.79). Limited data revealed that associations may exist between ADIPOQ rs2241766 and knee OA in Asians (dominant: OR = 1.35, 95% CI 1.03–1.78), between VISFATIN rs4730153 and knee OA in Asians (allele: OR = 0.58, 95% CI 0.41–0.83; dominant: OR = 0.57, 95% CI 0.39–0.83), and between VISFATIN rs16872158 and knee OA in Asians (allele: OR = 1.84, 95% CI 1.26–2.68; dominant: OR = 1.94, 95% CI 1.31–2.89).

Conclusions

Adipokines gene polymorphisms may be associated with knee OA. The association was observed in LEPR rs1137101 in the present study. In addition, limited data revealed that associations may also exist in ADIPOQ rs2241766, VISFATIN rs4730153 and VISFATIN rs16872158.

Prospero registration

CRD42020187664.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12891-022-05111-4.

Methylation quantitative trait locus analysis of chronic postsurgical pain uncovers epigenetic mediators of genetic risk

Background: Overlap of pathways enriched by single nucleotide polymorphisms and DNA-methylation underlying chronic postsurgical pain (CPSP), prompted pilot study of CPSP-associated methylation quantitative trait loci (meQTL). Materials & methods: Children undergoing spine-fusion were recruited prospectively. Logistic-regression for genome- and epigenome-wide CPSP association and DNA-methylation-single nucleotide polymorphism association/mediation analyses to identify meQTLs were followed by functional genomics analyses. Results: CPSP (n = 20/58) and non-CPSP groups differed in pain-measures. Of 2753 meQTLs, DNA-methylation at 127 cytosine-guanine dinucleotides mediated association of 470 meQTLs with CPSP (p < 0.05). At PARK16 locus, CPSP risk meQTLs were associated with decreased DNA-methylation at RAB7L1 and increased DNA-methylation at PM20D1. Corresponding RAB7L1/PM20D1 blood eQTLs (GTEx) and cytosine-guanine dinucleotide-loci enrichment for histone marks, transcription factor binding sites and ATAC-seq peaks suggest altered transcription factor-binding. Conclusion: CPSP-associated meQTLs indicate epigenetic mechanisms mediate genetic risk. Clinical trial registration: NCT01839461, NCT01731873 (ClinicalTrials.gov).

Genetic and Epigenetic Interplay Within a COLGALT2 Enhancer Associated With Osteoarthritis

OBJECTIVE

The osteoarthritis (OA)-associated single-nucleotide polymorphism (SNP) rs11583641 is located in COLGALT2, encoding a posttranslational modifier of collagen. In cartilage, the SNP genotype correlates with DNA methylation in a putative enhancer. This study was undertaken to characterize the mechanistic relationship between rs11583641, the putative enhancer, and COLGALT2 expression using cartilage samples from human patients and a chondrocyte cell model.

METHODS

Nucleic acids were extracted from articular cartilage samples obtained from patients with OA (n = 137). Samples were genotyped, and DNA methylation was quantified at 12 CpGs using pyrosequencing. The putative enhancer was deleted in Tc28a2 chondrocytes using clustered regularly interspaced short palindromic repeat/Cas9, and the impact on nearby gene expression was determined using real-time quantitative polymerase chain reaction. Targeted modulation of the epigenome using catalytically dead Cas9 (dCas9) constructs fused to DNA methyltransferase 3a or ten-eleven translocase 1 allowed for the investigation of a causal relationship between DNA methylation and enhancer activity.

RESULTS

The genotype at rs11583641 correlated with DNA methylation at 3 CpGs, and the presence of the OA risk allele, C, corresponded to reduced levels of methylation. Deletion of the enhancer resulted in a 2.7-fold reduction in COLGALT2 expression. Targeted methylation and demethylation of the CpGs had antagonistic effects on COLGALT2 expression. An allelic imbalance in the expression of COLGALT2 was identified in the cartilage from patients with OA, with relative overexpression of the OA risk allele. Allelic expression ratios correlated with DNA methylation at 4 CpGs.

CONCLUSION

COLGALT2 is a target of OA genetic risk at this locus. The genotype at rs11583641 impacts DNA methylation in a gene enhancer, which, in turn, modulates COLGALT2 expression. COLGALT2 encodes an enzyme that initiates posttranslational glycosylation of collagens and is therefore a compelling OA susceptibility target.

Genetic and Epigenetic Fine-Tuning of TGFB1 Expression Within the Human Osteoarthritic Joint

OBJECTIVE

Osteoarthritis (OA) is an age-related disease characterized by articular cartilage degeneration. It is largely heritable, and genetic screening has identified single-nucleotide polymorphisms (SNPs) marking genomic risk loci. One such locus is marked by the G>A SNP rs75621460, downstream of TGFB1. This gene encodes transforming growth factor β1, the correct expression of which is essential for cartilage maintenance. This study investigated the regulatory activity of rs75621460 to characterize its impact on TGFB1 expression in disease-relevant patient samples (n = 319) and in Tc28a2 immortalized chondrocytes.

METHODS

Articular cartilage samples from human patients were genotyped, and DNA methylation levels were quantified using pyrosequencing. Gene reporter and electrophoretic mobility shift assays were used to determine differential nuclear protein binding to the region. The functional impact of DNA methylation on TGFB1 expression was tested using targeted epigenome editing.

RESULTS

The analyses showed that SNP rs75621460 was located within a TGFB1 enhancer region, and the OA risk allele A altered transcription factor binding, with decreased enhancer activity. Protein complexes binding to A (but not G) induced DNA methylation at flanking CG dinucleotides. Strong correlations between patient DNA methylation levels and TGFB1 expression were observed, with directly opposing effects in the cartilage and the synovium at this locus. This demonstrated biologic pleiotropy in the impact of the SNP within different tissues of the articulating joint.

CONCLUSION

The OA risk SNP rs75621460 impacts TGFB1 expression by modulating the function of a gene enhancer. We propose a mechanism by which the SNP impacts enhancer function, providing novel biologic insight into one mechanism of OA genetic risk, which may facilitate the development of future pharmacologic therapies.

Interplay between genetics and epigenetics in osteoarthritis

Research into the molecular genetics of osteoarthritis (OA) has been substantially bolstered in the past few years by the implementation of powerful genome-wide scans that have revealed a large number of novel risk loci associated with the disease. This refreshing wave of discovery has occurred concurrently with epigenetic studies of joint tissues that have examined DNA methylation, histone modifications and regulatory RNAs. These epigenetic analyses have involved investigations of joint development, homeostasis and disease and have used both human samples and animal models. What has become apparent from a comparison of these two complementary approaches is that many OA genetic risk signals interact with, map to or correlate with epigenetic mediators. This discovery implies that epigenetic mechanisms, and their effect on gene expression, are a major conduit through which OA genetic risk polymorphisms exert their functional effects. This observation is particularly exciting as it provides mechanistic insight into OA susceptibility. Furthermore, this knowledge reveals avenues for attenuating the negative effect of risk-conferring alleles by exposing the epigenome as an exploitable target for therapeutic intervention in OA.

Increased WISP1 expression in human osteoarthritic articular cartilage is epigenetically regulated and decreases cartilage matrix production

OBJECTIVES

Previously, we have shown the involvement of Wnt-activated protein Wnt-1-induced signaling protein 1 (WISP1) in the development of OA in mice. Here, we aimed to characterize the relation between WISP1 expression and human OA and its regulatory epigenetic determinants.

METHODS

Preserved and lesioned articular cartilage from end-stage OA patients and non-OA-diagnosed individuals was collected. WISP1 expression was determined using immunohistochemistry and damage was classified using Mankin scoring. RNA expression and DNA methylation were assessed in silico from genome-wide datasets (microarray analysis and RNA sequencing, and 450 k-methylationarrays, respectively). Effects of WISP1 were tested in pellet cultures of primary human chondrocytes.

RESULTS

WISP1 expression in cartilage of OA patients was increased compared with non-OA-diagnosed controls and, within OA patients, WISP1 was even higher in lesioned compared with preserved regions, with expression strongly correlating with Mankin score. In early symptomatic OA patients with disease progression, higher synovial WISP1 expression was observed as compared with non-progressors. Notably, increased WISP1 expression was inversely correlated with methylation levels of a positional CpG-dinucleotide (cg10191240), with lesioned areas showing strong hypomethylation for this CpG as compared with preserved cartilage. Additionally, we observed that methylation levels were allele-dependent for an intronic single-nucleotide polymorphism nearby cg10191240. Finally, addition of recombinant WISP1 to pellets of primary chondrocytes strongly inhibited deposition of extracellular matrix as reflected by decreased pellet circumference, proteoglycan content and decreased expression of matrix components.

CONCLUSION

Increased WISP1 expression is found in lesioned human articular cartilage, and appears epigenetically regulated via DNA methylation. In vitro assays suggest that increased WISP1 is detrimental for cartilage integrity.

Lymphocyte DNA methylation mediates genetic risk at shared immune-mediated disease loci

Background

Defining regulatory mechanisms through which noncoding risk variants influence the cell-mediated pathogenesis of immune-mediated disease (IMD) has emerged as a priority in the post–genome-wide association study era.

Objectives

With a focus on rheumatoid arthritis, we sought new insight into genetic mechanisms of adaptive immune dysregulation to help prioritize molecular pathways for targeting in this and related immune pathologies.

Methods

Whole-genome methylation and transcriptional data from isolated CD4⁺ T cells and B cells of more than 100 genotyped and phenotyped patients with inflammatory arthritis, all of whom were naive to immunomodulatory treatments, were obtained. Analysis integrated these comprehensive data with genome-wide association study findings across IMDs and other publicly available resources.

Results

We provide strong evidence that disease-associated DNA variants regulate cis-CpG methylation in CD4⁺ T and/or B cells at 37% RA loci. Using paired, cell-specific transcriptomic data and causal inference testing, we identify examples where site-specific DNA methylation in turn mediates gene expression, including FCRL3 in both cell types and ORMDL3/GSDMB, IL6ST/ANKRD55, and JAZF1 in CD4⁺ T cells. A number of genes regulated in this way highlight mechanisms common to RA and other IMDs including multiple sclerosis and asthma, in turn distinguishing them from osteoarthritis, a primarily degenerative disease. Finally, we corroborate the observed effects experimentally.

Conclusions

Our observations highlight important mechanisms of genetic risk in RA and the wider context of immune dysregulation. They confirm the utility of DNA methylation profiling as a tool for causal gene prioritization and, potentially, therapeutic targeting in complex IMD.

Genome-wide analysis of aberrant methylation of enhancer DNA in human osteoarthritis

Background

Osteoarthritis is a chronic musculoskeletal disease characterized by age-related gradual thinning and a high risk in females. Recent studies have shown that DNA methylation plays important roles in osteoarthritis. However, the genome-wide pattern of methylation in enhancers in osteoarthritis remains unclear.

Methods

To explore the function of enhancers in osteoarthritis, we quantified CpG methylation in human enhancers based on a public dataset that included methylation profiles of 470,870 CpG probes in 108 samples from patients with hip and knee osteoarthritis and hip tissues from healthy individuals. Combining various bioinformatics analysis tools, we systematically analyzed aberrant DNA methylation of the enhancers throughout the genome in knee osteoarthritis and hip osteoarthritis.

Results

We identified 16,816 differentially methylated CpGs, and nearly half (8111) of them were from enhancers, suggesting major DNA methylation changes in both types of osteoarthritis in the enhancer regions. A detailed analysis of hip osteoarthritis identified 2426 differentially methylated CpGs in enhancers between male and female patients, and 84.5% of them were hypomethylated in female patients and enriched in phenotypes related to hip osteoarthritis in females. Next, we explored the enhancer methylation dynamics among patients with knee osteoarthritis and identified 280 differentially methylated enhancer CpGs that were enriched in the human phenotypes and disease ontologies related to osteoarthritis. Finally, a comparison of enhancer methylation between knee osteoarthritis and hip osteoarthritis revealed organ source-dependent differences in enhancer methylation.

Conclusion

Our findings indicate that aberrant methylation of enhancers is related to osteoarthritis phenotypes, and a comprehensive atlas of enhancer methylation is useful for further analysis of the epigenetic regulation of osteoarthritis and the development of clinical drugs for treatment of osteoarthritis.

Runx2 plays a central role in Osteoarthritis development

Osteoarthritis (OA) is the most common form of arthritis, is the leading cause of impaired mobility in the elderly, and accounts for more than a third of chronic moderate to severe pain. As a degenerative joint disorder, OA affects the whole joint and results in synovial hyperplasia, degradation of articular cartilage, subchondral sclerosis, osteophyte formation, and chronic pain. Currently, there is no effective drug to decelerate OA progression and molecular targets for drug development have been insufficiently investigated. Anti-OA drug development can benefit from more and precise knowledge of molecular targets for drug development. Runt-related transcription factor 2 (Runx2) is a key transcription factor controlling osteoblast and chondrocyte differentiation and is among the most promising potential therapeutic targets. Notably, Runx2 expression is upregulated in several murine OA models, suggesting a role in disease pathogenesis. In this review article, we summarized recent findings on Runx2 related to OA development and evaluated its potential as a therapeutic target.

The translational potential of this article

A better understanding of the role of Runx2 in osteoarthritis pathogenesis will contribute to the development of novel intervention of osteoarthritis disease.

Discovery and analysis of methylation quantitative trait loci (mQTLs) mapping to novel osteoarthritis genetic risk signals

OBJECTIVE

Osteoarthritis (OA) is polygenic with over 90 independent genome-wide association loci so far reported. A key next step is the identification of target genes and the molecular mechanisms through which this genetic risk operates. The majority of OA risk-conferring alleles are predicted to act by modulating gene expression. DNA methylation at CpG dinucleotides may be a functional conduit through which this occurs and is detectable by mapping methylation quantitative trait loci, or mQTLs. This approach can therefore provide functional insight into OA risk and will prioritize genes for subsequent investigation. That was our goal, with a focus on the largest set of OA loci yet to be reported.

METHOD

We investigated DNA methylation, genotype and RNA sequencing data derived from the cartilage of patients who had undergone arthroplasty and combined this with in silico analyses of expression quantitative trait loci, epigenomes and chromatin interactions.

RESULTS

We investigated 42 OA risk loci and in ten of these we identified 24 CpGs in which methylation correlated with genotype (false discovery rate (FDR) P-values ranging from 0.049 to 1.73x10^-25). In silico analyses of these mQTLs prioritised genes and regulatory elements at the majority of the ten loci, with COLGALT2 (encoding a collagen galactosyltransferase), COL11A2 (encoding a polypeptide chain of type XI collagen) and WWP2 (encoding a ubiquitin ligase active during chondrogenesis) emerging as particularly compelling target genes.

CONCLUSION

We have highlighted the pivotal role of DNA methylation as a link between genetic risk and OA and prioritized genes for further investigation.

Prioritization of PLEC and GRINA as Osteoarthritis Risk Genes Through the Identification and Characterization of Novel Methylation Quantitative Trait Loci

OBJECTIVE

To identify methylation quantitative trait loci (mQTLs) correlating with osteoarthritis (OA) risk alleles and to undertake mechanistic characterization as a means of target gene prioritization.

METHODS

We used genome-wide genotyping and cartilage DNA methylation array data in a discovery screen of novel OA risk loci. This was followed by methylation, gene expression analysis, and genotyping studies in additional cartilage samples, accompanied by in silico analyses.

RESULTS

We identified 4 novel OA mQTLs. The most significant mQTL contained 9 CpG sites where methylation correlated with OA risk genotype, with 5 of the CpG sites having P values <1 × 10^-10 . The 9 CpG sites reside in an interval of only 7.7 kb within the PLEC gene and form 2 distinct clusters. We were able to prioritize PLEC and the adjacent gene GRINA as independent targets of the OA risk. We identified PLEC and GRINA expression QTLs operating in cartilage, as well as methylation-expression QTLs operating on the 2 genes. GRINA and PLEC also demonstrated differential expression between OA hip and non-OA hip cartilage.

CONCLUSION

PLEC encodes plectin, a cytoskeletal protein that maintains tissue integrity by regulating intracellular signaling in response to mechanical stimuli. GRINA encodes the ionotropic glutamate receptor TMBIM3 (transmembrane BAX inhibitor 1 motif-containing protein family member 3), which regulates cell survival. Based on our results, we hypothesize that in a joint predisposed to OA, expression of these genes alters in order to combat aberrant biomechanics, and that this is epigenetically regulated. However, carriage of the OA risk-conferring allele at this locus hinders this response and contributes to disease development.

Multi-omics integrative analysis identified SNP-methylation-mRNA: Interaction in peripheral blood mononuclear cells

Genetic variants have potential influence on DNA methylation and thereby regulate mRNA expression. This study aimed to comprehensively reveal the relationships among SNP, methylation and mRNA, and identify methylation-mediated regulation patterns in human peripheral blood mononuclear cells (PBMCs). Based on in-house multi-omics datasets from 43 Chinese Han female subjects, genome-wide association trios were constructed by simultaneously testing the following three association pairs: SNP-methylation, methylation-mRNA and SNP-mRNA. Causal inference test (CIT) was used to identify methylation-mediated genetic effects on mRNA. A total of 64,184 significant cis-methylation quantitative trait loci (meQTLs) were identified (FDR < 0.05). Among the 745 constructed trios, 464 trios formed SNP-methylation-mRNA regulation chains (CIT). Network analysis (Cytoscape 3.3.0) constructed multiple complex regulation networks among SNP, methylation and mRNA (eg a total of 43 SNPs simultaneously connected to cg22517527 and further to PRMT2, DIP2A and YBEY). The regulation chains were supported by the evidence from 4DGenome database, relevant to immune or inflammatory related diseases/traits, and overlapped with previous eQTLs from dbGaP and GTEx. The results provide new insights into the regulation patterns among SNP, DNA methylation and mRNA expression, especially for the methylation-mediated effects, and also increase our understanding of functional mechanisms underlying the established associations.

Identification of a novel, methylation-dependent, RUNX2 regulatory region associated with osteoarthritis risk

Osteoarthritis (OA) is a common, multifactorial and polygenic skeletal disease that, in its severest form, requires joint replacement surgery to restore mobility and to relieve chronic pain. Using tissues from the articulating joints of 260 patients with OA and a range of in vitro experiments, including CRISPR-Cas9, we have characterized an intergenic regulatory element. Here, genotype at an OA risk locus correlates with differential DNA methylation, with altered gene expression of both a transcriptional regulator (RUNX2), and a chromatin remodelling protein (SUPT3H). RUNX2 is a strong candidate for OA susceptibility, with its encoded protein being essential for skeletogenesis and healthy joint function. The OA risk locus includes single nucleotide polymorphisms (SNPs) located within and flanking the differentially methylated region (DMR). The OA association SNP, rs10948172, demonstrates particularly strong correlation with methylation, and two intergenic SNPs falling within the DMR (rs62435998 and rs62435999) demonstrate genetic and epigenetic effects on the regulatory activity of this region. We therefore posit that the OA signal mediates its effect by modulating the methylation of the regulatory element, which then impacts on gene expression, with RUNX2 being the principal target. Our study highlights the interplay between DNA methylation, OA genetic risk and the downstream regulation of genes critical to normal joint function.

Role of Epigenomics in Bone and Cartilage Disease

Phenotypic variation in skeletal traits and diseases is the product of genetic and environmental factors. Epigenetic mechanisms include information-containing factors, other than DNA sequence, that cause stable changes in gene expression and are maintained during cell divisions. They represent a link between environmental influences, genome features, and the resulting phenotype. The main epigenetic factors are DNA methylation, posttranslational changes of histones, and higher-order chromatin structure. Sometimes non-coding RNAs, such as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), are also included in the broad term of epigenetic factors. There is rapidly expanding experimental evidence for a role of epigenetic factors in the differentiation of bone cells and the pathogenesis of skeletal disorders, such as osteoporosis and osteoarthritis. However, different from genetic factors, epigenetic signatures are cell- and tissue-specific and can change with time. Thus, elucidating their role has particular difficulties, especially in human studies. Nevertheless, epigenomewide association studies are beginning to disclose some disease-specific patterns that help to understand skeletal cell biology and may lead to development of new epigenetic-based biomarkers, as well as new drug targets useful for treating diffuse and localized disorders. Here we provide an overview and update of recent advances on the role of epigenomics in bone and cartilage diseases. © 2019 American Society for Bone and Mineral Research.

Involvement of epigenetics in osteoarthritis

Osteoarthritis (OA) is the most prevalent chronic age-related arthritic disease that mainly affects the diarthrodial joints. Nevertheless, there is no treatment currently available that can effectively reduce symptoms or slow down or stop disease progression. The lack of disease-modifying therapies could be explained by the complex pathogenesis of OA, which is still not completely understood. Intertwined epigenetic mechanisms such as DNA methylation, histone modifications, and noncoding RNAs (ncRNAs) have been indicated as important cellular tools to maintain tissue homeostasis upon environmental challenges. The current review illustrates that dysfunctional epigenetic control mechanisms in the articular cartilage likely play an important role in driving OA pathophysiology.

Functional Characterization of the Osteoarthritis Genetic Risk Residing at ALDH1A2 Identifies rs12915901 as a Key Target Variant

OBJECTIVE

To identify the functional single-nucleotide polymorphisms (SNPs) and mechanisms conferring increased risk of hand osteoarthritis (OA) at the ALDH1A2 locus, which is a retinoic acid regulatory gene.

METHODS

Tissue samples from 247 patients with knee, hip, or hand OA who had undergone joint surgery were included. RNA-sequencing analysis was used to investigate differential expression of ALDH1A2 and other retinoic acid signaling pathway genes in cartilage. Expression of ALDH1A2 in joint tissues obtained from multiple sites was quantified using quantitative reverse transcription-polymerase chain reaction. Allelic expression imbalance (AEI) was measured by pyrosequencing. The consequences of ALDH1A2 depletion by RNA interference were assessed in primary human chondrocytes. In silico and in vitro analyses were used to pinpoint which, among 62 highly correlated SNPs, could account for the association at the locus.

RESULTS

ALDH1A2 expression was observed across multiple joint tissue samples, including osteochondral tissue from the hand. The expression of ALDH1A2 and of several retinoic acid signaling genes was different in diseased cartilage compared to non-diseased cartilage, with ALDH1A2 showing lower levels in OA cartilage. Experimental depletion of ALDH1A2 resulted in changes in the expression levels of a number of chondrogenic markers, including SOX9. In addition, reduced expression of the OA risk-conferring allele was witnessed in a number of joint tissues, with the strongest effect in cartilage. The intronic SNP rs12915901 recapitulated the AEI observed in patient tissues, while the Ets transcription factors were identified as potential mediators of this effect.

CONCLUSION

The ALDH1A2 locus seems to increase the risk of hand OA through decreased expression of ALDH1A2 in joint tissues, with the effect dependent on rs12915901. These findings indicate a mechanism that may now be targeted to modulate OA risk.

Breast cancer family history and allele-specific DNA methylation in the legacy girls study

Family history, a well-established risk factor for breast cancer, can have both genetic and environmental contributions. Shared environment in families as well as epigenetic changes that also may be influenced by shared genetics and environment may also explain familial clustering of cancers. Epigenetic regulation, such as DNA methylation, can change the activity of a DNA segment without a change in the sequence; environmental exposures experienced across the life course can induce such changes. However, genetic-epigenetic interactions, detected as methylation quantitative trait loci (mQTLs; a.k.a. meQTLs) and haplotype-dependent allele-specific methylation (hap-ASM), can also contribute to inter-individual differences in DNA methylation patterns. To identify differentially methylated regions (DMRs) associated with breast cancer susceptibility, we examined differences in white blood cell DNA methylation in 29 candidate genes in 426 girls (ages 6-13 years) from the LEGACY Girls Study, 239 with and 187 without a breast cancer family history (BCFH). We measured methylation by targeted massively parallel bisulfite sequencing (bis-seq) and observed BCFH DMRs in two genes: ESR1 (Δ4.9%, P = 0.003) and SEC16B (Δ3.6%, P = 0.026), each of which has been previously implicated in breast cancer susceptibility and pubertal development. These DMRs showed high inter-individual variability in methylation, suggesting the presence of mQTLs/hap-ASM. Using single nucleotide polymorphisms data in the bis-seq amplicon, we found strong hap-ASM in SEC16B (with allele specific-differences ranging from 42% to 74%). These findings suggest that differential methylation in genes relevant to breast cancer susceptibility may be present early in life, and that inherited genetic factors underlie some of these epigenetic differences.

The Epigenomic Landscape in Osteoarthritis

PURPOSE OF REVIEW

Epigenomics has emerged as a key player in our rapidly evolving understanding of osteoarthritis. Historical studies implicated epigenetic alterations, particularly DNA methylation, in OA pathogenesis; however, recent technological advances have resulted in numerous epigenome-wide studies examining in detail epigenetic modifications in OA. The purpose of this article is to introduce basic concepts in epigenetics and their recent applications to the study of osteoarthritis development and progression.

RECENT FINDINGS

Epigenetics describes three major phenomena: DNA modification via methylation, histone sidechain modifications, and short noncoding RNA sequences which work in concert to regulate gene transcription in a heritable fashion. Cartilage has been the most widely studied tissue in OA, and differential methylation of genes involved in inflammation, cell cycle, TGFβ, and HOX genes have been confirmed several times. Bone studies suggest similar findings, and the intriguing possibility of epigenetic changes in subchondral bone during many OA processes. Multiple studies have demonstrated the involvement of certain noncoding RNAs, particularly miR-140, in OA development via modulation of key catabolic factors. Although much work has been done, much is still unknown. Future epigenomic studies will no doubt continue to widen our understanding of extraarticular tissues and OA pathogenesis, and studies in animal models may offer glimpses into epigenome alterations in the earliest stages of OA.

Genetic polymorphism in the asporin gene is not a key risk factor for osteoarthritis: Evidence based on an updated cumulative meta-analysis

To provide an evidence-based medical basis for the treatment of osteoarthritis, a meta-analysis was performed to assess the association between asporin (ASPN) gene polymorphism and susceptibility to osteoarthritis (OA). The current study searched the literature from January 1st, 1915 through February 1st, 2017 using the Cochrane Library, PubMed, the Excerpta Medica database (EMBASE) and three main Chinese databases (VIP, CNKI and Wan Fang). Cohort and case-control studies that explored the association between different types of ASPN alleles and OA susceptibility were evaluated. The K/L grading system, clinical and radiological diagnoses were used for OA diagnosis. A random-effects model was used in a pooled analysis to adjust for heterogeneity of the included studies, and the differences between treatment groups were reported as odds ratio (OR), 95% confidence intervals (CIs) and P-values. Begg's funnel plots and Egger's tests were used to assess publication bias in the present meta-analysis. Following document retrieval and screening, a total of 10 studies were deemed eligible, including 4,842 patients and 3,661 healthy subjects. Results of the multivariate meta-regression analysis revealed that the study sample size was a source of heterogeneity between studies. The D17 allele was a risk factor for the development of OA (OR=1.33, 95% CI: 1.02–1.73, P<0.05). The other alleles were not considered as risk factors for development of OA (P>0.05). The results of the meta-analysis verified that ASPN polymorphisms were not significantly relevant to an increased OA risk. However, the mechanisms contributing to the association between ASPN polymorphisms and OA risk still require further study.

MeQTL analysis of childhood obesity links epigenetics with a risk SNP rs17782313 near MC4R from meta-analysis

Earlier GWAS has identified that rs17782313 near MC4R was associated with obesity. However, subsequent studies showed conflicting results, especially among childhood. Besides, the mechanisms underlying the association between rs17782313 and childhood obesity remain largely unexplored, and genetic and epigenetic may interact and together affect the development of childhood obesity. We conducted a comprehensive meta-analysis to assess the association between rs17782313 and childhood obesity. MeQTL and eQTL analysis was applied to explore the effect of rs17782313 on DNA methylation and MC4R expression. We found that rs17782313 near MC4R was associated with increased childhood obesity risk and BMI z-score in several inheritable models (P < 0.05). Additionally, the similar trend was observed among subgroups of Asians, Caucasian. Furthermore, meQTL and eQTL analysis indicated that individuals carrying rs17782313 TT genotype were significantly associated with increased methylation level of cg10097150 located in MC4R promoter and decreased expression of MC4R than those with CT/CC genotype (P = 1.7 × 10⁻⁴ and P = 1.9 × 10⁻³ respectively). Our results strongly confirmed that rs17782313 was associated with increased risk of childhood obesity. Furthermore, rs17782313 T allele was correlated with promoter hypermethylation and decreased expression of MC4R, thus involved in the development of childhood obesity.

Comprehensive characterization, annotation and innovative use of Infinium DNA methylation BeadChip probes

Illumina Infinium DNA Methylation BeadChips represent the most widely used genome-scale DNA methylation assays. Existing strategies for masking Infinium probes overlapping repeats or single nucleotide polymorphisms (SNPs) are based largely on ad hoc assumptions and subjective criteria. In addition, the recently introduced MethylationEPIC (EPIC) array expands on the utility of this platform, but has not yet been well characterized. We present in this paper an extensive characterization of probes on the EPIC and HM450 microarrays, including mappability to the latest genome build, genomic copy number of the 3΄ nested subsequence and influence of polymorphisms including a previously unrecognized color channel switch for Type I probes. We show empirical evidence for exclusion criteria for underperforming probes, providing a sounder basis than current ad hoc criteria for exclusion. In addition, we describe novel probe uses, exemplified by the addition of a total of 1052 SNP probes to the existing 59 explicit SNP probes on the EPIC array and the use of these probes to predict ethnicity. Finally, we present an innovative out-of-band color channel application for the dual use of 62 371 probes as internal bisulfite conversion controls.

Integrative epigenomics, transcriptomics and proteomics of patient chondrocytes reveal genes and pathways involved in osteoarthritis

Osteoarthritis (OA) is a common disease characterized by cartilage degeneration and joint remodeling. The underlying molecular changes underpinning disease progression are incompletely understood. We investigated genes and pathways that mark OA progression in isolated primary chondrocytes taken from paired intact versus degraded articular cartilage samples across 38 patients undergoing joint replacement surgery (discovery cohort: 12 knee OA, replication cohorts: 17 knee OA, 9 hip OA patients). We combined genome-wide DNA methylation, RNA sequencing, and quantitative proteomics data. We identified 49 genes differentially regulated between intact and degraded cartilage in at least two -omics levels, 16 of which have not previously been implicated in OA progression. Integrated pathway analysis implicated the involvement of extracellular matrix degradation, collagen catabolism and angiogenesis in disease progression. Using independent replication datasets, we showed that the direction of change is consistent for over 90% of differentially expressed genes and differentially methylated CpG probes. AQP1, COL1A1 and CLEC3B were significantly differentially regulated across all three -omics levels, confirming their differential expression in human disease. Through integration of genome-wide methylation, gene and protein expression data in human primary chondrocytes, we identified consistent molecular players in OA progression that replicated across independent datasets and that have translational potential.

DNA methyltransferase 3b regulates articular cartilage homeostasis by altering metabolism

Osteoarthritis (OA) is the most common form of arthritis worldwide. It is a complex disease affecting the whole joint but is generally characterized by progressive degradation of articular cartilage. Recent genome-wide association screens have implicated distinct DNA methylation signatures in OA patients. We show that the de novo DNA methyltransferase (Dnmt) 3b, but not Dnmt3a, is present in healthy murine and human articular chondrocytes and its expression decreases in OA mouse models and in chondrocytes from human OA patients. Targeted deletion of Dnmt3b in murine articular chondrocytes results in an early-onset and progressive postnatal OA-like pathology. RNA-Seq and methylC-Seq analyses of Dnmt3b loss-of-function chondrocytes show that cellular metabolic processes are affected. Specifically, TCA metabolites and mitochondrial respiration are elevated. Importantly, a chondroprotective effect was found following Dnmt3b gain of function in murine articular chondrocytes in vitro and in vivo. This study shows that Dnmt3b plays a significant role in regulating postnatal articular cartilage homeostasis. Cellular pathways regulated by Dnmt3b in chondrocytes may provide novel targets for therapeutic approaches to treat OA.

Osteoarthritis year in review 2016: genetics, genomics and epigenetics

The purpose of this narrative review is to provide an overview of last year's publications in the field of genetics, genomics and epigenetics in the osteoarthritis (OA) field. Major themes arising from a Pubmed search on (epi)genetics in OA were identified. In addition, general developments in the fast evolving field of (epi)genetics are reviewed and relevance for the OA field is summarized. In the last 5 years, a number of genome-wide association studies have identified a modest number of genetic loci associated to OA. Continued functional research into these DNA variants is showing putative biological mechanisms underlying these associations. Over the last year, no additional large genome-wide association studies were published, but there clearly remains much to be discovered in the OA genetic field. A lot of research has been done into the epigenetics of OA over the last year. Several genome-wide screens examining the methylome of osteoarthritic cartilage were done. Pathway analysis confirmed deregulation of developmental and extracellular pathways in OA cartilage. Over the last year many microRNAs (miRNAs) have been identified that potentially play important roles in cartilage homeostasis and/or OA process. Continued research will learn whether these identified miRNAs are truly causal and can be used in clinical applications. Many of the epigenetic findings need further confirmation, but they highlight potential novel pathways involved in cartilage biology and OA.

The first international workshop on the epigenetics of osteoarthritis

Osteoarthritis (OA) is a major clinical problem across the world, in part due to the lack of disease-modifying drugs resulting, to a significant degree, from our incomplete understanding of the underlying molecular mechanisms of the disease. Emerging evidence points to a role of epigenetics in the pathogenesis of OA, but research in this area is still in its early stages. In order to summarize current knowledge and to facilitate the potential coordination of future research activities, the first international workshop on the epigenetics of OA was held in Amsterdam in October 2015. Recent findings on DNA methylation and hydroxymethylation, histone modifications, noncoding RNAs, and other epigenetic mechanisms were presented and discussed. The workshop demonstrated the advantage of bringing together those working in this nascent field and highlights from the event are summarized in this report in the form of summaries from invited speakers and organizers.

The impact of methylation quantitative trait loci (mQTLs) on active smoking-related DNA methylation changes

BACKGROUND

Methylation quantitative trait loci (mQTLs) are the genetic variants that may affect the DNA methylation patterns of CpG sites. However, their roles in influencing the disturbances of smoking-related epigenetic changes have not been well established. This study was conducted to address whether mQTLs exist in the vicinity of smoking-related CpG sites (± 50 kb) and to examine their associations with smoking exposure and all-cause mortality in older adults.

RESULTS

We obtained DNA methylation profiles in whole blood samples by Illumina Infinium Human Methylation 450 BeadChip array of two independent subsamples of the ESTHER study (discovery set, n = 581; validation set, n = 368) and their corresponding genotyping data using the Illumina Infinium OncoArray BeadChip. After correction for multiple testing (FDR), we successfully identified that 70 out of 151 previously reported smoking-related CpG sites were significantly associated with 192 SNPs within the 50 kb search window of each locus. The 192 mQTLs significantly influenced the active smoking-related DNA methylation changes, with percentage changes ranging from 0.01 to 18.96%, especially for the weakly/moderately smoking-related CpG sites. However, these identified mQTLs were not directly associated with active smoking exposure or all-cause mortality.

CONCLUSIONS

Our findings clearly demonstrated that if not dealt with properly, the mQTLs might impair the power of epigenetic-based models of smoking exposure to a certain extent. In addition, such genetic variants could be the key factor to distinguish between the heritable and smoking-induced impact on epigenome disparities. These mQTLs are of special importance when DNA methylation markers measured by Illumina Infinium assay are used for any comparative population studies related to smoking-related cancers and chronic diseases.

Functional genomics in osteoarthritis: Past, present, and future

Osteoarthritis (OA) is a common complex disease of high public health burden. OA is characterized by the degeneration of affected joints leading to pain and reduced mobility. Over the last few years, several studies have focused on the genomic changes underpinning OA. Here, we provide a comprehensive overview of genome-wide, non-hypothesis-driven functional genomics (methylation, gene, and protein expression) studies of knee and hip OA in humans. Individual studies have generally been limited in sample size and hence power, and have differed in their approaches; nonetheless, some common themes have started to emerge, notably the role played by biological processes related to the extracellular matrix, immune response, the WNT pathway, angiogenesis, and skeletal development. Larger-scale studies and streamlined, robust methodologies will be needed to further elucidate the biological etiology of OA going forward. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 34:1105-1110, 2016.

Analysis of genetics and DNA methylation in osteoarthritis: What have we learnt about the disease?

Osteoarthritis (OA) is a chronic musculoskeletal disease characterised by the destruction of articular cartilage, synovial inflammation and bone remodelling. Disease aetiology is complex and highly heritable, with genetic variation estimated to contribute to 50% of OA occurrence. Epigenetic alterations, including DNA methylation changes, have also been implicated in OA pathophysiology. This review examines what genetic and DNA methylation studies have taught us about the genes and pathways involved in OA pathology. The influence of DNA methylation on the molecular mechanisms underlying OA genetic risk and the consequence of this interaction on disease susceptibility and penetrance are also discussed.