Genome-wide analyses using UK Biobank data provide insights into the genetic architecture of osteoarthritis

Human gene and disease associations for clinical-genomics and precision medicine research

We are entering the era of personalized medicine in which an individual's genetic makeup will eventually determine how a doctor can tailor his or her therapy. Therefore, it is becoming critical to understand the genetic basis of common diseases, for example, which genes predispose and rare genetic variants contribute to diseases, and so on. Our study focuses on helping researchers, medical practitioners, and pharmacists in having a broad view of genetic variants that may be implicated in the likelihood of developing certain diseases. Our focus here is to create a comprehensive database with mobile access to all available, authentic and actionable genes, SNPs, and classified diseases and drugs collected from different clinical and genomics databases worldwide, including Ensembl, GenCode, ClinVar, GeneCards, DISEASES, HGMD, OMIM, GTR, CNVD, Novoseek, Swiss-Prot, LncRNADisease, Orphanet, GWAS Catalog, SwissVar, COSMIC, WHO, and FDA. We present a new cutting-edge gene-SNP-disease-drug mobile database with a smart phone application, integrating information about classified diseases and related genes, germline and somatic mutations, and drugs. Its database includes over 59 000 protein-coding and noncoding genes; over 67 000 germline SNPs and over a million somatic mutations reported for over 19 000 protein-coding genes located in over 1000 regions, published with over 3000 articles in over 415 journals available at the PUBMED; over 80 000 ICDs; over 123 000 NDCs; and over 100 000 classified gene-SNP-disease associations. We present an application that can provide new insights into the information about genetic basis of human complex diseases and contribute to assimilating genomic with phenotypic data for the availability of gene-based designer drugs, precise targeting of molecular fingerprints for tumor, appropriate drug therapy, predicting individual susceptibility to disease, diagnosis, and treatment of rare illnesses are all a few of the many transformations expected in the decade to come.

Unpicking observational relationships between hip shape and osteoarthritis: hype or hope?

PURPOSE OF REVIEW

To review recent findings concerning the observational relationship between hip shape and hip osteoarthritis (HOA) and their shared genetic influences, and the potential for clinical application.

RECENT FINDINGS

Recent observational studies have strengthened the evidence that specific shape deformities, such as cam and acetabular dysplasia, are related to HOA. Statistical shape modelling has emerged as a method to measure hip shape holistically, with the added advantage that this can be applied to dual X-ray absorptiometry scan images. This has led to several additional aspects of hip shape variation being identified, such as a wider femoral neck and larger lesser trochanter, in association with HOA. Furthermore, this method has formed the basis of genetic studies identifying novel genetic influences on hip shape, several of which are shared with known genetic risk factors for HOA.

SUMMARY

Shared genetic influences of hip shape and HOA raise the possibility that hip shape plays a casual role in the development of HOA, justifying preventive approaches aiming to combat these adverse consequences.

Discovery and analysis of methylation quantitative trait loci (mQTLs) mapping to novel osteoarthritis genetic risk signals

OBJECTIVE

Osteoarthritis (OA) is polygenic with over 90 independent genome-wide association loci so far reported. A key next step is the identification of target genes and the molecular mechanisms through which this genetic risk operates. The majority of OA risk-conferring alleles are predicted to act by modulating gene expression. DNA methylation at CpG dinucleotides may be a functional conduit through which this occurs and is detectable by mapping methylation quantitative trait loci, or mQTLs. This approach can therefore provide functional insight into OA risk and will prioritize genes for subsequent investigation. That was our goal, with a focus on the largest set of OA loci yet to be reported.

METHOD

We investigated DNA methylation, genotype and RNA sequencing data derived from the cartilage of patients who had undergone arthroplasty and combined this with in silico analyses of expression quantitative trait loci, epigenomes and chromatin interactions.

RESULTS

We investigated 42 OA risk loci and in ten of these we identified 24 CpGs in which methylation correlated with genotype (false discovery rate (FDR) P-values ranging from 0.049 to 1.73x10^-25). In silico analyses of these mQTLs prioritised genes and regulatory elements at the majority of the ten loci, with COLGALT2 (encoding a collagen galactosyltransferase), COL11A2 (encoding a polypeptide chain of type XI collagen) and WWP2 (encoding a ubiquitin ligase active during chondrogenesis) emerging as particularly compelling target genes.

CONCLUSION

We have highlighted the pivotal role of DNA methylation as a link between genetic risk and OA and prioritized genes for further investigation.

Urate and osteoarthritis: Evidence for a reciprocal relationship

Hyperuricemia is a common condition, and in a subset of patients leads to gout, the most common inflammatory arthritis. Osteoarthritis is the most common form of arthritis overall, and gout and osteoarthritis frequently coexist in the same patient. However, the relationship between the two remains poorly defined. More particularly, the impact of osteoarthritis on the development of gout, and the impact of gout on the development of osteoarthritis, remain to be determined. Additionally, whether hyperuricemia mediates osteoarthritis in the absence of gout is uncertain. Here, we review the evidence linking gout and osteoarthritis, with a special focus on the role of hyperuricemia in the presence or absence of gout. Since disease modifying agents are currently available for hyperuricemia and gout but not for osteoarthritis, a contributory role for urate in the pathogenesis of osteoarthritis could have important clinical implications.

Integrating Genome-Wide Association Studies With Pathway Analysis and Gene Expression Analysis Highlights Novel Osteoarthritis Risk Pathways and Genes

Osteoarthritis (OA) is the most common degenerative joint disorder worldwide. To identify more genetic signals, genome-wide association study (GWAS) has been widely used and elucidated some OA susceptibility genes. However, these susceptibility genes could only explain only a small part of heritability of OA. It is suggested that the identification of disease-related pathways may contribute to understand the genomic etiology of OA. Here, we integrated the GWAS into pathway analysis to identify novel OA risk pathways. In this study, we first selected 187 independent genetic variants identified by GWAS (P < 1.00E−05) and found that most of these genetic variants are noncoding mutations. We then conducted an expression quantitative trait loci analysis and found that 165 of these 187 genetic variants could significantly regulate the expression of nearby genes. Third, we identified OA susceptibility genes corresponding to these genetic variants, conducted a pathway analysis, and identified novel OA-related KEGG pathways, GO biological processes, GO molecular functions, and GO cellular components. In KEGG database, transforming growth factor β signaling pathway is the most significant signal (P = 5.98E−05) and is the only pathway after the BH multiple-test adjustment with false discovery rate (FDR) = 0.02. In GO database, we identified 24 statistically significant GO biological processes, one statistically significant GO molecular function, and five statistically significant GO cellular components (FDR < 0.05). These signals are related with chondrocyte differentiation and development, which are all known biological pathways associated with OA. Finally, we conducted an OA case–control gene expression analysis to evaluate the differential expression of these OA risk genes. Using an OA case–control gene expression analysis, we showed that 44 risk genes were suggestively differentially expressed in OA cases compared with controls (P < 0.05). Three genes, WWP2, COG5, and MAPT, were statistically differentially expressed in OA cases compared with controls (P < 0.05/122 = 4.10E−04). Hence, our findings may contribute to understanding the genomic etiology of OA.

Causal Factors for Knee, Hip, and Hand Osteoarthritis: A Mendelian Randomization Study in the UK Biobank

OBJECTIVE

There is no curative treatment for osteoarthritis (OA), which is the most common form of arthritis. This study was undertaken to identify causal risk factors of knee, hip, and hand OA.

METHODS

Individual-level data from 384,838 unrelated participants in the UK Biobank study were analyzed. Mendelian randomization (MR) analyses were performed to test for causality for body mass index (BMI), bone mineral density (BMD), serum high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, and triglyceride levels, type 2 diabetes, systolic blood pressure (BP), and C-reactive protein (CRP) levels. The primary outcome measure was OA determined using hospital diagnoses (all sites, n = 48,431; knee, n = 19,727; hip, n = 11,875; hand, n = 2,330). Odds ratios (ORs) with 95% confidence intervals (95% CIs) were calculated.

RESULTS

MR analyses demonstrated a robust causal association of genetically determined BMI with all OA (OR per SD increase 1.57 [95% CI 1.44-1.71]), and with knee OA and hip OA, but not with hand OA. Increased genetically determined femoral neck BMD was causally associated with all OA (OR per SD increase 1.14 [95% CI 1.06-1.22]), knee OA, and hip OA. Low systolic BP was causally associated with all OA (OR per SD decrease 1.55 [95% CI 1.29-1.87]), knee OA, and hip OA. There was no evidence of causality for the other tested metabolic factors or CRP level.

CONCLUSION

Our findings indicate that BMI exerts a major causal effect on the risk of OA at weight-bearing joints, but not at the hand. Evidence of causality of all OA, knee OA, and hip OA was also observed for high femoral neck BMD and low systolic BP. However, we found no evidence of causality for other metabolic factors or CRP level.

Genome-wide association study of knee pain identifies associations with GDF5 and COL27A1 in UK Biobank

Knee pain is one of the most common musculoskeletal complaints that brings people to medical attention. Approximately 50% of individuals over the age of 50 report an experience of knee pain within the past 12 months. We sought to identify the genetic variants associated with knee pain in 171,516 subjects from the UK Biobank cohort and seek supporting evidence in cohorts from 23andMe, the Osteoarthritis Initiative, and the Johnston County Osteoarthritis Project. We identified two loci that reached genome-wide significance in the UK Biobank: rs143384, located in GDF5 (P = 1.32 × 10-12), a gene previously implicated in osteoarthritis; and rs2808772, located near COL27A1 (P = 1.49 × 10-8). These findings were supported in cohorts with self-reported osteoarthritis/radiographic knee osteoarthritis without pain information. In this report on genome-wide association of knee pain, we identified two loci in or near GDF5 and COL27A1 that are associated with knee pain.

Genome Engineering for Osteoarthritis: From Designer Cells to Disease-Modifying Drugs

BACKGROUND

Osteoarthritis (OA) is a highly prevalent degenerative joint disease involving joint cartilage and its surrounding tissues. OA is the leading cause of pain and disability worldwide. At present, there are no disease-modifying OA drugs, and the primary therapies include exercise and nonsteroidal anti-inflammatory drugs until total joint replacement at the end-stage of the disease.

METHODS

In this review, we summarized the current state of knowledge in genetic and epigenetic associations and risk factors for OA and their potential diagnostic and therapeutic applications.

RESULTS

Genome-wide association studies and analysis of epigenetic modifications (such as miRNA expression, DNA methylation and histone modifications) conducted across various populations support the notion that there is a genetic basis for certain subsets of OA pathogenesis.

CONCLUSION

With recent advances in the development of genome editing technologies such as the CRISPR-Cas9 system, these genetic and epigenetic alternations in OA can be used as platforms from which potential biomarkers for the diagnosis, prognosis, drug response, and development of potential personalized therapeutic targets for OA can be approached. Furthermore, genome editing has allowed the development of "designer" cells, whereby the receptors, gene regulatory networks, or transgenes can be modified as a basis for new cell-based therapies.

Prioritization of PLEC and GRINA as Osteoarthritis Risk Genes Through the Identification and Characterization of Novel Methylation Quantitative Trait Loci

OBJECTIVE

To identify methylation quantitative trait loci (mQTLs) correlating with osteoarthritis (OA) risk alleles and to undertake mechanistic characterization as a means of target gene prioritization.

METHODS

We used genome-wide genotyping and cartilage DNA methylation array data in a discovery screen of novel OA risk loci. This was followed by methylation, gene expression analysis, and genotyping studies in additional cartilage samples, accompanied by in silico analyses.

RESULTS

We identified 4 novel OA mQTLs. The most significant mQTL contained 9 CpG sites where methylation correlated with OA risk genotype, with 5 of the CpG sites having P values <1 × 10^-10 . The 9 CpG sites reside in an interval of only 7.7 kb within the PLEC gene and form 2 distinct clusters. We were able to prioritize PLEC and the adjacent gene GRINA as independent targets of the OA risk. We identified PLEC and GRINA expression QTLs operating in cartilage, as well as methylation-expression QTLs operating on the 2 genes. GRINA and PLEC also demonstrated differential expression between OA hip and non-OA hip cartilage.

CONCLUSION

PLEC encodes plectin, a cytoskeletal protein that maintains tissue integrity by regulating intracellular signaling in response to mechanical stimuli. GRINA encodes the ionotropic glutamate receptor TMBIM3 (transmembrane BAX inhibitor 1 motif-containing protein family member 3), which regulates cell survival. Based on our results, we hypothesize that in a joint predisposed to OA, expression of these genes alters in order to combat aberrant biomechanics, and that this is epigenetically regulated. However, carriage of the OA risk-conferring allele at this locus hinders this response and contributes to disease development.

Chondroprotective Factors in Osteoarthritis: a Joint Affair

PURPOSE OF THE REVIEW

Osteoarthritis is widely regarded as a spectrum of conditions that affect all joint tissues, typified by a common entity: cartilage loss. Here, we review recent progress and challenges in chondroprotection and discuss new strategies to prevent cartilage loss in osteoarthritis.

RECENT FINDINGS

Advances in clinical, molecular, and cellular characterization are enabling improved stratification of osteoarthritis subtypes. Integration of next-generation sequencing and "omics" approaches with clinically relevant readouts shows promise in delineating both subtypes of disease and meaningful trial end points. Novel delivery strategies are enabling joint-specific delivery. Chondroprotection requires a whole joint approach, stratification of patient groups, and use of patient-relevant end points. Drug development should continue to explore new targets, while using modern technologies and recent knowledge to re-visit unsuccessful therapeutics from the past. The overarching goal for chondroprotection is to provide the right treatment(s) for the right patient at the right time.

Expression analysis of the osteoarthritis genetic susceptibility mapping to the matrix Gla protein gene MGP

BACKGROUND

Osteoarthritis (OA) is a common disease of older individuals that impacts detrimentally on the quality and the length of life. It is characterised by the painful loss of articular cartilage and is polygenic and multifactorial. Genome-wide association scans have highlighted over 90 osteoarthritis genetic signals, some of which reside within or close to highly plausible candidate genes. An example is an association to polymorphisms within and adjacent to the matrix Gla protein gene MGP. We set out to undertake a functional study of this gene.

METHODS

Nucleic acid was extracted from cartilage, infrapatellar fat pad, synovium, trabecular bone, trapezium and peripheral whole blood from OA patients and also from mesenchymal stem cells (MSCs) subjected to chondrogenesis. Expression of MGP was measured by quantitative PCR (qPCR), RNA-sequencing and allelic expression imbalance (AEI) analysis. Matrix Gla protein was depleted in chondrocytes by knocking down MGP expression using RNA interference (RNAi) and the effect on a range of genes assessed by qPCR.

RESULTS

MGP is expressed in joint tissues, blood and chondrocytes cultured from MSCs. There is a higher expression in diseased versus non-diseased cartilage. Polymorphisms that are associated with OA also correlate with the expression of MGP, with the OA risk-conferring allele showing significantly reduced expression in cartilage, fat pad and synovium but increased expression in blood. Depletion of Matrix Gla protein had a significant effect on the majority of genes tested, with an increased expression of catabolic genes that encode enzymes that degrade cartilage.

CONCLUSIONS

MGP expression is subject to cis-acting regulators that correlate with the OA association signal. These are active in a range of joint tissues but have effects which are particularly strong in cartilage. An opposite effect is observed in blood, highlighting the context-specific nature of the regulation of this gene's expression. Recapitulation of the genetic deficit in cartilage chondrocytes is pro-catabolic.

MOB (Mps one Binder) Proteins in the Hippo Pathway and Cancer

The family of MOBs (monopolar spindle-one-binder proteins) is highly conserved in the eukaryotic kingdom. MOBs represent globular scaffold proteins without any known enzymatic activities. They can act as signal transducers in essential intracellular pathways. MOBs have diverse cancer-associated cellular functions through regulatory interactions with members of the NDR/LATS kinase family. By forming additional complexes with serine/threonine protein kinases of the germinal centre kinase families, other enzymes and scaffolding factors, MOBs appear to be linked to an even broader disease spectrum. Here, we review our current understanding of this emerging protein family, with emphases on post-translational modifications, protein-protein interactions, and cellular processes that are possibly linked to cancer and other diseases. In particular, we summarise the roles of MOBs as core components of the Hippo tissue growth and regeneration pathway.

Functional testing of thousands of osteoarthritis-associated variants for regulatory activity

To date, genome-wide association studies have implicated at least 35 loci in osteoarthritis but, due to linkage disequilibrium, the specific variants underlying these associations and the mechanisms by which they contribute to disease risk have yet to be pinpointed. Here, we functionally test 1,605 single nucleotide variants associated with osteoarthritis for regulatory activity using a massively parallel reporter assay. We identify six single nucleotide polymorphisms (SNPs) with differential regulatory activity between the major and minor alleles. We show that the most significant SNP, rs4730222, exhibits differential nuclear protein binding in electrophoretic mobility shift assays and drives increased expression of an alternative isoform of HBP1 in a heterozygote chondrosarcoma cell line, in a CRISPR-edited osteosarcoma cell line, and in chondrocytes derived from osteoarthritis patients. This study provides a framework for prioritization of GWAS variants and highlights a role of HBP1 and Wnt signaling in osteoarthritis pathogenesis.

Role of SREBP2 gene polymorphism on knee osteoarthritis in the South Indian Hyderabad Population: A hospital based study with G595C variant

Introduction

Osteoarthritis (OA) is a multifactorial disease with genetic factors playing a crucial role, and it has been associated with a family history of obesity. G595C polymorphism in the sterol regulatory element-binding protein 2 (SREBP2) gene has demonstrated an association with knee osteoarthritis (KOA) patients. However, this polymorphism has been never explored in an Indian population. Hence, the current study aimed to examine whether G595C (rs2228314) polymorphism in SREBP2 gene was associated with KOA susceptibility in the South Indian Hyderabad population.

Methods

G595C polymorphism was genotyped with 200 KOA cases and 200 healthy controls using polymerase chain reaction-restriction fragment length polymorphism analysis.

Results

A significant association was observed between age, body mass index (BMI), and family histories in KOA cases and controls (p < 0.05). The current allele (C vs G; OR-2.8 [95%CI = 2.1-3.7]; p < 0.0001) and genotype analysis confirms the significant association with (GC + CC vs GG; OR-3.5 [95%CI = 2.3-5.3]; p < 0.0001 & GC vs GG + CC; OR-1.7 [95%CI = 1.0-2.9]; p = 0.02) KOA vs. control subjects. On stratification analysis, genotype CC and C allele were associated with KOA. Gender association failed to demonstrate positive genotype frequencies (p > 0.05). Multifactor-dimensionality reduction (MDR) analysis showed a positive association with BMI and G595C genotypes (p < 0.05); 51% of the homozygous variant CC genotypes were present in obesity subjects.

Conclusion

In conclusion, our findings suggest that G595C polymorphism in SREBP2 gene is associated with KOA in the South Indian Hyderabad population and presents scope for further investigation of the gene's function in KOA.

eQTLs Weighted Genetic Correlation Analysis Detected Brain Region Differences in Genetic Correlations for Complex Psychiatric Disorders

BACKGROUND

Psychiatric disorders are usually caused by the dysfunction of various brain regions. Incorporating the genetic information of brain regions into correlation analysis can provide novel clues for pathogenetic and therapeutic studies of psychiatric disorders.

METHODS

The latest genome-wide association study (GWAS) summary data of schizophrenia (SCZ), bipolar disorder (BIP), autism spectrum disorder (AUT), major depression disorder (MDD), and attention-deficit/hyperactivity disorder (ADHD) were obtained from the Psychiatric GWAS Consortium (PGC). The expression quantitative trait loci (eQTLs) datasets of 10 brain regions were driven from the genotype-tissue expression (GTEx) database. The PGC GWAS summaries were first weighted by the GTEx eQTLs summaries for each brain region. Linkage disequilibrium score regression was applied to the weighted GWAS summary data to detect genetic correlation for each pair of 5 disorders.

RESULTS

Without considering brain region difference, significant genetic correlations were observed for BIP vs SCZ (P = 1.68 × 10-63), MDD vs SCZ (P = 5.08 × 10-45), ADHD vs MDD (P = 1.93 × 10-44), BIP vs MDD (P = 6.39 × 10-9), AUT vs SCZ (P = .0002), and ADHD vs SCZ (P = .0002). Utilizing brain region related eQTLs weighted LD score regression, different strengths of genetic correlations were observed within various brain regions for BIP vs SCZ, MDD vs SCZ, ADHD vs MDD, and SCZ vs ADHD. For example, the most significant genetic correlations were observed at anterior cingulate cortex (P = 1.13 × 10-34) for BIP vs SCZ.

CONCLUSIONS

This study provides new clues for elucidating the mechanism of genetic correlations among various psychiatric disorders.

Cardiometabolic Polygenic Risk Scores and Osteoarthritis Outcomes: A Mendelian Randomization Study Using Data From the Malmö Diet and Cancer Study and the UK Biobank

Objective

To investigate the causal role of cardiometabolic risk factors in osteoarthritis (OA) using associated genetic variants.

Methods

We studied 27,691 adults from the Malmö Diet and Cancer Study (MDCS) and replicated novel findings among 376,435 adults from the UK Biobank. Trait‐specific polygenic risk scores for low‐density lipoprotein (LDL) and high‐density lipoprotein (HDL) cholesterol levels, triglyceride levels, body mass index (BMI), fasting plasma glucose (FPG) levels, and systolic blood pressure (BP) were used to test the associations of genetically predicted elevations in each trait with incident OA diagnosis (n = 3,559), OA joint replacement (n = 2,780), or both (total OA; n = 4,226) in Mendelian randomization (MR) analyses in the MDCS, and with self‐reported and/or hospital‐diagnosed OA (n = 65,213) in the UK Biobank. Multivariable MR, MR‐Egger, and weighted median MR were used to adjust for potential pleiotropic biases.

Results

In the MDCS, genetically predicted elevation in LDL cholesterol level was associated with a lower risk of OA diagnosis (odds ratio [OR] 0.83 [95% confidence interval (95% CI) 0.73–0.95] per 1SD increase) and total OA (OR 0.87 [95% CI 0.78–0.98]), which was supported by multivariable MR for OA diagnosis (OR 0.84 [95% CI 0.75–0.95]) and total OA (0.87 [95% CI 0.78–0.97]), and by conventional 2‐sample MR for OA diagnosis (OR 0.86 [95% CI 0.75–0.98]). MR‐Egger indicated no pleiotropic bias. Genetically predicted elevation in BMI was associated with an increased risk of OA diagnosis (OR 1.65 [95% CI 1.14–2.41]), while MR‐Egger indicated pleiotropic bias and a larger association with OA diagnosis (OR 3.25 [1.26–8.39]), OA joint replacement (OR 3.81 [95% CI 1.39–10.4]), and total OA (OR 3.41 [95% CI 1.43–8.15]). No associations were observed between genetically predicted HDL cholesterol level, triglyceride level, FPG level, or systolic BP and OA outcomes. The associations with LDL cholesterol levels were replicated in the UK Biobank (OR 0.95 [95% CI 0.93–0.98]).

Conclusion

Our MR study provides evidence of a causal role of lower LDL cholesterol level and higher BMI in OA.

Nature vs nurture in knee osteoarthritis - the importance of age, sex and body mass index

OBJECTIVE

(1) To estimate the life-time genetic contribution for knee osteoarthritis (OA) surgery and (2) to explore any differences in the genetic contribution across age, sex and body mass index (BMI).

METHODS

We studied the sex-specific genetic contribution to knee OA surgery in a prospective cohort study of 62,490 twins aged 35 years or older with a follow-up period of up to 47 years (10,092 identical and 21,153 non-identical twin pairs, 54% women). To study interactions with age, we graphed the heritabilities over the lifespan for men and women. We also studied the sex-specific heritability across strata of the median BMI to explore any interactions with BMI.

RESULTS

The overall heritability of knee OA surgery was 0.53 (95% confidence intervals [CI] = 0.31-0.75), with higher heritability among women (H² = 0.80 (95% CI = 0.73-0.87)) than men (H² = 0.39 (95% CI = 0.10-0.69)). For men, the heritability started to rise after age 68. The genetic contribution was particularly low in men above median BMI (H²_{≥23.7 kg/m2} = 0.08, 95% CI = -0.32-0.48). For women, the heritability was consistently high from age 50 to death, independently of BMI (H²_{≥22.5 kg/m2} = 0.77, 95% CI = 0.66-0.87).

CONCLUSION

There is a higher and more consistent genetic contribution for knee OA surgery in women than men. In men the genetic contribution was relatively low and varied with age and BMI.

IL-10 Could Play a Role in the Interrelation between Diabetes Mellitus and Osteoarthritis

The association between osteoarthritis (OA), obesity and metabolic syndrome suggests an interrelation between OA and diabetes mellitus (DM). Little is known about the role of anti-inflammatory cytokine interleukin (IL)-10 in the interrelation between OA and DM. Hence, the effects of IL-10 under hyperglycemia (HG) and hyperinsulinemia (HI) in human articular chondrocytes (hAC) and chondrosarcoma cell line Okayama University Medical School (OUMS)-27 were examined. HAC and OUMS-27, cultured in normoglycemic (NG) and HG conditions were stimulated with insulin and/or IL-10. Cell survival, metabolic activity, proliferation and extracellular matrix (ECM) synthesis were immunocytochemically examined. No significant differences in vitality of hAC neither in pure NG (NG_w/o) nor HG (HG_w/o) conditions were found. Applying HI and/or IL-10 in both conditions reduced significantly the vitality of hAC but not of OUMS-27. HG impaired significantly hAC metabolism. When combined with HI + IL-10 or IL-10 alone it decreased also significantly hAC proliferation compared to NG_w/o. In OUMS-27 it induced only a trend of impaired proliferation compared to NG_w/o. hAC but not OUMS-27 reduced significantly their collagen type (col) I, SOX9 and proteoglycan (PG) synthesis in HG combined with HI +/− IL-10 compared to NG_w/o. IL-10 could not moderate HI and HG effects. In contrast to hAC OUMS-27 showed limited sensitivity as DM model.

The Korea Biobank Array: Design and Identification of Coding Variants Associated with Blood Biochemical Traits

We introduce the design and implementation of a new array, the Korea Biobank Array (referred to as KoreanChip), optimized for the Korean population and demonstrate findings from GWAS of blood biochemical traits. KoreanChip comprised >833,000 markers including >247,000 rare-frequency or functional variants estimated from >2,500 sequencing data in Koreans. Of the 833 K markers, 208 K functional markers were directly genotyped. Particularly, >89 K markers were presented in East Asians. KoreanChip achieved higher imputation performance owing to the excellent genomic coverage of 95.38% for common and 73.65% for low-frequency variants. From GWAS (Genome-wide association study) using 6,949 individuals, 28 associations were successfully recapitulated. Moreover, 9 missense variants were newly identified, of which we identified new associations between a common population-specific missense variant, rs671 (p.Glu457Lys) of ALDH2, and two traits including aspartate aminotransferase (P = 5.20 × 10⁻¹³) and alanine aminotransferase (P = 4.98 × 10⁻⁸). Furthermore, two novel missense variants of GPT with rare frequency in East Asians but extreme rarity in other populations were associated with alanine aminotransferase (rs200088103; p.Arg133Trp, P = 2.02 × 10⁻⁹ and rs748547625; p.Arg143Cys, P = 1.41 × 10⁻⁶). These variants were successfully replicated in 6,000 individuals (P = 5.30 × 10⁻⁸ and P = 1.24 × 10⁻⁶). GWAS results suggest the promising utility of KoreanChip with a substantial number of damaging variants to identify new population-specific disease-associated rare/functional variants.

Identification of Novel Loci Associated With Hip Shape: A Meta-Analysis of Genomewide Association Studies

We aimed to report the first genomewide association study (GWAS) meta-analysis of dual-energy X-ray absorptiometry (DXA)-derived hip shape, which is thought to be related to the risk of both hip osteoarthritis and hip fracture. Ten hip shape modes (HSMs) were derived by statistical shape modeling using SHAPE software, from hip DXA scans in the Avon Longitudinal Study of Parents and Children (ALSPAC; adult females), TwinsUK (mixed sex), Framingham Osteoporosis Study (FOS; mixed), Osteoporotic Fractures in Men study (MrOS), and Study of Osteoporotic Fractures (SOF; females) (total N = 15,934). Associations were adjusted for age, sex, and ancestry. Five genomewide significant (p < 5 × 10-9 , adjusted for 10 independent outcomes) single-nucleotide polymorphisms (SNPs) were associated with HSM1, and three SNPs with HSM2. One SNP, in high linkage disequilibrium with rs2158915 associated with HSM1, was associated with HSM5 at genomewide significance. In a look-up of previous GWASs, three of the identified SNPs were associated with hip osteoarthritis, one with hip fracture, and five with height. Seven SNPs were within 200 kb of genes involved in endochondral bone formation, namely SOX9, PTHrP, RUNX1, NKX3-2, FGFR4, DICER1, and HHIP. The SNP adjacent to DICER1 also showed osteoblast cis-regulatory activity of GSC, in which mutations have previously been reported to cause hip dysplasia. For three of the lead SNPs, SNPs in high LD (r2  > 0.5) were identified, which intersected with open chromatin sites as detected by ATAC-seq performed on embryonic mouse proximal femora. In conclusion, we identified eight SNPs independently associated with hip shape, most of which were associated with height and/or mapped close to endochondral bone formation genes, consistent with a contribution of processes involved in limb growth to hip shape and pathological sequelae. These findings raise the possibility that genetic studies of hip shape might help in understanding potential pathways involved in hip osteoarthritis and hip fracture. © 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals, Inc.

Role of Epigenomics in Bone and Cartilage Disease

Phenotypic variation in skeletal traits and diseases is the product of genetic and environmental factors. Epigenetic mechanisms include information-containing factors, other than DNA sequence, that cause stable changes in gene expression and are maintained during cell divisions. They represent a link between environmental influences, genome features, and the resulting phenotype. The main epigenetic factors are DNA methylation, posttranslational changes of histones, and higher-order chromatin structure. Sometimes non-coding RNAs, such as microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), are also included in the broad term of epigenetic factors. There is rapidly expanding experimental evidence for a role of epigenetic factors in the differentiation of bone cells and the pathogenesis of skeletal disorders, such as osteoporosis and osteoarthritis. However, different from genetic factors, epigenetic signatures are cell- and tissue-specific and can change with time. Thus, elucidating their role has particular difficulties, especially in human studies. Nevertheless, epigenomewide association studies are beginning to disclose some disease-specific patterns that help to understand skeletal cell biology and may lead to development of new epigenetic-based biomarkers, as well as new drug targets useful for treating diffuse and localized disorders. Here we provide an overview and update of recent advances on the role of epigenomics in bone and cartilage diseases. © 2019 American Society for Bone and Mineral Research.

EGFR Signaling: Friend or Foe for Cartilage?

Recent studies using genetically modified mice, pharmacological approaches, and human samples have highlighted an important role for the epidermal growth factor receptor (EGFR), selected ligands, and downstream components in endochondral bone formation and joint homeostasis. Although most data demonstrate an important function of this pathway in endochondral ossification and articular cartilage growth, conflicting results on its role in osteoarthritis have been reported. In some contexts, inactivation of EGFR signaling has been shown to protect joints from surgically induced osteoarthritis, whereas in others, similar manipulations worsened joint pathology. The current review summarizes recent studies of cartilage EGFR signaling in long bone development and diseases, provides potential explanations for the reported discrepancies, and suggests directions for future work to clarify the potential of this pathway as target for osteoarthritis treatment. © 2019 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

The 2018 Otto Aufranc Award: How Does Genome-wide Variation Affect Osteolysis Risk After THA?

BACKGROUND

Periprosthetic osteolysis resulting in aseptic loosening is a leading cause of THA revision. Individuals vary in their susceptibility to osteolysis and heritable factors may contribute to this variation. However, the overall contribution that such variation makes to osteolysis risk is unknown.

QUESTIONS/PURPOSES

We conducted two genome-wide association studies to (1) identify genetic risk loci associated with susceptibility to osteolysis; and (2) identify genetic risk loci associated with time to prosthesis revision for osteolysis.

METHODS

The Norway cohort comprised 2624 patients after THA recruited from the Norwegian Arthroplasty Registry, of whom 779 had undergone revision surgery for osteolysis. The UK cohort included 890 patients previously recruited from hospitals in the north of England, 317 who either had radiographic evidence of and/or had undergone revision surgery for osteolysis. All participants had received a fully cemented or hybrid THA using a small-diameter metal or ceramic-on-conventional polyethylene bearing. Osteolysis susceptibility case-control analyses and quantitative trait analyses for time to prosthesis revision (a proxy measure of the speed of osteolysis onset) in those patients with osteolysis were undertaken in each cohort separately after genome-wide genotyping. Finally, a meta-analysis of the two independent cohort association analysis results was undertaken.

RESULTS

Genome-wide association analysis identified four independent suggestive genetic signals for osteolysis case-control status in the Norwegian cohort and 11 in the UK cohort (p ≤ 5 x 10). After meta-analysis, five independent genetic signals showed a suggestive association with osteolysis case-control status at p ≤ 5 x 10 with the strongest comprising 18 correlated variants on chromosome 7 (lead signal rs850092, p = 1.13 x 10). Genome-wide quantitative trait analysis in cases only showed a total of five and nine independent genetic signals for time to revision at p ≤ 5 x 10, respectively. After meta-analysis, 11 independent genetic signals showed suggestive evidence of an association with time to revision at p ≤ 5 x 10 with the largest association block comprising 174 correlated variants in chromosome 15 (lead signal rs10507055, p = 1.40 x 10).

CONCLUSIONS

We explored the heritable biology of osteolysis at the whole genome level and identify several genetic loci that associate with susceptibility to osteolysis or with premature revision surgery. However, further studies are required to determine a causal association between the identified signals and osteolysis and their functional role in the disease.

CLINICAL RELEVANCE

The identification of novel genetic risk loci for osteolysis enables new investigative avenues for clinical biomarker discovery and therapeutic intervention in this disease.

Evidence for Genetic Contribution to Variation in Posttraumatic Osteoarthritis in Mice

OBJECTIVE

Recombinant inbred mouse strains generated from an LG/J and SM/J intercross offer a unique resource to study complex genetic traits such as osteoarthritis (OA). We undertook this study to determine the susceptibility of 14 strains to various phenotypes characteristic of posttraumatic OA. We hypothesized that phenotypic variability is associated with genetic variability.

METHODS

Ten-week-old male mice underwent surgical destabilization of the medial meniscus (DMM) to induce posttraumatic OA. Mice were killed 8 weeks after surgery, and knee joints were processed for histology to score cartilage degeneration and synovitis. Micro-computed tomography was used to analyze trabecular bone parameters including subchondral bone plate thickness and synovial ectopic calcifications. Gene expression in the knees was assessed using a QuantiGene Plex assay.

RESULTS

Broad-sense heritability ranged from 0.18 to 0.58, which suggested that the responses to surgery were moderately heritable. The LGXSM-33, LGXSM-5, LGXSM-46, and SM/J strains were highly susceptible to OA, while the LGXSM-131b, LGXSM-163, LGXSM-35, LGXSM-128a, LGXSM-6, and LG/J strains were relatively OA resistant. This study was the first to accomplish measurement of genetic correlations of phenotypes that are characteristic of posttraumatic OA. Cartilage degeneration was significantly positively associated with synovitis (r = 0.83-0.92), and subchondral bone plate thickness was negatively correlated with ectopic calcifications (r = -0.59). Moreover, we showed that 40 of the 78 genes tested were significantly correlated with various OA phenotypes. However, unlike the OA phenotypes, there was no evidence for genetic variation in differences in gene expression levels between DMM-operated and sham-operated knees.

CONCLUSION

For these mouse strains, various characteristics of posttraumatic OA varied with genetic composition, which demonstrated a genetic basis for susceptibility to posttraumatic OA. The heritability of posttraumatic OA was established. Phenotypes exhibited various degrees of correlations; cartilage degeneration was positively correlated with synovitis, but not with the formation of ectopic calcifications. Further investigation of the genome regions that contain genes implicated in OA, as well as further investigation of gene expression data, will be useful for studying mechanisms of OA and identifying therapeutic targets.

Identification of new therapeutic targets for osteoarthritis through genome-wide analyses of UK Biobank data

Osteoarthritis is the most common musculoskeletal disease and the leading cause of disability globally. Here, we perform a genome-wide association study for osteoarthritis (77,052 cases and 378,169 controls), analysing 4 phenotypes: knee osteoarthritis, hip osteoarthritis, knee and/or hip osteoarthritis, and any osteoarthritis. We discover 64 signals, 52 of them novel, more than doubling the number of established disease loci. Six signals fine map to a single variant. We identify putative effector genes by integrating eQTL colocalization, fine-mapping, human rare disease, animal model, and osteoarthritis tissue expression data. We find enrichment for genes underlying monogenic forms of bone development diseases, and for the collagen formation and extracellular matrix organisation biological pathways. Ten of the likely effector genes, including TGFB1, FGF18, CTSK and IL11 have therapeutics approved or in clinical trials, with mechanisms of action supportive of evaluation for efficacy in osteoarthritis.

Targeting of chondrocyte plasticity via connexin43 modulation attenuates cellular senescence and fosters a pro-regenerative environment in osteoarthritis

Osteoarthritis (OA), a chronic disease characterized by articular cartilage degeneration, is a leading cause of disability and pain worldwide. In OA, chondrocytes in cartilage undergo phenotypic changes and senescence, restricting cartilage regeneration and favouring disease progression. Similar to other wound-healing disorders, chondrocytes from OA patients show a chronic increase in the gap junction channel protein connexin43 (Cx43), which regulates signal transduction through the exchange of elements or recruitment/release of signalling factors. Although immature or stem-like cells are present in cartilage from OA patients, their origin and role in disease progression are unknown. In this study, we found that Cx43 acts as a positive regulator of chondrocyte-mesenchymal transition. Overactive Cx43 largely maintains the immature phenotype by increasing nuclear translocation of Twist-1 and tissue remodelling and proinflammatory agents, such as MMPs and IL-1β, which in turn cause cellular senescence through upregulation of p53, p16^INK4a and NF-κB, contributing to the senescence-associated secretory phenotype (SASP). Downregulation of either Cx43 by CRISPR/Cas9 or Cx43-mediated gap junctional intercellular communication (GJIC) by carbenoxolone treatment triggered rediferentiation of osteoarthritic chondrocytes into a more differentiated state, associated with decreased synthesis of MMPs and proinflammatory factors, and reduced senescence. We have identified causal Cx43-sensitive circuit in chondrocytes that regulates dedifferentiation, redifferentiation and senescence. We propose that chondrocytes undergo chondrocyte-mesenchymal transition where increased Cx43-mediated GJIC during OA facilitates Twist-1 nuclear translocation as a novel mechanism involved in OA progression. These findings support the use of Cx43 as an appropriate therapeutic target to halt OA progression and to promote cartilage regeneration.

Update on the epidemiology, risk factors and disease outcomes of osteoarthritis

Osteoarthritis (OA) is the most frequent form of arthritis and a leading cause of pain and disability worldwide. OA can affect any synovial joint, although the hip, knee, hand, foot and spine are the most commonly affected sites. Knowledge about the occurrence and risk factors for OA is important to define the clinical and public health burden of the disease to understand mechanisms of disease occurrence and may also help to inform the development of population-wide prevention strategies. In this article, we review the occurrence and risk factors for OA and also consider patient-reported outcome measures that have been used for the assessment of the disease.

Meta-analysis of Icelandic and UK data sets identifies missense variants in SMO, IL11, COL11A1 and 13 more new loci associated with osteoarthritis

Osteoarthritis has a highly negative impact on quality of life because of the associated pain and loss of joint function. Here we describe the largest meta-analysis so far of osteoarthritis of the hip and the knee in samples from Iceland and the UK Biobank (including 17,151 hip osteoarthritis patients, 23,877 knee osteoarthritis patients, and more than 562,000 controls). We found 23 independent associations at 22 loci in the additive meta-analyses, of which 16 of the loci were novel: 12 for hip and 4 for knee osteoarthritis. Two associations are between rare or low-frequency missense variants and hip osteoarthritis, affecting the genes SMO (rs143083812, frequency 0.11%, odds ratio (OR) = 2.8, P = 7.9 × 10^-12, p.Arg173Cys) and IL11 (rs4252548, frequency 2.08%, OR = 1.30, P = 2.1 × 10^-11, p.Arg112His). A common missense variant in the COL11A1 gene also associates with hip osteoarthritis (rs3753841, frequency 61%, P = 5.2 × 10^-10, OR = 1.08, p.Pro1284Leu). In addition, using a recessive model, we confirm an association between hip osteoarthritis and a variant of CHADL¹ (rs117018441, P = 1.8 × 10^-25, OR = 5.9). Furthermore, we observe a complex relationship between height and risk of osteoarthritis.

Chromatin accessibility landscape of articular knee cartilage reveals aberrant enhancer regulation in osteoarthritis

Osteoarthritis (OA) is a common joint disorder with increasing impact in an aging society. While genetic and transcriptomic analyses have revealed some genes and non-coding loci associated to OA, the pathogenesis remains incompletely understood. Chromatin profiling, which provides insight into gene regulation, has not been reported in OA mainly due to technical difficulties. Here, we employed Assay for Transposase-Accessible Chromatin with high throughput sequencing (ATAC-seq) to map the accessible chromatin landscape in articular knee cartilage of OA patients. We identified 109,215 accessible chromatin regions for cartilages, of which 71% were annotated as enhancers. By overlaying them with genetic and DNA methylation data, we have determined potential OA-relevant enhancers and their putative target genes. Furthermore, through integration with RNA-seq data, we characterized genes that are altered both at epigenomic and transcriptomic levels in OA. These genes are enriched in pathways regulating ossification and mesenchymal stem cell (MSC) differentiation. Consistently, the differentially accessible regions in OA are enriched for MSC-specific enhancers and motifs of transcription factor families involved in osteoblast differentiation. In conclusion, we demonstrate how direct chromatin profiling of clinical tissues can provide comprehensive epigenetic information for a disease and suggest candidate genes and enhancers of translational potential.

Biobank integration of large-scale clinical and histopathology melanoma studies within the European Cancer Moonshot Lund Center

We present the Cancer Moonshot clinical project located at the European center in Lund. Here, tissue and blood samples have been collected and stored in a large-scale biobank. Multiple clinical centers around the world are participating and tissue and blood samples are sent to the European Cancer Moonshot Lund Center that acts as the clinical hub. Our center has been developed to generate and build large-scale biostorage archives of patient melanoma samples, which is then combined with a histopathological capability to characterize the patient tumours. Such a large-scale clinical sample processing initiative has begun with the aim of creating high-end histopathology indexing with database computational power and including proteogenomic analysis. The biobank at Lund has become an important resource in clinical research worldwide. Following suite, several national health programs are being initiated with the aim of also building large-scale biobank storages with a wealth of high-quality patient samples. In our Cancer Moonshot R&D activities, samples in the biobanks and the data derived from these samples are being used to build an understanding of disease presentation and using this information to move towards 'Big Data' proteogenomic and mass spectrometry imaging studies. Additionally, we report here a sample processing workflow that has been adapted to a fully-automated biobank processing strategy for large-scale studies.

TLR10 and NFKBIA contributed to the risk of hip osteoarthritis: systematic evaluation based on Han Chinese population

Multiple lines of evidence have confirmed the importance of genetic factors for hip osteoarthritis (HOA). Our study aimed to investigate the associations of TLR10 and NFKBIA with respect to the HOA risk in Han Chinese individuals. A total of 1,043 HOA patients and 2,664 controls were recruited. Then, 23 tag single-nucleotide polymorphisms (SNPs) in the TLR10 and NFKBIA genes were selected for genotyping. Genetic association analyses were conducted in both single-marker and haplotype-based ways. Gene by gene, two-way interactions were analysed using a case-only method. Multiple bioinformatics tools were utilised to examine the potential functional significance of the SNPs. Two significant SNPs, rs11096957 (OR = 1.26, P = 1.35 × 10⁻⁵) and rs2273650 (OR = 1.2, P = 1.57 × 10⁻³), were significantly associated with HOA risk. Rs11096957 was also associated with the severity of the HOA. Bioinformatics analysis indicated that the allele T of rs2273650 would create new miRNA/SNP target duplexes, which suggests that rs2273650 could alter the NFKBIA expression by affecting the miRNA/SNP target duplexes. Our study identified significant association signals from NFKBIA with HOA for the first time, and it also confirmed the contribution of TLR10 to the HOA risk. These findings would provide clues for identifying individuals at high risk of HOA.

Medical relevance of protein-truncating variants across 337,205 individuals in the UK Biobank study

Protein-truncating variants can have profound effects on gene function and are critical for clinical genome interpretation and generating therapeutic hypotheses, but their relevance to medical phenotypes has not been systematically assessed. Here, we characterize the effect of 18,228 protein-truncating variants across 135 phenotypes from the UK Biobank and find 27 associations between medical phenotypes and protein-truncating variants in genes outside the major histocompatibility complex. We perform phenome-wide analyses and directly measure the effect in homozygous carriers, commonly referred to as “human knockouts,” across medical phenotypes for genes implicated as being protective against disease or associated with at least one phenotype in our study. We find several genes with strong pleiotropic or non-additive effects. Our results illustrate the importance of protein-truncating variants in a variety of diseases.

Protein-truncating variants (PTVs) are predicted to significantly affect a gene’s function and, thus, human traits. Here, DeBoever et al. systematically analyze PTVs in more than 300,000 individuals across 135 phenotypes and identify 27 associations between PTVs and medical conditions.