Cartilage Acidic Protein 1 in Plasma Associates With Prevalent Osteoarthritis and Predicts Future Risk as Well as Progression to Joint Replacements : Results From the UK Biobank Resource

Funding Information: Supported by deCODE genetics/Amgen Inc. Publisher Copyright: © 2022 The Authors. Arthritis & Rheumatology published by Wiley Periodicals LLC on behalf of American College of Rheumatology.OBJECTIVE: The level of cartilage acidic protein 1 (CRTAC1) in plasma was recently discovered to be associated with osteoarthritis (OA) risk and progression to joint replacement in Iceland. This study was undertaken to validate these findings in an independent population. METHODS: In this study, 1,462 plasma proteins were measured in 54,265 participants from the UK Biobank on the Olink Explore platform. We analyzed the association of plasma proteins with prevalent OA, incident OA, and progression to joint replacement. We assessed the specificity of OA association through comparison of associations with inflammatory joint diseases and with previous joint replacement. RESULTS: The CRTAC1 protein showed the strongest association with prevalent knee OA (odds ratio [OR] 1.34 [95% confidence interval (95% CI) 1.27, 1.41]) and was associated with hip OA (OR 1.19 [95% CI 1.11, 1.28]). It predicted incident diagnosis of OA in the knee (hazard ratio [HR] 1.40 [95% CI 1.35, 1.46]) and hip (HR 1.25 [95% CI 1.19, 1.31]), as well as progression to joint replacement (HR 1.20 [95% CI 1.08, 1.33] for the knee and HR 1.22 [95% CI 1.08, 1.38] for the hip), while no association was found with inflammatory joint diseases. Individuals in the highest quintile of risk based on CRTAC1 level, age, sex, and body mass index had a 10-fold risk of knee or hip OA within 5 years compared to those in the lowest quintile. Adding aggrecan core protein (ACAN) and neurocan core protein (NCAN) to the model improved the prediction of OA but not joint replacement. Furthermore, we replicated the association of CUB domain-containing protein 1 with prior joint replacement. CONCLUSION: Plasma CRTAC1 is a specific biomarker for OA and a predictor of OA risk and progression to joint replacement. Adding ACAN and NCAN protein levels to the CRTAC1 model improved the prediction of OA.Peer reviewe

Multiple genetic loci for bone mineral density and fractures

To access publisher full text version of this article. Please click on the hyperlink in Additional Links fieldBACKGROUND: Bone mineral density influences the risk of osteoporosis later in life and is useful in the evaluation of the risk of fracture. We aimed to identify sequence variants associated with bone mineral density and fracture. METHODS: We performed a quantitative trait analysis of data from 5861 Icelandic subjects (the discovery set), testing for an association between 301,019 single-nucleotide polymorphisms (SNPs) and bone mineral density of the hip and lumbar spine. We then tested for an association between 74 SNPs (most of which were implicated in the discovery set) at 32 loci in replication sets of Icelandic, Danish, and Australian subjects (4165, 2269, and 1491 subjects, respectively). RESULTS: Sequence variants in five genomic regions were significantly associated with bone mineral density in the discovery set and were confirmed in the replication sets (combined P values, 1.2x10(-7) to 2.0x10(-21)). Three regions are close to or within genes previously shown to be important to the biologic characteristics of bone: the receptor activator of nuclear factor-kappaB ligand gene (RANKL) (chromosomal location, 13q14), the osteoprotegerin gene (OPG) (8q24), and the estrogen receptor 1 gene (ESR1) (6q25). The two other regions are close to the zinc finger and BTB domain containing 40 gene (ZBTB40) (1p36) and the major histocompatibility complex region (6p21). The 1p36, 8q24, and 6p21 loci were also associated with osteoporotic fractures, as were loci at 18q21, close to the receptor activator of the nuclear factor-kappaB gene (RANK), and loci at 2p16 and 11p11. CONCLUSIONS: We have discovered common sequence variants that are consistently associated with bone mineral density and with low-trauma fractures in three populations of European descent. Although these variants alone are not clinically useful in the prediction of risk to the individual person, they provide insight into the biochemical pathways underlying osteoporosis

Meta-analysis of erosive hand osteoarthritis identifies four common variants that associate with relatively large effect

[Abstract] Objectives: Erosive hand osteoarthritis (EHOA) is a severe subset of hand osteoarthritis (OA). It is unclear if EHOA is genetically different from other forms of OA. Sequence variants at ten loci have been associated with hand OA but none with EHOA. Methods: We performed meta-analysis of EHOA in 1484 cases and 550 680 controls, from 5 populations. To identify causal genes, we performed eQTL and plasma pQTL analyses, and developed one zebrafish mutant. We analysed associations of variants with other traits and estimated shared genetics between EHOA and other traits. Results: Four common sequence variants associated with EHOA, all with relatively high effect. Rs17013495 (SPP1/MEPE, OR=1.40, p=8.4×10-14) and rs11243284 (6p24.3, OR=1.35, p=4.2×10-11) have not been associated with OA, whereas rs11631127 (ALDH1A2, OR=1.46, p=7.1×10-18), and rs1800801 (MGP, OR=1.37, p=3.6×10-13) have previously been associated with hand OA. The association of rs1800801 (MGP) was consistent with a recessive mode of inheritance in contrast to its additive association with hand OA (OR homozygotes vs non-carriers=2.01, 95% CI 1.71 to 2.37). All four variants associated nominally with finger OA, although with substantially lower effect. We found shared genetic components between EHOA and other OA measures, grip strength, urate levels and gout, but not rheumatoid arthritis. We identified ALDH1A2, MGP and BMP6 as causal genes for EHOA, with loss-of-function Bmp6 zebrafish mutants displaying EHOA-like phenotypes. Conclusions: We report on significant genetic associations with EHOA. The results support the view of EHOA as a form of severe hand OA and partly separate it from OA in larger joints

Sequence variant at 8q24.21 associates with sciatica caused by lumbar disc herniation.

Efst á síðunni er hægt að nálgast greinina í heild sinni með því að smella á hlekkinn To access publisher's full text version of this article, please click on the hyperlink in Additional Links field or click on the hyperlink at the top of the page marked FilesLumbar disc herniation (LDH) is common and often debilitating. Microdiscectomy of herniated lumbar discs (LDHsurg) is performed on the most severe cases to resolve the resulting sciatica. Here we perform a genome-wide association study on 4,748 LDHsurg cases and 282,590 population controls and discover 37 highly correlated markers associating with LDHsurg at 8q24.21 (between CCDC26 and GSDMC), represented by rs6651255[C] (OR=0.81; P=5.6 × 10(-12)) with a stronger effect among younger patients than older. As rs6651255[C] also associates with height, we performed a Mendelian randomization analysis using height polygenic risk scores as instruments to estimate the effect of height on LDHsurg risk, and found that the marker's association with LDHsurg is much greater than predicted by its effect on height. In light of presented findings, we speculate that the effect of rs6651255 on LDHsurg is driven by susceptibility to developing severe and persistent sciatica upon LDH.European Commission National Institutes of Healt

Causal mechanisms and balancing selection inferred from genetic associations with polycystic ovary syndrome.

Polycystic ovary syndrome (PCOS) is the most common reproductive disorder in women, yet there is little consensus regarding its aetiology. Here we perform a genome-wide association study of PCOS in up to 5,184 self-reported cases of White European ancestry and 82,759 controls, with follow-up in a further ∼2,000 clinically validated cases and ∼100,000 controls. We identify six signals for PCOS at genome-wide statistical significance (P<5 × 10(-8)), in/near genes ERBB4/HER4, YAP1, THADA, FSHB, RAD50 and KRR1. Variants in/near three of the four epidermal growth factor receptor genes (ERBB2/HER2, ERBB3/HER3 and ERBB4/HER4) are associated with PCOS at or near genome-wide significance. Mendelian randomization analyses indicate causal roles in PCOS aetiology for higher BMI (P=2.5 × 10(-9)), higher insulin resistance (P=6 × 10(-4)) and lower serum sex hormone binding globulin concentrations (P=5 × 10(-4)). Furthermore, genetic susceptibility to later menopause is associated with higher PCOS risk (P=1.6 × 10(-8)) and PCOS-susceptibility alleles are associated with higher serum anti-Müllerian hormone concentrations in girls (P=8.9 × 10(-5)). This large-scale study implicates an aetiological role of the epidermal growth factor receptors, infers causal mechanisms relevant to clinical management and prevention, and suggests balancing selection mechanisms involved in PCOS risk.This work was supported by the Medical Research Council [U106179472; MC_U106179472; U106179471; MC_U106179471] and the National Human Genome Research Institute of the National Institutes of Health (grant number R44HG006981 to 23andMe). The UK Medical Research Council and Wellcome Trust (092731), together with the University of Bristol, provide core support for the ALSPAC study. AMH assays in ALSPAC were funded with a grant from the US National Institute of Health (R01 DK077659). DAL works in a unit that receives funding from the University of Bristol and the UK Medical Research Council (MC_UU_12013/5).This is the final version of the article. It first appeared from Nature Publishing Group via http://dx.doi.org/10.1038/ncomms946

European bone mineral density loci are also associated with BMD in East-Asian populations

To access publisher full text version of this article. Please click on the hyperlink in Additional Links fieldMost genome-wide association (GWA) studies have focused on populations of European ancestry with limited assessment of the influence of the sequence variants on populations of other ethnicities. To determine whether markers that we have recently shown to associate with Bone Mineral Density (BMD) in Europeans also associate with BMD in East-Asians we analysed 50 markers from 23 genomic loci in samples from Korea (n = 1,397) and two Chinese Hong Kong sample sets (n = 3,869 and n = 785). Through this effort we identified fourteen loci that associated with BMD in East-Asian samples using a false discovery rate (FDR) of 0.05; 1p36 (ZBTB40, P = 4.3×10(-9)), 1p31 (GPR177, P = 0.00012), 3p22 (CTNNB1, P = 0.00013), 4q22 (MEPE, P = 0.0026), 5q14 (MEF2C, P = 1.3×10(-5)), 6q25 (ESR1, P = 0.0011), 7p14 (STARD3NL, P = 0.00025), 7q21 (FLJ42280, P = 0.00017), 8q24 (TNFRSF11B, P = 3.4×10(-5)), 11p15 (SOX6, P = 0.00033), 11q13 (LRP5, P = 0.0033), 13q14 (TNFSF11, P = 7.5×10(-5)), 16q24 (FOXL1, P = 0.0010) and 17q21 (SOST, P = 0.015). Our study marks an early effort towards the challenge of cataloguing bone density variants shared by many ethnicities by testing BMD variants that have been established in Europeans, in East-Asians

Collaborative Meta-analysis: Associations of 150 Candidate Genes With Osteoporosis and Osteoporotic Fracture

To access publisher full text version of this article. Please click on the hyperlink in Additional Links fieldBACKGROUND: Osteoporosis is a highly heritable trait. Many candidate genes have been proposed as being involved in regulating bone mineral density (BMD). Few of these findings have been replicated in independent studies. OBJECTIVE: To assess the relationship between BMD and fracture and all common single-nucleotide polymorphisms (SNPs) in previously proposed osteoporosis candidate genes. DESIGN: Large-scale meta-analysis of genome-wide association data. SETTING: 5 international, multicenter, population-based studies. PARTICIPANTS: Data on BMD were obtained from 19 195 participants (14 277 women) from 5 populations of European origin. Data on fracture were obtained from a prospective cohort (n = 5974) from the Netherlands. MEASUREMENTS: Systematic literature review using the Human Genome Epidemiology Navigator identified autosomal genes previously evaluated for association with osteoporosis. We explored the common SNPs arising from the haplotype map of the human genome (HapMap) across all these genes. BMD at the femoral neck and lumbar spine was measured by dual-energy x-ray absorptiometry. Fractures were defined as clinically apparent, site-specific, validated nonvertebral and vertebral low-energy fractures. RESULTS: 150 candidate genes were identified and 36 016 SNPs in these loci were assessed. SNPs from 9 gene loci (ESR1, LRP4, ITGA1, LRP5, SOST, SPP1, TNFRSF11A, TNFRSF11B, and TNFSF11) were associated with BMD at either site. For most genes, no SNP was statistically significant. For statistically significant SNPs (n = 241), effect sizes ranged from 0.04 to 0.18 SD per allele. SNPs from the LRP5, SOST, SPP1, and TNFRSF11A loci were significantly associated with fracture risk; odds ratios ranged from 1.13 to 1.43 per allele. These effects on fracture were partially independent of BMD at SPP1 and SOST. Limitation: Only common polymorphisms in linkage disequilibrium with SNPs in HapMap could be assessed, and previously reported associations for SNPs in some candidate genes could not be excluded. CONCLUSION: In this large-scale collaborative genome-wide meta-analysis, 9 of 150 candidate genes were associated with regulation of BMD, 4 of which also significantly affected risk for fracture. However, most candidate genes had no consistent association with BMD

The sequences of 150,119 genomes in the UK Biobank

Detailed knowledge of how diversity in the sequence of the human genome affects phenotypic diversity depends on a comprehensive and reliable characterization of both sequences and phenotypic variation. Over the past decade, insights into this relationship have been obtained from whole-exome sequencing or whole-genome sequencing of large cohorts with rich phenotypic data(1,2). Here we describe the analysis of whole-genome sequencing of 150,119 individuals from the UK Biobank(3). This constitutes a set of high-quality variants, including 585,040,410 single-nucleotide polymorphisms, representing 7.0% of all possible human single-nucleotide polymorphisms, and 58,707,036 indels. This large set of variants allows us to characterize selection based on sequence variation within a population through a depletion rank score of windows along the genome. Depletion rank analysis shows that coding exons represent a small fraction of regions in the genome subject to strong sequence conservation. We define three cohorts within the UK Biobank: a large British Irish cohort, a smaller African cohort and a South Asian cohort. A haplotype reference panel is provided that allows reliable imputation of most variants carried by three or more sequenced individuals. We identified 895,055 structural variants and 2,536,688 microsatellites, groups of variants typically excluded from large-scale whole-genome sequencing studies. Using this formidable new resource, we provide several examples of trait associations for rare variants with large effects not found previously through studies based on whole-exome sequencing and/or imputation

Genome-wide association identifies seven loci for pelvic organ prolapse in Iceland and the UK Biobank.

To access publisher's full text version of this article, please click on the hyperlink in Additional Links field or click on the hyperlink at the top of the page marked DownloadPelvic organ prolapse (POP) is a downward descent of one or more of the pelvic organs, resulting in a protrusion of the vaginal wall and/or uterus. We performed a genome-wide association study of POP using data from Iceland and the UK Biobank, a total of 15,010 cases with hospital-based diagnosis code and 340,734 female controls, and found eight sequence variants at seven loci associating with POP (P 5%) and one with minor allele frequency of 4.87%. Some of the variants associating with POP also associated with traits of similar pathophysiology. Of these, rs3820282, which may alter the estrogen-based regulation of WNT4, also associates with leiomyoma of uterus, gestational duration and endometriosis. Rs3791675 at EFEMP1, a gene involved in connective tissue homeostasis, also associates with hernias and carpal tunnel syndrome. Our results highlight the role of connective tissue metabolism and estrogen exposure in the etiology of POP.UCL Hospitals NIHR Biomedical Research Centr

Novel Genetic Variants for Cartilage Thickness and Hip Osteoarthritis

Osteoarthritis is one of the most frequent and disabling diseases of the elderly. Only few genetic variants have been identified for osteoarthritis, which is partly due to large phenotype heterogeneity. To reduce heterogeneity, we here examined cartilage thickness, one of the structural components of joint health. We conducted a genome-wide association study of minimal joint space width (mJSW), a proxy for cartilage thickness, in a discovery set of 13,013 participants from five different cohorts and replication in 8,227 individuals from seven independent cohorts. We identified five genome-wide significant (GWS, P≤5·0×10−8) SNPs annotated to four distinct loci. In addition, we found two additional loci that were significantly replicated, but results of combined meta-analysis fell just below the genome wide significance threshold. The four novel associated genetic loci were located in/near TGFA (rs2862851), PIK3R1 (rs10471753), SLBP/FGFR3 (rs2236995), and TREH/DDX6 (rs496547), while the other two (DOT1L and SUPT3H/RUNX2) were previously identified. A systematic prioritization for underlying causal genes was performed using diverse lines of evidence. Exome sequencing data (n = 2,050 individuals) indicated that there were no rare exonic variants that could explain the identified associations. In addition, TGFA, FGFR3 and PIK3R1 were differentially expressed in OA cartilage lesions versus non-lesioned cartilage in the same individuals. In conclusion, we identified four novel loci (TGFA, PIK3R1, FGFR3 and TREH) and confirmed two loci known to be associated with cartilage thickness.The identified associations were not caused by rare exonic variants. This is the first report linking TGFA to human OA, which may serve as a new target for future therapies