Late-onset Alzheimer's disease (LOAD) is the most common form of dementia in the elderly. The National Institute of Aging-Late Onset Alzheimer's Disease Family Study and the National Cell Repository for Alzheimer's Disease conducted a joint genome-wide association study (GWAS) of multiplex LOAD families (3,839 affected and unaffected individuals from 992 families plus additional unrelated neurologically evaluated normal subjects) using the 610 IlluminaQuad panel. This cohort represents the largest family-based GWAS of LOAD to date, with analyses limited here to the European-American subjects. SNPs near APOE gave highly significant results (e.g., rs2075650, p = 3.2×10(-81)), but no other genome-wide significant evidence for association was obtained in the full sample. Analyses that stratified on APOE genotypes identified SNPs on chromosome 10p14 in CUGBP2 with genome-wide significant evidence for association within APOE ε4 homozygotes (e.g., rs201119, p = 1.5×10(-8)). Association in this gene was replicated in an independent sample consisting of three cohorts. There was evidence of association for recently-reported LOAD risk loci, including BIN1 (rs7561528, p = 0.009 with, and p = 0.03 without, APOE adjustment) and CLU (rs11136000, p = 0.023 with, and p = 0.008 without, APOE adjustment), with weaker support for CR1. However, our results provide strong evidence that association with PICALM (rs3851179, p = 0.69 with, and p = 0.039 without, APOE adjustment) and EXOC3L2 is affected by correlation with APOE, and thus may represent spurious association. Our results indicate that genetic structure coupled with ascertainment bias resulting from the strong APOE association affect genome-wide results and interpretation of some recently reported associations. We show that a locus such as APOE, with large effects and strong association with disease, can lead to samples that require appropriate adjustment for this locus to avoid both false positive and false negative evidence of association. We suggest that similar adjustments may also be needed for many other large multi-site studies.
The PR interval on the electrocardiogram reflects atrial and atrioventricular nodal conduction time. The PR interval is heritable, provides important information about arrhythmia risk, and has been suggested to differ among human races. Genome-wide association (GWA) studies have identified common genetic determinants of the PR interval in individuals of European and Asian ancestry, but there is a general paucity of GWA studies in individuals of African ancestry. We performed GWA studies in African American individuals from four cohorts (n = 6,247) to identify genetic variants associated with PR interval duration. Genotyping was performed using the Affymetrix 6.0 microarray. Imputation was performed for 2.8 million single nucleotide polymorphisms (SNPs) using combined YRI and CEU HapMap phase II panels. We observed a strong signal (rs3922844) within the gene encoding the cardiac sodium channel (SCN5A) with genome-wide significant association (p<2.5 x 10⁻⁸) in two of the four cohorts and in the meta-analysis. The signal explained 2% of PR interval variability in African Americans (beta = 5.1 msec per minor allele, 95% CI = 4.1-6.1, p = 3 x 10⁻²³). This SNP was also associated with PR interval (beta = 2.4 msec per minor allele, 95% CI = 1.8-3.0, p = 3 x 10⁻¹⁶) in individuals of European ancestry (n = 14,042), but with a smaller effect size (p for heterogeneity <0.001) and variability explained (0.5%). Further meta-analysis of the four cohorts identified genome-wide significant associations with SNPs in SCN10A (rs6798015), MEIS1 (rs10865355), and TBX5 (rs7312625) that were highly correlated with SNPs identified in European and Asian GWA studies. African ancestry was associated with increased PR duration (13.3 msec, p = 0.009) in one but not the other three cohorts. Our findings demonstrate the relevance of common variants to African Americans at four loci previously associated with PR interval in European and Asian samples and identify an association signal at one of these loci that is more strongly associated with PR interval in African Americans than in Europeans.
Damage to the optic nerve (e.g. from glaucoma) has an adverse and often irreversible impact on vision. Earlier studies have suggested that the size of the optic nerve head could be governed by hereditary factors. We conducted a genome-wide association study (GWAS) on 4445 Singaporean individuals (n = 2132 of Indian and n = 2313 of Malay ancestry, respectively), with replication in Rotterdam, the Netherlands (n = 9326 individuals of Caucasian ancestry) using the most widely reported parameter for optic disc traits, the optic disc area. We identified a novel locus on chromosome 22q13.1, CARD10, which strongly associates with optic disc area in both Singaporean cohorts as well as in the Rotterdam Study (RS; rs9607469, per-allele change in optic disc area = 0.051 mm(2); P(meta) = 2.73×10(-12)) and confirmed the association between CDC7/TGFBR3 (lead single nucleotide polymorphism (SNP) rs1192415, P(meta) = 7.57×10(-17)) and ATOH7 (lead SNP rs7916697, P(meta) = 2.00 × 10(-15)) and optic disc area in Asians. This is the first Asian-based GWAS on optic disc area, identifying a novel locus for the optic disc area, but also confirming the results found in Caucasian persons suggesting that there are general genetic determinants applicable to the size of the optic disc across different ethnicities.
Genome-wide association studies and candidate gene studies in ulcerative colitis have identified 18 susceptibility loci. We conducted a meta-analysis of six ulcerative colitis genome-wide association study datasets, comprising 6,687 cases and 19,718 controls, and followed up the top association signals in 9,628 cases and 12,917 controls. We identified 29 additional risk loci (P < 5 × 10(-8)), increasing the number of ulcerative colitis-associated loci to 47. After annotating associated regions using GRAIL, expression quantitative trait loci data and correlations with non-synonymous SNPs, we identified many candidate genes that provide potentially important insights into disease pathogenesis, including IL1R2, IL8RA-IL8RB, IL7R, IL12B, DAP, PRDM1, JAK2, IRF5, GNA12 and LSP1. The total number of confirmed inflammatory bowel disease risk loci is now 99, including a minimum of 28 shared association signals between Crohn's disease and ulcerative colitis.
A three-stage approach was undertaken using genome-wide, case-control, and case-only association studies to identify genetic variants associated with heart failure mortality. In an Amish founder population (n = 851), cardiac hypertrophy, a trait integral to the adaptive response to failure, was found to be heritable (h² = 0.28, p = 0.0002) and GWAS revealed 21 candidate hypertrophy SNPs. In a case (n = 1,610)-control (n = 463) study in unrelated Caucasians, one of the SNPs associated with hypertrophy (rs2207418, p = 8 × 10⁻⁶), was associated with heart failure, RR = 1.85(1.25-2.73, p = 0.0019). In heart failure cases rs2207418 was associated with increased mortality, HR = 1.51(1.20-1.97, p = 0.0004). There was consistency between studies, with the GG allele being associated with increased ventricular mass (~13 g/m²) in the Amish, heart failure risk, and heart failure mortality. This SNP is in a gene desert of chromosome 20p12. Five genes are within 2.0 mbp of rs2207418 but with low LD between their SNPs and rs2207418. A region near this SNP is highly conserved in multiple vertebrates (lod score = 1,208). This conservation and the internal consistency across studies suggests that this region has biologic importance in heart failure, potentially acting as an enhancer or repressor element. rs2207418 may be useful for predicting a more progressive form of heart failure that may require aggressive therapy.
C-reactive protein (CRP) is a heritable marker of chronic inflammation that is strongly associated with cardiovascular disease. We sought to identify genetic variants that are associated with CRP levels.
Genome-wide association studies (GWAS) have identified several susceptibility loci for bipolar disorder (BP), most notably ANK3. However, most of the inherited risk for BP remains unexplained. One reason for the limited success may be the genetic heterogeneity of BP. Clinical sub-phenotypes of BP may identify more etiologically homogeneous subsets of patients, which can be studied with increased power to detect genetic variation. Here, we report on a mega-analysis of two widely studied sub-phenotypes of BP, age at onset and psychotic symptoms, which are familial and clinically significant. We combined data from three GWAS: NIMH Bipolar Disorder Genetic Association Information Network (GAIN-BP), NIMH Bipolar Disorder Genome Study (BiGS), and a German sample. The combined sample consisted of 2,836 BP cases with information on sub-phenotypes and 2,744 controls. Imputation was performed, resulting in 2.3 million SNPs available for analysis. No SNP reached genome-wide significance for either sub-phenotype. In addition, no SNP reached genome-wide significance in a meta-analysis with an independent replication sample. We had 80% power to detect associations with a common SNP at an OR of 1.6 for psychotic symptoms and a mean difference of 1.8 years in age at onset. Age at onset and psychotic symptoms in BP may be influenced by many genes of smaller effect sizes or other variants not measured well by SNP arrays, such as rare alleles.
Diabetic retinopathy (DR) is a microvascular complication of diabetes with a complex multifactorial pathogenesis. The aim of this study was to identify the susceptibility genes that increase the risk of DR in type 2 diabetes (T2D) and to further elucidate the underlying mechanism of DR pathogenesis.
High-percent mammographic density adjusted for age and body mass index is one of the strongest risk factors for breast cancer. We conducted a meta analysis of five genome-wide association studies of percent mammographic density and report an association with rs10995190 in ZNF365 (combined P = 9.6 × 10(-10)). Common variants in ZNF365 have also recently been associated with susceptibility to breast cancer.
Uric acid is the primary byproduct of purine metabolism. Hyperuricemia is associated with body mass index (BMI), sex, and multiple complex diseases including gout, hypertension (HTN), renal disease, and type 2 diabetes (T2D). Multiple genome-wide association studies (GWAS) in individuals of European ancestry (EA) have reported associations between serum uric acid levels (SUAL) and specific genomic loci. The purposes of this study were: 1) to replicate major signals reported in EA populations; and 2) to use the weak LD pattern in African ancestry population to better localize (fine-map) reported loci and 3) to explore the identification of novel findings cognizant of the moderate sample size.
Kawasaki disease (KD) is an acute systemic vasculitis syndrome that primarily affects infants and young children. Its etiology is unknown; however, epidemiological findings suggest that genetic predisposition underlies disease susceptibility. Taiwan has the third-highest incidence of KD in the world, after Japan and Korea. To investigate novel mechanisms that might predispose individuals to KD, we conducted a genome-wide association study (GWAS) in 250 KD patients and 446 controls in a Han Chinese population residing in Taiwan, and further validated our findings in an independent Han Chinese cohort of 208 cases and 366 controls. The most strongly associated single-nucleotide polymorphisms (SNPs) detected in the joint analysis corresponded to three novel loci. Among these KD-associated SNPs three were close to the COPB2 (coatomer protein complex beta-2 subunit) gene: rs1873668 (p = 9.52×10⁻⁵), rs4243399 (p = 9.93×10⁻⁵), and rs16849083 (p = 9.93×10⁻⁵). We also identified a SNP in the intronic region of the ERAP1 (endoplasmic reticulum amino peptidase 1) gene (rs149481, p(best) = 4.61×10⁻⁵). Six SNPs (rs17113284, rs8005468, rs10129255, rs2007467, rs10150241, and rs12590667) clustered in an area containing immunoglobulin heavy chain variable regions genes, with p(best)-values between 2.08×10⁻⁵ and 8.93×10⁻⁶, were also identified. This is the first KD GWAS performed in a Han Chinese population. The novel KD candidates we identified have been implicated in T cell receptor signaling, regulation of proinflammatory cytokines, as well as antibody-mediated immune responses. These findings may lead to a better understanding of the underlying molecular pathogenesis of KD.
Crohn's disease (CD) and celiac disease (CelD) are chronic intestinal inflammatory diseases, involving genetic and environmental factors in their pathogenesis. The two diseases can co-occur within families, and studies suggest that CelD patients have a higher risk to develop CD than the general population. These observations suggest that CD and CelD may share common genetic risk loci. Two such shared loci, IL18RAP and PTPN2, have already been identified independently in these two diseases. The aim of our study was to explicitly identify shared risk loci for these diseases by combining results from genome-wide association study (GWAS) datasets of CD and CelD. Specifically, GWAS results from CelD (768 cases, 1,422 controls) and CD (3,230 cases, 4,829 controls) were combined in a meta-analysis. Nine independent regions had nominal association p-value <1.0 x 10⁻⁵ in this meta-analysis and showed evidence of association to the individual diseases in the original scans (p-value < 1 x 10⁻² in CelD and < 1 x 10⁻³ in CD). These include the two previously reported shared loci, IL18RAP and PTPN2, with p-values of 3.37 x 10⁻⁸ and 6.39 x 10⁻⁹, respectively, in the meta-analysis. The other seven had not been reported as shared loci and thus were tested in additional CelD (3,149 cases and 4,714 controls) and CD (1,835 cases and 1,669 controls) cohorts. Two of these loci, TAGAP and PUS10, showed significant evidence of replication (Bonferroni corrected p-values <0.0071) in the combined CelD and CD replication cohorts and were firmly established as shared risk loci of genome-wide significance, with overall combined p-values of 1.55 x 10⁻¹⁰ and 1.38 x 10⁻¹¹ respectively. Through a meta-analysis of GWAS data from CD and CelD, we have identified four shared risk loci: PTPN2, IL18RAP, TAGAP, and PUS10. The combined analysis of the two datasets provided the power, lacking in the individual GWAS for single diseases, to detect shared loci with a relatively small effect.
Aortic root diameter is a clinically relevant trait due to its known relationship with the pathogenesis of aortic regurgitation and risk for aortic dissection. African Americans are an understudied population despite a particularly high burden of cardiovascular diseases. We report a genome-wide association study on aortic root diameter among African Americans enrolled in the HyperGEN study. We invoked a two-stage, mixed model procedure to jointly identify SNP allele and copy number variation effects.
High blood concentration of the N-terminal cleavage product of the B-type natriuretic peptide (NT-proBNP) is strongly associated with cardiac dysfunction and is increasingly used for heart failure diagnosis. To identify genetic variants associated with NT-proBNP level, we performed a genome-wide association analysis in 1325 individuals from South Tyrol, Italy, and followed up the most significant results in 1746 individuals from two German population-based studies. A genome-wide significant signal in the MTHFR-CLCN6-NPPA-NPPB gene cluster was replicated, after correction for multiple testing (replication one-sided P-value = 8.4 × 10(-10)). A conditional regression analysis of 128 single-nucleotide polymorphisms in the region of interest identified novel variants in the CLCN6 gene as independently associated with NT-proBNP. In this locus, four haplotypes were associated with increased NT-proBNP levels (haplotype-specific combined P-values from 8.3 × 10(-03) to 9.3 × 10(-11)). The observed increase in the NT-proBNP level was proportional to the number of haplotype copies present (i.e. dosage effect), with an increase associated with two copies that varied between 20 and 100 pg/ml across populations. The identification of novel variants in the MTHFR-CLCN6-NPPA-NPPB cluster provides new insights into the biological mechanisms of cardiac dysfunction.
Genome-wide association studies have identified several common genetic variants associated with breast cancer risk. It is likely, however, that a substantial proportion of such loci have not yet been discovered.
Multiple sclerosis (MS) is the most common cause of chronic neurologic disability beginning in early to middle adult life. Results from recent genome-wide association studies (GWAS) have substantially lengthened the list of disease loci and provide convincing evidence supporting a multifactorial and polygenic model of inheritance. Nevertheless, the knowledge of MS genetics remains incomplete, with many risk alleles still to be revealed.
Left ventricular (LV) mass and related phenotypes are heritable, important predictors of cardiovascular disease, particularly in hypertensive individuals.
Left ventricular mass, intraventricular septal thickness
Insulin-like growth factor-I (IGF-I) and insulin-like growth factor-binding protein-3 (IGFBP-3) are involved in cell replication, proliferation, differentiation, protein synthesis, carbohydrate homeostasis and bone metabolism. Circulating IGF-I and IGFBP-3 concentrations predict anthropometric traits and risk of cancer and cardiovascular disease. In a genome-wide association study of 10 280 middle-aged and older men and women from four community-based cohort studies, we confirmed a known association of single nucleotide polymorphisms in the IGFBP3 gene region on chromosome 7p12.3 with IGFBP-3 concentrations using a significance threshold of P < 5 × 10(-8) (P = 3.3 × 10(-101)). Furthermore, the same IGFBP3 gene locus (e.g. rs11977526) that was associated with IGFBP-3 concentrations was also associated with the opposite direction of effect, with IGF-I concentration after adjustment for IGFBP-3 concentration (P = 1.9 × 10(-26)). A novel and independent locus on chromosome 7p12.3 (rs700752) had genome-wide significant associations with higher IGFBP-3 (P = 4.4 × 10(-21)) and higher IGF-I (P = 4.9 × 10(-9)) concentrations; when the two measurements were adjusted for one another, the IGF-I association was attenuated but the IGFBP-3 association was not. Two additional loci demonstrated genome-wide significant associations with IGFBP-3 concentration (rs1065656, chromosome 16p13.3, P = 1.2 × 10(-11), IGFALS, a confirmatory finding; and rs4234798, chromosome 4p16.1, P = 4.5 × 10(-10), SORCS2, a novel finding). Together, the four genome-wide significant loci explained 6.5% of the population variation in IGFBP-3 concentration. Furthermore, we observed a borderline statistically significant association between IGF-I concentration and FOXO3 (rs2153960, chromosome 6q21, P = 5.1 × 10(-7)), a locus associated with longevity. These genetic loci deserve further investigation to elucidate the biological basis for the observed associations and clarify their possible role in IGF-mediated regulation of cell growth and metabolism.
Type I interferon (IFN) is used for the treatment of chronic hepatitis C virus (HCV) infection. Despite advances in antiviral therapy, a large proportion of patients remain infected following current therapies. Through a genome-wide scan, we found two variants (rs8099917 and rs12979860) in the IL-28B locus that affect the outcome of PEG-IFN and ribavirin combination therapy, consistent with recent studies (P = 6.52×10(-8); odds ratio 2.46 and P = 8.63×10(-8), odds ratio 2.40, respectively). Significant associations were also observed in the case of IFN monotherapy for HCV genotypes 1b and 2a. With rs8099917, HCV genotype 1b patients had a significantly lower frequency of the favourable genotype (86.6 %) compared with healthy controls (91.7 %), and HCV genotype 2a patients had an intermediate frequency (89.9 %). Similar results were found for rs12979860. Fine-mapping analysis revealed that rs8099917 had the strongest association with treatment outcome and 14 others, including four novel single nucleotide polymorphisms, had comparable associations. Haplotype analysis revealed that none of the haplotypes showed stronger association than any single marker. Early non-responders who could not achieve 2 log viral decline during the first 12 weeks of treatment had higher odds ratios for these two variants. The favourable allele of rs8099917 is also associated with initial viral decline at 2 and 4 weeks following the start of therapy. Multivariate analysis of PEG-IFN and ribavirin-treated patients showed that rs8099917 genotype, viral load, fibrosis and age were significant predictors of response to therapy. Common variation at the IL-28B locus is predictive of various IFN-based therapies for HCV independent of regimen or HCV genotype.
GWAS of Alzheimer's disease (late onset)
