Annotating Transcriptional Effects of Genetic Variants in Disease-Relevant Tissue: Transcriptome-Wide Allelic Imbalance in Osteoarthritic Cartilage

Osteoarthritis as an Enhanceropathy: Gene Regulation in Complex Musculoskeletal Disease

PURPOSE OF REVIEW

Osteoarthritis is a complex and highly polygenic disease. Over 100 reported osteoarthritis risk variants fall in non-coding regions of the genome, ostensibly conferring functional effects through the disruption of regulatory elements impacting target gene expression. In this review, we summarise the progress that has advanced our knowledge of gene enhancers both within the field of osteoarthritis and more broadly in complex diseases.

RECENT FINDINGS

Advances in technologies such as ATAC-seq have facilitated our understanding of chromatin states in specific cell types, bolstering the interpretation of GWAS and the identification of effector genes. Their application to osteoarthritis research has revealed enhancers as the principal regulatory element driving disease-associated changes in gene expression. However, tissue-specific effects in gene regulatory mechanisms can contribute added complexity to biological interpretation. Understanding gene enhancers and their altered activity in specific cell and tissue types is the key to unlocking the genetic complexity of osteoarthritis. The use of single-cell technologies in osteoarthritis research is still in its infancy. However, such tools offer great promise in improving our functional interpretation of osteoarthritis GWAS and the identification of druggable targets. Large-scale collaborative efforts will be imperative to understand tissue and cell-type specific molecular mechanisms underlying enhancer function in disease.

Risk and Resilience Variants in the Retinoic Acid Metabolic and Developmental Pathways Associated with Risk of FASD Outcomes

Fetal Alcohol Spectrum Disorder (FASD) is a common neurodevelopmental disorder that affects an estimated 2-5% of North Americans. FASD is induced by prenatal alcohol exposure (PAE) during pregnancy and while there is a clear genetic contribution, few genetic factors are currently identified or understood. In this study, using a candidate gene approach, we performed a genetic variant analysis of retinoic acid (RA) metabolic and developmental signaling pathway genes on whole exome sequencing data of 23 FASD-diagnosed individuals. We found risk and resilience alleles in ADH and ALDH genes known to normally be involved in alcohol detoxification at the expense of RA production, causing RA deficiency, following PAE. Risk and resilience variants were also identified in RA-regulated developmental pathway genes, especially in SHH and WNT pathways. Notably, we also identified significant variants in the causative genes of rare neurodevelopmental disorders sharing comorbidities with FASD, including STRA6 (Matthew-Wood), SOX9 (Campomelic Dysplasia), FDG1 (Aarskog), and 22q11.2 deletion syndrome (TBX1). Although this is a small exploratory study, the findings support PAE-induced RA deficiency as a major etiology underlying FASD and suggest risk and resilience variants may be suitable biomarkers to determine the risk of FASD outcomes following PAE.

Thyroid hormone induces ossification and terminal maturation in a preserved OA cartilage biomimetic model

OBJECTIVE

To characterize aspects of triiodothyronine (T3) induced chondrocyte terminal maturation within the molecular osteoarthritis pathophysiology using the previously established T3 human ex vivo osteochondral explant model.

DESIGNS

RNA-sequencing was performed on explant cartilage obtained from OA patients (n = 8), that was cultured ex vivo with or without T3 (10 ng/ml), and main findings were validated using RT-qPCR in an independent sample set (n = 22). Enrichment analysis was used for functional clustering and comparisons with available OA patient RNA-sequencing and GWAS datasets were used to establish relevance for OA pathophysiology by linking to OA patient genomic profiles.

RESULTS

Besides the upregulation of known hypertrophic genes EPAS1 and ANKH, T3 treatment resulted in differential expression of 247 genes with main pathways linked to extracellular matrix and ossification. CCDC80, CDON, ANKH and ATOH8 were among the genes found to consistently mark early, ongoing and terminal maturational OA processes in patients. Furthermore, among the 37 OA risk genes that were significantly affected in cartilage by T3 were COL12A1, TNC, SPARC and PAPPA.

CONCLUSIONS

RNA-sequencing results show that metabolic activation and recuperation of growth plate morphology are induced by T3 in OA chondrocytes, indicating terminal maturation is accelerated. The molecular mechanisms involved in hypertrophy were linked to all stages of OA pathophysiology and will be used to validate disease models for drug testing.

Specific isoforms of the ubiquitin ligase gene WWP2 are targets of osteoarthritis genetic risk via a differentially methylated DNA sequence

BACKGROUND

Transitioning from a genetic association signal to an effector gene and a targetable molecular mechanism requires the application of functional fine-mapping tools such as reporter assays and genome editing. In this report, we undertook such studies on the osteoarthritis (OA) risk that is marked by single nucleotide polymorphism (SNP) rs34195470 (A > G). The OA risk-conferring G allele of this SNP associates with increased DNA methylation (DNAm) at two CpG dinucleotides within WWP2. This gene encodes a ubiquitin ligase and is the host gene of microRNA-140 (miR-140). WWP2 and miR-140 are both regulators of TGFβ signaling.

METHODS

Nucleic acids were extracted from adult OA (arthroplasty) and foetal cartilage. Samples were genotyped and DNAm quantified by pyrosequencing at the two CpGs plus 14 flanking CpGs. CpGs were tested for transcriptional regulatory effects using a chondrocyte cell line and reporter gene assay. DNAm was altered using epigenetic editing, with the impact on gene expression determined using RT-qPCR. In silico analysis complemented laboratory experiments.

RESULTS

rs34195470 genotype associates with differential methylation at 14 of the 16 CpGs in OA cartilage, forming a methylation quantitative trait locus (mQTL). The mQTL is less pronounced in foetal cartilage (5/16 CpGs). The reporter assay revealed that the CpGs reside within a transcriptional regulator. Epigenetic editing to increase their DNAm resulted in altered expression of the full-length and N-terminal transcript isoforms of WWP2. No changes in expression were observed for the C-terminal isoform of WWP2 or for miR-140.

CONCLUSIONS

As far as we are aware, this is the first experimental demonstration of an OA association signal targeting specific transcript isoforms of a gene. The WWP2 isoforms encode proteins with varying substrate specificities for the components of the TGFβ signaling pathway. Future analysis should focus on the substrates regulated by the two WWP2 isoforms that are the targets of this genetic risk.

E3 ubiquitin ligase WWP2 as a promising therapeutic target for diverse human diseases

Mammalian E3 ubiquitin ligases have emerged in recent years as critical regulators of cellular homeostasis due to their roles in targeting substrate proteins for ubiquitination and triggering subsequent downstream signals. In this review, we describe the multiple roles of WWP2, an E3 ubiquitin ligase with unique and important functions in regulating a wide range of biological processes, including DNA repair, gene expression, signal transduction, and cell-fate decisions. As such, WWP2 has evolved to play a key role in normal physiology and diseases, such as tumorigenesis, skeletal development and diseases, immune regulation, cardiovascular disease, and others. We attempt to provide an overview of the biochemical, physiological, and pathophysiological roles of WWP2, as well as open questions for future research, particularly in the context of putative therapeutic opportunities.

The role of WWP1 and WWP2 in bone/cartilage development and diseases

Bone and cartilage diseases are often associated with trauma and senescence, manifested as pain and limited mobility. The repair of bone and cartilage lesion by mesenchymal stem cells is regulated by various transcription factors. WW domain-containing protein 1 (WWP1) and WW domain-containing protein 2 (WWP2) are named for WW domain which recognizes PPXY (phono Ser Pro and Pro Arg) motifs of substrate. WWP1and WWP2 are prominent components of the homologous to the E6-AP carboxyl terminus (HECT) subfamily, a group of the ubiquitin ligase. Recently, some studies have found that WWP1 and WWP2 play an important role in the pathogenesis of bone and cartilage diseases and regulate the level and the transactivation of various transcription factors through ubiquitination. Therefore, this review summarizes the distribution and effects of WWP1 and WWP2 in the development of bone and cartilage, discusses the potential mechanism and therapeutic drugs in bone and cartilage diseases such as osteoarthritis, fracture, and osteoporosis.

Mediation of the Same Epigenetic and Transcriptional Effect by Independent Osteoarthritis Risk-Conferring Alleles on a Shared Target Gene, COLGALT2

OBJECTIVE

Over 100 DNA variants have been associated with osteoarthritis (OA), including rs1046934, located within a linkage disequilibrium block encompassing part of COLGALT2 and TSEN15. The present study was undertaken to determine the target gene(s) and the mechanism of action of the OA locus using human fetal cartilage, cartilage from OA and femoral neck fracture arthroplasty patients, and a chondrocyte cell model.

METHODS

Genotyping and methylation array data of DNA from human OA cartilage samples (n = 87) were used to determine whether the rs1046934 genotype is associated with differential DNA methylation at proximal CpGs. Results were replicated in DNA from human arthroplasty (n = 132) and fetal (n = 77) cartilage samples using pyrosequencing. Allelic expression imbalance (AEI) measured the effects of genotype on COLGALT2 and TSEN15 expression. Reporter gene assays and epigenetic editing determined the functional role of regions harboring differentially methylated CpGs. In silico analyses complemented these experiments.

RESULTS

Three differentially methylated CpGs residing within regulatory regions were detected in the human OA cartilage array data, and 2 of these were replicated in human arthroplasty and fetal cartilage. AEI was detected for COLGALT2 and TSEN15, with associations between expression and methylation for COLGALT2. Reporter gene assays confirmed that the CpGs are in chondrocyte enhancers, with epigenetic editing results directly linking methylation with COLGALT2 expression.

CONCLUSION

COLGALT2 is a target of this OA locus. We previously characterized another OA locus, marked by rs11583641, that independently targets COLGALT2. The genotype of rs1046934, like rs11583641, mediates its effect by modulating expression of COLGALT2 via methylation changes to CpGs located in enhancers. Although the single-nucleotide polymorphisms, CpGs, and enhancers are distinct between the 2 independent OA risk loci, their effect on COLGALT2 is the same. COLGALT2 is the target of independent OA risk loci sharing a common mechanism of action.

WWP2 confers risk to osteoarthritis by affecting cartilage matrix deposition via hypoxia associated genes

OBJECTIVE

To explore the co-expression network of the osteoarthritis (OA) risk gene WWP2 in articular cartilage and study cartilage characteristics when mimicking the effect of OA risk allele rs1052429-A on WWP2 expression in a human 3D in vitro model of cartilage.

METHOD

Co-expression behavior of WWP2 with genes expressed in lesioned OA articular cartilage (N = 35 samples) was explored. By applying lentiviral particle mediated WWP2 upregulation in 3D in vitro pellet cultures of human primary chondrocytes (N = 8 donors) the effects of upregulation on cartilage matrix deposition was evaluated. Finally, we transfected primary chondrocytes with miR-140 mimics to evaluate whether miR-140 and WWP2 are involved in similar pathways.

RESULTS

Upon performing Spearman correlations in lesioned OA cartilage, 98 highly correlating genes (|ρ| > 0.7) were identified. Among these genes, we identified GJA1, GDF10, STC2, WDR1, and WNK4. Subsequent upregulation of WWP2 on 3D chondrocyte pellet cultures resulted in a decreased expression of COL2A1 and ACAN and an increase in EPAS1 expression. Additionally, we observed a decreased expression of GDF10, STC2, and GJA1. Proteomics analysis identified 42 proteins being differentially expressed with WWP2 upregulation, which were enriched for ubiquitin conjugating enzyme activity. Finally, upregulation of miR-140 in 2D chondrocytes resulted in significant upregulation of WWP2 and WDR1.

CONCLUSIONS

Mimicking the effect of OA risk allele rs1052429-A on WWP2 expression initiates detrimental processes in the cartilage shown by a response in hypoxia associated genes EPAS1, GDF10, and GJA1 and a decrease in anabolic markers, COL2A1 and ACAN.

Allelic expression imbalance in articular cartilage and subchondral bone refined genome-wide association signals in osteoarthritis

OBJECTIVES

To present an unbiased approach to identify positional transcript single nucleotide polymorphisms (SNPs) of osteoarthritis (OA) risk loci by allelic expression imbalance (AEI) analyses using RNA sequencing of articular cartilage and subchondral bone from OA patients.

METHODS

RNA sequencing from 65 articular cartilage and 24 subchondral bone from OA patients was used for AEI analysis. AEI was determined for all genes present in the 100 regions reported by the GWAS catalog that were also expressed in cartilage or bone. The count fraction of the alternative allele (φ) was calculated for each heterozygous individual with the risk-SNP or with the SNP in linkage disequilibrium (LD) with it (r2 > 0.6). Furthermore, a meta-analysis was performed to generate a meta-φ (null hypothesis median φ = 0.49) and P-value for each SNP.

RESULTS

We identified 30 transcript SNPs (28 in cartilage and 2 in subchondral bone) subject to AEI in 29 genes. Notably, 10 transcript SNPs were located in genes not previously reported in the GWAS catalog, including two long intergenic non-coding RNAs (lincRNAs), MALAT1 (meta-φ = 0.54, FDR = 1.7x10-4) and ILF3-DT (meta-φ = 0.6, FDR = 1.75x10-5). Moreover, 12 drugs were interacting with 7 genes displaying AEI, of which 7 drugs have been already approved.

CONCLUSIONS

By prioritizing proxy transcript SNPs that mark AEI in cartilage and/or subchondral bone at loci harboring GWAS signals, we present an unbiased approach to identify the most likely functional OA risk-SNP and gene. We identified 10 new potential OA risk genes ready for further, translation towards underlying biological mechanisms.

Developmental principles informing human pluripotent stem cell differentiation to cartilage and bone

Human pluripotent stem cells can differentiate into any cell type given appropriate signals and hence have been used to research early human development of many tissues and diseases. Here, we review the major biological factors that regulate cartilage and bone development through the three main routes of neural crest, lateral plate mesoderm and paraxial mesoderm. We examine how these routes have been used in differentiation protocols that replicate skeletal development using human pluripotent stem cells and how these methods have been refined and improved over time. Finally, we discuss how pluripotent stem cells can be employed to understand human skeletal genetic diseases with a developmental origin and phenotype, and how developmental protocols have been applied to gain a better understanding of these conditions.

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.

Genetic Variants and Protein Alterations of Selenium- and T-2 Toxin-Responsive Genes Are Associated With Chondrocytic Damage in Endemic Osteoarthropathy

The mechanism of environmental factors in Kashin–Beck disease (KBD) remains unknown. We aimed to identify single nucleotide polymorphisms (SNPs) and protein alterations of selenium- and T-2 toxin–responsive genes to provide new evidence of chondrocytic damage in KBD. This study sampled the cubital venous blood of 258 subjects including 129 sex-matched KBD patients and 129 healthy controls for SNP detection. We applied an additive model, a dominant model, and a recessive model to identify significant SNPs. We then used the Comparative Toxicogenomics Database (CTD) to select selenium- and T-2 toxin–responsive genes with the candidate SNP loci. Finally, immunohistochemistry was applied to verify the protein expression of candidate genes in knee cartilage obtained from 15 subjects including 5 KBD, 5 osteoarthritis (OA), and 5 healthy controls. Forty-nine SNPs were genotyped in the current study. The C allele of rs6494629 was less frequent in KBD than in the controls (OR = 0.63, p = 0.011). Based on the CTD database, PPARG, ADAM12, IL6, SMAD3, and TIMP2 were identified to interact with selenium, sodium selenite, and T-2 toxin. KBD was found to be significantly associated with rs12629751 of PPARG (additive model: OR = 0.46, p = 0.012; dominant model: OR = 0.45, p = 0.049; recessive model: OR = 0.18, p = 0.018), rs1871054 of ADAM12 (dominant model: OR = 2.19, p = 0.022), rs1800796 of IL6 (dominant model: OR = 0.30, p = 0.003), rs6494629 of SMAD3 (additive model: OR = 0.65, p = 0.019; dominant model: OR = 0.52, p = 0.012), and rs4789936 of TIMP2 (recessive model: OR = 5.90, p = 0.024). Immunohistochemistry verified significantly upregulated PPARG, ADAM12, SMAD3, and TIMP2 in KBD compared with OA and normal controls (p < 0.05). Genetic polymorphisms of PPARG, ADAM12, SMAD3, and TIMP2 may contribute to the risk of KBD. These genes could promote the pathogenesis of KBD by disturbing ECM homeostasis.

CRLF1 and CLCF1 in Development, Health and Disease

Cytokines and their receptors have a vital function in regulating various processes such as immune function, inflammation, haematopoiesis, cell growth and differentiation. The interaction between a cytokine and its specific receptor triggers intracellular signalling cascades that lead to altered gene expression in the target cell and consequent changes in its proliferation, differentiation, or activation. In this review, we highlight the role of the soluble type I cytokine receptor CRLF1 (cytokine receptor-like factor-1) and the Interleukin (IL)-6 cytokine CLCF1 (cardiotrophin-like cytokine factor 1) during development in physiological and pathological conditions with particular emphasis on Crisponi/cold-induced sweating syndrome (CS/CISS) and discuss new insights, challenges and possibilities arising from recent studies.

Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation.

Elucidating mechano-pathology of osteoarthritis: transcriptome-wide differences in mechanically stressed aged human cartilage explants

Background

Failing of intrinsic chondrocyte repair after mechanical stress is known as one of the most important initiators of osteoarthritis. Nonetheless, insight into these early mechano-pathophysiological processes in age-related human articular cartilage is still lacking. Such insights are needed to advance clinical development. To highlight important molecular processes of osteoarthritis mechano-pathology, the transcriptome-wide changes following injurious mechanical stress on human aged osteochondral explants were characterized.

Methods

Following mechanical stress at a strain of 65% (65%MS) on human osteochondral explants (n_65%MS = 14 versus n_control = 14), RNA sequencing was performed. Differential expression analysis between control and 65%MS was performed to determine mechanical stress-specific changes. Enrichment for pathways and protein-protein interactions was analyzed with Enrichr and STRING.

Results

We identified 156 genes significantly differentially expressed between control and 65%MS human osteochondral explants. Of note, IGFBP5 (FC = 6.01; FDR = 7.81 × 10⁻³) and MMP13 (FC = 5.19; FDR = 4.84 × 10⁻²) were the highest upregulated genes, while IGFBP6 (FC = 0.19; FDR = 3.07 × 10⁻⁴) was the most downregulated gene. Protein-protein interactions were significantly higher than expected by chance (P = 1.44 × 10⁻¹⁵ with connections between 116 out of 156 genes). Pathway analysis showed, among others, enrichment for cellular senescence, insulin-like growth factor (IGF) I and II binding, and focal adhesion.

Conclusions

Our results faithfully represent transcriptomic wide consequences of mechanical stress in human aged articular cartilage with MMP13, IGF binding proteins, and cellular senescence as the most notable results. Acquired knowledge on the as such identified initial, osteoarthritis-related, detrimental responses of chondrocytes may eventually contribute to the development of effective disease-modifying osteoarthritis treatments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13075-021-02595-8.

Characterization of dynamic changes in Matrix Gla Protein (MGP) gene expression as function of genetic risk alleles, osteoarthritis relevant stimuli, and the vitamin K inhibitor warfarin

OBJECTIVE

We here aimed to characterize changes of Matrix Gla Protein (MGP) expression in relation to its recently identified OA risk allele rs1800801-T in OA cartilage, subchondral bone and human ex vivo osteochondral explants subjected to OA related stimuli. Given that MGP function depends on vitamin K bioavailability, we studied the effect of frequently prescribed vitamin K antagonist warfarin.

METHODS

Differential (allelic) mRNA expression of MGP was analyzed using RNA-sequencing data of human OA cartilage and subchondral bone. Human osteochondral explants were used to study exposures to interleukin one beta (IL-1β; inflammation), triiodothyronine (T3; Hypertrophy), warfarin, or 65% mechanical stress (65%MS) as function of rs1800801 genotypes.

RESULTS

We confirmed that the MGP risk allele rs1800801-T was associated with lower expression and that MGP was significantly upregulated in lesioned as compared to preserved OA tissues, mainly in risk allele carriers, in both cartilage and subchondral bone. Moreover, MGP expression was downregulated in response to OA like triggers in cartilage and subchondral bone and this effect might be reduced in carriers of the rs1800801-T risk allele. Finally, warfarin treatment in cartilage increased COL10A1 and reduced SOX9 and MMP3 expression and in subchondral bone reduced COL1A1 and POSTN expression.

DISCUSSION & CONCLUSIONS

Our data highlights that the genetic risk allele lowers MGP expression and upon OA relevant triggers may hamper adequate dynamic changes in MGP expression, mainly in cartilage. The determined direct negative effect of warfarin on human explant cultures functionally underscores the previously found association between vitamin K deficiency and OA.

Association between miRNA Target Sites and Incidence of Primary Osteoarthritis in Women from Volga-Ural Region of Russia: A Case-Control Study

Over the past decades, numerous studies on the genetic markers of osteoarthritis (OA) have been conducted. MiRNA targets sites are a promising new area of research. In this study, we analyzed the polymorphic variants in 3′ UTR regions of COL1A1, COL11A1, ADAMTS5, MMP1, MMP13, SOX9, GDF5, FGF2, FGFR1, and FGFRL1 genes to examine the association between miRNA target site alteration and the incidence of OA in women from the Volga-Ural region of Russia using competitive allele-specific PCR. The T allele of the rs9659030 was associated with generalized OA (OR = 2.0), whereas the C allele of the rs229069 was associated with total OA (OR = 1.43). The T allele of the rs13317 was associated with the total OA (OR = 1.67). After Benjamini-Hochberg correction, only rs13317 remained statistically significant. According to ethnic heterogeneity, associations between the T allele (rs1061237) with OA in women of Russian descent (OR = 1.77), the G allele (rs6854081) in women of Tatar descent (OR = 4.78), the C allele (rs229069) and the T allele (rs73611720) in women of mixed descent and other ethnic groups (OR = 2.25 and OR = 3.02, respectively) were identified. All associations remained statistically significant after Benjamini-Hochberg correction. Together, this study identified miRNA target sites as a genetic marker for the development of OA in various ethnic groups.

Multi-tissue epigenetic analysis of the osteoarthritis susceptibility locus mapping to the plectin gene PLEC

OBJECTIVE

In cartilage, the osteoarthritis (OA) associated single nucleotide polymorphism (SNP) rs11780978 correlates with differential expression of PLEC, and with differential methylation of PLEC CpG dinucleotides, forming eQTLs and mQTLs respectively. This implies that methylation links chondrocyte genotype and phenotype, thus driving the functional effect of this genetic risk signal. PLEC encodes plectin, a cytoskeletal protein that enables tissues to respond to mechanical forces. We sought to assess whether these PLEC functional effects were cartilage specific.

METHOD

Cartilage, fat pad, synovium and peripheral blood were collected from patients undergoing arthroplasty. PLEC CpGs were analysed for mQTLs and allelic expression imbalance (AEI) was performed to test for eQTLs. Plectin was knocked down in a mesenchymal stem cell (MSC) line using CRISPR/Cas9 and cells phenotyped by RNA-sequencing.

RESULTS

mQTLs were discovered in fat pad, synovium and blood. Their effects were however stronger in the joint tissues and of comparable effect between these tissues. We observed AEI in synovium in the same direction as for cartilage and correlations between methylation and PLEC expression. Knocking-down plectin impacted on pathways reported to have a role in OA, including Wnt signalling, glycosaminoglycan biosynthesis and immune regulation.

CONCLUSIONS

Synovium is also a target of the rs11780978 OA association functionally operating on PLEC. In fat pad, mQTLs were identified but these did not correlate with PLEC expression, suggesting the functional effect is not joint-wide. Our study highlights interplay between genetic risk, DNA methylation and gene expression in OA, and reveals clear differences between tissues from the same diseased joint.

Genome-wide association of phenotypes based on clustering patterns of hand osteoarthritis identify WNT9A as novel osteoarthritis gene

Background

Despite recent advances in the understanding of the genetic architecture of osteoarthritis (OA), only two genetic loci have been identified for OA of the hand, in part explained by the complexity of the different hand joints and heterogeneity of OA pathology.

Methods

We used data from the Rotterdam Study (RSI, RSII and RSIII) to create three hand OA phenotypes based on clustering patterns of radiographic OA severity to increase power in our modest discovery genome-wide association studies in the RS (n=8700), and sought replication in an independent cohort, the Framingham Heart Study (n=1203). We used multiple approaches that leverage different levels of information and functional data to further investigate the underlying biological mechanisms and candidate genes for replicated loci. We also attempted to replicate known OA loci at other joint sites, including the hips and knees.

Results

We found two novel genome-wide significant loci for OA in the thumb joints. We identified WNT9A as a possible novel causal gene involved in OA pathogenesis. Furthermore, several previously identified genetic loci for OA seem to confer risk for OA across multiple joints: TGFa, RUNX2, COL27A1, ASTN2, IL11 and GDF5 loci.

Conclusions

We identified a robust novel genetic locus for hand OA on chromosome 1, of which WNT9A is the most likely causal gene. In addition, multiple genetic loci were identified to be associated with OA across multiple joints. Our study confirms the potential for novel insight into the genetic architecture of OA by using biologically meaningful stratified phenotypes.

A single cell transcriptional atlas of early synovial joint development

Synovial joint development begins with the formation of the interzone, a region of condensed mesenchymal cells at the site of the prospective joint. Recently, lineage-tracing strategies have revealed that Gdf5-lineage cells native to and from outside the interzone contribute to most, if not all, of the major joint components. However, there is limited knowledge of the specific transcriptional and signaling programs that regulate interzone formation and fate diversification of synovial joint constituents. To address this, we have performed single cell RNA-Seq analysis of 7329 synovial joint progenitor cells from the developing murine knee joint from E12.5 to E15.5. By using a combination of computational analytics, in situ hybridization and in vitro characterization of prospectively isolated populations, we have identified the transcriptional profiles of the major developmental paths for joint progenitors. Our freely available single cell transcriptional atlas will serve as a resource for the community to uncover transcriptional programs and cell interactions that regulate synovial joint development.

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.

Osteoarthritis year in review 2019: genetics, genomics and epigenetics

Although osteoarthritis (OA) aetiology is complex, genetic, genomic and epigenetic studies published within the last decade have advanced our understanding of the molecular processes underlying this common musculoskeletal disease. The purpose of this narrative review is to highlight the key research articles within the OA genetics, genomics and epigenetics fields that were published between April 2018 and April 2019. The review focuses on the identification of new OA genetic risk loci, genomics techniques that have been used for the first time in human cartilage and new publicly available databases, and datasets that will aid OA functional studies. Fifty-six new OA susceptibility loci were identified by two large scale genome wide association study meta-analyses, increasing the number of genome-wide significant risk loci to 90. OA risk variants are enriched near genes involved in skeletal development and morphology, and show genetic overlap with height, hip shape, bone area and developmental dysplasia of the hip. Several functional studies of OA loci were published, including a genome-wide analysis of genetic variation on cartilage gene expression. A specialised data portal for exploring cross-species skeletal transcriptomic datasets has been developed, and the first use of cartilage single cell RNAseq analysis reported. This year also saw the systematic identification of all microRNAs, long non-coding RNAs and circular RNAs expressed in human OA cartilage. Putative transcriptional regulatory regions have been mapped in human chondrocytes genome-wide, providing a dataset that will facilitate the prioritisation and characterisation of OA genetic and epigenetic loci.

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.