METAL: fast and efficient meta-analysis of genomewide association scans

Association of genetic loci with sleep apnea in European Americans and African-Americans: the Candidate Gene Association Resource (CARe)

Although obstructive sleep apnea (OSA) is known to have a strong familial basis, no genetic polymorphisms influencing apnea risk have been identified in cross-cohort analyses. We utilized the National Heart, Lung, and Blood Institute (NHLBI) Candidate Gene Association Resource (CARe) to identify sleep apnea susceptibility loci. Using a panel of 46,449 polymorphisms from roughly 2,100 candidate genes on a customized Illumina iSelect chip, we tested for association with the apnea hypopnea index (AHI) as well as moderate to severe OSA (AHI≥15) in 3,551 participants of the Cleveland Family Study and two cohorts participating in the Sleep Heart Health Study.

Among 647 African-Americans, rs11126184 in the pleckstrin (PLEK) gene was associated with OSA while rs7030789 in the lysophosphatidic acid receptor 1 (LPAR1) gene was associated with AHI using a chip-wide significance threshold of p-value<2×10⁻⁶. Among 2,904 individuals of European ancestry, rs1409986 in the prostaglandin E2 receptor (PTGER3) gene was significantly associated with OSA. Consistency of effects between rs7030789 and rs1409986 in LPAR1 and PTGER3 and apnea phenotypes were observed in independent clinic-based cohorts.

Novel genetic loci for apnea phenotypes were identified through the use of customized gene chips and meta-analyses of cohort data with replication in clinic-based samples. The identified SNPs all lie in genes associated with inflammation suggesting inflammation may play a role in OSA pathogenesis.

A meta-analysis of genome-wide association studies for serum total IgE in diverse study populations

BACKGROUND

IgE is both a marker and mediator of allergic inflammation. Despite reported differences in serum total IgE levels by race-ethnicity, African American and Latino subjects have not been well represented in genetic studies of total IgE.

OBJECTIVE

We sought to identify the genetic predictors of serum total IgE levels.

METHODS

We used genome-wide association data from 4292 subjects (2469 African Americans, 1564 European Americans, and 259 Latinos) in the EVE Asthma Genetics Consortium. Tests for association were performed within each cohort by race-ethnic group (ie, African American, Latino, and European American) and asthma status. The resulting P values were meta-analyzed, accounting for sample size and direction of effect. Top single nucleotide polymorphism associations from the meta-analysis were reassessed in 6 additional cohorts comprising 5767 subjects.

RESULTS

We identified 10 unique regions in which the combined association statistic was associated with total serum IgE levels (P<5.0×10(-6)) and the minor allele frequency was 5% or greater in 2 or more population groups. Variant rs9469220, corresponding to HLA-DQB1, was the single nucleotide polymorphism most significantly associated with serum total IgE levels when assessed in both the replication cohorts and the discovery and replication sets combined (P=.007 and 2.45×10(-7), respectively). In addition, findings from earlier genome-wide association studies were also validated in the current meta-analysis.

CONCLUSION

This meta-analysis independently identified a variant near HLA-DQB1 as a predictor of total serum IgE levels in multiple race-ethnic groups. This study also extends and confirms the findings of earlier genome-wide association analyses in African American and Latino subjects.

Population-based case-control association studies

This unit provides an overview of the design and analysis of population-based case-control studies of genetic risk factors for complex disease. Considerations specific to genetic studies are emphasized. The unit reviews basic study designs differentiating case-control studies from others, presents different genetic association strategies (candidate gene, genome-wide association, and high-throughput sequencing), introduces basic methods of statistical analysis for case-control data and approaches to combining case-control studies, and discusses measures of association and impact. Admixed populations, controlling for confounding (including population stratification), consideration of multiple loci and environmental risk factors, and complementary analyses of haplotypes, genes, and pathways are briefly discussed. Readers are referred to basic texts on epidemiology for more details on general conduct of case-control studies.

Imputation of exome sequence variants into population- based samples and blood-cell-trait-associated loci in African Americans: NHLBI GO Exome Sequencing Project

Researchers have successfully applied exome sequencing to discover causal variants in selected individuals with familial, highly penetrant disorders. We demonstrate the utility of exome sequencing followed by imputation for discovering low-frequency variants associated with complex quantitative traits. We performed exome sequencing in a reference panel of 761 African Americans and then imputed newly discovered variants into a larger sample of more than 13,000 African Americans for association testing with the blood cell traits hemoglobin, hematocrit, white blood count, and platelet count. First, we illustrate the feasibility of our approach by demonstrating genome-wide-significant associations for variants that are not covered by conventional genotyping arrays; for example, one such association is that between higher platelet count and an MPL c.117G>T (p.Lys39Asn) variant encoding a p.Lys39Asn amino acid substitution of the thrombopoietin receptor gene (p = 1.5 × 10(-11)). Second, we identified an association between missense variants of LCT and higher white blood count (p = 4 × 10(-13)). Third, we identified low-frequency coding variants that might account for allelic heterogeneity at several known blood cell-associated loci: MPL c.754T>C (p.Tyr252His) was associated with higher platelet count; CD36 c.975T>G (p.Tyr325(∗)) was associated with lower platelet count; and several missense variants at the α-globin gene locus were associated with lower hemoglobin. By identifying low-frequency missense variants associated with blood cell traits not previously reported by genome-wide association studies, we establish that exome sequencing followed by imputation is a powerful approach to dissecting complex, genetically heterogeneous traits in large population-based studies.

Large-scale gene-centric meta-analysis across 32 studies identifies multiple lipid loci

Genome-wide association studies (GWASs) have identified many SNPs underlying variations in plasma-lipid levels. We explore whether additional loci associated with plasma-lipid phenotypes, such as high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), total cholesterol (TC), and triglycerides (TGs), can be identified by a dense gene-centric approach. Our meta-analysis of 32 studies in 66,240 individuals of European ancestry was based on the custom ∼50,000 SNP genotyping array (the ITMAT-Broad-CARe array) covering ∼2,000 candidate genes. SNP-lipid associations were replicated either in a cohort comprising an additional 24,736 samples or within the Global Lipid Genetic Consortium. We identified four, six, ten, and four unreported SNPs in established lipid genes for HDL-C, LDL-C, TC, and TGs, respectively. We also identified several lipid-related SNPs in previously unreported genes: DGAT2, HCAR2, GPIHBP1, PPARG, and FTO for HDL-C; SOCS3, APOH, SPTY2D1, BRCA2, and VLDLR for LDL-C; SOCS3, UGT1A1, BRCA2, UBE3B, FCGR2A, CHUK, and INSIG2 for TC; and SERPINF2, C4B, GCK, GATA4, INSR, and LPAL2 for TGs. The proportion of explained phenotypic variance in the subset of studies providing individual-level data was 9.9% for HDL-C, 9.5% for LDL-C, 10.3% for TC, and 8.0% for TGs. This large meta-analysis of lipid phenotypes with the use of a dense gene-centric approach identified multiple SNPs not previously described in established lipid genes and several previously unknown loci. The explained phenotypic variance from this approach was comparable to that from a meta-analysis of GWAS data, suggesting that a focused genotyping approach can further increase the understanding of heritability of plasma lipids.

Genetic risk factors for ischaemic stroke and its subtypes (the METASTROKE collaboration): a meta-analysis of genome-wide association studies

BACKGROUND

Various genome-wide association studies (GWAS) have been done in ischaemic stroke, identifying a few loci associated with the disease, but sample sizes have been 3500 cases or less. We established the METASTROKE collaboration with the aim of validating associations from previous GWAS and identifying novel genetic associations through meta-analysis of GWAS datasets for ischaemic stroke and its subtypes.

METHODS

We meta-analysed data from 15 ischaemic stroke cohorts with a total of 12 389 individuals with ischaemic stroke and 62 004 controls, all of European ancestry. For the associations reaching genome-wide significance in METASTROKE, we did a further analysis, conditioning on the lead single nucleotide polymorphism in every associated region. Replication of novel suggestive signals was done in 13 347 cases and 29 083 controls.

FINDINGS

We verified previous associations for cardioembolic stroke near PITX2 (p=2·8×10(-16)) and ZFHX3 (p=2·28×10(-8)), and for large-vessel stroke at a 9p21 locus (p=3·32×10(-5)) and HDAC9 (p=2·03×10(-12)). Additionally, we verified that all associations were subtype specific. Conditional analysis in the three regions for which the associations reached genome-wide significance (PITX2, ZFHX3, and HDAC9) indicated that all the signal in each region could be attributed to one risk haplotype. We also identified 12 potentially novel loci at p<5×10(-6). However, we were unable to replicate any of these novel associations in the replication cohort.

INTERPRETATION

Our results show that, although genetic variants can be detected in patients with ischaemic stroke when compared with controls, all associations we were able to confirm are specific to a stroke subtype. This finding has two implications. First, to maximise success of genetic studies in ischaemic stroke, detailed stroke subtyping is required. Second, different genetic pathophysiological mechanisms seem to be associated with different stroke subtypes.

FUNDING

Wellcome Trust, UK Medical Research Council (MRC), Australian National and Medical Health Research Council, National Institutes of Health (NIH) including National Heart, Lung and Blood Institute (NHLBI), the National Institute on Aging (NIA), the National Human Genome Research Institute (NHGRI), and the National Institute of Neurological Disorders and Stroke (NINDS).

Transferability and fine-mapping of glucose and insulin quantitative trait loci across populations: CARe, the Candidate Gene Association Resource

AIMS/HYPOTHESIS

Hyperglycaemia disproportionately affects African-Americans (AfAs). We tested the transferability of 18 single-nucleotide polymorphisms (SNPs) associated with glycaemic traits identified in European ancestry (EuA) populations in 5,984 non-diabetic AfAs.

METHODS

We meta-analysed SNP associations with fasting glucose (FG) or insulin (FI) in AfAs from five cohorts in the Candidate Gene Association Resource. We: (1) calculated allele frequency differences, variations in linkage disequilibrium (LD), fixation indices (F(st)s) and integrated haplotype scores (iHSs); (2) tested EuA SNPs in AfAs; and (3) interrogated within ± 250 kb around each EuA SNP in AfAs.

RESULTS

Allele frequency differences ranged from 0.6% to 54%. F(st) exceeded 0.15 at 6/16 loci, indicating modest population differentiation. All iHSs were <2, suggesting no recent positive selection. For 18 SNPs, all directions of effect were the same and 95% CIs of association overlapped when comparing EuA with AfA. For 17 of 18 loci, at least one SNP was nominally associated with FG in AfAs. Four loci were significantly associated with FG (GCK, p = 5.8 × 10(-8); MTNR1B, p = 8.5 × 10(-9); and FADS1, p = 2.2 × 10(-4)) or FI (GCKR, p = 5.9 × 10(-4)). At GCK and MTNR1B the EuA and AfA SNPs represented the same signal, while at FADS1, and GCKR, the EuA and best AfA SNPs were weakly correlated (r(2) <0.2), suggesting allelic heterogeneity for association with FG at these loci.

CONCLUSIONS/INTERPRETATION

Few glycaemic SNPs showed strict evidence of transferability from EuA to AfAs. Four loci were significantly associated in both AfAs and those with EuA after accounting for varying LD across ancestral groups, with new signals emerging to aid fine-mapping.

Genetic association analysis of complex diseases incorporating intermediate phenotype information

Genetic researchers often collect disease related quantitative traits in addition to disease status because they are interested in understanding the pathophysiology of disease processes. In genome-wide association (GWA) studies, these quantitative phenotypes may be relevant to disease development and serve as intermediate phenotypes or they could be behavioral or other risk factors that predict disease risk. Statistical tests combining both disease status and quantitative risk factors should be more powerful than case-control studies, as the former incorporates more information about the disease. In this paper, we proposed a modified inverse-variance weighted meta-analysis method to combine disease status and quantitative intermediate phenotype information. The simulation results showed that when an intermediate phenotype was available, the inverse-variance weighted method had more power than did a case-control study of complex diseases, especially in identifying susceptibility loci having minor effects. We further applied this modified meta-analysis to a study of imputed lung cancer genotypes with smoking data in 1154 cases and 1137 matched controls. The most significant SNPs came from the CHRNA3-CHRNA5-CHRNB4 region on chromosome 15q24–25.1, which has been replicated in many other studies. Our results confirm that this CHRNA region is associated with both lung cancer development and smoking behavior. We also detected three significant SNPs—rs1800469, rs1982072, and rs2241714—in the promoter region of the TGFB1 gene on chromosome 19 (p = 1.46×10⁻⁵, 1.18×10⁻⁵, and 6.57×10⁻⁶, respectively). The SNP rs1800469 is reported to be associated with chronic obstructive pulmonary disease and lung cancer in cigarette smokers. The present study is the first GWA study to replicate this result. Signals in the 3q26 region were also identified in the meta-analysis. We demonstrate the intermediate phenotype can potentially enhance the power of complex disease association analysis and the modified meta-analysis method is robust to incorporate intermediate phenotype or other quantitative risk factor in the analysis.

Exome sequencing and complex disease: practical aspects of rare variant association studies

Genetic association and linkage studies can provide insights into complex disease biology, guiding the development of new diagnostic and therapeutic strategies. Over the past decade, genetic association studies have largely focused on common, easy to measure genetic variants shared between many individuals. These common variants typically have subtle functional consequence and translating the resulting association signals into biological insights can be challenging. In the last few years, exome sequencing has emerged as a cost-effective strategy for extending these studies to include rare coding variants, which often have more marked functional consequences. Here, we provide practical guidance in the design and analysis of complex trait association studies focused on rare, coding variants.

FTO genotype is associated with phenotypic variability of body mass index

There is evidence across several species for genetic control of phenotypic variation of complex traits, such that the variance among phenotypes is genotype dependent. Understanding genetic control of variability is important in evolutionary biology, agricultural selection programmes and human medicine, yet for complex traits, no individual genetic variants associated with variance, as opposed to the mean, have been identified. Here we perform a meta-analysis of genome-wide association studies of phenotypic variation using ∼170,000 samples on height and body mass index (BMI) in human populations. We report evidence that the single nucleotide polymorphism (SNP) rs7202116 at the FTO gene locus, which is known to be associated with obesity (as measured by mean BMI for each rs7202116 genotype), is also associated with phenotypic variability. We show that the results are not due to scale effects or other artefacts, and find no other experiment-wise significant evidence for effects on variability, either at loci other than FTO for BMI or at any locus for height. The difference in variance for BMI among individuals with opposite homozygous genotypes at the FTO locus is approximately 7%, corresponding to a difference of ∼0.5 kilograms in the standard deviation of weight. Our results indicate that genetic variants can be discovered that are associated with variability, and that between-person variability in obesity can partly be explained by the genotype at the FTO locus. The results are consistent with reported FTO by environment interactions for BMI, possibly mediated by DNA methylation. Our BMI results for other SNPs and our height results for all SNPs suggest that most genetic variants, including those that influence mean height or mean BMI, are not associated with phenotypic variance, or that their effects on variability are too small to detect even with samples sizes greater than 100,000.

Genetic variants near PDGFRA are associated with corneal curvature in Australians

PURPOSE

Irregularity in the corneal curvature (CC) is highly associated with various eye disorders such as keratoconus and myopia. The sample had limited power to find genomewide significant (5 × 10(-8)) hits but good power for replication. Thus, an attempt was made to test whether alleles in the FRAP1 and PDGFRA genes, recently found to be associated with CC in Asian populations, also influence CC in Australians of North European ancestry. Results of initial genomewide association studies (GWAS) for CC in Australians were also reported.

METHODS

Two population-based cohorts of 1788 Australian twins and their families, as well as 1013 individuals from a birth cohort from Western Australia, were genotyped using genomewide arrays. Following separate individual analysis and quality control, the results from each cohort underwent meta-analysis.

RESULTS

Meta-analysis revealed significant replication of association between rs2114039 and corneal curvature (P = 0.0045). The SNP rs2114039 near PDGFRA has been previously implicated in Asians. No SNP at the FRAP1 locus was found to be associated in our Australian samples. No SNP surpassed the genomewide significance threshold of 5 × 10(-8). The SNP with strongest association was rs2444240 (P = 3.658 × 10(-7)), which is 31 kb upstream to the TRIM29 gene.

CONCLUSIONS

A significant role of the PDGFRA gene in determining corneal curvature in the Australian population was confirmed in this study, also highlighting the putative association of the TRIM29 locus with CC.

Whole-genome studies identify solute carrier transporters in cellular susceptibility to paclitaxel

OBJECTIVE

The clinical use of paclitaxel is limited by variable responses and the potential for significant toxicity. To date, studies of associations between variants in candidate genes and paclitaxel effects have yielded conflicting results. We aimed to evaluate the relationships between global gene expression and paclitaxel sensitivity.

METHODS

We utilized well-genotyped lymphoblastoid cell lines derived from the International HapMap Project to evaluate the relationships between cellular susceptibility to paclitaxel and global gene expression. Cells were exposed to varying concentrations of paclitaxel to evaluate paclitaxel-induced cytotoxicity and apoptosis. Among the top genes, we identified solute carrier (SLC) genes associated with paclitaxel sensitivity and narrowed down the list to those that had single nucleotide polymorphisms associated with both the expression level of the SLC gene and also with paclitaxel sensitivity. We performed an independent validation in an independent set of cell lines and also conducted functional studies using RNA interference.

RESULTS

Of all genes associated with paclitaxel-induced cytotoxicity at P less than 0.05 (1713 genes), there was a significant enrichment in SLC genes (31 genes). A subset of SLC genes, namely SLC31A2, SLC43A1, SLC35A5, and SLC41A2, was associated with paclitaxel sensitivity and had regulating single nucleotide polymorphisms that were also associated with paclitaxel-induced cytotoxicity. Multivariate modeling demonstrated that those four SLC genes together explained 20% of the observed variability in paclitaxel susceptibility. Using RNA interference, we demonstrated increased paclitaxel susceptibility with knockdown of three SLC genes, SLC31A2, SLC35A5, and SLC41A2.

CONCLUSION

Our findings are novel and lend further support to the role of transporters, specifically solute carriers, in mediating cellular susceptibility to paclitaxel.

CSK regulatory polymorphism is associated with systemic lupus erythematosus and influences B-cell signaling and activation

C-src tyrosine kinase, Csk, physically interacts with the intracellular phosphatase Lyp (PTPN22) and can modify the activation state of downstream Src kinases, such as Lyn, in lymphocytes. We identified an association of Csk with systemic lupus erythematosus (SLE) and refined its location to an intronic polymorphism rs34933034 (OR 1.32, p = 1.04 × 10⁻⁹). The risk allele is associated with increased CSK expression and augments inhibitory phosphorylation of Lyn. In carriers of the risk allele, B cell receptor (BCR)-mediated activation of mature B cells, as well as plasma IgM, are increased. Moreover, the fraction of transitional B cells is doubled in the cord blood of carriers of the risk allele compared to non-risk haplotypes due to an expansion of the late transitional cells, a stage targeted by selection mechanisms. This suggests that the Lyp-Csk complex increases susceptibility to lupus at multiple maturation and activation points of B cells.

Insights into pancreatic cancer etiology from pathway analysis of genome-wide association study data

BACKGROUND

Pancreatic cancer is the fourth leading cause of cancer death in the U.S. and the etiology of this highly lethal disease has not been well defined. To identify genetic susceptibility factors for pancreatic cancer, we conducted pathway analysis of genome-wide association study (GWAS) data in 3,141 pancreatic cancer patients and 3,367 controls with European ancestry.

METHODS

Using the gene set ridge regression in association studies (GRASS) method, we analyzed 197 pathways identified from the Kyoto Encyclopedia of Genes and Genomes database. We used the logistic kernel machine (LKM) test to identify major contributing genes to each pathway. We conducted functional enrichment analysis of the most significant genes (P<0.01) using the Database for Annotation, Visualization, and Integrated Discovery (DAVID).

RESULTS

Two pathways were significantly associated with risk of pancreatic cancer after adjusting for multiple comparisons (P<0.00025) and in replication testing: neuroactive ligand-receptor interaction, (Ps<0.00002), and the olfactory transduction pathway (P = 0.0001). LKM test identified four genes that were significantly associated with risk of pancreatic cancer after Bonferroni correction (P<1×10(-5)): ABO, HNF1A, OR13C4, and SHH. Functional enrichment analysis using DAVID consistently found the G protein-coupled receptor signaling pathway (including both neuroactive ligand-receptor interaction and olfactory transduction pathways) to be the most significant pathway for pancreatic cancer risk in this study population.

CONCLUSION

These novel findings provide new perspectives on genetic susceptibility to and molecular mechanisms of pancreatic cancer.

Genome-wide association studies identify CHRNA5/3 and HTR4 in the development of airflow obstruction

RATIONALE

Genome-wide association studies (GWAS) have identified loci influencing lung function, but fewer genes influencing chronic obstructive pulmonary disease (COPD) are known.

OBJECTIVES

Perform meta-analyses of GWAS for airflow obstruction, a key pathophysiologic characteristic of COPD assessed by spirometry, in population-based cohorts examining all participants, ever smokers, never smokers, asthma-free participants, and more severe cases.

METHODS

Fifteen cohorts were studied for discovery (3,368 affected; 29,507 unaffected), and a population-based family study and a meta-analysis of case-control studies were used for replication and regional follow-up (3,837 cases; 4,479 control subjects). Airflow obstruction was defined as FEV(1) and its ratio to FVC (FEV(1)/FVC) both less than their respective lower limits of normal as determined by published reference equations.

MEASUREMENTS AND MAIN RESULTS

The discovery meta-analyses identified one region on chromosome 15q25.1 meeting genome-wide significance in ever smokers that includes AGPHD1, IREB2, and CHRNA5/CHRNA3 genes. The region was also modestly associated among never smokers. Gene expression studies confirmed the presence of CHRNA5/3 in lung, airway smooth muscle, and bronchial epithelial cells. A single-nucleotide polymorphism in HTR4, a gene previously related to FEV(1)/FVC, achieved genome-wide statistical significance in combined meta-analysis. Top single-nucleotide polymorphisms in ADAM19, RARB, PPAP2B, and ADAMTS19 were nominally replicated in the COPD meta-analysis.

CONCLUSIONS

These results suggest an important role for the CHRNA5/3 region as a genetic risk factor for airflow obstruction that may be independent of smoking and implicate the HTR4 gene in the etiology of airflow obstruction.

Genome-wide association studies of asthma in population-based cohorts confirm known and suggested loci and identify an additional association near HLA

RATIONALE

Asthma has substantial morbidity and mortality and a strong genetic component, but identification of genetic risk factors is limited by availability of suitable studies.

OBJECTIVES

To test if population-based cohorts with self-reported physician-diagnosed asthma and genome-wide association (GWA) data could be used to validate known associations with asthma and identify novel associations.

METHODS

The APCAT (Analysis in Population-based Cohorts of Asthma Traits) consortium consists of 1,716 individuals with asthma and 16,888 healthy controls from six European-descent population-based cohorts. We examined associations in APCAT of thirteen variants previously reported as genome-wide significant (P<5 x 10(-8)) and three variants reported as suggestive (P<5× 10(-7)). We also searched for novel associations in APCAT (Stage 1) and followed-up the most promising variants in 4,035 asthmatics and 11,251 healthy controls (Stage 2). Finally, we conducted the first genome-wide screen for interactions with smoking or hay fever.

MAIN RESULTS

We observed association in the same direction for all thirteen previously reported variants and nominally replicated ten of them. One variant that was previously suggestive, rs11071559 in RORA, now reaches genome-wide significance when combined with our data (P = 2.4 × 10(-9)). We also identified two genome-wide significant associations: rs13408661 near IL1RL1/IL18R1 (P(Stage1+Stage2) = 1.1x10(-9)), which is correlated with a variant recently shown to be associated with asthma (rs3771180), and rs9268516 in the HLA region (P(Stage1+Stage2) = 1.1x10(-8)), which appears to be independent of previously reported associations in this locus. Finally, we found no strong evidence for gene-environment interactions with smoking or hay fever status.

CONCLUSIONS

Population-based cohorts with simple asthma phenotypes represent a valuable and largely untapped resource for genetic studies of asthma.

Estimation of genetic effects on BMI during adolescence in an ethnically diverse cohort: The National Longitudinal Study of Adolescent Health

OBJECTIVE

The contribution of genetic variants to body mass index (BMI) during adolescence across multiethnic samples is largely unknown. We selected genetic loci associated with BMI or obesity in European-descent samples and examined them in a multiethnic adolescent sample.

DESIGN AND SAMPLE

In 5103 European American (EA), 1748 African American (AfA), 1304 Hispanic American (HA) and 439 Asian American (AsA) participants of the National Longitudinal Study of Adolescent Health (Add Health; ages 12-21 years, 47.5% male), we assessed the association between 41 established obesity-related single-nucleotide polymorphisms (SNPs) with BMI using additive genetic models, stratified by race/ethnicity, and in a pooled meta-analysis sample. We also compared the magnitude of effect for BMI-SNP associations in EA and AfA adolescents to comparable effect estimates from 11 861 EA and AfA adults in the Atherosclerosis Risk in Communities study (ages 45-64 years, 43.2% male).

RESULTS

Thirty-five of 41 BMI-SNP associations were directionally consistent with published studies in European populations, 18 achieved nominal significance (P<0.05; effect sizes from 0.19 to 0.71 kg m(-2) increase in BMI per effect allele), while 4 (FTO, TMEM18, TFAP2B, MC4R) remained significant after Bonferroni correction (P<0.0015). Of 41 BMI-SNP associations in AfA, HA and AsA adolescents, nine, three and five, respectively, were directionally consistent and nominally significant. In the pooled meta-analysis, 36 of 41 effect estimates were directionally consistent and 21 of 36 were nominally significant. In EA adolescents, BMI effect estimates were larger (P<0.05) for variants near TMEM18, PTER and MC4R and smaller for variants near MTIF3 and NRXN3 compared with EA adults.

CONCLUSION

Our findings suggest that obesity susceptibility loci may have a comparatively stronger role during adolescence than during adulthood, with variation across race/ethnic subpopulation.

Salt sensitivity of blood pressure is associated with polymorphisms in the sodium-bicarbonate cotransporter

Previous studies have demonstrated that single nucleotide polymorphisms (SNPs) of the sodium-bicarbonate co-transporter gene (SLC4A5) are associated with hypertension. We tested the hypothesis that SNPs in SLC4A5 are associated with salt sensitivity of blood pressure in 185 whites consuming an isocaloric constant diet with a randomized order of 7 days of low Na(+) (10 mmol/d) and 7 days of high Na(+) (300 mmol/d) intake. Salt sensitivity was defined as a ≥ 7-mm Hg increase in mean arterial pressure during a randomized transition between high and low Na(+) diet. A total of 35 polymorphisms in 17 candidate genes were assayed, 25 of which were tested for association. Association analyses with salt sensitivity revealed 3 variants that associated with salt sensitivity, 2 in SLC4A5 (P<0.001) and 1 in GRK4 (P=0.020). Of these, 2 SNPs in SLC4A5 (rs7571842 and rs10177833) demonstrated highly significant results and large effects sizes, using logistic regression. These 2 SNPs had P values of 1.0 × 10(-4) and 3.1 × 10(-4) with odds ratios of 0.221 and 0.221 in unadjusted regression models, respectively, with the G allele at both sites conferring protection. These SNPs remained significant after adjusting for body mass index and age (P=8.9 × 10(-5) and 2.6 × 10(-4) and odds ratios 0.210 and 0.286, respectively). Furthermore, the association of these SNPs with salt sensitivity was replicated in a second hypertensive population. Meta-analysis demonstrated significant associations of both SNPs with salt sensitivity (rs7571842 [P=1.2 × 10(-5)]; rs1017783 [P=1.1 × 10(-4)]). In conclusion, SLC4A5 variants are strongly associated with salt sensitivity of blood pressure in 2 separate white populations.

Variation in genes related to cochlear biology is strongly associated with adult-onset deafness in border collies

Domestic dogs can suffer from hearing losses that can have profound impacts on working ability and quality of life. We have identified a type of adult-onset hearing loss in Border Collies that appears to have a genetic cause, with an earlier age of onset (3–5 years) than typically expected for aging dogs (8–10 years). Studying this complex trait within pure breeds of dog may greatly increase our ability to identify genomic regions associated with risk of hearing impairment in dogs and in humans. We performed a genome-wide association study (GWAS) to detect loci underlying adult-onset deafness in a sample of 20 affected and 28 control Border Collies. We identified a region on canine chromosome 6 that demonstrates extended support for association surrounding SNP Chr6.25819273 (p-value = 1.09×10⁻¹³). To further localize disease-associated variants, targeted next-generation sequencing (NGS) of one affected and two unaffected dogs was performed. Through additional validation based on targeted genotyping of additional cases (n = 23 total) and controls (n = 101 total) and an independent replication cohort of 16 cases and 265 controls, we identified variants in USP31 that were strongly associated with adult-onset deafness in Border Collies, suggesting the involvement of the NF-κB pathway. We found additional support for involvement of RBBP6, which is critical for cochlear development. These findings highlight the utility of GWAS–guided fine-mapping of genetic loci using targeted NGS to study hereditary disorders of the domestic dog that may be analogous to human disorders.

The domestic dog offers a unique opportunity to study complex disorders similar to those seen in humans, but within the context of the much simpler genetic backgrounds of pure breeds, which represent closed populations. We performed a whole-genome search for genetic risk factors of adult-onset deafness in the Border Collie, a breed of herding dog that relies on acute hearing to perceive and respond to commands while working. Adult-onset deafness in Border Collies typically begins in early adulthood and is similar to age-related hearing loss in humans. This earlier onset has particular impact on the utility of working Border Collies and the livelihoods of their owners, and it appears to have a genetic cause. We identified three genetic variants that were strongly associated with adult-onset deafness in a sample of 405 Border Collies. These variants are located in two genes that have previously been linked to deafness, one involved in ear development and another that appears to mitigate tissue damage in the ear. These results provide new insight regarding genetic risk factors for age-related hearing loss in both dogs and humans.

Quick, "imputation-free" meta-analysis with proxy-SNPs

Background

Meta-analysis (MA) is widely used to pool genome-wide association studies (GWASes) in order to a) increase the power to detect strong or weak genotype effects or b) as a result verification method. As a consequence of differing SNP panels among genotyping chips, imputation is the method of choice within GWAS consortia to avoid losing too many SNPs in a MA. YAMAS (Yet Another Meta Analysis Software), however, enables cross-GWAS conclusions prior to finished and polished imputation runs, which eventually are time-consuming.

Results

Here we present a fast method to avoid forfeiting SNPs present in only a subset of studies, without relying on imputation. This is accomplished by using reference linkage disequilibrium data from 1,000 Genomes/HapMap projects to find proxy-SNPs together with in-phase alleles for SNPs missing in at least one study. MA is conducted by combining association effect estimates of a SNP and those of its proxy-SNPs. Our algorithm is implemented in the MA software YAMAS. Association results from GWAS analysis applications can be used as input files for MA, tremendously speeding up MA compared to the conventional imputation approach. We show that our proxy algorithm is well-powered and yields valuable ad hoc results, possibly providing an incentive for follow-up studies. We propose our method as a quick screening step prior to imputation-based MA, as well as an additional main approach for studies without available reference data matching the ethnicities of study participants. As a proof of principle, we analyzed six dbGaP Type II Diabetes GWAS and found that the proxy algorithm clearly outperforms naïve MA on the p-value level: for 17 out of 23 we observe an improvement on the p-value level by a factor of more than two, and a maximum improvement by a factor of 2127.

Conclusions

YAMAS is an efficient and fast meta-analysis program which offers various methods, including conventional MA as well as inserting proxy-SNPs for missing markers to avoid unnecessary power loss. MA with YAMAS can be readily conducted as YAMAS provides a generic parser for heterogeneous tabulated file formats within the GWAS field and avoids cumbersome setups. In this way, it supplements the meta-analysis process.

Integration of genome-wide association studies with biological knowledge identifies six novel genes related to kidney function

In conducting genome-wide association studies (GWAS), analytical approaches leveraging biological information may further understanding of the pathophysiology of clinical traits. To discover novel associations with estimated glomerular filtration rate (eGFR), a measure of kidney function, we developed a strategy for integrating prior biological knowledge into the existing GWAS data for eGFR from the CKDGen Consortium. Our strategy focuses on single nucleotide polymorphism (SNPs) in genes that are connected by functional evidence, determined by literature mining and gene ontology (GO) hierarchies, to genes near previously validated eGFR associations. It then requires association thresholds consistent with multiple testing, and finally evaluates novel candidates by independent replication. Among the samples of European ancestry, we identified a genome-wide significant SNP in FBXL20 (P = 5.6 × 10(-9)) in meta-analysis of all available data, and additional SNPs at the INHBC, LRP2, PLEKHA1, SLC3A2 and SLC7A6 genes meeting multiple-testing corrected significance for replication and overall P-values of 4.5 × 10(-4)-2.2 × 10(-7). Neither the novel PLEKHA1 nor FBXL20 associations, both further supported by association with eGFR among African Americans and with transcript abundance, would have been implicated by eGFR candidate gene approaches. LRP2, encoding the megalin receptor, was identified through connection with the previously known eGFR gene DAB2 and extends understanding of the megalin system in kidney function. These findings highlight integration of existing genome-wide association data with independent biological knowledge to uncover novel candidate eGFR associations, including candidates lacking known connections to kidney-specific pathways. The strategy may also be applicable to other clinical phenotypes, although more testing will be needed to assess its potential for discovery in general.

Genome-wide association and population genetic analysis of C-reactive protein in African American and Hispanic American women

C-reactive protein (CRP) is a systemic inflammation marker that predicts future cardiovascular risk. CRP levels are higher in African Americans and Hispanic Americans than in European Americans, but the genetic determinants of CRP in these admixed United States minority populations are largely unknown. We performed genome-wide association studies (GWASs) of 8,280 African American (AA) and 3,548 Hispanic American (HA) postmenopausal women from the Women's Health Initiative SNP Health Association Resource. We discovered and validated a CRP-associated variant of triggering receptors expressed by myeloid cells 2 (TREM2) in chromosomal region 6p21 (p = 10(-10)). The TREM2 variant associated with higher CRP is common in Africa but rare in other ancestral populations. In AA women, the CRP region in 1q23 contained a strong admixture association signal (p = 10(-17)), which appears to be related to several independent CRP-associated alleles; the strongest of these is present only in African ancestral populations and is associated with higher CRP. Of the other genomic loci previously associated with CRP through GWASs of European populations, most loci (LEPR, IL1RN, IL6R, GCKR, NLRP3, HNF1A, HNF4A, and APOC1) showed consistent patterns of association with CRP in AA and HA women. In summary, we have identified a common TREM2 variant associated with CRP in United States minority populations. The genetic architecture underlying the CRP phenotype in AA women is complex and involves genetic variants shared across populations, as well as variants specific to populations of African descent.

Genome-wide association study of intracranial aneurysms confirms role of Anril and SOX17 in disease risk

BACKGROUND

Genomewide association studies have identified novel genetic factors that contribute to intracranial aneurysm (IA) susceptibility. We sought to confirm previously reported loci, to identify novel risk factors, and to evaluate the contribution of these factors to familial and sporadic IA.

METHOD

We utilized 2 complementary samples, one recruited on the basis of a dense family history of IA (discovery sample 1: 388 IA cases and 397 controls) and the other without regard to family history (discovery sample 2: 1095 IA cases and 1286 controls). Imputation was used to generate a common set of single nucleotide polymorphisms (SNP) across samples, and a logistic regression model was used to test for association in each sample. Results from each sample were then combined in a metaanalysis.

RESULTS

There was only modest overlap in the association results obtained in the 2 samples. In neither sample did results reach genomewide significance. However, the metaanalysis yielded genomewide significance for SNP on chromosome 9p (CDKN2BAS; rs6475606; P=3.6×10(-8)) and provided further evidence to support the previously reported association of IA with SNP in SOX17 on chromosome 8q (rs1072737; P=8.7×10(-5)). Analyses suggest that the effect of smoking acts multiplicatively with the SNP genotype, and smoking has a greater effect on risk than SNP genotype.

CONCLUSIONS

In addition to replicating several previously reported loci, we provide further evidence that the association on chromosome 9p is attributable to variants in CDKN2BAS (also known as ANRIL, an antisense noncoding RNA).

Common variation at 6q16 within HACE1 and LIN28B influences susceptibility to neuroblastoma

Neuroblastoma is a cancer of the sympathetic nervous system that accounts for approximately 10% of all pediatric oncology deaths. Here, we report a genome-wide association study of 2,817 neuroblastoma cases and 7,473 controls. We identified two new associations at 6q16, the first within HACE1 (rs4336470; combined P=2.7×10(-11); odds ratio 1.26, 95% confidence interval (CI) 1.18-1.35) and the second within LIN28B (rs17065417; combined P=1.2×10(-8); odds ratio 1.38, 95% CI 1.23-1.54). Expression of LIN28B and let-7 miRNA correlated with rs17065417 genotype in neuroblastoma cell lines, and we observed significant growth inhibition upon depletion of LIN28B, specifically in neuroblastoma cells that were homozygous for the risk allele. Low HACE1 and high LIN28B expression in diagnostic primary neuroblastomas were associated with worse overall survival (P=0.008 and 0.014, respectively). Taken together, these data show that common variants in HACE1 and LIN28B influence neuroblastoma susceptibility and indicate that both genes likely have a role in disease progression.

The role of ECE1 variants in cognitive ability in old age and Alzheimer's disease risk

The β-amyloid peptide may play a central role in Alzheimer's disease (AD) pathogenesis. We have evaluated variants in seven Aβ-degrading genes (ACE, ECE1, ECE2, IDE, MME, PLAU, and TF) for association with AD risk in the Genetic and Environmental Risk in Alzheimer's Disease Consortium 1 (GERAD1) cohort, and with three cognitive phenotypes in the Lothian Birth Cohort 1936 (LBC1936), using 128 and 121 SNPs, respectively. In GERAD1, we identified a significant association between a four-SNP intragenic ECE1 haplotype and risk of AD in individuals that carried at least one APOE ε4 allele (P = 0.00035, odds ratio = 1.61). In LBC1936, we identified a significant association between a different two-SNP ECE1 intragenic haplotype and non-verbal reasoning in individuals lacking the APOE ε4 allele (P = 0.00036, β = -0.19). Both results showed a trend towards significance after permutation (0.05 < P < 0.10). A follow-up cognitive genetic study evaluated the association of ECE1 SNPs in three additional cohorts of non-demented older people. Meta-analysis of the four cohorts identified the significant association (Z < 0.05) of SNPs in the ECE-1b promoter with non-verbal reasoning scores, particularly in individuals lacking the APOE ε4 allele. Our genetic findings are not wholly consistent. Nonetheless, the AD associated intronic haplotype is linked to the 338A variant of known ECE1b promoter variant, 338C>A (rs213045). We observed significantly less expression from the 338A variant in two human neuroblastoma cell lines and speculate that this promoter may be subject to tissue-specific regulation.

Large-scale association analyses identify new loci influencing glycemic traits and provide insight into the underlying biological pathways

Through genome-wide association meta-analyses of up to 133,010 individuals of European ancestry without diabetes, including individuals newly genotyped using the Metabochip, we have increased the number of confirmed loci influencing glycemic traits to 53, of which 33 also increase type 2 diabetes risk (q < 0.05). Loci influencing fasting insulin concentration showed association with lipid levels and fat distribution, suggesting impact on insulin resistance. Gene-based analyses identified further biologically plausible loci, suggesting that additional loci beyond those reaching genome-wide significance are likely to represent real associations. This conclusion is supported by an excess of directionally consistent and nominally significant signals between discovery and follow-up studies. Functional analysis of these newly discovered loci will further improve our understanding of glycemic control.

Genome-wide association uncovers shared genetic effects among personality traits and mood states

Measures of personality and psychological distress are correlated and exhibit genetic covariance. We conducted univariate genome-wide SNP (~2.5 million) and gene-based association analyses of these traits and examined the overlap in results across traits, including a prediction analysis of mood states using genetic polygenic scores for personality. Measures of neuroticism, extraversion, and symptoms of anxiety, depression, and general psychological distress were collected in eight European cohorts (n ranged 546-1,338; maximum total n = 6,268) whose mean age ranged from 55 to 79 years. Meta-analysis of the cohort results was performed, with follow-up associations of the top SNPs and genes investigated in independent cohorts (n = 527-6,032). Suggestive association (P = 8 × 10(-8)) of rs1079196 in the FHIT gene was observed with symptoms of anxiety. Other notable associations (P < 6.09 × 10(-6)) included SNPs in five genes for neuroticism (LCE3C, POLR3A, LMAN1L, ULK3, SCAMP2), KIAA0802 for extraversion, and NOS1 for general psychological distress. An association between symptoms of depression and rs7582472 (near to MGAT5 and NCKAP5) was replicated in two independent samples, but other replication findings were less consistent. Gene-based tests identified a significant locus on chromosome 15 (spanning five genes) associated with neuroticism which replicated (P < 0.05) in an independent cohort. Support for common genetic effects among personality and mood (particularly neuroticism and depressive symptoms) was found in terms of SNP association overlap and polygenic score prediction. The variance explained by individual SNPs was very small (up to 1%) confirming that there are no moderate/large effects of common SNPs on personality and related traits.

New susceptibility loci associated with kidney disease in type 1 diabetes

Diabetic kidney disease, or diabetic nephropathy (DN), is a major complication of diabetes and the leading cause of end-stage renal disease (ESRD) that requires dialysis treatment or kidney transplantation. In addition to the decrease in the quality of life, DN accounts for a large proportion of the excess mortality associated with type 1 diabetes (T1D). Whereas the degree of glycemia plays a pivotal role in DN, a subset of individuals with poorly controlled T1D do not develop DN. Furthermore, strong familial aggregation supports genetic susceptibility to DN. However, the genes and the molecular mechanisms behind the disease remain poorly understood, and current therapeutic strategies rarely result in reversal of DN. In the GEnetics of Nephropathy: an International Effort (GENIE) consortium, we have undertaken a meta-analysis of genome-wide association studies (GWAS) of T1D DN comprising ~2.4 million single nucleotide polymorphisms (SNPs) imputed in 6,691 individuals. After additional genotyping of 41 top ranked SNPs representing 24 independent signals in 5,873 individuals, combined meta-analysis revealed association of two SNPs with ESRD: rs7583877 in the AFF3 gene (P = 1.2 × 10(-8)) and an intergenic SNP on chromosome 15q26 between the genes RGMA and MCTP2, rs12437854 (P = 2.0 × 10(-9)). Functional data suggest that AFF3 influences renal tubule fibrosis via the transforming growth factor-beta (TGF-β1) pathway. The strongest association with DN as a primary phenotype was seen for an intronic SNP in the ERBB4 gene (rs7588550, P = 2.1 × 10(-7)), a gene with type 2 diabetes DN differential expression and in the same intron as a variant with cis-eQTL expression of ERBB4. All these detected associations represent new signals in the pathogenesis of DN.

Genome-wide association study of African and European Americans implicates multiple shared and ethnic specific loci in sarcoidosis susceptibility

Sarcoidosis is a systemic inflammatory disease characterized by the formation of granulomas in affected organs. Genome-wide association studies (GWASs) of this disease have been conducted only in European population. We present the first sarcoidosis GWAS in African Americans (AAs, 818 cases and 1,088 related controls) followed by replication in independent sets of AAs (455 cases and 557 controls) and European Americans (EAs, 442 cases and 2,284 controls). We evaluated >6 million SNPs either genotyped using the Illumina Omni1-Quad array or imputed from the 1000 Genomes Project data. We identified a novel sarcoidosis-associated locus, NOTCH4, that reached genome-wide significance in the combined AA samples (rs715299, P(AA-meta) = 6.51 × 10(-10)) and demonstrated the independence of this locus from others in the MHC region in the same sample. We replicated previous European GWAS associations within HLA-DRA, HLA-DRB5, HLA-DRB1, BTNL2, and ANXA11 in both our AA and EA datasets. We also confirmed significant associations to the previously reported HLA-C and HLA-B regions in the EA but not AA samples. We further identified suggestive associations with several other genes previously reported in lung or inflammatory diseases.

Evaluating genetic risk for prostate cancer among Japanese and Latinos

BACKGROUND

There have been few genome-wide association studies (GWAS) of prostate cancer among diverse populations. To search for novel prostate cancer risk variants, we conducted GWAS of prostate cancer in Japanese and Latinos. In addition, we tested prostate cancer risk variants and developed genetic risk models of prostate cancer for Japanese and Latinos.

METHODS

Our first-stage GWAS of prostate cancer included Japanese (cases/controls = 1,033/1,042) and Latino (cases/controls = 1,043/1,057) from the Multiethnic Cohort (MEC). Significant associations from stage I (P < 1.0 × 10(-4)) were examined in silico in GWAS of prostate cancer (stage II) in Japanese (cases/controls = 1,583/3,386) and Europeans (cases/controls = 1,854/1,894).

RESULTS

No novel stage I single-nucleotide polymorphism (SNP) outside of known risk regions reached genome-wide significance. For Japanese, in stage I, the most notable putative novel association was seen with 10 SNPs (P ≤ 8.0 × 10(-6)) at chromosome 2q33; however, this was not replicated in stage II. For Latinos, the most significant association was observed with rs17023900 at the known 3p12 risk locus (stage I: OR = 1.45; P = 7.01 × 10(-5) and stage II: OR = 1.58; P = 3.05 × 10(-7)). The majority of the established risk variants for prostate cancer, 79% and 88%, were positively associated with prostate cancer in Japanese and Latinos (stage I), respectively. The cumulative effects of these variants significantly influence prostate cancer risk (OR per allele = 1.10; P = 2.71 × 10(-25) and OR = 1.07; P = 1.02 × 10(-16) for Japanese and Latinos, respectively).

CONCLUSION AND IMPACT

Our GWAS of prostate cancer did not identify novel genome-wide significant variants. However, our findings show that established risk variants for prostate cancer significantly contribute to risk among Japanese and Latinos.

Novel association approach for variable number tandem repeats (VNTRs) identifies DOCK5 as a susceptibility gene for severe obesity

Variable number tandem repeats (VNTRs) constitute a relatively under-examined class of genomic variants in the context of complex disease because of their sequence complexity and the challenges in assaying them. Recent large-scale genome-wide copy number variant mapping and association efforts have highlighted the need for improved methodology for association studies using these complex polymorphisms. Here we describe the in-depth investigation of a complex region on chromosome 8p21.2 encompassing the dedicator of cytokinesis 5 (DOCK5) gene. The region includes two VNTRs of complex sequence composition which flank a common 3975 bp deletion, all three of which were genotyped by polymerase chain reaction and fragment analysis in a total of 2744 subjects. We have developed a novel VNTR association method named VNTRtest, suitable for association analysis of multi-allelic loci with binary and quantitative outcomes, and have used this approach to show significant association of the DOCK5 VNTRs with childhood and adult severe obesity (P(empirical)= 8.9 × 10(-8) and P= 3.1 × 10(-3), respectively) which we estimate explains ~0.8% of the phenotypic variance. We also identified an independent association between the 3975 base pair (bp) deletion and obesity, explaining a further 0.46% of the variance (P(combined)= 1.6 × 10(-3)). Evidence for association between DOCK5 transcript levels and the 3975 bp deletion (P= 0.027) and both VNTRs (P(empirical)= 0.015) was also identified in adipose tissue from a Swedish family sample, providing support for a functional effect of the DOCK5 deletion and VNTRs. These findings highlight the potential role of DOCK5 in human obesity and illustrate a novel approach for analysis of the contribution of VNTRs to disease susceptibility through association studies.

Fine-mapping and initial characterization of QT interval loci in African Americans

The QT interval (QT) is heritable and its prolongation is a risk factor for ventricular tachyarrhythmias and sudden death. Most genetic studies of QT have examined European ancestral populations; however, the increased genetic diversity in African Americans provides opportunities to narrow association signals and identify population-specific variants. We therefore evaluated 6,670 SNPs spanning eleven previously identified QT loci in 8,644 African American participants from two Population Architecture using Genomics and Epidemiology (PAGE) studies: the Atherosclerosis Risk in Communities study and Women's Health Initiative Clinical Trial. Of the fifteen known independent QT variants at the eleven previously identified loci, six were significantly associated with QT in African American populations (P≤1.20×10(-4)): ATP1B1, PLN1, KCNQ1, NDRG4, and two NOS1AP independent signals. We also identified three population-specific signals significantly associated with QT in African Americans (P≤1.37×10(-5)): one at NOS1AP and two at ATP1B1. Linkage disequilibrium (LD) patterns in African Americans assisted in narrowing the region likely to contain the functional variants for several loci. For example, African American LD patterns showed that 0 SNPs were in LD with NOS1AP signal rs12143842, compared with European LD patterns that indicated 87 SNPs, which spanned 114.2 Kb, were in LD with rs12143842. Finally, bioinformatic-based characterization of the nine African American signals pointed to functional candidates located exclusively within non-coding regions, including predicted binding sites for transcription factors such as TBX5, which has been implicated in cardiac structure and conductance. In this detailed evaluation of QT loci, we identified several African Americans SNPs that better define the association with QT and successfully narrowed intervals surrounding established loci. These results demonstrate that the same loci influence variation in QT across multiple populations, that novel signals exist in African Americans, and that the SNPs identified as strong candidates for functional evaluation implicate gene regulatory dysfunction in QT prolongation.

Linkage analysis for plasma amyloid beta levels in persons with hypertension implicates Aβ-40 levels to presenilin 2

Plasma concentrations of Aβ40 and Aβ42 rise with age and are increased in people with mutations that cause early-onset Alzheimer's disease (AD). Amyloid beta (Aβ) plasma levels were successfully used as an (endo)phenotype for gene discovery using a linkage approach in families with dominant forms of disease. Here, we searched for loci involved in Aβ plasma levels in a series of non-demented patients with hypertension in the Erasmus Rucphen Family study. Aβ40 and Aβ42 levels were determined in 125 subjects with severe hypertension. All patients were genotyped with a 6,000 single nucleotide polymorphisms (SNPs) illumina array designed for linkage analysis. We conducted linkage analysis of plasma Aβ levels. None of the linkage analyses yielded genome-wide significant logarithm of odds (LOD) score over 3.3, but there was suggestive evidence for linkage (LOD > 1.9) for two regions: 1q41 (LOD = 2.07) and 11q14.3 (LOD = 2.97), both for Aβ40. These regions were followed up with association analysis in the study subjects and in 320 subjects from a population-based cohort. For the Aβ40 region on chromosome 1, association of several SNPs was observed at the presenilin 2 gene (PSEN2) (p = 2.58 × 10(-4) for rs6703170). On chromosome 11q14-21, we found some association (p = 3.1 × 10(-3) for rs2514299). This linkage study of plasma concentrations of Aβ40 and Aβ42 yielded two suggestive regions, of which one points toward a known locus for familial AD.

Insulin clearance: confirmation as a highly heritable trait, and genome-wide linkage analysis

AIMS/HYPOTHESIS

We have previously documented a high heritability of insulin clearance in a Hispanic cohort. Here, our goal was to confirm the high heritability in a second cohort and search for genetic loci contributing to insulin clearance.

METHODS

Hyperinsulinaemic-euglycaemic clamps were performed in 513 participants from 140 Hispanic families. Heritability was estimated for clamp-derived insulin clearance and a two-phase genome-wide linkage scan was conducted using a variance components approach. Linkage peaks were further investigated by candidate gene association analysis in two cohorts.

RESULTS

The covariate-adjusted heritability of insulin clearance was 73%, indicating that the majority of the phenotypic variance is due to genetic factors. In the Phase 1 linkage scan, no signals with a logarithm of odds (LOD) score >2 were detected. In the Phase 2 scan, two linkage peaks with an LOD >2 for insulin clearance were identified on chromosomes 15 (LOD 3.62) and 20 (LOD 2.43). These loci harbour several promising candidate genes for insulin clearance, with 12 single nucleotide polymorphisms (SNPs) on chromosome 15 and six SNPs on chromosome 20 being associated with insulin clearance in both Hispanic cohorts.

CONCLUSIONS/INTERPRETATION

In a second Hispanic cohort, we confirmed that insulin clearance is a highly heritable trait and identified chromosomal loci that harbour genes regulating insulin clearance. The identification of such genes may improve our understanding of how the body clears insulin, thus leading to improved risk assessment, diagnosis, prevention and therapy of diabetes, as well as of other hyperinsulinaemic disorders, such as the metabolic syndrome and polycystic ovary syndrome.

Extended analysis of a genome-wide association study in primary sclerosing cholangitis detects multiple novel risk loci

BACKGROUND & AIMS

A limited number of genetic risk factors have been reported in primary sclerosing cholangitis (PSC). To discover further genetic susceptibility factors for PSC, we followed up on a second tier of single nucleotide polymorphisms (SNPs) from a genome-wide association study (GWAS).

METHODS

We analyzed 45 SNPs in 1221 PSC cases and 3508 controls. The association results from the replication analysis and the original GWAS (715 PSC cases and 2962 controls) were combined in a meta-analysis comprising 1936 PSC cases and 6470 controls. We performed an analysis of bile microbial community composition in 39 PSC patients by 16S rRNA sequencing.

RESULTS

Seventeen SNPs representing 12 distinct genetic loci achieved nominal significance (p(replication) <0.05) in the replication. The most robust novel association was detected at chromosome 1p36 (rs3748816; p(combined)=2.1 × 10(-8)) where the MMEL1 and TNFRSF14 genes represent potential disease genes. Eight additional novel loci showed suggestive evidence of association (p(repl) <0.05). FUT2 at chromosome 19q13 (rs602662; p(comb)=1.9 × 10(-6), rs281377; p(comb)=2.1 × 10(-6) and rs601338; p(comb)=2.7 × 10(-6)) is notable due to its implication in altered susceptibility to infectious agents. We found that FUT2 secretor status and genotype defined by rs601338 significantly influence biliary microbial community composition in PSC patients.

CONCLUSIONS

We identify multiple new PSC risk loci by extended analysis of a PSC GWAS. FUT2 genotype needs to be taken into account when assessing the influence of microbiota on biliary pathology in PSC.

Impact of common variation in bone-related genes on type 2 diabetes and related traits

Exploring genetic pleiotropy can provide clues to a mechanism underlying the observed epidemiological association between type 2 diabetes and heightened fracture risk. We examined genetic variants associated with bone mineral density (BMD) for association with type 2 diabetes and glycemic traits in large well-phenotyped and -genotyped consortia. We undertook follow-up analysis in ∼19,000 individuals and assessed gene expression. We queried single nucleotide polymorphisms (SNPs) associated with BMD at levels of genome-wide significance, variants in linkage disequilibrium (r(2) > 0.5), and BMD candidate genes. SNP rs6867040, at the ITGA1 locus, was associated with a 0.0166 mmol/L (0.004) increase in fasting glucose per C allele in the combined analysis. Genetic variants in the ITGA1 locus were associated with its expression in the liver but not in adipose tissue. ITGA1 variants appeared among the top loci associated with type 2 diabetes, fasting insulin, β-cell function by homeostasis model assessment, and 2-h post-oral glucose tolerance test glucose and insulin levels. ITGA1 has demonstrated genetic pleiotropy in prior studies, and its suggested role in liver fibrosis, insulin secretion, and bone healing lends credence to its contribution to both osteoporosis and type 2 diabetes. These findings further underscore the link between skeletal and glucose metabolism and highlight a locus to direct future investigations.

Association testing of previously reported variants in a large case-control meta-analysis of diabetic nephropathy

We formed the GEnetics of Nephropathy-an International Effort (GENIE) consortium to examine previously reported genetic associations with diabetic nephropathy (DN) in type 1 diabetes. GENIE consists of 6,366 similarly ascertained participants of European ancestry with type 1 diabetes, with and without DN, from the All Ireland-Warren 3-Genetics of Kidneys in Diabetes U.K. and Republic of Ireland (U.K.-R.O.I.) collection and the Finnish Diabetic Nephropathy Study (FinnDiane), combined with reanalyzed data from the Genetics of Kidneys in Diabetes U.S. Study (U.S. GoKinD). We found little evidence for the association of the EPO promoter polymorphism, rs161740, with the combined phenotype of proliferative retinopathy and end-stage renal disease in U.K.-R.O.I. (odds ratio [OR] 1.14, P = 0.19) or FinnDiane (OR 1.06, P = 0.60). However, a fixed-effects meta-analysis that included the previously reported cohorts retained a genome-wide significant association with that phenotype (OR 1.31, P = 2 × 10(-9)). An expanded investigation of the ELMO1 locus and genetic regions reported to be associated with DN in the U.S. GoKinD yielded only nominal statistical significance for these loci. Finally, top candidates identified in a recent meta-analysis failed to reach genome-wide significance. In conclusion, we were unable to replicate most of the previously reported genetic associations for DN, and significance for the EPO promoter association was attenuated.

Functional and genetic analysis of haplotypic sequence variation at the nicastrin genomic locus

Nicastrin (NCSTN) is a component of the γ-secretase complex and therefore potentially a candidate risk gene for Alzheimer's disease. Here, we have developed a novel functional genomics methodology to express common locus haplotypes to assess functional differences. DNA recombination was used to engineer 5 bacterial artificial chromosomes (BACs) to each express a different haplotype of the NCSTN locus. Each NCSTN-BAC was delivered to knockout nicastrin (Ncstn(-/-)) cells and clonal NCSTN-BAC(+)/Ncstn(-/-) cell lines were created for functional analyses. We showed that all NCSTN-BAC haplotypes expressed nicastrin protein and rescued γ-secretase activity and amyloid beta (Aβ) production in NCSTN-BAC(+)/Ncstn(-/-) lines. We then showed that genetic variation at the NCSTN locus affected alternative splicing in human postmortem brain tissue. However, there was no robust functional difference between clonal cell lines rescued by each of the 5 different haplotypes. Finally, there was no statistically significant association of NCSTN with disease risk in the 4 cohorts. We therefore conclude that it is unlikely that common variation at the NCSTN locus is a risk factor for Alzheimer's disease.

Identification of a novel percent mammographic density locus at 12q24

Percent mammographic density adjusted for age and body mass index (BMI) is one of the strongest risk factors for breast cancer and has a heritable component that remains largely unidentified. We performed a three-stage genome-wide association study (GWAS) of percent mammographic density to identify novel genetic loci associated with this trait. In stage 1, we combined three GWASs of percent density comprised of 1241 women from studies at the Mayo Clinic and identified the top 48 loci (99 single nucleotide polymorphisms). We attempted replication of these loci in 7018 women from seven additional studies (stage 2). The meta-analysis of stage 1 and 2 data identified a novel locus, rs1265507 on 12q24, associated with percent density, adjusting for age and BMI (P = 4.43 × 10(-8)). We refined the 12q24 locus with 459 additional variants (stage 3) in a combined analysis of all three stages (n = 10 377) and confirmed that rs1265507 has the strongest association in the 12q24 region (P = 1.03 × 10(-8)). Rs1265507 is located between the genes TBX5 and TBX3, which are members of the phylogenetically conserved T-box gene family and encode transcription factors involved in developmental regulation. Understanding the mechanism underlying this association will provide insight into the genetics of breast tissue composition.

Genetic associations for activated partial thromboplastin time and prothrombin time, their gene expression profiles, and risk of coronary artery disease

Activated partial thromboplastin time (aPTT) and prothrombin time (PT) are clinical tests commonly used to screen for coagulation-factor deficiencies. One genome-wide association study (GWAS) has been reported previously for aPTT, but no GWAS has been reported for PT. We conducted a GWAS and meta-analysis to identify genetic loci for aPTT and PT. The GWAS for aPTT was conducted in 9,240 individuals of European ancestry from the Atherosclerosis Risk in Communities (ARIC) study, and the GWAS for PT was conducted in 2,583 participants from the Genetic Study of Three Population Microisolates in South Tyrol (MICROS) and the Lothian Birth Cohorts (LBC) of 1921 and 1936. Replication was assessed in 1,041 to 3,467 individuals. For aPTT, previously reported associations with KNG1, HRG, F11, F12, and ABO were confirmed. A second independent association in ABO was identified and replicated (rs8176704, p = 4.26 × 10(-24)). Pooling the ARIC and replication data yielded two additional loci in F5 (rs6028, p = 3.22 × 10(-9)) and AGBL1 (rs2469184, p = 3.61 × 10(-8)). For PT, significant associations were identified and confirmed in F7 (rs561241, p = 3.71 × 10(-56)) and PROCR/EDEM2 (rs2295888, p = 5.25 × 10(-13)). Assessment of existing gene expression and coronary artery disease (CAD) databases identified associations of five of the GWAS loci with altered gene expression and two with CAD. In summary, eight genetic loci that account for ∼29% of the variance in aPTT and two loci that account for ∼14% of the variance in PT were detected and supported by functional data.

A follow-up association study of two genetic variants for bone mineral density variation in Caucasians

SUMMARY

We tested whether two genetic variants were associated with BMD at multiple clinically relevant skeletal sites in Caucasians. We found that variant rs7776725 is consistently associated with hip, spine, wrist and whole-body BMD, which highlights the potential importance of this variant or linked variants for osteoporosis.

INTRODUCTION

A recent genome-wide association study identified two single nucleotide polymorphisms (SNPs), rs7776725 and rs1721400, that were associated with bone mineral density (BMD) variation at the radius, tibia and calcaneus in a Korean population. In this study, we aimed to test whether the association of these two genetic variants can be replicated in Caucasians and whether their association with BMD can be extended to other clinically relevant skeletal sites.

METHODS

We performed this study in two large cohorts of unrelated US Caucasians. Area BMD at the hip, spine, wrist (ultra-distal radius) and whole body were measured with Hologic dual-energy X-ray absorptiometer. SNPs were genotyped with Affymetrix human genome-wide genotyping arrays. Association analyses were performed using PLINK.

RESULTS

We detected highly significant association (combined p = 1.42 × 10(-16)) of rs7776725 with wrist BMD but only borderline association signal (combined p = 0.017) for rs1721400 with wrist BMD. In addition, we found that rs7776725 was associated with BMD at the hip, spine and whole body. At the FAM3C gene locus where rs7776725 was located, we identified several other SNPs (rs4727922, rs1803389, rs718766 and rs7793554) that were also associated with BMD.

CONCLUSIONS

This is the first follow-up association study of rs7776725 and rs1721400 with BMD. The rs7776725 showed consistent association with BMD at multiple clinically important skeletal sites, which highlighted the potential importance of rs7776725 or linked SNPs for risk of osteoporosis. Further in-depth re-sequencing studies and functional assays are necessary to elucidate the underlying mechanisms.

Genome-wide association study in Han Chinese identifies four new susceptibility loci for coronary artery disease

We performed a meta-analysis of 2 genome-wide association studies of coronary artery disease comprising 1,515 cases and 5,019 controls followed by replication studies in 15,460 cases and 11,472 controls, all of Chinese Han ancestry. We identify four new loci for coronary artery disease that reached the threshold of genome-wide significance (P < 5 × 10(-8)). These loci mapped in or near TTC32-WDR35, GUCY1A3, C6orf10-BTNL2 and ATP2B1. We also replicated four loci previously identified in European populations (in or near PHACTR1, TCF21, CDKN2A-CDKN2B and C12orf51). These findings provide new insights into pathways contributing to the susceptibility for coronary artery disease in the Chinese Han population.

C-reactive protein (CRP) promoter polymorphisms influence circulating CRP levels in a genome-wide association study of African Americans

C-reactive protein (CRP) is an acute phase reactant protein produced primarily by the liver. Circulating CRP levels are influenced by genetic and non-genetic factors, including infection and obesity. Genome-wide association studies (GWAS) provide an unbiased approach towards identifying loci influencing CRP levels. None of the six GWAS for CRP levels has been conducted in an African ancestry population. The present study aims to: (i) identify genetic variants that influence serum CRP in African Americans (AA) using a genome-wide association approach and replicate these findings in West Africans (WA), (ii) assess transferability of major signals for CRP reported in European ancestry populations (EA) to AA and (iii) use the weak linkage disequilibrium (LD) structure characteristic of African ancestry populations to fine-map the previously reported CRP locus. The discovery cohort comprised 837 unrelated AA, with the replication of significant single-nucleotide polymorphisms (SNPs) assessed in 486 WA. The association analysis was conducted with 2 366 856 genotyped and imputed SNPs under an additive genetic model with adjustment for appropriate covariates. Genome-wide and replication significances were set at P < 5 × 10(-8) and P < 0.05, respectively. Ten SNPs in (CRP pseudogene-1) CRPP1 and CRP genes were associated with serum CRP (P = 2.4 × 10(-09) to 4.3 × 10(-11)). All but one of the top-scoring SNPs associated with CRP in AA were successfully replicated in WA. CRP signals previously identified in EA samples were transferable to AAs, and we were able to fine-map this signal, reducing the region of interest from the 25 kb of LD around the locus in the HapMap CEU sample to only 8 kb in our AA sample.

A functional haplotype of UBE2L3 confers risk for systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is an autoimmune disease with diverse clinical manifestations characterized by the development of pathogenic autoantibodies manifesting in inflammation of target organs such as the kidneys, skin and joints. Genome-wide association studies have identified genetic variants in the UBE2L3 region that are associated with SLE in subjects of European and Asian ancestry. UBE2L3 encodes an ubiquitin-conjugating enzyme, UBCH7, involved in cell proliferation and immune function. In this study, we sought to further characterize the genetic association in the region of UBE2L3 and use molecular methods to determine the functional effect of the risk haplotype. We identified significant associations between variants in the region of UBE2L3 and SLE in individuals of European and Asian ancestry that exceeded a Bonferroni-corrected threshold (P<1 × 10(-4)). A single risk haplotype was observed in all associated populations. Individuals harboring the risk haplotype display a significant increase in both UBE2L3 mRNA expression (P=0.0004) and UBCH7 protein expression (P=0.0068). The results suggest that variants carried on the SLE-associated UBE2L3 risk haplotype influence autoimmunity by modulating UBCH7 expression.

Multiple loci influencing hippocampal degeneration identified by genome scan

OBJECTIVE

Large genome-wide association studies (GWASs) have identified many novel genes influencing Alzheimer disease (AD) risk, but most of the genetic variance remains unexplained. We conducted a 2-stage GWAS for AD-related quantitative measures of hippocampal volume (HV), total cerebral volume (TCV), and white matter hyperintensities (WMH).

METHODS

Brain magnetic resonance imaging measures of HV, TCV, and WMH were obtained from 981 Caucasian and 419 African American AD cases and their cognitively normal siblings in the MIRAGE (Multi Institutional Research in Alzheimer's Genetic Epidemiology) Study, and from 168 AD cases, 336 individuals with mild cognitive impairment, and 188 controls in the Alzheimer's Disease Neuroimaging Initiative Study. A GWAS for each trait was conducted in the 2 Caucasian data sets in stage 1. Results from the 2 data sets were combined by meta-analysis. In stage 2, 1 single nucleotide polymorphism (SNP) from each region that was nominally significant in each data set (p < 0.05) and strongly associated in both data sets (p < 1.0 × 10(-5)) was evaluated in the African American data set.

RESULTS

Twenty-two markers (14 for HV, 3 for TCV, and 5 for WMH) from distinct regions met criteria for evaluation in stage 2. Novel genome-wide significant associations (p < 5.0 × 10(-8)) were attained for HV with SNPs in the APOE, F5/SELP, LHFP, and GCFC2 gene regions. All of these associations were supported by evidence in each data set. Associations with different SNPs in the same gene (p < 1 × 10(-5) in Caucasians and p < 2.2 × 10(-4) in African Americans) were also observed for PICALM with HV, SYNPR with TCV, and TTC27 with WMH.

INTERPRETATION

Our study demonstrates the efficacy of endophenotypes for broadening our understanding of the genetic basis of AD.

Novel late-onset Alzheimer disease loci variants associate with brain gene expression

OBJECTIVE

Recent genome-wide association studies (GWAS) of late-onset Alzheimer disease (LOAD) identified 9 novel risk loci. Discovery of functional variants within genes at these loci is required to confirm their role in Alzheimer disease (AD). Single nucleotide polymorphisms that influence gene expression (eSNPs) constitute an important class of functional variants. We therefore investigated the influence of the novel LOAD risk loci on human brain gene expression.

METHODS

We measured gene expression levels in the cerebellum and temporal cortex of autopsied AD subjects and those with other brain pathologies (∼400 total subjects). To determine whether any of the novel LOAD risk variants are eSNPs, we tested their cis-association with expression of 6 nearby LOAD candidate genes detectable in human brain (ABCA7, BIN1, CLU, MS4A4A, MS4A6A, PICALM) and an additional 13 genes ±100 kb of these SNPs. To identify additional eSNPs that influence brain gene expression levels of the novel candidate LOAD genes, we identified SNPs ±100 kb of their location and tested for cis-associations.

RESULTS

CLU rs11136000 (p = 7.81 × 10(-4)) and MS4A4A rs2304933/rs2304935 (p = 1.48 × 10(-4)-1.86 × 10(-4)) significantly influence temporal cortex expression levels of these genes. The LOAD-protective CLU and risky MS4A4A locus alleles associate with higher brain levels of these genes. There are other cis-variants that significantly influence brain expression of CLU and ABCA7 (p = 4.01 × 10(-5)-9.09 × 10(-9)), some of which also associate with AD risk (p = 2.64 × 10(-2)-6.25 × 10(-5)).

CONCLUSIONS

CLU and MS4A4A eSNPs may at least partly explain the LOAD risk association at these loci. CLU and ABCA7 may harbor additional strong eSNPs. These results have implications in the search for functional variants at the novel LOAD risk loci.

Common genetic variants in the microRNA biogenesis pathway are not associated with breast cancer risk in Asian women

BACKGROUND

Although the role of miRNA in cancer development and progression has been well established, the association between genetic variants in miRNA biogenesis pathway genes and breast cancer risk has been yet unclear.

METHODS

We analyzed data from two genome-wide association studies conducted in East Asian women including 5,066 cases and 4,337 controls. Among the single-nucleotide polymorphisms (SNP), which were directly genotyped or imputed, we selected 237 SNPs in 32 genes involved in miRNA biogenesis pathway and its regulation.

RESULTS

Although eight SNPs were nominally associated with breast cancer risk in combined samples (P < 0.05), none of them were significant after adjustment for multiple comparisons.

CONCLUSIONS

The common genetic variants in miRNA biogenesis pathway genes may not be associated with breast cancer risk.

IMPACT

This study suggests no association between the polymorphisms in miRNA biogenesis pathway genes and breast cancer risk. Studies with large sample size and more genetic variants should be warranted to adequately evaluate the potential association.

Brain expression genome-wide association study (eGWAS) identifies human disease-associated variants

Genetic variants that modify brain gene expression may also influence risk for human diseases. We measured expression levels of 24,526 transcripts in brain samples from the cerebellum and temporal cortex of autopsied subjects with Alzheimer's disease (AD, cerebellar n = 197, temporal cortex n = 202) and with other brain pathologies (non–AD, cerebellar n = 177, temporal cortex n = 197). We conducted an expression genome-wide association study (eGWAS) using 213,528 cisSNPs within ±100 kb of the tested transcripts. We identified 2,980 cerebellar cisSNP/transcript level associations (2,596 unique cisSNPs) significant in both ADs and non–ADs (q<0.05, p = 7.70×10⁻⁵–1.67×10⁻⁸²). Of these, 2,089 were also significant in the temporal cortex (p = 1.85×10⁻⁵–1.70×10⁻¹⁴¹). The top cerebellar cisSNPs had 2.4-fold enrichment for human disease-associated variants (p<10⁻⁶). We identified novel cisSNP/transcript associations for human disease-associated variants, including progressive supranuclear palsy SLCO1A2/rs11568563, Parkinson's disease (PD) MMRN1/rs6532197, Paget's disease OPTN/rs1561570; and we confirmed others, including PD MAPT/rs242557, systemic lupus erythematosus and ulcerative colitis IRF5/rs4728142, and type 1 diabetes mellitus RPS26/rs1701704. In our eGWAS, there was 2.9–3.3 fold enrichment (p<10⁻⁶) of significant cisSNPs with suggestive AD–risk association (p<10⁻³) in the Alzheimer's Disease Genetics Consortium GWAS. These results demonstrate the significant contributions of genetic factors to human brain gene expression, which are reliably detected across different brain regions and pathologies. The significant enrichment of brain cisSNPs among disease-associated variants advocates gene expression changes as a mechanism for many central nervous system (CNS) and non–CNS diseases. Combined assessment of expression and disease GWAS may provide complementary information in discovery of human disease variants with functional implications. Our findings have implications for the design and interpretation of eGWAS in general and the use of brain expression quantitative trait loci in the study of human disease genetics.

Genetic variants that regulate gene expression levels can also influence human disease risk. Discovery of genomic loci that alter brain gene expression levels (brain expression quantitative trait loci = eQTLs) can be instrumental in the identification of genetic risk underlying both central nervous system (CNS) and non–CNS diseases. To systematically assess the role of brain eQTLs in human disease and to evaluate the influence of brain region and pathology in eQTL mapping, we performed an expression genome-wide association study (eGWAS) in 773 brain samples from the cerebellum and temporal cortex of ∼200 autopsied subjects with Alzheimer's disease (AD) and ∼200 with other brain pathologies (non–AD). We identified ∼3,000 significant associations between cisSNPs near ∼700 genes and their cerebellar transcript levels, which replicate in ADs and non–ADs. More than 2,000 of these associations were reproducible in the temporal cortex. The top cisSNPs are enriched for both CNS and non–CNS disease-associated variants. We identified novel and confirmed previous cisSNP/transcript associations for many disease loci, suggesting gene expression regulation as their mechanism of action. These findings demonstrate the reproducibility of the eQTL approach across different brain regions and pathologies, and advocate the combined use of gene expression and disease GWAS for identification and functional characterization of human disease-associated variants.

Genome-wide association study of circulating estradiol, testosterone, and sex hormone-binding globulin in postmenopausal women

Genome-wide association studies (GWAS) have successfully identified common genetic variants that contribute to breast cancer risk. Discovering additional variants has become difficult, as power to detect variants of weaker effect with present sample sizes is limited. An alternative approach is to look for variants associated with quantitative traits that in turn affect disease risk. As exposure to high circulating estradiol and testosterone, and low sex hormone-binding globulin (SHBG) levels is implicated in breast cancer etiology, we conducted GWAS analyses of plasma estradiol, testosterone, and SHBG to identify new susceptibility alleles. Cancer Genetic Markers of Susceptibility (CGEMS) data from the Nurses’ Health Study (NHS), and Sisters in Breast Cancer Screening data were used to carry out primary meta-analyses among ∼1600 postmenopausal women who were not taking postmenopausal hormones at blood draw. We observed a genome-wide significant association between SHBG levels and rs727428 (joint β = -0.126; joint P = 2.09×10^–16), downstream of the SHBG gene. No genome-wide significant associations were observed with estradiol or testosterone levels. Among variants that were suggestively associated with estradiol (P<10^–5), several were located at the CYP19A1 gene locus. Overall results were similar in secondary meta-analyses that included ∼900 NHS current postmenopausal hormone users. No variant associated with estradiol, testosterone, or SHBG at P<10^–5 was associated with postmenopausal breast cancer risk among CGEMS participants. Our results suggest that the small magnitude of difference in hormone levels associated with common genetic variants is likely insufficient to detectably contribute to breast cancer risk.

Caution in interpreting results from imputation analysis when linkage disequilibrium extends over a large distance: a case study on venous thrombosis

By applying an imputation strategy based on the 1000 Genomes project to two genome-wide association studies (GWAS), we detected a susceptibility locus for venous thrombosis on chromosome 11p11.2 that was missed by previous GWAS analyses that had been conducted on the same datasets. A comprehensive linkage disequilibrium and haplotype analysis of the whole locus where twelve SNPs exhibited association p-values lower than 2.23 10(-11) and the use of independent case-control samples demonstrated that the culprit variant was a rare variant located ~1 Mb away from the original hits, not tagged by current genome-wide genotyping arrays and even not well imputed in the original GWAS samples. This variant was in fact the rs1799963, also known as the FII G20210A prothrombin mutation. This work may be of major interest not only for its scientific impact but also for its methodological findings.