Hereditary axonopathies are frequently caused by mutations in proteins that reside in the endoplasmic reticulum (ER). Which of the many ER functions are pathologically relevant, however, remains to be determined. REEP1 is an ER protein mutated in hereditary spastic paraplegia (HSP) and hereditary motor neuropathy (HMN). We found that HSP-associated missense variants at the N-terminus of REEP1 abolish ER targeting, whereas two more central variants are either rare benign SNPs or confer pathogenicity via a different mechanism. The mis-targeted variants accumulate at lipid droplets (LDs). N-terminal tagging, deletion of the N-terminus, and expression of a minor REEP1 isoform had the same effect. We also confirmed an increase in LD size upon cooverexpression of atlastins and REEP1. Neither wild-type REEP1, LD-targeted HSP variants, nor a non-LD-targeted HMN variant reproduced this effect when expressed alone. We conclude that the N-terminus of REEP1 is necessary for proper targeting to and/or retention in the ER. The protein's potential to also associate with LDs corroborates a synergistic effect with atlastins on LD size. Interestingly, LD size is also altered upon knockdown of seipin, mutations of which also cause HSP and HMN. Regulation of LDs may thus be an ER function critical for long-term axonal maintenance.

Osteoarthritis (OA) is a polygenic disease of older people resulting in the breakdown of cartilage within articular joints. Although a leading cause of disability, there are no disease-modifying therapies. Evidence is emerging to support the origins of OA in skeletogenesis. Whilst methylation QTLs (mQTLs) co-localizing with OA GWAS 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 foetal limb (FL) and cartilage (FC) tissues using a range of molecular genetic techniques. We compared our results to aged cartilage samples (AC) and identified significant OA-mQTLs at 14 CpGs and 29 CpGs in FL and FC tissues, respectively. Differential methylation was observed at 26 sites between foetal and aged cartilage, with the majority becoming actively hypermethylated in old age. Notably, 6/9 OA effector genes showed allelic expression imbalances during foetal 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 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.

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.

Regulation of transcription occurs in a cell type specific manner orchestrated by epigenetic mechanisms including DNA methylation. Methylation changes may also play a key role in lineage specification during stem cell differentiation. To further our understanding of epigenetic regulation in chondrocytes we characterised the DnA methylation changes during chondrogenesis of mesenchymal stem cells (MSCs) by Infinium 450 K methylation array. Significant DNA hypomethylation was identified during chondrogenic differentiation including changes at many key cartilage gene loci. Integration with chondrogenesis gene expression data revealed an enrichment of significant CpGs in upregulated genes, while characterisation of significant CpG loci indicated their predominant localisation to enhancer regions. Comparison with methylation profiles of other tissues, including healthy and diseased adult cartilage, identified chondrocyte-specific regions of hypomethylation and the overlap with differentially methylated cpGs in osteoarthritis. taken together we have associated DnA methylation levels with the chondrocyte phenotype. the consequences of which has potential to improve cartilage generation for tissue engineering purposes and also to provide context for observed methylation changes in cartilage diseases such as osteoarthritis.

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

Dynamic modifications of chromatin allow rapid access of the gene regulatory machinery to condensed genomic regions facilitating subsequent gene expression. Inflammatory cytokine stimulation of cells can cause rapid gene expression changes through direct signalling pathwaymediated transcription factor activation and regulatory element binding. Here we used the Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) to assess regions of the genome that are differentially accessible following treatment of cells with interleukin-1 (IL-1). We identified 126,483 open chromatin regions, with 241 regions significantly differentially accessible following stimulation, with 64 and 177 more or less accessible, respectively. These differentially accessible regions predominantly correspond to regions of the genome marked as enhancers. Motif searching identified an overrepresentation of a number of transcription factors, most notably RelA, in the regions becoming more accessible, with analysis of ChIP-seq data confirmed RelA binding to these regions. A significant correlation in differential chromatin accessibility and gene expression was also observed. Functionality in regulating gene expression was confirmed using CRISPR/Cas9 genome-editing to delete regions that became more accessible following stimulation in the genes MMP13, IKBKE and C1QTNF1. These same regions were also accessible for activation using a dCas9-transcriptional activator and showed enhancer activity in a cellular model. Together, these data describe and functionally validate a number of dynamically accessible chromatin regions involved in inflammatory signalling.

Osteoarthritis (OA) is a polygenic disease of older people resulting in the breakdown of cartilage within articular joints. Although a leading cause of disability, there are no disease-modifying therapies. Evidence is emerging to support the origins of OA in skeletogenesis. Whilst methylation QTLs (mQTLs) co-localizing with OA GWAS 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 compared our results to aged cartilage samples (AC) and identified significant OA-mQTLs at 14 CpGs and 29 CpGs in FL and FC tissues, respectively. Differential methylation was observed at 26 sites between fetal and aged cartilage, 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 skeletogenesis. We show striking differences in the spatiotemporal function of these loci, contributing to our understanding of OA etiology, with implications for the timing and strategy of pharmacological interventions.

Human mesenchymal stem cells are able to differentiate into chondrocytes, the cell type found in cartilage, making them an accessible system to study gene regulation during this process. Epigenetic mechanisms such as histone modifications and DNA methylation together with transcription factor binding play a role in activating and repressing gene expression. In this study, we investigated the genome-wide histone modification changes during chondrocyte differentiation. Integration of this data with DNA methylation and SOX9 transcription factor ChIP-seq revealed epigenetic changes at gene enhancer elements.Regions of the genome that transition from non-enhancers to enhancers in chondrocytes are enriched for SOX9 transcription factor binding sites. Luciferase reporter assays revealed that enhancer activity may be modulated by manipulating DNA methylation and SOX9 expression.This study has defined important regulatory elements in chondrocytes which could serve as targets for future mechanistic studies. AbstractEpigenetic mechanisms are known to regulate gene expression during chondrogenesis. In this study, we have characterised the epigenome during in vitro differentiation of human mesenchymal stem cells (hMSCs) into chondrocytes. Chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) was used to assess a range of N-terminal post-transcriptional modifications (marks) to histone H3 lysines (H3K4me3, H3K4me1, H3K27ac, H3K27me3 and H3K36me3) in both hMSCs and differentiated chondrocytes.Chromatin states were characterised using histone ChIP-seq and cis-regulatory elements were identified in chondrocytes. Chondrocyte enhancers were associated with chondrogenesis related gene ontology (GO) terms. In silico analysis and integration of DNA methylation data with chondrogenesis chromatin states revealed that enhancers marked by histone marks H3K4me1 and H3K27ac were de-methylated during in vitro chondrogenesis. Similarity analysis between hMSC and chondrocyte chromatin states defined in this studywith epigenomes of cell-types defined by the Roadmap Epigenomics project revealed that enhancers are more distinct between cell-types compared to other chromatin states. Motif analysis revealed that the transcription factor SOX9 is enriched in chondrocyte enhancers.Luciferase reporter assays confirmed that chondrocyte enhancers characterised in this study exhibited enhancer activity which may be modulated by inducing DNA methylation and SOX9 overexpression. Altogether, these integrated data illustrate the cross-talk between different epigenetic mechanisms during chondrocyte differentiation.