ProbABEL package for genome-wide association analysis of imputed data

Genome-wide association study identifies a common variant associated with risk of endometrial cancer

Endometrial cancer is the most common malignancy of the female genital tract in developed countries. To identify genetic variants associated with endometrial cancer risk, we performed a genome-wide association study involving 1,265 individuals with endometrial cancer (cases) from Australia and the UK and 5,190 controls from the Wellcome Trust Case Control Consortium. We compared genotype frequencies in cases and controls for 519,655 SNPs. Forty seven SNPs that showed evidence of association with endometrial cancer in stage 1 were genotyped in 3,957 additional cases and 6,886 controls. We identified an endometrial cancer susceptibility locus close to HNF1B at 17q12 (rs4430796, P = 7.1 × 10(-10)) that is also associated with risk of prostate cancer and is inversely associated with risk of type 2 diabetes.

PredictABEL: an R package for the assessment of risk prediction models

The rapid identification of genetic markers for multifactorial diseases from genome-wide association studies is fuelling interest in investigating the predictive ability and health care utility of genetic risk models. Various measures are available for the assessment of risk prediction models, each addressing a different aspect of performance and utility. We developed PredictABEL, a package in R that covers descriptive tables, measures and figures that are used in the analysis of risk prediction studies such as measures of model fit, predictive ability and clinical utility, and risk distributions, calibration plot and the receiver operating characteristic plot. Tables and figures are saved as separate files in a user-specified format, which include publication-quality EPS and TIFF formats. All figures are available in a ready-made layout, but they can be customized to the preferences of the user. The package has been developed for the analysis of genetic risk prediction studies, but can also be used for studies that only include non-genetic risk factors. PredictABEL is freely available at the websites of GenABEL ( http://www.genabel.org ) and CRAN ( http://cran.r-project.org/).

Genome-wide meta-analysis identifies regions on 7p21 (AHR) and 15q24 (CYP1A2) as determinants of habitual caffeine consumption

We report the first genome-wide association study of habitual caffeine intake. We included 47,341 individuals of European descent based on five population-based studies within the United States. In a meta-analysis adjusted for age, sex, smoking, and eigenvectors of population variation, two loci achieved genome-wide significance: 7p21 (P = 2.4 × 10(-19)), near AHR, and 15q24 (P = 5.2 × 10(-14)), between CYP1A1 and CYP1A2. Both the AHR and CYP1A2 genes are biologically plausible candidates as CYP1A2 metabolizes caffeine and AHR regulates CYP1A2.

Genome-wide association study identifies genetic variants influencing F-cell levels in sickle-cell patients

Fetal hemoglobin (HbF) level has emerged as an important prognostic factor in sickle-cell disease (SCD) and can be measured by the proportion of HbF-containing erythrocytes (F-cells). Recently, BCL11A (zinc-finger protein) was identified as a regulator of HbF, and the strongest association signals were observed either directly for rs766432 or for correlated single-nucleotide polymorphisms (SNPs). To identify additional independently associated genetic variants, we performed a genome-wide association study (GWAS) on the proportion of F-cells in individuals of African ancestry with SCD from the Silent Infarct Transfusion (SIT) Trial cohort. Our study not only confirms the association of rs766432 (P-value <3.32 × 10(-13)), but also identifies an independent novel intronic SNP, rs7606173, associated with F-cells (P-value <1.81 × 10(-15)). The F-cell variances explained independently by these two SNPs are ∼13% (rs7606173) and ∼11% (rs766432), whereas, together they explain ∼16%. Additionally, in men, we identify a novel locus on chromosome 17, glucagon-like peptide-2 receptor (GLP2R), associated with F-cell regulation (rs12103880; P-value <3.41 × 10(-8)). GLP2R encodes a G protein-coupled receptor and involved in proliferative and anti-apoptotic cellular responses. These findings highlight the importance of denser genetic screens and suggest further exploration of the BCL11A and GLP2R loci to gain additional insight into HbF/F-cell regulation.

A genome-wide association study of upper aerodigestive tract cancers conducted within the INHANCE consortium

Genome-wide association studies (GWAS) have been successful in identifying common genetic variation involved in susceptibility to etiologically complex disease. We conducted a GWAS to identify common genetic variation involved in susceptibility to upper aero-digestive tract (UADT) cancers. Genome-wide genotyping was carried out using the Illumina HumanHap300 beadchips in 2,091 UADT cancer cases and 3,513 controls from two large European multi-centre UADT cancer studies, as well as 4,821 generic controls. The 19 top-ranked variants were investigated further in an additional 6,514 UADT cancer cases and 7,892 controls of European descent from an additional 13 UADT cancer studies participating in the INHANCE consortium. Five common variants presented evidence for significant association in the combined analysis (p ≤ 5 × 10⁻⁷). Two novel variants were identified, a 4q21 variant (rs1494961, p = 1×10⁻⁸) located near DNA repair related genes HEL308 and FAM175A (or Abraxas) and a 12q24 variant (rs4767364, p =2 × 10⁻⁸) located in an extended linkage disequilibrium region that contains multiple genes including the aldehyde dehydrogenase 2 (ALDH2) gene. Three remaining variants are located in the ADH gene cluster and were identified previously in a candidate gene study involving some of these samples. The association between these three variants and UADT cancers was independently replicated in 5,092 UADT cancer cases and 6,794 controls non-overlapping samples presented here (rs1573496-ADH7, p = 5 × 10⁻⁸); rs1229984-ADH1B, p = 7 × 10⁻⁹; and rs698-ADH1C, p = 0.02). These results implicate two variants at 4q21 and 12q24 and further highlight three ADH variants in UADT cancer susceptibility.

Statistical Optimization of Pharmacogenomics Association Studies: Key Considerations from Study Design to Analysis

Research in human genetics and genetic epidemiology has grown significantly over the previous decade, particularly in the field of pharmacogenomics. Pharmacogenomics presents an opportunity for rapid translation of associated genetic polymorphisms into diagnostic measures or tests to guide therapy as part of a move towards personalized medicine. Expansion in genotyping technology has cleared the way for widespread use of whole-genome genotyping in the effort to identify novel biology and new genetic markers associated with pharmacokinetic and pharmacodynamic endpoints. With new technology and methodology regularly becoming available for use in genetic studies, a discussion on the application of such tools becomes necessary. In particular, quality control criteria have evolved with the use of GWAS as we have come to understand potential systematic errors which can be introduced into the data during genotyping. There have been several replicated pharmacogenomic associations, some of which have moved to the clinic to enact change in treatment decisions. These examples of translation illustrate the strength of evidence necessary to successfully and effectively translate a genetic discovery. In this review, the design of pharmacogenomic association studies is examined with the goal of optimizing the impact and utility of this research. Issues of ascertainment, genotyping, quality control, analysis and interpretation are considered.

Genome-wide meta-analysis of joint tests for genetic and gene-environment interaction effects

BACKGROUND

There is growing interest in the study of gene-environment interactions in the context of genome-wide association studies (GWASs). These studies will likely require meta-analytic approaches to have sufficient power.

METHODS

We describe an approach for meta-analysis of a joint test for genetic main effects and gene-environment interaction effects. Using simulation studies based on a meta-analysis of five studies (total n = 10,161), we compare the power of this test to the meta-analysis of marginal test of genetic association and the meta-analysis of standard 1 d.f. interaction tests across a broad range of genetic main effects and gene-environment interaction effects.

RESULTS

We show that the joint meta-analysis is valid and can be more powerful than classical meta-analytic approaches, with a potential gain of power over 50% compared to the marginal test. The standard interaction test had less than 1% power in almost all the situations we considered. We also show that regardless of the test used, sample sizes far exceeding those of a typical individual GWAS will be needed to reliably detect genes with subtle gene-environment interaction patterns.

CONCLUSION

The joint meta-analysis is an attractive approach to discover markers which may have been missed by initial GWASs focusing on marginal marker-trait associations.

Association of HSP70 and its co-chaperones with Alzheimer's disease

The heat shock protein (HSP) 70 family has been implicated in the pathology of Alzheimer's disease (AD). In this study, we examined common genetic variations in the 80 genes encoding HSP70 and its co-chaperones. We conducted a study in a series of 462 patients and 5238 unaffected participants derived from the Rotterdam Study, a population-based study including 7983 persons aged 55 years and older. We genotyped a total of 12,053 Single Nucleotide Polymorphisms (SNPs) using the HumanHap550K Genotyping BeadChip from Illumina. Replication was performed in two independent cohort studies, the Framingham Heart study (FHS; n = 806) and Cardiovascular Health Study (CHS; n = 2150). When adjusting for multiple testing, we found a small but consistent, though not significant effect of rs12118313 located 32 kb from PFDN2, with an OR of 1.19 (p-value from meta-analysis = 0.003). However this SNP was in the intron of another gene, suggesting it is unlikely this SNP reflects the effect of PFDN2. In a formal pathway analysis we found nominally significant evidence for an association of BAG, DNAJA and prefoldin with AD. These findings corroborate with those of a study of 2032 AD patients and 5328 controls, in which several members of the prefoldin family showed evidence for association to AD. Our study did not reveal evidence for a genetic variant if the HSP70 family with a major effect on AD. However, our findings of the single SNP analysis and pathway analysis suggest that multiple genetic variants in prefoldin are associated with AD.

Practical Consideration of Genotype Imputation: Sample Size, Window Size, Reference Choice, and Untyped Rate

Imputation offers a promising way to infer the missing and/or untyped genotypes in genetic studies. In practice, however, many factors may affect the quality of imputation. In this study, we evaluated the influence of untyped rate, sizes of the study sample and the reference sample, window size, and reference choice (for admixed population), as the factors affecting the quality of imputation. The results show that in order to obtain good imputation quality, it is necessary to have an untyped rate less than 50%, a reference sample size greater than 50, and a window size of greater than 500 SNPs (roughly 1 MB in base pairs). Compared with the whole-region imputation, piecewise imputation with large-enough window sizes provides improved efficacy. For an admixed study sample, if only an external reference panel is used, it should include samples from the ancestral populations that represent the admixed population under investigation. Internal references are strongly recommended. When internal references are limited, however, augmentation by external references should be used carefully. More specifically, augmentation with samples from the major source populations of the admixture can lower the quality of imputation; augmentation with seemingly genetically unrelated cohorts may improve the quality of imputation.

Genome-wide association study of renal cell carcinoma identifies two susceptibility loci on 2p21 and 11q13.3

We conducted a two-stage genome-wide association study of renal cell carcinoma (RCC) in 3,772 affected individuals (cases) and 8,505 controls of European background from 11 studies and followed up 6 SNPs in 3 replication studies of 2,198 cases and 4,918 controls. Two loci on the regions of 2p21 and 11q13.3 were associated with RCC susceptibility below genome-wide significance. Two correlated variants (r² = 0.99 in controls), rs11894252 (P = 1.8 × 10⁻⁸) and rs7579899 (P = 2.3 × 10⁻⁹), map to EPAS1 on 2p21, which encodes hypoxia-inducible-factor-2 alpha, a transcription factor previously implicated in RCC. The second locus, rs7105934, at 11q13.3, contains no characterized genes (P = 7.8 × 10⁻¹⁴). In addition, we observed a promising association on 12q24.31 for rs4765623, which maps to SCARB1, the scavenger receptor class B, member 1 gene (P = 2.6 × 10⁻⁸). Our study reports previously unidentified genomic regions associated with RCC risk that may lead to new etiological insights.

Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index

Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and ∼ 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 × 10⁻⁸), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation.

Identification of Diabetic Retinopathy Genes through a Genome-Wide Association Study among Mexican-Americans from Starr County, Texas

To identify genetic loci for severe diabetic retinopathy, 286 Mexican-Americans with type 2 diabetes from Starr County, Texas, completed physical examinations including fundus photography for diabetic retinopathy grading. Individuals with moderate-to-severe non-proliferative and proliferative diabetic retinopathy were defined as cases. Direct genotyping was performed using the Affymetrix GeneChip Human Mapping 100 K Set, and SNPs passing quality control criteria were used to impute markers available in HapMap Phase III Mexican population (MXL) in Los Angeles, California. Two directly genotyped markers were associated with severe diabetic retinopathy at a P-value less than .0001: SNP rs2300782 (P = 6.04 × 10(-5)) mapped to an intron region of CAMK4 (calcium/calmodulin-dependent protein kinase IV) on chromosome 5, and SNP rs10519765 (P = 6.21 × 10(-5)) on chromosomal 15q13 in the FMN1 (formin 1) gene. Using well-imputed markers based on the HapMap III Mexican population, we identified an additional 32 SNPs located in 11 chromosomal regions with nominal association with severe diabetic retinopathy at P-value less than .0001. None of these markers were located in traditional candidate genes for diabetic retinopathy or diabetes itself. However, these signals implicate genes involved in inflammation, oxidative stress and cell adhesion for the development and progression of diabetic retinopathy.

A large-scale candidate gene association study of age at menarche and age at natural menopause

Recent genome-wide association (GWA) studies have identified several novel genetic loci associated with age at menarche and age at natural menopause. However, the stringent significance threshold used in GWA studies potentially led to false negatives and true associations may have been overlooked. Incorporating biologically relevant information, we examined whether common genetic polymorphisms in candidate genes of nine groups of biologically plausible pathways and related phenotypes are associated with age at menarche and age at natural menopause. A total of 18,862 genotyped and imputed single nucleotide polymorphisms (SNPs) in 278 genes were assessed for their associations with these two traits among a total of 24,341 women from the Nurses' Health Study (NHS, N = 2,287) and the Women's Genome Health Study (WGHS, N = 22,054). Linear regression was used to assess the marginal association of each SNP with each phenotype. We adjusted for multiple testing within each gene to identify statistically significant SNP associations at the gene level. To evaluate the overall evidence for an excess of statistically significant gene associations over the proportion expected by chance, we applied a one-sample test of proportion to each group of candidate genes. The steroid-hormone metabolism and biosynthesis pathway was found significantly associated with both age at menarche and age at natural menopause (P = 0.040 and 0.011, respectively). In addition, the group of genes associated with precocious or delayed puberty was found significantly associated with age at menarche (P = 0.013), and the group of genes involved in premature ovarian failure with age at menopause (P = 0.025).

New loci associated with central cornea thickness include COL5A1, AKAP13 and AVGR8

Central corneal thickness (CCT) is a highly heritable trait, which has been proposed to influence disorders of the anterior segment of the eye. A genome-wide association study (GWAS) of CCT was performed in 2269 individuals from three Croatian and one Scottish population. In the discovery set (1445 individuals), two genome-wide significant associations were identified for single nucleotide polymorphisms rs12447690 (β = 0.23 SD, P = 4.4 × 10(-9)) and rs1536482 (β = 0.22 SD, P = 7.1 × 10(-8)) for which the closest candidate genes (although ≥90 kb away) were zinc finger 469 (ZNF469) on 16q24.2 and collagen 5 alpha 1 (COL5A1) on 9q34.2, respectively. Only the ZNF469 association was confirmed in our replication set (824 individuals, P = 8.0 × 10(-4)) but COL5A1 remained a suggestive association in the combined sample (β = 0.16 SD, P = 1.1 × 10(-6)). Following a larger meta-analysis including recently published CCT GWAS summary data, COL5A1 was genome-wide significant (β = 0.13 SD, P = 5.1 × 10(-8)), together with two additional novel loci. The second new locus (defined by rs1034200) was 5 kb from the AVGR8 gene, encoding a putative transcription factor with typical ZNF and KRAB domains, in chromosomal region 13q12.11 (β = 0.14 SD, P = 3.5 × 10(-9)). The third new locus (rs6496932), on 15q25.3 (β = 0.13, P = 1.4 × 10(-8)), was within a wide linkage disequilibrium block extending into the 5' end of the AKAP13 gene, encoding a scaffold protein concerned with signal transduction from the cell surface. These associations offer mechanistic insights into the regulation of CCT and offer new candidate genes for susceptibility to common disorders in which CCT has been implicated, including primary open-angle glaucoma and keratoconus.

Common variants in the calcium-sensing receptor gene are associated with total serum calcium levels

Serum calcium levels are tightly regulated. We performed genome-wide association studies (GWAS) in population-based studies participating in the CHARGE Consortium to uncover common genetic variations associated with total serum calcium levels. GWAS of serum calcium concentrations was performed in 20 611 individuals of European ancestry for ∼2.5 million genotyped and imputed single-nucleotide polymorphisms (SNPs). The SNP with the lowest P-value was rs17251221 (P = 2.4 * 10(-22), minor allele frequency 14%) in the calcium-sensing receptor gene (CASR). This lead SNP was associated with higher serum calcium levels [0.06 mg/dl (0.015 mmol/l) per copy of the minor G allele] and accounted for 0.54% of the variance in serum calcium concentrations. The identification of variation in CASR that influences serum calcium concentration confirms the results of earlier candidate gene studies. The G allele of rs17251221 was also associated with higher serum magnesium levels (P = 1.2 * 10(-3)), lower serum phosphate levels (P = 2.8 * 10(-7)) and lower bone mineral density at the lumbar spine (P = 0.038), but not the femoral neck. No additional genomic loci contained SNPs associated at genome-wide significance (P < 5 * 10(-8)). These associations resemble clinical characteristics of patients with familial hypocalciuric hypercalcemia, an autosomal-dominant disease arising from rare inactivating mutations in the CASR gene. We conclude that common genetic variation in the CASR gene is associated with similar but milder features in the general population.