PINK1 and Parkin target Miro for phosphorylation and degradation to arrest mitochondrial motility.

Cells keep their energy balance and avoid oxidative stress by regulating mitochondrial movement, distribution, and clearance. We report here that two Parkinson's disease proteins, the Ser/Thr kinase PINK1 and ubiquitin ligase Parkin, participate in this regulation by arresting mitochondrial movement. PINK1 phosphorylates Miro, a component of the primary motor/adaptor complex that anchors kinesin to the mitochondrial surface. The phosphorylation of Miro activates proteasomal degradation of Miro in a Parkin-dependent manner. Removal of Miro from the mitochondrion also detaches kinesin from its surface. By preventing mitochondrial movement, the PINK1/Parkin pathway may quarantine damaged mitochondria prior to their clearance. PINK1 has been shown to act upstream of Parkin, but the mechanism corresponding to this relationship has not been known. We propose that PINK1 phosphorylation of substrates triggers the subsequent action of Parkin and the proteasome.

Genetic risk of osteoarthritis operates during human skeletogenesis.

Osteoarthritis (OA) is a polygenic disease of older people resulting in the breakdown of cartilage within articular joints. Although it is a leading cause of disability, there are no disease-modifying therapies. Evidence is emerging to support the origins of OA in skeletogenesis. Whereas methylation quantitative trait loci (mQTLs) co-localizing with OA genome-wide association study signals have been identified in aged human cartilage and used to identify effector genes and variants, such analyses have never been conducted during human development. Here, for the first time, we have investigated the developmental origins of OA genetic risk at seven well-characterized OA risk loci, comprising 39 OA-mQTL CpGs, in human fetal limb (FL) and cartilage (FC) tissues using a range of molecular genetic techniques. We identified significant OA-mQTLs at 14 and 29 CpGs in FL and FC tissues, respectively, and compared our results with aged cartilage samples (AC). Differential methylation was observed at 26 sites between FC and AC, with the majority becoming actively hypermethylated in old age. Notably, 6/9 OA effector genes showed allelic expression imbalances during fetal development. Finally, we conducted ATAC-sequencing in cartilage from the developing and aged hip and knee to identify accessible chromatin regions and found enrichment for transcription factor binding motifs including SOX9 and FOS/JUN. For the first time, we have demonstrated the activity of OA-mQTLs and expression imbalance of OA effector genes during human skeletogenesis. We show striking differences in the spatiotemporal function of these loci, contributing to our understanding of OA aetiology, with implications for the timing and strategy of pharmacological interventions.

Progressive liver, kidney, and heart degeneration in children and adults affected by TULP3 mutations.

Organ fibrosis is a shared endpoint of many diseases, yet underlying mechanisms are not well understood. Several pathways governed by the primary cilium, a sensory antenna present on most vertebrate cells, have been linked with fibrosis. Ciliopathies usually start early in life and represent a considerable disease burden. We performed massively parallel sequencing by using cohorts of genetically unsolved individuals with unexplained liver and kidney failure and correlated this with clinical, imaging, and histopathological analyses. Mechanistic studies were conducted with a vertebrate model and primary cells. We detected bi-allelic deleterious variants in TULP3, encoding a critical adaptor protein for ciliary trafficking, in a total of 15 mostly adult individuals, originating from eight unrelated families, with progressive degenerative liver fibrosis, fibrocystic kidney disease, and hypertrophic cardiomyopathy with atypical fibrotic patterns on histopathology. We recapitulated the human phenotype in adult zebrafish and confirmed disruption of critical ciliary cargo composition in several primary cell lines derived from affected individuals. Further, we show interaction between TULP3 and the nuclear deacetylase SIRT1, with roles in DNA damage repair and fibrosis, and report increased DNA damage ex vivo. Transcriptomic studies demonstrated upregulation of profibrotic pathways with gene clusters for hypertrophic cardiomyopathy and WNT and TGF-ß signaling. These findings identify variants in TULP3 as a monogenic cause for progressive degenerative disease of major organs in which affected individuals benefit from early detection and improved clinical management. Elucidation of mechanisms crucial for DNA damage repair and tissue maintenance will guide novel therapeutic avenues for this and similar genetic and non-genomic diseases.

Three decades of advancements in osteoarthritis research: insights from transcriptomic, proteomic, and metabolomic studies.

OBJECTIVE: Osteoarthritis (OA) is a complex disease involving contributions from both local joint tissues and systemic sources. Patient characteristics, encompassing sociodemographic and clinical variables, are intricately linked with OA rendering its understanding challenging. Technological advancements have allowed for a comprehensive analysis of transcripts, proteomes and metabolomes in OA tissues/fluids through omic analyses. The objective of this review is to highlight the advancements achieved by omic studies in enhancing our understanding of OA pathogenesis over the last three decades. DESIGN: We conducted an extensive literature search focusing on transcriptomics, proteomics and metabolomics within the context of OA. Specifically, we explore how these technologies have identified individual transcripts, proteins, and metabolites, as well as distinctive endotype signatures from various body tissues or fluids of OA patients, including insights at the single-cell level, to advance our understanding of this highly complex disease. RESULTS: Omic studies reveal the description of numerous individual molecules and molecular patterns within OA-associated tissues and fluids. This includes the identification of specific cell (sub)types and associated pathways that contribute to disease mechanisms. However, there remains a necessity to further advance these technologies to delineate the spatial organization of cellular subtypes and molecular patterns within OA-afflicted tissues. CONCLUSIONS: Leveraging a multi-omics approach that integrates datasets from diverse molecular detection technologies, combined with patients' clinical and sociodemographic features, and molecular and regulatory networks, holds promise for identifying unique patient endophenotypes. This holistic approach can illuminate the heterogeneity among OA patients and, in turn, facilitate the development of tailored therapeutic interventions.

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.

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.

ARL3 Mutations Cause Joubert Syndrome by Disrupting Ciliary Protein Composition.

Joubert syndrome (JBTS) is a genetically heterogeneous autosomal-recessive neurodevelopmental ciliopathy. We investigated further the underlying genetic etiology of Joubert syndrome by studying two unrelated families in whom JBTS was not associated with pathogenic variants in known JBTS-associated genes. Combined autozygosity mapping of both families highlighted a candidate locus on chromosome 10 (chr10: 101569997-109106128, UCSC Genome Browser hg 19), and exome sequencing revealed two missense variants in ARL3 within the candidate locus. The encoded protein, ADP ribosylation factor-like GTPase 3 (ARL3), is a small GTP-binding protein that is involved in directing lipid-modified proteins into the cilium in a GTP-dependent manner. Both missense variants replace the highly conserved Arg149 residue, which we show to be necessary for the interaction with its guanine nucleotide exchange factor ARL13B, such that the mutant protein is associated with reduced INPP5E and NPHP3 localization in cilia. We propose that ARL3 provides a potential hub in the network of proteins implicated in ciliopathies, whereby perturbation of ARL3 leads to the mislocalization of multiple ciliary proteins as a result of abnormal displacement of lipidated protein cargo.

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.

Insight into human Miro1/2 domain organization based on the structure of its N-terminal GTPase.

Dysfunction in mitochondrial dynamics is believed to contribute to a host of neurological disorders and has recently been implicated in cancer metastasis. The outer mitochondrial membrane adapter protein Miro functions in the regulation of mitochondrial mobility and degradation, however, the structural basis for its roles in mitochondrial regulation remain unknown. Here, we report a 1.7Å crystal structure of N-terminal GTPase domain (nGTPase) of human Miro1 bound unexpectedly to GTP, thereby revealing a non-catalytic configuration of the putative GTPase active site. We identify two conserved surfaces of the nGTPase, the "SELFYY" and "ITIP" motifs, that are potentially positioned to mediate dimerization or interaction with binding partners. Additionally, we report small angle X-ray scattering (SAXS) data obtained from the intact soluble HsMiro1 and its paralog HsMiro2. Taken together, the data allow modeling of a crescent-shaped assembly of the soluble domain of HsMiro1/2. PDB RSEFERENCE: Crystal structure of the human Miro1 N-terminal GTPase bound to GTP, 6D71.

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.

Red Blood Cell Metabolic Responses to Torpor and Arousal in the Hibernator Arctic Ground Squirrel.

Arctic ground squirrels provide a unique model to investigate metabolic responses to hibernation in mammals. During winter months these rodents are exposed to severe hypothermia, prolonged fasting, and hypoxemia. In the light of their role in oxygen transport/off-loading and owing to the absence of nuclei and organelles (and thus de novo protein synthesis capacity), mature red blood cells have evolved metabolic programs to counteract physiological or pathological hypoxemia. However, red blood cell metabolism in hibernation has not yet been investigated. Here we employed targeted and untargeted metabolomics approaches to investigate erythrocyte metabolism during entrance to torpor to arousal, with a high resolution of the intermediate time points. We report that torpor and arousal promote metabolism through glycolysis and pentose phosphate pathway, respectively, consistent with previous models of oxygen-dependent metabolic modulation in mature erythrocytes. Erythrocytes from hibernating squirrels showed up to 100-fold lower levels of biomarkers of reperfusion injury, such as the pro-inflammatory dicarboxylate succinate. Altered tryptophan metabolism during torpor was here correlated to the accumulation of potentially neurotoxic catabolites kynurenine, quinolinate, and picolinate. Arousal was accompanied by alterations of sulfur metabolism, including sudden spikes in a metabolite putatively identified as thiorphan (level 1 confidence)-a potent inhibitor of several metalloproteases that play a crucial role in nociception and inflammatory complication to reperfusion secondary to ischemia or hemorrhage. Preliminary studies in rats showed that intravenous injection of thiorphan prior to resuscitation mitigates metabolic and cytokine markers of reperfusion injury, etiological contributors to inflammatory complications after shock.

Bialleleic PKD1 mutations underlie early-onset autosomal dominant polycystic kidney disease in Saudi Arabian families.

BACKGROUND: Polycystic kidney disease (PKD) is one of the most common genetic renal diseases and may be inherited in an autosomal dominant or autosomal recessive pattern. Pathogenic variants in two major genes, PKD1 and PKD2, and two rarer genes, GANAB and DNAJB11, cause autosomal dominant PKD (ADPKD). Early onset and severe PKD can occur with PKD1 and PKD2 pathogenic variants and such phenotypes may be modified by second alleles inherited in trans. Homozygous or compound heterozygous hypomorphic PKD1 variants may also cause a moderate to severe disease PKD phenotype. METHODS: Targeted renal gene panel followed by Sanger sequencing of PKD1 gene were employed to investigate molecular causes in early onset PKD patients. RESULTS: In this study, we report four consanguineous Saudi Arabian families with early onset PKD which were associated with biallelic variants in PKD1 gene. CONCLUSIONS: Our findings confirm that PKD1 alleles may combine to produce severe paediatric onset PKD mimicking the more severe autosomal recessive ciliopathy syndromes associated with PKD. Screening of parents of such children may also reveal subclinical PKD phenotypes.

UbMES and UbFluor: Novel probes for ring-between-ring (RBR) E3 ubiquitin ligase PARKIN.

Ring-between-ring (RBR) E3 ligases have been implicated in autoimmune disorders and neurodegenerative diseases. The functions of many RBR E3s are poorly defined, and their regulation is complex, involving post-translational modifications and allosteric regulation with other protein partners. The functional complexity of RBRs, coupled with the complexity of the native ubiquitination reaction that requires ATP and E1 and E2 enzymes, makes it difficult to study these ligases for basic research and therapeutic purposes. To address this challenge, we developed novel chemical probes, ubiquitin C-terminal fluorescein thioesters UbMES and UbFluor, to qualitatively and quantitatively assess the activity of the RBR E3 ligase PARKIN in a simple experimental setup and in real time using fluorescence polarization. First, we confirmed that PARKIN does not require an E2 enzyme for substrate ubiquitination, lysine selection, and polyubiquitin chain formation. Second, we confirmed that UbFluor quantitatively detects naturally occurring activation states of PARKIN caused by Ser65 phosphorylation (pPARKIN) and phosphorylated ubiquitin (pUb). Third, we showed that both pUb and the ubiquitin-accepting substrate contribute to maximal pPARKIN ubiquitin conjugation turnover. pUb enhances the transthiolation step, whereas the substrate clears the pPARKIN∼Ub thioester intermediate. Finally, we established that UbFluor can quantify activation or inhibition of PARKIN by structural mutations. These results demonstrate the feasibility of using UbFluor for quantitative studies of the biochemistry of RBR E3s and for high-throughput screening of small-molecule activators or inhibitors of PARKIN and other RBR E3 ligases.

Neural Signaling Metabolites May Modulate Energy Use in Hibernation.

Despite an epidemic in obesity and metabolic syndrome limited means exist to effect adiposity or metabolic rate other than life style changes. Here we review evidence that neural signaling metabolites may modulate thermoregulatory pathways and offer novel means to fine tune energy use. We extend prior reviews on mechanisms that regulate thermogenesis and energy use in hibernation by focusing primarily on the neural signaling metabolites adenosine, AMP and glutamate.

A service evaluation of an integrated model of palliative care of cystic fibrosis.

BACKGROUND: Patients with advanced cystic fibrosis have severe symptoms with a complex trajectory of exacerbations and recovery. They are often awaiting lung transplantation, and many die without receiving specialist palliative care. AIM: We introduced an integrated model whereby palliative specialists joined the cystic fibrosis team to provide palliative care in parallel with standard care. DESIGN: A service evaluation of this model of care was undertaken in a prospective case series documenting symptoms and outcomes, the views of the cystic fibrosis team and the experience of the palliative specialists. SETTING/PARTICIPANTS: Over 3 years, 28 (10%) of 282 patients attending the cystic fibrosis centre had specialist palliative care. RESULTS: They had advanced lung disease (mean forced expiratory volume in 1 s (FEV1) = 0.86 L (25% predicted)), and 17 died: 6 were on a transplant waiting list at death; 10 were unsuitable and 1 died post transplantation. All who died over these 3 years had specialist palliative care. Four patients had successful transplants. Assessment showed a high prevalence of breathlessness, cough, pain, vomiting and fatigue, with a significant impact on daily life. The cystic fibrosis team rated this model of care highly, felt that palliative care should be members of the team, and thought that patients had found it helpful. The palliative specialists gained knowledge of cystic fibrosis, found it beneficial to meet patients earlier in the disease, and identified unmet needs in managing bereavement and the effects of deaths on other patients with cystic fibrosis. CONCLUSION: This model has been successful in overcoming the difficulties in access to specialist palliative care for patients with cystic fibrosis.

Fourteen monogenic genes account for 15% of nephrolithiasis/nephrocalcinosis.

Nephrolithiasis is a prevalent condition with a high morbidity. Although dozens of monogenic causes have been identified, the fraction of single-gene disease has not been well studied. To determine the percentage of cases that can be molecularly explained by mutations in 1 of 30 known kidney stone genes, we conducted a high-throughput mutation analysis in a cohort of consecutively recruited patients from typical kidney stone clinics. The cohort comprised 272 genetically unresolved individuals (106 children and 166 adults) from 268 families with nephrolithiasis (n=256) or isolated nephrocalcinosis (n=16). We detected 50 likely causative mutations in 14 of 30 analyzed genes, leading to a molecular diagnosis in 14.9% (40 of 268) of all cases; 20 of 50 detected mutations were novel (40%). The cystinuria gene SLC7A9 (n=19) was most frequently mutated. The percentage of monogenic cases was notably high in both the adult (11.4%) and pediatric cohorts (20.8%). Recessive causes were more frequent among children, whereas dominant disease occurred more abundantly in adults. Our study provides an in-depth analysis of monogenic causes of kidney stone disease. We suggest that knowledge of the molecular cause of nephrolithiasis and nephrocalcinosis may have practical implications and might facilitate personalized treatment.

Survey of phosphorylation near drug binding sites in the Protein Data Bank (PDB) and their effects.

While it is currently estimated that 40 to 50% of eukaryotic proteins are phosphorylated, little is known about the frequency and local effects of phosphorylation near pharmaceutical inhibitor binding sites. In this study, we investigated how frequently phosphorylation may affect the binding of drug inhibitors to target proteins. We examined the 453 non-redundant structures of soluble mammalian drug target proteins bound to inhibitors currently available in the Protein Data Bank (PDB). We cross-referenced these structures with phosphorylation data available from the PhosphoSitePlus database. Three hundred twenty-two of 453 (71%) of drug targets have evidence of phosphorylation that has been validated by multiple methods or labs. For 132 of 453 (29%) of those, the phosphorylation site is within 12 Å of the small molecule-binding site, where it would likely alter small molecule binding affinity. We propose a framework for distinguishing between drug-phosphorylation site interactions that are likely to alter the efficacy of drugs versus those that are not. In addition we highlight examples of well-established drug targets, such as estrogen receptor alpha, for which phosphorylation may affect drug affinity and clinical efficacy. Our data suggest that phosphorylation may affect drug binding and efficacy for a significant fraction of drug target proteins.

Structural coupling of the EF hand and C-terminal GTPase domains in the mitochondrial protein Miro.

Miro is a highly conserved calcium-binding GTPase at the regulatory nexus of mitochondrial transport and autophagy. Here we present crystal structures comprising the tandem EF hand and carboxy terminal GTPase (cGTPase) domains of Drosophila Miro. The structures reveal two previously unidentified 'hidden' EF hands, each paired with a canonical EF hand. Each EF hand pair is bound to a helix that structurally mimics an EF hand ligand. A key nucleotide-sensing element and a Pink1 phosphorylation site both lie within an extensive EF hand-cGTPase interface. Our results indicate structural mechanisms for calcium, nucleotide and phosphorylation-dependent regulation of mitochondrial function by Miro.

Mechanism and regulation of kinesin-5, an essential motor for the mitotic spindle.

Mitotic cell division is the most fundamental task of all living cells. Cells have intricate and tightly regulated machinery to ensure that mitosis occurs with appropriate frequency and high fidelity. A core element of this machinery is the kinesin-5 motor protein, which plays essential roles in spindle formation and maintenance. In this review, we discuss how the structural and mechanical properties of kinesin-5 motors uniquely suit them to their mitotic role. We describe some of the small molecule inhibitors and regulatory proteins that act on kinesin-5, and discuss how these regulators may influence the process of cell division. Finally, we touch on some more recently described functions of kinesin-5 motors in non-dividing cells. Throughout, we highlight a number of open questions that impede our understanding of both this motor's function and the potential utility of kinesin-5 inhibitors.

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.