Evaluation and extension of a kernel-based method for gene-gene interaction tests of common variants

Interaction is likely to play a significant role in complex diseases, and various methods are available for identifying interactions between variants in genome-wide association studies (GWAS). Kernel-based variance component methods such as SKAT are flexible and computationally efficient methods for identifying marginal associations. A kernel-based variance component method, called the Gene-centric Gene-Gene Interaction with Smoothing-sPline ANOVA model (SPA3G) was proposed to identify gene-gene interactions for a quantitative trait. For interaction testing, the SPA3G method performs better than some SNP-based approaches under many scenarios. In this thesis, we evaluate the properties of the SPA3G method and extend SPA3G using alternative p-value approximations and interaction kernels. This thesis focuses on common variants only. Our simulation results show that the allele matching interaction kernel, combined with the method of moments p-value approximation, leads to inflated type I error in small samples. For small samples, we propose a Principal Component (PC)-based interaction kernel and computing p-values with a 3-moment adjustment that yield more appropriate type I error. We also propose a weighted PC kernel that has higher power than competing approaches when interaction effects are sparse. By combining the two proposed kernels, we develop omnibus methods that obtain near-optimal power in most settings. Finally, we illustrate how to analyze the interaction between selected gene pairs on the age at natural menopause (ANM) from the Framingham Heart Study

Evaluation and extension of a kernel-based method for gene-gene interaction tests of common variants

Interaction is likely to play a significant role in complex diseases, and various methods are available for identifying interactions between variants in genome-wide association studies (GWAS). Kernel-based variance component methods such as SKAT are flexible and computationally efficient methods for identifying marginal associations. A kernel-based variance component method, called the Gene-centric Gene-Gene Interaction with Smoothing-sPline ANOVA model (SPA3G) was proposed to identify gene-gene interactions for a quantitative trait. For interaction testing, the SPA3G method performs better than some SNP-based approaches under many scenarios. In this thesis, we evaluate the properties of the SPA3G method and extend SPA3G using alternative p-value approximations and interaction kernels. This thesis focuses on common variants only. Our simulation results show that the allele matching interaction kernel, combined with the method of moments p-value approximation, leads to inflated type I error in small samples. For small samples, we propose a Principal Component (PC)-based interaction kernel and computing p-values with a 3-moment adjustment that yield more appropriate type I error. We also propose a weighted PC kernel that has higher power than competing approaches when interaction effects are sparse. By combining the two proposed kernels, we develop omnibus methods that obtain near-optimal power in most settings. Finally, we illustrate how to analyze the interaction between selected gene pairs on the age at natural menopause (ANM) from the Framingham Heart Study

Human longevity is influenced by many genetic variants: evidence from 75,000 UK Biobank participants

This is the final version of the article. Available from the publisher via the DOI in this record.Variation in human lifespan is 20 to 30% heritable in twins but few genetic variants have been identified. We undertook a Genome Wide Association Study (GWAS) using age at death of parents of middle-aged UK Biobank participants of European decent (n=75,244 with father's and/or mother's data, excluding early deaths). Genetic risk scores for 19 phenotypes (n=777 proven variants) were also tested. In GWAS, a nicotine receptor locus(CHRNA3, previously associated with increased smoking and lung cancer) was associated with fathers' survival. Less common variants requiring further confirmation were also identified. Offspring of longer lived parents had more protective alleles for coronary artery disease, systolic blood pressure, body mass index, cholesterol and triglyceride levels, type-1 diabetes, inflammatory bowel disease and Alzheimer's disease. In candidate analyses, variants in the TOMM40/APOE locus were associated with longevity, but FOXO variants were not. Associations between extreme longevity (mother >=98 years, fathers >=95 years, n=1,339) and disease alleles were similar, with an additional association with HDL cholesterol (p=5.7x10-3). These results support a multiple protective factors model influencing lifespan and longevity (top 1% survival) in humans, with prominent roles for cardiovascular-related pathways. Several of these genetically influenced risks, including blood pressure and tobacco exposure, are potentially modifiable.This work was generously funded by an award to DM, TF, AM, LH and CB by the Medical Research Council MR/M023095/1. 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. S.E.J. is funded by the Medical Research Council (grant: MR/M005070/1). J.T. is funded by a Diabetes Research and Wellness Foundation Fellowship. R.B. is funded by the Wellcome Trust and Royal Society grant: 104150/Z/14/Z. M.A.T., M.N.W. and A.M. are supported by the Wellcome Trust Institutional Strategic Support Award (WT097835MF). R.M.F. is a Sir Henry Dale Fellow (Wellcome Trust and Royal Society grant: 104150/Z/14/Z). A.R.W. H.Y., and T.M.F. are supported by the European Research Council grant: 323195:GLUCOSEGENES-FP7-IDEAS-ERC. The funders had no influence on study design, data collection and analysis, decision to publish, or preparation of the manuscript. The Framingham Heart Study is supported by Contract No. N01-HC-25195 and HHSN268201500001I and its contract with Affymetrix, Inc for genotyping services (Contract No. N02-HL-6-4278). The phenotypegenotype association analyses were supported by National Institute of Aging R01AG29451. This work has made use of the resources provided by the University of Exeter Science Strategy and resulting Systems Biology initiative. Primarily these include high-performance computing facilities managed by Konrad Paszkiewicz of the College of Environmental and Life Sciences and Pete Leggett of the University of Exeter Academics services unit

GWAS of epigenetic aging rates in blood reveals a critical role for TERT.

DNA methylation age is an accurate biomarker of chronological age and predicts lifespan, but its underlying molecular mechanisms are unknown. In this genome-wide association study of 9907 individuals, we find gene variants mapping to five loci associated with intrinsic epigenetic age acceleration (IEAA) and gene variants in three loci associated with extrinsic epigenetic age acceleration (EEAA). Mendelian randomization analysis suggests causal influences of menarche and menopause on IEAA and lipoproteins on IEAA and EEAA. Variants associated with longer leukocyte telomere length (LTL) in the telomerase reverse transcriptase gene (TERT) paradoxically confer higher IEAA (P < 2.7 × 10-11). Causal modeling indicates TERT-specific and independent effects on LTL and IEAA. Experimental hTERT-expression in primary human fibroblasts engenders a linear increase in DNA methylation age with cell population doubling number. Together, these findings indicate a critical role for hTERT in regulating the epigenetic clock, in addition to its established role of compensating for cell replication-dependent telomere shortening

Genome-wide association meta-analysis of fish and EPA plus DHA consumption in 17 US and European cohorts

Background Regular fish and omega-3 consumption may have several health benefits and are recommended by major dietary guidelines. Yet, their intakes remain remarkably variable both within and across populations, which could partly owe to genetic influences. Objective To identify common genetic variants that influence fish and dietary eicosapentaenoic acid plus docosahexaenoic acid (EPA+DHA) consumption. Design We conducted genome-wide association (GWA) meta-analysis of fish (n = 86,467) and EPA+DHA (n = 62,265) consumption in 17 cohorts of European descent from the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) Consortium Nutrition Working Group. Results from cohort-specific GWA analyses (additive model) for fish and EPA+DHA consumption were adjusted for age, sex, energy intake, and population stratification, and meta-analyzed separately using fixed-effect meta-analysis with inverse variance weights (METAL software). Additionally, heritability was estimated in 2 cohorts. Results Heritability estimates for fish and EPA+DHA consumption ranged from 0.13-0.24 and 0.12-0.22, respectively. A significant GWA for fish intake was observed for rs9502823 on chromosome 6: each copy of the minor allele (Freq(A) = 0.015) was associated with 0.029 servings/day (similar to 1 serving/month) lower fish consumption (P = 1.96x10(-8)). No significant association was observed for EPA+DHA, although rs7206790 in the obesity-associated FTO gene was among top hits (P = 8.18x10(-7)). Post-hoc calculations demonstrated 95% statistical power to detect a genetic variant associated with effect size of 0.05% for fish and 0.08% for EPA+DHA. Conclusions These novel findings suggest that non-genetic personal and environmental factors are principal determinants of the remarkable variation in fish consumption, representing modifiable targets for increasing intakes among all individuals. Genes underlying the signal at rs72838923 and mechanisms for the association warrant further investigation.Peer reviewe

Genome-wide association meta-analysis of fish and EPA+DHA consumption in 17 US and European cohorts

Background: Regular fish and omega-3 consumption may have several health benefits and are recommended by major dietary guidelines. Yet, their intakes remain remarkably variable both within and across populations, which could partly owe to genetic influences. Objective: To identify common genetic variants that influence fish and dietary eicosapentaenoic acid plus docosahexaenoic acid (EPA+DHA) consumption. Design: We conducted genome-wide association (GWA) meta-analysis of fish (n = 86, 467) and EPA +DHA (n = 62, 265) consumption in 17 cohorts of European descent from the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) Consortium Nutrition Working Group. Results from cohort-specific GWA analyses (additive model) for fish and EPA+DHA consumption were adjusted for age, sex, energy intake, and population stratification, and meta-analyzed separately using fixed-effect meta-analysis with inverse variance weights (METAL software). Additionally, heritability was estimated in 2 cohorts. Results: Heritability estimates for fish and EPA+DHA consumption ranged from 0.13

Genome-Wide Association of Body Fat Distribution in African Ancestry Populations Suggests New Loci

Central obesity, measured by waist circumference (WC) or waist-hip ratio (WHR), is a marker of body fat distribution. Although obesity disproportionately affects minority populations, few studies have conducted genome-wide association study (GWAS) of fat distribution among those of predominantly African ancestry (AA). We performed GWAS of WC and WHR, adjusted and unadjusted for BMI, in up to 33,591 and 27,350 AA individuals, respectively. We identified loci associated with fat distribution in AA individuals using meta-analyses of GWA results for WC and WHR (stage 1). Overall, 25 SNPs with single genomic control (GC)-corrected p-values<5.0×10−6 were followed-up (stage 2) in AA with WC and with WHR. Additionally, we interrogated genomic regions of previously identified European ancestry (EA) WHR loci among AA. In joint analysis of association results including both Stage 1 and 2 cohorts, 2 SNPs demonstrated association, rs2075064 at LHX2, p = 2.24×10−8 for WC-adjusted-for-BMI, and rs6931262 at RREB1, p = 2.48×10−8 for WHR-adjusted-for-BMI. However, neither signal was genome-wide significant after double GC-correction (LHX2: p = 6.5×10−8; RREB1: p = 5.7×10−8). Six of fourteen previously reported loci for waist in EA populations were significant (p<0.05 divided by the number of independent SNPs within the region) in AA studied here (TBX15-WARS2, GRB14, ADAMTS9, LY86, RSPO3, ITPR2-SSPN). Further, we observed associations with metabolic traits: rs13389219 at GRB14 associated with HDL-cholesterol, triglycerides, and fasting insulin, and rs13060013 at ADAMTS9 with HDL-cholesterol and fasting insulin. Finally, we observed nominal evidence for sexual dimorphism, with stronger results in AA women at the GRB14 locus (p for interaction = 0.02). In conclusion, we identified two suggestive loci associated with fat distribution in AA populations in addition to confirming 6 loci previously identified in populations of EA. These findings reinforce the concept that there are fat distribution loci that are independent of generalized adiposity

Genome-Wide Association of Body Fat Distribution in African Ancestry Populations Suggests New Loci

Central obesity, measured by waist circumference (WC) or waist-hip ratio (WHR), is a marker of body fat distribution. Although obesity disproportionately affects minority populations, few studies have conducted genome-wide association study (GWAS) of fat distribution among those of predominantly African ancestry (AA). We performed GWAS of WC and WHR, adjusted and unadjusted for BMI, in up to 33,591 and 27,350 AA individuals, respectively. We identified loci associated with fat distribution in AA individuals using meta-analyses of GWA results for WC and WHR (stage 1). Overall, 25 SNPs with single genomic control (GC)-corrected p-values<5.0×10−6 were followed-up (stage 2) in AA with WC and with WHR. Additionally, we interrogated genomic regions of previously identified European ancestry (EA) WHR loci among AA. In joint analysis of association results including both Stage 1 and 2 cohorts, 2 SNPs demonstrated association, rs2075064 at LHX2, p = 2.24×10−8 for WC-adjusted-for-BMI, and rs6931262 at RREB1, p = 2.48×10−8 for WHR-adjusted-for-BMI. However, neither signal was genome-wide significant after double GC-correction (LHX2: p = 6.5×10−8; RREB1: p = 5.7×10−8). Six of fourteen previously reported loci for waist in EA populations were significant (p<0.05 divided by the number of independent SNPs within the region) in AA studied here (TBX15-WARS2, GRB14, ADAMTS9, LY86, RSPO3, ITPR2-SSPN). Further, we observed associations with metabolic traits: rs13389219 at GRB14 associated with HDL-cholesterol, triglycerides, and fasting insulin, and rs13060013 at ADAMTS9 with HDL-cholesterol and fasting insulin. Finally, we observed nominal evidence for sexual dimorphism, with stronger results in AA women at the GRB14 locus (p for interaction = 0.02). In conclusion, we identified two suggestive loci associated with fat distribution in AA populations in addition to confirming 6 loci previously identified in populations of EA. These findings reinforce the concept that there are fat distribution loci that are independent of generalized adiposity

Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.Peer reviewe