Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations.

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation.

Evaluation of shared genetic aetiology between osteoarthritis and bone mineral density identifies SMAD3 as a novel osteoarthritis risk locus.

Osteoarthritis (OA) is a common complex disease with high public health burden and no curative therapy. High bone mineral density (BMD) is associated with an increased risk of developing OA, suggesting a shared underlying biology. Here, we performed the first systematic overlap analysis of OA and BMD on a genome wide scale. We used summary statistics from the GEFOS consortium for lumbar spine (n = 31,800) and femoral neck (n = 32,961) BMD, and from the arcOGEN consortium for three OA phenotypes (hip, ncases=3,498; knee, ncases=3,266; hip and/or knee, ncases=7,410; ncontrols=11,009). Performing LD score regression we found a significant genetic correlation between the combined OA phenotype (hip and/or knee) and lumbar spine BMD (rg=0.18, P = 2.23 × 10-2), which may be driven by the presence of spinal osteophytes. We identified 143 variants with evidence for cross-phenotype association which we took forward for replication in independent large-scale OA datasets, and subsequent meta-analysis with arcOGEN for a total sample size of up to 23,425 cases and 236,814 controls. We found robustly replicating evidence for association with OA at rs12901071 (OR 1.08 95% CI 1.05-1.11, Pmeta=3.12 × 10-10), an intronic variant in the SMAD3 gene, which is known to play a role in bone remodeling and cartilage maintenance. We were able to confirm expression of SMAD3 in intact and degraded cartilage of the knee and hip. Our findings provide the first systematic evaluation of pleiotropy between OA and BMD, highlight genes with biological relevance to both traits, and establish a robust new OA genetic risk locus at SMAD3.

A novel variant in GLIS3 is associated with osteoarthritis.

OBJECTIVES: Osteoarthritis (OA) is a complex disease, but its genetic aetiology remains poorly characterised. To identify novel susceptibility loci for OA, we carried out a genome-wide association study (GWAS) in individuals from the largest UK-based OA collections to date. METHODS: We carried out a discovery GWAS in 5414 OA individuals with knee and/or hip total joint replacement (TJR) and 9939 population-based controls. We followed-up prioritised variants in OA subjects from the interim release of the UK Biobank resource (up to 12 658 cases and 50 898 controls) and our lead finding in operated OA subjects from the full release of UK Biobank (17 894 cases and 89 470 controls). We investigated its functional implications in methylation, gene expression and proteomics data in primary chondrocytes from 12 pairs of intact and degraded cartilage samples from patients undergoing TJR. RESULTS: We detect a genome-wide significant association at rs10116772 with TJR (P=3.7×10-8; for allele A: OR (95% CI) 0.97 (0.96 to 0.98)), an intronic variant in GLIS3, which is expressed in cartilage. Variants in strong correlation with rs10116772 have been associated with elevated plasma glucose levels and diabetes. CONCLUSIONS: We identify a novel susceptibility locus for OA that has been previously implicated in diabetes and glycaemic traits.

Widespread epigenomic, transcriptomic and proteomic differences between hip osteophytic and articular chondrocytes in osteoarthritis.

Objectives: To identify molecular differences between chondrocytes from osteophytic and articular cartilage tissue from OA patients. Methods: We investigated genes and pathways by combining genome-wide DNA methylation, RNA sequencing and quantitative proteomics in isolated primary chondrocytes from the cartilaginous layer of osteophytes and matched areas of low- and high-grade articular cartilage across nine patients with OA undergoing hip replacement surgery. Results: Chondrocytes from osteophytic cartilage showed widespread differences to low-grade articular cartilage chondrocytes. These differences were similar to, but more pronounced than, differences between chondrocytes from osteophytic and high-grade articular cartilage, and more pronounced than differences between high- and low-grade articular cartilage. We identified 56 genes with significant differences between osteophytic chondrocytes and low-grade articular cartilage chondrocytes on all three omics levels. Several of these genes have known roles in OA, including ALDH1A2 and cartilage oligomeric matrix protein, which have functional genetic variants associated with OA from genome-wide association studies. An integrative gene ontology enrichment analysis showed that differences between osteophytic and low-grade articular cartilage chondrocytes are associated with extracellular matrix organization, skeletal system development, platelet aggregation and regulation of ERK1 and ERK2 cascade. Conclusion: We present a first comprehensive view of the molecular landscape of chondrocytes from osteophytic cartilage as compared with articular cartilage chondrocytes from the same joints in OA. We found robust changes at genes relevant to chondrocyte function, providing insight into biological processes involved in osteophyte development and thus OA progression.

Combination therapy as a potential risk factor for the development of type 2 diabetes in patients with schizophrenia: the GOMAP study.

BACKGROUND: Schizophrenia (SCZ) is associated with increased risk of type 2 diabetes (T2D). The potential diabetogenic effect of concomitant application of psychotropic treatment classes in patients with SCZ has not yet been evaluated. The overarching goal of the Genetic Overlap between Metabolic and Psychiatric disease (GOMAP) study is to assess the effect of pharmacological, anthropometric, lifestyle and clinical measurements, helping elucidate the mechanisms underlying the aetiology of T2D. METHODS: The GOMAP case-control study (Genetic Overlap between Metabolic and Psychiatric disease) includes hospitalized patients with SCZ, some of whom have T2D. We enrolled 1653 patients with SCZ; 611 with T2D and 1042 patients without T2D. This is the first study of SCZ and T2D comorbidity at this scale in the Greek population. We retrieved detailed information on first- and second-generation antipsychotics (FGA, SGA), antidepressants and mood stabilizers, applied as monotherapy, 2-drug combination, or as 3- or more drug combination. We assessed the effects of psychotropic medication, body mass index, duration of schizophrenia, number of hospitalizations and physical activity on risk of T2D. Using logistic regression, we calculated crude and adjusted odds ratios (OR) to identify associations between demographic factors and the psychiatric medications. RESULTS: Patients with SCZ on a combination of at least three different classes of psychiatric drugs had a higher risk of T2D [OR 1.81 (95% CI 1.22-2.69); p = 0.003] compared to FGA alone therapy, after adjustment for age, BMI, sex, duration of SCZ and number of hospitalizations. We did not find evidence for an association of SGA use or the combination of drugs belonging to two different classes of psychiatric medications with increased risk of T2D [1.27 (0.84-1.93), p = 0.259 and 0.98 (0.71-1.35), p = 0.885, respectively] compared to FGA use. CONCLUSIONS: We find an increased risk of T2D in patients with SCZ who take a combination of at least three different psychotropic medication classes compared to patients whose medication consists only of one or two classes of drugs.

Radiographic endophenotyping in hip osteoarthritis improves the precision of genetic association analysis.

OBJECTIVE: Osteoarthritis (OA) has a strong genetic component but the success of previous genome-wide association studies (GWAS) has been restricted due to insufficient sample sizes and phenotype heterogeneity. Our aim was to examine the effect of clinically relevant endophenotyping according to site of maximal joint space narrowing (maxJSN) and bone remodelling response on GWAS signal detection in hip OA. METHODS: A stratified GWAS meta-analysis was conducted in 2118 radiographically defined hip OA cases and 6500 population-based controls. Signals were followed up by analysing differential expression of proximal genes for bone remodelling endophenotypes in 33 pairs of macroscopically intact and OA-affected cartilage. RESULTS: We report suggestive evidence (p<5×10-6) of association at 6 variants with OA endophenotypes that would have been missed by using presence of hip OA as the disease end point. For example, in the analysis of hip OA cases with superior maxJSN versus cases with non-superior maxJSN we detected association with a variant in the LRCH1 gene (rs754106, p=1.49×10-7, OR (95% CIs) 0.70 (0.61 to 0.80)). In the comparison of hypertrophic with non-hypertrophic OA the most significant variant was located between STT3B and GADL1 (rs6766414, p=3.13×10-6, OR (95% CIs) 1.45 (1.24 to 1.69)). Both of these associations were fully attenuated in non-stratified analyses of all hip OA cases versus population controls (p>0.05). STT3B was significantly upregulated in OA-affected versus intact cartilage, particularly in the analysis of hypertrophic and normotrophic compared with atrophic bone remodelling pattern (p=4.2×10-4). CONCLUSIONS: Our findings demonstrate that stratification of OA cases into more homogeneous endophenotypes can identify genes of potential functional importance otherwise obscured by disease heterogeneity.

A Dietary Pattern with High Sugar Content Is Associated with Cardiometabolic Risk Factors in the Pomak Population.

The present study describes the geographically isolated Pomak population and its particular dietary patterns in relationship to cardiovascular risk factors. We collected a population-based cohort in a cross-sectional study, with detailed anthropometric, biochemical, clinical, and lifestyle parameter information. Dietary patterns were derived through principal component analysis based on a validated food-frequency questionnaire, administered to 1702 adult inhabitants of the Pomak villages on the Rhodope mountain range in Greece. A total of 69.9% of the participants were female with a population mean age of 44.9 years; 67% of the population were overweight or obese with a significantly different prevalence for obesity between men and women (17.5% vs. 37.5%, respectively, p < 0.001). Smoking was more prevalent in men (45.8% vs. 2.2%, p < 0.001), as 97.3% of women had never smoked. Four dietary patterns emerged as characteristic of the population, and were termed "high in sugars", "quick choices", "balanced", and "homemade". Higher adherence to the "high in sugars" dietary pattern was associated with increased glucose levels (p < 0.001) and increased risk of hypertension (OR (95% CI) 2.61 (1.55, 4.39), p < 0.001) and nominally associated with high blood glucose levels (OR (95% CI) 1.85 (1.11, 3.08), p = 0.018), compared to lower adherence. Overall, we characterize the dietary patterns of the Pomak population and describe associations with cardiovascular risk factors.

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 performed a genome-wide association study for osteoarthritis (77,052 cases and 378,169 controls), analyzing four phenotypes: knee osteoarthritis, hip osteoarthritis, knee and/or hip osteoarthritis, and any osteoarthritis. We discovered 64 signals, 52 of them novel, more than doubling the number of established disease loci. Six signals fine-mapped to a single variant. We identified putative effector genes by integrating expression quantitative trait loci (eQTL) colocalization, fine-mapping, and human rare-disease, animal-model, and osteoarthritis tissue expression data. We found enrichment for genes underlying monogenic forms of bone development diseases, and for the collagen formation and extracellular matrix organization biological pathways. Ten of the likely effector genes, including TGFB1 (transforming growth factor beta 1), FGF18 (fibroblast growth factor 18), CTSK (cathepsin K), and IL11 (interleukin 11), have therapeutics approved or in clinical trials, with mechanisms of action supportive of evaluation for efficacy in osteoarthritis.

Evidence for genetic contribution to the increased risk of type 2 diabetes in schizophrenia.

The epidemiologic link between schizophrenia (SCZ) and type 2 diabetes (T2D) remains poorly understood. Here, we investigate the presence and extent of a shared genetic background between SCZ and T2D using genome-wide approaches. We performed a genome-wide association study (GWAS) and polygenic risk score analysis in a Greek sample collection (GOMAP) comprising three patient groups: SCZ only (n = 924), T2D only (n = 822), comorbid SCZ and T2D (n = 505); samples from two separate Greek cohorts were used as population-based controls (n = 1,125). We used genome-wide summary statistics from two large-scale GWAS of SCZ and T2D from the PGC and DIAGRAM consortia, respectively, to perform genetic overlap analyses, including a regional colocalisation test. We show for the first time that patients with comorbid SCZ and T2D have a higher genetic predisposition to both disorders compared to controls. We identify five genomic regions with evidence of colocalising SCZ and T2D signals, three of which contain known loci for both diseases. We also observe a significant excess of shared association signals between SCZ and T2D at nine out of ten investigated p value thresholds. Finally, we identify 29 genes associated with both T2D and SCZ, several of which have been implicated in biological processes relevant to these disorders. Together our results demonstrate that the observed comorbidity between SCZ and T2D is at least in part due to shared genetic mechanisms.

Genome-wide association study of developmental dysplasia of the hip identifies an association with GDF5.

Developmental dysplasia of the hip (DDH) is the most common skeletal developmental disease. However, its genetic architecture is poorly understood. We conduct the largest DDH genome-wide association study to date and replicate our findings in independent cohorts. We find the heritable component of DDH attributable to common genetic variants to be 55% and distributed equally across the autosomal and X-chromosomes. We identify replicating evidence for association between GDF5 promoter variation and DDH (rs143384, effect allele A, odds ratio 1.44, 95% confidence interval 1.34-1.56, P = 3.55 × 10-22). Gene-based analysis implicates GDF5 (P = 9.24 × 10-12), UQCC1 (P = 1.86 × 10- 10), MMP24 (P = 3.18 × 10-9), RETSAT (P = 3.70 × 10- 8) and PDRG1 (P = 1.06 × 10- 7) in DDH susceptibility. We find shared genetic architecture between DDH and hip osteoarthritis, but no predictive power of osteoarthritis polygenic risk score on DDH status, underscoring the complex nature of the two traits. We report a scalable, time-efficient recruitment strategy and establish for the first time to our knowledge a robust DDH genetic association locus at GDF5.

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

Osteoarthritis is a common complex disease imposing a large public-health burden. Here, we performed a genome-wide association study for osteoarthritis, using data across 16.5 million variants from the UK Biobank resource. After performing replication and meta-analysis in up to 30,727 cases and 297,191 controls, we identified nine new osteoarthritis loci, in all of which the most likely causal variant was noncoding. For three loci, we detected association with biologically relevant radiographic endophenotypes, and in five signals we identified genes that were differentially expressed in degraded compared with intact articular cartilage from patients with osteoarthritis. We established causal effects on osteoarthritis for higher body mass index but not for triglyceride levels or genetic predisposition to type 2 diabetes.

Evaluating the glucose raising effect of established loci via a genetic risk score.

Recent genome-wide association studies have identified several single nucleotide polymorphisms (SNPs) associated with glucose levels. We tested the hypothesis here whether the cumulative effect of glucose raising SNPs, assessed via a score, is associated with glucose levels. A total of 1,434 participants of Greek descent from the THISEAS study and 1,160 participants form the GOMAP study were included in this analysis. We developed a genetic risk score (GRS), based on the known glucose-raising loci, in order to investigate the cumulative effect of known glucose loci on glucose levels. In the THISEAS study, the GRS score was significantly associated with increased glucose levels (mmol/L) (ß ± SE: 0.024 ± 0.004, P = 8.27e-07). The effect of the genetic risk score was also significant in the GOMAP study (ß ± SE: 0.011 ± 0.005, P = 0.031). In the meta-analysis of the two studies both scores were significantly associated with higher glucose levels GRS: ß ± SE: 0.019 ± 0.003, P = 1.41e-09. Also, variants at the SLC30A8, PROX1, MTNR1B, ADRA2A, G6PC2, LPIN3 loci indicated nominal evidence for association with glucose levels (p < 0.05). We replicate associations of the established glucose raising variants in the Greek population and confirm directional consistency of effects (binomial sign test p = 6.96e-05). We also demonstrate that the cumulative effect of the established glucose loci yielded a significant association with increasing glucose levels.

The Genetic Epidemiological Landscape of Hip and Knee Osteoarthritis: Where Are We Now and Where Are We Going?

Osteoarthritis (OA) is a complex disease that affects the whole joint, with multiple biological and environmental factors contributing to its development. The heritable component for primary OA accounts for ∼50% of susceptibility. So far, candidate gene studies and genome-wide association scans have established 18 OA-associated loci. These findings account for 11% of the heritability, explaining a rather small fraction of the genetic component. To further unravel the genetic architecture of OA, the field needs to facilitate more precise phenotypic definitions, high genome coverage, and large sample metaanalyses, expecting the identification of rare and low frequency variants with potentially higher penetrance, and more accurate methods for calculating phenotype-genotype correlation. Expression analysis, epigenetics, and investigation of interactions can also help clarify the implicated transcriptional regulatory pathways and provide insights into further novel pathogenic OA mechanisms leading to diagnostic biomarker identification and new, more focused therapeutic disease approaches.

Genetic characterization of Greek population isolates reveals strong genetic drift at missense and trait-associated variants.

Isolated populations are emerging as a powerful study design in the search for low-frequency and rare variant associations with complex phenotypes. Here we genotype 2,296 samples from two isolated Greek populations, the Pomak villages (HELIC-Pomak) in the North of Greece and the Mylopotamos villages (HELIC-MANOLIS) in Crete. We compare their genomic characteristics to the general Greek population and establish them as genetic isolates. In the MANOLIS cohort, we observe an enrichment of missense variants among the variants that have drifted up in frequency by more than fivefold. In the Pomak cohort, we find novel associations at variants on chr11p15.4 showing large allele frequency increases (from 0.2% in the general Greek population to 4.6% in the isolate) with haematological traits, for example, with mean corpuscular volume (rs7116019, P=2.3 × 10(-26)). We replicate this association in a second set of Pomak samples (combined P=2.0 × 10(-36)). We demonstrate significant power gains in detecting medical trait associations.