The effect of FTO variation on increased osteoarthritis risk is mediated through body mass index : a mendelian randomisation study

Objective: Variation in the fat mass and obesity-associated (FTO) gene influences susceptibility to obesity. A variant in the FTO gene has been implicated in genetic risk to osteoarthritis (OA). We examined the role of the FTO polymorphism rs8044769 in risk of knee and hip OA in cases and controls incorporating body mass index (BMI) information. Methods: 5409 knee OA patients, 4355 hip OA patients and up to 5362 healthy controls from 7 independent cohorts from the UK and Australia were genotyped for rs8044769. The association of the FTO variant with OA was investigated in case/control analyses with and without BMI adjustment and in analyses matched for BMI category. A mendelian randomisation approach was employed using the FTO variant as the instrumental variable to evaluate the role of overweight on OA. Results: In the meta-analysis of all overweight (BMI≥25) samples versus normal-weight controls irrespective of OA status the association of rs8044769 with overweight is highly significant (OR[CIs] for allele G=1.14 [01.08 to 1.19], p=7.5×10−7). A significant association with knee OA is present in the analysis without BMI adjustment (OR[CIs]=1.08[1.02 to 1.14], p=0.009) but the signal fully attenuates after BMI adjustment (OR[CIs]=0.99[0.93 to 1.05], p=0.666). We observe no evidence for association in the BMI-matched meta-analyses. Using mendelian randomisation approaches we confirm the causal role of overweight on OA. Conclusions: Our data highlight the contribution of genetic risk to overweight in defining risk to OA but the association is exclusively mediated by the effect on BMI. This is consistent with what is known of the biology of the FTO gene and supports the causative role of high BMI in OA

Allelic expression analysis of the osteoarthritis susceptibility locus that maps to MICAL3

<p>Abstract</p> <p>Background</p> <p>A genome-wide association scan with subsequent replication study that involved over 67,000 individuals of European ancestry has produced evidence of association of single nucleotide polymorphism rs2277831 to primary osteoarthritis (OA) with a P-value of 2.9 × 10^-5. rs2277831, an A/G transition, is located in an intron of <it>MICAL3</it>. This gene is located on chromosome 22q11.21 and the association signal encompasses two additional genes, <it>BCL2L13 </it>and <it>BID</it>. It is becoming increasingly apparent that many common complex traits are mediated by <it>cis</it>-acting regulatory polymorphisms that influence, in a tissue-specific manner, gene expression or transcript stability.</p> <p>Methods</p> <p>We used total and allelic expression analysis to assess whether the OA association to rs2277831 is mediated by an influence on MICAL3, BCL2L13 or BID expression. Using RNA extracted from joint tissues of 60 patients who had undergone elective joint replacement surgery, we assessed whether rs2277831 correlated with allelic expression of either of the three genes by: 1) measuring the expression of each gene by quantitative PCR and then stratifying the data by genotype at rs2277831 and 2) accurately discriminating and quantifying the mRNA synthesised from the alleles of OA patients using allelic-quantitative PCR.</p> <p>Results</p> <p>We found no evidence for a correlation between gene expression and genotype at rs2277831, with P-values of 0.09 for <it>BCL2L13</it>, 0.07 for <it>BID </it>and 0.33 for <it>MICAL3</it>. In the allelic expression analysis we observed several examples of significant (p < 0.05) allelic imbalances, with an allelic expression ratio of 2.82 observed in <it>BCL2L13 </it>(P = 0.004), 2.09 at <it>BID </it>(P = 0.001) and the most extreme case being at <it>MICAL3</it>, with an allelic expression ratio of 5.47 (P = 0.001). However, there was no correlation observed between the pattern of allelic expression and the genotype at rs2277831.</p> <p>Conclusions</p> <p>In the tissues that we have studied, our data do not support our hypothesis that the association between rs2277831 and OA is due to the effect this SNP has on <it>MICAL3, BCL2L13 </it>or <it>BID </it>gene expression. Instead, our data point towards other functional effects accounting for the OA associated signal.</p

Identification of new susceptibility loci for osteoarthritis (arcOGEN):a genome-wide association study

To access publisher's full text version of this article. Please click on the hyperlink in Additional Links field.Osteoarthritis is the most common form of arthritis worldwide and is a major cause of pain and disability in elderly people. The health economic burden of osteoarthritis is increasing commensurate with obesity prevalence and longevity. Osteoarthritis has a strong genetic component but the success of previous genetic studies has been restricted due to insufficient sample sizes and phenotype heterogeneity. We undertook a large genome-wide association study (GWAS) in 7410 unrelated and retrospectively and prospectively selected patients with severe osteoarthritis in the arcOGEN study, 80% of whom had undergone total joint replacement, and 11,009 unrelated controls from the UK. We replicated the most promising signals in an independent set of up to 7473 cases and 42,938 controls, from studies in Iceland, Estonia, the Netherlands, and the UK. All patients and controls were of European descent. We identified five genome-wide significant loci (binomial test p≤5·0×10(-8)) for association with osteoarthritis and three loci just below this threshold. The strongest association was on chromosome 3 with rs6976 (odds ratio 1·12 [95% CI 1·08-1·16]; p=7·24×10(-11)), which is in perfect linkage disequilibrium with rs11177. This SNP encodes a missense polymorphism within the nucleostemin-encoding gene GNL3. Levels of nucleostemin were raised in chondrocytes from patients with osteoarthritis in functional studies. Other significant loci were on chromosome 9 close to ASTN2, chromosome 6 between FILIP1 and SENP6, chromosome 12 close to KLHDC5 and PTHLH, and in another region of chromosome 12 close to CHST11. One of the signals close to genome-wide significance was within the FTO gene, which is involved in regulation of bodyweight-a strong risk factor for osteoarthritis. All risk variants were common in frequency and exerted small effects. Our findings provide insight into the genetics of arthritis and identify new pathways that might be amenable to future therapeutic intervention.Arthritis Research UK 1803

Rare variation at the TNFAIP3 locus and susceptibility to rheumatoid arthritis

Genome-wide association studies (GWAS) conducted using commercial single nucleotide polymorphisms (SNP) arrays have proven to be a powerful tool for the detection of common disease susceptibility variants. However, their utility for the detection of lower frequency variants is yet to be practically investigated. Here we describe the application of a rare variant collapsing method to a large genome-wide SNP dataset, the Wellcome Trust Case Control Consortium rheumatoid arthritis (RA) GWAS. We partitioned the data into gene-centric bins and collapsed genotypes of low frequency variants (defined here as MAF ≤0.05) into a single count coupled with univariate analysis. We then prioritised gene regions for further investigation in an independent cohort of 3,355 cases and 2,427 controls based on rare variant signal p value and prior evidence to support involvement in RA. A total of 14,536 gene bins were investigated in the primary analysis and signals mapping to the TNFAIP3 and chr17q24 loci were selected for further investigation. We detected replicating association to low frequency variants in the TNFAIP3 gene (combined p = 6.6 × 10−6). Even though rare variants are not well-represented and can be difficult to genotype in GWAS, our study supports the application of low frequency variant collapsing methods to genome-wide SNP datasets as a means of exploiting data that are routinely ignored

Whole-Genome Sequencing Coupled to Imputation Discovers Genetic Signals for Anthropometric Traits

Deep sequence-based imputation can enhance the discovery power of genome-wide association studies by assessing previously unexplored variation across the common-and low-frequency spectra. We applied a hybrid whole-genome sequencing (WGS) and deep imputation approach to examine the broader allelic architecture of 12 anthropometric traits associated with height, body mass, and fat distribution in up to 267,616 individuals. We report 106 genome-wide significant signals that have not been previously identified, including 9 low-frequency variants pointing to functional candidates. Of the 106 signals, 6 are in genomic regions that have not been implicated with related traits before, 28 are independent signals at previously reported regions, and 72 represent previously reported signals for a different anthropometric trait. 71% of signals reside within genes and fine mapping resolves 23 signals to one or two likely causal variants. We confirm genetic overlap between human monogenic and polygenic anthropometric traits and find signal enrichment in cis expression QTLs in relevant tissues. Our results highlight the potential of WGS strategies to enhance biologically relevant discoveries across the frequency spectrum.Peer reviewe

Age-dependent alteration of TGF-β signalling in osteoarthritis

Osteoarthritis (OA) is a disease of articular cartilage, with aging as the main risk factor. In OA, changes in chondrocytes lead to the autolytic destruction of cartilage. Transforming growth factor-β has recently been demonstrated to signal not only via activin receptor-like kinase 5 (ALK5)-induced Smad2/3 phosphorylation, but also via ALK1-induced Smad1/5/8 phosphorylation in articular cartilage. In aging cartilage and experimental OA, the ratio ALK1/ALK5 has been found to be increased, and the expression of ALK1 is correlated with matrix metalloproteinase-13 expression. The age-dependent shift towards Smad1/5/8 signalling might trigger the differentiation of articular chondrocytes with an autolytic phenotype

Evaluation of the genetic overlap between osteoarthritis with body mass index and height using genome-wide association scan data.

OBJECTIVES: Obesity as measured by body mass index (BMI) is one of the major risk factors for osteoarthritis. In addition, genetic overlap has been reported between osteoarthritis and normal adult height variation. We investigated whether this relationship is due to a shared genetic aetiology on a genome-wide scale. METHODS: We compared genetic association summary statistics (effect size, p value) for BMI and height from the GIANT consortium genome-wide association study (GWAS) with genetic association summary statistics from the arcOGEN consortium osteoarthritis GWAS. Significance was evaluated by permutation. Replication of osteoarthritis association of the highlighted signals was investigated in an independent dataset. Phenotypic information of height and BMI was accounted for in a separate analysis using osteoarthritis-free controls. RESULTS: We found significant overlap between osteoarthritis and height (p=3.3×10(-5) for signals with p≤0.05) when the GIANT and arcOGEN GWAS were compared. For signals with p≤0.001 we found 17 shared signals between osteoarthritis and height and four between osteoarthritis and BMI. However, only one of the height or BMI signals that had shown evidence of association with osteoarthritis in the arcOGEN GWAS was also associated with osteoarthritis in the independent dataset: rs12149832, within the FTO gene (combined p=2.3×10(-5)). As expected, this signal was attenuated when we adjusted for BMI. CONCLUSIONS: We found a significant excess of shared signals between both osteoarthritis and height and osteoarthritis and BMI, suggestive of a common genetic aetiology. However, only one signal showed association with osteoarthritis when followed up in a new dataset

Cardiometabolic effects of genetic upregulation of the interleukin 1 receptor antagonist: a Mendelian randomisation analysis

Background To investigate potential cardiovascular and other effects of long-term pharmacological interleukin 1 (IL-1) inhibition, we studied genetic variants that produce inhibition of IL-1, a master regulator of inflammation. Methods We created a genetic score combining the effects of alleles of two common variants (rs6743376 and rs1542176) that are located upstream of IL1RN, the gene encoding the IL-1 receptor antagonist (IL-1Ra; an endogenous inhibitor of both IL-1 alpha and IL-1 beta); both alleles increase soluble IL-1Ra protein concentration. We compared effects on inflammation biomarkers of this genetic score with those of anakinra, the recombinant form of IL-1Ra, which has previously been studied in randomised trials of rheumatoid arthritis and other inflammatory disorders. In primary analyses, we investigated the score in relation to rheumatoid arthritis and four cardiometabolic diseases (type 2 diabetes, coronary heart disease, ischaemic stroke, and abdominal aortic aneurysm; 453 411 total participants). In exploratory analyses, we studied the relation of the score to many disease traits and to 24 other disorders of proposed relevance to IL-1 signalling (746 171 total participants). Findings For each IL1RN minor allele inherited, serum concentrations of IL-1Ra increased by 0.22 SD (95% CI 0.18-0.25; 12.5%; p=9.3 x 10(-33)), concentrations of interleukin 6 decreased by 0.02 SD (-0.04 to -0.01; -1,7%; p=3.5 x 10(-3)), and concentrations of C-reactive protein decreased by 0.03 SD (-0.04 to -0.02; -3.4%; p=7.7 x 10(-14)). We noted the effects of the genetic score on these inflammation biomarkers to be directionally concordant with those of anakinra. The allele count of the genetic score had roughly log-linear, dose-dependent associations with both IL-1Ra concentration and risk of coronary heart disease. For people who carried four IL-1Ra-raising alleles, the odds ratio for coronary heart disease was 1.15 (1.08-1.22; p=1.8 x 10(-6)) compared with people who carried no IL-1Ra-raising alleles; the per-allele odds ratio for coronary heart disease was 1.03 (1.02-1.04; p=3.9 x 10(-10)). Perallele odds ratios were 0.97 (0.95-0.99; p=9.9 x 10(-4)) for rheumatoid arthritis, 0.99 (0.97-1.01; p=0.47) for type 2 diabetes, 1.00 (0.98-1.02; p=0.92) for ischaemic stroke, and 1.08 (1.04-1.12; p=1.8 x 10(-5)) for abdominal aortic aneurysm. In exploratory analyses, we observed per-allele increases in concentrations of proatherogenic lipids, including LDL-cholesterol, but no clear evidence of association for blood pressure, glycaemic traits, or any of the 24 other disorders studied. Modelling suggested that the observed increase in LDL-cholesterol could account for about a third of the association observed between the genetic score and increased coronary risk. Interpretation Human genetic data suggest that long-term dual IL-1 alpha/beta inhibition could increase cardiovascular risk and, conversely, reduce the risk of development of rheumatoid arthritis. The cardiovascular risk might, in part, be mediated through an increase in proatherogenic lipid concentrations. Copyright (C) The Interleukin 1 Genetics Consortium. Open Access article distributed under the terms of CC-BY-NC-ND

The genetic epidemiology of joint shape and the development of osteoarthritis

Congruent, low-friction relative movement between the articulating elements of a synovial joint is an essential pre-requisite for sustained, efficient, function. Where disorders of joint formation or maintenance exist, mechanical overloading and osteoarthritis (OA) follow. The heritable component of OA accounts for ~ 50% of susceptible risk. Although almost 100 genetic risk loci for OA have now been identified, and the epidemiological relationship between joint development, joint shape and osteoarthritis is well established, we still have only a limited understanding of the contribution that genetic variation makes to joint shape and how this modulates OA risk. In this article, a brief overview of synovial joint development and its genetic regulation is followed by a review of current knowledge on the genetic epidemiology of established joint shape disorders and common shape variation. A summary of current genetic epidemiology of OA is also given, together with current evidence on the genetic overlap between shape variation and OA. Finally, the established genetic risk loci for both joint shape and osteoarthritis are discussed

No evidence of an association between mitochondrial DNA variants and osteoarthritis in 7393 cases and 5122 controls.

OBJECTIVES: Osteoarthritis (OA) has a complex aetiology with a strong genetic component. Genome-wide association studies implicate several nuclear genes in the aetiology, but a major component of the heritability has yet to be defined at the molecular level. Initial studies implicate maternally inherited variants of mitochondrial DNA (mtDNA) in subgroups of patients with OA based on gender and specific joint involvement, but these findings have not been replicated. METHODS: The authors studied 138 maternally inherited mtDNA variants genotyped in a two cohort genetic association study across a total of 7393 OA cases from the arcOGEN consortium and 5122 controls genotyped in the Wellcome Trust Case Control consortium 2 study. RESULTS: Following data quality control we examined 48 mtDNA variants that were common in cohort 1 and cohort 2, and found no association with OA. None of the phenotypic subgroups previously associated with mtDNA haplogroups were associated in this study. CONCLUSIONS: We were not able to replicate previously published findings in the largest mtDNA association study to date. The evidence linking OA to mtDNA is not compelling at present