ProbABEL package for genome-wide association analysis of imputed data

Common genetic variants do not associate with CAD in familial hypercholesterolemia

In recent years, multiple loci dispersed on the genome have been shown to be associated with coronary artery disease (CAD). We investigated whether these common genetic variants also hold value for CAD prediction in a large cohort of patients with familial hypercholesterolemia (FH). We genotyped a total of 41 single-nucleotide polymorphisms (SNPs) in 1701 FH patients, of whom 482 patients (28.3%) had at least one coronary event during an average follow up of 66 years. The association of each SNP with event-free survival time was calculated with a Cox proportional hazard model. In the cardiovascular disease risk factor adjusted analysis, the most significant SNP was rs1122608:G>T in the SMARCA4 gene near the LDL-receptor (LDLR) gene, with a hazard ratio for CAD risk of 0.74 (95% CI 0.49-0.99; P-value 0.021). However, none of the SNPs reached the Bonferroni threshold. Of all the known CAD loci analyzed, the SMARCA4 locus near the LDLR had the strongest negative association with CAD in this high-risk FH cohort. The effect is contrary to what was expected. None of the other loci showed association with CAD.

Global analysis of DNA methylation variation in adipose tissue from twins reveals links to disease-associated variants in distal regulatory elements

Epigenetic modifications such as DNA methylation play a key role in gene regulation and disease susceptibility. However, little is known about the genome-wide frequency, localization, and function of methylation variation and how it is regulated by genetic and environmental factors. We utilized the Multiple Tissue Human Expression Resource (MuTHER) and generated Illumina 450K adipose methylome data from 648 twins. We found that individual CpGs had low variance and that variability was suppressed in promoters. We noted that DNA methylation variation was highly heritable (h(2)median = 0.34) and that shared environmental effects correlated with metabolic phenotype-associated CpGs. Analysis of methylation quantitative-trait loci (metQTL) revealed that 28% of CpGs were associated with nearby SNPs, and when overlapping them with adipose expression quantitative-trait loci (eQTL) from the same individuals, we found that 6% of the loci played a role in regulating both gene expression and DNA methylation. These associations were bidirectional, but there were pronounced negative associations for promoter CpGs. Integration of metQTL with adipose reference epigenomes and disease associations revealed significant enrichment of metQTL overlapping metabolic-trait or disease loci in enhancers (the strongest effects were for high-density lipoprotein cholesterol and body mass index [BMI]). We followed up with the BMI SNP rs713586, a cg01884057 metQTL that overlaps an enhancer upstream of ADCY3, and used bisulphite sequencing to refine this region. Our results showed widespread population invariability yet sequence dependence on adipose DNA methylation but that incorporating maps of regulatory elements aid in linking CpG variation to gene regulation and disease risk in a tissue-dependent manner.

Gene-environment interactions in cancer epidemiology: a National Cancer Institute Think Tank report

Cancer risk is determined by a complex interplay of genetic and environmental factors. Genome-wide association studies (GWAS) have identified hundreds of common (minor allele frequency [MAF] > 0.05) and less common (0.01 < MAF < 0.05) genetic variants associated with cancer. The marginal effects of most of these variants have been small (odds ratios: 1.1-1.4). There remain unanswered questions on how best to incorporate the joint effects of genes and environment, including gene-environment (G × E) interactions, into epidemiologic studies of cancer. To help address these questions, and to better inform research priorities and allocation of resources, the National Cancer Institute sponsored a "Gene-Environment Think Tank" on January 10-11, 2012. The objective of the Think Tank was to facilitate discussions on (1) the state of the science, (2) the goals of G × E interaction studies in cancer epidemiology, and (3) opportunities for developing novel study designs and analysis tools. This report summarizes the Think Tank discussion, with a focus on contemporary approaches to the analysis of G × E interactions. Selecting the appropriate methods requires first identifying the relevant scientific question and rationale, with an important distinction made between analyses aiming to characterize the joint effects of putative or established genetic and environmental factors and analyses aiming to discover novel risk factors or novel interaction effects. Other discussion items include measurement error, statistical power, significance, and replication. Additional designs, exposure assessments, and analytical approaches need to be considered as we move from the current small number of success stories to a fuller understanding of the interplay of genetic and environmental factors.

Genome-wide association study on dimethylarginines reveals novel AGXT2 variants associated with heart rate variability but not with overall mortality

AIMS

The purpose of this study was to identify novel genetic variants influencing circulating asymmetric dimethylarginine (ADMA) and symmetric dimethylarginine (SDMA) levels and to evaluate whether they have a prognostic value on cardiovascular mortality.

METHODS AND RESULTS

We conducted a genome-wide association study on the methylarginine traits and investigated the predictive value of the new discovered variants on mortality. Our meta-analyses replicated the previously known locus for ADMA levels in DDAH1 (rs997251; P = 1.4 × 10(-40)), identified two non-synomyous polymorphisms for SDMA levels in AGXT2 (rs37369; P = 1.4 × 10(-40) and rs16899974; P = 1.5 × 10(-38)) and one in SLC25A45 (rs34400381; P = 2.5 × 10(-10)). We also fine-mapped the AGXT2 locus for further independent association signals. The two non-synonymous AGXT2 variants independently associated with SDMA levels were also significantly related with short-term heart rate variability (HRV) indices in young adults. The major allele (C) of the novel non-synonymous rs16899974 (V498L) variant associated with decreased SDMA levels and an increase in the ratio between the low- and high-frequency spectral components of HRV (P = 0.00047). Furthermore, the SDMA decreasing allele (G) of the non-synomyous SLC25A45 (R285C) variant was associated with a lower resting mean heart rate during the HRV measurements (P = 0.0046), but not with the HRV indices. None of the studied genome-wide significant variants had any major effect on cardiovascular or total mortality in patients referred for coronary angiography.

CONCLUSIONS

AGXT2 has an important role in SDMA metabolism in humans. AGXT2 may additionally have an unanticipated role in the autonomic nervous system regulation of cardiac function.

Genome-wide association study of autistic-like traits in a general population study of young adults

Lay abstract: It has been proposed that autistic-like traits in the general population lie on a continuum, with clinical Autism Spectrum Disorder (ASD), representing the extreme end of this distribution. The current study undertook a genome-wide association (GWA) scan of 965 young Western Australian adults to identify novel risk variants associated with autistic-like traits. No associations reached genome-wide significance; however, a review of nominally associated single nucleotide polymorphisms (SNPs) indicated two positional candidate loci that have been previously implicated in autistic-like trait etiology.

Scientific abstract: Research has proposed that autistic-like traits in the general population lie on a continuum, with clinical ASD representing the extreme end of this distribution. Inherent in this proposal is that biological mechanisms associated with clinical ASD may also underpin variation in autistic-like traits within the general population. A GWA study using 2,462,046 SNPs was undertaken for ASD in 965 individuals from the Western Australian Pregnancy Cohort (Raine) Study. No SNP associations reached genome-wide significance (p < 5.0 × 10⁻⁸). However, investigations into nominal observed SNP associations (p < 1.0 × 10⁻⁵) add support to two positional candidate genes previously implicated in ASD etiology, PRKCB1, and CBLN1. The rs198198 SNP (p = 9.587 × 10⁻⁶), is located within an intron of the protein kinase C, beta 1 (PRKCB1) gene on chromosome 16p11. The PRKCB1 gene has been previously reported in linkage and association studies for ASD, and its mRNA expression has been shown to be significantly down regulated in ASD cases compared with controls. The rs16946931 SNP (p = 1.78 × 10⁻⁶) is located in a region flanking the Cerebellin 1 (CBLN1) gene on chromosome 16q12.1. The CBLN1 gene is involved with synaptogenesis and is part of a gene family previously implicated in ASD. This GWA study is only the second to examine SNPs associated with autistic-like traits in the general population, and provides evidence to support roles for the PRKCB1 and CBLN1 genes in risk of clinical ASD.

Genome-wide association study on serum alkaline phosphatase levels in a Chinese population

BACKGROUND

Serum alkaline phosphatase (ALP) is a complex phenotype influenced by both genetic and environmental factors. Recent Genome-Wide Association Studies (GWAS) have identified several loci affecting ALP levels; however, such studies in Chinese populations are limited. We performed a GWAS analyzing the association between 658,288 autosomal SNPs and serum ALP in 1,461 subjects, and replicated the top SNPs in an additional 8,830 healthy Chinese Han individuals. The interactions between significant locus and environmental factors on serum ALP levels were further investigated.

RESULTS

The association between ABO locus and serum ALP levels was replicated (P = 2.50 × 10⁻²¹, 1.12 × 10⁻⁵⁶ and 2.82 × 10⁻²⁷ for SNP rs8176720, rs651007 and rs7025162 on ABO locus, respectively). SNP rs651007 accounted for 2.15% of the total variance of serum ALP levels independently of the other 2 SNPs. When comparing our findings with previously published studies, ethnic differences were observed across populations. A significant interaction between ABO rs651007 and overweight and obesity was observed (FDR for interaction was 0.036); for individuals with GG genotype, those with normal weight and those who were overweight or obese have similar serum ALP concentrations; minor allele A of rs651007 remarkably reduced serum ALP levels, but this effect was attenuated in overweight and obese individuals.

CONCLUSIONS

Our findings indicate that ABO locus is a major determinant for serum ALP levels in Chinese Han population. Overweight and obesity modifies the effect of ABO locus on serum ALP concentrations.

Genome-Wide Association Study Identifies 3 Genomic Loci Significantly Associated With Serum Levels of Homoarginine

Background—

Low serum levels of the amino acid derivative, homoarginine, have been associated with increased risk of total and cardiovascular mortality. Homoarginine deficiency may be related to renal and heart diseases, but the pathophysiologic role of homoarginine and the genetic regulation of its serum levels are largely unknown.

Methods and Results—

In 3041 patients of the Ludwigshafen Risk and Cardiovascular Health (LURIC) study referred for coronary angiography and 2102 participants of the Young Finns Study (YFS), we performed a genome-wide association study to identify genomic loci associated with homoarginine serum levels and tested for associations of identified single-nucleotide polymorphisms with mortality in LURIC. We found genome-wide significant associations with homoarginine serum levels on chromosome 2 at the carbamoyl phosphate synthetase I locus, on chromosome 5 at the alanine-glyoxylate aminotransferase 2 locus, and on chromosome 15 at the glycine amidinotransferase locus, as well as a suggestive association on chromosome 6 at the Homo sapiens mediator complex subunit 23 gene/arginase I locus. All loci harbor enzymes located in the mitochondrium are involved in arginine metabolism. The strongest association was observed for rs1153858 at the glycine amidinotransferase locus with a P value of 1.25E-45 in the combined analysis and has been replicated in both the Die Deutsche Diabetes Dialyse Studie (4D study) and the Graz Endocrine Causes of Hypertension (GECOH) study.

Conclusions—

In our genome-wide association study, we identified 3 chromosomal regions significantly associated with serum homoarginine and another region with suggestive association, providing novel insights into the genetic regulation of homoarginine.

Genome-wide meta-analysis of systolic blood pressure in children with sickle cell disease

In pediatric sickle cell disease (SCD) patients, it has been reported that higher systolic blood pressure (SBP) is associated with increased risk of a silent cerebral infarction (SCI). SCI is a major cause of neurologic morbidity in children with SCD, and blood pressure is a potential modulator of clinical manifestations of SCD; however, the risk factors underlying these complications are not well characterized. The aim of this study was to identify genetic variants that influence SBP in an African American population in the setting of SCD, and explore the use of SBP as an endo-phenotype for SCI. We conducted a genome-wide meta-analysis for SBP using two SCD cohorts, as well as a candidate screen based on published SBP loci. A total of 1,617 patients were analyzed, and while no SNP reached genome-wide significance (P-value<5.0 x 10(-8)), a number of suggestive candidate loci were identified. The most significant SNP, rs7952106 (P-value=8.57 x 10(-7)), was in the DRD2 locus on chromosome 11. In a gene-based association analysis, MIR4301 (micro-RNA4301), which resides in an intron of DRD2, was the most significant gene (P-value=5.2 x 10(-5)). Examining 27 of the previously reported SBP associated SNPs, 4 SNPs were nominally significant. A genetic risk score was constructed to assess the aggregated genetic effect of the published SBP variants, demonstrating a significant association (P=0.05). In addition, we also assessed whether these variants are associated with SCI, validating the use of SBP as an endo-phenotype for SCI. Three SNPs were nominally associated, and only rs2357790 (5' CACNB2) was significant for both SBP and SCI. None of these SNPs retained significance after Bonferroni correction. Taken together, our results suggest the importance of DRD2 genetic variation in the modulation of SBP, and extend the aggregated importance of previously reported SNPs in the modulation of SBP in an African American cohort, more specifically in children with SCD.

Association between a genetic variant related to glutamic acid metabolism and coronary heart disease in individuals with type 2 diabetes

IMPORTANCE

Diabetes is associated with an elevated risk of coronary heart disease (CHD). Previous studies have suggested that the genetic factors predisposing to excess cardiovascular risk may be different in diabetic and nondiabetic individuals.

OBJECTIVE

To identify genetic determinants of CHD that are specific to patients with diabetes.

DESIGN, SETTING, AND PARTICIPANTS

We studied 5 independent sets of CHD cases and CHD-negative controls from the Nurses' Health Study (enrolled in 1976 and followed up through 2008), Health Professionals Follow-up Study (enrolled in 1986 and followed up through 2008), Joslin Heart Study (enrolled in 2001-2008), Gargano Heart Study (enrolled in 2001-2008), and Catanzaro Study (enrolled in 2004-2010). Included were a total of 1517 CHD cases and 2671 CHD-negative controls, all with type 2 diabetes. Results in diabetic patients were compared with those in 737 nondiabetic CHD cases and 1637 nondiabetic CHD-negative controls from the Nurses' Health Study and Health Professionals Follow-up Study cohorts. Exposures included 2,543,016 common genetic variants occurring throughout the genome.

MAIN OUTCOMES AND MEASURES

Coronary heart disease--defined as fatal or nonfatal myocardial infarction, coronary artery bypass grafting, percutaneous transluminal coronary angioplasty, or angiographic evidence of significant stenosis of the coronary arteries.

RESULTS

A variant on chromosome 1q25 (rs10911021) was consistently associated with CHD risk among diabetic participants, with risk allele frequencies of 0.733 in cases vs 0.679 in controls (odds ratio, 1.36 [95% CI, 1.22-1.51]; P = 2 × 10(-8)). No association between this variant and CHD was detected among nondiabetic participants, with risk allele frequencies of 0.697 in cases vs 0.696 in controls (odds ratio, 0.99 [95% CI, 0.87-1.13]; P = .89), consistent with a significant gene × diabetes interaction on CHD risk (P = 2 × 10(-4)). Compared with protective allele homozygotes, rs10911021 risk allele homozygotes were characterized by a 32% decrease in the expression of the neighboring glutamate-ammonia ligase (GLUL) gene in human endothelial cells (P = .0048). A decreased ratio between plasma levels of γ-glutamyl cycle intermediates pyroglutamic and glutamic acid was also shown in risk allele homozygotes (P = .029).

CONCLUSION AND RELEVANCE

A single-nucleotide polymorphism (rs10911021) was identified that was significantly associated with CHD among persons with diabetes but not in those without diabetes and was functionally related to glutamic acid metabolism, suggesting a mechanistic link.

Causal and synthetic associations of variants in the SERPINA gene cluster with alpha1-antitrypsin serum levels

Several infrequent genetic polymorphisms in the SERPINA1 gene are known to substantially reduce concentration of alpha1-antitrypsin (AAT) in the blood. Since low AAT serum levels fail to protect pulmonary tissue from enzymatic degradation, these polymorphisms also increase the risk for early onset chronic obstructive pulmonary disease (COPD). The role of more common SERPINA1 single nucleotide polymorphisms (SNPs) in respiratory health remains poorly understood.

We present here an agnostic investigation of genetic determinants of circulating AAT levels in a general population sample by performing a genome-wide association study (GWAS) in 1392 individuals of the SAPALDIA cohort.

Five common SNPs, defined by showing minor allele frequencies (MAFs) >5%, reached genome-wide significance, all located in the SERPINA gene cluster at 14q32.13. The top-ranking genotyped SNP rs4905179 was associated with an estimated effect of β = −0.068 g/L per minor allele (P = 1.20*10⁻¹²). But denser SERPINA1 locus genotyping in 5569 participants with subsequent stepwise conditional analysis, as well as exon-sequencing in a subsample (N = 410), suggested that AAT serum level is causally determined at this locus by rare (MAF<1%) and low-frequent (MAF 1–5%) variants only, in particular by the well-documented protein inhibitor S and Z (PI S, PI Z) variants. Replication of the association of rs4905179 with AAT serum levels in the Copenhagen City Heart Study (N = 8273) was successful (P<0.0001), as was the replication of its synthetic nature (the effect disappeared after adjusting for PI S and Z, P = 0.57). Extending the analysis to lung function revealed a more complex situation. Only in individuals with severely compromised pulmonary health (N = 397), associations of common SNPs at this locus with lung function were driven by rarer PI S or Z variants. Overall, our meta-analysis of lung function in ever-smokers does not support a functional role of common SNPs in the SERPINA gene cluster in the general population.

Low levels of alpha1-antitrypsin (AAT) in the blood are a well-established risk factor for accelerated loss in lung function and chronic obstructive pulmonary disease. While a few infrequent genetic polymorphisms are known to influence the serum levels of this enzyme, the role of common genetic variants has not been examined so far. The present genome-wide scan for associated variants in approximately 1400 Swiss inhabitants revealed a chromosomal locus containing the functionally established variants of AAT deficiency and variants previously associated with lung function and emphysema. We used dense genotyping of this genetic region in more than 5500 individuals and subsequent conditional analyses to unravel which of these associated variants contribute independently to the phenotype's variability. All associations of common variants could be attributed to the rarer functionally established variants, a result which was then replicated in an independent population-based Danish cohort. Hence, this locus represents a textbook example of how a large part of a trait's heritability can be hidden in infrequent genetic polymorphisms. The attempt to transfer these results to lung function furthermore suggests that effects of common variants in this genetic region in ever-smokers may also be explained by rarer variants, but only in individuals with hampered pulmonary health.

Genome-wide association study of a heart failure related metabolomic profile among African Americans in the Atherosclerosis Risk in Communities (ARIC) study

Both the prevalence and incidence of heart failure (HF) are increasing, especially among African Americans, but no large-scale, genome-wide association study (GWAS) of HF-related metabolites has been reported. We sought to identify novel genetic variants that are associated with metabolites previously reported to relate to HF incidence. GWASs of three metabolites identified previously as risk factors for incident HF (pyroglutamine, dihydroxy docosatrienoic acid, and X-11787, being either hydroxy-leucine or hydroxy-isoleucine) were performed in 1,260 African Americans free of HF at the baseline examination of the Atherosclerosis Risk in Communities (ARIC) study. A significant association on chromosome 5q33 (rs10463316, MAF = 0.358, P-value = 1.92 × 10(-10) ) was identified for pyroglutamine. One region on chromosome 2p13 contained a nonsynonymous substitution in N-acetyltransferase 8 (NAT8) was associated with X-11787 (rs13538, MAF = 0.481, P-value = 1.71 × 10(-23) ). The smallest P-value for dihydroxy docosatrienoic acid was rs4006531 on chromosome 8q24 (MAF = 0.400, P-value = 6.98 × 10(-7) ). None of the above SNPs were individually associated with incident HF, but a genetic risk score (GRS) created by summing the most significant risk alleles from each metabolite detected 11% greater risk of HF per allele. In summary, we identified three loci associated with previously reported HF-related metabolites. Further use of metabolomics technology will facilitate replication of these findings in independent samples.

Association of heat shock proteins with all-cause mortality

Experimental mild heat shock is widely known as an intervention that results in extended longevity in various models along the evolutionary lineage. Heat shock proteins (HSPs) are highly upregulated immediately after a heat shock. The elevation in HSP levels was shown to inhibit stress-mediated cell death, and recent experiments indicate a highly versatile role for these proteins as inhibitors of programmed cell death. In this study, we examined common genetic variations in 31 genes encoding all members of the HSP70, small HSP, and heat shock factor (HSF) families for their association with all-cause mortality. Our discovery cohort was the Rotterdam study (RS1) containing 5,974 participants aged 55 years and older (3,174 deaths). We assessed 4,430 single nucleotide polymorphisms (SNPs) using the HumanHap550K Genotyping BeadChip from Illumina. After adjusting for multiple testing by permutation analysis, three SNPs showed evidence for association with all-cause mortality in RS1. These findings were followed in eight independent population-based cohorts, leading to a total of 25,007 participants (8,444 deaths). In the replication phase, only HSF2 (rs1416733) remained significantly associated with all-cause mortality. Rs1416733 is a known cis-eQTL for HSF2. Our findings suggest a role of HSF2 in all-cause mortality.

Dense mapping of IL2RA shows no association with Graves' disease in Chinese Han population

OBJECTIVE

Associations between IL2RA and various autoimmune diseases have been reported in Caucasians. We investigated whether genetic polymorphisms at the IL2RA locus were associated with Graves' disease (GD) in the Chinese Han population.

DESIGN

We performed a genome-wide association study (GWAS) in 1 536 GD patients and 1 516 controls. The 1000 Genomes Project data were adopted as references for imputation analysis. After forward and conditional logistic regressions, we found that rs11256313 was the major risk variant in the CD25/IL2RA region. Thus, we further genotyped rs11256313 in a replication cohort with 3 694 GD patients and 3 510 controls using ABI 7900HT TaqMan Real-Time PCR System.

RESULTS

Nine single nucleotide polymorphisms (SNPs) in the IL2RA block were nominally associated with GD in our GWAS (0·01 < P < 0·05). After imputation analysis, 13 imputed SNPs in the IL2RA block were weakly associated with GD (P ≤ 0·05). Logistic regression analysis suggested that the imputed rs11256313 could represent the IL2RA block (P = 0·003). However, we failed to replicate the association of rs11256313 in a larger cohort (P = 0·145). A subphenotype analysis of rs11256313 on thyroid hormone receptor antibody (TRAb) and gender showed that there was no association in any of the subphenotype groups (P > 0·05).

CONCLUSIONS

The results suggested that common genetic polymorphisms at IL2RA do not exert a significant genetic effect on the development of GD in the Chinese Han population. Previously reported associations between CD25/IL2RA and autoimmune diseases including GD in Caucasians again imply that heterogeneity exists in different ethnic populations.

A polymorphism affecting MYB binding within the promoter of the PDCD4 gene is associated with severe asthma in children

A previous genome-wide association study in asthma revealed putative associations that merit further investigation. In this study, the genome-wide significant associations of SNPs at the 5% false discovery rate were examined in independent groups of severe asthmatics. The panel consisted of 397 severe asthmatic adults, 116 severe asthmatic children, and a collection of 207 family-trios with an asthmatic proband. Three SNPs in the PDCD4 gene (rs6585018:G>A, rs1322997:C>A, and rs34104444:G>A) were significantly associated with severe childhood asthma (P values: 0.003, 0.002, 0.004) and total immunoglobulin E (IgE) levels (P values: 0.034, 0.041, 0.052). In an independent group of 234 asthmatic children and 652 controls, PDCD4 SNPs rs1407696:T>G and rs11195360:T>C were associated with total IgE levels (P values: 0.006, 0.014). In silico analysis of PDCD4 locus showed that rs6585018:G>A had the potential to affect MYB transcription factor binding, shown to act as a PDCD4-transcription inducer. Electromobility shift assays and reporter assays revealed that rs6585018:G>A alters MYB binding thereby influencing the expression of PDCD4. SNPs within MYB itself confer susceptibility to eosinophilia and asthma. Our association between a variant MYB binding site in PDCD4 and the severest form of childhood asthma therefore suggests that PDCD4 is a novel molecule of importance to asthmatic inflammatory responses.

Fast association tests for genes with FAST

Gene-based tests of association can increase the power of a genome-wide association study by aggregating multiple independent effects across a gene or locus into a single stronger signal. Recent gene-based tests have distinct approaches to selecting which variants to aggregate within a locus, modeling the effects of linkage disequilibrium, representing fractional allele counts from imputation, and managing permutation tests for p-values. Implementing these tests in a single, efficient framework has great practical value. Fast ASsociation Tests (Fast) addresses this need by implementing leading gene-based association tests together with conventional SNP-based univariate tests and providing a consolidated, easily interpreted report. Fast scales readily to genome-wide SNP data with millions of SNPs and tens of thousands of individuals, provides implementations that are orders of magnitude faster than original literature reports, and provides a unified framework for performing several gene based association tests concurrently and efficiently on the same data.

AVAILABILITY

https://bitbucket.org/baderlab/fast/downloads/FAST.tar.gz, with documentation at https://bitbucket.org/baderlab/fast/wiki/Home.

Genetic association analysis and meta-analysis of imputed SNPs in longitudinal studies

In this paper we propose a new method to analyze time-to-event data in longitudinal genetic studies. This method address the fundamental problem of incorporating uncertainty when analyzing survival data and imputed single-nucleotide polymorphisms (SNPs) from genome-wide association studies (GWAS). Our method incorporates uncertainty in the likelihood function, the opposite of existing methods that incorporate the uncertainty in the design matrix. Through simulation studies and real data analyses, we show that our proposed method is unbiased and provides powerful results. We also show how combining results from different GWAS (meta-analysis) may lead to wrong results when effects are not estimated using our approach. The model is implemented in an R package that is designed to analyze uncertainty not only arising from imputed SNPs, but also from copy number variants.

Novel genetic variants associated with lumbar disc degeneration in northern Europeans: a meta-analysis of 4600 subjects

OBJECTIVE

Lumbar disc degeneration (LDD) is an important cause of low back pain, which is a common and costly problem. LDD is characterised by disc space narrowing and osteophyte growth at the circumference of the disc. To date, the agnostic search of the genome by genome-wide association (GWA) to identify common variants associated with LDD has not been fruitful. This study is the first GWA meta-analysis of LDD.

METHODS

We have developed a continuous trait based on disc space narrowing and osteophytes growth which is measurable on all forms of imaging (plain radiograph, CT scan and MRI) and performed a meta-analysis of five cohorts of Northern European extraction each having GWA data imputed to HapMap V.2.

RESULTS

This study of 4600 individuals identified four single nucleotide polymorphisms with p<5×10(-8), the threshold set for genome-wide significance. We identified a variant in the PARK2 gene (p=2.8×10(-8)) associated with LDD. Differential methylation at one CpG island of the PARK2 promoter was observed in a small subset of subjects (β=8.74×10(-4), p=0.006).

CONCLUSIONS

LDD accounts for a considerable proportion of low back pain and the pathogenesis of LDD is poorly understood. This work provides evidence of association of the PARK2 gene and suggests that methylation of the PARK2 promoter may influence degeneration of the intervertebral disc. This gene has not previously been considered a candidate in LDD and further functional work is needed on this hitherto unsuspected pathway.

Genome-wide association and longitudinal analyses reveal genetic loci linking pubertal height growth, pubertal timing and childhood adiposity

The pubertal height growth spurt is a distinctive feature of childhood growth reflecting both the central onset of puberty and local growth factors. Although little is known about the underlying genetics, growth variability during puberty correlates with adult risks for hormone-dependent cancer and adverse cardiometabolic health. The only gene so far associated with pubertal height growth, LIN28B, pleiotropically influences childhood growth, puberty and cancer progression, pointing to shared underlying mechanisms. To discover genetic loci influencing pubertal height and growth and to place them in context of overall growth and maturation, we performed genome-wide association meta-analyses in 18 737 European samples utilizing longitudinally collected height measurements. We found significant associations (P < 1.67 × 10(-8)) at 10 loci, including LIN28B. Five loci associated with pubertal timing, all impacting multiple aspects of growth. In particular, a novel variant correlated with expression of MAPK3, and associated both with increased prepubertal growth and earlier menarche. Another variant near ADCY3-POMC associated with increased body mass index, reduced pubertal growth and earlier puberty. Whereas epidemiological correlations suggest that early puberty marks a pathway from rapid prepubertal growth to reduced final height and adult obesity, our study shows that individual loci associating with pubertal growth have variable longitudinal growth patterns that may differ from epidemiological observations. Overall, this study uncovers part of the complex genetic architecture linking pubertal height growth, the timing of puberty and childhood obesity and provides new information to pinpoint processes linking these traits.

A genome-wide association study of sleep habits and insomnia

Several aspects of sleep behavior such as timing, duration and quality have been demonstrated to be heritable. To identify common variants that influence sleep traits in the population, we conducted a genome-wide association study of six sleep phenotypes assessed by questionnaire in a sample of 2,323 individuals from the Australian Twin Registry. Genotyping was performed on the Illumina 317, 370, and 610K arrays and the SNPs in common between platforms were used to impute non-genotyped SNPs. We tested for association with more than 2,000,000 common polymorphisms across the genome. While no SNPs reached the genome-wide significance threshold, we identified a number of associations in plausible candidate genes. Most notably, a group of SNPs in the third intron of the CACNA1C gene ranked as most significant in the analysis of sleep latency (P = 1.3 × 10⁻⁶). We attempted to replicate this association in an independent sample from the Chronogen Consortium (n = 2,034), but found no evidence of association (P = 0.73). We have identified several other suggestive associations that await replication in an independent sample. We did not replicate the results from previous genome-wide analyses of self-reported sleep phenotypes after correction for multiple testing.

Genetic loci for retinal arteriolar microcirculation

Narrow arterioles in the retina have been shown to predict hypertension as well as other vascular diseases, likely through an increase in the peripheral resistance of the microcirculatory flow. In this study, we performed a genome-wide association study in 18,722 unrelated individuals of European ancestry from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium and the Blue Mountain Eye Study, to identify genetic determinants associated with variations in retinal arteriolar caliber. Retinal vascular calibers were measured on digitized retinal photographs using a standardized protocol. One variant (rs2194025 on chromosome 5q14 near the myocyte enhancer factor 2C MEF2C gene) was associated with retinal arteriolar caliber in the meta-analysis of the discovery cohorts at genome-wide significance of P-value <5×10⁻⁸. This variant was replicated in an additional 3,939 individuals of European ancestry from the Australian Twins Study and Multi-Ethnic Study of Atherosclerosis (rs2194025, P-value = 2.11×10⁻¹² in combined meta-analysis of discovery and replication cohorts). In independent studies of modest sample sizes, no significant association was found between this variant and clinical outcomes including coronary artery disease, stroke, myocardial infarction or hypertension. In conclusion, we found one novel loci which underlie genetic variation in microvasculature which may be relevant to vascular disease. The relevance of these findings to clinical outcomes remains to be determined.

Genomic association analysis of common variants influencing antihypertensive response to hydrochlorothiazide

To identify novel genes influencing blood pressure response to thiazide diuretic therapy for hypertension, we conducted genome-wide association meta-analyses of ≈1.1 million single-nucleotide polymorphisms in a combined sample of 424 European Americans with primary hypertension treated with hydrochlorothiazide from the Pharmacogenomic Evaluation of Antihypertensive Responses study (n=228) and the Genetic Epidemiology of Responses to Antihypertensive study (n=196). Polymorphisms associated with blood pressure response at P<10(-5) were tested for replication of the associations in independent samples of hydrochlorothiazide-treated European hypertensives. The rs16960228 polymorphism in protein kinase C, α replicated for same-direction association with diastolic blood pressure response in the Nordic Diltiazem study (n=420) and the Genetics of Drug Responsiveness in Essential Hypertension study (n=206), and the combined 4-study meta-analysis P value achieved genome-wide significance (P=3.3 × 10(-8)). Systolic or diastolic blood pressure responses were consistently greater in carriers of the rs16960228 A allele than in GG homozygotes (>4/4 mm Hg) across study samples. The rs2273359 polymorphism in the GNAS-EDN3 region also replicated for same-direction association with systolic blood pressure response in the Nordic Diltiazem study, and the combined 3-study meta-analysis P value approached genome-wide significance (P=5.5 × 10(-8)). The findings document clinically important effects of genetic variation at novel loci on blood pressure response to a thiazide diuretic, which may be a basis for individualization of antihypertensive drug therapy and identification of new drug targets.

Sex-stratified genome-wide association studies including 270,000 individuals show sexual dimorphism in genetic loci for anthropometric traits

Given the anthropometric differences between men and women and previous evidence of sex-difference in genetic effects, we conducted a genome-wide search for sexually dimorphic associations with height, weight, body mass index, waist circumference, hip circumference, and waist-to-hip-ratio (133,723 individuals) and took forward 348 SNPs into follow-up (additional 137,052 individuals) in a total of 94 studies. Seven loci displayed significant sex-difference (FDR<5%), including four previously established (near GRB14/COBLL1, LYPLAL1/SLC30A10, VEGFA, ADAMTS9) and three novel anthropometric trait loci (near MAP3K1, HSD17B4, PPARG), all of which were genome-wide significant in women (P<5×10(-8)), but not in men. Sex-differences were apparent only for waist phenotypes, not for height, weight, BMI, or hip circumference. Moreover, we found no evidence for genetic effects with opposite directions in men versus women. The PPARG locus is of specific interest due to its role in diabetes genetics and therapy. Our results demonstrate the value of sex-specific GWAS to unravel the sexually dimorphic genetic underpinning of complex traits.

GWAS on your notebook: fast semi-parallel linear and logistic regression for genome-wide association studies

Background

Genome-wide association studies have become very popular in identifying genetic contributions to phenotypes. Millions of SNPs are being tested for their association with diseases and traits using linear or logistic regression models. This conceptually simple strategy encounters the following computational issues: a large number of tests and very large genotype files (many Gigabytes) which cannot be directly loaded into the software memory. One of the solutions applied on a grand scale is cluster computing involving large-scale resources. We show how to speed up the computations using matrix operations in pure R code.

Results

We improve speed: computation time from 6 hours is reduced to 10-15 minutes. Our approach can handle essentially an unlimited amount of covariates efficiently, using projections. Data files in GWAS are vast and reading them into computer memory becomes an important issue. However, much improvement can be made if the data is structured beforehand in a way allowing for easy access to blocks of SNPs. We propose several solutions based on the R packages ff and ncdf.

We adapted the semi-parallel computations for logistic regression. We show that in a typical GWAS setting, where SNP effects are very small, we do not lose any precision and our computations are few hundreds times faster than standard procedures.

Conclusions

We provide very fast algorithms for GWAS written in pure R code. We also show how to rearrange SNP data for fast access.

17q25 Locus is associated with white matter hyperintensity volume in ischemic stroke, but not with lacunar stroke status

BACKGROUND AND PURPOSE

Recently, a novel locus at 17q25 was associated with white matter hyperintensities (WMH) on MRI in stroke-free individuals. We aimed to replicate the association with WMH volume (WMHV) in patients with ischemic stroke. If the association acts by promoting a small vessel arteriopathy, it might be expected to also associate with lacunar stroke.

METHODS

We quantified WMH on MRI in the stroke-free hemisphere of 2588 ischemic stroke cases. Association between WMHV and 6 single-nucleotide polymorphisms at chromosome 17q25 was assessed by linear regression. These single-nucleotide polymorphisms were also investigated for association with lacunar stroke in 1854 cases and 51 939 stroke-free controls from METASTROKE. Meta-analyses with previous reports and a genetic risk score approach were applied to identify other novel WMHV risk variants and uncover shared genetic contributions to WMHV in community participants without stroke and ischemic stroke.

RESULTS

Single-nucleotide polymorphisms at 17q25 were associated with WMHV in ischemic stroke, the most significant being rs9894383 (P=0.0006). In contrast, there was no association between any single-nucleotide polymorphism and lacunar stroke. A genetic risk score analysis revealed further genetic components to WMHV shared between community participants without stroke and ischemic stroke.

CONCLUSIONS

This study provides support for an association between the 17q25 locus and WMH. In contrast, it is not associated with lacunar stroke, suggesting that the association does not act by promoting small-vessel arteriopathy or the same arteriopathy responsible for lacunar infarction.

Genome-wide analysis of BMI in adolescents and young adults reveals additional insight into the effects of genetic loci over the life course

Genetic loci for body mass index (BMI) in adolescence and young adulthood, a period of high risk for weight gain, are understudied, yet may yield important insight into the etiology of obesity and early intervention. To identify novel genetic loci and examine the influence of known loci on BMI during this critical time period in late adolescence and early adulthood, we performed a two-stage meta-analysis using 14 genome-wide association studies in populations of European ancestry with data on BMI between ages 16 and 25 in up to 29 880 individuals. We identified seven independent loci (P < 5.0 × 10⁻⁸) near FTO (P = 3.72 × 10⁻²³), TMEM18 (P = 3.24 × 10⁻¹⁷), MC4R (P = 4.41 × 10⁻¹⁷), TNNI3K (P = 4.32 × 10⁻¹¹), SEC16B (P = 6.24 × 10⁻⁹), GNPDA2 (P = 1.11 × 10⁻⁸) and POMC (P = 4.94 × 10⁻⁸) as well as a potential secondary signal at the POMC locus (rs2118404, P = 2.4 × 10⁻⁵ after conditioning on the established single-nucleotide polymorphism at this locus) in adolescents and young adults. To evaluate the impact of the established genetic loci on BMI at these young ages, we examined differences between the effect sizes of 32 published BMI loci in European adult populations (aged 18-90) and those observed in our adolescent and young adult meta-analysis. Four loci (near PRKD1, TNNI3K, SEC16B and CADM2) had larger effects and one locus (near SH2B1) had a smaller effect on BMI during adolescence and young adulthood compared with older adults (P < 0.05). These results suggest that genetic loci for BMI can vary in their effects across the life course, underlying the importance of evaluating BMI at different ages.

One thousand genomes imputation in the National Cancer Institute Breast and Prostate Cancer Cohort Consortium aggressive prostate cancer genome-wide association study

BACKGROUND

Genotype imputation substantially increases available markers for analysis in genome-wide association studies (GWAS) by leveraging linkage disequilibrium from a reference panel. We sought to (i) investigate the performance of imputation from the August 2010 release of the 1000 Genomes Project (1000GP) in an existing GWAS of prostate cancer, (ii) look for novel associations with prostate cancer risk, (iii) fine-map known prostate cancer susceptibility regions using an approximate Bayesian framework and stepwise regression, and (iv) compare power and efficiency of imputation and de novo sequencing.

METHODS

We used 2,782 aggressive prostate cancer cases and 4,458 controls from the NCI Breast and Prostate Cancer Cohort Consortium aggressive prostate cancer GWAS to infer 5.8 million well-imputed autosomal single nucleotide polymorphisms (SNPs).

RESULTS

Imputation quality, as measured by correlation between imputed and true allele counts, was higher among common variants than rare variants. We found no novel prostate cancer associations among a subset of 1.2 million well-imputed low-frequency variants. At a genome-wide sequencing cost of $2,500, imputation from SNP arrays is a more powerful strategy than sequencing for detecting disease associations of SNPs with minor allele frequencies (MAF) above 1%.

CONCLUSIONS

1000GP imputation provided dense coverage of previously identified prostate cancer susceptibility regions, highlighting its potential as an inexpensive first-pass approach to fine mapping in regions such as 5p15 and 8q24. Our study shows 1000GP imputation can accurately identify low-frequency variants and stresses the importance of large sample size when studying these variants.

Gene variants in angiogenesis and lymphangiogenesis and cutaneous melanoma progression

BACKGROUND

Angiogenesis and lymphangiogenesis are important in the progression of melanoma. We investigated associations between genetic variants in these pathways with sentinel lymph node (SLN) metastasis and mortality in 2 independent series of patients with melanoma.

METHODS

Participants at Moffitt Cancer Center were 552 patients, all Caucasian, with primary cutaneous melanoma referred for SLN biopsy. A total of 177 patients had SLN metastasis, among whom 60 died from melanoma. Associations between 238 single-nucleotide polymorphisms (SNP) in 26 genes and SLN metastasis were estimated as ORs and 95% confidence intervals (CI) using logistic regression. Competing risk regression was used to estimate HRs and 95% CI for each SNP and melanoma-specific mortality. We attempted to replicate significant findings using data from a genome-wide association study comprising 1,115 patients with melanoma who were referred for SLN biopsy from MD Anderson Cancer Center (MDACC), among whom 189 patients had SLN metastasis and 92 patients died from melanoma.

RESULTS

In the Moffitt dataset, we observed significant associations in 18 SNPs with SLN metastasis and 17 SNPs with mortality. Multiple SNPs in COL18A1, EGF receptor (EGFR), FLT1, interleukin (IL)-10, platelet-derived growth factor D (PDGFD), PIK3CA, and toll-like receptor (TLR)-3 were associated with the risk of SLN metastasis and/or patient mortality. The MDACC data set replicated an association between mortality and rs2220377 in PDGFD. Furthermore, in a meta-analysis, 3 additional SNPs were significantly associated with SLN metastasis (EGFR rs723526 and TLR3 rs3775292) and melanoma-specific death (TLR3 rs7668666).

CONCLUSIONS

These findings suggest that genetic variation in angiogenesis and lymphangiogenesis contributes to regional nodal metastasis and progression of melanoma.

IMPACT

Additional research attempting to replicate these results is warranted.

Genome-wide association study of age at menarche in African-American women

African-American (AA) women have earlier menarche on average than women of European ancestry (EA), and earlier menarche is a risk factor for obesity and type 2 diabetes among other chronic diseases. Identification of common genetic variants associated with age at menarche has a potential value in pointing to the genetic pathways underlying chronic disease risk, yet comprehensive genome-wide studies of age at menarche are lacking for AA women. In this study, we tested the genome-wide association of self-reported age at menarche with common single-nucleotide polymorphisms (SNPs) in a total of 18 089 AA women in 15 studies using an additive genetic linear regression model, adjusting for year of birth and population stratification, followed by inverse-variance weighted meta-analysis (Stage 1). Top meta-analysis results were then tested in an independent sample of 2850 women (Stage 2). First, while no SNP passed the pre-specified P < 5 × 10(-8) threshold for significance in Stage 1, suggestive associations were found for variants near FLRT2 and PIK3R1, and conditional analysis identified two independent SNPs (rs339978 and rs980000) in or near RORA, strengthening the support for this suggestive locus identified in EA women. Secondly, an investigation of SNPs in 42 previously identified menarche loci in EA women demonstrated that 25 (60%) of them contained variants significantly associated with menarche in AA women. The findings provide the first evidence of cross-ethnic generalization of menarche loci identified to date, and suggest a number of novel biological links to menarche timing in AA women.

Identification of heart rate-associated loci and their effects on cardiac conduction and rhythm disorders

Elevated resting heart rate is associated with greater risk of cardiovascular disease and mortality. In a 2-stage meta-analysis of genome-wide association studies in up to 181,171 individuals, we identified 14 new loci associated with heart rate and confirmed associations with all 7 previously established loci. Experimental downregulation of gene expression in Drosophila melanogaster and Danio rerio identified 20 genes at 11 loci that are relevant for heart rate regulation and highlight a role for genes involved in signal transmission, embryonic cardiac development and the pathophysiology of dilated cardiomyopathy, congenital heart failure and/or sudden cardiac death. In addition, genetic susceptibility to increased heart rate is associated with altered cardiac conduction and reduced risk of sick sinus syndrome, and both heart rate-increasing and heart rate-decreasing variants associate with risk of atrial fibrillation. Our findings provide fresh insights into the mechanisms regulating heart rate and identify new therapeutic targets.

Identification of 23 new prostate cancer susceptibility loci using the iCOGS custom genotyping array

Prostate cancer is the most frequently diagnosed cancer in males in developed countries. To identify common prostate cancer susceptibility alleles, we genotyped 211,155 SNPs on a custom Illumina array (iCOGS) in blood DNA from 25,074 prostate cancer cases and 24,272 controls from the international PRACTICAL Consortium. Twenty-three new prostate cancer susceptibility loci were identified at genome-wide significance (P < 5 × 10(-8)). More than 70 prostate cancer susceptibility loci, explaining ∼30% of the familial risk for this disease, have now been identified. On the basis of combined risks conferred by the new and previously known risk loci, the top 1% of the risk distribution has a 4.7-fold higher risk than the average of the population being profiled. These results will facilitate population risk stratification for clinical studies.

Large-scale genotyping identifies 41 new loci associated with breast cancer risk

Breast cancer is the most common cancer among women. Common variants at 27 loci have been identified as associated with susceptibility to breast cancer, and these account for ∼9% of the familial risk of the disease. We report here a meta-analysis of 9 genome-wide association studies, including 10,052 breast cancer cases and 12,575 controls of European ancestry, from which we selected 29,807 SNPs for further genotyping. These SNPs were genotyped in 45,290 cases and 41,880 controls of European ancestry from 41 studies in the Breast Cancer Association Consortium (BCAC). The SNPs were genotyped as part of a collaborative genotyping experiment involving four consortia (Collaborative Oncological Gene-environment Study, COGS) and used a custom Illumina iSelect genotyping array, iCOGS, comprising more than 200,000 SNPs. We identified SNPs at 41 new breast cancer susceptibility loci at genome-wide significance (P < 5 × 10(-8)). Further analyses suggest that more than 1,000 additional loci are involved in breast cancer susceptibility.

Genome-wide association studies identify four ER negative-specific breast cancer risk loci

Estrogen receptor (ER)-negative tumors represent 20-30% of all breast cancers, with a higher proportion occurring in younger women and women of African ancestry. The etiology and clinical behavior of ER-negative tumors are different from those of tumors expressing ER (ER positive), including differences in genetic predisposition. To identify susceptibility loci specific to ER-negative disease, we combined in a meta-analysis 3 genome-wide association studies of 4,193 ER-negative breast cancer cases and 35,194 controls with a series of 40 follow-up studies (6,514 cases and 41,455 controls), genotyped using a custom Illumina array, iCOGS, developed by the Collaborative Oncological Gene-environment Study (COGS). SNPs at four loci, 1q32.1 (MDM4, P = 2.1 × 10(-12) and LGR6, P = 1.4 × 10(-8)), 2p24.1 (P = 4.6 × 10(-8)) and 16q12.2 (FTO, P = 4.0 × 10(-8)), were associated with ER-negative but not ER-positive breast cancer (P > 0.05). These findings provide further evidence for distinct etiological pathways associated with invasive ER-positive and ER-negative breast cancers.

Using multiple measures for quantitative trait association analyses: application to estimated glomerular filtration rate

Studies of multiple measures of a quantitative trait can have greater precision and thus statistical power compared with single-measure studies, but this has rarely been studied in the relation to quantitative trait measurement error models in genetic association studies. Using estimated glomerular filtration rate (eGFR), a quantitative measure of kidney function, as an example we constructed measurement error models of a quantitative trait with systematic and random error components. We then examined the effects on precision of the parameter estimate between genetic loci and eGFR, resulting from varying the correlation and contribution of the error components. We also compared the empirical results from three genome-wide association studies (GWAS) of kidney function in 9049 European Americans: a single measure model, a three-measure model of the same biomarker of kidney function and a six-measure model of different biomarkers of kidney function. Simulations showed that given the same amount of overall errors, inclusion of measures with less correlated systematic errors led to greater gain in precision. The empirical GWAS results confirmed that both the three- and six-measure models detected more eGFR-associated genomic loci with stronger statistical association than the single-measure model despite some heterogeneity among the measures. Multiple measures of a quantitative trait can increase the statistical power of a study without additional participant recruitment. However, careful attention must be paid to the correlation of systematic errors and inconsistent associations when different biomarkers or methods are used to measure the quantitative trait.

Pleiotropy and pathway analyses of genetic variants associated with both type 2 diabetes and prostate cancer

AIMS

Epidemiological evidence shows that diabetes is associated with a reduced risk of prostate cancer. The objective of this study was to identify genes that may contribute to both type 2 diabetes and prostate cancer outcomes and the biological pathways these diseases may share.

METHODS

The Atherosclerosis Risk in Communities (ARIC) Study is a population-based prospective cohort study in four U.S. communities that included a baseline examination in 1987-89 and three follow-up exams at three year intervals. Participants were 45-64 years old at baseline. We conducted a genomewide association (GWA) study of incident type 2 diabetes in males, summarized variation across genetic loci into a polygenic risk score, and determined if that diabetes risk score was also associated with incident prostate cancer in the same study population. Secondarily we conducted a separate GWA study of prostate cancer, performed a pathway analysis of both type 2 diabetes and prostate cancer, and qualitatively determined if any of the biochemical pathways identified were shared between the two outcomes.

RESULTS

We found that the polygenic risk score for type 2 diabetes was not statistically significantly associated with prostate cancer. The pathway analysis also found no overlap between pathways associated with type 2 diabetes and prostate cancer. However, it did find that the growth hormone signaling pathway was statistically significantly associated with type 2 diabetes (p=0.0001).

CONCLUSION

The inability of this study to find an association between type 2 diabetes polygenic risk scores with prostate cancer or biological pathways in common suggests that shared genetic variants may not contribute significantly to explaining shared etiology.

Exploring the genetic basis of chronic periodontitis: a genome-wide association study

Chronic periodontitis (CP) is a common oral disease that confers substantial systemic inflammatory and microbial burden and is a major cause of tooth loss. Here, we present the results of a genome-wide association study of CP that was carried out in a cohort of 4504 European Americans (EA) participating in the Atherosclerosis Risk in Communities (ARIC) Study (mean age—62 years, moderate CP—43% and severe CP—17%). We detected no genome-wide significant association signals for CP; however, we found suggestive evidence of association (P < 5 × 10⁻⁶) for six loci, including NIN, NPY, WNT5A for severe CP and NCR2, EMR1, 10p15 for moderate CP. Three of these loci had concordant effect size and direction in an independent sample of 656 adult EA participants of the Health, Aging, and Body Composition (Health ABC) Study. Meta-analysis pooled estimates were severe CP (n = 958 versus health: n = 1909)—NPY, rs2521634 [G]: odds ratio [OR = 1.49 (95% confidence interval (CI = 1.28–1.73, P = 3.5 × 10⁻⁷))]; moderate CP (n = 2293)—NCR2, rs7762544 [G]: OR = 1.40 (95% CI = 1.24–1.59, P = 7.5 × 10⁻⁸), EMR1, rs3826782 [A]: OR = 2.01 (95% CI = 1.52–2.65, P = 8.2 × 10⁻⁷). Canonical pathway analysis indicated significant enrichment of nervous system signaling, cellular immune response and cytokine signaling pathways. A significant interaction of NUAK1 (rs11112872, interaction P = 2.9 × 10⁻⁹) with smoking in ARIC was not replicated in Health ABC, although estimates of heritable variance in severe CP explained by all single nucleotide polymorphisms increased from 18 to 52% with the inclusion of a genome-wide interaction term with smoking. These genome-wide association results provide information on multiple candidate regions and pathways for interrogation in future genetic studies of CP.

Seven new loci associated with age-related macular degeneration

Age-related macular degeneration (AMD) is a common cause of blindness in older individuals. To accelerate the understanding of AMD biology and help design new therapies, we executed a collaborative genome-wide association study, including >17,100 advanced AMD cases and >60,000 controls of European and Asian ancestry. We identified 19 loci associated at P < 5 × 10(-8). These loci show enrichment for genes involved in the regulation of complement activity, lipid metabolism, extracellular matrix remodeling and angiogenesis. Our results include seven loci with associations reaching P < 5 × 10(-8) for the first time, near the genes COL8A1-FILIP1L, IER3-DDR1, SLC16A8, TGFBR1, RAD51B, ADAMTS9 and B3GALTL. A genetic risk score combining SNP genotypes from all loci showed similar ability to distinguish cases and controls in all samples examined. Our findings provide new directions for biological, genetic and therapeutic studies of AMD.

Genome-wide meta-analyses of multiancestry cohorts identify multiple new susceptibility loci for refractive error and myopia

Refractive error is the most common eye disorder worldwide and is a prominent cause of blindness. Myopia affects over 30% of Western populations and up to 80% of Asians. The CREAM consortium conducted genome-wide meta-analyses, including 37,382 individuals from 27 studies of European ancestry and 8,376 from 5 Asian cohorts. We identified 16 new loci for refractive error in individuals of European ancestry, of which 8 were shared with Asians. Combined analysis identified 8 additional associated loci. The new loci include candidate genes with functions in neurotransmission (GRIA4), ion transport (KCNQ5), retinoic acid metabolism (RDH5), extracellular matrix remodeling (LAMA2 and BMP2) and eye development (SIX6 and PRSS56). We also confirmed previously reported associations with GJD2 and RASGRF1. Risk score analysis using associated SNPs showed a tenfold increased risk of myopia for individuals carrying the highest genetic load. Our results, based on a large meta-analysis across independent multiancestry studies, considerably advance understanding of the mechanisms involved in refractive error and myopia.

Genetic determinants of trabecular and cortical volumetric bone mineral densities and bone microstructure

Most previous genetic epidemiology studies within the field of osteoporosis have focused on the genetics of the complex trait areal bone mineral density (aBMD), not being able to differentiate genetic determinants of cortical volumetric BMD (vBMD), trabecular vBMD, and bone microstructural traits. The objective of this study was to separately identify genetic determinants of these bone traits as analysed by peripheral quantitative computed tomography (pQCT). Separate GWA meta-analyses for cortical and trabecular vBMDs were performed. The cortical vBMD GWA meta-analysis (n = 5,878) followed by replication (n = 1,052) identified genetic variants in four separate loci reaching genome-wide significance (RANKL, rs1021188, p = 3.6×10⁻¹⁴; LOC285735, rs271170, p = 2.7×10⁻¹²; OPG, rs7839059, p = 1.2×10⁻¹⁰; and ESR1/C6orf97, rs6909279, p = 1.1×10⁻⁹). The trabecular vBMD GWA meta-analysis (n = 2,500) followed by replication (n = 1,022) identified one locus reaching genome-wide significance (FMN2/GREM2, rs9287237, p = 1.9×10⁻⁹). High-resolution pQCT analyses, giving information about bone microstructure, were available in a subset of the GOOD cohort (n = 729). rs1021188 was significantly associated with cortical porosity while rs9287237 was significantly associated with trabecular bone fraction. The genetic variant in the FMN2/GREM2 locus was associated with fracture risk in the MrOS Sweden cohort (HR per extra T allele 0.75, 95% confidence interval 0.60–0.93) and GREM2 expression in human osteoblasts. In conclusion, five genetic loci associated with trabecular or cortical vBMD were identified. Two of these (FMN2/GREM2 and LOC285735) are novel bone-related loci, while the other three have previously been reported to be associated with aBMD. The genetic variants associated with cortical and trabecular bone parameters differed, underscoring the complexity of the genetics of bone parameters. We propose that a genetic variant in the RANKL locus influences cortical vBMD, at least partly, via effects on cortical porosity, and that a genetic variant in the FMN2/GREM2 locus influences GREM2 expression in osteoblasts and thereby trabecular number and thickness as well as fracture risk.

Osteoporosis is a common highly heritable skeletal disease characterized by reduced bone mineral density (BMD) and deteriorated bone microstructure, resulting in an increased risk of fracture. Most previous genetic epidemiology studies have focused on the genetics of the complex trait BMD, not being able to separate genetic determinants of the trabecular and cortical bone compartments and bone microstructure. The trabecular and cortical BMDs can be analysed separately by computed tomography. Therefore, we performed separate genome-wide association studies for trabecular and cortical BMDs, demonstrating that the genetic determinants of cortical and trabecular BMDs differ. Genetic variants in the RANKL, LOC285735, OPG, and ESR1 loci were associated with cortical BMD, while a genetic variant in the FMN2/GREM2 locus was associated with trabecular BMD. Two of these are novel bone-related loci. Follow-up analyses of bone microstructure demonstrated that a genetic variant in the RANKL locus is associated with cortical porosity and that the FMN2/GREM2 locus is associated with trabecular number and thickness. We propose that a genetic variant in the RANKL locus influences cortical BMD via effects on cortical porosity, and that a genetic variant in the FMN2/GREM2 locus influences trabecular BMD and fracture risk via effects on both trabecular number and thickness.

Genome-wide association study identified the human leukocyte antigen region as a novel locus for plasma beta-2 microglobulin

Beta-2 microglobulin (B2M) is a component of the major histocompatibility complex (MHC) class I molecule and has been studied as a biomarker of kidney function, cardiovascular diseases and mortality. Little is known about the genes influencing its levels directly or through glomerular filtration rate (GFR). We conducted a genome-wide association study of plasma B2M levels in 6738 European Americans from the Atherosclerosis Risk in Communities study to identify novel loci for B2M and assessed its association with known estimated GFR (eGFR) loci. We identified 2 genome-wide significant loci. One was in the human leukocyte antigen (HLA) region on chromosome 6 (lowest p value = 1.8 × 10(-23) for rs9264638). At this locus, 6 index SNPs accounted for 3.2 % of log(B2M) variance, and their association with B2M could largely be explained by imputed classical alleles of the MHC class I genes: HLA-A, HLA-B, or HLA-C. The index SNPs at this locus were not associated with eGFR based on serum creatinine (eGFRcr). The other locus of B2M was on chromosome 12 (rs3184504 at SH2B3, beta = 0.02, p value = 3.1 × 10(-8)), which was previously implicated as an eGFR locus. In conclusion, although B2M is known to be a component of MHC class I molecule, the association between HLA class I alleles and plasma B2M levels in a community-based population is novel. The identification of the two novel loci for B2M extends our understanding of its metabolism and informs its use as a kidney filtration biomarker.

Genome-wide association study identified the human leukocyte antigen region as a novel locus for plasma beta-2 microglobulin

Beta-2 microglobulin (B2M) is a component of the major histocompatibility complex (MHC) class I molecule and has been studied as a biomarker of kidney function, cardiovascular diseases and mortality. Little is known about the genes influencing its levels directly or through glomerular filtration rate (GFR). We conducted a genome-wide association study of plasma B2M levels in 6738 European Americans from the Atherosclerosis Risk in Communities study to identify novel loci for B2M and assessed its association with known estimated GFR (eGFR) loci. We identified 2 genome-wide significant loci. One was in the human leukocyte antigen (HLA) region on chromosome 6 (lowest p value = 1.8 × 10(-23) for rs9264638). At this locus, 6 index SNPs accounted for 3.2 % of log(B2M) variance, and their association with B2M could largely be explained by imputed classical alleles of the MHC class I genes: HLA-A, HLA-B, or HLA-C. The index SNPs at this locus were not associated with eGFR based on serum creatinine (eGFRcr). The other locus of B2M was on chromosome 12 (rs3184504 at SH2B3, beta = 0.02, p value = 3.1 × 10(-8)), which was previously implicated as an eGFR locus. In conclusion, although B2M is known to be a component of MHC class I molecule, the association between HLA class I alleles and plasma B2M levels in a community-based population is novel. The identification of the two novel loci for B2M extends our understanding of its metabolism and informs its use as a kidney filtration biomarker.

Genome-wide association analyses suggest NELL1 influences adverse metabolic response to HCTZ in African Americans

Hydrochlorothiazide (HCTZ) is one of the most widely prescribed antihypertensive medications. Although it is well known that HCTZ is associated with hyperglycemia and hypertriglyceridemia, the mechanisms underlying these adverse effects are not well understood. We performed a genome-wide association study and meta-analysis of the change in fasting plasma glucose and triglycerides in response to HCTZ from two different clinical trials: the Pharmacogenomic Evaluation of Antihypertensive Responses and the Genetic Epidemiology of Responses to Antihypertensive studies. Two single-nucleotide polymorphisms (rs12279250 and rs4319515 (r²=0.73)), located at 11p15.1 in the NELL1 gene, achieved genome-wide significance for association with change in fasting plasma triglycerides in African Americans, whereby each variant allele was associated with a 28 mg dl⁻¹ increase in the change in triglycerides. NELL1 encodes a cytoplasmic protein that contains epidermal growth factor-like repeats and has been shown to represses adipogenic differentiation. These findings may represent a novel mechanism underlying HCTZ-induced adverse metabolic effects.

Genome-wide association study of retinopathy in individuals without diabetes

Background

Mild retinopathy (microaneurysms or dot-blot hemorrhages) is observed in persons without diabetes or hypertension and may reflect microvascular disease in other organs. We conducted a genome-wide association study (GWAS) of mild retinopathy in persons without diabetes.

Methods

A working group agreed on phenotype harmonization, covariate selection and analytic plans for within-cohort GWAS. An inverse-variance weighted fixed effects meta-analysis was performed with GWAS results from six cohorts of 19,411 Caucasians. The primary analysis included individuals without diabetes and secondary analyses were stratified by hypertension status. We also singled out the results from single nucleotide polymorphisms (SNPs) previously shown to be associated with diabetes and hypertension, the two most common causes of retinopathy.

Results

No SNPs reached genome-wide significance in the primary analysis or the secondary analysis of participants with hypertension. SNP, rs12155400, in the histone deacetylase 9 gene (HDAC9) on chromosome 7, was associated with retinopathy in analysis of participants without hypertension, −1.3±0.23 (beta ± standard error), p = 6.6×10⁻⁹. Evidence suggests this was a false positive finding. The minor allele frequency was low (∼2%), the quality of the imputation was moderate (r² ∼0.7), and no other common variants in the HDAC9 gene were associated with the outcome. SNPs found to be associated with diabetes and hypertension in other GWAS were not associated with retinopathy in persons without diabetes or in subgroups with or without hypertension.

Conclusions

This GWAS of retinopathy in individuals without diabetes showed little evidence of genetic associations. Further studies are needed to identify genes associated with these signs in order to help unravel novel pathways and determinants of microvascular diseases.

Identification of a candidate gene for astigmatism

PURPOSE

Astigmatism is a common refractive error that reduces vision, where the curvature and refractive power of the cornea in one meridian are less than those of the perpendicular axis. It is a complex trait likely to be influenced by both genetic and environmental factors. Twin studies of astigmatism have found approximately 60% of phenotypic variance is explained by genetic factors. This study aimed to identify susceptibility loci for astigmatism.

METHODS

We performed a meta-analysis of seven genome-wide association studies that included 22,100 individuals of European descent, where astigmatism was defined as the number of diopters of cylinder prescription, using fixed effect inverse variance-weighted methods.

RESULTS

A susceptibility locus was identified with lead single nucleotide polymorphism rs3771395 on chromosome 2p13.3 (meta-analysis, P = 1.97 × 10(-7)) in the VAX2 gene. VAX2 plays an important role in the development of the dorsoventral axis of the eye. Animal studies have shown a gradient in astigmatism along the vertical plane, with corresponding changes in refraction, particularly in the ventral field.

CONCLUSIONS

This finding advances the understanding of refractive error, and provides new potential pathways to be evaluated with regard to the development of astigmatism.

Association of genome-wide variation with highly sensitive cardiac troponin-T levels in European Americans and Blacks: a meta-analysis from atherosclerosis risk in communities and cardiovascular health studies

BACKGROUND

High levels of cardiac troponin T, measured by a highly sensitive assay (hs-cTnT), are strongly associated with incident coronary heart disease and heart failure. To date, no large-scale genome-wide association study of hs-cTnT has been reported. We sought to identify novel genetic variants that are associated with hs-cTnT levels.

METHODS AND RESULTS

We performed a genome-wide association in 9491 European Americans and 2053 blacks free of coronary heart disease and heart failure from 2 prospective cohorts: the Atherosclerosis Risk in Communities Study and the Cardiovascular Health Study. Genome-wide association studies were conducted in each study and race stratum. Fixed-effect meta-analyses combined the results of linear regression from 2 cohorts within each race stratum and then across race strata to produce overall estimates and probability values. The meta-analysis identified a significant association at chromosome 8q13 (rs10091374; P=9.06×10(-9)) near the nuclear receptor coactivator 2 (NCOA2) gene. Overexpression of NCOA2 can be detected in myoblasts. An additional analysis using logistic regression and the clinically motivated 99th percentile cut point detected a significant association at 1q32 (rs12564445; P=4.73×10(-8)) in the gene TNNT2, which encodes the cardiac troponin T protein itself. The hs-cTnT-associated single-nucleotide polymorphisms were not associated with coronary heart disease in a large case-control study, but rs12564445 was significantly associated with incident heart failure in Atherosclerosis Risk in Communities Study European Americans (hazard ratio=1.16; P=0.004).

CONCLUSIONS

We identified 2 loci, near NCOA2 and in the TNNT2 gene, at which variation was significantly associated with hs-cTnT levels. Further use of the new assay should enable replication of these results.

Novel locus including FGF21 is associated with dietary macronutrient intake

Dietary intake of macronutrients (carbohydrate, protein, and fat) has been associated with risk of chronic conditions such as obesity and diabetes. Family studies have reported a moderate contribution of genetics to variation in macronutrient intake. In a genome-wide meta-analysis of a population-based discovery cohort (n = 33 533), rs838133 in FGF21 (19q13.33), rs197273 near TRAF family member-associated NF-kappa-B activator (TANK) (2p24.2), and rs10163409 in FTO (16q12.2) were among the top associations (P < 10(-5)) for percentage of total caloric intake from protein and carbohydrate. rs838133 was replicated in silico in an independent sample from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium (CHARGE) Nutrition Working Group (n = 38 360) and attained genome-wide significance in combined analysis (Pjoint = 7.9 × 10(-9)). A cytokine involved in cellular metabolism, FGF21 is a potential susceptibility gene for obesity and type 2 diabetes. Our results highlight the potential of genetic variation for determining dietary macronutrient intake.

TATES: efficient multivariate genotype-phenotype analysis for genome-wide association studies

To date, the genome-wide association study (GWAS) is the primary tool to identify genetic variants that cause phenotypic variation. As GWAS analyses are generally univariate in nature, multivariate phenotypic information is usually reduced to a single composite score. This practice often results in loss of statistical power to detect causal variants. Multivariate genotype–phenotype methods do exist but attain maximal power only in special circumstances. Here, we present a new multivariate method that we refer to as TATES (Trait-based Association Test that uses Extended Simes procedure), inspired by the GATES procedure proposed by Li et al (2011). For each component of a multivariate trait, TATES combines p-values obtained in standard univariate GWAS to acquire one trait-based p-value, while correcting for correlations between components. Extensive simulations, probing a wide variety of genotype–phenotype models, show that TATES's false positive rate is correct, and that TATES's statistical power to detect causal variants explaining 0.5% of the variance can be 2.5–9 times higher than the power of univariate tests based on composite scores and 1.5–2 times higher than the power of the standard MANOVA. Unlike other multivariate methods, TATES detects both genetic variants that are common to multiple phenotypes and genetic variants that are specific to a single phenotype, i.e. TATES provides a more complete view of the genetic architecture of complex traits. As the actual causal genotype–phenotype model is usually unknown and probably phenotypically and genetically complex, TATES, available as an open source program, constitutes a powerful new multivariate strategy that allows researchers to identify novel causal variants, while the complexity of traits is no longer a limiting factor.

The genome-wide association study (GWAS) is the primary tool to identify genetic variants that cause phenotypic variation. As GWAS methods are generally univariate in nature, multivariate phenotypic information is usually reduced to a single composite score, which frequently results in a considerable loss of statistical power to detect causal variants. Multivariate genotype–phenotype methods do exist but attain maximal power only in special circumstances. We present a new multivariate method called TATES (Trait-based Association Test that uses Extended Simes procedure). Extensive simulations show that TATES's false positive rate is correct, and that TATES's statistical power to detect causal variants explaining 0.5% of the variance can be 2.5–9 times higher than the power of univariate tests of composite scores and 1.5–2 times higher than the power of the standard MANOVA. Unlike other multivariate methods, TATES uncovers both genetic variants that are common to multiple phenotypes as well as phenotype specific variants. TATES thus provides a more complete view of the genetic architecture of complex traits and constitutes a powerful new multivariate strategy that allows researchers to identify novel causal variants.