Genome-wide meta-analysis of muscle weakness identifies 15 susceptibility loci in older men and women
=====================================================================================================

* Garan Jones
* Katerina Trajanoska
* Adam J Santanasto
* Najada Stringa
* Chia-Ling Kuo
* Janice L Atkins
* Joshua R Lewis
* ThuyVy Duong
* Shengjun Hong
* Mary L Biggs
* Jian’an Luan
* Chloe Sarnowski
* Kathryn L Lunetta
* Toshiko Tanaka
* Mary K Wojczynski
* Ryan Cvejkus
* Maria Nethander
* Sahar Ghasemi
* Jingyun Yang
* M. Carola Zillikens
* Stefan Walter
* Kamil Sicinski
* Erika Kague
* Cheryl L Ackert-Bicknell
* Dan E Arking
* B Gwen Windham
* Eric Boerwinkle
* Megan L Grove
* Misa Graff
* Dominik Spira
* Ilja Demuth
* Nathalie van der Velde
* Lisette C P G M de Groot
* Bruce M Psaty
* Michelle C Odden
* Alison E Fohner
* Claudia Langenberg
* Nicholas J Wareham
* Stefania Bandinelli
* Natasja M van Schoor
* Martijn Huisman
* Qihua Tan
* Joseph Zmuda
* Dan Mellström
* Magnus Karlsson
* David A Bennett
* Aron S Buchman
* Philip L De Jager
* Andre G Uitterlinden
* Uwe Völker
* Thomas Kocher
* Alexander Teumer
* Leocadio Rodriguéz-Mañas
* Francisco J García García
* José A Carnicero
* Pamela Herd
* Lars Bertram
* Claes Ohlsson
* Joanne M Murabito
* George A Kuchel
* Luigi Ferrucci
* David Melzer
* David Karasik
* Fernando Rivadeneira
* Douglas P Kiel
* Luke C Pilling

## Abstract

Low muscle strength is an important heritable indicator of poor health linked to morbidity and mortality in older people. In a genome-wide association study meta-analysis of 256,523 Europeans aged 60 years and over from 22 cohorts we identified 15 loci associated with muscle weakness (European Working Group on Sarcopenia in Older People definition: n=48,596 cases, 18.9% of total), including 12 loci not implicated in previous analyses of continuous measures of grip strength. Loci include genes reportedly involved in autoimmune disease (*HLA-DQA1* p=4*10−17), arthritis (*GDF5* p=4*10−13), cell cycle control and cancer protection, regulation of transcription, and others involved in the development and maintenance of the musculoskeletal system. Using Mendelian randomization we report possible overlapping causal pathways, including diabetes susceptibility, hematological parameters, and the immune system. We conclude that muscle weakness in older adults has distinct mechanisms from continuous strength, including several pathways considered to be hallmarks of ageing.

Keywords
*   Muscle weakness
*   strength
*   ageing
*   genetics
*   genome wide association study

## Introduction

Age-associated loss of muscle strength (termed dynapenia)1 is one of the characteristic changes occurring with advancing age, and muscle weakness is considered a fundamental component of frailty and sarcopenia2. Individuals over 70 years old typically demonstrate up to 20% lost muscle mass compared with individual in their twenties3. Although definitions of reduced muscle function in older people have focused on loss of muscle mass (“sarcopenia”) evidence now shows that muscle weakness itself is often more predictive of negative health outcomes4. Muscle weakness causes difficulties in daily functioning (i.e. disability) and low muscle strength (measured as hand grip strength, considered a biomarker of general dynapenia) is predictive of future morbidity and mortality3 over the long term5. Despite intensive research, causes of and contributors to muscle weakness in later life remain to be fully elucidated6. Importantly, muscle strength is heritable, and can thus be studied using genetic analysis.

Previously a genome-wide association study (GWAS) by the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) consortium identified two loci associated with maximum hand grip strength (as a quantitative trait) in 27,581 Europeans aged 65 and over7. Another study on maximum hand grip strength (divided by weight) in mostly middle-aged UK Biobank participants (334,925 people aged 40 to 70, mean aged 56) identified and replicated 64 loci, many of which are known to have a role in determining anthropometric measures of body size8,9. These previous studies that considered grip strength as a continuous phenotype across young and old individuals may not provide insights into the age related loss of muscle strength that leads to a magnitude of weakness sufficient to call it a disease.

Given the limited data on genetic contributions to a clinically meaningful level of muscle weakness in older adults, we aimed to determine the genetic variants and investigate causal pathways associated with low measured grip strength in 256,523 older adults (aged 60+ years) of European ancestries from the CHARGE consortium. The primary analysis was based on the established 2010 European Working Group on Sarcopenia in Older People (EWGSOP) definition of low grip strength, and results were compared to an analysis of the alternative Foundations of the National Institutes of Health (FNIH) definition based on its association with functional outcomes, with additional analyses stratified by sex.

## Results

### Study description

The meta-analysis comprised 256,523 individuals of European descent aged 60 years or older at assessment from 22 independent cohorts with maximum hand grip strength recorded – including the UK Biobank, the US Health and Retirement Study (HRS), the Framingham Heart Study (FHS), and others. 46,596 (18.9%) of all participants had muscle weakness (dynapenia) based on hand grip strength (EWGSOP definition: grip strength < 30 Kg Male; < 20 Kg Female). Individual study characteristics are described in the Supplementary Information and in Supplementary Table 1.

Our primary analysis of EWGSOP definition low grip strength will be described first, with subsequent additional analyses described in later sections: these include analysis of males and females separately and use of the alternative low grip strength criteria provided by the FNIH data.

### GWAS of low muscle strength identifies 15 loci

We found 15 genomic risk loci to be associated (p<5*10−8; 8 loci p<5*10−9) with EWGSOP definition low hand grip strength in our GWAS meta-analysis of 22 cohorts (N=256,523, n=48,596 cases), adjusted for age, sex, and technical covariates (Figure 1, Table 1). The strongest associations were with variants close to *HLA-DQA1* (rs34415150, beta/log-OR per G allele=0.0833, p=4.4*10−17), *GDF5* (rs143384, beta per A allele=0.0545, p=4.5*10−13) and *DYM* (rs62102286, beta per T allele=0.0487, p=5.5*10−11). Twelve of the fifteen lead SNPs from the GWAS have not previously been identified in studies of continuous grip strength in all ages 10(Table 2). This included the three most strongly associated variants near *HLA-DQA1* (previously implicated in rheumatoid arthritis: see Table 2), *GDF5* (‘Growth differentiation factor 5’: previously implicated in height, waist hip ratio, muscle mass, and osteoarthritis) and *DYM* (‘Dymeclin’: implicated in in height). Six other variants were previously linked to height and four to osteoarthritis. None were significantly (p<5*10−8) associated with lean muscle mass, although rs10952289 near AOC1 is nominally associated with appendicular lean muscle mass (p=6*10−4)11. The test of cohort heterogeneity in METAL for all 15 lead SNPs was not statistically significant (nominal het p>0.05). Full summary statistics for the meta-analysis are available for download ([http://www.kp4cd.org/dataset_downloads/mskkp](http://www.kp4cd.org/dataset_downloads/mskkp)).

View this table:
[Table 1:](http://medrxiv.org/content/early/2020/05/21/2020.05.14.20100354/T1)

Table 1: Fifteen genomic risk loci associated with low grip strength in 256,523 older men and women

![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/05/21/2020.05.14.20100354/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2020/05/21/2020.05.14.20100354/F1)

Figure 1: Manhattan plot of low grip strength genome-wide association study -log10 p-values
The p-values are from a fixed-effects meta-analysis of 256,523 Europeans aged 60 or older from 22 cohorts. The outcome was low hand grip strength grip strength cutoff (males <30 Kg, females <20 Kg). The x-axis is the chromosomal location, and the y-axis is the –log10 p-value for each genetic variant. The horizontal red line is the threshold for genome-wide significance (*p*<5*10−8). Fifteen genomic loci cross the threshold, and the lead variant (most significantly associated with low strength) is described in Table 1. The nearest gene is displayed for each locus.

View this table:
[Table 2:](http://medrxiv.org/content/early/2020/05/21/2020.05.14.20100354/T2)

Table 2: Low grip loci that appear in the GWAS catalog of published variant-trait associations

Overall, three of the fifteen identified lead variants (or proxies) have not previously been implicated in anthropometric or musculoskeletal phenotypes in the GWAS catalogue (see Table 2). This included: *RBBP6* (‘RB binding protein 6, ubiquitin ligase’: a tumor suppressor gene), *ALDH1A2* (‘Aldehyde Dehydrogenase 1 Family Member A2’: involved in the synthesis of retinoic acid), and a variant near *FTO* (‘FTO Alpha-Ketoglutarate Dependent Dioxygenase’: involved in the oxidative demethylation of different forms of RNA). Although the lead *ALDH1A2* SNP itself has not been identified in previous GWAS, other independent variants (R2<0.6) at the same locus (e.g. rs3204689) have been found to be associated with osteoarthritis (Supplementary Table 2).

The LDSC regression intercept was not significantly greater than 1, suggesting no bias due to population stratification or cryptic relatedness12 (see Supplementary Results for more information). The single nucleotide polymorphism (SNP) based heritability (*h2*) of low grip strength was 0.044 (SE 0.0027), i.e. 4.4%, by LD Score Regression.

In sex-stratified analysis there were eight significant genomic risk loci associated with EWGSOP low grip strength in females only (total N=132,443 with n=33,548 cases, 25.3%; see Supplementary Table 3). Seven of the eight loci were either present in the main analysis or were correlated with corresponding variants, for example rs201754 (chr13:51078446, Female only) and rs3118903 (chr13: 51099577, both sexes) with D’=0.97 and r2=0.91. rs7185040 (chr16: 2145787), mapped to gene *PKD1*, was only significant in the analysis of females only, although the association is borderline in the analysis of males and females together (females p = 3*10−8; combined analysis p = 5.5*10−8).

The analysis of males only (total N=118,371 with 13,327 cases, 11.3%) identified three genomic loci associated with the EWGSOP low grip strength definition. One of these variants, rs35225200, was in linkage disequilibrium with a lead SNP from the combined analysis (rs13107325 R2=0.89, D’=1.0). While the two remaining lead variants appeared to be distinct signals (see Supplementary Table 4 for details): rs774787160 mapped to gene *DSCAM* and was not associated with low grip strength in females (males p = 1*10−8; females p = 0.9) and rs145933237 mapped to mir-466, which was only nominally associated in females (males p = 2*10−8; females p=0.01). The reduced number of cases in males could contribute to the fewer associations observed, and the borderline significance (p-values close to 5*10−8). The results from sex-stratified analyses need to be interpreted with caution, as a recent pre-print on bioRxiv has shown that variants are associated with sex in cohorts such as the UK Biobank, likely due to their effect on participation13.

Analysis of 116 mitochondrial genetic variants (MAF>0.01) available in the UK Biobank directly genotyped microarray data identified two associated with EWGSOP-defined low hand grip strength (p<0.00043, i.e. Bonferroni-adjustment for mitochondrial variants). rs41518645 is a missense variant (p.Asp171Asn) in *MT-CYB*, identified in Plink logistic regression analysis (p=0.0003). rs201950015 is intronic, located between genes *CO1* and *ATP6/8* (p=0.00042). No variants achieved “genome-wide” significant (p<5*10−8). See Supplementary Table 5.

### Subset of overall grip strength genetic variants are associated with low grip strength

After reviewing the 64 published variants associated with continuous measures of grip strength in the UK Biobank cohort8, for all ages (40 to 70 years), we found that only three of these SNPs were significant genome-wide (p<5*10−8) in the EWGSOP low grip strength GWAS (Supplementary Table 6). The top association was rs13107325 (*SLC39A8* linear grip p=4.4*10−23, low grip strength meta-analysis p=7.4*10−11), then rs2430740 (*C12orf60* linear grip p=6*10−12, low grip strength meta-analysis p=2.6*10−9) and finally rs11236203 (*POLD3*, linear grip p=8.4*10−10, low grip strength meta-analysis p=2.8*10−8). Two less-significant associations (< 1*10−6) were seen with rs3821269 (*TGFA* linear grip p=3.5E*10−15, low grip strength meta-analysis p=8.0*10−8) and rs1556659 (*ENSG00000232985* linear grip p=1.1*10−11, low grip strength meta-analysis p=3.5*10−7).

A previous GWAS by the CHARGE consortium identified two loci associated with maximum hand grip strength recorded in 27,581 Europeans aged 65 or older7. We found that neither locus was associated with low grip strength in our analysis (rs752045 EWGSOP p=0.67, rs3121278 EWGSOP p=0.15).

### GWAS of low grip strength based on FNIH criteria

In secondary analysis we performed GWAS using the low grip strength definition published by the FNIH14. This criterion uses lower grip strength cut-offs (<26Kg for males and <16Kg for females) than the EWGSOP definition15, resulting in fewer cases (n=19,345, 7.6% of total). Five loci were significant in the analysis (p<5*10−8), only one of which was not identified in the EWGSOP low grip strength analysis described previously (see Supplementary Table 7, either the same SNP or in high LD with the EWGSOP lead SNP at that loci, for example rs3771501 and rs958685 R2=0.90). This single base-pair deletion (rs1403785912 – chr9:4284961:T:-) mapped to *GLIS3* (“GLIS family zinc finger 3” – a repressor and activator of transcription), and may be specifically associated with more “extreme” weakness (EWGSOP p-value = 1.2*10−3).

### Gene Expression and Pathways

We used data from the Genotype-Tissue Expression project (GTEx) v8 to identify whether the variants associated with low grip strength affect expression of genes (Table 1; Supplementary Table 8). Of the top 15 EWGSOP-associated variants 12 are eQTLs for at least one gene. For 8 of these, the nearest gene to the variant by chromosomal location is known to have altered expression, but other genes in the locus may also be affected. This is consistent with a recent study showing that the “nearest gene” is often a good candidate for being a causal pathway16. For the top two loci (*HLA-DQA1* and *GDF5*) the variants are eQTLs for these nearest genes, however for the *SLC39A8* locus the lead SNP (rs13107325) is not an eQTL for *SLC39A8*, but is an eQTL for *UBE2D3* (‘Ubiquitin Conjugating Enzyme E2 D3’) in the aorta.

In MAGMA analysis we found 80 GO processes enriched in low grip strength-associated genes (see Supplementary Table 9 for details), mainly involved in the immune system and antigen presentation.

Using MetaXcan17 we identified 24 genes with expression in skeletal muscle significantly enriched in the low grip strength GWAS, after Benjamini-Hochberg adjustment for multiple testing (Supplementary Table 10). We used the International Mouse Phenotyping Consortium database ([http://www.mousephenotype.org](http://www.mousephenotype.org)) to investigate the possible phenotypes associated with the genes highlighted by MetaXcan, with many having clear effects on relevant phenotypes such as “growth”, “lean mass”, “body weight”, “angiogenesis”, “thymus involution”, and “lipid metabolism” (Supplementary Table 11).

We used LDSC-SEG to determine tissue-specific gene expression and chromatin modification enrichment in the low grip strength GWAS results18. We found no significant enrichment for the genetic determinants of low grip strength in expression profiles and epigenetic changes after adjustment for multiple testing, although genes expressed in the tibial nerve were borderline enriched (Benjamini-Hochberg-adjusted p-value 0.1). See Supplementary Table 12 for details.

### Genetic correlations and Mendelian randomization

We assessed ten common age-related diseases for their genetic correlation with low grip strength (Figure 2) using published genome-wide summary statistics and LDSC12. The largest genetic overlap was with osteoarthritis (29.7% genetic correlation, SE 6.3%), but also strong positive correlations with coronary artery disease (19.5%, SE 3.0%), type-2 diabetes (15.8%, SE 3.1%) and rheumatoid arthritis (12.7%, SE 7.9%). We observed no significant genetic correlation after multiple-testing correction with the remaining diseases examined (osteoporotic fracture risk, Alzheimer’s disease, stroke, chronic kidney disease, breast cancer and colorectal cancer). We also determined genetic correlations with five anthropometric traits (Figure 2), and found significant positive correlations with waist:hip ratio (13.0%, SE 2.7%) and BMI (9.1%, SE 2.5%), i.e. greater adiposity correlated with weakness in 60+ year olds. Significant negative correlations were observed with lean muscle mass (whole body: –30.9%, SE 6.1%; and appendicular: –26.5%, SE 5.8%) and with height (−37.4%, SE 2.3%). See Supplementary Table 13 for details.

![Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/05/21/2020.05.14.20100354/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2020/05/21/2020.05.14.20100354/F2)

Figure 2: Low grip strength genetic correlations with ten common diseases and five anthropometric traits
Genome-wide genetic correlations between low muscle strength and published summary statistics for common age-related diseases and low muscle strength risk factors. Full results available in Supplementary Table 13.

We examined 83 traits in Mendelian randomization analysis to find evidence for shared causal pathways with weakness (low grip EWGSOP) at older ages: primary results presented are betas from inverse variance-weighted regression using the ‘TwoSampleMR’ R package19 (Figure 3; Supplementary Table 14). We found significantly increased likelihood of weakness (multiple testing-adjusted p-values<0.05) with genetically predicted rheumatoid arthritis (Odds ratio=1.03, Benjamini-Hochberg adjusted p-value=5.3*10−4), presence of type-2 diabetes (OR=1.05, BH p = 2.5*10−3), or asthma and allergic disease (OR=1.07, BH p=9.4*10−3) (Supplementary Table 14). Genetic predisposition to greater age of menarche (OR=0.92, BH p=5.3*10−4), birth weight (OR=0.80, BH p = 5.3*10−4) and waist-hip ratio (WHR) adjusted for BMI in women only (OR=0.85, BH p=5.3*10−4) were protective of low grip strength as defined by the EWGSOP definition in both sexes (after adjustment for multiple testing). For each significant analysis we also examined the results from the weighted-median and MR-Egger tests to check consistency and for horizontal pleiotropy. Only birth weight had an MR-Egger beta that was inconsistent with the main effect (−0.03 compared to –0.2), although the intercept did not significantly deviate from 0 (p=0.1) and the MR-Egger confidence intervals overlap the IVW effect (95% CIs –0.33 to 0.27). Additionally, the WHR association should be interpreted with caution, as the analysis of WHR variants associated with both sexes were not statistically significant (nominal IVW p=0.01).

![Figure 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/05/21/2020.05.14.20100354/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2020/05/21/2020.05.14.20100354/F3)

Figure 3: Traits sharing causal pathways with low grip strength at older agesidentified in Mendelian Randomization analysis
In Mendelian randomization analyses we estimated whether 83 exposures may share causal pathways with low grip strength in people aged 60 and older. Those identified as significant (multiple testing-adjusted p<0.05) in at least one analysis (all participants, or males or females separately) are included in the figure. *WHR (adjBMI Women) = Waist-Hip Ratio SNPs identified in GWAS analysis of females only, adjusted for Body Mass Index. Full results of Mendelian randomization available in Supplementary Table 14.

The analysis of females only also identified depression (OR=1.21, BH p=2.75*10−2), and colorectal cancer (OR=1.06, BH p=2.75*10–2) (Supplementary Table 15). Although the MR-Egger intercept for depression was not statistically different from 0, there is potential for horizontal pleiotropy confounding this result (seen on Supplementary Figure 3).

In the males, type-2 diabetes significantly increased odds of EWGSOP low grip (OR=1.08, BH p=3.0*10−3) whereas greater plateletcrit (the volume occupied by platelets in the blood as a percentage) (OR=0.91, p=3.0*10−3) and other hematological parameters appear to be protective (Supplementary Table 16).

To explore the effect of genetic predisposition to low grip strength at older ages we created an unweighted genetic risk score (GRS) in the UK Biobank European sample by summing the number of low grip strength-associated alleles (15 genetic variants so 30 alleles, mean number of alleles = 12.6, SD=2.3). We first confirmed the association with low grip strength in UK Biobank participants (OR per allele 1.036: 95% CIs 1.030 to 1.041, p=6*10−40). The low grip GRS was also associated with increased Frailty Index20 (increase in points per allele = 0.013: 0.006 to 0.021, p=4*10−4).

### Low grip strength loci independence from musculoskeletal traits and diseases

To determine whether the genetic variants associated with low grip strength identified in the GWAS were independent of anthropometric traits or musculoskeletal co-morbidities we performed regression analyses in the UK Biobank cohort with adjustment for the following covariates: height, weight, skeletal muscle mass (determined using bioimpedance analysis), osteoarthritis, Rheumatoid arthritis, osteoporosis, Dupytrens contracture (one or more fingers permanently bent), and rhizarthrosis (arthritis of the thumb). The association between eight of the fifteen EWGSOP loci and low grip strength was attenuated after adjusting for height, including rs143384 (initial UKB p=3.7*10−11; adjusted UKB p=3.8*10−2) and rs7624084 (initial UKB p=9.3*10−7; adjusted UKB p=4.9*10−1). Adjustment for weight or BMI did not substantially attenuate any of the associations. Overall, the associations were not attenuated by adjustment for osteoarthritis, Rheumatoid arthritis, osteoporosis, Dupytrens, or rhizarthrosis. See Supplementary Information for diagnostic codes, and Supplementary Table 17 for detailed results.

## Discussion

In this study of 256,523 Europeans aged 60 years and over we identified 15 genetic loci – 3 novel –associated with the EWGSOP definition of low grip strength (dynapenia), plus two additional loci for the FNIH definition that used a more extreme definition for muscle weakness. Only three of these were previously implicated in previous analyses of continuous strength measures, suggesting that the genetic causes of clinically meaningful weakness at older ages are partly distinct. We found prominent overlaps with osteoarthritis and Rheumatoid arthritis, but also with cardiovascular disease and type-2 diabetes. Additional links to asthma and allergy were also found. The novel pathways implicated appear to include hallmark mechanisms of ageing21, for example cell cycle control related to the cancer control retinoblastoma pathway. However other ageing pathways such as telomere length22, and many lifespan-associated loci including *APOE*23, were not associated. Three of the low grip strength-associated genetic signals identified have not appeared in GWAS prior to the time of analysis, further demonstrating that low strength in older people may have distinct genetic underpinnings.

The strongest association found was rs34415150, near the *HLA-DQA1* gene. Genetic variants at this locus have been implicated in a wide range of conditions, including autoimmune diseases such as Rheumatoid arthritis24, and continuous grip strength8. HLA haplotypes *HLA-DQA1**03:01 and *HLA-DRB1**04:01, have been previously linked to sarcopenia in older UK Biobank participants25. HLA-DQA1 is associated with chronic inflammation in muscle of untreated children with juvenile dermatomyositis (inflammatory myopathies in children, which one of the characteristics is muscle weakness)26. Additionally, in a multi-trait analysis of age-related diseases *HLA-DQA1* was confirmed and it may therefore contribute to underlying ageing mechanisms as a “geroscience locus”27.

Overall five of the fifteen genomic risk loci for EWGSOP low grip strength have been previously associated, or in linkage disequilibrium with variants associated with osteoarthritis (rs143384 – *GDF5*, rs13107325 – *SLC39A8*, rs34464763 – *C12orf60*, rs958685 – *TGFA* and rs79723785 – *BRSK1*), of which two are also linked to adiposity, a known risk factor for arthritis28. We found that rs143384 in the 5’ untranslated region of growth/differentiation factor 5 (*GDF5*) was the second most strongly associated variant with low grip strength. *GDF5* is a protein in the transforming growth factor beta (TGF-β) family, with key roles in bone and joint development29,30. It was the first locus identified for osteoarthritis31, with a reported odds ratio of 1.79, as well as one of the first identified for height32. GDF5 is known to reduce expression of cartilage extracellular matrix-degrading enzymes in human primary chondrocytes33, thereby may be an potential intervention target for avoiding weakness at older ages, although work is needed to determine when intervention would be most effective. In follow-up analyses we found that the association between rs143384 and low grip strength was independent of prevalent osteoarthritis in UK Biobank participants, although we cannot rule out an effect of sub-clinical osteoarthritis. However, the association was completely attenuated after adjustment for standing height, suggesting the effect of the variant is mediated by developmental traits such as bone length. Mice with loss of *Gdf5* function exhibited severely impaired knee development, and the regulatory region pinpointed to mediate this effect in humans includes osteoarthritis-associated genetic variants34. Although the observed association between variants mapping to *GDF5* and low strength might be mediated by hand osteoarthritis pain compromising grip strength, there is evidence of a direct effect of GDF5 on muscle35.

The locus on Chr 18 was near to the DYM gene. Loss of function mutations in this gene are associated with Dyggve-Melchior-Clausen syndrome and Smith-McCort dysplasia respectively. Mice lacking this gene present with chondrodysplasia resulting from impaired endochondrial bone formation and abnormalities of the growth plate that begin to manifest shortly birth36. In homozygous mutant mice, and in patients with loss of function mutations in this gene, both the axial and the appendicular skeleton are affected. As is noted in Table 2, this gene is also associated with height in the GWAS catalogue. Interestingly, this gene is highly expressed in skeletal muscle in humans37, however its function in muscle is not completely understood.

*SLC39A8* encodes for the metal ion transporter ZIP8 which has been shown to be upregulated in chondrocytes present in osteoarthritic cartilage38. The lead SNP rs13107325 is a missense variant within *SLC39A8* which has been previously associated with osteoarthritis39.

On chromosome 12 genetic variants known to affect expression of *MGP* (Matrix Gla Protein) in the tibial nerve (among other tissues) are associated with osteoarthritis of the hand but not of the hip or knee40. MGP is an inhibitor of arterial and soft tissue calcifications, with links to atherosclerosis41. Consistent with this, older women with severe abdominal aortic calcification have greater decline in grip strength over 5 years42. More recently a study suggested MGP may also regulate muscle development and atrophy43.

We also identified variants known to affect Transforming Growth Factor Alpha (*TGFA*) expression (increased in the testis, decreased in the brain), which is implicated in cell proliferation, differentiation and development. This locus has previously been identified for overall strength9, and suggestive evidence from a study of 1,323 participants linked rs2862851, a variant in linkage disequilibrium with the lead SNP rs958685 (R2=0.90, D’=1.0), with increased risk of osteoarthritis in the knee (OR=1.4, p=3.1*10−4)44.

The low grip strength locus on chromosome 15 is near *ALDH1A2*, which has a key role in the pathogenesis of osteoarthritis45. Low grip strength-associated variants at this locus have previously been identified for severe osteoarthritis of the hand45. The lack of this gene in mice is perinatally lethal, however at embryonic day 18.5 mice lacking this gene do present with numerous cartilage gene defects46. We also identified variants that affect expression of *CHRDL2* (Chordin Like 2) in the thyroid, known to interact with mouse Gdf5, which is upregulated in human osteoarthritic joint cartilage cell line47. In addition, down-regulation of *CHRDL2* expression has been linked to the progression of severe osteoarthritis in the knee joint48, suggesting a role for CHRDL2 in cartilage repair.

Two of the identified loci (*RBBP6* ‘RB Binding Protein 6, Ubiquitin Ligase’ and *ZBTB38* ‘Zinc Finger And BTB Domain Containing 38’) form an axis involved in DNA replication and chromosomal stability49. RBBP6 ubiquitinates the transcriptional repressor ZBTB38, destabilizing it and reducing its action on the replication factor MCM10. In mice, *Zbtb38* is highly expressed in skeletal muscle, loss of this methyl-CpG-binding protein (which is also known as *Cibz*), promotes myogenic differentiation. Conversely, expression of *Zbtb38* is decreased in satellite cells during muscle regeneration50, suggesting that like other members of this gene family, this gene is involved in cellular differentiation. Variants associated with low grip strength in our analysis are known to decrease ZBTB38 expression in whole blood and skeletal muscle (among others) and increase expression in the skin.

Mitochondrial dysfunction is a hallmark of aging, yet we found no variants associated with low strength at genome-wide significance levels. Variants in *MT-CYB* were nominally significant (p=0.0003): *MT-CYB* (mitochondrial cytochrome b) is part of the mitochondrial respiratory chain, and essential for Complex III formation. Monogenic diseases associated with *MT-CYB* include exercise intolerance and “Additional features include lactic acidosis, muscle weakness and/or myoglobinuria”51. Recent work by Cohen *et al*. have highlighted the importance of mitochondrial peptides such as humanin in many age-related diseases52, however we found no variants in these genes associated with low grip strength passing multiple testing correction.

We observed minimal overlap between loci associated with low grip strength and general anthropometric traits such as height and continuous measures of strength. The SNP-based heritability estimate for EWGSOP low grip strength in older adults was 4.4% (SE 0.3%). This compares to 13% (SE 0.4%) SNP-based heritability for continuous grip strength in UK Biobank participants aged 40 to 708. This may be partly explained by our using a binary cut-off for low grip compared to the quantitative grip analysis. Additionally, we observed only limited negative genetic correlation between low grip strength and lean muscle mass (∼30%): taken together, these results emphasize that the genetics of weakness and overall strength are distinct.

Despite seeing little overlap at the individual locus level between low grip strength at older ages and other diseases we did observe significant genetic correlations, especially with osteoarthritis (30% overlap). However in Mendelian Randomization analysis we did not observe a causal relationship between osteoarthritis and low grip strength: taken together, this suggests that osteoarthritis shares causal risk factors and biological pathways with low grip strength at older ages, such as obesity, but may not cause it. Additionally, osteoarthritis has diverse loci associated with different joints, for instance osteoarthritis in the hand appears to have a distinct genetic signal: the low grip strength locus associated with *MGP* expression has been previously linked to osteoarthritis in the fingers and hand, but not the hip or knee40. Although our results were robust to adjustment for osteoarthritis (including rhizarthrosis – arthritis of the thumb) this may suggest that arthritis in the hand needs to be accounted for in measures of muscle weakness, as hand grip strength may not always reflect muscle strength elsewhere, e.g. lower-extremity strength.

Our Mendelian Randomization analyses highlighted specific traits and diseases which may share causal pathways with weakness at older ages. This included growth and development traits such as birth weight, waist:hip ratio, and pubertal timing (age at menarche) in women (highly genetically correlated – 75% – with age at voice breaking in men53), where greater values were protective. Although puberty timing is highly polygenic, it is strongly genetically correlated with BMI (−35%)53. These results are consistent with the observation that growth and development traits are associated with strength trajectories in later life54. We chose to not adjust for body size in our primary analysis in case interesting or novel effects were masked, but in follow-up analysis determined that all except two of the loci identified were predominantly independent of participant height and weight. Raised red blood cell parameters – especially plateletcrit, the proportion of blood occupied by platelets –appear to be protective in males but associations were attenuated or non-significant in females. A number of studies have recently reported a link between raised platelet counts, inflammation, and sarcopenia cross-sectionally55. However our results suggest that plateletcrit across the lifecourse (rather than after sarcopenia onset) may be different. Lastly, only four of the conditions we investigated (which included coronary artery disease, and some common cancers) appear to causally increase risk of weakness in older people: depression, asthma/allergic diseases, Rheumatoid arthritis, and type-2 diabetes. These are diverse conditions and further underlines the multifactorial causes of weakness in older people.

There are a number of limitations to our analyses. The data included are predominantly from subjects at the younger end of the age 60 plus demographic, and is not enriched for frail individuals. The sample size for sex-specific analyses is limited, especially for men. Our analysis was limited to relatively common variants (prevalence > 1%) in subjects with European ancestries only. The analyses of the FNIH strength cutpoints also have limited power, given the low prevalence of the phenotype, although studying the extremes of a continuous trait can provide increased power if stronger associations are uncovered. We did not include additional analysis of the revised EWGSOP2 low grip cut-points6 as these are almost identical to the FNIH criteria, which many cohorts had already analyzed; we future analyses should include this. Analysis of rare and structural variants, and analyses in other ancestral groups will give a more complete picture of the genetic landscape for low grip strength or dynapenia. Some Mendelian Randomization analyses have limited power due to the lack of strong instruments, and therefore null results for these analyses should be interpreted with caution.

To conclude, genetic variation in 15 loci are related to muscle weakness in people aged 60 plus, of European descent, with limited overlap with loci associated with the full range of muscle strength in 40 to 70 year olds. Novel loci implicated may be involved in hallmark pathways of ageing including cell cycle control and inflammation, along with loci implicated in arthritis and pathways involved in the development and maintenance of the musculoskeletal system.

## Methods

### GWAS of low grip strength in older people

We conducted a GWAS meta-analysis of low grip strength in participants aged 60 years or older of European ancestry from 22 studies yielding a combined sample of 254,894 individuals. Individual studies used different genotyping platforms and imputation was predominantly performed using the Haplotype reference consortium (HRC) v1.1 panel. See Supplementary Information for details on individual study methods.

Two definitions of low muscle hand grip strength were utilized at the time of analysis. The primary analysis was of the 2010 EWGSOP criteria for sarcopenic grip strength (Grip strength < 30 Kg Male; < 20 Kg Female). In secondary analysis we considered a more data-driven definition with more extreme thresholds by the Foundation for the National Institutes of Health (FNIH) sarcopenia project 2014 (Grip strength < 26 Kg Male; < 16 Kg Female) for comparison. Given the known differences in strength between males and females (on average) we also performed sex stratified analyses.

GWAS was performed by each cohort individually (see Supplementary Methods) using regression models, adjusted for age, sex (except in sex-specific models), and population substructure, accounting for relatedness and technical covariates as required by the individual study. No adjustment for anthropometric measures was made in the primary analysis, but the effects were explored in sensitivity analyses (see below). Fixed-effects inverse variance weighted meta-analysis was performed using METAL56 using the GWAS summary statistics generated by each cohort, with genomic control for population structure (see Supplementary Methods for details). The following quality control filters were applied: minor allele frequency (MAF) > 0.01, imputation info score of > 0.4, and the variant present in at least two studies (UK Biobank – the largest included cohort – plus at least one other). The final analysis therefore included 9,678,524 genetic variants. Associations that achieved a p<5*10−8 were considered statistically significant, with those reaching the more stringent threshold of p<5*10−9 highlighted.

Distinct loci were initially defined as two significant variants separated by > 500kb. To identify independent signals at each locus we used FUMA (Functional Mapping and Annotation of Genome-Wide Association Studies)57, which uses Linkage Disequilibrium (LD) information to determine independence (r2 threshold = 0.1 for independent significant single nucleotide polymorphisms (SNP)). We used Linkage Disequilibrium Score Regression (LDSC, v1.0.0) to estimate the level of bias (i.e. from population stratification and cryptic relatedness) in the GWAS, and the SNP-based heritability of low grip strength12.

### Locus overlap with diseases and anthropometric traits

The GWAS catalog of published loci-trait associations10 was searched to identify whether low grip strength-associated loci determined from our meta-analysis are known to influence other traits or diseases. In addition, we performed sensitivity analyses in the UK Biobank sample to determine whether associations between variants and low grip strength identified in the meta-analysis were robust to adjustment for the following traits or diseases: height, weight, body mass index, osteoarthritis, rheumatoid arthritis, and osteoporosis (prevalent diseases were from the self-reported data at baseline or in the Hospital Episode Statistics). Analyses were performed in STATA (v 15) using logistic regression models adjusted for age, sex, principal components 1 to 10, assessment center, and genotyping array (the UK Biobank using two different Affymetrix microarrays that shared > 95% of sites – see Supplementary Methods).

### Gene Ontology Pathways, Tissue Enrichment, and eQTL analyses

We utilized FUMA to perform functional interpretation of the GWAS results57. In particular, FUMA performs gene-set analysis (using Multi-marker Analysis of GenoMic Annotation (MAGMA)58) to identify pathways enriched amongst the significant genes (weighted by the SNP-associations in proximity to them), in addition to searching eQTL databases to identify SNPs that significantly alter the expression of genes in various tissues.

We used MetaXcan to determine whether gene-level transcriptomic associations in GTEx v7 skeletal muscle data were enriched in the GWAS summary statistics for low grip strength17. The analysis included 7,512 genes with measured expression in the dataset; we applied Benjamini-Hochberg multiple-testing correction, with adjusted p-values < 0.05 deemed to be significant.

LD Score Regression applied to specifically expressed genes (LDSC-SEG) allows the identification of enriched tissue activity associated with GWAS results18. We applied LDSC-SEG (v1.0.0) to the GWAS summary statistics using the datasets ‘Multi\_tissue_gene_expr’ and ‘Multi_tissue_chromatin’ provided by the authors. We applied Benjamini-Hochberg multiple-testing correction for the 703 tests (n gene expression = 205, n chromatin = 498), with adjusted p-values < 0.05 deemed to be significant.

### Genetic correlations and Mendelian randomization

We investigated the genetic correlation between the low grip strength trait and 10 diseases – chosen because they are common, chronic diseases of aging – using LDSC (v1.0.0)12 and published GWAS summary statistics for the following diseases: Alzheimer’s disease59, breast cancer60, chronic kidney disease61, coronary artery disease62, osteoporotic fracture risk63, osteoarthritis64, prostate cancer65, rheumatoid arthritis66, stroke67, and type-2 diabetes68. We also calculated genetic correlations with the following anthropometric traits: height69, body mass index (BMI)69, waist:hip ratio (WHR)70, whole-body lean mass71, and appendicular lean mass71.

We also undertook a phenotype-wide Mendelian randomization (MR) association study to examine the causal effect of 83 traits on low hand grip strength. We used the ‘TwoSampleMR’ (v0.4.23) package in R19 to perform the analysis of genetic instruments from the 83 traits, which including those traits with clear biological rationale (for example, adiposity) and others that are more exploratory for hypothesis generation (for example, puberty timing). Follow up sensitivity analysis of the identified traits was by the MR-Egger and using weighted median estimation methods provided in the package.

## Data Availability

The GWAS summary statistics and supporting information on low grip strength in older people are available on the Musculoskeletal Knowledge Portal (http://www.kp4cd.org/dataset_downloads/mskkp). All relevant additional data is available on request from the authors or individual studies.

[http://www.kp4cd.org/dataset\_downloads/mskkp](http://www.kp4cd.org/dataset_downloads/mskkp) 

## Data Availability

The GWAS summary statistics and supporting information on low grip strength in older people are available on the Musculoskeletal Knowledge Portal ([http://www.kp4cd.org/dataset_downloads/mskkp](http://www.kp4cd.org/dataset_downloads/mskkp)). All relevant additional data is available on request from the authors or individual studies.

## Author Information

Garan Jones, Katerina Trajanoska, Adam J. Santanasto, Najada Stringa, Chia-Ling Kuo, and Janice L. Atkins contributed equally to this work.

Joanne M. Murabito, George A. Kuchel, Luigi Ferrucci, David Melzer, David Karasik, Fernando Rivadeneira, Douglas K. Kiel, and Luke C. Pilling jointly supervised this work.

## Contributions

D.E.A., B.G.W., E.B., M.L.G., M.G., I.D., N.v.d.V, L.C.P.G.M.d.G, B.M.P., C.L., N.J.W., S.B., N.M.v.S, M.H, J.M.Z., D.M., M.K., D.A.B., A.S.B., P.L.D.J., A.G.U., A.T., L.R.M., F.G.G., J.C., P.H., L.B., C.O., J.M.M, L.F., D.K., F.R., and D.P.K. designed and managed individual studies. B.G.W., N.M.v.S, M.H, D.A.B., A.S.B., P.L.D.J., L.R.M., F.G.G., J.C., and P.H. collected data. G.J., K.T., A.J.S., N.S, J.R.L, J.M.M, G.A.K, L.F., D.M., D.K., F.R., D.P.K. and L.C.P. reviewed the analysis plan. G.J., K.T., A.J.S., N.S, C.K., T.V.D., S.H., M.L.B., J.L., C.S., K. L. L., T.T, M.K.W., R.C., M.N., S.G., J.Y., M.C., S.W., K.S., A.T., and L.C.P. analyzed the data. G.J. and L.C.P. performed the meta-analysis. G.J., C.A-B and L.C.P. performed the pathway and other analyses. J.M.M, G.A.K, L.F., D.M., D.K., F.R., D.P.K., and L.C.P. supervised the overall study design. G.J., K.T., A.J.S., N.S, C.K., J.L.A., J.M.M, G.A.K, L.F., D.M., D.K., F.R., D.P.K. and L.C.P. wrote the manuscript. G.J., K.T., A.J.S., N.S, C.K., J.L.A., J.R.L, T.V.D., S.H., M.L.B., J.L., C.S., K. L. L., T.T, M.K.W., R.C., M.N., S.G., J.Y., M.C., S.W., K.S., E.K., C.A-B., D.E.A., B.G.W., E.B., M.L.G., M.G., D.S., I.D., N.v.d.V, L.C.P.G.M.d.G, B.M.P., M.C.O., A.E.F., C.L., N.J.W., S.B., N.M.v.S, M.H, Q.T., J.M.Z., D.M., M.K., D.A.B., A.S.B., P.L.D.J., A.G.U., U.V., T.K., A.T., L.R.M., F.G.G., J.C., P.H., L.B., C.O., J.M.M, G.A.K, L.F., D.M., D.K., F.R., D.P.K. and L.C.P. reviewed the manuscript.

## Competing Interests

Dr Kiel is a consultant for Solarea Bio, received institutional grants from Amgen and Radius Health, and received royalties for publication Wolters Kluwer. Dr. Odden serves as a consultant for Cricket Health, a kidney care company. Dr. Psaty serves on the Steering Committee of the Yale Open Data Access Project funded by Johnson & Johnson. All other authors declare no competing financial interests.

## Acknowledgements

This study was part funded by an award to DM by the UK Medical Research Council (MR/M023095/1). A full list of acknowledgements and grant support can be found in the Supplementary Information. The authors would like to acknowledge the use of the University of Exeter High-Performance Computing (HPC) facility in carrying out this work.

*   Received May 14, 2020.
*   Revision received May 14, 2020.
*   Accepted May 19, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.Clark, B. C. & Manini, T. M. What is dynapenia? Nutrition 28, 495–503 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.nut.2011.12.002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22469110&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000303095200003&link_type=ISI) 

2.  2.Manini, T. M. & Clark, B. C. Dynapenia and Aging: An Update Journals Gerontol. Ser. A 67A, 28–40 (2012).
    
    

3.  3.Mitchell, W. K. et al. Sarcopenia, Dynapenia, and the Impact of Advancing Age on Human Skeletal Muscle Size and Strength; a Quantitative Review. Front. Physiol. 3, 260 (2012).
    
    

4.  4.Cawthon, P. M. et al. Establishing the Link Between Lean Mass and Grip Strength Cut Points With Mobility Disability and Other Health Outcomes: Proceedings of the Sarcopenia Definition and Outcomes Consortium Conference. Journals Gerontol. Ser. A (2019). doi:10.1093/gerona/glz081
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/gerona/glz081&link_type=DOI) 

5.  5.Rantanen, T. et al. Midlife hand grip strength as a predictor of old age disability. J. Am. Med. Assoc. 281, 558–560 (1999).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.281.6.558&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10022113&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000078452900035&link_type=ISI) 

6.  6.Cruz-Jentoft, A. J. et al. Sarcopenia: revised European consensus on definition and diagnosis. Age Ageing 48, 16–31 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ageing/afy169&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

7.  7.Matteini, A. M. et al. GWAS analysis of handgrip and lower body strength in older adults in the CHARGE consortium. Aging Cell 15, 792–800 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/acel.12468&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27325353&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

8.  8.Tikkanen, E. et al. Biological Insights Into Muscular Strength: Genetic Findings in the UK Biobank. Sci. Rep. 8, 6451 (2018).
    
    

9.  9.Willems, S. M. et al. Large-scale GWAS identifies multiple loci for hand grip strength providing biological insights into muscular fitness. Nat. Commun. 8, 16015 (2017).
    
    

10. 10.Buniello, A. et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 47, D1005–D1012 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gky1120&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30445434&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

11. 11.Zillikens, M. C. et al. Large meta-analysis of genome-wide association studies identifies five loci for lean body mass. Nat. Commun. 8, 80 (2017).
    
    

12. 12.Bulik-Sullivan, B. K. et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat Genet advance on, 291–295 (2015).
    
    

13. 13.Pirastu, N. et al. Genetic analyses identify widespread sex-differential participation bias. bioRxiv (unpublished Prepr. (2020). doi:10.1101/2020.03.22.001453
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wMy4yMi4wMDE0NTN2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA1LzIxLzIwMjAuMDUuMTQuMjAxMDAzNTQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

14. 14.Studenski, S. a et al. The FNIH sarcopenia project: rationale, study description, conference recommendations, and final estimates. J. Gerontol. A. Biol. Sci. Med. Sci. 69, 547–58 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/gerona/glu010&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24737557&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000334738000007&link_type=ISI) 

15. 15.Cruz-Jentoft, A. J. et al. Sarcopenia: European consensus on definition and diagnosis: Report of the European Working Group on Sarcopenia in Older People. Age Ageing 39, 412–423 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ageing/afq034&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20392703&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000278968300003&link_type=ISI) 

16. 16.Stacey, D. et al. ProGeM: a framework for the prioritization of candidate causal genes at molecular quantitative trait loci. Nucleic Acids Res. 47, e3–e3 (2019).
    
    

17. 17.Barbeira, A. N. et al. Exploring the phenotypic consequences of tissue specific gene expression variation inferred from GWAS summary statistics. Nat. Commun. 9, 1–20 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-02974-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29317637&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

18. 18.Finucane, H. K. et al. Heritability enrichment of specifically expressed genes identifies disease-relevant tissues and cell types. Nat. Genet. 50, 621–629 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0081-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29632380&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

19. 19.Hemani, G. et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife 7, (2018).
    
    

20. 20.Williams, D. M., Jylhava, J., Pedersen, N. L. & Hagg, S. A frailty index for UK Biobank participants. Journals Gerontol. Med. Sci. 74, 582–587 (2019).
    
    

21. 21.López-Otín, C., Blasco, M. A., Partridge, L., Serrano, M. & Kroemer, G. The hallmarks of aging. Cell 153, 1194–217 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2013.05.039&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23746838&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000319979200007&link_type=ISI) 

22. 22.Kuo, C.-L., Pilling, L. C., Kuchel, G. A., Ferrucci, L. & Melzer, D. Telomere length and aging-related outcomes in humans: A Mendelian randomization study in 261,000 older participants. Aging Cell e13017 (2019). doi:10.1111/acel.13017
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/acel.13017&link_type=DOI) 

23. 23.Timmers, P. R. et al. Genomics of 1 million parent lifespans implicates novel pathways and common diseases and distinguishes survival chances. Elife 8, 1–71 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/elife.39324&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

24. 24.Cortes, A., Albers, P. K., Dendrou, C. A., Fugger, L. & McVean, G. Identifying cross-disease components of genetic risk across hospital data in the UK Biobank. Nat. Genet. (2019). doi:10.1038/s41588-019-0550-4
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0550-4&link_type=DOI) 

25. 25.Jones, G. et al. Sarcopenia and variation in the Human Leukocyte Antigen complex. J. Gerontol. A. Biol. Sci. Med. Sci. (2019). doi:10.1093/gerona/glz042
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/gerona/glz042&link_type=DOI) 

26. 26.Chen, Y.-W. et al. Duration of chronic inflammation alters gene expression in muscle from untreated girls with juvenile dermatomyositis. BMC Immunol. 9, 43 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2172-9-43&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18671865&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

27. 27.Melzer, D., Pilling, L. C. & Ferrucci, L. The genetics of human ageing. Nat. Rev. Genet. (2019). doi:10.1038/s41576-019-0183-6
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41576-019-0183-6&link_type=DOI) 

28. 28.Murphy, L. et al. Lifetime risk of symptomatic knee osteoarthritis. Arthritis Care Res. 59, 1207–1213 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/art.24021&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000259669700002&link_type=ISI) 

29. 29.Francis-West, P. H. et al. Mechanisms of GDF-5 action during skeletal development. Development 126, 1305–15 (1999).
    
    [Abstract](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGV2ZWxvcCI7czo1OiJyZXNpZCI7czoxMDoiMTI2LzYvMTMwNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA1LzIxLzIwMjAuMDUuMTQuMjAxMDAzNTQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

30. 30.Capellini, T. D. et al. Ancient selection for derived alleles at a GDF5 enhancer influencing human growth and osteoarthritis risk. Nat. Genet. 49, 1202–1210 (2017).
    
    

31. 31.Miyamoto, Y. et al. A functional polymorphism in the 5’ UTR of GDF5 is associated with susceptibility to osteoarthritis. Nat. Genet. 39, 529–33 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/2005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17384641&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000245271200024&link_type=ISI) 

32. 32.Sanna, S. et al. Common variants in the GDF5-UQCC region are associated with variation in human height. Nat. Genet. 40, 198–203 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.74&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18193045&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000252732900018&link_type=ISI) 

33. 33.Uhalte, E. C., Wilkinson, J. M., Southam, L. & Zeggini, E. Pathways to understanding the genomic aetiology of osteoarthritis. Hum. Mol. Genet. 26, R193–R201 (2017).
    
    

34. 34.Pregizer, S. K. et al. Impact of broad regulatory regions on Gdf5 expression and function in knee development and susceptibility to osteoarthritis. Ann. Rheum. Dis. 77, 450 (2018).
    
    

35. 35.Traoré, M. et al. An embryonic CaVβ1 isoform promotes muscle mass maintenance via GDF5 signaling in adult mouse. Sci. Transl. Med. 11, (2019).
    
    

36. 36.Osipovich, A. B., Jennings, J. L., Lin, Q., Link, A. J. & Ruley, H. E. Dyggve-Melchior-Clausen syndrome: chondrodysplasia resulting from defects in intracellular vesicle traffic. Proc. Natl. Acad. Sci. U. S. A. 105, 16171–6 (2008).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA1LzQyLzE2MTcxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMjEvMjAyMC4wNS4xNC4yMDEwMDM1NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

37. 37.Paupe, V. et al. Recent advances in Dyggve-Melchior-Clausen syndrome. Mol. Genet. Metab. 83, 51–9
    
    

38. 38.Kim, J.-H. et al. Regulation of the catabolic cascade in osteoarthritis by the zinc-ZIP8-MTF1 axis. Cell 156, 730–43 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2014.01.007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24529376&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

39. 39.Tachmazidou, I. et al. Identification of new therapeutic targets for osteoarthritis through genome-wide analyses of UK Biobank data. Nat. Genet. 51, 230–236 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0327-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30664745&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

40. 40.den Hollander, W. et al. Genome-wide association and functional studies identify a role for matrix Gla protein in osteoarthritis of the hand. Ann. Rheum. Dis. 76, 2046–2053 (2017).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6ImFubnJoZXVtZGlzIjtzOjU6InJlc2lkIjtzOjEwOiI3Ni8xMi8yMDQ2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMjEvMjAyMC4wNS4xNC4yMDEwMDM1NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

41. 41.Herrmann, S. M. et al. Polymorphisms of the human matrix Gla protein (MGP) gene, vascular calcification, and myocardial infarction. Arterioscler. Thromb. Vasc. Biol. 20, 2386–2393 (2000).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYXR2YmFoYSI7czo1OiJyZXNpZCI7czoxMDoiMjAvMTEvMjM4NiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA1LzIxLzIwMjAuMDUuMTQuMjAxMDAzNTQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

42. 42.Rodríguez, A. J. et al. Aortic Calcification is Associated with Five-Year Decline in Handgrip Strength in Older Women. Calcif. Tissue Int. 103, 589–598 (2018).
    
    

43. 43.Ahmad, S., Jan, A. T., Baig, M. H., Lee, E. J. & Choi, I. Matrix gla protein: An extracellular matrix protein regulates myostatin expression in the muscle developmental program. Life Sci. 172, 55–63 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.lfs.2016.12.011&link_type=DOI) 

44. 44.Cui, G. et al. Association of Common Variants in *TGFA* with Increased Risk of Knee Osteoarthritis Susceptibility. Genet. Test. Mol. Biomarkers 21, 586–591 (2017).
    
    

45. 45.Styrkarsdottir, U. et al. Severe osteoarthritis of the hand associates with common variants within the ALDH1A2 gene and with rare variants at 1p31. Nat. Genet. 46, 498–502 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2957&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24728293&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

46. 46.Vermot, J., Niederreither, K., Garnier, J.-M., Chambon, P. & Dollé, P. Decreased embryonic retinoic acid synthesis results in a DiGeorge syndrome phenotype in newborn mice. Proc. Natl. Acad. Sci. U. S. A. 100, 1763–8 (2003).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMDoiMTAwLzQvMTc2MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA1LzIxLzIwMjAuMDUuMTQuMjAxMDAzNTQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

47. 47.Nakayama, N. et al. A novel chordin-like BMP inhibitor, CHL2, expressed preferentially in chondrocytes of developing cartilage and osteoarthritic joint cartilage. Development 131, 229–240 (2004).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGV2ZWxvcCI7czo1OiJyZXNpZCI7czo5OiIxMzEvMS8yMjkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNS8yMS8yMDIwLjA1LjE0LjIwMTAwMzU0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

48. 48.Chou, C.-H. et al. Insights into osteoarthritis progression revealed by analyses of both knee tibiofemoral compartments. Osteoarthr. Cartil. 23, 571–580 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.joca.2014.12.020&link_type=DOI) 

49. 49.Miotto, B. et al. The RBBP6/ZBTB38/MCM10 Axis Regulates DNA Replication and Common Fragile Site Stability. Cell Rep. 7, 575–587 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.celrep.2014.03.030&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24726359&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

50. 50.Oikawa, Y. et al. The methyl-CpG-binding protein CIBZ suppresses myogenic differentiation by directly inhibiting myogenin expression. Cell Res. 21, 1578–90 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/cr.2011.90&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21625269&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000296706700009&link_type=ISI) 

51. 51.UniProt. MT-CYB – Cytochrome b – Homo sapiens (Human) – MT-CYB gene & protein. Available at: [https://www.uniprot.org/uniprot/P00156](https://www.uniprot.org/uniprot/P00156). (Accessed: 20th April 2020)
    
    

52. 52.Yen, K., Lee, C., Mehta, H. & Cohen, P. The emerging role of the mitochondrial-derived peptide humanin in stress resistance. J. Mol. Endocrinol. 50, R11–9 (2013).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoiam1lIjtzOjU6InJlc2lkIjtzOjg6IjUwLzEvUjExIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMjEvMjAyMC4wNS4xNC4yMDEwMDM1NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

53. 53.Day, F. R. et al. Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk. Nat. Genet. 49, 834–841 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3841.&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

54. 54.Kuh, D., Hardy, R., Blodgett, J. M. & Cooper, R. Developmental factors associated with decline in grip strength from midlife to old age: A British birth cohort study. BMJ Open 9, 1–12 (2019).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYm1qb3BlbiI7czo1OiJyZXNpZCI7czoxMToiOS9TdXBwbF8zLzEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNS8yMS8yMDIwLjA1LjE0LjIwMTAwMzU0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

55. 55.Liaw, F. Y. et al. Higher Platelet-to-Lymphocyte Ratio Increased the Risk of Sarcopenia in the Community-Dwelling Older Adults. Sci. Rep. 7, (2017).
    
    

56. 56.Willer, C. J., Li, Y. & Abecasis, G. R. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics 26, 2190–1 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq340&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20616382&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000281738900017&link_type=ISI) 

57. 57.Watanabe, K., Taskesen, E., van Bochoven, A. & Posthuma, D. Functional mapping and annotation of genetic associations with FUMA. Nat. Commun. 8, 1826 (2017).
    
    

58. 58.de Leeuw, C. A., Mooij, J. M., Heskes, T. & Posthuma, D. MAGMA: Generalized Gene-Set Analysis of GWAS Data. PLoS Comput. Biol. 11, (2015).
    
    

59. 59.Jansen, I. E. et al. Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer’s disease risk. Nat. Genet. 2585331 (2019). doi:10.1038/s41588-018-0311-9
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0311-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30617256&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

60. 60.Michailidou, K. et al. Association analysis identifies 65 new breast cancer risk loci. Nature 551, 92–94 (2017).
    
    

61. 61.Pattaro, C. et al. Genetic associations at 53 loci highlight cell types and biological pathways relevant for kidney function. Nat. Commun. 7, 10023 (2016).
    
    

62. 62.van der Harst, P. & Verweij, N. Identification of 64 Novel Genetic Loci Provides an Expanded View on the Genetic Architecture of Coronary Artery Disease. Circ. Res. 122, 433–443 (2018).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTA6ImNpcmNyZXNhaGEiO3M6NToicmVzaWQiO3M6OToiMTIyLzMvNDMzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMjEvMjAyMC4wNS4xNC4yMDEwMDM1NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

63. 63.Trajanoska, K. et al. Assessment of the genetic and clinical determinants of fracture risk: genome wide association and mendelian randomisation study. BMJ 362, k3225 (2018).
    
    

64. 64.Zengini, E. et al. Genome-wide analyses using UK Biobank data provide insights into the genetic architecture of osteoarthritis. Nat. Genet. 50, 549–558 (2018).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

65. 65.Schumacher, F. R. et al. Association analyses of more than 140,000 men identify 63 new prostate cancer susceptibility loci. Nat. Genet. 50, 928–936 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0142-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

66. 66.Okada, Y. et al. Genetics of rheumatoid arthritis contributes to biology and drug discovery. Nature 506, 376–381 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature12873&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24390342&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000331477800043&link_type=ISI) 

67. 67.Malik, R. et al. Multiancestry genome-wide association study of 520,000 subjects identifies 32 loci associated with stroke and stroke subtypes. Nat. Genet. 50, 524–537 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0058-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29531354&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

68. 68.Mahajan, A. et al. Fine-mapping type 2 diabetes loci to single-variant resolution using high-density imputation and islet-specific epigenome maps. Nat. Genet. (2018). doi:10.1038/s41588-018-0241-6
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0241-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30297969&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

69. 69.Yengo, L. et al. Meta-analysis of genome-wide association studies for height and body mass index in ∼700000 individuals of European ancestry. Hum. Mol. Genet. 27, 3641–3649 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddy271&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30124842&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom) 

70. 70.Pulit, S. L. et al. Meta-analysis of genome-wide association studies for body fat distribution in 694 649 individuals of European ancestry. Hum. Mol. Genet. 28, 166–174 (2019).
    
    

71. 71.Zillikens, M. C. et al. Large meta-analysis of genome-wide association studies identifies five loci for lean body mass. Nat. Commun. 8, 80 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-017-00031-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28724990&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F21%2F2020.05.14.20100354.atom)