Candidate lung tumor susceptibility genes identified through whole-genome association analyses in inbred mice

We performed a whole-genome association analysis of lung tumor susceptibility using dense SNP maps ( approximately 1 SNP per 20 kb) in inbred mice. We reproduced the pulmonary adenoma susceptibility 1 (Pas1) locus identified in previous linkage studies and further narrowed this quantitative trait locus (QTL) to a region of less than 0.5 Mb in which at least two genes, Kras2 (Kirsten rat sarcoma oncogene 2) and Casc1 (cancer susceptibility candidate 1; also known as Las1), are strong candidates. Casc1 knockout mouse tumor bioassays showed that Casc1-deficient mice were susceptible to chemical induction of lung tumors. We also found three more genetic loci for lung adenoma development. Analysis of one of these candidate loci identified a previously uncharacterized gene Lasc1, bearing a nonsynonymous substitution (D102E). We found that the Lasc1 Glu102 allele preferentially promotes lung tumor cell growth. Our findings demonstrate the prospects for using dense SNP maps in laboratory mice to refine previous QTL regions and identify genetic determinants of complex traits.     

Evaluating and improving power in whole-genome association studies using fixed marker sets

Emerging technologies make it possible for the first time to genotype hundreds of thousands of SNPs simultaneously, enabling whole-genome association studies. Using empirical genotype data from the International HapMap Project, we evaluate the extent to which the sets of SNPs contained on three whole-genome genotyping arrays capture common SNPs across the genome, and we find that the majority of common SNPs are well captured by these products either directly or through linkage disequilibrium. We explore analytical strategies that use HapMap data to improve power of association studies conducted with these fixed sets of markers and show that limited inclusion of specific haplotype tests in association analysis can increase the fraction of common variants captured by 25-100%. Finally, we introduce a Bayesian approach to association analysis by weighting the likelihood of each statistical test to reflect the number of putative causal alleles to which it is correlated.     

Transferability of tag SNPs in genetic association studies in multiple populations

A general question for linkage disequilibrium-based association studies is how power to detect an association is compromised when tag SNPs are chosen from data in one population sample and then deployed in another sample. Specifically, it is important to know how well tags picked from the HapMap DNA samples capture the variation in other samples. To address this, we collected dense data uniformly across the four HapMap population samples and eleven other population samples. We picked tag SNPs using genotype data we collected in the HapMap samples and then evaluated the effective coverage of these tags in comparison to the entire set of common variants observed in the other samples. We simulated case-control association studies in the non-HapMap samples under a disease model of modest risk, and we observed little loss in power. These results demonstrate that the HapMap DNA samples can be used to select tags for genome-wide association studies in many samples around the world.     

Evaluating potential for whole-genome studies in Kosrae, an isolated population in Micronesia

Whole-genome association studies are predicted to be especially powerful in isolated populations owing to increased linkage disequilibrium (LD) and decreased allelic diversity, but this possibility has not been empirically tested. We compared genome-wide data on 113,240 SNPs typed on 30 trios from the Pacific island of Kosrae to the same markers typed in the 270 samples from the International HapMap Project. The extent of LD is longer and haplotype diversity is lower in Kosrae than in the HapMap populations. More than 98% of Kosraen haplotypes are present in HapMap populations, indicating that HapMap will be useful for genetic studies on Kosrae. The long-range LD around common alleles and limited diversity result in improved efficiency in genetic studies in this population and augments the power to detect association of 'hidden SNPs'.     

Common deletions and SNPs are in linkage disequilibrium in the human genome

Humans show great variation in phenotypic traits such as height, eye color and susceptibility to disease. Genomic DNA sequence differences among individuals are responsible for the inherited components of these complex traits. Reports suggest that intermediate and large-scale DNA copy number and structural variations are prevalent enough to be an important source of genetic variation between individuals. Because association studies to identify genomic loci associated with particular phenotypic traits have focused primarily on genotyping SNPs, it is important to determine whether common structural polymorphisms are in linkage disequilibrium with common SNPs, and thus can be assessed indirectly in SNP-based studies. Here we examine 100 deletion polymorphisms ranging from 70 bp to 7 kb. We show that common deletions and SNPs ascertained with similar criteria have essentially the same distribution of linkage disequilibrium with surrounding SNPs, indicating that these polymorphisms may share evolutionary history and that most deletion polymorphisms are effectively assayed by proxy in SNP-based association studies.     

Next-generation association studies for complex traits

A new study successfully applies complementary whole-genome sequencing and imputation approaches to establish robust disease associations in an isolated population. This strategy is poised to help elucidate the role of variants at the low end of the allele frequency spectrum in the genetic architecture of complex traits.     

Genome-wide efficient mixed-model analysis for association studies

Linear mixed models have attracted considerable attention recently as a powerful and effective tool for accounting for population stratification and relatedness in genetic association tests. However, existing methods for exact computation of standard test statistics are computationally impractical for even moderate-sized genome-wide association studies. To address this issue, several approximate methods have been proposed. Here, we present an efficient exact method, which we refer to as genome-wide efficient mixed-model association (GEMMA), that makes approximations unnecessary in many contexts. This method is approximately n times faster than the widely used exact method known as efficient mixed-model association (EMMA), where n is the sample size, making exact genome-wide association analysis computationally practical for large numbers of individuals.     

Efficiency and power in genetic association studies

We investigated selection and analysis of tag SNPs for genome-wide association studies by specifically examining the relationship between investment in genotyping and statistical power. Do pairwise or multimarker methods maximize efficiency and power? To what extent is power compromised when tags are selected from an incomplete resource such as HapMap? We addressed these questions using genotype data from the HapMap ENCODE project, association studies simulated under a realistic disease model, and empirical correction for multiple hypothesis testing. We demonstrate a haplotype-based tagging method that uniformly outperforms single-marker tests and methods for prioritization that markedly increase tagging efficiency. Examining all observed haplotypes for association, rather than just those that are proxies for known SNPs, increases power to detect rare causal alleles, at the cost of reduced power to detect common causal alleles. Power is robust to the completeness of the reference panel from which tags are selected. These findings have implications for prioritizing tag SNPs and interpreting association studies.     

Genome-wide association analyses identify 13 new susceptibility loci for generalized vitiligo

We previously reported a genome-wide association study (GWAS) identifying 14 susceptibility loci for generalized vitiligo. We report here a second GWAS (450 individuals with vitiligo (cases) and 3,182 controls), an independent replication study (1,440 cases and 1,316 controls) and a meta-analysis (3,187 cases and 6,723 controls) identifying 13 additional vitiligo-associated loci. These include OCA2-HERC2 (combined P = 3.80 × 10(-8)), MC1R (P = 1.82 × 10(-13)), a region near TYR (P = 1.57 × 10(-13)), IFIH1 (P = 4.91 × 10(-15)), CD80 (P = 3.78 × 10(-10)), CLNK (P = 1.56 × 10(-8)), BACH2 (P = 2.53 × 10(-8)), SLA (P = 1.58 × 10(-8)), CASP7 (P = 3.56 × 10(-8)), CD44 (P = 1.78 × 10(-9)), IKZF4 (P = 2.75 × 10(-14)), SH2B3 (P = 3.54 × 10(-18)) and TOB2 (P = 6.81 × 10(-10)). Most vitiligo susceptibility loci encode immunoregulatory proteins or melanocyte components that likely mediate immune targeting and the relationships among vitiligo, melanoma, and eye, skin and hair coloration.     

Principal components analysis corrects for stratification in genome-wide association studies

Population stratification--allele frequency differences between cases and controls due to systematic ancestry differences-can cause spurious associations in disease studies. We describe a method that enables explicit detection and correction of population stratification on a genome-wide scale. Our method uses principal components analysis to explicitly model ancestry differences between cases and controls. The resulting correction is specific to a candidate marker's variation in frequency across ancestral populations, minimizing spurious associations while maximizing power to detect true associations. Our simple, efficient approach can easily be applied to disease studies with hundreds of thousands of markers.     

Genome-wide association scan of tag SNPs identifies a susceptibility locus for lung cancer at 15q25.1

To identify risk variants for lung cancer, we conducted a multistage genome-wide association study. In the discovery phase, we analyzed 315,450 tagging SNPs in 1,154 current and former (ever) smoking cases of European ancestry and 1,137 frequency-matched, ever-smoking controls from Houston, Texas. For replication, we evaluated the ten SNPs most significantly associated with lung cancer in an additional 711 cases and 632 controls from Texas and 2,013 cases and 3,062 controls from the UK. Two SNPs, rs1051730 and rs8034191, mapping to a region of strong linkage disequilibrium within 15q25.1 containing PSMA4 and the nicotinic acetylcholine receptor subunit genes CHRNA3 and CHRNA5, were significantly associated with risk in both replication sets. Combined analysis yielded odds ratios of 1.32 (P < 1 x 10(-17)) for both SNPs. Haplotype analysis was consistent with there being a single risk variant in this region. We conclude that variation in a region of 15q25.1 containing nicotinic acetylcholine receptors genes contributes to lung cancer risk.     

A genome-wide association scan of nonsynonymous SNPs identifies a susceptibility variant for Crohn disease in ATG16L1

We performed a genome-wide association study of 19,779 nonsynonymous SNPs in 735 individuals with Crohn disease and 368 controls. A total of 7,159 of these SNPs were informative. We followed up on all 72 SNPs with P 0.4), these data suggest that the underlying biological process may be specific to Crohn disease.     

Integrated genotype calling and association analysis of SNPs, common copy number polymorphisms and rare CNVs

Accurate and complete measurement of single nucleotide (SNP) and copy number (CNV) variants, both common and rare, will be required to understand the role of genetic variation in disease. We present Birdsuite, a four-stage analytical framework instantiated in software for deriving integrated and mutually consistent copy number and SNP genotypes. The method sequentially assigns copy number across regions of common copy number polymorphisms (CNPs), calls genotypes of SNPs, identifies rare CNVs via a hidden Markov model (HMM), and generates an integrated sequence and copy number genotype at every locus (for example, including genotypes such as A-null, AAB and BBB in addition to AA, AB and BB calls). Such genotypes more accurately depict the underlying sequence of each individual, reducing the rate of apparent mendelian inconsistencies. The Birdsuite software is applied here to data from the Affymetrix SNP 6.0 array. Additionally, we describe a method, implemented in PLINK, to utilize these combined SNP and CNV genotypes for association testing with a phenotype.     

Extremely low-coverage sequencing and imputation increases power for genome-wide association studies

Genome-wide association studies (GWAS) have proven to be a powerful method to identify common genetic variants contributing to susceptibility to common diseases. Here, we show that extremely low-coverage sequencing (0.1-0.5×) captures almost as much of the common (>5%) and low-frequency (1-5%) variation across the genome as SNP arrays. As an empirical demonstration, we show that genome-wide SNP genotypes can be inferred at a mean r(2) of 0.71 using off-target data (0.24× average coverage) in a whole-exome study of 909 samples. Using both simulated and real exome-sequencing data sets, we show that association statistics obtained using extremely low-coverage sequencing data attain similar P values at known associated variants as data from genotyping arrays, without an excess of false positives. Within the context of reductions in sample preparation and sequencing costs, funds invested in extremely low-coverage sequencing can yield several times the effective sample size of GWAS based on SNP array data and a commensurate increase in statistical power.     

Parameters for reliable results in genetic association studies in common disease.

It is increasingly apparent that the identification of true genetic associations in common multifactorial disease will require studies comprising thousands rather than the hundreds of individuals employed to date. Using 2,873 families, we were unable to confirm a recently published association of the interleukin 12B gene in 422 type I diabetic families. These results emphasize the need for large datasets, small P values and independent replication if results are to be reliable.     

A genome-wide association scan of tag SNPs identifies a susceptibility variant for colorectal cancer at 8q24.21

Much of the variation in inherited risk of colorectal cancer (CRC) is probably due to combinations of common low risk variants. We conducted a genome-wide association study of 550,000 tag SNPs in 930 familial colorectal tumor cases and 960 controls. The most strongly associated SNP (P = 1.72 x 10(-7), allelic test) was rs6983267 at 8q24.21. To validate this finding, we genotyped rs6983267 in three additional CRC case-control series (4,361 affected individuals and 3,752 controls; 1,901 affected individuals and 1,079 controls; 1,072 affected individuals and 415 controls) and replicated the association, providing P = 1.27 x 10(-14) (allelic test) overall, with odds ratios (ORs) of 1.27 (95% confidence interval (c.i.): 1.16-1.39) and 1.47 (95% c.i.: 1.34-1.62) for heterozygotes and rare homozygotes, respectively. Analyses based on 1,477 individuals with colorectal adenoma and 2,136 controls suggest that susceptibility to CRC is mediated through development of adenomas (OR = 1.21, 95% c.i.: 1.10-1.34; P = 6.89 x 10(-5)). These data show that common, low-penetrance susceptibility alleles predispose to colorectal neoplasia.     

TLR4 mutations are associated with endotoxin hyporesponsiveness in humans

There is much variability between individuals in the response to inhaled toxins, but it is not known why certain people develop disease when challenged with environmental agents and others remain healthy. To address this, we investigated whether TLR4 (encoding the toll-like receptor-4), which has been shown to affect lipopolysaccharide (LPS) responsiveness in mice, underlies the variability in airway responsiveness to inhaled LPS in humans. Here we show that common, co-segregating missense mutations (Asp299Gly and Thr399Ile) affecting the extracellular domain of the TLR4 receptor are associated with a blunted response to inhaled LPS in humans. Transfection of THP-1 cells demonstrates that the Asp299Gly mutation (but not the Thr399Ile mutation) interrupts TLR4-mediated LPS signalling. Moreover, the wild-type allele of TLR4 rescues the LPS hyporesponsive phenotype in either primary airway epithelial cells or alveolar macrophages obtained from individuals with the TLR4 mutations. Our findings provide the first genetic evidence that common mutations in TLR4 are associated with differences in LPS responsiveness in humans, and demonstrate that gene-sequence changes can alter the ability of the host to respond to environmental stress.