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

A coding variant in CR1 interacts with APOE-ε4 to influence cognitive decline

Complement receptor 1 (CR1) is an Alzheimer's disease (AD) susceptibility locus that also influences AD-related traits such as episodic memory decline and neuritic amyloid plaque deposition. We implemented a functional fine-mapping approach, leveraging intermediate phenotypes to identify functional variant(s) within the CR1 locus. Using 1709 subjects (697 deceased) from the Religious Orders Study and the Rush Memory and Aging Project, we tested 41 single-nucleotide polymorphisms (SNPs) within the linkage disequilibrium block containing the published CR1 AD SNP (rs6656401) for associations with episodic memory decline, and then examined the functional consequences of the top result. We report that a coding variant in the LHR-D (long homologous repeat D) region of the CR1 gene, rs4844609 (Ser1610Thr, minor allele frequency = 0.02), is associated with episodic memory decline and accounts for the known effect of the index SNP rs6656401 (D' = 1, r(2)= 0.084) on this trait. Further, we demonstrate that the coding variant's effect is largely dependent on an interaction with APOE-ε4 and mediated by an increased burden of AD-related neuropathology. Finally, in our data, this coding variant is also associated with AD susceptibility (joint odds ratio = 1.4). Taken together, our analyses identify a CR1 coding variant that influences episodic memory decline; it is a variant known to alter the conformation of CR1 and points to LHR-D as the functional domain within the CR1 protein that mediates the effect on memory decline. We thus implicate C1q and MBL, which bind to LHR-D, as likely targets of the variant's effect and suggest that CR1 may be an important intermediate in the clearance of Aβ42 particles by C1q.

Genome-wide association for abdominal subcutaneous and visceral adipose reveals a novel locus for visceral fat in women

Body fat distribution, particularly centralized obesity, is associated with metabolic risk above and beyond total adiposity. We performed genome-wide association of abdominal adipose depots quantified using computed tomography (CT) to uncover novel loci for body fat distribution among participants of European ancestry. Subcutaneous and visceral fat were quantified in 5,560 women and 4,997 men from 4 population-based studies. Genome-wide genotyping was performed using standard arrays and imputed to ~2.5 million Hapmap SNPs. Each study performed a genome-wide association analysis of subcutaneous adipose tissue (SAT), visceral adipose tissue (VAT), VAT adjusted for body mass index, and VAT/SAT ratio (a metric of the propensity to store fat viscerally as compared to subcutaneously) in the overall sample and in women and men separately. A weighted z-score meta-analysis was conducted. For the VAT/SAT ratio, our most significant p-value was rs11118316 at LYPLAL1 gene (p = 3.1 × 10E-09), previously identified in association with waist-hip ratio. For SAT, the most significant SNP was in the FTO gene (p = 5.9 × 10E-08). Given the known gender differences in body fat distribution, we performed sex-specific analyses. Our most significant finding was for VAT in women, rs1659258 near THNSL2 (p = 1.6 × 10-08), but not men (p = 0.75). Validation of this SNP in the GIANT consortium data demonstrated a similar sex-specific pattern, with observed significance in women (p = 0.006) but not men (p = 0.24) for BMI and waist circumference (p = 0.04 [women], p = 0.49 [men]). Finally, we interrogated our data for the 14 recently published loci for body fat distribution (measured by waist-hip ratio adjusted for BMI); associations were observed at 7 of these loci. In contrast, we observed associations at only 7/32 loci previously identified in association with BMI; the majority of overlap was observed with SAT. Genome-wide association for visceral and subcutaneous fat revealed a SNP for VAT in women. More refined phenotypes for body composition and fat distribution can detect new loci not previously uncovered in large-scale GWAS of anthropometric traits.

Genome-wide association of pericardial fat identifies a unique locus for ectopic fat

Pericardial fat is a localized fat depot associated with coronary artery calcium and myocardial infarction. We hypothesized that genetic loci would be associated with pericardial fat independent of other body fat depots. Pericardial fat was quantified in 5,487 individuals of European ancestry from the Framingham Heart Study (FHS) and the Multi-Ethnic Study of Atherosclerosis (MESA). Genotyping was performed using standard arrays and imputed to ∼2.5 million Hapmap SNPs. Each study performed a genome-wide association analysis of pericardial fat adjusted for age, sex, weight, and height. A weighted z-score meta-analysis was conducted, and validation was obtained in an additional 3,602 multi-ethnic individuals from the MESA study. We identified a genome-wide significant signal in our primary meta-analysis at rs10198628 near TRIB2 (MAF 0.49, p = 2.7×10^-08). This SNP was not associated with visceral fat (p = 0.17) or body mass index (p = 0.38), although we observed direction-consistent, nominal significance with visceral fat adjusted for BMI (p = 0.01) in the Framingham Heart Study. Our findings were robust among African ancestry (n = 1,442, p = 0.001), Hispanic (n = 1,399, p = 0.004), and Chinese (n = 761, p = 0.007) participants from the MESA study, with a combined p-value of 5.4E-14. We observed TRIB2 gene expression in the pericardial fat of mice. rs10198628 near TRIB2 is associated with pericardial fat but not measures of generalized or visceral adiposity, reinforcing the concept that there are unique genetic underpinnings to ectopic fat distribution.

Pericardial fat is a localized fat depot associated with coronary artery calcium and myocardial infarction. To test whether genetic loci are associated with pericardial fat independent of other body fat depots, we measured pericardial fat in 5,487 individuals of European ancestry. After performing an unbiased screen using genome-wide association, we identified a genome-wide significant signal in our primary meta-analysis at rs10198628 near TRIB2 (MAF 0.49, p = 2.7×10^-08). This SNP was not associated with visceral fat (p = 0.17) or body mass index (p = 0.38). Our findings were robust among multi-ethnic participants from the MESA study, with a combined p-value of 5.4E-14. We observed TRIB2 gene expression in the pericardial fat of mice. rs10198628 near TRIB2 is associated with pericardial fat but not measures of generalized or visceral adiposity, reinforcing the concept that there are unique genetic underpinnings to ectopic fat distribution.

A study of prostaglandin pathway genes and interactions with current nonsteroidal anti-inflammatory drug use in colorectal adenoma

Colorectal cancer (CRC) is the second leading cause of cancer-related death and usually arises from colorectal polyps. Screening and removal of polyps reduce mortality from CRC. Colorectal polyps are known to aggregate in families; however the genetic determinants for risk of polyps are unknown. In addition, it has been shown that nonsteroidal anti-inflammatory drug (NSAID) use decreases the risk of CRC and the incidence and size of polyps. In this study, we used data from the Tennessee Colorectal Polyp Study and the Tennessee-Indiana Adenoma Recurrence Study to evaluate selected genes from the prostaglandin (PG) metabolism and signaling pathways for association with risk of polyps and for interactions with NSAIDs. Our design consisted of discovery and replication phases for a total of 2,551 Caucasian polyp cases and 3,285 Caucasian controls. We carried out multivariable logistic regression to test for association in both the discovery and replication phase and further examined the results with meta-analysis. We detected association signals in the genes PGE receptor 3 (PTGER3) and 15-hydroxyprostaglandin dehydrogenase (HPGD), both strong biologic candidates for influence on polyp risk. We did not observe the previously reported effects and effect modification in PG-endoperoxide synthase 2 (PTGS2), PGE receptor 2 (PTGER2), or PGE receptor 4 (PTGER4), although we did observe a single nucleotide polymorphism in PTGER2 associated with risk of multiple adenomas. We also observed effect modification of the HPGD signal by NSAID exposure.

Genome-wide associated loci influencing interleukin (IL)-10, IL-1Ra, and IL-6 levels in African Americans

Interleukins (ILs) are key mediators of the immune response and inflammatory process. Plasma levels of IL-10, IL-1Ra, and IL-6 are associated with metabolic conditions, show large inter-individual variations, and are under strong genetic control. Therefore, elucidation of the genetic variants that influence levels of these ILs provides useful insights into mechanisms of immune response and pathogenesis of diseases. We conducted a genome-wide association study (GWAS) of IL-10, IL-1Ra, and IL-6 levels in 707 non-diabetic African Americans using 5,396,780 imputed and directly genotyped single nucleotide polymorphisms (SNPs) with adjustment for gender, age, and body mass index. IL-10 levels showed genome-wide significant associations (p  < 5 × 10(-8)) with eight SNPs, the most significant of which was rs5743185 in the PMS1 gene (p = 2.30 × 10(-10)). We tested replication of SNPs that showed genome-wide significance in 425 non-diabetic individuals from West Africa, and successfully replicated rs17365948 in the YWHAZ gene (p = 0.02). IL-1Ra levels showed suggestive associations with two SNPs in the ASB3 gene (p = 2.55 × 10(-7)), ten SNPs in the IL-1 gene family (IL1F5, IL1F8, IL1F10, and IL1Ra, p = 1.04 × 10(-6) to 1.75 × 10(-6)), and 23 SNPs near the IL1A gene (p = 1.22 × 10(-6) to 1.63 × 10(-6)). We also successfully replicated rs4251961 (p = 0.009); this SNP was reported to be associated with IL-1Ra levels in a candidate gene study of Europeans. IL-6 levels showed genome-wide significant association with one SNP (RP11-314E23.1; chr6:133397598; p = 8.63 × 10(-9)). To our knowledge, this is the first GWAS on IL-10, IL-1Ra, and IL-6 levels. Follow-up of these findings may provide valuable insight into the pathobiology of IL actions and dysregulations in inflammation and human diseases.

Comprehensive literature review and statistical considerations for GWAS meta-analysis

Over the last decade, genome-wide association studies (GWAS) have become the standard tool for gene discovery in human disease research. While debate continues about how to get the most out of these studies and on occasion about how much value these studies really provide, it is clear that many of the strongest results have come from large-scale mega-consortia and/or meta-analyses that combine data from up to dozens of studies and tens of thousands of subjects. While such analyses are becoming more and more common, statistical methods have lagged somewhat behind. There are good meta-analysis methods available, but even when they are carefully and optimally applied there remain some unresolved statistical issues. This article systematically reviews the GWAS meta-analysis literature, highlighting methodology and software options and reviewing methods that have been used in real studies. We illustrate differences among methods using a case study. We also discuss some of the unresolved issues and potential future directions.

Replication and meta-analysis of the gene-environment interaction between body mass index and the interleukin-6 promoter polymorphism with higher insulin resistance

Insulin resistance (IR) is a complex disorder caused by an interplay of both genetic and environmental factors. Recent studies identified a significant interaction between body mass index (BMI) and the rs1800795 polymorphism of the interleukin-6 gene that influences both IR and onset of type 2 diabetes mellitus, with obese individuals homozygous for the C allele demonstrating the highest level of IR and greatest risk for type 2 diabetes mellitus. Replication of a gene-environment interaction is important to confirm the validity of the initial finding and extend the generalizability of the results to other populations. Thus, the objective of this study was to replicate this gene-environment interaction on IR in a hypertensive population and perform a meta-analysis with prior published results. The replication analysis was performed using white individuals with hypertension from the Hypertensive Pathotype cohort (N = 311), genotyped for rs1800795. Phenotype studies were conducted after participants consumed 2 diets--high sodium (200 mmol/d) and low sodium (10 mmol/d)--for 7 days each. Measurements for plasma glucose, insulin, and interleukin-6 were obtained after 8 hours of fasting. Insulin resistance was characterized by the homeostatic model assessment (HOMA-IR). In Hypertensive Pathotype, BMI was a significant effect modifier of the relationship between rs1800795 and HOMA-IR; higher BMI was associated with higher HOMA-IR among homozygote CC individuals when compared with major allele G carriers (P = .003). Furthermore, the meta-analysis in 1028 individuals confirmed the result, demonstrating the same significant interaction between rs1800795 and BMI on HOMA-IR (P = 1.05 × 10(-6)). This rare replication of a gene-environment interaction extends the generalizability of the results to hypertension while highlighting this polymorphism as a marker of IR in obese individuals.

Common variants at 9p21 and 8q22 are associated with increased susceptibility to optic nerve degeneration in glaucoma

Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63-0.75], p = 1.86×10⁻¹⁸), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21-1.43], p = 3.87×10⁻¹¹). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50-0.67], p = 1.17×10⁻¹²) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53-0.72], p = 8.88×10⁻¹⁰). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41-0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54-1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma.

Evaluation of the metabochip genotyping array in African Americans and implications for fine mapping of GWAS-identified loci: the PAGE study

The Metabochip is a custom genotyping array designed for replication and fine mapping of metabolic, cardiovascular, and anthropometric trait loci and includes low frequency variation content identified from the 1000 Genomes Project. It has 196,725 SNPs concentrated in 257 genomic regions. We evaluated the Metabochip in 5,863 African Americans; 89% of all SNPs passed rigorous quality control with a call rate of 99.9%. Two examples illustrate the value of fine mapping with the Metabochip in African-ancestry populations. At CELSR2/PSRC1/SORT1, we found the strongest associated SNP for LDL-C to be rs12740374 (p = 3.5 × 10(-11)), a SNP indistinguishable from multiple SNPs in European ancestry samples due to high correlation. Its distinct signal supports functional studies elsewhere suggesting a causal role in LDL-C. At CETP we found rs17231520, with risk allele frequency 0.07 in African Americans, to be associated with HDL-C (p = 7.2 × 10(-36)). This variant is very rare in Europeans and not tagged in common GWAS arrays, but was identified as associated with HDL-C in African Americans in a single-gene study. Our results, one narrowing the risk interval and the other revealing an associated variant not found in Europeans, demonstrate the advantages of high-density genotyping of common and rare variation for fine mapping of trait loci in African American samples.

Functional variant in the autophagy-related 5 gene promotor is associated with childhood asthma

RATIONALE AND OBJECTIVE

Autophagy is a cellular process directed at eliminating or recycling cellular proteins. Recently, the autophagy pathway has been implicated in immune dysfunction, the pathogenesis of inflammatory disorders, and response to viral infection. Associations between two genes in the autophagy pathway, ATG5 and ATG7, with childhood asthma were investigated.

METHODS

Using genetic and experimental approaches, we examined the association of 13 HapMap-derived tagging SNPs in ATG5 and ATG7 with childhood asthma in 312 asthmatic and 246 non-allergic control children. We confirmed our findings by using independent cohorts and imputation analysis. Finally, we evaluated the functional relevance of a disease associated SNP.

MEASUREMENTS AND MAIN RESULTS

We demonstrated that ATG5 single nucleotide polymorphisms rs12201458 and rs510432 were associated with asthma (p = 0.00085 and 0.0025, respectively). In three independent cohorts, additional variants in ATG5 in the same LD block were associated with asthma (p<0.05). We found that rs510432 was functionally relevant and conferred significantly increased promotor activity. Furthermore, Atg5 expression was increased in nasal epithelium of acute asthmatics compared to stable asthmatics and non-asthmatic controls.

CONCLUSION

Genetic variants in ATG5, including a functional promotor variant, are associated with childhood asthma. These results provide novel evidence for a role for ATG5 in childhood asthma.

Evidence of differential allelic effects between adolescents and adults for plasma high-density lipoprotein

A recent meta-analysis of genome-wide association (GWA) studies identified 95 loci that influence lipid traits in the adult population and found that collectively these explained about 25–30% of heritability for each trait. Little is known about how these loci affect lipid levels in early life, but there is evidence that genetic effects on HDL- and LDL-cholesterol (HDL-C, LDL-C) and triglycerides vary with age. We studied Australian adults (N = 10,151) and adolescents (N = 2,363) who participated in twin and family studies and for whom we have lipid phenotypes and genotype information for 91 of the 95 genetic variants. Heterogeneity tests between effect sizes in adult and adolescent cohorts showed an excess of heterogeneity for HDL-C (p_Het<0.05 at 5 out of 37 loci), but no more than expected by chance for LDL-C (1 out of 14 loci), or trigycerides (0 out 24). There were 2 (out of 5) with opposite direction of effect in adolescents compared to adults for HDL-C, but none for LDL-C. The biggest difference in effect size was for LDL-C at rs6511720 near LDLR, adolescents (0.021±0.033 mmol/L) and adults (0.157±0.023 mmol/L), p_Het = 0.013; followed by ZNF664 (p_Het = 0.018) and PABPC4 (p_Het = 0.034) for HDL-C. Our findings suggest that some of the previously identified variants associate differently with lipid traits in adolescents compared to adults, either because of developmental changes or because of greater interactions with environmental differences in adults.

Translational studies of lipoprotein-associated phospholipase A₂ in inflammation and atherosclerosis

OBJECTIVES

This study sought to examine the role of lipoprotein-associated phospholipase A₂ (Lp-PLA₂/PLA2G7) in human inflammation and coronary atherosclerosis.

BACKGROUND

Lp-PLA₂ has emerged as a potential therapeutic target in coronary heart disease. Data supporting Lp-PLA₂ are indirect and confounded by species differences; whether Lp-PLA₂ is causal in coronary heart disease remains in question.

METHODS

We examined inflammatory regulation of Lp-PLA₂ during experimental endotoxemia in humans, probed the source of Lp-PLA₂ in human leukocytes under inflammatory conditions, and assessed the relationship of variation in PLA2G7, the gene encoding Lp-PLA₂, with coronary artery calcification.

RESULTS

In contrast to circulating tumor necrosis factor-alpha and C-reactive protein, blood and monocyte Lp-PLA₂ messenger ribonucleic acid decreased transiently, and plasma Lp-PLA₂ mass declined modestly during endotoxemia. In vitro, Lp-PLA₂ expression increased dramatically during human monocyte to macrophage differentiation and further in inflammatory macrophages and foamlike cells. Despite only a marginal association of single nucleotide polymorphisms in PLA2G7 with Lp-PLA₂ activity or mass, numerous PLA2G7 single nucleotide polymorphisms were associated with coronary artery calcification. In contrast, several single nucleotide polymorphisms in CRP were significantly associated with plasma C-reactive protein levels but had no relation with coronary artery calcification.

CONCLUSIONS

Circulating Lp-PLA₂ did not increase during acute phase response in humans, whereas inflammatory macrophages and foam cells, but not circulating monocytes, are major leukocyte sources of Lp-PLA₂. Common genetic variation in PLA2G7 is associated with subclinical coronary atherosclerosis. These data link Lp-PLA₂ to atherosclerosis in humans while highlighting the challenge in using circulating Lp-PLA₂ as a biomarker of Lp-PLA₂ actions in the vasculature.

Fasting glucose GWAS candidate region analysis across ethnic groups in the Multiethnic Study of Atherosclerosis (MESA)

Genetic variants associated with fasting glucose in European ancestry populations are increasingly well understood. However, the nature of the associations between these single nucleotide polymorphisms (SNPs) and fasting glucose in other racial and ethnic groups is unclear. We sought to examine regions previously identified to be associated with fasting glucose in Caucasian genome-wide association studies (GWAS) across multiple ethnicities in the Multiethnic Study of Atherosclerosis (MESA). Nondiabetic MESA participants with fasting glucose measured at the baseline exam and with GWAS genotyping were included; 2,349 Caucasians, 664 individuals of Chinese descent, 1,366 African Americans, and 1,171 Hispanics. Genotype data were generated from the Affymetrix 6.0 array and imputation in IMPUTE. Fasting glucose was regressed on SNP dosage data in each ethnic group adjusting for age, gender, MESA study center, and ethnic-specific principal components. SNPs from the three gene regions with the strongest associations to fasting glucose in previous Caucasian GWAS (MTNR1B / GCK / G6PC2) were examined in depth. There was limited power to replicate associations in other ethnic groups due to smaller allele frequencies and limited sample size; SNP associations may also have differed across ethnic groups due to differing linkage disequilibrium patterns with causal variants. rs10830963 in MTNR1B and rs4607517 in GCK demonstrated consistent magnitude and direction of association with fasting glucose across ethnic groups, although the associations were often not nominally significant. In conclusion, certain SNPs in MTNR1B and GCK demonstrate consistent effects across four racial and ethnic groups, narrowing the putative region for these causal variants.

Identification of common variants associated with human hippocampal and intracranial volumes

Identifying genetic variants influencing human brain structures may reveal new biological mechanisms underlying cognition and neuropsychiatric illness. The volume of the hippocampus is a biomarker of incipient Alzheimer's disease and is reduced in schizophrenia, major depression and mesial temporal lobe epilepsy. Whereas many brain imaging phenotypes are highly heritable, identifying and replicating genetic influences has been difficult, as small effects and the high costs of magnetic resonance imaging (MRI) have led to underpowered studies. Here we report genome-wide association meta-analyses and replication for mean bilateral hippocampal, total brain and intracranial volumes from a large multinational consortium. The intergenic variant rs7294919 was associated with hippocampal volume (12q24.22; N = 21,151; P = 6.70 × 10(-16)) and the expression levels of the positional candidate gene TESC in brain tissue. Additionally, rs10784502, located within HMGA2, was associated with intracranial volume (12q14.3; N = 15,782; P = 1.12 × 10(-12)). We also identified a suggestive association with total brain volume at rs10494373 within DDR2 (1q23.3; N = 6,500; P = 5.81 × 10(-7)).

Common variants at 12q14 and 12q24 are associated with hippocampal volume

Aging is associated with reductions in hippocampal volume (HV) that are accelerated by Alzheimer’s disease and vascular risk factors. Our genome-wide association study of dementia-free persons (n=9,232) identified 46 SNPs at four loci with p-values <4.0×10^-7. Two additional samples (n=2,318) replicated associations at 12q24 within MSRB3/WIF1 (discovery + replication, rs17178006; p=5.3×10^-11) and at 12q14 near HRK/FBXW8 (rs7294919; p=2.9×10^-11). Remaining associations included one 2q24 SNP within DPP4 (rs6741949; p=2.9×10^-7) and nine 9p33 SNPs within ASTN2 (rs7852872; p=1.0×10^-7) that were also associated with HV (p<0.05) in a third younger, more heterogeneous sample (n=7,794). The ASTN2 SNP was also associated with decline in cognition in a largely independent sample (n=1,563). These associations implicate genes related to apoptosis (HRK), development (WIF1), oxidative stress (MSR3B), ubiquitination (FBXW8), enzymes targeted by new diabetes medications (DPP4), and neuronal migration (ASTN2), indicating novel genetic influences that influence hippocampal size and possibly the risk of cognitive decline and dementia.

Methods for meta-analyses of genome-wide association studies: critical assessment of empirical evidence

There has been a steep increase in the number of meta-analyses of genome-wide association (GWA) studies aimed at identifying genetic variants with increasingly smaller effects, but pressure to publish findings of new genetic associations has limited the time available for careful consideration of all of their methodological aspects. The authors surveyed the literature (2007-2010) to provide empirical evidence on the methods used in GWA meta-analyses, including their organization, requirements about the uniformity of methods used in primary studies, methods for data pooling, investigation of between-study heterogeneity, and quality of reporting. This review showed that a great variety of methods are being used, but the rationale for their choice is often unclear. It also highlights how important methodological aspects have received insufficient attention, potentially leading to missed opportunities for improving gene discovery and characterization. Evaluation of power to replicate findings was inadequate, and the number of variants selected for replication was not associated with replication sample size. A low proportion of GWA meta-analyses investigated the presence and magnitude of heterogeneity, even when there was little uniformity in methods used in primary studies. More methodological work is required before clear guidance can be offered as to optimal methods or tradeoffs between alternative methods. However, there is a clear need for guidelines for reporting the results of GWA meta-analyses.

Population structure of Hispanics in the United States: the multi-ethnic study of atherosclerosis

Using ~60,000 SNPs selected for minimal linkage disequilibrium, we perform population structure analysis of 1,374 unrelated Hispanic individuals from the Multi-Ethnic Study of Atherosclerosis (MESA), with self-identification corresponding to Central America (n = 93), Cuba (n = 50), the Dominican Republic (n = 203), Mexico (n = 708), Puerto Rico (n = 192), and South America (n = 111). By projection of principal components (PCs) of ancestry to samples from the HapMap phase III and the Human Genome Diversity Panel (HGDP), we show the first two PCs quantify the Caucasian, African, and Native American origins, while the third and fourth PCs bring out an axis that aligns with known South-to-North geographic location of HGDP Native American samples and further separates MESA Mexican versus Central/South American samples along the same axis. Using k-means clustering computed from the first four PCs, we define four subgroups of the MESA Hispanic cohort that show close agreement with self-identification, labeling the clusters as primarily Dominican/Cuban, Mexican, Central/South American, and Puerto Rican. To demonstrate our recommendations for genetic analysis in the MESA Hispanic cohort, we present pooled and stratified association analysis of triglycerides for selected SNPs in the LPL and TRIB1 gene regions, previously reported in GWAS of triglycerides in Caucasians but as yet unconfirmed in Hispanic populations. We report statistically significant evidence for genetic association in both genes, and we further demonstrate the importance of considering population substructure and genetic heterogeneity in genetic association studies performed in the United States Hispanic population.

Alarcón-Riquelme for the BIOLUPUS and GENLES Networks, Alarcón G, Anaya JM, Bae SC, Boackle S, Brown E, Chang DM, Cho SK, Criswell L, Edberg J, Freedman B, Gilkeson G, Jacob C, James J, Kamen D, Kimberly R, Kim JH, Martin J, Merrill J, Niewold T, Park SY, Petri M, Pons-Estel B, Ramsey-Goldman R, Reveille J, Scofield R, Song Y, Stevens A, Tsao B, Vila L, Vyse T, Yu CY, Guthridge J, Kaufman K, Harley J, Wakeland E, Langefeld C, Gaffney P, Montgomery C, Moser K. Identification of IRF8, TMEM39A, and IKZF3-ZPBP2 as susceptibility loci for systemic lupus erythematosus in a large-scale multiracial replication study

Systemic lupus erythematosus (SLE) is a chronic heterogeneous autoimmune disorder characterized by the loss of tolerance to self-antigens and dysregulated interferon responses. The etiology of SLE is complex, involving both heritable and environmental factors. Candidate-gene studies and genome-wide association (GWA) scans have been successful in identifying new loci that contribute to disease susceptibility; however, much of the heritable risk has yet to be identified. In this study, we sought to replicate 1,580 variants showing suggestive association with SLE in a previously published GWA scan of European Americans; we tested a multiethnic population consisting of 7,998 SLE cases and 7,492 controls of European, African American, Asian, Hispanic, Gullah, and Amerindian ancestry to find association with the disease. Several genes relevant to immunological pathways showed association with SLE. Three loci exceeded the genome-wide significance threshold: interferon regulatory factor 8 (IRF8; rs11644034; p(meta-Euro) = 2.08 × 10(-10)), transmembrane protein 39A (TMEM39A; rs1132200; p(meta-all) = 8.62 × 10(-9)), and 17q21 (rs1453560; p(meta-all) = 3.48 × 10(-10)) between IKAROS family of zinc finger 3 (AIOLOS; IKZF3) and zona pellucida binding protein 2 (ZPBP2). Fine mapping, resequencing, imputation, and haplotype analysis of IRF8 indicated that three independent effects tagged by rs8046526, rs450443, and rs4843869, respectively, were required for risk in individuals of European ancestry. Eleven additional replicated effects (5 × 10(-8) < p(meta-Euro) < 9.99 × 10(-5)) were observed with CFHR1, CADM2, LOC730109/IL12A, LPP, LOC63920, SLU7, ADAMTSL1, C10orf64, OR8D4, FAM19A2, and STXBP6. The results of this study increase the number of confirmed SLE risk loci and identify others warranting further investigation.

Role of common and rare APP DNA sequence variants in Alzheimer disease

OBJECTIVES

More than 30 different rare mutations, including copy number variants (CNVs), in the amyloid precursor protein gene (APP) cause early-onset familial Alzheimer disease (EOFAD), whereas the contribution of common APP variants to disease risk remains controversial. In this study we systematically assessed the role of both rare and common APP DNA variants in Alzheimer disease (AD) families.

METHODS

Families with EOFAD genetically linked to the APP region were screened for missense mutations and locus duplications of APP. Further, using genome-wide DNA microarray data, we examined the APP locus for CNVs in a total of 797 additional early- and late-onset AD pedigrees. Finally, 423 single nucleotide polymorphisms (SNPs) in the APP locus, including 2 promoter polymorphisms previously associated with AD risk, were tested in up to 4,200 individuals from multiplex AD families.

RESULTS

Analyses of 8 21q21-linked families revealed one family carrying a nonsynonymous mutation in exon 17 (Val717Leu) and another family with a partially penetrant 3.5-Mb locus duplication encompassing APP. CNV analysis in the APP locus revealed an additional family carrying a fully penetrant 380-kb duplication, merely spanning APP. Last, contrary to previous reports, association analyses of more than 400 different SNPs in or near APP failed to show significant effects on AD risk.

CONCLUSION

Our study shows that APP mutations and locus duplications are a very rare cause of EOFAD and that the contribution of common APP variants to AD susceptibility is insignificant. Furthermore, duplications of APP may not be fully penetrant, possibly indicating the existence of hitherto unknown protective genetic factors.

Finding genes and variants for lipid levels after genome-wide association analysis

PURPOSE OF REVIEW

We review the main findings from genome-wide association studies (GWAS) for levels of HDL-cholesterol, LDL-cholesterol and triglycerides, including approaches to identify the functional variant(s) or gene(s). We discuss study design and challenges related to whole genome or exome sequencing to identify novel genes and variants.

RECENT FINDINGS

GWAS have detected approximately 100 loci associated with one or more lipid trait. Fine mapping of several loci for LDL-cholesterol demonstrated that the trait variance explained may double when the functional variants responsible for the association signals are identified. Experimental follow-up of three loci identified by GWAS has identified functional genes GALNT2, TRIB1, and SORT1, and a functional variant at SORT1.

SUMMARY

The goal of genetic studies for lipid levels is to improve treatment and ultimately reduce the prevalence of heart disease. Many signals identified by GWAS have modest effect sizes, useful for identifying novel biologically relevant genes, but less useful for personalized medicine. Whole genome or exome sequencing studies may fill this gap by identifying rare variants of larger effect associated with lipid levels and heart disease.

Endophenotyping reveals differential phenotype-genotype correlations between myopia-associated polymorphisms and eye biometric parameters

PURPOSE

To investigate the association with ocular biometric parameters in myopia-associated single nucleotide polymorphisms (SNPs) of the gap junction protein delta 2 (GJD2), insulin-like growth factor-1 (IGF1) and hepatocyte growth factor (HGF) genes in two geographically different Chinese cohorts.

METHODS

In 814 unrelated Han Chinese individuals aged above 50 years including 362 inland residents and 432 island dwellers, comprehensive ophthalmic examinations were performed. Three SNPs, including GJD2 rs634990, IGF1 rs6214, and HGF rs3735520, were genotyped. Genetic association with ocular biometric parameters was analyzed in individual cohorts, using linear regression controlled for sex and age. Common associations shared by the two cohorts were revealed by meta-analysis.

RESULTS

Meta-analysis showed that GJD2 rs634990 alone was not associated with any biometric parameters (adjusted p>0.645). The T allele of IGF1 rs6214 was specifically associated with thicker lens (β±SE=0.055±0.022, adjusted p=0.034). The A allele of HGF rs3735520 was associated with longer vitreous chamber depth (β±SE=0.143±0.060, adjusted p=0.050). Significant interaction between HGF rs3735520 and GJD2 rs634990 was found in association with axial length and vitreous chamber depth (adjusted p=0.003 and 0.033, respectively), and possibly with spherical error (adjusted p=0.056).

CONCLUSIONS

Our endophenotyping analysis showed differential association between selected myopia-associated genes and ocular biometric parameters in our Chinese cohorts, which may underline substantial but diversified effects of these genes and their interaction on the development of eye structure and etiology of myopia.

Genome-wide association and functional follow-up reveals new loci for kidney function

Chronic kidney disease (CKD) is an important public health problem with a genetic component. We performed genome-wide association studies in up to 130,600 European ancestry participants overall, and stratified for key CKD risk factors. We uncovered 6 new loci in association with estimated glomerular filtration rate (eGFR), the primary clinical measure of CKD, in or near MPPED2, DDX1, SLC47A1, CDK12, CASP9, and INO80. Morpholino knockdown of mpped2 and casp9 in zebrafish embryos revealed podocyte and tubular abnormalities with altered dextran clearance, suggesting a role for these genes in renal function. By providing new insights into genes that regulate renal function, these results could further our understanding of the pathogenesis of CKD.

Associations between incident ischemic stroke events and stroke and cardiovascular disease-related genome-wide association studies single nucleotide polymorphisms in the Population Architecture Using Genomics and Epidemiology study

BACKGROUND

Genome-wide association studies (GWAS) have identified loci associated with ischemic stroke (IS) and cardiovascular disease (CVD) in European-descent individuals, but their replication in different populations has been largely unexplored.

METHODS AND RESULTS

Nine single nucleotide polymorphisms (SNPs) selected from GWAS and meta-analyses of stroke, and 86 SNPs previously associated with myocardial infarction and CVD risk factors, including blood lipids (high density lipoprotein [HDL], low density lipoprotein [LDL], and triglycerides), type 2 diabetes, and body mass index (BMI), were investigated for associations with incident IS in European Americans (EA) N=26 276, African-Americans (AA) N=8970, and American Indians (AI) N=3570 from the Population Architecture using Genomics and Epidemiology Study. Ancestry-specific fixed effects meta-analysis with inverse variance weighting was used to combine study-specific log hazard ratios from Cox proportional hazards models. Two of 9 stroke SNPs (rs783396 and rs1804689) were significantly associated with [corrected] IS hazard in AA; none were significant in this large EA cohort. Of 73 CVD risk factor SNPs tested in EA, 2 (HDL and triglycerides SNPs) were associated with IS. In AA, SNPs associated with LDL, HDL, and BMI were significantly associated with IS (3 of 86 SNPs tested). Out of 58 SNPs tested in AI, 1 LDL SNP was significantly associated with IS.

CONCLUSIONS

Our analyses showing lack of replication in spite of reasonable power for many stroke SNPs and differing results by ancestry highlight the need to follow up on GWAS findings and conduct genetic association studies in diverse populations. We found modest IS associations with BMI and lipids SNPs, though these findings require confirmation.

A meta-analysis and genome-wide association study of platelet count and mean platelet volume in african americans

Several genetic variants associated with platelet count and mean platelet volume (MPV) were recently reported in people of European ancestry. In this meta-analysis of 7 genome-wide association studies (GWAS) enrolling African Americans, our aim was to identify novel genetic variants associated with platelet count and MPV. For all cohorts, GWAS analysis was performed using additive models after adjusting for age, sex, and population stratification. For both platelet phenotypes, meta-analyses were conducted using inverse-variance weighted fixed-effect models. Platelet aggregation assays in whole blood were performed in the participants of the GeneSTAR cohort. Genetic variants in ten independent regions were associated with platelet count (N = 16,388) with p<5×10(-8) of which 5 have not been associated with platelet count in previous GWAS. The novel genetic variants associated with platelet count were in the following regions (the most significant SNP, closest gene, and p-value): 6p22 (rs12526480, LRRC16A, p = 9.1×10(-9)), 7q11 (rs13236689, CD36, p = 2.8×10(-9)), 10q21 (rs7896518, JMJD1C, p = 2.3×10(-12)), 11q13 (rs477895, BAD, p = 4.9×10(-8)), and 20q13 (rs151361, SLMO2, p = 9.4×10(-9)). Three of these loci (10q21, 11q13, and 20q13) were replicated in European Americans (N = 14,909) and one (11q13) in Hispanic Americans (N = 3,462). For MPV (N = 4,531), genetic variants in 3 regions were significant at p<5×10(-8), two of which were also associated with platelet count. Previously reported regions that were also significant in this study were 6p21, 6q23, 7q22, 12q24, and 19p13 for platelet count and 7q22, 17q11, and 19p13 for MPV. The most significant SNP in 1 region was also associated with ADP-induced maximal platelet aggregation in whole blood (12q24). Thus through a meta-analysis of GWAS enrolling African Americans, we have identified 5 novel regions associated with platelet count of which 3 were replicated in other ethnic groups. In addition, we also found one region associated with platelet aggregation that may play a potential role in atherothrombosis.

A genome-wide scan of Ashkenazi Jewish Crohn's disease suggests novel susceptibility loci

Crohn's disease (CD) is a complex disorder resulting from the interaction of intestinal microbiota with the host immune system in genetically susceptible individuals. The largest meta-analysis of genome-wide association to date identified 71 CD-susceptibility loci in individuals of European ancestry. An important epidemiological feature of CD is that it is 2-4 times more prevalent among individuals of Ashkenazi Jewish (AJ) descent compared to non-Jewish Europeans (NJ). To explore genetic variation associated with CD in AJs, we conducted a genome-wide association study (GWAS) by combining raw genotype data across 10 AJ cohorts consisting of 907 cases and 2,345 controls in the discovery stage, followed up by a replication study in 971 cases and 2,124 controls. We confirmed genome-wide significant associations of 9 known CD loci in AJs and replicated 3 additional loci with strong signal (p<5×10⁻⁶). Novel signals detected among AJs were mapped to chromosomes 5q21.1 (rs7705924, combined p = 2×10⁻⁸; combined odds ratio OR = 1.48), 2p15 (rs6545946, p = 7×10⁻⁹; OR = 1.16), 8q21.11 (rs12677663, p = 2×10⁻⁸; OR = 1.15), 10q26.3 (rs10734105, p = 3×10⁻⁸; OR = 1.27), and 11q12.1 (rs11229030, p = 8×10⁻⁹; OR = 1.15), implicating biologically plausible candidate genes, including RPL7, CPAMD8, PRG2, and PRG3. In all, the 16 replicated and newly discovered loci, in addition to the three coding NOD2 variants, accounted for 11.2% of the total genetic variance for CD risk in the AJ population. This study demonstrates the complementary value of genetic studies in the Ashkenazim.

Evolutionary conserved longevity genes and human cognitive abilities in elderly cohorts

Genetic influences have an important role in the ageing process. The genetic factors that influence success in bodily ageing may also contribute to the successful ageing of cognitive abilities. A comparative genomics approach found longevity genes conserved between yeast Saccharomyces cerevisiae and nematode Caenorhabditis elegans. We hypothesised that these longevity genes influence variance in cognitive ability and age-related cognitive decline in humans. Here, we investigated six of these genes that have human orthologs and show expression in the brain. We tested AFG3L2 (MIM: 604581, AFG3 ATPase family gene 3-like 2 (yeast)), FRAP1 (MIM: 601231, a FK506 binding protein 12-rapamycin associated protein), MAT1A, MAT2A (MIM: 610550 and 601468, methionine adenosyltransferases I alpha and II alpha, respectively), SYNJ1 and SYNJ2 (MIM: 604297 and 609410, synaptojanin-1 and synaptojanin-2, respectively) in approximately 1000 healthy older Scots: the Lothian Birth Cohort 1936 (LBC1936). They were tested on general cognitive ability at age 11 years. At a mean age of 70 years, they re-sat the same general cognitive ability test and underwent an additional battery of diverse cognitive tests. In all, 70 tag and functional SNPs in the six longevity genes were genotyped and tested for association with cognition and cognitive ageing in LBC1936. Suggestive associations were detected between SNPs in SYNJ2, MAT1A, AFG3L2 and SYNJ1 and a general memory factor and general cognitive ability at age 11 and 70 years. Replication studies for cognitive ability associations were performed in 2506 samples from the Cognitive Ageing Genetics in England and Scotland consortium. A meta-analysis replicated the SYNJ2 association with cognitive abilities (lowest P=0.00077). SYNJ2 is a novel gene in which variation is potentially associated with cognitive abilities.

Longevity candidate genes and their association with personality traits in the elderly

Human longevity and personality traits are both heritable and are consistently linked at the phenotypic level. We test the hypothesis that candidate genes influencing longevity in lower organisms are associated with variance in the five major dimensions of human personality (measured by the NEO-FFI and IPIP inventories) plus related mood states of anxiety and depression. Seventy single nucleotide polymorphisms (SNPs) in six brain expressed, longevity candidate genes (AFG3L2, FRAP1, MAT1A, MAT2A, SYNJ1, and SYNJ2) were typed in over 1,000 70-year old participants from the Lothian Birth Cohort of 1936 (LBC1936). No SNPs were associated with the personality and psychological distress traits at a Bonferroni corrected level of significance (P < 0.0002), but there was an over-representation of nominally significant (P < 0.05) SNPs in the synaptojanin-2 (SYNJ2) gene associated with agreeableness and symptoms of depression. Eight SNPs which showed nominally significant association across personality measurement instruments were tested in an extremely large replication sample of 17,106 participants. SNP rs350292, in SYNJ2, was significant: the minor allele was associated with an average decrease in NEO agreeableness scale scores of 0.25 points, and 0.67 points in the restricted analysis of elderly cohorts (most aged >60 years). Because we selected a specific set of longevity genes based on functional genomics findings, further research on other longevity gene candidates is warranted to discover whether they are relevant candidates for personality and psychological distress traits.

Genome-wide association study identifies novel loci associated with circulating phospho- and sphingolipid concentrations

Phospho- and sphingolipids are crucial cellular and intracellular compounds. These lipids are required for active transport, a number of enzymatic processes, membrane formation, and cell signalling. Disruption of their metabolism leads to several diseases, with diverse neurological, psychiatric, and metabolic consequences. A large number of phospholipid and sphingolipid species can be detected and measured in human plasma. We conducted a meta-analysis of five European family-based genome-wide association studies (N = 4034) on plasma levels of 24 sphingomyelins (SPM), 9 ceramides (CER), 57 phosphatidylcholines (PC), 20 lysophosphatidylcholines (LPC), 27 phosphatidylethanolamines (PE), and 16 PE-based plasmalogens (PLPE), as well as their proportions in each major class. This effort yielded 25 genome-wide significant loci for phospholipids (smallest P-value = 9.88×10(-204)) and 10 loci for sphingolipids (smallest P-value = 3.10×10(-57)). After a correction for multiple comparisons (P-value<2.2×10(-9)), we observed four novel loci significantly associated with phospholipids (PAQR9, AGPAT1, PKD2L1, PDXDC1) and two with sphingolipids (PLD2 and APOE) explaining up to 3.1% of the variance. Further analysis of the top findings with respect to within class molar proportions uncovered three additional loci for phospholipids (PNLIPRP2, PCDH20, and ABDH3) suggesting their involvement in either fatty acid elongation/saturation processes or fatty acid specific turnover mechanisms. Among those, 14 loci (KCNH7, AGPAT1, PNLIPRP2, SYT9, FADS1-2-3, DLG2, APOA1, ELOVL2, CDK17, LIPC, PDXDC1, PLD2, LASS4, and APOE) mapped into the glycerophospholipid and 12 loci (ILKAP, ITGA9, AGPAT1, FADS1-2-3, APOA1, PCDH20, LIPC, PDXDC1, SGPP1, APOE, LASS4, and PLD2) to the sphingolipid pathways. In large meta-analyses, associations between FADS1-2-3 and carotid intima media thickness, AGPAT1 and type 2 diabetes, and APOA1 and coronary artery disease were observed. In conclusion, our study identified nine novel phospho- and sphingolipid loci, substantially increasing our knowledge of the genetic basis for these traits.

Meta-analysis of untargeted metabolomic data from multiple profiling experiments

metaXCMS is a software program for the analysis of liquid chromatography/mass spectrometry-based untargeted metabolomic data. It is designed to identify the differences between metabolic profiles across multiple sample groups (e.g., 'healthy' versus 'active disease' versus 'inactive disease'). Although performing pairwise comparisons alone can provide physiologically relevant data, these experiments often result in hundreds of differences, and comparison with additional biologically meaningful sample groups can allow for substantial data reduction. By performing second-order (meta-) analysis, metaXCMS facilitates the prioritization of interesting metabolite features from large untargeted metabolomic data sets before the rate-limiting step of structural identification. Here we provide a detailed step-by-step protocol for going from raw mass spectrometry data to metaXCMS results, visualized as Venn diagrams and exported Microsoft Excel spreadsheets. There is no upper limit to the number of sample groups or individual samples that can be compared with the software, and data from most commercial mass spectrometers are supported. The speed of the analysis depends on computational resources and data volume, but will generally be less than 1 d for most users. metaXCMS is freely available at http://metlin.scripps.edu/metaxcms/.

A genome-wide association study of COPD identifies a susceptibility locus on chromosome 19q13

The genetic risk factors for chronic obstructive pulmonary disease (COPD) are still largely unknown. To date, genome-wide association studies (GWASs) of limited size have identified several novel risk loci for COPD at CHRNA3/CHRNA5/IREB2, HHIP and FAM13A; additional loci may be identified through larger studies. We performed a GWAS using a total of 3499 cases and 1922 control subjects from four cohorts: the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE); the Normative Aging Study (NAS) and National Emphysema Treatment Trial (NETT); Bergen, Norway (GenKOLS); and the COPDGene study. Genotyping was performed on Illumina platforms with additional markers imputed using 1000 Genomes data; results were summarized using fixed-effect meta-analysis. We identified a new genome-wide significant locus on chromosome 19q13 (rs7937, OR = 0.74, P = 2.9 × 10(-9)). Genotyping this single nucleotide polymorphism (SNP) and another nearby SNP in linkage disequilibrium (rs2604894) in 2859 subjects from the family-based International COPD Genetics Network study (ICGN) demonstrated supportive evidence for association for COPD (P = 0.28 and 0.11 for rs7937 and rs2604894), pre-bronchodilator FEV(1) (P = 0.08 and 0.04) and severe (GOLD 3&4) COPD (P = 0.09 and 0.017). This region includes RAB4B, EGLN2, MIA and CYP2A6, and has previously been identified in association with cigarette smoking behavior.

Common variants near MBNL1 and NKX2-5 are associated with infantile hypertrophic pyloric stenosis

Infantile hypertrophic pyloric stenosis (IHPS) is a severe condition characterized by hypertrophy of the pyloric sphincter muscle. We conducted a genome-wide association study (GWAS) on 1,001 surgery-confirmed cases and 2,401 controls from Denmark. The six most strongly associated loci were tested in a replication set of 796 cases and 876 controls. Three SNPs reached genome-wide significance. One of these SNPs, rs11712066 (odds ratio (OR) = 1.61; P = 1.5 × 10(-17)) at 3p25.1, is located 150 kb upstream of MBNL1, which encodes a factor that regulates splicing transitions occurring shortly after birth. The second SNP, rs573872 (OR = 1.41; P = 4.3 × 10(-12)), maps to an intergenic region at 3p25.2 approximately 1.3 Mb downstream of MBNL1. The third SNP, rs29784 (OR = 1.42; P = 1.5 × 10(-15)) at 5q35.2, is 64 kb downstream of NKX2-5, which is involved in development of cardiac muscle tissue and embryonic gut development.

SNP rs11190870 near LBX1 is associated with adolescent idiopathic scoliosis in southern Chinese

A study was conducted to validate the most significant single nucleotide polymorphism (SNP) from a genome-wide association study of Japanese adolescent idiopathic scoliosis (AIS) patients in an independent southern Chinese population. In total, 300 AIS patients fulfilled the clinical criteria and 788 controls with MRI scans of the spine were included in the replication study. We employed case-control analysis to study the association of SNP rs11190870 near LBX1 (ladybird homeobox 1) with AIS in a southern Chinese population. The results suggest that SNP rs11190870 is significantly associated with AIS (P=9.1 × 10(-10); odds ratio=1.85; 95% confidence interval=1.52-2.25). The results of this study confirm that SNP rs11190870 is associated with AIS.

A comprehensive genetic association study of Alzheimer disease in African Americans

OBJECTIVES

To evaluate the association of genetic variation with late-onset Alzheimer disease (AD) in African Americans, including genes implicated in recent genome-wide association studies of whites.

DESIGN

We analyzed a genome-wide set of 2.5 million imputed markers to evaluate the genetic basis of AD in an African American population.

SUBJECTS

Five hundred thirteen well-characterized African American AD cases and 496 cognitively normal African American control subjects.

SETTING

Data were collected from multiple sites as part of the Multi-Institutional Research on Alzheimer Genetic Epidemiology (MIRAGE) Study and the Henry Ford Health System as part of the Genetic and Environmental Risk Factors for Alzheimer Disease Among African Americans (GenerAAtions) Study.

RESULTS

Several significant single-nucleotide polymorphisms (SNPs) were observed in the region of the apolipoprotein E gene (APOE). After adjusting for the confounding effects of APOE genotype, one of these SNPs, rs6859 in PVRL2, remained significantly associated with AD (P = .0087). Association was also observed with SNPs in CLU, PICALM, BIN1, EPHA1, MS4A, ABCA7, and CD33, although the effect direction for some SNPs and the most significant SNPs differed from findings in data sets consisting of whites. Finally, using the African American genome-wide association study data set as a discovery sample, we obtained suggestive evidence of association with SNPs for several novel candidate genes.

CONCLUSIONS

Some genes contribute to AD pathogenesis in both white and African American cohorts, although it is unclear whether the causal variants are the same. A larger African American sample will be needed to confirm novel gene associations, which may be population specific.

Analysis of family- and population-based samples in cohort genome-wide association studies

Cohort studies typically sample unrelated individuals from a population, although family members of index cases may also be recruited to investigate shared familial risk factors. Recruitment of family members may be incomplete or ancillary to the main cohort, resulting in a mixed sample of independent family units, including unrelated singletons and multiplex families. Multiple methods are available to perform genome-wide association (GWA) analysis of binary or continuous traits in families, but it is unclear whether methods known to perform well on ascertained pedigrees, sibships, or trios are appropriate in analysis of a mixed unrelated cohort and family sample. We present simulation studies based on Multi-Ethnic Study of Atherosclerosis (MESA) pedigree structures to compare the performance of several popular methods of GWA analysis for both quantitative and dichotomous traits in cohort studies. We evaluate approaches suitable for analysis of families, and combined the best performing methods with population-based samples either by meta-analysis, or by pooled analysis of family- and population-based samples (mega-analysis), comparing type 1 error and power. We further assess practical considerations, such as availability of software and ability to incorporate covariates in statistical modeling, and demonstrate our recommended approaches through quantitative and binary trait analysis of HDL cholesterol (HDL-C) in 2,553 MESA family- and population-based African-American samples. Our results suggest linear modeling approaches that accommodate family-induced phenotypic correlation (e.g., variance-component model for quantitative traits or generalized estimating equations for dichotomous traits) perform best in the context of combined family- and population-based cohort GWAS.

A genome-wide association study of gestational diabetes mellitus in Korean women

Knowledge regarding the genetic risk loci for gestational diabetes mellitus (GDM) is still limited. In this study, we performed a two-stage genome-wide association analysis in Korean women. In the stage 1 genome scan, 468 women with GDM and 1,242 nondiabetic control women were compared using 2.19 million genotyped or imputed markers. We selected 11 loci for further genotyping in stage 2 samples of 931 case and 783 control subjects. The joint effect of stage 1 plus stage 2 studies was analyzed by meta-analysis. We also investigated the effect of known type 2 diabetes variants in GDM. Two loci known to be associated with type 2 diabetes had a genome-wide significant association with GDM in the joint analysis. rs7754840, a variant in CDKAL1, had the strongest association with GDM (odds ratio 1.518; P=6.65×10(-16)). A variant near MTNR1B, rs10830962, was also significantly associated with the risk of GDM (1.454; P=2.49×10(-13)). We found that there is an excess of association between known type 2 diabetes variants and GDM above what is expected under the null hypothesis. In conclusion, we have confirmed that genetic variants in CDKAL1 and near MTNR1B are strongly associated with GDM in Korean women. There seems to be a shared genetic basis between GDM and type 2 diabetes.

Genome-wide association study of body height in African Americans: the Women's Health Initiative SNP Health Association Resource (SHARe)

Height is a complex trait under strong genetic influence. To date, numerous genetic loci have been associated with height in individuals of European ancestry. However, few large-scale discovery genome-wide association studies (GWAS) of height in minority populations have been conducted and thus information about population-specific height regulation is limited. We conducted a GWA analysis of height in 8149 African-American (AA) women from the Women's Health Initiative. Genetic variants with P< 5 × 10(-5) (n = 169) were followed up in a replication data set (n = 20 809) and meta-analyzed in a total of 28 958 AAs and African-descent individuals. Twelve single-nucleotide polymorphisms (SNPs) representing 7 independent loci were significantly associated with height at P < 5 × 10(-8). We identified novel SNPs in 17q23 (TMEM100/PCTP) and Xp22.3 (ARSE) reflecting population-specific regulation of height in AAs and replicated five loci previously reported in European-descent populations [4p15/LCORL, 11q13/SERPINH1, 12q14/HMGA2, 17q23/MAP3K3 (mitogen-activated protein kinase3) and 18q21/DYM]. In addition, we performed an admixture mapping analysis of height which is both complementary and supportive to the GWA analysis and suggests potential associations between ancestry and height on chromosomes 4 (4q21), 15 (15q26) and 17 (17q23). Our findings provide insight into the genetic architecture of height and support the investigation of non-European-descent populations for identifying genetic factors associated with complex traits. Specifically, we identify new loci that may reflect population-specific regulation of height and report several known height loci that are important in determining height in African-descent populations.

A functional haplotype in EIF2AK3, an ER stress sensor, is associated with lower bone mineral density

EIF2AK3 is a type I transmembrane protein that functions as an endoplasmic reticulum (ER) stress sensor to regulate global protein synthesis. Rare mutations in EIF2AK3 cause Wolcott-Rallison syndrome (OMIM 226980), an autosomal recessive disorder characterized by diabetes, epiphyseal dysplasia, osteoporosis, and growth retardation. To investigate the role of common genetic variation in EIF2AK3 as a determinant of bone mineral density (BMD) and osteoporosis, we sequenced all exons and flanking regions, then genotyped six potentially functional single nucleotide polymorphisms (SNPs) in this gene in 997 Amish subjects for association analysis, and attempted replication in 887 Mexican Americans. We found that the minor allele of a nonsynonymous SNP rs13045 had borderline associations with decreased forearm BMD in both discovery and replication cohorts (unadjusted p = 0.036 and β = -0.007 for the Amish; unadjusted p = 0.031 and β = -0.008 for Mexican Americans). A meta-analysis indicated this association achieved statistical significance in the combined sample (unadjusted p = 0.003; Bonferroni corrected p = 0.009). Rs13045 and three other potentially functional SNPs, a promoter SNP (rs6547787) and two nonsynonymous SNPs (rs867529 and rs1805165), formed two haplotypes: a low-BMD associated haplotype, denoted haplotype B [minor allele frequency (MAF) = 0.311] and a common haplotype A (MAF = 0.676). There were no differences in mRNA expression in lymphoblastoid cell lines between the two haplotypes. However, after treating lymphoblastoid cell lines with thapsigargin to induce ER stress, cell lines with haplotype B showed increased sensitivity to ER stress (p = 0.014) compared with cell lines with haplotype A. Taken together, our results suggest that common nonsynonymous sequence variants in EIF2AK3 have a modest effect on ER stress response and may contribute to the risk for low BMD through this mechanism.

Association between chromosome 9p21 variants and the ankle-brachial index identified by a meta-analysis of 21 genome-wide association studies

BACKGROUND

Genetic determinants of peripheral arterial disease (PAD) remain largely unknown. To identify genetic variants associated with the ankle-brachial index (ABI), a noninvasive measure of PAD, we conducted a meta-analysis of genome-wide association study data from 21 population-based cohorts.

METHODS AND RESULTS

Continuous ABI and PAD (ABI ≤0.9) phenotypes adjusted for age and sex were examined. Each study conducted genotyping and imputed data to the ≈2.5 million single nucleotide polymorphisms (SNPs) in HapMap. Linear and logistic regression models were used to test each SNP for association with ABI and PAD using additive genetic models. Study-specific data were combined using fixed effects inverse variance weighted meta-analyses. There were a total of 41 692 participants of European ancestry (≈60% women, mean ABI 1.02 to 1.19), including 3409 participants with PAD and with genome-wide association study data available. In the discovery meta-analysis, rs10757269 on chromosome 9 near CDKN2B had the strongest association with ABI (β=-0.006, P=2.46×10(-8)). We sought replication of the 6 strongest SNP associations in 5 population-based studies and 3 clinical samples (n=16 717). The association for rs10757269 strengthened in the combined discovery and replication analysis (P=2.65×10(-9)). No other SNP associations for ABI or PAD achieved genome-wide significance. However, 2 previously reported candidate genes for PAD and 1 SNP associated with coronary artery disease were associated with ABI: DAB21P (rs13290547, P=3.6×10(-5)), CYBA (rs3794624, P=6.3×10(-5)), and rs1122608 (LDLR, P=0.0026).

CONCLUSIONS

Genome-wide association studies in more than 40 000 individuals identified 1 genome wide significant association on chromosome 9p21 with ABI. Two candidate genes for PAD and 1 SNP for coronary artery disease are associated with ABI.

BioSMACK: a linux live CD for genome-wide association analyses

Recent advances in high-throughput genotyping technologies have enabled us to conduct a genome-wide association study (GWAS) on a large cohort. However, analyzing millions of single nucleotide polymorphisms (SNPs) is still a difficult task for researchers conducting a GWAS. Several difficulties such as compatibilities and dependencies are often encountered by researchers using analytical tools, during the installation of software. This is a huge obstacle to any research institute without computing facilities and specialists. Therefore, a proper research environment is an urgent need for researchers working on GWAS. We developed BioSMACK to provide a research environment for GWAS that requires no configuration and is easy to use. BioSMACK is based on the Ubuntu Live CD that offers a complete Linux-based operating system environment without installation. Moreover, we provide users with a GWAS manual consisting of a series of guidelines for GWAS and useful examples. BioSMACK is freely available at http://ksnp.cdc. go.kr/biosmack.

Meta-analyses identify 13 loci associated with age at menopause and highlight DNA repair and immune pathways

To newly identify loci for age at natural menopause, we carried out a meta-analysis of 22 genome-wide association studies (GWAS) in 38,968 women of European descent, with replication in up to 14,435 women. In addition to four known loci, we identified 13 loci newly associated with age at natural menopause (at P < 5 × 10(-8)). Candidate genes located at these newly associated loci include genes implicated in DNA repair (EXO1, HELQ, UIMC1, FAM175A, FANCI, TLK1, POLG and PRIM1) and immune function (IL11, NLRP11 and PRRC2A (also known as BAT2)). Gene-set enrichment pathway analyses using the full GWAS data set identified exoDNase, NF-κB signaling and mitochondrial dysfunction as biological processes related to timing of menopause.

Identification of novel germline polymorphisms governing capecitabine sensitivity

BACKGROUND

Capecitabine, an oral 5-fluorouracil (5-FU) prodrug, is widely used in the treatment of breast, colorectal, and gastric cancers. To guide the selection of patients with potentially the greatest benefit of experiencing antitumor efficacy, or, alternatively, of developing toxicities, identifying genomic predictors of capecitabine sensitivity could permit its more informed use.

METHODS

The objective of this study was to perform capecitabine sensitivity genome-wide association studies (GWAS) using 503 well genotyped human cell lines from individuals representing multiple different world populations. A meta-analysis that included all ethnic populations then enabled the identification of novel germline determinants (single nucleotide polymorphisms [SNPs]) of capecitabine susceptibility.

RESULTS

First, an intrapopulation GWAS of Caucasian individuals identified reference SNP 4702484 (rs4702484) (within adenylate cyclase 2 [ADCY2]) at a level reaching genome-wide significance (P = 5.2 × 10(-8) ). This SNP is located upstream of the 5 methyltetrahydrofolate-homocysteine methyltransferase reductase (MTRR) gene, and it is known that the enzyme for MTRR is involved in the methionine-folate biosynthesis and metabolism pathway, which is the primary target of 5-FU-related compounds, although the authors were unable to identify a direct relation between rs4702484 and MTRR expression in a tested subset of cells. In the meta-analysis, 4 SNPs comprised the top hits, which, again, included rs4702484 and 3 additional SNPs (rs8101143, rs576523, and rs361433) that approached genome-wide significance (P values from 1.9 × 10(-7) to 8.8 × 10(-7) ). The meta-analysis also identified 1 missense variant (rs11722476; serine to asparagine) within switch/sucrose nonfermentable-related, matrix-associated, actin-dependent regulator of chromatin (SMARCAD1), a novel gene for association with capecitabine/5-FU susceptibility.

CONCLUSIONS

Toward the goal of individualizing cancer chemotherapy, the current study identified novel SNPs and genes associated with capecitabine sensitivity that are potentially informative and testable in any patient regardless of ethnicity.

Absence of functional LIN28B mutations in a large cohort of patients with idiopathic central precocious puberty

AIM

To investigate LIN28B gene variants in children with idiopathic central precocious puberty (CPP).

PATIENTS AND METHODS

We studied 178 Brazilian children with CPP (171 girls, 16.8% familial cases). A large multiethnic group (1,599 subjects; Multiethnic Cohort, MEC) was used as control. DNA analysis and biochemical in vitro studies were performed.

RESULTS

A heterozygous LIN28B variant, p.H199R, was identified in a girl who developed CPP at 5.2 years. This variant was absent in 310 Brazilian control individuals, but it was found in the same allele frequency in women from the MEC cohort, independent of the age of menarche. Functional studies revealed that when ectopically expressed in cells, the mutant protein was capable of binding pre-let-7 microRNA and inhibiting let-7 expression to the same extent as wild-type Lin28B protein. Other rare LIN28B variants (p.P173P, c.198+ 32_33delCT, g.9575731A>C and c.-11C>T) were identified in CPP patients and controls. Therefore, no functional mutation was identified.

CONCLUSION

In vitro studies revealed that the rare LIN28B p.H199R variant identified in a girl with CPP does not affect the Lin28B function in the regulation of let-7 expression. Although LIN28B SNPs were associated with normal pubertal timing, rare variations in this gene do not seem to be commonly involved in the molecular pathogenesis of CPP.

Chromosome 9 ALS and FTD locus is probably derived from a single founder

We and others have recently reported an association between amyotrophic lateral sclerosis (ALS) and single nucleotide polymorphisms on chromosome 9p21 in several populations. Here we show that the associated haplotype is the same in all populations and that several families previously shown to have genetic linkage to this region also share this haplotype. The most parsimonious explanation of these data are that there is a single founder for this form of disease.

Eight genetic loci associated with variation in lipoprotein-associated phospholipase A2 mass and activity and coronary heart disease: meta-analysis of genome-wide association studies from five community-based studies

AIMS

Lipoprotein-associated phospholipase A2 (Lp-PLA2) generates proinflammatory and proatherogenic compounds in the arterial vascular wall and is a potential therapeutic target in coronary heart disease (CHD). We searched for genetic loci related to Lp-PLA2 mass or activity by a genome-wide association study as part of the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium.

METHODS AND RESULTS

In meta-analyses of findings from five population-based studies, comprising 13 664 subjects, variants at two loci (PLA2G7, CETP) were associated with Lp-PLA2 mass. The strongest signal was at rs1805017 in PLA2G7 [P = 2.4 × 10(-23), log Lp-PLA2 difference per allele (beta): 0.043]. Variants at six loci were associated with Lp-PLA2 activity (PLA2G7, APOC1, CELSR2, LDL, ZNF259, SCARB1), among which the strongest signals were at rs4420638, near the APOE-APOC1-APOC4-APOC2 cluster [P = 4.9 × 10(-30); log Lp-PLA2 difference per allele (beta): -0.054]. There were no significant gene-environment interactions between these eight polymorphisms associated with Lp-PLA2 mass or activity and age, sex, body mass index, or smoking status. Four of the polymorphisms (in APOC1, CELSR2, SCARB1, ZNF259), but not PLA2G7, were significantly associated with CHD in a second study.

CONCLUSION

Levels of Lp-PLA2 mass and activity were associated with PLA2G7, the gene coding for this protein. Lipoprotein-associated phospholipase A2 activity was also strongly associated with genetic variants related to low-density lipoprotein cholesterol levels.

Genomewide association study using a high-density single nucleotide polymorphism array and case-control design identifies a novel essential hypertension susceptibility locus in the promoter region of endothelial NO synthase

Essential hypertension is a multifactorial disorder and is the main risk factor for renal and cardiovascular complications. The research on the genetics of hypertension has been frustrated by the small predictive value of the discovered genetic variants. The HYPERGENES Project investigated associations between genetic variants and essential hypertension pursuing a 2-stage study by recruiting cases and controls from extensively characterized cohorts recruited over many years in different European regions. The discovery phase consisted of 1865 cases and 1750 controls genotyped with 1M Illumina array. Best hits were followed up in a validation panel of 1385 cases and 1246 controls that were genotyped with a custom array of 14 055 markers. We identified a new hypertension susceptibility locus (rs3918226) in the promoter region of the endothelial NO synthase gene (odds ratio: 1.54 [95% CI: 1.37-1.73]; combined P=2.58 · 10(-13)). A meta-analysis, using other in silico/de novo genotyping data for a total of 21 714 subjects, resulted in an overall odds ratio of 1.34 (95% CI: 1.25-1.44; P=1.032 · 10(-14)). The quantitative analysis on a population-based sample revealed an effect size of 1.91 (95% CI: 0.16-3.66) for systolic and 1.40 (95% CI: 0.25-2.55) for diastolic blood pressure. We identified in silico a potential binding site for ETS transcription factors directly next to rs3918226, suggesting a potential modulation of endothelial NO synthase expression. Biological evidence links endothelial NO synthase with hypertension, because it is a critical mediator of cardiovascular homeostasis and blood pressure control via vascular tone regulation. This finding supports the hypothesis that there may be a causal genetic variation at this locus.

Genome-wide meta-analysis of psoriatic arthritis identifies susceptibility locus at REL

Psoriatic arthritis (PsA) is a chronic inflammatory musculoskeletal disease affecting up to 30% of psoriasis vulgaris (PsV) cases and approximately 0.25% to 1% of the general population. To identify common susceptibility loci, we performed a meta-analysis of three imputed genome-wide association studies (GWAS) on psoriasis, stratified for PsA. A total of 1,160,703 SNPs were analyzed in the discovery set consisting of 535 PsA cases and 3,432 controls from Germany, the United States and Canada. We followed up two SNPs in 1,931 PsA cases and 6,785 controls comprising six independent replication panels from Germany, Estonia, the United States and Canada. In the combined analysis, a genome-wide significant association was detected at 2p16 near the REL locus encoding c-Rel (rs13017599, P=1.18×10⁻⁸, OR=1.27, 95% CI=1.18–1.35). The rs13017599 polymorphism is known to associate with rheumatoid arthritis (RA), and another SNP near REL (rs702873) was recently implicated in PsV susceptibility. However, conditional analysis indicated that rs13017599, rather than rs702873, accounts for the PsA association at REL. We hypothesize that c-Rel, as a member of the Rel/NF-κB family, is associated with PsA in the context of disease pathways that involve other identified PsA and PsV susceptibility genes including TNIP1, TNFAIP3 and NFκBIA.

Genome-wide association of copy-number variation reveals an association between short stature and the presence of low-frequency genomic deletions

Height is a model polygenic trait that is highly heritable. Genome-wide association studies have identified hundreds of single-nucleotide polymorphisms associated with stature, but the role of structural variation in determining height is largely unknown. We performed a genome-wide association study of copy-number variation and stature in a clinical cohort of children who had undergone comparative genomic hybridization (CGH) microarray analysis for clinical indications. We found that subjects with short stature had a greater global burden of copy-number variants (CNVs) and a greater average CNV length than did controls (p < 0.002). These associations were present for lower-frequency (<5%) and rare (<1%) deletions, but there were no significant associations seen for duplications. Known gene-deletion syndromes did not account for our findings, and we saw no significant associations with tall stature. We then extended our findings into a population-based cohort and found that, in agreement with the clinical cohort study, an increased burden of lower-frequency deletions was associated with shorter stature (p = 0.015). Our results suggest that in individuals undergoing copy-number analysis for clinical indications, short stature increases the odds that a low-frequency deletion will be found. Additionally, copy-number variation might contribute to genetic variation in stature in the general population.

GWAtoolbox: an R package for fast quality control and handling of genome-wide association studies meta-analysis data

SUMMARY

The GWAtoolbox is an R package that standardizes and accelerates the handling of data from genome-wide association studies (GWAS), particularly in the context of large-scale GWAS meta-analyses. A key feature of GWAtoolbox is its ability to perform quality control (QC) of any number of files in a matter of minutes. The implemented workflow has been structured to check three particular data quality aspects: (i) data formatting, (ii) quality of the GWAS results and (iii) data consistency across studies. Output consists of an extensive list of quality statistics and plots which allow inspection of individual files and between-study comparison to identify systematic bias.

AVAILABILITY

http://www.eurac.edu/GWAtoolbox

CONTACT

cfuchsb@umich.edu; daniel.taliun@eurac.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Common variation at 3p22.1 and 7p15.3 influences multiple myeloma risk

To identify risk variants for multiple myeloma, we conducted a genome-wide association study of 1,675 individuals with multiple myeloma and 5,903 control subjects. We identified risk loci for multiple myeloma at 3p22.1 (rs1052501 in ULK4; odds ratio (OR) = 1.32; P = 7.47 × 10(-9)) and 7p15.3 (rs4487645, OR = 1.38; P = 3.33 × 10(-15)). In addition, we observed a promising association at 2p23.3 (rs6746082, OR = 1.29; P = 1.22 × 10(-7)). Our study identifies new genomic regions associated with multiple myeloma risk that may lead to new etiological insights.