The genetic architecture of osteoarthritis: insights from UK Biobank

Osteoarthritis is a common complex disease with huge public health burden. Here we perform a genome-wide association study for osteoarthritis using data across 16.5 million variants from the UK Biobank resource. Following replication and meta-analysis in up to 30,727 cases and 297,191 controls, we report 9 new osteoarthritis loci, in all of which the most likely causal variant is non-coding. For three loci, we detect association with biologically-relevant radiographic endophenotypes, and in five signals we identify genes that are differentially expressed in degraded compared to intact articular cartilage from osteoarthritis patients. We establish causal effects for higher body mass index, but not for triglyceride levels or type 2 diabetes liability, on osteoarthritis.

[1]  P. Donnelly,et al.  A new multipoint method for genome-wide association studies by imputation of genotypes , 2007, Nature Genetics.

[2]  J. Mort,et al.  Rho/ROCK and MEK/ERK activation by transforming growth factor-α induces articular cartilage degradation , 2010, Laboratory Investigation.

[3]  Wj Gauderman,et al.  QUANTO 1.1: A computer program for power and sample size calculations for genetic-epidemiology studies , 2006 .

[4]  M. Daly,et al.  Proteins Encoded in Genomic Regions Associated with Immune-Mediated Disease Physically Interact and Suggest Underlying Biology , 2011, PLoS genetics.

[5]  E. Zeggini,et al.  Advances in osteoarthritis genetics , 2013, Journal of Medical Genetics.

[6]  Michael Boehnke,et al.  Recommended Joint and Meta‐Analysis Strategies for Case‐Control Association Testing of Single Low‐Count Variants , 2013, Genetic epidemiology.

[7]  Beatriz Gámez,et al.  p38 MAPK Signaling in Osteoblast Differentiation , 2016, Front. Cell Dev. Biol..

[8]  German Tischler,et al.  biobambam: tools for read pair collation based algorithms on BAM files , 2013, Source Code for Biology and Medicine.

[9]  R. Pearson,et al.  Histopathology grading systems for characterisation of human knee osteoarthritis--reproducibility, variability, reliability, correlation, and validity. , 2011, Osteoarthritis and cartilage.

[10]  Yin Xiao,et al.  Inhibition of p38 pathway leads to OA-like changes in a rat animal model. , 2012, Rheumatology.

[11]  E. Zeggini,et al.  Pathways to understanding the genomic aetiology of osteoarthritis , 2017, Human molecular genetics.

[12]  S. Ochiai,et al.  Oxidative stress reaction in the meniscus of Bach 1 deficient mice: Potential prevention of meniscal degeneration , 2008, Journal of orthopaedic research : official publication of the Orthopaedic Research Society.

[13]  M. Englund,et al.  Association between occupation and knee and hip replacement due to osteoarthritis: a case-control study , 2010, Arthritis research & therapy.

[14]  Shijie Wang,et al.  [Junctophilin, a molecule that couples endo (sarco) plasma membrane to cell membrane]. , 2009, Sheng li ke xue jin zhan [Progress in physiology].

[15]  G. Wiche,et al.  The many faces of plectin and plectinopathies: pathology and mechanisms , 2012, Acta Neuropathologica.

[16]  Eleni Zengini,et al.  The Genetic Epidemiological Landscape of Hip and Knee Osteoarthritis: Where Are We Now and Where Are We Going? , 2016, The Journal of Rheumatology.

[17]  Stephen Burgess,et al.  Sensitivity Analyses for Robust Causal Inference from Mendelian Randomization Analyses with Multiple Genetic Variants , 2016, Epidemiology.

[18]  Jun Li,et al.  Prevalence and associated factors of knee osteoarthritis in a rural Chinese adult population: an epidemiological survey , 2015, BMC Public Health.

[19]  D. A. van der Windt,et al.  Subsets of symptomatic hand osteoarthritis in community-dwelling older adults in the United Kingdom: prevalence, inter-relationships, risk factor profiles and clinical characteristics at baseline and 3-years , 2013, Osteoarthritis and cartilage.

[20]  C. Hiller,et al.  Self-reported knee pain and disability among healthy individuals: reference data and factors associated with the Knee injury and Osteoarthritis Outcome Score (KOOS) and KOOS-Child. , 2017, Osteoarthritis and Cartilage.

[21]  F. Cicuttini,et al.  How Are Obesity and Body Composition Related to Patellar Cartilage? A Systematic Review , 2017, The Journal of Rheumatology.

[22]  E. Zeggini,et al.  Integrative epigenomics, transcriptomics and proteomics of patient chondrocytes reveal genes and pathways involved in osteoarthritis , 2016, Scientific Reports.

[23]  M. Robinson,et al.  Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. , 2015, F1000Research.

[24]  Yusuke Nakamura,et al.  A functional polymorphism in the 5′ UTR of GDF5 is associated with susceptibility to osteoarthritis , 2007, Nature Genetics.

[25]  Kari Stefansson,et al.  Severe osteoarthritis of the hand associates with common variants within the ALDH1A2 gene and with rare variants at 1p31 , 2014, Nature Genetics.

[26]  Yu Sun,et al.  Bioinformatics analysis of gene expression profiles of osteoarthritis. , 2015, Acta histochemica.

[27]  Andres Metspalu,et al.  Identification of new susceptibility loci for osteoarthritis (arcOGEN): a genome-wide association study , 2012, The Lancet.

[28]  A. Butterworth,et al.  Mendelian Randomization Analysis With Multiple Genetic Variants Using Summarized Data , 2013, Genetic epidemiology.

[29]  A. Iagnocco,et al.  Diabetes Is an Independent Predictor for Severe Osteoarthritis , 2013, Diabetes Care.

[30]  C. Cooper,et al.  An increased rate of falling leads to a rise in fracture risk in postmenopausal women with self-reported osteoarthritis: a prospective multinational cohort study (GLOW) , 2012, Annals of the rheumatic diseases.

[31]  A. Uitterlinden,et al.  The DOT1L rs12982744 polymorphism is associated with osteoarthritis of the hip with genome-wide statistical significance in males , 2013, Annals of the rheumatic diseases.

[32]  John Loughlin,et al.  An SNP in the 5'-UTR of GDF5 is associated with osteoarthritis susceptibility in Europeans and with in vivo differences in allelic expression in articular cartilage. , 2007, Human molecular genetics.

[33]  J. Loughlin.  The genetic epidemiology of human primary osteoarthritis: current status , 2005, Expert Reviews in Molecular Medicine.

[34]  Yusuke Nakamura,et al.  New Sequence Variants in HLA Class II/III Region Associated with Susceptibility to Knee Osteoarthritis Identified by Genome-Wide Association Study , 2010, PloS one.

[35]  Francisco S. Roque,et al.  A large-scale analysis of tissue-specific pathology and gene expression of human disease genes and complexes , 2008, Proceedings of the National Academy of Sciences.

[36]  Bjarni V. Halldórsson,et al.  A genome-wide association study identifies an osteoarthritis susceptibility locus on chromosome 7q22. , 2010, Arthritis and rheumatism.

[37]  Tom R. Gaunt,et al.  The UK10K project identifies rare variants in health and disease , 2016 .

[38]  Ricardo Pong-Wong,et al.  Evaluating the contribution of genetic and familial shared environment to common disease using the UK Biobank , 2016, Nature Genetics.

[39]  A. Hofman,et al.  A variant in MCF2L is associated with osteoarthritis. , 2011, American journal of human genetics.

[40]  Lincoln Stein,et al.  Reactome: a knowledgebase of biological pathways , 2004, Nucleic Acids Res..

[41]  R. Derynck,et al.  Epidermal growth factor receptor function is necessary for normal craniofacial development and palate closure , 1999, Nature Genetics.

[42]  J. Buxbaum,et al.  A SPECTRAL APPROACH INTEGRATING FUNCTIONAL GENOMIC ANNOTATIONS FOR CODING AND NONCODING VARIANTS , 2015, Nature Genetics.

[43]  Tom R. Gaunt,et al.  LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis , 2016, bioRxiv.

[44]  Gunwoo Kim,et al.  Transforming growth factor alpha controls the transition from hypertrophic cartilage to bone during endochondral bone growth. , 2012, Bone.

[45]  K. Muir,et al.  Home based exercise programme for knee pain and knee osteoarthritis: randomised controlled trial , 2002, BMJ : British Medical Journal.

[46]  Gabor T. Marth,et al.  A global reference for human genetic variation , 2015, Nature.

[47]  Matthew E. Ritchie,et al.  limma powers differential expression analyses for RNA-sequencing and microarray studies , 2015, Nucleic acids research.

[48]  G. Hawker,et al.  Impact of Osteoarthritis on Difficulty Walking: A Population‐Based Study , 2018, Arthritis care & research.

[49]  Hiroyuki Ogata,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 1999, Nucleic Acids Res..

[50]  Zoltán Kutalik,et al.  Quality control and conduct of genome-wide association meta-analyses , 2014, Nature Protocols.

[51]  L. Li,et al.  Junctophilin 3 expresses in pancreatic beta cells and is required for glucose-stimulated insulin secretion , 2016, Cell Death and Disease.

[52]  H. Dorfman,et al.  Biochemical and metabolic abnormalities in articular cartilage from osteo-arthritic human hips. II. Correlation of morphology with biochemical and metabolic data. , 1971, The Journal of bone and joint surgery. American volume.

[53]  N. Sheehan,et al.  Assessing the suitability of summary data for two-sample Mendelian randomization analyses using MR-Egger regression: the role of the I2 statistic , 2016, International journal of epidemiology.

[54]  A. Hofman,et al.  A meta-analysis of genome-wide association studies identifies novel variants associated with osteoarthritis of the hip , 2013, Annals of the rheumatic diseases.

[55]  Ying Zhang,et al.  A review of crosstalk between MAPK and Wnt signals and its impact on cartilage regeneration , 2014, Cell and Tissue Research.

[56]  Pall I. Olason,et al.  A human phenome-interactome network of protein complexes implicated in genetic disorders , 2007, Nature Biotechnology.

[57]  D. Lawlor,et al.  Re: Estimation of bias in nongenetic observational studies using "Mendelian triangulation" by Bautista et al. , 2007, Annals of epidemiology.

[58]  Qian Chen,et al.  CD95-induced osteoarthritic chondrocyte apoptosis and necrosis: dependency on p38 mitogen-activated protein kinase , 2006, Arthritis research & therapy.

[59]  D. Shi,et al.  A meta-analysis of European and Asian cohorts reveals a global role of a functional SNP in the 5' UTR of GDF5 with osteoarthritis susceptibility. , 2008, Human molecular genetics.

[60]  George Davey Smith,et al.  Mendelian randomization: Using genes as instruments for making causal inferences in epidemiology , 2008, Statistics in medicine.

[61]  E. Tai,et al.  Genome-wide association studies in East Asians identify new loci for waist-hip ratio and waist circumference , 2016, Scientific Reports.

[62]  Tamara S. Roman,et al.  New genetic loci link adipose and insulin biology to body fat distribution , 2014, Nature.

[63]  G. Davey Smith,et al.  Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator , 2016, Genetic epidemiology.

[64]  E. Zeggini,et al.  Evaluation of shared genetic aetiology between osteoarthritis and bone mineral density identifies SMAD3 as a novel osteoarthritis risk locus , 2017, Human molecular genetics.

[65]  F. Hartwig,et al.  Why internal weights should be avoided (not only) in MR-Egger regression. , 2016, International journal of epidemiology.

[66]  David T. Broome,et al.  Mitogen-activated protein kinase phosphatase-1: function and regulation in bone and related tissues , 2016, Connective tissue research.

[67]  Qian Chen,et al.  Mitogen-activated protein kinase p38 induces HDAC4 degradation in hypertrophic chondrocytes. , 2015, Biochimica et biophysica acta.

[68]  G. Davey Smith,et al.  Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression , 2015, International journal of epidemiology.

[69]  D. Moore,et al.  Reduced limb length and worsened osteoarthritis in adult mice after genetic inhibition of p38 MAP kinase activity in cartilage. , 2008, Arthritis and rheumatism.

[70]  Debbie A Lawlor,et al.  Triangulation in aetiological epidemiology , 2016, International journal of epidemiology.

[71]  N. Sheehan,et al.  A framework for the investigation of pleiotropy in two‐sample summary data Mendelian randomization , 2017, Statistics in medicine.

[72]  Tom R. Gaunt,et al.  Improved imputation of low-frequency and rare variants using the UK10K haplotype reference panel , 2015, Nature Communications.

[73]  Joris M. Mooij,et al.  MAGMA: Generalized Gene-Set Analysis of GWAS Data , 2015, PLoS Comput. Biol..

[74]  Yusuke Nakamura,et al.  Common variants in DVWA on chromosome 3p24.3 are associated with susceptibility to knee osteoarthritis , 2008, Nature Genetics.

[75]  Yun Li,et al.  METAL: fast and efficient meta-analysis of genomewide association scans , 2010, Bioinform..

[76]  Bjarni V. Halldórsson,et al.  Meta-analysis of genome-wide association studies confirms a susceptibility locus for knee osteoarthritis on chromosome 7q22 , 2010, Annals of the rheumatic diseases.

[77]  M. Doherty,et al.  Assessment of a genetic contribution to osteoarthritis of the hip: sibling study , 2000, BMJ : British Medical Journal.

[78]  T. Moilanen,et al.  Adiponectin associates with markers of cartilage degradation in osteoarthritis and induces production of proinflammatory and catabolic factors through mitogen-activated protein kinase pathways , 2011, Arthritis research & therapy.

[79]  L. Nordsletten,et al.  Distribution of osteoarthritis in a Norwegian population-based cohort: associations to risk factor profiles and health-related quality of life , 2017, Rheumatology International.

[80]  Hashem A. Shihab,et al.  MR-Base: a platform for systematic causal inference across the phenome using billions of genetic associations , 2016, bioRxiv.

[81]  T. Spector,et al.  Osteoarthritis: New Insights. Part 1: The Disease and Its Risk Factors , 2000, Annals of Internal Medicine.

[82]  Shankung Lin,et al.  Zinc finger protein ZFP36L1 promotes osteoblastic differentiation but represses adipogenic differentiation of mouse multipotent cells , 2017, Oncotarget.

[83]  Brian V. Jenkins,et al.  A mouse forward genetics screen identifies LISTERIN as an E3 ubiquitin ligase involved in neurodegeneration , 2009, Proceedings of the National Academy of Sciences.

[84]  A. Rosenthal,et al.  The adverse effects of diabetes on osteoarthritis: update on clinical evidence and molecular mechanisms. , 2015, Osteoarthritis and cartilage.

[85]  Ross M. Fraser,et al.  Genetic studies of body mass index yield new insights for obesity biology , 2015, Nature.

[86]  Q. Jin,et al.  Association between the +104T/C polymorphism in the 5'UTR of GDF5 and susceptibility to knee osteoarthritis: a meta-analysis. , 2013, Molecular medicine reports.

[87]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[88]  M. Ashburner,et al.  Gene Ontology: tool for the unification of biology , 2000, Nature Genetics.

[89]  A. Hofman,et al.  Novel Genetic Variants for Cartilage Thickness and Hip Osteoarthritis , 2016, PLoS genetics.

[90]  Susumu Goto,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 2000, Nucleic Acids Res..

[91]  U. Sambamoorthi,et al.  The Association of Perceived Memory Loss with Osteoarthritis and Related Joint Pain in a Large Appalachian Population , 2018, Pain medicine.

[92]  M. Daly,et al.  LD Score regression distinguishes confounding from polygenicity in genome-wide association studies , 2014, Nature Genetics.

[93]  Yin Xiao,et al.  ERK-1/2 and p38 in the regulation of hypertrophic changes of normal articular cartilage chondrocytes induced by osteoarthritic subchondral osteoblasts. , 2010, Arthritis and rheumatism.

[94]  Richard C. Laughlin,et al.  Cellular and molecular mechanisms that underlie Entamoeba histolytica pathogenesis: prospects for intervention , 2005, Expert Reviews in Molecular Medicine.

[95]  R. Coffey,et al.  From wavy hair to naked proteins: the role of transforming growth factor alpha in health and disease. , 2014, Seminars in cell & developmental biology.

[96]  M. Robinson,et al.  Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences , 2015, F1000Research.

[97]  S. Murakami,et al.  Constitutive activation of MKK6 in chondrocytes of transgenic mice inhibits proliferation and delays endochondral bone formation. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[98]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[99]  T. Aigner,et al.  Transforming growth factor alpha suppression of articular chondrocyte phenotype and Sox9 expression in a rat model of osteoarthritis. , 2007, Arthritis and rheumatism.

[100]  Paul Dieppe,et al.  Osteoarthritis: New Insights. Part 2: Treatment Approaches , 2000, Annals of Internal Medicine.

[101]  A. Hofman,et al.  Cam Deformity and Acetabular Dysplasia as Risk Factors for Hip Osteoarthritis , 2017, Arthritis & rheumatology.

[102]  Ross M. Fraser,et al.  Defining the role of common variation in the genomic and biological architecture of adult human height , 2014, Nature Genetics.

[103]  P. Casey,et al.  Protein prenylation: unique fats make their mark on biology , 2016, Nature Reviews Molecular Cell Biology.

[104]  P. Elliott,et al.  UK Biobank: An Open Access Resource for Identifying the Causes of a Wide Range of Complex Diseases of Middle and Old Age , 2015, PLoS medicine.

[105]  G. Davey Smith,et al.  Epidemiology--is it time to call it a day? , 2001, International journal of epidemiology.

[106]  J. Shendure,et al.  A general framework for estimating the relative pathogenicity of human genetic variants , 2014, Nature Genetics.

[107]  V. Rosen,et al.  Human Migratory Meniscus Progenitor Cells Are Controlled via the TGF-β Pathway , 2014, Stem cell reports.

[108]  J. Horwood.  UK Biobank Data: Come and Get It , 2014 .

[109]  Chang-jie Zhang,et al.  Explore the variation of MMP3, JNK, p38 MAPKs, and autophagy at the early stage of osteoarthritis , 2016, IUBMB life.

[110]  H. Stefánsson,et al.  Whole-genome sequencing identifies rare genotypes in COMP and CHADL associated with high risk of hip osteoarthritis , 2017, Nature Genetics.

[111]  Nuala A Sheehan,et al.  Detecting pleiotropy in Mendelian randomisation studies with summary data and a continuous outcome , 2015, Statistics in medicine.

[112]  Jonathan A C Sterne,et al.  The impact of residual and unmeasured confounding in epidemiologic studies: a simulation study. , 2007, American journal of epidemiology.

[113]  Ming D. Li,et al.  Genome-wide meta-analyses identify multiple loci associated with smoking behavior , 2010, Nature Genetics.

[114]  G. Walko,et al.  Plectin–intermediate filament partnership in skin, skeletal muscle, and peripheral nerve , 2013, Histochemistry and Cell Biology.

[115]  Alessandro Vullo,et al.  Ensembl 2017 , 2016, Nucleic Acids Res..

[116]  Michael A Mont,et al.  What is the Impact of Body Mass Index on Cardiovascular and Musculoskeletal Health? , 2017, Surgical technology international.

[117]  Henning Hermjakob,et al.  The Reactome pathway knowledgebase , 2013, Nucleic Acids Res..

[118]  A. Palmieri,et al.  Comparison between osteoblasts derived from human dental pulp stem cells and osteosarcoma cell lines , 2008, Cell biology international.

[119]  S. Goldring,et al.  Osteoarthritis: a disease of the joint as an organ. , 2012, Arthritis and rheumatism.

[120]  I. Welch,et al.  Context-specific protection of TGFα null mice from osteoarthritis , 2016, Scientific Reports.

[121]  P. Korkusuz,et al.  Mutation in exon 1f of PLEC, leading to disruption of plectin isoform 1f, causes autosomal-recessive limb-girdle muscular dystrophy. , 2010, American journal of human genetics.

[122]  Chi Ma,et al.  Pioglitazone inhibits advanced glycation end product-induced matrix metalloproteinases and apoptosis by suppressing the activation of MAPK and NF-κB , 2016, Apoptosis.

[123]  Jie Huang,et al.  Discovery and refinement of genetic loci associated with cardiometabolic risk using dense imputation maps , 2016, Nature Genetics.

[124]  S. Ebrahim,et al.  'Mendelian randomization': can genetic epidemiology contribute to understanding environmental determinants of disease? , 2003, International journal of epidemiology.

[125]  Charity W. Law,et al.  voom: precision weights unlock linear model analysis tools for RNA-seq read counts , 2014, Genome Biology.

[126]  T. Spector,et al.  Genetic epidemiology of hip and knee osteoarthritis , 2011, Nature Reviews Rheumatology.

[127]  R. Altman.  Early management of osteoarthritis. , 2010, The American journal of managed care.

[128]  Gregory A. Poland,et al.  Fine Mapping Causal Variants with an Approximate Bayesian Method Using Marginal Test Statistics , 2015, Genetics.

[129]  Yin Xiao,et al.  Osteoarthritic cartilage chondrocytes alter subchondral bone osteoblast differentiation via MAPK signalling pathway involving ERK1/2. , 2010, Bone.

[130]  N. Eriksson,et al.  Genome-wide analysis identifies 12 loci influencing human reproductive behavior , 2016 .

[131]  J. Jordan,et al.  Associations of educational attainment, occupation and community poverty with knee osteoarthritis in the Johnston County (North Carolina) osteoarthritis project , 2011, Arthritis research & therapy.

[132]  Fernando Pires Hartwig,et al.  Robust inference in two-sample Mendelian randomisation via the zero modal pleiotropy assumption , 2017, bioRxiv.

[133]  J. Protheroe,et al.  Current evidence on risk factors for knee osteoarthritis in older adults: a systematic review and meta-analysis. , 2015, Osteoarthritis and cartilage.

[134]  Joris M. Mooij,et al.  MAGMA: Generalized Gene-Set Analysis of GWAS Data , 2015, PLoS Comput. Biol..

[135]  M. Lotz,et al.  Bach1 deficiency reduces severity of osteoarthritis through upregulation of heme oxygenase-1 , 2015, Arthritis Research & Therapy.

[136]  Rob Patro,et al.  Salmon provides fast and bias-aware quantification of transcript expression , 2017, Nature Methods.

[137]  F. Berenbaum,et al.  Association between diabetes mellitus and osteoarthritis: systematic literature review and meta-analysis , 2015, RMD Open.