Research Paper Volume 9, Issue 12 pp 2504—2520

Human longevity: 25 genetic loci associated in 389,166 UK biobank participants

Luke C. Pilling1, , Chia-Ling Kuo2, , Kamil Sicinski3, , Jone Tamosauskaite1, , George A. Kuchel4, , Lorna W. Harries5, , Pamela Herd6, , Robert Wallace7, , Luigi Ferrucci8, , David Melzer1,4, ,

  • 1 Epidemiology and Public Health Group, University of Exeter Medical School, RILD Level 3, Royal Devon & Exeter Hospital, Exeter, EX2 5DW, UK
  • 2 Department of Community Medicine and Health Care, Connecticut Institute for Clinical and Translational Science, Institute for Systems Genomics, University of Connecticut Health Center, CT 06269, USA
  • 3 Center for Demography of Health and Aging, University of Wisconsin, Madison, WI 53706, USA
  • 4 UConn Center on Aging, University of Connecticut, Farmington, CT 06030, USA
  • 5 Institute of Biomedical and Clinical Sciences, University of Exeter Medical School, RILD Level 3, Royal Devon & Exeter Hospital, Exeter, UK
  • 6 La Follette School of Public Affairs and the Department of Sociology, University of Wisconsin, Madison, WI 53706, USA
  • 7 College of Public Health, University of Iowa, Iowa City, IA 52242, USA
  • 8 National Institute on Aging, Baltimore, MD 21224, USA

Received: September 26, 2017       Accepted: November 26, 2017       Published: December 6, 2017      

https://doi.org/10.18632/aging.101334
How to Cite

Copyright: Pilling et al. This is an open‐access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

We undertook a genome-wide association study (GWAS) of parental longevity in European descent UK Biobank participants. For combined mothers' and fathers' attained age, 10 loci were associated (p<5*10-8), including 8 previously identified for traits including survival, Alzheimer’s and cardiovascular disease. Of these, 4 were also associated with longest 10% survival (mothers age ≥90 years, fathers ≥87 years), with 2 additional associations including MC2R intronic variants (coding for the adrenocorticotropic hormone receptor). Mother’s age at death was associated with 3 additional loci (2 linked to autoimmune conditions), and 8 for fathers only. An attained age genetic risk score associated with parental survival in the US Health and Retirement Study and the Wisconsin Longitudinal Study and with having a centenarian parent (n=1,181) in UK Biobank. The results suggest that human longevity is highly polygenic with prominent roles for loci likely involved in cellular senescence and inflammation, plus lipid metabolism and cardiovascular conditions. There may also be gender specific routes to longevity.

Introduction

For centuries there has been great interest in factors enabling some individuals to live to very old ages. Extreme longevity can be achieved in animal models by knocking out many specific genes [1]. In humans there is a growing number of inherited (germline) genetic variants linked to overall [2] and longer human lifespan [3,4] in conventional genome wide association studies (GWAS), and in studies weighted for the large numbers of previously reported variant associations with related diseases and traits [5,6].

Human longevity or lifespan has moderate heritability in twin studies, estimated at 20-30% [7], with health related behaviors, environmental exposures and chance likely explaining much of the remaining variability. GWAS have mostly studied general survival or compared longer lived individuals (e.g. aged 80+, 90+ or centenarians) to younger controls. An alternative approach to studying older individuals themselves is to study parental longevity. Biological offspring of longer lived parents inherit a combination of parental genetic variants, and tend to remain healthier and live longer than the offspring of shorter lived parents [8]. In the US Health and Retirement Study (HRS), we found that mortality in offspring declined progressively with later parental ages of death [9]. In addition, offspring of longer lived parents had progressively lower incidence of cardiovascular disease and cancers, plus reduced rates of cognitive impairment [10]. These findings were replicated in 186,151 non-adopted UK Biobank participants [11]. Statistical associations of mothers’ and fathers’ attained ages with offspring health outcomes in HRS were additive. These findings suggest that combined parental longevity is a quantitative trait of slower aging, likely with a significant genetic component.

We [12] and separately Joshi et al. [13] previously undertook a GWAS of parental longevity in the first one third UK Biobank (UKB) participants aged 55-70 [12]. We determined the percent of variance of parental age at death explained by the directly genotyped variants on the UKB arrays (n=845,997): the variance attributed was greatest for the sum of normalized mother's and father's age at death (10.2%, SD=1.26%) and smaller for mothers (6.08% SD=1.21%) and fathers (5.79% SD=1.23%) separately.

In the current analysis we aimed to extend our previous GWAS to include data from the newly available fully genotyped UKB cohort, essentially tripling our previous sample size. We aimed to identify all common genetic variants (prevalence ≥0.1%) associated with longer parental lifespan. We studied combined parents’ attained ages but also tested associations with both parents having longest 10% survival, to establish whether variants were associated with long life rather than early death. We also examined mothers’ and fathers’ age of death separately. We also tested a combined genetic risk score of our findings in two US aging cohorts, the HRS and the Wisconsin Longitudinal Study (WLS).

Results

Subjects meeting inclusion criteria (Table 1) were 389,166 UK Biobank participants of European descent with data recorded on parents' current ages or parents' ages of death (n=181,048 had one or more living parents). GWAS analyses had >99% power to detect an allele of 1% minor allele frequency, accounting for 0.1% of the variance in the phenotype, after correction for multiple statistical testing (alpha=5x10-8). Variants with minor allele frequencies of >0.1 were included.

Table 1. Summary statistics for 389,166 UK Biobank participants of European descent with genetic and parental age data.

nmin - maxmean (SD)
Age of study participant (years)389,16640 - 7356.66 (8)
BMI (Kg/m2)387,95112.12 - 74.6827.31 (4.73)
Mother's age at death (years)246,94156 - 10777.39 (9.77)
Mother's age if alive (years)165,99660 - 10578.58 (8.15)
Father's age at death (years)317,65246 - 10672.18 (11.05)
Father's age if alive (years)97,65960 - 10378.02 (7.37)
Combined parents age at death
Years (average)208,11851 - 10175.41 (7.79)
Sum of z-scores208,118-3.12 - 3.310.01 (1)
Combined parents age (alive or dead)
Years (total)389,166102 - 202151.57 (14.98)
n%
Sex
Females212,67254.65
Males176,49445.35
Smoking status
Never212,83354.87
Former136,64735.23
Current38,4149.9
Highest education level attained
None62,51916.21
Secondary65,58517
College-level71,04818.42
Professional / University186,57248.37
Both parents top 10% survival
No (both died <80)79,76791.74
Yes (≥90 mothers, ≥87 fathers)7,1828.26
Combined parents age (alive or dead)
Either alive181,04846.52
Both dead208,11853.48

GWAS of parental attained age

We examined associations with both parents’ combined attained ages (including ages of living parents), using Martingale residuals from Cox’s proportional hazards regression models, following the approach used by Joshi et al. [13]. In this analysis we found 10 loci associated at genome wide significance (p<5x10-8) (Table 2; Supplementary Table 1 for details of all significant variants; full results available to download https://doi.org/10.6084/m9.figshare.5439382.v1).

Table 2. GWAS associations with parents' attained age phenotypes.

P-values for phenotype:Effect for the analysis in bold
Lead SNPCHRBPA1A0Parents' attained ageParents' age at deathBoth parents top 10%Mothers' attained ageFathers' attained ageA1FBetaSEImplicated gene(s)
Combined parents' age results
Common & in GWAS catalogue of findings
rs6026331109821511TG2.7E-082.9E-032.4E-032.1E-033.6E-090.217-0.01500.0027CLESR2 ... PSRC1
rs28383322632592796CT5.3E-119.8E-032.0E-043.5E-092.0E-040.7830.01820.0028HLA-DRB1… HLA-DQA1
rs557304996161005610CT1.7E-185.4E-112.7E-085.5E-105.9E-170.920-0.03610.0041LPA
rs1556516922100176GC4.7E-161.5E-102.5E-088.7E-051.3E-200.502-0.01810.0022CDKN2B-AS1 (ANRIL)
rs713782812111932800CT3.4E-143.1E-083.2E-062.2E-072.7E-120.4830.01700.0022SH2B3 / ATXN2
rs13172861578896129AG1.2E-261.8E-114.0E-092.1E-071.5E-310.668-0.02540.0024CHRNA3
rs175148461591416550CA7.1E-106.6E-072.0E-047.8E-041.2E-070.530-0.01390.0023FURIN
rs4293581945411941TC1.4E-745.0E-227.1E-461.3E-681.9E-280.845-0.05660.0031APOE / APOC1
Common
rs16278046107400428CA1.1E-044.0E-084.5E-065.7E-032.4E-040.6840.01970.0035BEND3
rs7844965827442064GA7.7E-091.5E-046.0E-044.9E-042.2E-070.7680.01540.0027EPHX2
rs619789281475321714TC2.0E-087.7E-053.3E-071.0E-041.9E-070.6880.01360.0024PROX2
rs289261731813886719GA6.4E-062.4E-052.3E-082.4E-042.5E-020.953-0.02640.0047MC2R
Rare (<2%)
rs146254978174867799TC2.1E-024.6E-081.9E-034.3E-013.0E-040.980-0.06920.0126FPGT / TNNI3K
rs13913745911119269958GA5.3E-031.1E-012.7E-083.0E-023.1E-020.996-0.09060.0162USP2-AS1
Loci only emerging in mothers' analysis
rs3130507631147476GA3.5E-066.1E-028.9E-012.1E-102.4E-010.720-0.01550.0024PSORS1C3 … [x]
13:318715141331871514TG5.6E-044.5E-017.0E-014.7E-083.4E-010.9950.09790.0180B3GALTL
rs619496501364836488TC4.2E-051.5E-031.5E-032.9E-086.7E-010.9380.02510.0045(intergenic)
Loci only emerging in fathers' analysis
rs3131621631425499AG8.9E-081.9E-024.8E-013.5E-053.6E-080.629-0.01240.0022MICA … MICB
rs13262617859838133AG9.6E-037.3E-031.3E-018.1E-013.1E-080.9700.03600.0064TOX
rs6190574711113639842AC3.6E-053.5E-027.0E-022.2E-015.5E-090.809-0.01600.0028ZW10
rs740114151547660194GA7.0E-063.9E-022.6E-025.0E-021.4E-080.881-0.01920.0034SEMA6D
rs124619641941341229AG1.8E-041.5E-017.8E-014.8E-018.2E-090.504-0.01260.0022EGLN2 … CYP2A6
rs744449831945745607TC4.1E-059.1E-031.1E-016.9E-029.1E-090.741-0.01430.0025EXOC3L2 … MARK4
rs61087842010964366TC7.7E-044.3E-015.7E-027.8E-011.2E-100.586-0.01430.0022C20orf187
rs22735002061986949TC4.4E-078.9E-021.8E-013.8E-031.5E-080.856-0.01760.0031CHRNA4
Lead SNP = SNP with lowest p-value in 1Mb region
BP = Genomic position, build 37 (hg19), A1F = frequency of A1, the effect allele, P = BOLT LMM p-value (A1F, BETA, and SE are for p-value in bold), SE = Standard Error
Beta = A1 effect on outcome. For attained age (Martingale residuals) a negative BETA = reduced hazard, i.e. increased attained age
Implicated gene(s) = Variants in locus usually intronic, or exonic (indicated by underlining). If two genes separated by dots, the SNP is intergenic

Associations with APOE related rs429358 (p=5x10-22) variant were confirmed, as were associations with the CHRNA3/5 (nicotinic acetylcholine receptor) locus (rs1317286, 33.2% minor allele frequency, p=1.2x10-26). Associations were also found for six other loci already implicated in disease specific GWAS; the strongest association was with rs55730499 (8% frequency, p=1.7x10-18) in an intron of LPA gene, associated with lipoprotein A and LDL cholesterol levels, plus CAD [14]. Variants in the 9p21 region, which includes CDKN2B-AS1 (ANRIL) and CDKN2A/B, were also associated (lead SNP rs1556516, 49.8% frequency, p=4.7x10-16); this locus is involved in cellular senescence and associated with myocardial infarction and peripheral vascular disease [14]. Genetic variants across the 9p21.3 locus have previously also been associated with type-2 diabetes and cancers [14] and in longevity/survival analyses (see below), but only the CAD-associated loci were associated (p<5x10-8) with parent’s lifespan (Figure 2; Supplementary Table 2).

Manhattan plots for GWAS of parent’s lifespan. (A) combined parents' attained age GWAS including living parents, (B) combined age at death of mother and father, (C) both parents reached top 10% of survival (mothers reached ≥90 years, father reached ≥87 years). Red points are variants (p-5) previously associated with other traits. The blue line indicates p=5x10-8. The Y-axis has been capped at 20 (p=1x10-20) to aid in visualising the results.

Figure 1. Manhattan plots for GWAS of parent’s lifespan. (A) combined parents' attained age GWAS including living parents, (B) combined age at death of mother and father, (C) both parents reached top 10% of survival (mothers reached ≥90 years, father reached ≥87 years). Red points are variants (p<1x10-5) previously associated with other traits. The blue line indicates p=5x10-8. The Y-axis has been capped at 20 (p=1x10-20) to aid in visualising the results.

LocusZoom plot of the CDKN2A/B/ANRIL gene cluster and known disease variants. Known disease-associated genetic variants have been highlighted in the 9p21.3 region, in addition to the lead SNP identified in our analysis. Association p-values indicate the strength of association with parents' attained age. The variants known to influence coronary artery disease (CAD) and related-traits are those most-associated with parents’ attained age. Variants previously associated with cancer are markedly less associated with parents’ attained age, with type-2 diabetes variants not associated (p>0.05). Details, including references, are included in Supplementary Table 2.

Figure 2. LocusZoom plot of the CDKN2A/B/ANRIL gene cluster and known disease variants. Known disease-associated genetic variants have been highlighted in the 9p21.3 region, in addition to the lead SNP identified in our analysis. Association p-values indicate the strength of association with parents' attained age. The variants known to influence coronary artery disease (CAD) and related-traits are those most-associated with parents’ attained age. Variants previously associated with cancer are markedly less associated with parents’ attained age, with type-2 diabetes variants not associated (p>0.05). Details, including references, are included in Supplementary Table 2.

Associations with parents' attained age were also present for rs7137828 (48.3% frequency, p=3.4x10-14) near ATXN2, coding for Ataxin-2 which is involved in endocytosis and is associated with blood traits and autoimmune disease [14]. SNPs in LD (R2=0.91) with rs7137828 include rs3184504 (48% frequency, p=4x10-8), a missense variant in the SH2B3 gene, associated with blood pressure, maternal birth weight, TNF-α levels, colorectal cancer, Inflammatory Bowel Disease, LDL cholesterol levels and stroke [14]. Other known loci included rs28383322 (located between HLA-DRB1 and HLA-DQA1, p=5.3x10-11), previously associated with atopic dermatitis; rs17514846 (located in an intron of FURIN, p=7.1x10-10), previously associated with blood pressure; and rs602633 (located between CLESR2 and PSRC1, p=2.7x10-8), previously associated with CAD and LDL [14].

Two associations were discovered for novel variants not previously linked in GWAS to other traits or diseases; rs7844965, located in an intron of EPHX2 (23.2% frequency, p=7.7x10-9), and rs61978928, located in the 3’ UTR of gene PROX2 (31.2% frequency, p=2x10-8).

A parents' attained age genetic risk score in independent cohorts

We calculated a genetic risk score (GRS) combining the 10 attained parental age associated variants, weighted by the estimated effect sizes. Genotyped samples of European descent from the Health and Retirement Study (HRS, n=12,940 in genetic analyses) and the Wisconsin Longitudinal Study (WLS, n=8,952 from 6,683 families) were too small to test individual variant associations: the sample size needed for 80% power to detect an association between the APOE variant and parents' lifespan at p<0.05 (R2=0.0005 in UK Biobank) is 15,598, but much larger samples were required for the smaller-effect variants discovered in UK Biobank. We therefore tested the parents’ attained age GRS; the sample size needed for an 80% chance of detecting an association of this magnitude (R2=0.002 in UKB) with p<0.05 was 3,924. In Cox’s proportional hazards survival models against both dead and alive parents attained ages in HRS (n=12,345, number of events=9,809), increasing GRS was associated with longer survival (mortality Hazard Ratio=0.961, 95% CIs=0.942 to 0.980, p=6x10-5, adjusted for offspring age, sex, and genetic principal components 1-5). In the WLS (n=7,437, number of events=6,683, errors clustered by family), increasing GRS was associated with longer survival (mortality Hazard Ratio=0.963, 95% CIs=0.938 to 0.988, p=4x10-3, adjusted for age at saliva collection, sex, and genetic principal components 1-5).

Functional implications and pathways analysis

We investigated effects of the parents' attained age SNPs on expression of genes in multiple tissues, such as whole blood [15], using the FUMA tools [16]. Expression of 53 genes were associated, including 8 of the 10 nearest genes to the variants identified above (Supplementary Table 5 and Supplementary Table 6).

In MAGENTA gene-set analysis, only two pathways were significantly enriched after correction for multiple testing, both related to nicotinic acetylcholine receptor binding and downstream events (Supplementary Table 4).

Sensitivity analyses of combined parents' attained age analyses

We generated QQ plots for all GWAS (Supplementary Figure 2) and observed some genomic inflation (Lambda GC for parents' attained age = 1.199), consistent with human longevity being highly polygenic, but perhaps also representing some effect of relatedness in inflating results. The BOLT-LMM software used for analysis accounts for relatedness in the methods, but we investigated further; the 10 variants associated with attained parental age were tested in Cox proportional hazards (CoxPH) models against parents’ attained age, to assess whether the associations observed with rank-normalized Martingale residuals were robust. These models also excluded related (≥3rd-degree) participants, leaving a reduced sample size of n=326,864, including n=153,101 with censored data due to at least one parent still being alive. Five of the 10 loci remained genome wide significant (p<5x10-8) in these models; APOE, CHRNA3, LPA, ANRIL, and ATAXN2/SH2B3 (Supplementary Table 3) with the remaining five directionally consistent and associated at p≤3.1*10-6.

We performed sensitivity analyses only including parents who had died (Figure 1B; Table 2), as the ages of living parents at the UKB study baseline are clearly an incomplete measure of their longevity. This analysis included a smaller sample size of n=208,118 UKB participants with both parental ages of death known. The same five of the ten loci reported above remained genome-wide significant (APOE, CHRNA3, LPA, ANRIL, and SH2B3) with the remaining loci remaining nominally significant. Two additional variants not identified in the previous analysis emerged as candidates for follow-up: rs1627804 (32% frequency, p=4x10-8), in an intron of BEND3, a transcriptional repressor which associates with the NoRC (nucleolar remodeling complex) complex and playing a key role in repressing ribosomal DNA transcription; and rs146254978 (2% frequency, p=4.6x10-8), located in an intron of both FPGT and TNNI3K.

GWAS of both parents achieving longest 10% survival

We next tested binary associations with having both parents surviving to the longest 10% of lifespans, using sex-specific age-at-death cut-offs (father ≥87, mother ≥90, including living parents aged over these age cut-offs), to assess whether associations were attributable to longer lives rather than early deaths. This analysis confirmed associations (Figure 1C; Table 1) with variants in the APOE, CHRNA3, LPA and CDKN2B-AS1, with a suggestive association with ATXN2 (p=3.4x10-6). Two further loci just reached significance in this analysis: rs28926173 (4.7% frequency, p=2.3x10-8) in an intron of gene MC2R, and rs139137459 (0.4% frequency, p=2.7x10-8) in an intron of USP2-AS1 (near USP2). The MC2R variant was associated at p=6.4x10-6 with parents’ attained age, but the USP2-AS1 variant had a low population frequency (0.4%) and was not significant in the continuous trait analyses.

Genetic risk score associations with centenarian status and offspring health status in UKB

There were 691 UKB participants with a centenarian mother (the longest lived 0.1% of maternal deaths in the sample), and 747 with a father aged ≥98 years (the corresponding age cut-off of longest lived 0.1%); only five participants had both parents in the longest lived 0.1%. Although this relatively small group (0.3% of the full attained age sample) were included in the discovery analysis, we tested the GRS for association with having a centenarian parent (in the un-related UK Biobank participants, n=1,181), against controls whose (both) parents died before age 80 (n=66,280). Increasing attained age GRS was associated with increased odds of having a centenarian parent (OR per standard deviation of GRS=1.193, 95% CIs=1.130 to 1.260, p=3.0x10-10).

GWAS of mothers' and fathers' age at death separately

We tested for parental gender differences by conducting separate GWAS for mothers' and fathers' attained ages (including living parents). Eleven additional genetic variants were identified, three only associated with mother’s age at death, eight associated with father’s age at death (Table 2). Results included variants previously associated with the autoimmune diseases psoriasis and idiopathic membranous nephropathy (for mothers), and an independent nicotine-dependence locus on chromosome 20 for fathers. Novel gene-longevity associations included genes MICA/B (a locus previously linked to auto-immune conditions rheumatoid arthritis and multiple sclerosis), TOX, ZW10, SEMA6D, EGLN2/CYP2A6, EXOC3L2/MARK4 and C20orf187.

Comparison with recently reported GWAS analyses of longevity

A recent large analysis by Joshi et al. [2], with up to 294,998 participants in their analysis of European parents' attained age, included the first third of UK Biobank genotyping data plus 26 other cohorts. For survival past age 40 years, associations were confirmed for APOE and CHRNA3/5, plus novel associations with overall longevity (deaths or survival to age 40 years and over) with HLA-DQA1/DRB1 and LPA, and with replication of CDKN2A/B and SH2B3. While these results overlap with our attained age analysis, they also reported an association with FOXO3A variants, but we found no significant associations for this locus (p>0.01) (see Supplementary Table 2 for details). No analyses of longer or extreme survival were reported.

McDaid et al. [6] undertook a Bayesian association scan of general lifespan (in the first third of UK Biobank), weighted for 58 disease related GWAS studies, and identified 16 individual variants at genome wide significance. Of these, four (LPA, CDKN2BAS, CHRNA5 and APOC1) reached conventional GWAS significance in our analyses, and CELSR2 only reached significance for mothers' age at death.

While the above studies focused mainly on overall lifespan, for aging research survival to extreme ages is of special interest. In a powerful study of extreme longevity, Fortney et al. used samples including 801 centenarians and 5406 subjects aged 90 plus. Disease GWAS results were used to weight genetic variants, using a false discovery rate of 10% [5], at which associations with APOE, SH2B3/ATXN2, ANRIL, LPA, and HLA were reported (although only APOE reached conventional genome wide significance). Fortney et al. also identified variants at KCNT2 and FADS1, which were not associated with parents’ attained age in our analysis, with the ABO SNP (rs514659) only just reaching nominal significance (p=0.015).

A previous meta-analysis comparing older individuals (age 80 or 90+) to younger controls reported a genome wide association at the 5q33.3 locus, but no significant associations were found for the other two loci (p>0.01).

Various other variants, including in mitochondrial DNA, have been suggested as associated with longevity. In our GWAS we tested associations with mitochondrial variants, plus variants on the X and Y chromosomes, using the directly genotyped microarray data: none were significantly associated (p>1x10-5) with the studied phenotypes.

Disease-associated genetic variants; candidates for future work

We investigated genetic variants associated with parents' attained age (p<1x10-5) that have previously been associated with other traits for promising candidates for future work and identified 12 independent signals (Table 3). These included SNPs associated with BMI (rs11126666 [17],), breast cancer (rs11552449 [18],), educational attainment (rs7593947 [19],), and kidney function (estimated Glomerular Filtration Rate, rs4014195 [20],), in addition to a variant in IL6R (rs6689306) associated with CAD.

Table 3. Variants implicated in disease and trait GWAS that were associated with parents' attained age at near genome wide significance (p>1x10-8 but p<1*10-5).

SNPCHRBPA1A0A1FREQBETASEPGenePrevious associations with locus
rs115524491114448389CT0.826-0.01340.00306.3E-06DCLRE1BBreast Cancer [18]
rs66893061154395946AG0.4190.01110.00231.5E-06IL6RCAD [40]
rs11126666226928811GA0.745-0.01170.00264.5E-06KCNK3BMI [17]
rs7593947260704933AT0.526-0.01100.00231.2E-06BCL11AEducational attainment [19]
4:77410318477410318CA0.5440.01000.00229.4E-06SHROOM3Red blood cell count [41]
rs1986734477420784CT0.4760.01010.00226.9E-06SHROOM3Creatinine [42]
rs806794626200677AG0.7280.01150.00255.0E-06HIST1H2BF … HIST1H4EHeight [43]
rs40141951165506822CG0.656-0.01080.00243.7E-06RNASEH2C … AP5B1eGFR [20]
rs63318511100593538GC0.285-0.01120.00258.9E-06ARHGAP42DBP, SBP [44]
rs562898211911188247GA0.8820.01550.00357.3E-06SMARCA4 … LDLRCAD [40]
rs65117201911202306GT0.8810.01650.00351.8E-06LDLRLDL [45]
rs177249921918454825AG0.7330.01240.00258.6E-07PGPEP1BMI [17]
BP = Genomic position, build 37 (hg19), BETA = Effect on outcome with respect to A1 (the effect allele), P = BOLT LMM p-value

Discussion

In this analysis we aimed to extend our previously published GWAS study of parental longevity to include the newly available full UKB data, achieving an unprecedented sample size for longevity studies. We have confirmed associations with APOE and CHRNA3 variants, and additionally identified 8 novel genome wide significant variants for our principal phenotype of combined parents’ attained age, plus several additional variants for our related longevity phenotypes. Overall, the findings suggest that human longevity is a polygenic trait, influenced by large numbers of modest effect variants. Many of the variants identified have been previously implicated in age-related diseases or traits, notably with coronary artery disease and lipid levels. Several loci are near genes involved in inflammation, an important pathway in the development of biological aging [21]. In additional analyses of mothers’ and fathers’ attained ages separately we identified 11 additional genetic variants, three only associated with mother’s age at death and eight associated only with father’s age at death.

Parental lifespan-associated variants included rs1556516, located in an intron of CDKN2B-AS1 (ANRIL), a long non-coding RNA; variants in this region, 9p21.3, have previously been associated with CAD, type-2 diabetes, and cancer [14]. ANRIL is known to regulate neighbor tumour suppressors CDKN2A (encoding proteins p14arf and p16) and CDKN2B (encoding protein p15) epigenetically through directing polycomb gene regulators to modify the chromatin state [22]. ANRIL therefore plays a key role in regulation of cell proliferation and senescence. It is interesting to note that removal of CDKN2a expressing senescent cells achieved significant rescue of aging phenotypes in mice [23]. Intriguingly, the variants in 9p21 associated with lifespan are predominantly related to arterial disease (coronary and peripheral arterial disease), rather than those influencing cancer or type II diabetes risk (Figure 2). Decreased expression of specific ANRIL isoforms is reported to affect expression of ADIPOR1 (adiponectin receptor with a role in fatty acid catabolism and inflammation), VAMP3 (membrane protein with a role in phagocytosis) and C11orf10 (in proximity to genes associated with metabolic syndrome) [24]. The lifespan-increasing allele for the lead SNP (rs1556516-G) is known to increase expression of ANRIL [25]; atherosclerosis-related phenotypes are associated with lower expression of ANRIL, suggesting a causal mechanism of action for this SNP on lifespan.

rs11065979 (chromosome 11) is located near the ATXN2 gene, which is known to be associated with autoimmune disease and blood traits [14]. The variant is also highly correlated (R2=0.91) with an exonic variant in SH2B3, previously linked to CAD, blood pressure, TNF-α, cell counts and colorectal cancer [14]. SH2B3 encodes the lymphocyte adapter protein (LNK), which plays a key role linking inflammation and hypertension, particularly in the kidney and blood vessels [26].

EPHX2 is a hydrolase involved in the degradation of vasoactive and anti-inflammatory acids, with possible function in LDL-cholesterol metabolism [27]. Exonic polymorphisms in EPHX2 are associated with coronary artery calcification [28]. PROX2 is a transcription factor of unknown function implicated in cardiac muscle and neuron differentiation by gene ontology [29]. Genetic variants in FURIN, a gene located in chromosome 15, were previously associated with blood pressure and are here associated with parents’ attained age; the protein is a widely-expressed membrane-bound protease, and is essential for regulatory and effector T-cell function [30].

Other novel genes implicated include MC2R, where the rare allele is associated with 15.7% increased likelihood of both parents reaching the top 10% of survival. This gene encodes the adrenocorticotropic hormone receptor, involved in the secretion of cortisol. This pathway is key to the stress response, which is regarded as a pillar of the aging process [31]. Lack of this functioning receptor can result in chronic hormone-imbalance and various conditions, including subclinical Cushing syndrome, plus impaired blood sugar regulation and immune system function [32].

For four of the novel parents’ attained age-associated variants it is the rare allele that is associated with increased lifespan; rs602633 near CLESR2 and PSRC1; rs28383322 near HLA-DRB1 and HRL-DQA1; rs7844965 in EPHX2; and rs61978928 in PROX2. This is unusual, as for the other genetic variants it is the common allele that confers the advantage (the rare variant tends to be associated with increased disease-risk, such as CAD).

We have shown that a genetic risk score based on the 10 variants associated with combined parents' age of death is associated with parental survival in two independent US aging cohorts, suggesting that our overall findings are robust. In addition to the more common variants identified (i.e. minor allele frequency >1%) there were several rarer variants, such as rs139137459 in USP2-AS1 (0.4% frequency): these associations are less certain, as the GWAS significance level of p<5*10-8 may not be sufficiently rigorous for testing low frequency variants. The identified SNPs previously linked to GWAS traits (such as a variant in IL6R) have a high prior probability and are especially likely to be robust.

This study is limited to European ancestry (white British) UK Biobank participants, and results may not be applicable to other populations. UK Biobank is a healthier volunteer study with a sample that was less likely to be obese, to smoke or to drink regularly, compared to the general population [33], perhaps reducing our power to detect variants associated with health behaviors, although the sample still contains many participants who were exposed to these factors. We were not able to verify reported parental ages or ages of death, perhaps introducing misclassification. Also, ages of birth of deceased parents were not recorded, so it is not possible to compare study and national cohort survival.

We are limited by the coverage of the available genotyping microarray, and have used data on imputed SNPs from the Haplotype Reference Consortium panel, with imputation using the UK10K reference panel to follow. As with all GWAS discovered loci, the identification of affected genes is not straightforward, but we have followed the convention of highlighting the nearest genes, which for 8 of the 10 attained age variants was confirmed by the eQTL data.

We studied the normal range of parental ages at death, as well as the top end of longevity (top 10% of survival), to establish whether associations were driven by early deaths. Our longest lived 10% analysis is similar in age (87+ years for fathers, 90+ for mothers) to those studied in many of the previous GWAS of longevity [3,4]. Analyses with bigger samples of those attaining extreme ages are needed, but unfortunately our sample with extreme surviving parents (e.g. centenarians) is too small for GWAS analyses.

In addition to these limitations, our analysis has several strengths. Our sample size yielded unprecedented power to detect variant associations: indeed the power (>99% to detect an allele of 1% minor allele frequency accounting for 0.1% of phenotype variance) is sufficient to suggest that we have identified all moderate to larger effect common genotyped or imputed variants in our studied population. Our analyses were undertaken in a single cohort: previous longevity GWAS combined data from diverse cohorts, which has been suggested as explaining difficulties in detecting and replicating longevity variants [34].

Future work

Replication and extension of the novel associations in other large cohorts and in different ancestries is needed, especially for variants narrowly achieving genome wide significance. The associations identified should be examined in future studies of age-related traits and disease. More work is needed to characterize the variants identified for biological functions; this could include epigenetic and gene expression studies, such as analysis of allele-specific expression imbalance in heterozygote individuals to determine the gene of action for genetic variants [35]. The variants identified in our analysis may provide useful information for designing targeted prevention or precision treatments to promote human longevity, perhaps particularly related to targets in inflammatory pathways and cardiovascular disease.

Conclusions

Our analysis of parental longevity involved the largest sample size thus far achieved, with ample power (>99%) to detect even modest associations. The results indicate that human longevity is a polygenetic trait, including likely prominent effects involving cell senescence, inflammation, autoimmunity and the stress response, plus lipid levels and vascular disease. The results suggest that the identified variants are involved in achieving long lifespans, and are not limited to causing early deaths. The results also suggest that there may be gender differences in routes to longevity, with autoimmune variants being particularly important in women. Overall these results are consistent with some previous analyses in humans and with certain key hallmarks or pillars of the biology of aging. The variants identified may be useful for identifying targeted prevention and precision treatments for extending human lifespans.

Methods

Parental age at death and longevity phenotypes

Participants were asked the age at which their parents had died (or their current age if still alive). Analyses were performed separately on mother’s age at death and father’s age at death, and also on a combined phenotype. To reduce the effect of higher ages at death of mothers (compared to the fathers) we first z-transformed the mothers and fathers age at deaths before combining the z-scores into a single summed phenotype. Offspring of parents who died prematurely were excluded because the cause of death of the participant’s parents was not asked, so we could not exclude accidental deaths explicitly. To determine the age at death considered to be too young to be part of the normal distribution of age-related mortality, we used previously described methods to define a normal distribution of ages at death for mothers and fathers separately, as described previously [9]). In brief, non-linear least square regression models were used to fit a normal curve on age at death data from the modal (M) age at death to the oldest ages; the younger (left) half of the curve is then imposed, to avoid the severe skewing of the distribution caused by the excess of deaths below the modal age. The analysis identified the following early death cut-points (based on modal age at death minus 1 Standard Deviation): mothers 57 years, fathers 46 years.

In addition, combined parents’ age (including alive parents) were derived from Cox’s proportional hazards regression models using Martingale residuals, as described by Joshi et al. [13]. Events were defined as the both parents having died. We also defined a binary phenotype by determining the top 10% of age at death for mothers (≥90 years) and fathers (≥87 years) separately. Controls were participants where both parents died before the age of 80. We were able to increase the number of “long lived” parents by including alive participants beyond the cut-off.

UK Biobank

Between 2006 and 2010, 503,325 volunteers (aged 45 to 69 years old) were recruited from across the United Kingdom to the UK Biobank study [36]. Of these, 389,166 participants met the inclusion criteria for at least one analysis: participants with complete genetics data and either date of death of their parents or their current age if alive. Participants were excluded if they reported themselves as adopted, or whose parents died more than 1 standard deviation below the modal age of death (mothers <57 years, fathers <46 years). Several longevity phenotypes are defined below based on the age at death of the participant’s parents.

Genetic data was available on 488,377 UK Biobank participants after genotype calling and quality control performed centrally by the UK Biobank team [37]. We selected 451,447 participants identified as ‘white European’ through self-report and verified through principal components analysis based on genotypes. Briefly, principal components were generated in the 1000 Genomes Cohort using high-confidence SNPs to obtain their individual loadings. These loadings were then used to project all of the UK Biobank samples into the same principal component space and individuals were then clustered using principal components 1 to 4. Related individuals were identified through kinship analysis, although these participants were included in genome-wide analysis using BOLT-LMM (see below), those related to the third-degree or closer were excluded in sensitivity analyses.

Imputation of 39,235,157 genetic variants from the Haplotype Reference Consortium panel was performed using IMPUTE4 centrally by the UK Biobank team [37]. After filtering for variants with MAF ≥0.1%, missingness <1.5%, imputation quality >0.1 and with Hardy-Weinberg equilibrium (HWE) P>1x10-6 within the white British participants 11,516,125 imputed autosomal variants were eligible for the analyses.

We also utilized data directly from the microarrays for variants on the X (n=19,381) and Y (n=284) chromosomes, and on the mitochondrial genome (n=135), which were unavailable in the imputed dataset.

Health and Retirement Study (HRS)

In HRS v2 data, phase 1-3 samples were combined (n=15,708) and SNPs were kept if no discordant genotype call among duplicate samples (2,315,518 SNPs). Related or duplicated participants and genotyping controls were removed, as well as those with a missing call rate equal to or greater than 2%. After the sample quality control, 15,454 samples were kept and 84% of them (n=12,940) were of European descent.

2,075,208 out of 2,315,518 SNPs passed the composite filters by the University of Washington (Hardy-Weinberg equilibrium p-value < 1e-4, missing call rate ≥ 2%, sex difference in allelic frequency ≥ 2% and sex difference in heterozygosity > 0.3) and had a minor allele frequency greater than zero. As recommended, SNPs with a minor allele frequency smaller than 2% further were removed. At the end, 1,437,608 SNPs were carried forward for SNP pruning. The purpose is to choose less correlated SNPs for principal components analysis.

The SNP pruning was performed following the procedure: a) consider a window of 50 SNPs, b) calculate LD between each pair of SNPs in the window, c) remove one of a pair of SNPs if the LD is greater than 0.5, d) shift the window 5 SNPs forward and repeat the procedure. 514,246 SNPs remained after the SNP pruning and were used to run principal components analysis. As a result, the top 10 principal components were used to adjust for ancestry in a regression model. Both SNP pruning and principal components analyses were performed in PLINK1.9

The samples in the genetic data were mapped to the phenotype data of HRS to retrieve parental longevity phenotypes. We used the phenotype data from the most recent wave or visit of HRS where parental data was missing if the parent died before the last visit. To recover the data, we sequentially looked back at the last visit until the parental data was available. Father’s and mother’s data therefore could be from different visits and the participant’s age varied with phenotypes. Mothers dying before 61 years and fathers dying before 46 years were excluded from the analysis, based on the same analysis of distribution of ages of death applied in UK Biobank. Parent’s age was set to missing if he or she was still alive and younger than the child’s age +10.

The genetic risk score analysis was conducted in a Cox regression model (against combined parents' attained age, alive or dead: derived by summing the z-transformed mothers and fathers' attained ages) with adjustment for age and sex of the participant and the top 10 principal components.

Wisconsin Longitudinal Study (WLS)

In the WLS saliva samples were collected by mail in 2006-7 and during in-person interviews in 2011 for those who did not mail back samples. Genotyping on the Illumina HumanOmniExpress array was completed at Johns Hopkins’ Center for Inherited Disease Research (CIDR) and data cleaning and imputation was performed in collaboration with the Genetic Analysis Center at the University of Washington. After quality control, there are a total of 9,012 respondents and 688,596 SNPs available for analysis. The sample has 4,601 singletons and 2,263 families with more than one member genotyped. Only 66 respondents are not of European descent.

Imputation of 33,549,349 genetic variants from the 1000 Genomes Project reference panel (Phase 3) was performed using IMPUTE2, with 686,143 study SNPs serving as the imputation basis. Of the 10 variants used to construct the GRS, two were directly genotyped and eight were imputed.

Parent’s age was constructed by subtracting the year of birth from the year of death. When the latter was not available, the year of the most recent interview was used instead and the parent was classified as not dead. Mothers dying before 56 years and fathers dying before 46 years were excluded from the analysis, based on the same modelling of distributions of ages of death applied in UK Biobank. Parent’s age was also set to missing if the disagreement between sibling reports exceeded 10 years. The analysis phenotype was the standardized sum of both parents' ages and was set to missing if either element was not available.

The genetic risk score analysis was conducted in a Cox regression with censored data, adjusting for age at saliva collection, sex of the participant and the top 5 principal components. Standard errors were clustered at the family level. Limiting the sample to the primary WLS respondents had a negligible effect on the results.

Power Calculations

R package `pwr` and function ‘pwr.f2.test’ was utilized for power calculations of regression analysis. The numerator degrees of freedom were set to 1, significance level to 0.05, and power to 0.8.

Genome-Wide Association Study

We used BOLT-LMM v2.2 to model the associations between imputed variants (dosages) and each phenotype [38], which uses a linear mixed effects model approach. We looked at the results for variants with imputation quality >0.1, HWE p-values >1x10-6 and minor allele frequencies >0.1% in the white/British subset used for all analyses. For variants on the X, Y and mitochondrial chromosomes only in the directly-genotyped data we used Plink (v1.9) [39] in linear (additive) or binary (fisher) models, as appropriate, adjusted for the same covariates as above including the first 5 principal components from FlashPCA.

Pathways Analysis and Functional Annotation

We utilized FUMA (Functional Mapping and Annotation of Genome-Wide Association Studies) to investigate the genetic variants associated with parents’ attained age [16]. In particular, FUMA performs gene-set analysis (using MAGMA) to identify pathways enriched amongst the significant genes (weighted by the SNP-associations in proximity to them). In addition, FUMA searched eQTL databases to identify SNPs that significant alter the expression of genes in various tissues.

Acknowledgments

This research has been conducted using the UK Biobank Resource, under application 1417. The authors wish to thank the UK Biobank participants and coordinators for this unique dataset.

Conflicts of Interest

All authors declare no conflicts of interest.

Funding

This work was generously funded by an award to DM and LH by the Medical Research Council MR/M023095/1. LF is supported by the Intramural Research Program of the National Institute on Aging, U.S. National Institutes of Health. Input from CK and GK was supported by the University of Connecticut Health Center. The Health and Retirement Study (HRS) is a longitudinal project sponsored by the National Institute on Aging (NIA U01AG009740) and the Social Security Administration. This research uses data from the Wisconsin Longitudinal Study (WLS) of the University of Wisconsin-Madison. Since 1991, the WLS has been supported principally by the National Institute on Aging (AG09775, AG21079 and AG33285), with additional support from the Vilas Estate Trust, the National Science Foundation, the Spencer Foundation and the Graduate School of the University of Wisconsin-Madison. A public use file of data from the WLS is available from the Wisconsin Longitudinal Study, University of Wisconsin-Madison, 1180 Observatory Drive, Madison, Wisconsin 53706 and at (http://www.ssc.wisc.edu/wlsresearch/data). This material is the result of work supported with resources at the William S. Middleton Memorial Veterans Hospital, Madison, Wisconsin.

References

  • 1. Kenyon CJ. The genetics of ageing. Nature. 2010; 464:504–12. https://doi.org/10.1038/nature08980 [PubMed]
  • 2. Joshi PK, Pirastu N, Kentistou KA, Fischer K, Hofer E, Schraut KE, Clark DW, Nutile T, Barnes CL, Timmers PR, Shen X, Gandin I, McDaid AF, et al. Genome-wide meta-analysis associates HLA-DQA1/DRB1 and LPA and lifestyle factors with human longevity. Nat Commun. 2017; 8:910. https://doi.org/10.1038/s41467-017-00934-5 [PubMed]
  • 3. Broer L, Buchman AS, Deelen J, Evans DS, Faul JD, Lunetta KL, Sebastiani P, Smith JA, Smith AV, Tanaka T, Yu L, Arnold AM, Aspelund T, et al. GWAS of longevity in CHARGE consortium confirms APOE and FOXO3 candidacy. J Gerontol A Biol Sci Med Sci. 2015; 70:110–18. https://doi.org/10.1093/gerona/glu166 [PubMed]
  • 4. Deelen J, Beekman M, Uh HW, Broer L, Ayers KL, Tan Q, Kamatani Y, Bennet AM, Tamm R, Trompet S, Guðbjartsson DF, Flachsbart F, Rose G, et al. Genome-wide association meta-analysis of human longevity identifies a novel locus conferring survival beyond 90 years of age. Hum Mol Genet. 2014; 23:4420–32. https://doi.org/10.1093/hmg/ddu139 [PubMed]
  • 5. Fortney K, Dobriban E, Garagnani P, Pirazzini C, Monti D, Mari D, Atzmon G, Barzilai N, Franceschi C, Owen AB, Kim SK. Genome-Wide Scan Informed by Age-Related Disease Identifies Loci for Exceptional Human Longevity. PLoS Genet. 2015; 11:e1005728. https://doi.org/10.1371/journal.pgen.1005728 [PubMed]
  • 6. McDaid AF, Joshi PK, Porcu E, Komljenovic A, Li H, Sorrentino V, Litovchenko M, Bevers RP, Rüeger S, Reymond A, Bochud M, Deplancke B, Williams RW, et al. Bayesian association scan reveals loci associated with human lifespan and linked biomarkers. Nat Commun. 2017; 8:15842. https://doi.org/10.1038/ncomms15842 [PubMed]
  • 7. Hjelmborg J, Iachine I, Skytthe A, Vaupel JW, McGue M, Koskenvuo M, Kaprio J, Pedersen NL, Christensen K. Genetic influence on human lifespan and longevity. Hum Genet. 2006; 119:312–21. https://doi.org/10.1007/s00439-006-0144-y [PubMed]
  • 8. Gudmundsson H, Gudbjartsson DF, Frigge M, Gulcher JR, Stefánsson K. Inheritance of human longevity in Iceland. Eur J Hum Genet. 2000; 8:743–49. https://doi.org/10.1038/sj.ejhg.5200527 [PubMed]
  • 9. Dutta A, Henley W, Robine JM, Langa KM, Wallace RB, Melzer D. Longer lived parents: protective associations with cancer incidence and overall mortality. J Gerontol A Biol Sci Med Sci. 2013; 68:1409–18. https://doi.org/10.1093/gerona/glt061 [PubMed]
  • 10. Dutta A, Henley W, Robine JM, Llewellyn D, Langa KM, Wallace RB, Melzer D. Aging children of long-lived parents experience slower cognitive decline. Alzheimers Dement. 2013; 10:S315–22. https://doi.org/10.1016/j.jalz.2013.07.002 [PubMed]
  • 11. Atkins JL, Pilling LC, Ble A, Dutta A, Harries LW, Murray A, Brayne C, Robine JM, Kuchel GA, Ferrucci L, Melzer D. Longer-Lived Parents and Cardiovascular Outcomes: 8-Year Follow-Up In 186,000 U.K. Biobank Participants. J Am Coll Cardiol. 2016; 68:874–75. https://doi.org/10.1016/j.jacc.2016.05.072 [PubMed]
  • 12. Pilling LC, Atkins JL, Bowman K, Jones SE, Tyrrell J, Beaumont RN, Ruth KS, Tuke MA, Yaghootkar H, Wood AR, Freathy RM, Murray A, Weedon MN, et al. Human longevity is influenced by many genetic variants: evidence from 75,000 UK Biobank participants. Aging (Albany NY). 2016; 8:547–60. https://doi.org/10.18632/aging.100930 [PubMed]
  • 13. Joshi PK, Fischer K, Schraut KE, Campbell H, Esko T, Wilson JF. Variants near CHRNA3/5 and APOE have age- and sex-related effects on human lifespan. Nat Commun. 2016; 7:11174. https://doi.org/10.1038/ncomms11174 [PubMed]
  • 14. Welter D, MacArthur J, Morales J, Burdett T, Hall P, Junkins H, Klemm A, Flicek P, Manolio T, Hindorff L, Parkinson H. The NHGRI GWAS Catalog, a curated resource of SNP-trait associations. Nucleic Acids Res. 2014; 42:D1001–06. https://doi.org/10.1093/nar/gkt1229 [PubMed]
  • 15. Zhernakova DV, Deelen P, Vermaat M, van Iterson M, van Galen M, Arindrarto W, van ’t Hof P, Mei H, van Dijk F, Westra HJ, Bonder MJ, van Rooij J, Verkerk M, et al. Identification of context-dependent expression quantitative trait loci in whole blood. Nat Genet. 2017; 49:139–45. https://doi.org/10.1038/ng.3737 [PubMed]
  • 16. Watanabe K, Taskesen E, van Bochoven A, Posthuma D. FUMA: functional mapping and annotation of genetic associations. bioRxiv. 2017; 110023. https://doi.org/10.1101/110023
  • 17. Locke AE, Kahali B, Berndt SI, Justice AE, Pers TH, Day FR, Powell C, Vedantam S, Buchkovich ML, Yang J, Croteau-Chonka DC, Esko T, Fall T, et al, and LifeLines Cohort Study, and ADIPOGen Consortium, and AGEN-BMI Working Group, and CARDIOGRAMplusC4D Consortium, and CKDGen Consortium, and GLGC, and ICBP, and MAGIC Investigators, and MuTHER Consortium, and MIGen Consortium, and PAGE Consortium, and ReproGen Consortium, and GENIE Consortium, and International Endogene Consortium. Genetic studies of body mass index yield new insights for obesity biology. Nature. 2015; 518:197–206. https://doi.org/10.1038/nature14177 [PubMed]
  • 18. Michailidou K, Beesley J, Lindstrom S, Canisius S, Dennis J, Lush MJ, Maranian MJ, Bolla MK, Wang Q, Shah M, Perkins BJ, Czene K, Eriksson M, et al, and BOCS, and kConFab Investigators, and AOCS Group, and NBCS, and GENICA Network. Genome-wide association analysis of more than 120,000 individuals identifies 15 new susceptibility loci for breast cancer. Nat Genet. 2015; 47:373–80. https://doi.org/10.1038/ng.3242 [PubMed]
  • 19. Okbay A, Beauchamp JP, Fontana MA, Lee JJ, Pers TH, Rietveld CA, Turley P, Chen GB, Emilsson V, Meddens SF, Oskarsson S, Pickrell JK, Thom K, et al, and LifeLines Cohort Study. Genome-wide association study identifies 74 loci associated with educational attainment. Nature. 2016; 533:539–42. https://doi.org/10.1038/nature17671 [PubMed]
  • 20. Pattaro C, Teumer A, Gorski M, Chu AY, Li M, Mijatovic V, Garnaas M, Tin A, Sorice R, Li Y, Taliun D, Olden M, Foster M, et al, and ICBP Consortium, and AGEN Consortium, and CARDIOGRAM, and CHARGe-Heart Failure Group, and ECHOGen Consortium. Genetic associations at 53 loci highlight cell types and biological pathways relevant for kidney function. Nat Commun. 2016; 7:10023. https://doi.org/10.1038/ncomms10023 [PubMed]
  • 21. Castelo-Branco C, Soveral I. The immune system and aging: a review. Gynecol Endocrinol. 2014; 30:16–22. https://doi.org/10.3109/09513590.2013.852531 [PubMed]
  • 22. Aguilo F, Di Cecilia S, Walsh MJ. Long Non-coding RNA ANRIL and Polycomb in Human Cancers and Cardiovascular Disease. Curr Top Microbiol Immunol. 2016; 394:29–39. https://doi.org/10.1007/82_2015_455 [PubMed]
  • 23. Baker DJ, Wijshake T, Tchkonia T, LeBrasseur NK, Childs BG, van de Sluis B, Kirkland JL, van Deursen JM. Clearance of p16Ink4a-positive senescent cells delays ageing-associated disorders. Nature. 2011; 479:232–36. https://doi.org/10.1038/nature10600 [PubMed]
  • 24. Bochenek G, Häsler R, El Mokhtari NE, König IR, Loos BG, Jepsen S, Rosenstiel P, Schreiber S, Schaefer AS. The large non-coding RNA ANRIL, which is associated with atherosclerosis, periodontitis and several forms of cancer, regulates ADIPOR1, VAMP3 and C11ORF10. Hum Mol Genet. 2013; 22:4516–27. https://doi.org/10.1093/hmg/ddt299 [PubMed]
  • 25. Congrains A, Kamide K, Oguro R, Yasuda O, Miyata K, Yamamoto E, Kawai T, Kusunoki H, Yamamoto H, Takeya Y, Yamamoto K, Onishi M, Sugimoto K, et al. Genetic variants at the 9p21 locus contribute to atherosclerosis through modulation of ANRIL and CDKN2A/B. Atherosclerosis. 2012; 220:449–55. https://doi.org/10.1016/j.atherosclerosis.2011.11.017 [PubMed]
  • 26. Dale BL, Madhur MS. Linking inflammation and hypertension via LNK/SH2B3. Curr Opin Nephrol Hypertens. 2016; 25:87–93. https://doi.org/10.1097/MNH.0000000000000196 [PubMed]
  • 27. Sato K, Emi M, Ezura Y, Fujita Y, Takada D, Ishigami T, Umemura S, Xin Y, Wu LL, Larrinaga-Shum S, Stephenson SH, Hunt SC, Hopkins PN. Soluble epoxide hydrolase variant (Glu287Arg) modifies plasma total cholesterol and triglyceride phenotype in familial hypercholesterolemia: intrafamilial association study in an eight-generation hyperlipidemic kindred. J Hum Genet. 2004; 49:29–34. https://doi.org/10.1007/s10038-003-0103-6 [PubMed]
  • 28. Wei Q, Doris PA, Pollizotto MV, Boerwinkle E, Jacobs DRJr, Siscovick DS, Fornage M. Sequence variation in the soluble epoxide hydrolase gene and subclinical coronary atherosclerosis: interaction with cigarette smoking. Atherosclerosis. 2007; 190:26–34. https://doi.org/10.1016/j.atherosclerosis.2006.02.021 [PubMed]
  • 29. Barrell D, Dimmer E, Huntley RP, Binns D, O’Donovan C, Apweiler R. The GOA database in 2009--an integrated Gene Ontology Annotation resource. Nucleic Acids Res. 2009; 37:D396–403. https://doi.org/10.1093/nar/gkn803 [PubMed]
  • 30. Pesu M, Watford WT, Wei L, Xu L, Fuss I, Strober W, Andersson J, Shevach EM, Quezado M, Bouladoux N, Roebroek A, Belkaid Y, Creemers J, O’Shea JJ. T-cell-expressed proprotein convertase furin is essential for maintenance of peripheral immune tolerance. Nature. 2008; 455:246–50. https://doi.org/10.1038/nature07210 [PubMed]
  • 31. Kennedy BK, Berger SL, Brunet A, Campisi J, Cuervo AM, Epel ES, Franceschi C, Lithgow GJ, Morimoto RI, Pessin JE, Rando TA, Richardson A, Schadt EE, et al. Geroscience: linking aging to chronic disease. Cell. 2014; 159:709–13. https://doi.org/10.1016/j.cell.2014.10.039 [PubMed]
  • 32. Meimaridou E, Hughes CR, Kowalczyk J, Chan LF, Clark AJ, Metherell LA. ACTH resistance: genes and mechanisms. Endocr Dev. 2013; 24:57–66. https://doi.org/10.1159/000342504 [PubMed]
  • 33. Fry A, Littlejohns TJ, Sudlow C, Doherty N, Adamska L, Sprosen T, Collins R, Allen NE. Comparison of Sociodemographic and Health-Related Characteristics of UK Biobank Participants With Those of the General Population. Am J Epidemiol. 2017; 186:1026–34. https://doi.org/10.1093/aje/kwx246 [PubMed]
  • 34. Sebastiani P, Bae H, Gurinovich A, Soerensen M, Puca A, Perls TT. Limitations and risks of meta-analyses of longevity studies. Mech Ageing Dev. 2017; 165:139–46. https://doi.org/10.1016/j.mad.2017.01.008 [PubMed]
  • 35. Locke JM, Hysenaj G, Wood AR, Weedon MN, Harries LW. Targeted allelic expression profiling in human islets identifies cis-regulatory effects for multiple variants identified by type 2 diabetes genome-wide association studies. Diabetes. 2015; 64:1484–91. https://doi.org/10.2337/db14-0957 [PubMed]
  • 36. UK Biobank. UK Biobank genotype collections and QC on 150k participants. 2014. .
  • 37. Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, Motyer A, Vukcevic D, Delaneau O, O'Connell J, Cortes A, Welsh S, McVean G, et al. Genome-wide genetic data on ~500,000 UK Biobank participants. bioRxiv. 2017. https://doi.org/10.1101/166298
  • 38. Loh PR, Tucker G, Bulik-Sullivan BK, Vilhjálmsson BJ, Finucane HK, Salem RM, Chasman DI, Ridker PM, Neale BM, Berger B, Patterson N, Price AL. Efficient Bayesian mixed-model analysis increases association power in large cohorts. Nat Genet. 2015; 47:284–90. https://doi.org/10.1038/ng.3190 [PubMed]
  • 39. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007; 81:559–75. https://doi.org/10.1086/519795 [PubMed]
  • 40. Nikpay M, Goel A, Won HH, Hall LM, Willenborg C, Kanoni S, Saleheen D, Kyriakou T, Nelson CP, Hopewell JC, Webb TR, Zeng L, Dehghan A, et al. A comprehensive 1,000 Genomes-based genome-wide association meta-analysis of coronary artery disease. Nat Genet. 2015; 47:1121–30. https://doi.org/10.1038/ng.3396 [PubMed]
  • 41. Astle WJ, Elding H, Jiang T, Allen D, Ruklisa D, Mann AL, Mead D, Bouman H, Riveros-Mckay F, Kostadima MA, Lambourne JJ, Sivapalaratnam S, Downes K, et al. The Allelic Landscape of Human Blood Cell Trait Variation and Links to Common Complex Disease. Cell. 2016; 167:1415–1429.e19. https://doi.org/10.1016/j.cell.2016.10.042 [PubMed]
  • 42. Kettunen J, Demirkan A, Würtz P, Draisma HH, Haller T, Rawal R, Vaarhorst A, Kangas AJ, Lyytikäinen LP, Pirinen M, Pool R, Sarin AP, Soininen P, et al. Genome-wide study for circulating metabolites identifies 62 loci and reveals novel systemic effects of LPA. Nat Commun. 2016; 7:11122. https://doi.org/10.1038/ncomms11122 [PubMed]
  • 43. Wood AR, Esko T, Yang J, Vedantam S, Pers TH, Gustafsson S, Chu AY, Estrada K, Luan J, Kutalik Z, Amin N, Buchkovich ML, Croteau-Chonka DC, et al, and Electronic Medical Records and Genomics (eMEMERGEGE) Consortium, and MIGen Consortium, and PAGEGE Consortium, and LifeLines Cohort Study. Defining the role of common variation in the genomic and biological architecture of adult human height. Nat Genet. 2014; 46:1173–86. https://doi.org/10.1038/ng.3097 [PubMed]
  • 44. Ehret GB, Munroe PB, Rice KM, Bochud M, Johnson AD, Chasman DI, Smith AV, Tobin MD, Verwoert GC, Hwang SJ, Pihur V, Vollenweider P, O’Reilly PF, et al, and CHARGE-HF consortium. Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk. Nature. 2011; 478:103–09. https://doi.org/10.1038/nature10405 [PubMed]
  • 45. Willer CJ, Schmidt EM, Sengupta S, Peloso GM, Gustafsson S, Kanoni S, Ganna A, Chen J, Buchkovich ML, Mora S, Beckmann JS, Bragg-Gresham JL, Chang HY, et al, and Global Lipids Genetics Consortium. Discovery and refinement of loci associated with lipid levels. Nat Genet. 2013; 45:1274–83. https://doi.org/10.1038/ng.2797 [PubMed]