Genome-wide association studies of 27 accelerometry-derived physical activity measurements identifies novel loci and genetic mechanisms
=======================================================================================================================================

* Guanghao Qi
* Diptavo Dutta
* Andrew Leroux
* Debashree Ray
* Ciprian Crainiceanu
* Nilanjan Chatterjee

## Abstract

Physical activity (PA) is an important risk factor for a wide range of diseases. Previous genome-wide association studies (GWAS), based on self-reported data or a small number of phenotypes derived from accelerometry, have identified a limited number of genetic loci associated with habitual PA and provided evidence for involvement of central nervous system in mediating genetic effects. In this study, we derived 27 PA phenotypes from wrist accelerometry data obtained from 93,745 UK Biobank study participants. Single-variant association analysis based on mixed-effects models and transcriptome-wide association studies (TWAS) together identified 6 novel loci that were not detected by previous studies. For both novel and previously known loci, we discovered associations with novel phenotypes including active-to-sedentary transition probability, light-intensity PA, activity during different times of the day and proxy phenotypes to sleep and circadian patterns. Follow-up studies indicated the role of the blood and immune system in modulating the genetic effects and a secondary role of the digestive and endocrine systems.

## Introduction

Regular physical activity (PA) is associated with lower risk of a wide range of diseases, including cancer, diabetes, cardiovascular disease1, Alzheimer’s disease2, as well as mortality3,4. However, studies have indicated that large majority of US adults and adolescents are insufficiently active5, and thus PA interventions have great potential to improve public health. PA was shown to have a substantial genetic component, and understanding its genetic mechanism can inform the design of individualized interventions6,7. For example, people who are genetically pre-disposed to low PA may benefit more from early and more frequent guidance.

A number of previous genome-wide association studies (GWAS) on physical activity have relied on self-reported phenotypes, which are subject to perception and recall error8–11. Recently, wearable devices have been used extensively to collect physical activity data objectively and continuously for multiple days. To date, there have been two GWAS based on acceleromtery-derived activity phenotypes. Both studies used data from the UK Biobank study12,13 but only focused on a few summaries of these high-density PA measurements. One study considered two accelerometry-derived phenotypes (average acceleration and fraction accelerations > 425 milli-gravities) and identified 3 loci associated with PA11. A second study used a machine learning approach to extract PA phenotypes, including overall activity, sleep duration, sedentary time, walking and moderate intensity activity14. This study identified 14 loci associated with PA and found that the central nervous system (CNS) plays an essential role in modulating the genetic effects on PA. However, both studies used a small number of phenotypes, which may not capture the complexity of PA patterns.

Recent studies suggest that in addition to the total volume of activity, other PA summaries may be strongly associated with human health and mortality risk. For example, the transition between active and sedentary states was strongly associated with measures of health and mortality4,15. PA relative amplitude, a proxy for sleep quality and circadian rhythm, was strongly associated with mental health16. Moderate-to-vigorous PA (MVPA) and light intensity PA (LIPA) have also been reported to be associated with health17,18. Thus, there is increasing evidence that objectively measured PA in the free-living environment is a highly complex phenotype that requires a large number of summaries that provide complementary information. Understanding the genetic mechanisms behind these summaries is critical for understanding the genetic regulation of activity behavior and informing targeted interventions.

In this paper, we conducted genome-wide association analysis using 27 accelerometry-derived PA measurements from UK Biobank data12,13. The phenotypes cover a wide range of features including volumes of activity, activity during different times of the day, active to sedentary transition probabilities and various principal components (**Supplementary Table 1**). We conducted GWAS using a mixed-model-based method, fastGWA19, to identify variants associated with the above phenotypes. We also conducted transcriptome-wide association studies (TWAS)20 across 48 tissues to identify genes and tissues harboring the associations. We further conducted tissue-specific heritability enrichment21,22, gene-set enrichment23 and genetic correlation24 analyses to further reveal the underlying biological mechanisms. We identified 6 novel loci associated with PA and showed that, in addition to the CNS, blood and immune related mechanisms could play an important role in modulating the genetic effects on activity, and digestive and endocrine tissues could play a secondary role.

## Results

### Genetic Loci Associated with Physical Activity

Single-variant genome-wide association analysis identified a total of 16 independent loci, including five novel ones compared to previous studies (**Table 1 and Fig. 1**). The locus indexed by single-nucleotide polymorphism (SNP) rs301799 on chromosome 1 was associated with total log acceleration (TLA) between 6pm to 8pm, representing early evening activity. Three novel loci were discovered on chromosome 3: the locus indexed by rs3836464 was associated with active-to-sedentary transition probability (ASTP); the locus indexed by rs9818758 was associated with relative amplitude, which is a proxy sleep behavior and circadian rhythm25; the locus indexed by insertion-deletion (INDEL) 3:131647162\_TA_T (no rsid available) was associated with TLA 2am-4am which is a proxy phenotype for activity during sleep. LIPA appeared to be associated with other SNPs near 3:131647162_TA_T but not the lead variant itself, indicating multiple independent signals at the same locus (**Fig. 1**). Another locus indexed by rs2138543 is associated with TLA 6am-8am which represents early morning activity, and the second principal component of log-acceleration (PC2), for which the interpretation is less clear (**Table 1**).

View this table:
[Table 1.](http://medrxiv.org/content/early/2021/02/18/2021.02.15.21251499/T1)

Table 1. Significant loci associated with physical activity in single-variant analysis.

![Fig. 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.15.21251499/F1.medium.gif)

[Fig. 1](http://medrxiv.org/content/early/2021/02/18/2021.02.15.21251499/F1)

Fig. 1 Manhattan plot for 18 traits that are significantly associated with at least one variant at *p* < 2. 63 ×10−9 in single-variant analysis.
The red dashed line is *p* = 2.63 ×10−9 accounting for the number of independent traits. The blue dashed line is the standard genome-wide significance threshold *p* = 5 ×10−8. Five novel loci that have not been discovered in previous GWAS of physical activity, sleep duration and circadian rhythm are circled out (see Table 1).

Our analysis also identified novel phenotypes for several known loci (**Table 1**). The strongest signal was seen for the locus indexed by rs113851554 which is associated with multiple sleep and circadian rhythm proxy phenotypes including TLA 12am-2am (*p* = 6.7 ×10−37), TLA 2am-4am (*p* = 7.9 ×10−39), average log acceleration during the least active 5 hours of the day (L5, *p* = 1.3 ×10−33), timing of L5 (*p* = 5.4 ×10−22) and PA relative amplitude (*p* = 6.9 ×10−15). This locus was previously identified to be associated with accelerometry-derived sleep duration in UK Biobank14. Among other known loci, 5 were only discovered in the GWAS of self-reported circadian rhythm26 but not in the other studies considered (**Table 1**, last column). In our analysis, the loci indexed by rs1144566, rs9369062 and rs12927162 were associated with sleep proxy phenotypes including timing of L5, TLA 12am-2am, TLA 2am-4am and TLA 10pm-12am. Two other loci, indexed by rs2909950 and rs12717867, were associated with TLA 6pm-8pm and LIPA, respectively.

### Transcriptome-Wide Association Study and Colocalization Analysis

We performed transcriptome-wide association studies (TWAS) 20,27 for each PA trait based on gene expression data across 48 tissues available through GTEx (version 7)28. Our analysis identified 15 loci (**Table 2, Fig. 2, Supplementary Table 2**) with significant association in at least one trait-tissue pair analysis after correcting for multiple testing (Benjamini-Hochberg corrected p-value < 2.5 ×10−6, see **Methods**). We identified a novel locus and an underlying index pseudogene *PDXDC2P (16q22*.*1)*, the expression of which in esophagus mucosa and EBV transformed lymphocytes appeared to be genetically associated with TLA 6am-8am (**Fig. 2**). The locus was not previously reported by any prior GWAS and were not close to any of the 5 novel regions detected by our single variants analysis (**Table 2**).

View this table:
[Table 2.](http://medrxiv.org/content/early/2021/02/18/2021.02.15.21251499/T2)

Table 2. Significant loci associated with physical activity identified in TWAS.

![Fig. 2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.15.21251499/F2.medium.gif)

[Fig. 2](http://medrxiv.org/content/early/2021/02/18/2021.02.15.21251499/F2)

Fig. 2 TWAS Manhattan plots for two tissues that harbor the novel TWAS locus represented by *PDXDC2P*.
The figure shows -log10 p-value for all the genes expressed in esophagus mucosa (digestive, upper panel) or EBV-transformed lymphocyte cells (blood/Immune, lower panel) as reported in GTEx v7. Significant genes with FDR corrected p-values< 2.5 ×10−6 are circled (see Methods for details). Only the PA traits that are significantly associated with at least one variant at *p* < 2.63 ×10−9 in single-variant analysis are shown.

The TWAS analysis also identified novel PA phenotypes, potential target genes and underlying tissues for many of the known loci or novel loci detected through single variant analysis (**Table 2**). Consistent with a previous study14, the TWAS analysis showed that genetic association for PA traits often points towards involvement of CNS (**Table 2**). Further, our analysis indicates consistent involvement of blood and immune, digestive and endocrine systems in modulating the genetic effects on PA. Among the 15 loci significant in TWAS analysis, the lead genes of 4 loci were significantly associated with PA phenotypes via the blood and immune tissues. For example, the genetically predicted expression of *PBX3* and *KANSL1* in the blood and immune tissues were each associated with 3 PA phenotypes. The genes associated with PA via blood and immune tissues were also associated via digestive (e.g., esophagus mucosa) and/or endocrine (e.g., thyroid) tissues, but only one of them overlapped with the 9 genes that were associated with PA phenotypes through the CNS (**Table 2**). Another locus, represented by *C3orf62*, were associated with PA relative amplitude only via the digestive and endocrine tissues but not the blood/immune tissues or CNS. These findings suggested that the genetic regulation of PA occurs via at least two different pathways: a primary pathway involving the CNS (brain in particular) and a secondary pathway involving the blood/immune system and, potentially, the digestive and endocrine systems.

Several genes that were found to be significantly associated to specific PA traits in our TWAS analysis, were also found to be highly overlapping with genes that were previously reported to be associated with various traits and diseases including but not limited to neuropsychiatric diseases, behavioral traits, anthropometric traits and autoimmune diseases (**Supplementary Fig. 1**). For example, we found that the genes associated with TLA across different tissues, are enriched for genes that have been associated with neuroticism, bipolar disorder, Parkinson’s disease, cognitive function and several others indicating the putative involvement of the CNS in the genetic mechanism of TLA. Additionally, the genes associated with relative amplitude overlapped highly with those associated with several autoimmune diseases like inflammatory bowel disease, ulcerative colitis in addition to different behavioral and cognitive traits (**Supplementary Fig. 1**). These results further supported the possible involvement of both CNS as well as the blood and immune system in the genetic mechanism of PA traits.

We performed a colocalization analysis to gain further insights on the tissue specific activity of the significant genetic loci. Among the 16 loci significantly associated with PA, 9 loci colocalized with the eQTL signals for at least one gene and one tissue with a colocalization probability (PP4) > 0.8 (**Supplementary Table 3**). Colocalization occurred in a similar set of tissues as those that harbored the TWAS associations (**Table 2 and Supplementary Table 3**), namely the CNS, blood and immune, digestive and endocrine tissues, and also in a number of cardiovascular tissues that were not highlighted by TWAS. Among the 15 lead genes for TWAS significant loci, the eQTL signal of 4 genes (*RERE, C3orf62, PBX3* and *RP11-396F22*.*1*) colocalized with PA GWAS signal in at least one tissue. Colocalization also occurred in two other secondary genes (*RP5-1115A15*.*1* and *CASC10*).

### Analysis of Heritability and Co-Heritability

Our fastGWA analysis estimated genome-wide heritability of PA phenotypes as an intermediate output. The estimates appeared to be dependent on the sparsity level of the genetic relationship matrix (**Supplementary Fig. 2**). We chose the results under the lower cutoff (0.02) since it captured more subtle relatedness and should give more accurate heritability estimates. The estimates of heritability varied across different PA phenotypes. A number of traits were estimated to have higher heritability than others, including TLA (0.15), TLA 6pm-8pm (0.15), MVPA (0.14) (**Supplementary Fig. 2**). Afternoon and pre-sleep evening activity (TLA 4pm to 12am) appeared to be more heritable than morning activity (TLA 2am to 12pm). As could be expected, phenotypes with higher heritability tend to have a higher average *X*2 statistic for genetic associations, and a QQ plot which deviate further from the null line (**Supplementary Fig. 3**).

We further used stratified LD-score regression for partitioning heritability by functional annotations of genome21,22. Consistent with TWAS findings, this analysis also indicated possible role for blood and immune system in addition to CNS for genetic regulation of PA (**Fig. 3**). In particular, heritabilities for both TLA and LIPA were enriched for DNase I hypersensitivity sites (DHS) in primary B cells from peripheral blood and that for TLA 12pm-2pm were enriched for H3K27ac in spleen. We also found potential enrichment in other traits, though they were not significant after FDR adjustment. For example, for TLA 8am-10am, MVPA, and ASTP the heritability enrichment in active chromatin regions of blood/immune tissues were all close to being statistically significant (**Supplementary Fig. 4**).

![Fig. 3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.15.21251499/F3.medium.gif)

[Fig. 3](http://medrxiv.org/content/early/2021/02/18/2021.02.15.21251499/F3)

Fig. 3 Tissue-specific heritability enrichment p-values for traits with significant enrichment at FDR < 0.05 in blood and immune tissues.
The analysis was conducted using tissue/cell type specific stratified LD score regression based on 6 chromatin-based annotations in 111 tissues and cell types described in Finucane et al, Nature Genetics 2018 (PMID: 29632380). Each dot corresponds to an annotation in a tissue or cell type. A complete list of tissue and cell types is provided in Supplementary Table 7 of the above paper. Black line corresponds to FDR < 0.05 (-log(p-value)=3.83) across all combinations of trait, tissue, and histone mark. Red line corresponds to p = 0.05. See Supplementary Fig. 4 for the enrichment p-values for the rest of the traits. CNS: central nervous system.

We further used LD score regression24,29 to explore genetic correlation between PA phenotypes and four broad groups of complex traits and diseases (**Supplementary Fig. 5 and Supplementary Table 4**). Genetic correlations were identified (FDR < 10%) between PA phenotypes and: (1) neurological, psychiatric and cognitive traits, including Alzheimer’s disease (AD), attention-deficit hyperactivity disorder (ADHD), depressive symptoms, intelligence, and neo-conscientiousness; (2) auto-immune diseases, with the strongest correlation for multiple sclerosis and weaker correlations for Crohn’s disease and primary billary cirrhosis; (3) obesity-related anthropometric traits and (4) cholesterol levels. These results broadly supported our previous results indicating the role of CNS and blood/immune related mechanisms in the genetics of PA traits.

## Discussion

In summary, our study provided novel insights to genetic architecture of physical activity through genome-wide association analysis of an extensive set of accelerometry based PA phenotypes, derived in the UK biobank study, and a series of follow-up genomic analyses. We identified a total of six novel loci, most of which were associated with PA phenotypes not considered in previous studies 11,14,26,30. Our analysis also identified novel phenotypes associated with the known loci. Further, we provided multiple independent lines of evidence that genetic mechanisms for association for PA involve the blood and immune system, which can be potential targets for developing future interventions.

Compared to the 15 loci identified by the two previous GWASs on accelerometry-based PA 11,14, the novel loci we discovered have increased the number of PA susceptibility loci by 40%. Most of the novel loci were connected to the expression of genes, pseudogenes or long non-coding RNAs (lncRNA, **Table 1 and 2**). The novel locus index by rs301799 overlapped with the TWAS locus indexed by *RERE* (**Table 1 and 2**). The *RERE* gene was shown to be important for early development of brain, eyes, inner ear, heart and kidneys, which could have complex effects on an individual’s ability to perform physical activity 31. The novel locus indexed by rs9818758 overlaps with the TWAS locus index by *C3orf62*. Though it was unclear how *C3orf62* is involved in PA, two secondary genes in the locus, *ARIH2* and *DAG1* (**Table 2**), appeared to be involved in relevant biological processes. *ARIH2* was found to be essential for embryogenesis by regulating the immune system32; *DAG1* was found to play a role in the regeneration of skeletal muscles33. Another two novel loci were connected to pseudogene *PDXDC2P* and lncRNA *RP11-396F22*.*1* of which the function is less clear and may be worth future lab investigation.

The novel phenotypes in this study provided important insights into the genetic architecture of PA, which may have been overlooked by previous GWASs on a small number of phenotypes. The accelerometry-based study by Doherty et al identified the genetic associations with overall activity, sleep duration and sedentary time14; the study by Klimentidis et al studied the average acceleration and the duration of active states11. Our results found that there can be different genetic architecture for PA during different times of the day, and there can be unique variants that only affect certain PA patterns, like ASTP, LIPA and relative amplitude, but not others (**Table 2**). The heritability and genetic correlation can also vary across different PA phenotypes (**Supplementary Fig. 2 and Supplementary Fig. 5**).

TWAS and tissue-specific heritability enrichment analysis suggested that in addition to the CNS, the blood and immune system could be also associated with PA. This finding was further supported by colocalization, gene-set enrichment and genetic correlation analyses. A previous study14, which explored enrichment of heritability for PA traits by tissue-specific gene expression patterns, identified potential modulating role of the CNS, adrenal/pancreatic and skeletal muscle tissues. Our study, which used a more extended set of phenotypes and chromatin-state-based annotations, confirmed previous findings and further highlighted the role of the blood and immune system. Previous medical literature has established the effect of PA on immune functions. A study showed that higher PA is associated with elevation of T-regulatory cells and lower risk for autoimmune diseases34. Multiple studies showed that regular moderately intense PA boost immune functions in older adults and protects against age-related inflammatory disorders35–37. Our analysis suggested the link between PA and immune functions in the genetic pathways and future studies are needed to better understand the underlying mechanisms and causal directions.

In addition to the blood and immune system, TWAS and enrichment analysis also suggested that the digestive system and endocrine system could be involved in modulating the genetic effects on PA. A previous study found that PA has complex effects on gastroinstestinal health38: acute strenuous activity may provoke gastrointestinal symptoms while low-intensity activity could have benefits. Interestingly, three TWAS loci that were significant in digestive tissues were associated with PA phenotypes that are proxies for meal-time activity: *PDXDC2P* with TLA 6am-8am, *RERE* with TLA 6pm-8PM and *KANSL1* with TLA 4pm-6pm (**Table 2**). It was also known that multiple organs in the endocrine system produce hormones that regulate physiological functions of the body, which can have complex bidirectional relationships with PA39–43. Our TWAS analysis indicated that the genes associated with PA via the blood and immune system tended to also be associated with the digestive and endocrine systems, but do not usually overlap with the genes associated with the CNS. This suggests that the blood and immune, digestive and endocrine systems may be involved in the same broad pathway that affects PA, which is different from that of CNS.

This study has a number of limitations. Though we derived a more extensive set of PA phenotypes than previous studies, information was still lost when collapsing a 7-day continuous times series of wrist accelerometry into 31 PA phenotypes. The ideal approach would be to conduct a GWAS utilizing all the information across the 7 days of accelerometry measurements. Results could outline genetic regulation of a continuous course of PA over time. The current analysis of TLA during 12 non-overlapping two-hour time intervals during the day, indicated that different genetic variants may affect PA during different times of the day (**Tables 1 and 2**). Another limitation is that some of the phenotypes are not directly interpretable. For example, the PCs of log acceleration are less interpretable than other phenotypes, such as TLA and ASTP.

However, they do reflect important features of physical activity and warrant further investigations. A potential solution is to obtain proxy measurements that are interpretable and highly correlated with PC scores.

In conclusion, we conducted association studies on a wide range of PA phenotypes and identified 6 novel loci associated with PA. We found that in addition to the CNS, the blood and immune system may also play an important role in the genetic mechanisms of PA, and the digestive and endocrine systems could also be involved in the blood and immune pathway.

## Materials and Methods

### Study Cohort and Physical Activity Phenotypes

The UK Biobank study consists of ∼500,000 individuals in the United Kingdom with comprehensive genotype and phenotype data13. We used a subset of the 103,712 individuals who were invited and agreed to participate in the accelerometry sub-study where participants wore a wrist-worn accelerometer for up to 7 days12. Accelerometry data from participants are available at multiple resolutions. Here, the individual-specific set of accelerometry-based phenotypes was derived from the five-second level acceleration data provided by the UK Biobank team. Individuals were screened for poor quality data using indictors provided by the UK Biobank. In addition, we required individuals to have at least 3 days (12am-12am) of sufficient wear time defined as estimated wear time greater than 95% of the day (>= 1368 minutes). Our inclusion criteria for this analysis closely mirrors that described in a related paper from our group44 with the exception that we did not exclude participants younger than 50 at the time of accelerometer wear or based on missing demographic and lifestyle data and instead excluded individuals based on ancestry and genotype data (see subsection *Genotype Data* below).

Physical activity phenotypes were all calculated at the day level and then averaged within study participants across days to obtain one measure for each phenotype and study participant. This led to in 31 PA phenotypes for 93,745 study participants that covered a wide spectrum of information: 1) total volume of activity (total acceleration (TA), total log acceleration (TLA)); 2) activity during 12 disjoint two-hour windows of the day (TLA 12am-2am, TLA 2am-4am, …, TLA 10pm-12am); 3) duration of sedentary state (ST), LIPA and MVPA; 4) PA principal components (PC1-6); 5) active-to-sedentary transition probability (ASTP) and sedentary-to-active transition probability (SATP); 6) proxy phenotypes for circadian patterns, including dynamic activity ratio estimate (DARE), activity during the most active 10 hours (M10) and least active 5 hours (L5) of the day, timing of M10 and L5, and PA relative amplitude. They included most of the phenotypes used in the previous PA association studies as well (See **Supplementary Table 1** for details)11,14. The exact procedure for deriving study participant-specific phenotypes is described in detail in the supplemental material of the related paper from our group44. The phenotypes were inverse-normal transformed ![Formula][1]</img>  where the transformed variables have mean 0 and variance 1.

### Removing Highly Correlated Phenotypes

Some of the initial 31 PA phenotypes were highly correlated (**Supplementary Fig. 6**). To avoid counting similar phenotypes multiple times, if two phenotypes had correlation > 0.8 we removed one of them. First, we removed total acceleration (TA), duration of sedentary state (ST), PC1 and M10 due to their high correlation with total log acceleration (TLA). TLA was retained as the main metric for the total volume of activity. Most previous studies used TA instead of TLA as the main metric for the volume of activity. However, the distribution of TA is highly skewed, which could lead to lower power for association testing and instability of findings due to dependence on extreme observations. In total, 4 phenotypes were removed and 27 PA phenotypes were retained for the association analysis.

### Genotype Data

The imputed genotype data for ∼93 million variants, using UK10K, 1000 Genomes (Phase 3) and Haplotype Reference Consortium as reference panel, provided by UK Biobank were used and merged with the PA phenotype data. We excluded study participants according to the following criteria: 1) non-white ancestry; 2) putative sex chromosome aneuploidy; 3) an excessive number of relatives (more than 10 putative third-degree relatives in the kinship table); 4) sample was not in the input for phasing of chr1-chr22. After applying these exclusion criteria the sample was further reduced to 88,411 study particiants for downstream analysis.

We conducted variant quality control to ensure that genetic variants with poor genotyping quality do not affect the results. Specifically, variants that satisfy any of the following criteria were removed: 1) imputation INFO score < 0.8; 2) MAF < 0.01; 3) Hardy-Weinberg Equilibrium (HWE) p-value < 1 ×10−6; 4) missing in more than 10% study participants. After the filtering, 8,951,705 variants remained for downstream analysis, of which 8,067,228 (90.1%) were single nucleiotide polymorphisms (SNPs) and the rest (9.9%) were insertion-deletions (INDELs).

### Association Analysis

We used a fast mixed-effects model method, fastGWA19, for genome-wide association analysis. Like other mixed-effects model methods, fastGWA allows the inclusion of related and unrelated individuals but improves computational efficiency by incorporating a sparse genetic relationship matrix (GRM). The GRM measures the genetic similarity between individuals and each element is the correlation of genotypes between a pair of individuals. We constructed the GRM using LD-pruned variants that had MAF > 5% and were present in HapMap3 (LD-pruning was done in PLINK using the following set up as recommended in Jiang et al 19: window size = 1000Kb, step-size = 100 and r2 = 0.9). We further computed a sparse-GRM at sparsity level 0.05 to capture the genetic relatedness between the closely related individuals only and reduced others to zero. We used the Haseman-Elston regression to estimate the variance of the random effects as an intermediate step of fastGWA. This approach is orders of magnitude faster than the previous state-of-the-art, BOLT-LMM45,46.

Models were adjusted for age, sex and the first 20 genetic principal components as covariates. Because the PA phenotypes are correlated, principal component analysis (PCA) was conducted on the phenotypes to estimate the number of independent phenotypes before setting the GWAS significance threshold. At least 19 phenotype PCs were needed to explain 99% percent of the PA phenotypic variance (**Supplementary Fig. 7**). Variants with p-value below the threshold 5 ×10−8/19 = 2.63 ×10−9 were declared to be statistically significant, which accounted for the number of independent phenotypes. LD clumping was conducted based on the minimum p-value across phenotypes. The requirements for the lead SNPs of different loci were to have r2 < 0.1 and be at least > 500kb apart.

A locus was defined as novel if its lead variant is > 500kb from the lead variant of any known loci discovered by the following GWASs on PA, sleep, and circadian rhythm: (1) Doherty et al study on a smaller set of accelerometry-derived PA phenotypes14; (2) Klimentidis et al study on self-reported and accelerometry-derived PA11; (3) Dashti et al study on self-reported sleep duration30; (4) Jones et al study on circadian rhythm26.

Transcriptome-wide association studies (TWAS) were conducted using the FUSION R program20 with reference models generated from 48 tissues of GTEx v728. TWAS analysis was limited to the 18 traits with at least one genome-wide significant variant (*p* < 2.63 ×10−9). Multiple testing due to the large number of tissue-trait combinations (48*18=864) was addressed by a two-stage adjustment approach: 1) for each variant, the Benjamini-Hochberg (BH) adjustment was applied across all tissue-trait pairs; 2) each variant with BH-adjusted p-value 2.5 × 10−6 was then identified (accounting for 20,000 protein-coding genes). Since there can be multiple genes in close proximity to each other, to identify independent loci detected by TWAS analysis, genes were clustered based on significant associations. A clumping approach was used, which selected the gene with the smallest minimum p-value across tissue-trait pairs and removed the other genes with a transcription start site (TSS) within 1Mb of the lead gene TSS. The process continued by identifying the gene with the next smallest minimum p-value and iterating. The only exception was when the lead gene of the cluster was not a protein-coding gene (e.g., pseudogene, lncRNA) and a protein-coding gene was in the cluster. In this case the protein-coding gene with the smallest minimum p-value was identified as the lead gene. This led to independent gene clusters at genomic loci which were least 1Mb apart, i.e., none of the lead gene TSS is within the cis region of another lead gene.

### Enrichment Analysis

Stratified LD score regression21,22 was used to identify the tissues and genomic annotations enriched by the heritability for PA. For tissue specific analysis, chromatin-based annotations were used as derived from the ENCODE and Roadmap data47,48 by Finucane et al21. The annotations were based on narrow peaks of DNase I hypersensitivity site (DHS) and five activating histone marks (H3K27ac, H3K4me3, H3K4me1, H3K9ac and H3K36me3) observed for 111 tissues or cell types, resulting in a total of 489 annotations. Stratified LD score regression computes the heritability attributed to each annotation and computes a coefficient and a p-value that characterize enrichment.

In a separate analysis, the enrichment of TWAS signals was evaluated among the genes that have been reported to be associated to different traits, using FUMA23,49. For a given PA trait, we defined a gene-set as the genes that were significant at an exome-wide level (*p* < 2.5 ×10−6) and investigated whether these genes overlapped with the genes that have been mapped to genome-wide significant variants for different traits as reported in GWAS catalog50. The collection of such genes have been detailed in Molecular Signatures Database (MSigDB)51. We used FUMA to compute the proportion of genes related to other diseases and traits that were also identified by our TWAS analysis and computed enrichment p-values using the Fisher’s exact test.

### Colocalization Analysis

For each susceptibility locus of PA (**Table 1**), colocalization analysis was conducted between its most significantly associated phenotype and eQTL effects on gene expression in 48 tissues in GTEx v728. SNPs within +-200kb radius of the lead SNP were used and genes that had at least one significant eQTL (q-value < 0.05) in the region were considered. Analysis was conducted using the R package COLOC 52 and GWAS and eQTL effects were identified as being colocalized if PP4 > 0.8.

### Heritability and Genetic Correlation Analysis

Heritability of activity phenotypes was estimated using Haseman-Elston regression as an intermediate output of fastGWA19. Our fastGWA analysis computed sparse GRM at sparsity level 0.05 as recommended by the fastGWA paper (see “Association analysis”). However, this cutoff may miss the subtle relatedness in the sample and affect heritability estimate. As a sensitivity analysis, we re-estimated the heritability using a lower sparsity threshold at 0.02 to capture more subtle relatedness. The genetic correlation between 18 PA traits and 238 complex traits and diseases was estimated using LD score regression24 implemented in LD Hub53. In particular, we focused on four broad groups of traits and diseases (A) cholesterol levels (B) anthropometric traits (C) autoimmune disease and (D) miscellaneous traits including psychiatric, neurological, cognitive and personality traits. For each trait and within each category, we applied a false discovery rate correction to the p-values corresponding to the genetic correlation estimated using LD score regression, to account for multiple testing. Any genetic correlation with FDR-adjusted p-value less that 10% were declared as significant.

## Supporting information

Supplementary Tables [[supplements/251499_file03.xlsx]](pending:yes)

## Data Availability

Data supporting the findings of this paper are available upon application to the UK Biobank study (https://www.ukbiobank.ac.uk/). The summary statistics will be made publicly available upon publication.

## Data Availability

Data supporting the findings of this paper are available upon application to the UK Biobank study ([https://www.ukbiobank.ac.uk/](https://www.ukbiobank.ac.uk/)). The summary statistics will be made publicly available upon publication.

## Competing Interests

Dr. Ciprian Crainiceanu is consulting with Bayer and Johnson and Johnson on methods development for wearable devices in clinical trials. The details of the contracts are disclosed through the Johns Hopkins University eDisclose system and have no direct or apparent relationship with the current paper. The other authors declare no conflict of interest.

## Ethics Statement

Approval was granted by the UK Biobank study under application ID 17712 to use the data in the present work. UK Biobank has ethical approval from the North West Multi-Centre Research Ethics Committee (approval number 16/NW/0274). All participants provided informed consent to participate.

## Web Resources

UK Biobank: [https://www.ukbiobank.ac.uk/](https://www.ukbiobank.ac.uk/)

fastGWA software: [https://cnsgenomics.com/software/gcta/#fastGWA](https://cnsgenomics.com/software/gcta/#fastGWA)

FUSION TWAS software: [http://gusevlab.org/projects/fusion/](http://gusevlab.org/projects/fusion/)

COLOC R package: [https://cran.r-project.org/web/packages/coloc/index.html](https://cran.r-project.org/web/packages/coloc/index.html)

LD score regression software: [https://github.com/bulik/ldsc](https://github.com/bulik/ldsc)

LD Hub: [http://ldsc.broadinstitute.org/ldhub/](http://ldsc.broadinstitute.org/ldhub/)

FUMA GWAS: [https://fuma.ctglab.nl/](https://fuma.ctglab.nl/)

PLINK: [https://www.cog-genomics.org/plink/2.0/](https://www.cog-genomics.org/plink/2.0/)

Molecular Signatures Database: [http://www.broadinstitute.org/msigdb](http://www.broadinstitute.org/msigdb)

## Supplementary Figures

![Supplementary Fig. 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.15.21251499/F4.medium.gif)

[Supplementary Fig. 1](http://medrxiv.org/content/early/2021/02/18/2021.02.15.21251499/F4)

Supplementary Fig. 1 Gene-set enrichment analysis for the significant genes in TWAS analysis for (A) TLA and (B) relative amplitude.
We use TLA and relative amplitude as examples because they capture two most important patterns of PA: total volume of activity and sleep quality. Genes previously reported to be associated to different traits and diseases as curated in the GWAS catalog were defined as gene-sets (see Methods). The genes associated to (A) TLA and (B) relative amplitude (p-value < 2.5 ×10−6) in transcriptome-wide analysis were checked for enriched overlap with the gene-sets curated from GWAS catalog for each trait (additionally reported in Molecular Signatures database).

![Supplementary Fig. 2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.15.21251499/F5.medium.gif)

[Supplementary Fig. 2](http://medrxiv.org/content/early/2021/02/18/2021.02.15.21251499/F5)

Supplementary Fig. 2 Estimates of heritability of 27 physical activity traits using Haseman-Elston regression and sparse genetic relationship matrix (GRM).
a) Cutoff at 0.05: correlations that are < 0.05 in the GRM are reduced to zero, as is in our fastGWA analysis and recommended by the fastGWA paper. b) Cutoff at 0.02: correlations that are < 0.02 in the GRM are reduced to zero.

![Supplementary Fig. 3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.15.21251499/F6.medium.gif)

[Supplementary Fig. 3](http://medrxiv.org/content/early/2021/02/18/2021.02.15.21251499/F6)

Supplementary Fig. 3 QQ plot for the p-values from fastGWA analysis for 27 activity phenotypes.
The average *χ*2 statistic is annotated at the top left corner of each panel.

![Supplementary Fig. 4](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.15.21251499/F7.medium.gif)

[Supplementary Fig. 4](http://medrxiv.org/content/early/2021/02/18/2021.02.15.21251499/F7)

Supplementary Fig. 4 Tissue-specific heritability enrichment for traits not included in Figure 3.
The analysis was conducted using tissue/cell type specific stratified LD score regression (Finucane et al, Nat Genet 2018, PMID 29632380) based on 6 chromatin-based annotations in 111 tissues and cell types. A complete list of tissue and cell types is provided in Supplementary Table 7 of the above paper. Black line corresponds to FDR < 0.05 (-log(p-value)=3.83) across all combinations of trait, tissue, and histone mark. Red line corresponds to p = 0.05. The analysis was conducted using tissue/cell type specific stratified LD score regression (Finucane et al, 2018).

![Supplementary Fig. 5](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.15.21251499/F8.medium.gif)

[Supplementary Fig. 5](http://medrxiv.org/content/early/2021/02/18/2021.02.15.21251499/F8)

Supplementary Fig. 5 Genetic correlation of 18 PA traits with other traits: (A) cholesterol levels (B) anthropometric traits (C) autoimmune diseases and (D) miscellaneous traits including psychiatric, neurological, cognitive and personality traits.
The PubMed ID’s for the corresponding for the GWAS summary statistics are in parentheses. The size of the circle and the darkness of the color represents the genetic correlation value. The asterisk sign represents if the corresponding genetic correlation was significant at an FDR correction threshold of 10% (corrected within the respective category). For genetic correlation value and p-values across all 237 traits analyzed, see Table S4.

![Supplementary Fig. 6](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.15.21251499/F9.medium.gif)

[Supplementary Fig. 6](http://medrxiv.org/content/early/2021/02/18/2021.02.15.21251499/F9)

Supplementary Fig. 6 Correlation matrix for 31 physical activity phenotypes.
Correlation coefficients > 0.8 or < −0.8 are marked with *.

![Supplementary Fig. 7](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/18/2021.02.15.21251499/F10.medium.gif)

[Supplementary Fig. 7](http://medrxiv.org/content/early/2021/02/18/2021.02.15.21251499/F10)

Supplementary Fig. 7 Number of phenotypic PCs and cumulative proportion of variance explained.
At least 19 PCs are needed to explain 99% of the variance (red dashed line).

## Acknowledgements

The UK Biobank data was accessed via application ID 17712. Research of Drs. Guanghao Qi, Diptavo Dutta and Nilanjan Chatterjee was supported by an R01 grant from the National Human Genome Research Institute [1 R01 HG010480-01].

*   Received February 15, 2021.
*   Revision received February 15, 2021.
*   Accepted February 18, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1.Kyu, H. H., Bachman, V. F., Alexander, L. T., Mumford, J. E., Afshin, A., Estep, K., Veerman, J. L., Delwiche, K., Iannarone, M. L., Moyer, M. L., Cercy, K., Vos, T., Murray, C. J. & Forouzanfar, M. H. Physical activity and risk of breast cancer, colon cancer, diabetes, ischemic heart disease, and ischemic stroke events: systematic review and dose-response meta-analysis for the Global Burden of Disease Study 2013. BMJ 354, i3857 (2016).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzNTQvYXVnMDhfNC9pMzg1NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzAyLzE4LzIwMjEuMDIuMTUuMjEyNTE0OTkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

2.  2.Rovio, S., Kåreholt, I., Helkala, E. L., Viitanen, M., Winblad, B., Tuomilehto, J., Soininen, H., Nissinen, A. & Kivipelto, M. Leisure-time physical activity at midlife and the risk of dementia and Alzheimer’s disease. Lancet Neurol 4, 705–711 (2005).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1474-4422(05)70198-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16239176&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233002400014&link_type=ISI) 

3.  3.Smirnova, E., Leroux, A., Cao, Q., Tabacu, L., Zipunnikov, V., Crainiceanu, C. & Urbanek, J. The predictive performance of objective measures of physical activity derived from accelerometry data for 5-year all-cause mortality in older adults: NHANES 2003-2006. J Gerontol A Biol Sci Med Sci (2019).
    
    

4.  4.Leroux, A., Di, J., Smirnova, E., Mcguffey, E. J., Cao, Q., Bayatmokhtari, E., Tabacu, L., Zipunnikov, V., Urbanek, J. K. & Crainiceanu, C. Organizing and analyzing the activity data in NHANES. Stat Biosci 11, 262–287 (2019).
    
    

5.  5.Piercy, K. L., Troiano, R. P., Ballard, R. M., Carlson, S. A., Fulton, J. E., Galuska, D. A., George, S. M. & Olson, R. D. The Physical Activity Guidelines for Americans. JAMA 320, 2020–2028 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2018.14854&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30418471&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

6.  6.Lightfoot, J. T., DE Geus, E. J. C., Booth, F. W., Bray, M. S., DEN Hoed, M., Kaprio, J., Kelly, S. A., Pomp, D., Saul, M. C., Thomis, M. A., Garland, T. & Bouchard, C. Biological/Genetic Regulation of Physical Activity Level: Consensus from GenBioPAC. Med Sci Sports Exerc 50, 863–873 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1249/MSS.0000000000001499&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29166322&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

7.  7.Moore-Harrison, T. & Lightfoot, J. T. Driven to be inactive? The genetics of physical activity. Prog Mol Biol Transl Sci 94, 271–290 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/B978-0-12-375003-7.00010-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21036329&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

8.  8.De Moor, M. H., Liu, Y. J., Boomsma, D. I., Li, J., Hamilton, J. J., Hottenga, J. J., Levy, S., Liu, X. G., Pei, Y. F., Posthuma, D., Recker, R. R., Sullivan, P. F., Wang, L., Willemsen, G., Yan, H., De Geus, E. J. & Deng, H. W. Genome-wide association study of exercise behavior in Dutch and American adults. Med Sci Sports Exerc 41, 1887–1895 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1249/MSS.0b013e3181a2f646&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19727025&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000270202700008&link_type=ISI) 

9.  9.Hara, M., Hachiya, T., Sutoh, Y., Matsuo, K., Nishida, Y., Shimanoe, C. et al. Genomewide Association Study of Leisure-Time Exercise Behavior in Japanese Adults. Med Sci Sports Exerc 50, 2433–2441 (2018).
    
    

10. 10.Kim, J., Min, H., Oh, S., Kim, Y., Lee, A. H. & Park, T. Joint identification of genetic variants for physical activity in Korean population. Int J Mol Sci 15, 12407–12421 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/ijms150712407&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25026172&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

11. 11.Klimentidis, Y. C., Raichlen, D. A., Bea, J., Garcia, D. O., Wineinger, N. E., Mandarino, L. J., Alexander, G. E., Chen, Z. & Going, S. B. Genome-wide association study of habitual physical activity in over 377,000 UK Biobank participants identifies multiple variants including CADM2 and APOE. Int J Obes (Lond) 42, 1161–1176 (2018).
    
    

12. 12.Doherty, A., Jackson, D., Hammerla, N., Plötz, T., Olivier, P., Granat, M. H., White, T., van Hees, V. T., Trenell, M. I., Owen, C. G., Preece, S. J., Gillions, R., Sheard, S., Peakman, T., Brage, S. & Wareham, N. J. Large Scale Population Assessment of Physical Activity Using Wrist Worn Accelerometers: The UK Biobank Study. PLoS One 12, e0169649 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=org/10.1371/journal.pone.0169649&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

13. 13.Bycroft, C., Freeman, C., Petkova, D., Band, G., Elliott, L. T., Sharp, K., Motyer, A., Vukcevic, D., Delaneau, O., O’Connell, J., Cortes, A., Welsh, S., Young, A., Effingham, M., McVean, G., Leslie, S., Allen, N., Donnelly, P. & Marchini, J. The UK Biobank resource with deep phenotyping and genomic data. Nature 562, 203–209 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-018-0579-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30305743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

14. 14.Doherty, A., Smith-Byrne, K., Ferreira, T., Holmes, M. V., Holmes, C., Pulit, S. L. & Lindgren, C. M. GWAS identifies 14 loci for device-measured physical activity and sleep duration. Nat. Commun. 9, 5257 (2018).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

15. 15.Schrack, J. A., Kuo, P. L., Wanigatunga, A. A., Di, J., Simonsick, E. M., Spira, A. P., Ferrucci, L. & Zipunnikov, V. Active-to-Sedentary Behavior Transitions, Fatigability, and Physical Functioning in Older Adults. J Gerontol A Biol Sci Med Sci 74, 560–567 (2019).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

16. 16.Rock, P., Goodwin, G., Harmer, C. & Wulff, K. Daily rest-activity patterns in the bipolar phenotype: A controlled actigraphy study. Chronobiol Int 31, 290–296 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3109/07420528.2013.843542&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24517177&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

17. 17.McGregor, D. E., Palarea-Albaladejo, J., Dall, P. M., Stamatakis, E. & Chastin, S. F. M. Differences in physical activity time-use composition associated with cardiometabolic risks. Prev Med Rep 13, 23–29 (2019).
    
    

18. 18.Young, D. R. & Haskell, W. L. Accumulation of Moderate-to-Vigorous Physical Activity and All-Cause Mortality. J Am Heart Assoc 7, (2018).
    
    

19. 19.Jiang, L., Zheng, Z., Qi, T., Kemper, K. E., Wray, N. R., Visscher, P. M. & Yang, J. A resource-efficient tool for mixed model association analysis of large-scale data. Nat. Genet. 51, 1749–1755 (2019).
    
    

20. 20.Gusev, A., Ko, A., Shi, H., Bhatia, G., Chung, W., Penninx, B. W. et al. Integrative approaches for large-scale transcriptome-wide association studies. Nat. Genet. 48, 245–252 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3506&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26854917&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

21. 21.Finucane, H. K., Reshef, Y. A., Anttila, V., Slowikowski, K., Gusev, A., Byrnes, A. 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%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

22. 22.Finucane, H. K., Bulik-Sullivan, B., Gusev, A., Trynka, G., Reshef, Y., Loh, P. R. et al. Partitioning heritability by functional annotation using genome-wide association summary statistics. Nat. Genet. 47, 1228–1235 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3404&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26414678&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

23. 23.Watanabe, K., Taskesen, E., van Bochoven, A. & Posthuma, D. Functional mapping and annotation of genetic associations with FUMA. Nat. Commun. 8, 1826 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-017-01261-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

24. 24.Bulik-Sullivan, B., Finucane, H. K., Anttila, V., Gusev, A., Day, F. R., Loh, P. R., ReproGen, C., Psychiatric, G. C., Genetic, C. F. A. N. O. T. W. T. C. C. C., Duncan, L., Perry, J. R., Patterson, N., Robinson, E. B., Daly, M. J., Price, A. L. & Neale, B. M. An atlas of genetic correlations across human diseases and traits. Nat. Genet. 47, 1236–1241 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3406&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26414676&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

25. 25.Rock, P., Goodwin, G., Harmer, C. & Wulff, K. Daily rest-activity patterns in the bipolar phenotype: A controlled actigraphy study. Chronobiol Int 31, 290–296 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3109/07420528.2013.843542&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24517177&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

26. 26.Jones, S. E., Lane, J. M., Wood, A. R., van Hees, V. T., Tyrrell, J., Beaumont, R. N. et al. Genome-wide association analyses of chronotype in 697,828 individuals provides insights into circadian rhythms. Nat. Commun. 10, 343 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-08259-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30696823&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

27. 27.Gamazon, E. R., Wheeler, H. E., Shah, K. P., Mozaffari, S. V., Aquino-Michaels, K., Carroll, R. J., Eyler, A. E., Denny, J. C., GTEx, C., Nicolae, D. L., Cox, N. J. & Im, H. K. A gene-based association method for mapping traits using reference transcriptome data. Nat. Genet. 47, 1091–1098 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3367&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26258848&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

28. 28.GTEx Consortium. Genetic effects on gene expression across human tissues. Nature 550, 204–213 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature24277&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29022597&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000412829500039&link_type=ISI) 

29. 29.Bulik-Sullivan, B. K., Loh, P. R., Finucane, H. K., Ripke, S., Yang, J., Schizophrenia, Working Group of the Psychiatric Genomics Consortium, Patterson, N., Daly, M. J., Price, A. L. & Neale, B. M. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat. Genet. 47, 291–295 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3211&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25642630&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

30. 30.Dashti, H. S., Jones, S. E., Wood, A. R., Lane, J. M., van Hees, V. T., Wang, H. et al. Genome-wide association study identifies genetic loci for self-reported habitual sleep duration supported by accelerometer-derived estimates. Nat. Commun. 10, 1100 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-08917-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30846698&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

31. 31.Fregeau, B., Kim, B. J., Hernández-García, A., Jordan, V. K., Cho, M. T., Schnur, R. E. et al. De Novo Mutations of RERE Cause a Genetic Syndrome with Features that Overlap Those Associated with Proximal 1p36 Deletions. Am. J. Hum. Genet. 98, 963–970 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2016.03.002&link_type=DOI) 

32. 32.Lin, A. E., Ebert, G., Ow, Y., Preston, S. P., Toe, J. G., Cooney, J. P., Scott, H. W., Sasaki, M., Saibil, S. D., Dissanayake, D., Kim, R. H., Wakeham, A., You-Ten, A., Shahinian, A., Duncan, G., Silvester, J., Ohashi, P. S., Mak, T. W. & Pellegrini, M. ARIH2 is essential for embryogenesis, and its hematopoietic deficiency causes lethal activation of the immune system. Nat. Immunol. 14, 27–33 (2013).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23179078&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

33. 33.Cohn, R. D., Henry, M. D., Michele, D. E., Barresi, R., Saito, F., Moore, S. A., Flanagan, J. D., Skwarchuk, M. W., Robbins, M. E., Mendell, J. R., Williamson, R. A. & Campbell, K. P. Disruption of DAG1 in differentiated skeletal muscle reveals a role for dystroglycan in muscle regeneration. Cell 110, 639–648 (2002).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0092-8674(02)00907-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12230980&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000177943700010&link_type=ISI) 

34. 34.Sharif, K., Watad, A., Bragazzi, N. L., Lichtbroun, M., Amital, H. & Shoenfeld, Y. Physical activity and autoimmune diseases: Get moving and manage the disease. Autoimmun Rev 17, 53–72 (2018).
    
    

35. 35.Dhalwani, N. N., O’Donovan, G., Zaccardi, F., Hamer, M., Yates, T., Davies, M. & Khunti, K. Long terms trends of multimorbidity and association with physical activity in older English population. Int J Behav Nutr Phys Act 13, 8 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12966-016-0330-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26785753&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

36. 36.Duggal, N. A., Niemiro, G., Harridge, S. D. R., Simpson, R. J. & Lord, J. M. Can physical activity ameliorate immunosenescence and thereby reduce age-related multi-morbidity. Nat Rev Immunol 19, 563–572 (2019).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

37. 37.Vancampfort, D., Koyanagi, A., Ward, P. B., Rosenbaum, S., Schuch, F. B., Mugisha, J., Richards, J., Firth, J. & Stubbs, B. Chronic physical conditions, multimorbidity and physical activity across 46 low-and middle-income countries. Int J Behav Nutr Phys Act 14, 6 (2017).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

38. 38.Peters, H. P., De Vries, W. R., Vanberge-Henegouwen, G. P. & Akkermans, L. M. Potential benefits and hazards of physical activity and exercise on the gastrointestinal tract. Gut 48, 435–439 (2001).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ3V0am5sIjtzOjU6InJlc2lkIjtzOjg6IjQ4LzMvNDM1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDIvMTgvMjAyMS4wMi4xNS4yMTI1MTQ5OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

39. 39.Ciloglu, F., Peker, I., Pehlivan, A., Karacabey, K., Ilhan, N., Saygin, O. & Ozmerdivenli, R. Exercise intensity and its effects on thyroid hormones. Neuro Endocrinol Lett 26, 830–834 (2005).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.diabet.2010.05.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16380698&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

40. 40.Hawkins, V. N., Foster-Schubert, K., Chubak, J., Sorensen, B., Ulrich, C. M., Stancyzk, F. Z., Plymate, S., Stanford, J., White, E., Potter, J. D. & McTiernan, A. Effect of exercise on serum sex hormones in men: a 12-month randomized clinical trial. Med Sci Sports Exerc 40, 223–233 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1249/mss.0b013e31815bbba9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18202581&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000252608300006&link_type=ISI) 

41. 41.Ennour-Idrissi, K., Maunsell, E. & Diorio, C. Effect of physical activity on sex hormones in women: a systematic review and meta-analysis of randomized controlled trials. Breast Cancer Res 17, 139 (2015).
    
    

42. 42.Alessa, H. B., Chomistek, A. K., Hankinson, S. E., Barnett, J. B., Rood, J., Matthews, C. E., Rimm, E. B., Willett, W. C., Hu, F. B. & Tobias, D. K. Objective Measures of Physical Activity and Cardiometabolic and Endocrine Biomarkers. Med Sci Sports Exerc 49, 1817–1825 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1249/MSS.0000000000001287&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28398945&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

43. 43.Hackney, A. C. & Saeidi, A. The thyroid axis, prolactin, and exercise in humans. Curr Opin Endocr Metab Res 9, 45–50 (2019).
    
    

44. 44.Leroux, A., Xu, S., Kundu, P., Muschelli, J., Smirnova, E., Chatterjee, N. & Crainiceanu, C. Quantifying the Predictive Performance of Objectively Measured Physical Activity on Mortality in the UK Biobank. J Gerontol A Biol Sci Med Sci (2020).
    
    

45. 45.Loh, P. R., Tucker, G., Bulik-Sullivan, B. K., Vilhjálmsson, B. J., Finucane, H. K., Salem, R. M., Chasman, D. I., Ridker, P. M., Neale, B. M., Berger, B., Patterson, N. & Price, A. L. Efficient Bayesian mixed-model analysis increases association power in large cohorts. Nat. Genet. 47, 284–290 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3190&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25642633&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

46. 46.Loh, P. R., Kichaev, G., Gazal, S., Schoech, A. P. & Price, A. L. Mixed-model association for biobank-scale datasets. Nat. Genet. 50, 906–908 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0144-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29892013&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

47. 47.ENCODE, Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature 489, 57–74 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature11247&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22955616&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000308347000039&link_type=ISI) 

48. 48.Roadmap Epigenomics Consortium, Kundaje, A., Meuleman, W., Ernst, J., Bilenky M., Yen, et al. Integrative analysis of 111 reference human epigenomes. Nature 518, 317–330 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature14248&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25693563&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

49. 49.Watanabe, K., Umićević Mirkov M., de Leeuw, C. A., van den Heuvel, M. P. & Posthuma, D. Genetic mapping of cell type specificity for complex traits. Nat. Commun. 10, 3222 (2019).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31324783&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

50. 50.Buniello, A., MacArthur, J. A. L., Cerezo, M., Harris, L. W., Hayhurst, J., Malangone, C. 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%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

51. 51.Liberzon, A., Subramanian, A., Pinchback, R., Thorvaldsdóttir, H., Tamayo, P. & Mesirov, J. P. Molecular signatures database (MSigDB) 3.0. Bioinformatics 27, 1739–1740 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btr260&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21546393&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000291261300036&link_type=ISI) 

52. 52.Giambartolomei, C., Vukcevic, D., Schadt, E. E., Franke, L., Hingorani, A. D., Wallace, C. & Plagnol, V. Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet. 10, e1004383 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1004383&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24830394&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom) 

53. 53.Zheng, J., Erzurumluoglu, A. M., Elsworth, B. L., Kemp, J. P., Howe, L., Haycock, P. C., Hemani, G., Tansey, K., Laurin, C., Early Genetics and Lifecourse Epidemiology (EAGLE) Eczema Consortium, Pourcain, B. S., Warrington, N. M., Finucane, H. K., Price, A. L., Bulik-Sullivan, B. K., Anttila, V., Paternoster, L., Gaunt, T. R., Evans, D. M. & Neale, B. M. LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis. Bioinformatics 33, 272–279 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btw613&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27663502&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F02%2F18%2F2021.02.15.21251499.atom)

 [1]: /embed/graphic-6.gif