Dual proteolytic pathways govern glycolysis and immune competence.

Proteasomes and lysosomes constitute the major cellular systems that catabolize proteins to recycle free amino acids for energy and new protein synthesis. Tripeptidyl peptidase II (TPPII) is a large cytosolic proteolytic complex that functions in tandem with the proteasome-ubiquitin protein degradation pathway. We found that autosomal recessive TPP2 mutations cause recurrent infections, autoimmunity, and neurodevelopmental delay in humans. We show that a major function of TPPII in mammalian cells is to maintain amino acid levels and that TPPII-deficient cells compensate by increasing lysosome number and proteolytic activity. However, the overabundant lysosomes derange cellular metabolism by consuming the key glycolytic enzyme hexokinase-2 through chaperone-mediated autophagy. This reduces glycolysis and impairs the production of effector cytokines, including IFN-γ and IL-1ß. Thus, TPPII controls the balance between intracellular amino acid availability, lysosome number, and glycolysis, which is vital for adaptive and innate immunity and neurodevelopmental health.

A Bayesian network approach incorporating imputation of missing data enables exploratory analysis of complex causal biological relationships.

Bayesian networks can be used to identify possible causal relationships between variables based on their conditional dependencies and independencies, which can be particularly useful in complex biological scenarios with many measured variables. Here we propose two improvements to an existing method for Bayesian network analysis, designed to increase the power to detect potential causal relationships between variables (including potentially a mixture of both discrete and continuous variables). Our first improvement relates to the treatment of missing data. When there is missing data, the standard approach is to remove every individual with any missing data before performing analysis. This can be wasteful and undesirable when there are many individuals with missing data, perhaps with only one or a few variables missing. This motivates the use of imputation. We present a new imputation method that uses a version of nearest neighbour imputation, whereby missing data from one individual is replaced with data from another individual, their nearest neighbour. For each individual with missing data, the subsets of variables to be used to select the nearest neighbour are chosen by sampling without replacement the complete data and estimating a best fit Bayesian network. We show that this approach leads to marked improvements in the recall and precision of directed edges in the final network identified, and we illustrate the approach through application to data from a recent study investigating the causal relationship between methylation and gene expression in early inflammatory arthritis patients. We also describe a second improvement in the form of a pseudo-Bayesian approach for upweighting certain network edges, which can be useful when there is prior evidence concerning their directions.

microRNA-seq of cartilage reveals an overabundance of miR-140-3p which contains functional isomiRs.

miR-140 is selectively expressed in cartilage. Deletion of the entire Mir140 locus in mice results in growth retardation and early-onset osteoarthritis-like pathology; however, the relative contribution of miR-140-5p or miR-140-3p to the phenotype remains to be determined. An unbiased small RNA sequencing approach identified miR-140-3p as significantly more abundant (>10-fold) than miR-140-5p in human cartilage. Analysis of these data identified multiple miR-140-3p isomiRs differing from the miRBase annotation at both the 5' and 3' end, with >99% having one of two seed sequences (5' bases 2-8). Canonical (miR-140-3p.2) and shifted (miR-140-3p.1) seed isomiRs were overexpressed in chondrocytes and transcriptomics performed to identify targets. miR-140-3p.1 and miR-140-3p.2 significantly down-regulated 694 and 238 genes, respectively, of which only 162 genes were commonly down-regulated. IsomiR targets were validated using 3'UTR luciferase assays. miR-140-3p.1 targets were enriched within up-regulated genes in rib chondrocytes of Mir140-null mice and within down-regulated genes during human chondrogenesis. Finally, through imputing the expression of miR-140 from the expression of the host gene WWP2 in 124 previously published data sets, an inverse correlation with miR-140-3p.1 predicted targets was identified. Together these data suggest the novel seed containing isomiR miR-140-3p.1 is more functional than original consensus miR-140-3p seed containing isomiR.

Interferon-α-mediated therapeutic resistance in early rheumatoid arthritis implicates epigenetic reprogramming.

OBJECTIVES: An interferon (IFN) gene signature (IGS) is present in approximately 50% of early, treatment naive rheumatoid arthritis (eRA) patients where it has been shown to negatively impact initial response to treatment. We wished to validate this effect and explore potential mechanisms of action. METHODS: In a multicentre inception cohort of eRA patients (n=191), we examined the whole blood IGS (MxA, IFI44L, OAS1, IFI6, ISG15) with reference to circulating IFN proteins, clinical outcomes and epigenetic influences on circulating CD19+ B and CD4+ T lymphocytes. RESULTS: We reproduced our previous findings demonstrating a raised baseline IGS. We additionally showed, for the first time, that the IGS in eRA reflects circulating IFN-α protein. Paired longitudinal analysis demonstrated a significant reduction between baseline and 6-month IGS and IFN-α levels (p<0.0001 for both). Despite this fall, a raised baseline IGS predicted worse 6-month clinical outcomes such as increased disease activity score (DAS-28, p=0.025) and lower likelihood of a good EULAR clinical response (p=0.034), which was independent of other conventional predictors of disease activity and clinical response. Molecular analysis of CD4+ T cells and CD19+ B cells demonstrated differentially methylated CPG sites and dysregulated expression of disease relevant genes, including PARP9, STAT1, and EPSTI1, associated with baseline IGS/IFNα levels. Differentially methylated CPG sites implicated altered transcription factor binding in B cells (GATA3, ETSI, NFATC2, EZH2) and T cells (p300, HIF1α). CONCLUSIONS: Our data suggest that, in eRA, IFN-α can cause a sustained, epigenetically mediated, pathogenic increase in lymphocyte activation and proliferation, and that the IGS is, therefore, a robust prognostic biomarker. Its persistent harmful effects provide a rationale for the initial therapeutic targeting of IFN-α in selected patients with eRA.

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

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

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.

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.

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.

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

Osteoarthritis (OA) is a common, painful and debilitating disease of articulating joints resulting from the age-associated loss of cartilage. Well-powered genetic studies have identified a number of DNA polymorphisms that are associated with OA susceptibility. Like most complex trait loci, these OA loci are thought to influence disease susceptibility through the regulation of gene expression, so-called expression quantitative loci, or eQTLs. One mechanism through which eQTLs act is epigenetic, by modulating DNA methylation. In such cases, there are quantitative differences in DNA methylation between the two alleles of the causal polymorphism, with the association signal referred to as a methylation quantitative trait locus, or meQTL. In this study, we aimed to investigate whether the OA susceptibility loci identified to date are functioning as meQTLs by integrating genotype data with whole genome methylation data of cartilage DNA. We investigated potential genotype-methylation correlations within a 1.0-1.5 Mb region surrounding each of 16 OA-associated single-nucleotide polymorphisms (SNPs) in 99 cartilage samples and identified four that function as meQTLs. Three of these replicated in an additional cohort of up to 62 OA patients. These observations suggest that OA susceptibility loci regulate the level of DNA methylation in cis and provide a mechanistic explanation as to how these loci impact upon OA susceptibility, further increasing our understanding of the role of genetics and epigenetics in this common disease.

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.

Defective Leukocyte Adhesion and Chemotaxis Contributes to Combined Immunodeficiency in Humans with Autosomal Recessive MST1 Deficiency.

PURPOSE: To investigate the clinical and functional aspects of MST1 (STK4) deficiency in a profoundly CD4-lymphopenic kindred with a novel homozygous nonsense mutation in STK4. Although recent studies have described the cellular effects of murine Mst1 deficiency, the phenotype of MST1-deficient human lymphocytes has yet to be fully explored. Patient lymphocytes were therefore investigated in the context of current knowledge of murine Mst1 deficiency. METHODS: Genetic etiology was identified by whole exome sequencing of genomic DNA from two siblings, combined with linkage analysis in the wider family. MST1 protein expression was assessed by immunoblotting. The ability of patient lymphocytes to adhere to ICAM-1 under flow conditions was measured, and transwell assays were used to assess chemotaxis. Chemokine receptor expression was examined by flow cytometry and receptor signalling by immunoblotting. RESULTS: A homozygous nonsense mutation in STK4 (c.442C > T, p.Arg148Stop) was found in the patients, leading to a lack of MST1 protein expression. Patient leukocytes exhibited deficient chemotaxis after stimulation with CXCL11, despite preserved expression of CXCR3. Patient lymphocytes were also unable to bind effectively to immobilised ICAM-1 under flow conditions, in keeping with a failure to develop high affinity binding. CONCLUSION: The observed abnormalities of adhesion and migration imply a profound trafficking defect among human MST1-deficient lymphocytes. By analogy with murine Mst1 deficiency and other defects of leucocyte trafficking, this is likely to contribute to immunodeficiency by impairing key aspects of T-cell development and function such as positive selection in the thymus, thymic egress and immune synapse formation in the periphery.

The identification of trans-acting factors that regulate the expression of GDF5 via the osteoarthritis susceptibility SNP rs143383.

rs143383 is a C to T transition SNP located in the 5'untranslated region (5'UTR) of the growth differentiation factor 5 gene GDF5. The T allele of the SNP is associated with increased risk of osteoarthritis (OA) in Europeans and in Asians. This susceptibility is mediated by the T allele producing less GDF5 transcript relative to the C allele, a phenomenon known as differential allelic expression (DAE). The aim of this study was to identify trans-acting factors that bind to rs143383 and which regulate this GDF5 DAE. Protein binding to the gene was investigated by two experimental approaches: 1) competition and supershift electrophoretic mobility shift assays (EMSAs) and 2) an oligonucleotide pull down assay followed by quantitative mass spectrometry. Binding was then confirmed in vivo by chromatin immunoprecipitation (ChIP), and the functional effects of candidate proteins investigated by RNA interference (RNAi) and over expression. Using these approaches the trans-acting factors Sp1, Sp3, P15, and DEAF-1 were identified as interacting with the GDF5 5'UTR. Knockdown and over expression of the factors demonstrated that Sp1, Sp3, and DEAF-1 are repressors of GDF5 expression. Depletion of DEAF-1 modulated the DAE of GDF5 and this differential allelic effect was confirmed following over expression, with the rs143383 T allele being repressed to a significantly greater extent than the rs143383 C allele. In combination, Sp1 and DEAF-1 had the greatest repressive activity. In conclusion, we have identified four trans-acting factors that are binding to GDF5, three of which are modulating GDF5 expression via the OA susceptibility locus rs143383.

Differential DNA methylation and expression of inflammatory and zinc transporter genes defines subgroups of osteoarthritic hip patients.

OBJECTIVES: We have previously shown that the cartilage DNA methylome delineates two clusters of osteoarthritic (OA) hip patients, characterised by differential methylation of inflammatory genes, while others have demonstrated a link between zinc homeostasis and inflammation in OA. We aimed to investigate these effects at the methylation and gene expression level. METHODS: We used our previously generated methylation data while quantitative PCR was used to measure gene expression using RNA from the hip cartilage of members of both clusters and from control individuals without hip OA. RESULTS: One of the OA clusters is characterised by the promoter hypomethylation and increased expression of inflammation-associated genes including IL1A and TNF. Furthermore, we show that the increase in expression of these genes is accompanied by increased expression of several zinc transporter genes. In addition, the zinc responsive transcription factor MTF1 is also upregulated, which is accompanied by an increase in the expression of its targets the metalloproteinases MMP13 and ADAMTS5. CONCLUSIONS: We have identified a subgroup of OA hip patients that are epigenetically and transcriptiomically characterised by a cartilage inflammatory phenotype with concurrent differential regulation of zinc regulators. The identification of subgroups enhances stratified phenotyping of OA patients and has important implications for future therapeutic applications.

Functional characterisation of the osteoarthritis susceptibility locus at chromosome 6q14.1 marked by the polymorphism rs9350591.

BACKGROUND: The arcOGEN genome-wide association study reported the rs9350591 C/T single nucleotide polymorphism (SNP) as marking a region on chromosome 6q14.1 that is associated with hip osteoarthritis (OA) in Europeans, with an odds ratio (OR) of 1.18 and a p-value of 2.42 × 10(-9). rs9350591 is an intergenic SNP surrounded by seven genes within 1 Mb. Six of the genes are expressed in cartilage. We sought to characterise this signal to assess whether the association of rs9350591 with OA is mediated by modulating gene expression. METHODS: Total RNA was extracted from hip or knee cartilage of 161 OA patients and from hip cartilage of 29 non-OA patients who had undergone hip replacements as a result of neck-of-femur (NOF) fractures. We used quantitative PCR (qPCR) to measure overall gene expression, and pyrosequencing to assess allelic expression of the genes. A mesenchymal stem cell (MSC) differentiation model was used to assess gene expression during chondrogenesis. RESULTS: We identified a significant decrease in the expression of SENP6 (p = 0.005) and MYO6 (p = 0.026) in OA hip cartilage relative to the non-OA hip control cartilage. However, we found no evidence for a correlation between gene expression and rs9350591 genotype for any of the six genes. In addition, we identified expression quantitative trait loci (eQTLs) operating on COL12A1, TMEM30A, SENP6 and MYO6, although these were not relevant to the OA associated signal. Finally, all genes were dynamically expressed during chondrogenesis. CONCLUSIONS: The regulation of gene expression at this locus is complex, highlighted by the down-regulation of SENP6 and MYO6 in OA hip cartilage and by eQTLs operating on four of the genes at the locus. However, modulation of gene expression in the end-stage OA cartilage that we have investigated is not the mechanism by which this association signal is operating. As implied by the dynamic patterns of gene expression throughout chondrogenesis, the association signal marked by rs9350591 could instead be exerting its effects during joint development.

Expression analysis of the osteoarthritis genetic susceptibility locus mapping to an intron of the MCF2L gene and marked by the polymorphism rs11842874.

BACKGROUND: Osteoarthritis (OA) is a painful, debilitating disease characterised by loss of articular cartilage with concurrent changes in other tissues of the synovial joint. Genetic association studies have shown that a number of common variants increase the risk of developing OA. Investigating their activity can uncover novel causal pathways and potentially highlight new treatment targets. One of the reported OA association signals is marked by the single nucleotide polymorphism (SNP) rs11842874 at chromosome 13q34. rs11842874 is positioned within a small linkage disequilibrium (LD) block within intron 4 of MCF2L, a gene encoding guanine-nucleotide exchange factor DBS. There are no non-synonymous SNPs that correlate with this association signal and we therefore set out to assess whether its effect on OA susceptibility is mediated by alteration of MCF2L expression. METHODS: Nucleic acid was extracted from cartilage, synovial membrane or infrapatellar fat pad tissues from OA patients. Expression of MCF2L was measured by quantitative PCR and RNA-sequencing whilst the presence of DBS was studied using immunohistochemistry. The functional effect of SNPs within the 13q34 locus was assessed using public databases and in vitro using luciferase reporter analysis. RESULTS: MCF2L gene and protein expression are detectable in joint tissues, with quantitative differences in the expression of the gene and in the transcript isoforms expressed between the tissues tested. There is an expression quantitative trait locus (eQTL) operating within synovial membrane tissue, with possession of the risk-conferring A allele of rs11842874 correlating with increased MCF2L expression. SNPs within the rs11842874 LD block reside within transcriptional regulatory elements and their direct analysis reveals that several show quantitative differences in regulatory activity at the allelic level. CONCLUSIONS: MCF2L is subject to a cis-acting eQTL in synovial membrane that correlates with the OA association signal. This signal contains several functional SNPs that could account for the susceptibility and which therefore merit further investigation. As far as we are aware, this is the first example of an OA susceptibility locus operating as an eQTL in synovial membrane tissue but not in cartilage.

CpG methylation regulates allelic expression of GDF5 by modulating binding of SP1 and SP3 repressor proteins to the osteoarthritis susceptibility SNP rs143383.

GDF5 encodes an extracellular signalling molecule that is essential for normal skeletal development. The rs144383 C to T SNP located in the 5'UTR of this gene is functional and has a pleiotropic effect on the musculoskeletal system, being a risk factor for knee-osteoarthritis (OA), congenital hip dysplasia, lumbar disc degeneration and Achilles tendon pathology. rs143383 exerts a joint-wide effect on GDF5 expression, with expression of the OA-associated T allele being significantly reduced relative to the C allele, termed allelic expression imbalance. We have previously reported that the GDF5 locus is subject to DNA methylation and that allelic imbalance of rs143383 is mediated by SP1, SP3 and DEAF1 transcriptional repressors. In this study, we have assayed GDF5 methylation in normal and osteoarthritic cartilage, and investigated the effect of methylation on the allelic imbalance of rs143383. We observed demethylation of the GDF5 5'UTR in OA knee cartilage relative to both OA (p = 0.009) and non-OA (p = 0.001) hip cartilage, with the most significant demethylation observed at the highly conserved +37 CpG site located 4 bp upstream of rs143383. Methylation modulates the level and direction of allelic imbalance of rs143383, with methylation of the +37 CpG dinucleotide within the SP1/SP3 binding site having an allele-specific effect on SP1 and SP3 binding. Furthermore, methylation attenuated the repressive effects of SP1, SP3 and DEAF1 on GDF5 promoter activity. This data suggest that the differential methylation of the +37 CpG site between osteoarthritic hip and knee cartilage may be responsible for the knee-specific effect of rs143383 on OA susceptibility.

Allelic expression analysis of the osteoarthritis susceptibility locus that maps to chromosome 3p21 reveals cis-acting eQTLs at GNL3 and SPCS1.

BACKGROUND: An osteoarthritis (OA) susceptibility locus has been mapped to chromosome 3p21, to a region of high linkage disequilibrium encompassing twelve genes. Six of these genes are expressed in joint tissues and we therefore assessed whether any of the six were subject to cis-acting regulatory polymorphisms active in these tissues and which could therefore account for the association signal. METHODS: We measured allelic expression using pyrosequencing assays that can distinguish mRNA output from each allele of a transcript single nucleotide polymorphism. We assessed RNA extracted from the cartilage and other joint tissues of OA patients who had undergone elective joint replacement surgery. A two-tailed Mann-Whitney exact test was used to test the significance of any allelic differences. RESULTS: GNL3 and SPCS1 demonstrated significant allelic expression imbalance (AEI) in OA cartilage (GNL3, mean AEI = 1.04, p = 0.0002; SPCS1, mean AEI = 1.07, p < 0.0001). Similar results were observed in other tissues. Expression of the OA-associated allele was lower than that of the non-associated allele for both genes. CONCLUSIONS: cis-acting regulatory polymorphisms acting on GNL3 and SPCS1 contribute to the OA association signal at chromosome 3p21, and these genes therefore merit further investigation.

Expression of the osteoarthritis-associated gene GDF5 is modulated epigenetically by DNA methylation.

GDF5 is involved in synovial joint development, maintenance and repair, and the rs143383 C/T single nucleotide polymorphism (SNP) located in the 5'UTR of GDF5 is associated, at the genome-wide significance level, with osteoarthritis susceptibility, and with other musculoskeletal phenotypes including height, congenital hip dysplasia and Achilles tendinopathy. There is a significant reduction in the expression of the disease-associated T allele relative to the C allele in synovial joint tissues, an effect influenced by a second SNP (rs143384, C/T) also within the 5'UTR. The differential allelic expression (DAE) imbalance of the C and T alleles of rs143383 varies intra- and inter-individually, suggesting that DAE may be modulated epigenetically. The C alleles of both SNPs form CpG dinucleotides that are potentially amenable to regulation by methylation. Here, we have examined whether DNA methylation regulates GDF5 expression and the allelic imbalance caused by rs143383. We observed methylation of the GDF5 promoter and 5'UTR in cell lines and joint tissues, with demethylation correlating with increased GDF5 expression. The CpG sites created by the C alleles at rs143383 and rs143384 were variably methylated, and treatment of a heterozygous cell line with a demethylating agent further increased the allelic expression imbalance between the C and T alleles. This demonstrates that the genetic effect of the rs143383 SNP on GDF5 expression is modulated epigenetically by DNA methylation. The variability in DAE of rs143383 is therefore partly accounted for by differences in DNA methylation that could influence the penetrance of this allele in susceptibility to common musculoskeletal diseases.

The genetics and functional analysis of primary osteoarthritis susceptibility.

Recent genome-wide association scans (GWASs) along with several adequately powered candidate gene studies have yielded a number of risk alleles for osteoarthritis (OA). This number is now sufficiently large to allow conclusions to be drawn regarding the nature of genetic susceptibility, including the fact that the risk alleles have variable effects depending on sex, ethnicity and on the skeletal site of the disease. Several of the alleles that have emerged from the GWASs are within or close to highly plausible candidate genes, including RUNX2 and CHST11. However, the majority of risk alleles do not map to genes previously reported to play a role in musculoskeletal biology, indicating that the GWAS datasets are telling us something new about the OA disease process. Functional studies have so far revealed that effects on gene expression are likely to be one of the main mechanisms through which OA susceptibility is acting. Epigenetic mechanisms such as DNA methylation also influence OA risk, and integration of genetic, transcriptomic and epigenetic data will allow us to use the genetic discoveries for informed development of new OA biological treatments.

Exome sequencing identifies GATA-2 mutation as the cause of dendritic cell, monocyte, B and NK lymphoid deficiency.

The human syndrome of dendritic cell, monocyte, B and natural killer lymphoid deficiency presents as a sporadic or autosomal dominant trait causing susceptibility to mycobacterial and other infections, predisposition to myelodysplasia and leukemia, and, in some cases, pulmonary alveolar proteinosis. Seeking a genetic cause, we sequenced the exomes of 4 unrelated persons, 3 with sporadic disease, looking for novel, heterozygous, and probably deleterious variants. A number of genes harbored novel variants in person, but only one gene, GATA2, was mutated in all 4 persons. Each person harbored a different mutation, but all were predicted to be highly deleterious and to cause loss or mutation of the C-terminal zinc finger domain. Because GATA2 is the only common mutated gene in 4 unrelated persons, it is highly probable to be the cause of dendritic cell, monocyte, B, and natural killer lymphoid deficiency. This disorder therefore constitutes a new genetic form of heritable immunodeficiency and leukemic transformation.