Assessing the impact of population stratification on genetic association studies

Population stratification refers to differences in allele frequencies between cases and controls due to systematic differences in ancestry rather than association of genes with disease. It has been proposed that false positive associations due to stratification can be controlled by genotyping a few dozen unlinked genetic markers. To assess stratification empirically, we analyzed data from 11 case-control and case-cohort association studies. We did not detect statistically significant evidence for stratification but did observe that assessments based on a few dozen markers lack power to rule out moderate levels of stratification that could cause false positive associations in studies designed to detect modest genetic risk factors. After increasing the number of markers and samples in a case-cohort study (the design most immune to stratification), we found that stratification was in fact present. Our results suggest that modest amounts of stratification can exist even in well designed studies.     

Genome-wide association study of prostate cancer in men of African ancestry identifies a susceptibility locus at 17q21

In search of common risk alleles for prostate cancer that could contribute to high rates of the disease in men of African ancestry, we conducted a genome-wide association study, with 1,047,986 SNP markers examined in 3,425 African-Americans with prostate cancer (cases) and 3,290 African-American male controls. We followed up the most significant 17 new associations from stage 1 in 1,844 cases and 3,269 controls of African ancestry. We identified a new risk variant on chromosome 17q21 (rs7210100, odds ratio per allele = 1.51, P = 3.4 × 10(-13)). The frequency of the risk allele is ～5% in men of African descent, whereas it is rare in other populations (<1%). Further studies are needed to investigate the biological contribution of this allele to prostate cancer risk. These findings emphasize the importance of conducting genome-wide association studies in diverse populations.     

Demonstrating stratification in a European American population

Population stratification occurs in case-control association studies when allele frequencies differ between cases and controls because of ancestry. Stratification may lead to false positive associations, although this issue remains controversial. Empirical studies have found little evidence of stratification in European-derived populations, but potentially significant levels of stratification could not be ruled out. We studied a European American panel discordant for height, a heritable trait that varies widely across Europe. Genotyping 178 SNPs and applying standard analytical methods yielded no evidence of stratification. But a SNP in the gene LCT that varies widely in frequency across Europe was strongly associated with height (P < 10(-6)). This apparent association was largely or completely due to stratification; rematching individuals on the basis of European ancestry greatly reduced the apparent association, and no association was observed in Polish or Scandinavian individuals. The failure of standard methods to detect this stratification indicates that new methods may be required.     

'Racial' differences in genetic effects for complex diseases

'Racial' differences are frequently debated in clinical, epidemiological and molecular research and beyond. In particular, there is considerable controversy regarding the existence and importance of 'racial' differences in genetic effects for complex diseases influenced by a large number of genes. An important question is whether ancestry influences the impact of each gene variant on the disease risk. Here, we addressed this question by examining the genetic effects for 43 validated gene-disease associations across 697 study populations of various descents. The frequencies of the genetic marker of interest in the control populations often (58%) showed large heterogeneity (statistical variability) between 'races'. Conversely, we saw large heterogeneity in the genetic effects (odds ratios) between 'races' in only 14% of cases. Genetic markers for proposed gene-disease associations vary in frequency across populations, but their biological impact on the risk for common diseases may usually be consistent across traditional 'racial' boundaries.     

Meta-analysis of genome-wide association studies of asthma in ethnically diverse North American populations

Asthma is a common disease with a complex risk architecture including both genetic and environmental factors. We performed a meta-analysis of North American genome-wide association studies of asthma in 5,416 individuals with asthma (cases) including individuals of European American, African American or African Caribbean, and Latino ancestry, with replication in an additional 12,649 individuals from the same ethnic groups. We identified five susceptibility loci. Four were at previously reported loci on 17q21, near IL1RL1, TSLP and IL33, but we report for the first time, to our knowledge, that these loci are associated with asthma risk in three ethnic groups. In addition, we identified a new asthma susceptibility locus at PYHIN1, with the association being specific to individuals of African descent (P = 3.9 × 10(-9)). These results suggest that some asthma susceptibility loci are robust to differences in ancestry when sufficiently large samples sizes are investigated, and that ancestry-specific associations also contribute to the complex genetic architecture of asthma.     

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

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

An Icelandic example of the impact of population structure on association studies

The impact of population structure on association studies undertaken to identify genetic variants underlying common human diseases is an issue of growing interest. Spurious associations of alleles with disease phenotypes may be obtained or true associations overlooked when allele frequencies differ notably among subpopulations that are not represented equally among cases and controls. Population structure influences even carefully designed studies and can affect the validity of association results. Most study designs address this problem by sampling cases and controls from groups that share the same nationality or self-reported ethnic background, with the implicit assumption that no substructure exists within such groups. We examined population structure in the Icelandic gene pool using extensive genealogical and genetic data. Our results indicate that sampling strategies need to take account of substructure even in a relatively homogenous genetic isolate. This will probably be even more important in larger populations.     

A model-based approach for analysis of spatial structure in genetic data

Characterizing genetic diversity within and between populations has broad applications in studies of human disease and evolution. We propose a new approach, spatial ancestry analysis, for the modeling of genotypes in two- or three-dimensional space. In spatial ancestry analysis (SPA), we explicitly model the spatial distribution of each SNP by assigning an allele frequency as a continuous function in geographic space. We show that the explicit modeling of the allele frequency allows individuals to be localized on the map on the basis of their genetic information alone. We apply our SPA method to a European and a worldwide population genetic variation data set and identify SNPs showing large gradients in allele frequency, and we suggest these as candidate regions under selection. These regions include SNPs in the well-characterized LCT region, as well as at loci including FOXP2, OCA2 and LRP1B.     

MYH9 is associated with nondiabetic end-stage renal disease in African Americans

As end-stage renal disease (ESRD) has a four times higher incidence in African Americans compared to European Americans, we hypothesized that susceptibility alleles for ESRD have a higher frequency in the West African than the European gene pool. We carried out a genome-wide admixture scan in 1,372 ESRD cases and 806 controls and found a highly significant association between excess African ancestry and nondiabetic ESRD (lod score = 5.70) but not diabetic ESRD (lod = 0.47) on chromosome 22q12. Each copy of the European ancestral allele conferred a relative risk of 0.50 (95% CI = 0.39-0.63) compared to African ancestry. Multiple common SNPs (allele frequencies ranging from 0.2 to 0.6) in the gene encoding nonmuscle myosin heavy chain type II isoform A (MYH9) were associated with two to four times greater risk of nondiabetic ESRD and accounted for a large proportion of the excess risk of ESRD observed in African compared to European Americans.     

Genome-wide association study in individuals of South Asian ancestry identifies six new type 2 diabetes susceptibility loci

We carried out a genome-wide association study of type-2 diabetes (T2D) in individuals of South Asian ancestry. Our discovery set included 5,561 individuals with T2D (cases) and 14,458 controls drawn from studies in London, Pakistan and Singapore. We identified 20 independent SNPs associated with T2D at P < 10(-4) for testing in a replication sample of 13,170 cases and 25,398 controls, also all of South Asian ancestry. In the combined analysis, we identified common genetic variants at six loci (GRB14, ST6GAL1, VPS26A, HMG20A, AP3S2 and HNF4A) newly associated with T2D (P = 4.1 × 10(-8) to P = 1.9 × 10(-11)). SNPs at GRB14 were also associated with insulin sensitivity (P = 5.0 × 10(-4)), and SNPs at ST6GAL1 and HNF4A were also associated with pancreatic beta-cell function (P = 0.02 and P = 0.001, respectively). Our findings provide additional insight into mechanisms underlying T2D and show the potential for new discovery from genetic association studies in South Asians, a population with increased susceptibility to T2D.     

Genome-wide copy number variation study associates metabotropic glutamate receptor gene networks with attention deficit hyperactivity disorder

Attention deficit hyperactivity disorder (ADHD) is a common, heritable neuropsychiatric disorder of unknown etiology. We performed a whole-genome copy number variation (CNV) study on 1,013 cases with ADHD and 4,105 healthy children of European ancestry using 550,000 SNPs. We evaluated statistically significant findings in multiple independent cohorts, with a total of 2,493 cases with ADHD and 9,222 controls of European ancestry, using matched platforms. CNVs affecting metabotropic glutamate receptor genes were enriched across all cohorts (P = 2.1 × 10(-9)). We saw GRM5 (encoding glutamate receptor, metabotropic 5) deletions in ten cases and one control (P = 1.36 × 10(-6)). We saw GRM7 deletions in six cases, and we saw GRM8 deletions in eight cases and no controls. GRM1 was duplicated in eight cases. We experimentally validated the observed variants using quantitative RT-PCR. A gene network analysis showed that genes interacting with the genes in the GRM family are enriched for CNVs in ～10% of the cases (P = 4.38 × 10(-10)) after correction for occurrence in the controls. We identified rare recurrent CNVs affecting glutamatergic neurotransmission genes that were overrepresented in multiple ADHD cohorts.     

A common variant at the TERT-CLPTM1L locus is associated with estrogen receptor-negative breast cancer

Estrogen receptor (ER)-negative breast cancer shows a higher incidence in women of African ancestry compared to women of European ancestry. In search of common risk alleles for ER-negative breast cancer, we combined genome-wide association study (GWAS) data from women of African ancestry (1,004 ER-negative cases and 2,745 controls) and European ancestry (1,718 ER-negative cases and 3,670 controls), with replication testing conducted in an additional 2,292 ER-negative cases and 16,901 controls of European ancestry. We identified a common risk variant for ER-negative breast cancer at the TERT-CLPTM1L locus on chromosome 5p15 (rs10069690: per-allele odds ratio (OR) = 1.18 per allele, P = 1.0 × 10(-10)). The variant was also significantly associated with triple-negative (ER-negative, progesterone receptor (PR)-negative and human epidermal growth factor-2 (HER2)-negative) breast cancer (OR = 1.25, P = 1.1 × 10(-9)), particularly in younger women (<50 years of age) (OR = 1.48, P = 1.9 × 10(-9)). Our results identify a genetic locus associated with estrogen receptor negative breast cancer subtypes in multiple populations.     

Genetic variation in DPP6 is associated with susceptibility to amyotrophic lateral sclerosis

We identified a SNP in the DPP6 gene that is consistently strongly associated with susceptibility to amyotrophic lateral sclerosis (ALS) in different populations of European ancestry, with an overall P value of 5.04 x 10(-8) in 1,767 cases and 1,916 healthy controls and with an odds ratio of 1.30 (95% confidence interval (CI) of 1.18-1.43). Our finding is the first report of a genome-wide significant association with sporadic ALS and may be a target for future functional studies.     

A genome-wide association meta-analysis identifies new childhood obesity loci

Multiple genetic variants have been associated with adult obesity and a few with severe obesity in childhood; however, less progress has been made in establishing genetic influences on common early-onset obesity. We performed a North American, Australian and European collaborative meta-analysis of 14 studies consisting of 5,530 cases (≥95th percentile of body mass index (BMI)) and 8,318 controls (<50th percentile of BMI) of European ancestry. Taking forward the eight newly discovered signals yielding association with P < 5 × 10(-6) in nine independent data sets (2,818 cases and 4,083 controls), we observed two loci that yielded genome-wide significant combined P values near OLFM4 at 13q14 (rs9568856; P = 1.82 × 10(-9); odds ratio (OR) = 1.22) and within HOXB5 at 17q21 (rs9299; P = 3.54 × 10(-9); OR = 1.14). Both loci continued to show association when two extreme childhood obesity cohorts were included (2,214 cases and 2,674 controls). These two loci also yielded directionally consistent associations in a previous meta-analysis of adult BMI(1).     

A unified mixed-model method for association mapping that accounts for multiple levels of relatedness

As population structure can result in spurious associations, it has constrained the use of association studies in human and plant genetics. Association mapping, however, holds great promise if true signals of functional association can be separated from the vast number of false signals generated by population structure. We have developed a unified mixed-model approach to account for multiple levels of relatedness simultaneously as detected by random genetic markers. We applied this new approach to two samples: a family-based sample of 14 human families, for quantitative gene expression dissection, and a sample of 277 diverse maize inbred lines with complex familial relationships and population structure, for quantitative trait dissection. Our method demonstrates improved control of both type I and type II error rates over other methods. As this new method crosses the boundary between family-based and structured association samples, it provides a powerful complement to currently available methods for association mapping.     

Genome-wide association study identifies susceptibility loci for IgA nephropathy

We carried out a genome-wide association study of IgA nephropathy, a major cause of kidney failure worldwide. We studied 1,194 cases and 902 controls of Chinese Han ancestry, with targeted follow up in Chinese and European cohorts comprising 1,950 cases and 1,920 controls. We identified three independent loci in the major histocompatibility complex, as well as a common deletion of CFHR1 and CFHR3 at chromosome 1q32 and a locus at chromosome 22q12 that each surpassed genome-wide significance (P values for association between 1.59 × 10?2? and 4.84 × 10?? and minor allele odds ratios of 0.63-0.80). These five loci explain 4-7% of the disease variance and up to a tenfold variation in interindividual risk. Many of the alleles that protect against IgA nephropathy impart increased risk for other autoimmune or infectious diseases, and IgA nephropathy risk allele frequencies closely parallel the variation in disease prevalence among Asian, European and African populations, suggesting complex selective pressures.     

Genome-wide association study identifies five loci associated with susceptibility to pancreatic cancer in Chinese populations

Pancreatic cancer has the lowest survival rate among human cancers, and there are no effective markers for its screening and early diagnosis. To identify genetic susceptibility markers for this cancer, we carried out a genome-wide association study on 981 individuals with pancreatic cancer (cases) and 1,991 cancer-free controls of Chinese descent using 666,141 autosomal SNPs. Promising associations were replicated in an additional 2,603 pancreatic cancer cases and 2,877 controls recruited from 25 hospitals in 16 provinces or cities in China. We identified five new susceptibility loci at chromosomes 21q21.3, 5p13.1, 21q22.3, 22q13.32 and 10q26.11 (P = 2.24 × 10(-13) to P = 4.18 × 10(-10)) in addition to 13q22.1 previously reported in populations of European ancestry. These results advance our understanding of the development of pancreatic cancer and highlight potential targets for the prevention or treatment of this cancer.     

A genome-wide association study identifies three new risk loci for Kawasaki disease

We performed a genome-wide association study (GWAS) of Kawasaki disease in Japanese subjects using data from 428 individuals with Kawasaki disease (cases) and 3,379 controls genotyped at 473,803 SNPs. We validated the association results in two independent replication panels totaling 754 cases and 947 controls. We observed significant associations in the FAM167A-BLK region at 8p22-23 (rs2254546, P = 8.2 × 10(-21)), in the human leukocyte antigen (HLA) region at 6p21.3 (rs2857151, P = 4.6 × 10(-11)) and in the CD40 region at 20q13 (rs4813003, P = 4.8 × 10(-8)). We also replicated the association of a functional SNP of FCGR2A (rs1801274, P = 1.6 × 10(-6)) identified in a recently reported GWAS of Kawasaki disease. Our findings provide new insights into the pathogenesis and pathophysiology of Kawasaki disease.     

Meta-analysis identifies multiple loci associated with kidney function-related traits in east Asian populations

Chronic kidney disease (CKD), impairment of kidney function, is a serious public health problem, and the assessment of genetic factors influencing kidney function has substantial clinical relevance. Here, we report a meta-analysis of genome-wide association studies for kidney function-related traits, including 71,149 east Asian individuals from 18 studies in 11 population-, hospital- or family-based cohorts, conducted as part of the Asian Genetic Epidemiology Network (AGEN). Our meta-analysis identified 17 loci newly associated with kidney function-related traits, including the concentrations of blood urea nitrogen, uric acid and serum creatinine and estimated glomerular filtration rate based on serum creatinine levels (eGFRcrea) (P < 5.0 × 10(-8)). We further examined these loci with in silico replication in individuals of European ancestry from the KidneyGen, CKDGen and GUGC consortia, including a combined total of ～110,347 individuals. We identify pleiotropic associations among these loci with kidney function-related traits and risk of CKD. These findings provide new insights into the genetics of kidney function.     

Meta-analysis identifies nine new loci associated with rheumatoid arthritis in the Japanese population

Rheumatoid arthritis is a common autoimmune disease characterized by chronic inflammation. We report a meta-analysis of genome-wide association studies (GWAS) in a Japanese population including 4,074 individuals with rheumatoid arthritis (cases) and 16,891 controls, followed by a replication in 5,277 rheumatoid arthritis cases and 21,684 controls. Our study identified nine loci newly associated with rheumatoid arthritis at a threshold of P < 5.0 × 10(-8), including B3GNT2, ANXA3, CSF2, CD83, NFKBIE, ARID5B, PDE2A-ARAP1, PLD4 and PTPN2. ANXA3 was also associated with susceptibility to systemic lupus erythematosus (P = 0.0040), and B3GNT2 and ARID5B were associated with Graves' disease (P = 3.5 × 10(-4) and 2.9 × 10(-4), respectively). We conducted a multi-ancestry comparative analysis with a previous meta-analysis in individuals of European descent (5,539 rheumatoid arthritis cases and 20,169 controls). This provided evidence of shared genetic risks of rheumatoid arthritis between the populations.