Showing papers by "Michael Boehnke published in 2011"

Rare-Variant Association Testing for Sequencing Data with the Sequence Kernel Association Test

Abstract: Sequencing studies are increasingly being conducted to identify rare variants associated with complex traits. The limited power of classical single-marker association analysis for rare variants poses a central challenge in such studies. We propose the sequence kernel association test (SKAT), a supervised, flexible, computationally efficient regression method to test for association between genetic variants (common and rare) in a region and a continuous or dichotomous trait while easily adjusting for covariates. As a score-based variance-component test, SKAT can quickly calculate p values analytically by fitting the null model containing only the covariates, and so can easily be applied to genome-wide data. Using SKAT to analyze a genome-wide sequencing study of 1000 individuals, by segmenting the whole genome into 30 kb regions, requires only 7 hr on a laptop. Through analysis of simulated data across a wide range of practical scenarios and triglyceride data from the Dallas Heart Study, we show that SKAT can substantially outperform several alternative rare-variant association tests. We also provide analytic power and sample-size calculations to help design candidate-gene, whole-exome, and whole-genome sequence association studies.

Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk

Abstract: Blood pressure is a heritable trait(1) influenced by several biological pathways and responsive to environmental stimuli. Over one billion people worldwide have hypertension (>= 140 mm Hg systolic blood pressure or >= 90 mm Hg diastolic blood pressure)(2). Even small increments in blood pressure are associated with an increased risk of cardiovascular events(3). This genome-wide association study of systolic and diastolic blood pressure, which used a multi-stage design in 200,000 individuals of European descent, identified sixteen novel loci: six of these loci contain genes previously known or suspected to regulate blood pressure (GUCY1A3-GUCY1B3, NPR3-C5orf23, ADM, FURIN-FES, GOSR2, GNAS-EDN3); the other ten provide new clues to blood pressure physiology. A genetic risk score based on 29 genome-wide significant variants was associated with hypertension, left ventricular wall thickness, stroke and coronary artery disease, but not kidney disease or kidney function. We also observed associations with blood pressure in East Asian, South Asian and African ancestry individuals. Our findings provide new insights into the genetics and biology of blood pressure, and suggest potential novel therapeutic pathways for cardiovascular disease prevention.

Large-scale genome-wide association analysis of bipolar disorder identifies a new susceptibility locus near ODZ4

Abstract: We conducted a combined genome-wide association study (GWAS) of 7,481 individuals with bipolar disorder (cases) and 9,250 controls as part of the Psychiatric GWAS Consortium. Our replication study tested 34 SNPs in 4,496 independent cases with bipolar disorder and 42,422 independent controls and found that 18 of 34 SNPs had P < 0.05, with 31 of 34 SNPs having signals with the same direction of effect (P = 3.8 × 10−7). An analysis of all 11,974 bipolar disorder cases and 51,792 controls confirmed genome-wide significant evidence of association for CACNA1C and identified a new intronic variant in ODZ4. We identified a pathway comprised of subunits of calcium channels enriched in bipolar disorder association intervals. Finally, a combined GWAS analysis of schizophrenia and bipolar disorder yielded strong association evidence for SNPs in CACNA1C and in the region of NEK4-ITIH1-ITIH3-ITIH4. Our replication results imply that increasing sample sizes in bipolar disorder will confirm many additional loci.

Physical activity attenuates the influence of FTO variants on obesity risk: A meta-analysis of 218,166 adults and 19,268 children

Abstract: Background: The FTO gene harbors the strongest known susceptibility locus for obesity. While many individual studies have suggested that physical activity (PA) may attenuate the effect of FTO on obesity risk, other studies have not been able to confirm this interaction. To confirm or refute unambiguously whether PA attenuates the association of FTO with obesity risk, we meta-analyzed data from 45 studies of adults (n=218,166) and nine studies of children and adolescents (n=19,268). Methods and Findings: All studies identified to have data on the FTO rs9939609 variant (or any proxy [r(2)>0.8]) and PA were invited to participate, regardless of ethnicity or age of the participants. PA was standardized by categorizing it into a dichotomous variable (physically inactive versus active) in each study. Overall, 25% of adults and 13% of children were categorized as inactive. Interaction analyses were performed within each study by including the FTOxPA interaction term in an additive model, adjusting for age and sex. Subsequently, random effects meta-analysis was used to pool the interaction terms. In adults, the minor (A-) allele of rs9939609 increased the odds of obesity by 1.23-fold/allele (95% CI 1.20-1.26), but PA attenuated this effect (p(interaction) = 0.001). More specifically, the minor allele of rs9939609 increased the odds of obesity less in the physically active group (odds ratio = 1.22/allele, 95% CI 1.19-1.25) than in the inactive group (odds ratio = 1.30/allele, 95% CI 1.24-1.36). No such interaction was found in children and adolescents. Conclusions: The association of the FTO risk allele with the odds of obesity is attenuated by 27% in physically active adults, highlighting the importance of PA in particular in those genetically predisposed to obesity.

Genome-wide association study identifies six new loci influencing pulse pressure and mean arterial pressure

Abstract: Numerous genetic loci have been associated with systolic blood pressure (SBP) and diastolic blood pressure (DBP) in Europeans(1-3). We now report genome-wide association studies of pulse pressure (PP) and mean arterial pressure (MAP). In discovery (N = 74,064) and follow-up studies (N = 48,607), we identified at genome-wide significance (P = 2.7 x 10(-8) to P = 2.3 x 10(-13)) four new PP loci (at 4q12 near CHIC2, 7q22.3 near PIK3CG, 8q24.12 in NOV and 11q24.3 near ADAMTS8), two new MAP loci (3p21.31 in MAP4 and 10q25.3 near ADRB1) and one locus associated with both of these traits (2q24.3 near FIGN) that has also recently been associated with SBP in east Asians. For three of the new PP loci, the estimated effect for SBP was opposite of that for DBP, in contrast to the majority of common SBP- and DBP-associated variants, which show concordant effects on both traits. These findings suggest new genetic pathways underlying blood pressure variation, some of which may differentially influence SBP and DBP.

Genome-Wide Association Identifies Nine Common Variants Associated With Fasting Proinsulin Levels and Provides New Insights Into the Pathophysiology of Type 2 Diabetes

Abstract: OBJECTIVE-Proinsulin is a precursor of mature insulin and C-peptide. Higher circulating proinsulin levels are associated with impaired beta-cell function, raised glucose levels, insulin resistance, and type 2 diabetes (T2D). Studies of the insulin processing pathway could provide new insights about T2D pathophysiology. RESEARCH DESIGN AND METHODS-We have conducted a meta-analysis of genome-wide association tests of similar to 2.5 million genotyped or imputed single nucleotide polymorphisms (SNPs) and fasting proinsulin levels in 10,701 nondiabetic adults of European ancestry, with follow-up of 23 loci in up to 16,378 individuals, using additive genetic models adjusted for age, sex, fasting insulin, and study-specific covariates. RESULTS-Nine SNPs at eight loci were associated with proinsulin levels (P < 5 x 10(-8)). Two loci (LARP6 and SGSM2) have not been previously related to metabolic traits, one (MADD) has been associated with fasting glucose, one (PCSK1) has been implicated in obesity, and four (TCF7L2, SLC3OA8, VPS13C/C2CD4A/B, and ARAP1, formerly CENTD2) increase T2D risk. The proinsulin-raising allele of ARAP1 was associated with a lower fasting glucose (P = 1.7 x 10(-4)), improved beta-cell function (P = 1.1 x 10(-5)), and lower risk of T2D (odds ratio 0.88; P = 7.8 x 10(-6)). Notably, PCSK1 encodes the protein prohormone convertase 1/3, the first enzyme in the insulin processing pathway. A genotype score composed of the nine proinsulin-raising alleles was not associated with coronary disease in two large case-control datasets. CONCLUSIONS-We have identified nine genetic variants associated with fasting proinsulin. Our findings illuminate the biology underlying glucose homeostasis and T2D development in humans and argue against a direct role of proinsulin in coronary artery disease pathogenesis. Diabetes 60:2624-2634, 2011

Low-coverage sequencing: Implications for design of complex trait association studies

Abstract: New sequencing technologies allow genomic variation to be surveyed in much greater detail than previously possible. While detailed analysis of a single individual typically requires deep sequencing, when many individuals are sequenced it is possible to combine shallow sequence data across individuals to generate accurate calls in shared stretches of chromosome. Here, we show that, as progressively larger numbers of individuals are sequenced, increasingly accurate genotype calls can be generated for a given sequence depth. We evaluate the implications of low-coverage sequencing for complex trait association studies. We systematically compare study designs based on genotyping of tagSNPs, sequencing of many individuals at depths ranging between 2× and 30×, and imputation of variants discovered by sequencing a subset of individuals into the remainder of the sample. We show that sequencing many individuals at low depth is an attractive strategy for studies of complex trait genetics. For example, for disease-associated variants with frequency >0.2%, sequencing 3000 individuals at 4× depth provides similar power to deep sequencing of >2000 individuals at 30× depth but requires only ~20% of the sequencing effort. We also show low-coverage sequencing can be used to build a reference panel that can drive imputation into additional samples to increase power further. We provide guidance for investigators wishing to combine results from sequenced, genotyped, and imputed samples.

Abdominal Aortic Aneurysm Is Associated with a Variant in Low-Density Lipoprotein Receptor-Related Protein 1

Abstract: Abdominal aortic aneurysm (AAA) is a common cause of morbidity and mortality and has a significant heritability. We carried out a genome-wide association discovery study of 1866 patients with AAA and 5435 controls and replication of promising signals (lead SNP with a p value < 1 × 10(-5)) in 2871 additional cases and 32,687 controls and performed further follow-up in 1491 AAA and 11,060 controls. In the discovery study, nine loci demonstrated association with AAA (p < 1 × 10(-5)). In the replication sample, the lead SNP at one of these loci, rs1466535, located within intron 1 of low-density-lipoprotein receptor-related protein 1 (LRP1) demonstrated significant association (p = 0.0042). We confirmed the association of rs1466535 and AAA in our follow-up study (p = 0.035). In a combined analysis (6228 AAA and 49182 controls), rs1466535 had a consistent effect size and direction in all sample sets (combined p = 4.52 × 10(-10), odds ratio 1.15 [1.10-1.21]). No associations were seen for either rs1466535 or the 12q13.3 locus in independent association studies of coronary artery disease, blood pressure, diabetes, or hyperlipidaemia, suggesting that this locus is specific to AAA. Gene-expression studies demonstrated a trend toward increased LRP1 expression for the rs1466535 CC genotype in arterial tissues; there was a significant (p = 0.029) 1.19-fold (1.04-1.36) increase in LRP1 expression in CC homozygotes compared to TT homozygotes in aortic adventitia. Functional studies demonstrated that rs1466535 might alter a SREBP-1 binding site and influence enhancer activity at the locus. In conclusion, this study has identified a biologically plausible genetic variant associated specifically with AAA, and we suggest that this variant has a possible functional role in LRP1 expression.

Fine Mapping of Five Loci Associated with Low-Density Lipoprotein Cholesterol Detects Variants That Double the Explained Heritability

Abstract: Complex trait genome-wide association studies (GWAS) provide an efficient strategy for evaluating large numbers of common variants in large numbers of individuals and for identifying trait-associated variants. Nevertheless, GWAS often leave much of the trait heritability unexplained. We hypothesized that some of this unexplained heritability might be due to common and rare variants that reside in GWAS identified loci but lack appropriate proxies in modern genotyping arrays. To assess this hypothesis, we re-examined 7 genes (APOE, APOC1, APOC2, SORT1, LDLR, APOB, and PCSK9) in 5 loci associated with low-density lipoprotein cholesterol (LDL-C) in multiple GWAS. For each gene, we first catalogued genetic variation by re-sequencing 256 Sardinian individuals with extreme LDL-C values. Next, we genotyped variants identified by us and by the 1000 Genomes Project (totaling 3,277 SNPs) in 5,524 volunteers. We found that in one locus (PCSK9) the GWAS signal could be explained by a previously described low-frequency variant and that in three loci (PCSK9, APOE, and LDLR) there were additional variants independently associated with LDL-C, including a novel and rare LDLR variant that seems specific to Sardinians. Overall, this more detailed assessment of SNP variation in these loci increased estimates of the heritability of LDL-C accounted for by these genes from 3.1% to 6.5%. All association signals and the heritability estimates were successfully confirmed in a sample of ,10,000 Finnish and Norwegian individuals. Our results thus suggest that focusing on variants accessible via GWAS can lead to clear underestimates of the trait heritability explained by a set of loci. Further, our results suggest that, as prelude to large-scale sequencing efforts, targeted re-sequencing efforts paired with large-scale genotyping will increase estimates of complex trait heritability explained by known loci.

Transferability of Type 2 Diabetes Implicated Loci in Multi-Ethnic Cohorts from Southeast Asia

Abstract: Recent large genome-wide association studies (GWAS) have identified multiple loci which harbor genetic variants associated with type 2 diabetes mellitus (T2D), many of which encode proteins not previously suspected to be involved in the pathogenesis of T2D. Most GWAS for T2D have focused on populations of European descent, and GWAS conducted in other populations with different ancestry offer a unique opportunity to study the genetic architecture of T2D. We performed genome-wide association scans for T2D in 3,955 Chinese (2,010 cases, 1,945 controls), 2,034 Malays (794 cases, 1,240 controls), and 2,146 Asian Indians (977 cases, 1,169 controls). In addition to the search for novel variants implicated in T2D, these multi-ethnic cohorts serve to assess the transferability and relevance of the previous findings from European descent populations in the three major ethnic populations of Asia, comprising half of the world's population. Of the SNPs associated with T2D in previous GWAS, only variants at CDKAL1 and HHEX/IDE/KIF11 showed the strongest association with T2D in the meta-analysis including all three ethnic groups. However, consistent direction of effect was observed for many of the other SNPs in our study and in those carried out in European populations. Close examination of the associations at both the CDKAL1 and HHEX/IDE/KIF11 loci provided some evidence of locus and allelic heterogeneity in relation to the associations with T2D. We also detected variation in linkage disequilibrium between populations for most of these loci that have been previously identified. These factors, combined with limited statistical power, may contribute to the failure to detect associations across populations of diverse ethnicity. These findings highlight the value of surveying across diverse racial/ethnic groups towards the fine-mapping efforts for the casual variants and also of the search for variants, which may be population-specific.

Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk

Abstract: Blood pressure (BP) is a heritable trait1 influenced by multiple biological pathways and is responsive to environmental stimuli. Over one billion people worldwide have hypertension (BP ≥140 mm Hg systolic [SBP] or ≥90 mm Hg diastolic [DBP])2. Even small increments in BP are associated with increased risk of cardiovascular events3. This genome-wide association study of SBP and DBP, which used a multi-stage design in 200,000 individuals of European descent, identified 16 novel loci: six of these loci contain genes previously known or suspected to regulate BP (GUCY1A3-GUCY1B3; NPR3-C5orf23; ADM; FURIN-FES; GOSR2; GNAS-EDN3); the other 10 provide new clues to BP physiology. A genetic risk score based on 29 genome-wide significant variants was associated with hypertension, left ventricular wall thickness, stroke, and coronary artery disease, but not kidney disease or kidney function. We also observed associations with BP in East Asian, South Asian, and African ancestry individuals. Our findings provide new insights into the genetics and biology of BP, and suggest novel potential therapeutic pathways for cardiovascular disease prevention. Genetic approaches have advanced the understanding of biological pathways underlying inter-individual variation in BP. For example, studies of rare Mendelian BP disorders have identified multiple defects in renal sodium handling pathways4. More recently two genomewide association studies (GWAS), each of >25,000 individuals of European-ancestry, identified 13 loci associated with SBP, DBP, and hypertension5,6. We now report results of a new meta-analysis of GWAS data that includes staged follow-up genotyping to identify additional BP loci. Primary analyses evaluated associations between 2.5 million genotyped or imputed single nucleotide polymorphisms (SNPs) and SBP and DBP in 69,395 individuals of European ancestry from 29 studies (Supplementary Materials Sections 1–3, Supplementary Tables 1– 2). Following GWAS meta-analysis, we conducted a three-stage validation experiment that made efficient use of available genotyping resources, to follow up top signals in up to 133,661 additional individuals of European descent (Supplementary Fig. 1 and Supplementary Materials Section 4). Twenty-nine independent SNPs at 28 loci were significantly associated with SBP, DBP, or both in the meta-analysis combining discovery and follow up data (Fig. 1, Table 1, Supplementary Figs 2–3, Supplementary Tables 3–5). All 29 SNPs attained association P <5×10−9, an order of magnitude beyond the standard genome-wide significance level for a single stage experiment (Table 1). Sixteen of these 29 associations were novel (Table 1). Two associations were near the FURIN and GOSR2 genes; prior targeted analyses of variants in these genes suggested they Note added in proof: Since this manuscript was submitted, Kato et al published a BP GWAS in East Asians that identified a SNP highly correlated to the SNP we report at the NPR3-c5orf23 locus28. Author contributions Full author contributions and roles are listed in the Supplementary Materials Section 19. NIH Public Access Author Manuscript Nature. Author manuscript; available in PMC 2012 May 01. Published in final edited form as: Nature. ; 478(7367): 103–109. doi:10.1038/nature10405. N IH PA Athor M anscript N IH PA Athor M anscript N IH PA Athor M anscript may be BP loci7,8. At the CACNB2 locus we validated association for a previously reported6 SNP rs4373814 and detected a novel independent association for rs1813353 (pairwise r2 =0.015 in HapMap CEU). Of our 13 previously reported associations5,6, only the association at PLCD3 was not supported by the current results (Supplementary Table 4). Some of the associations are in or near genes involved in pathways known to influence BP (NPR3, GUCY1A3-GUCY1B3, ADM, GNAS-EDN3, NPPA-NPPB, and CYP17A1; Supplementary Fig. 4). Twenty-two of the 28 loci did not contain genes that were a priori strong biological candidates. As expected from prior BP GWAS results, the effects of the novel variants on SBP and DBP were small (Fig. 1 and Table 1). For all variants, the observed directions of effects were concordant for SBP, DBP, and hypertension (Fig. 1, Table 1, Supplementary Fig. 3). Among the genes at the genome-wide significant loci, only CYP17A1, previously implicated in Mendelian congenital adrenal hyperplasia and hypertension, is known to harbour rare variants that have large effects on BP9. We performed several analyses to identify potential causal alleles and mechanisms. First, we looked up the 29 genome-wide significant index SNPs and their close proxies (r2>0.8) among cis-acting expression SNP (eSNP) results from multiple tissues (Supplementary Materials Section 5). For 13/29 index SNPs, we found association between nearby eSNP variants and expression level of at least one gene transcript (10−4 > p > 10−51, Supplementary Table 6). In 5 cases, the index BP SNP and the best eSNP from a genomewide survey were identical, highlighting potential mediators of the SNP-BP associations. Second, because changes in protein sequence are strong a priori candidates to be functional, we sought non-synonymous coding SNPs that were in high LD (r2 >0.8) with the 29 index SNPs. We identified such SNPsat 8 loci (Table 1, Supplementary Materials Section 6, Supplementary Table 7). In addition we performed analyses testing for differences in genetic effect according to body mass index (BMI) or sex, and analyses of copy number variants, pathway enrichment, and metabolomic data, but we did not find any statistically significant results (Supplementary Materials Sections 7–9, Supplementary Tables 8–10). We evaluated whether the BP variants we identified in Europeans were associated with BP in individuals of East Asian (N=29,719), South Asian (N=23,977), and African (N=19,775) ancestries (Table 1, Supplementary Tables 11–13). We found significant associations in individuals of East Asian ancestry for SNPs at 9 loci and in individuals of South Asian ancestry for SNPs at 6 loci; some have been reported previously (Supplementary Tables 12 and 15). The lack of significant association for individual SNPs may reflect small sample sizes, differences in allele frequencies or LD patterns, imprecise imputation for some ancestries using existing reference samples, or a genuinely different underlying genetic architecture. Because of limited power to detect effects of individual variants in the smaller non-European samples, we created genetic risk scores for SBP and DBP incorporating all 29 BP variants weighted according to effect sizes observed in the European samples. In each non-European ancestry group, risk scores were strongly associated with SBP (P=1.1×10−40 in East Asian, P=2.9×10−13 in South Asian, P=9.8×10−4 in African ancestry individuals) and DBP (P=2.9×10−48, P=9.5×10−15, and P=5.3×10−5, respectively; Supplementary Table 13). We also created a genetic risk score to assess association of the variants in aggregate with hypertension and with clinical measures of hypertensive complications including left ventricular mass, left ventricular wall thickness, incident heart failure, incident and prevalent stroke, prevalent coronary artery disease (CAD), kidney disease, and measures of kidney function, using results from other GWAS consortia (Table 2, Supplementary Materials Sections 10–11, Supplementary Table 14). The risk score was weighted using the average of Page 2 Nature. Author manuscript; available in PMC 2012 May 01. N IH PA Athor M anscript N IH PA Athor M anscript N IH PA Athor M anscript SBP and DBP effects for the 29 SNPs. In an independent sample of 23,294 women10, an increase of 1 standard deviation in the genetic risk score was associated with a 21% increase in the odds of hypertension (95% CI 19%–28%; Table 2, Supplementary Table 14). Among individuals in the top decile of the risk score, the prevalence of hypertension was 29% compared with 16% in the bottom decile (odds ratio 2.09, 95% CI 1.86–2.36). Similar results were observed in an independent hypertension case-control sample (Table 2). In our study, individuals in the top compared to bottom quintiles of genetic risk score differed by 4.6 mm Hg SBP and 3.0 mm Hg DBP, differences that approach population-averaged BP treatment effects for a single antihypertensive agent11. Epidemiologic data have shown that differences in SBP and DBP of this magnitude, across the population range of BP, are associated with an increase in cardiovascular disease risk3. Consistent with this and in line with findings from randomized trials of BP-lowering medication in hypertensive patients12,13, the genetic risk score was positively associated with left ventricular wall thickness (P=6.0×10−6), occurrence of stroke (P=3.3×10−5) and CAD (P=8.1×10−29). The same genetic risk score was not, however, significantly associated with chronic kidney disease or measures of kidney function, even though these renal outcomes were available in a similar sample size as for the other outcomes (Table 2). The absence of association with kidney phenotypes could be explained by a weaker causal relation of BP with kidney phenotypes than with CAD and stroke. This finding is consistent with the mismatch between observational data that show a positive association of BP with kidney disease, and clinical trial data that show inconsistent evidence of benefit of BP lowering on kidney disease prevention in patients with hypertension14. Thus, several lines of evidence converge to suggest that BP elevation may in part be a consequence rather than a cause of sub-clinical kidney disease. Our discovery meta-analysis (Supplementary Fig. 2) suggests an excess of modestly significant (10−5

Common variants show predicted polygenic effects on height in the tails of the distribution, except in extremely short individuals.

Abstract: Common genetic variants have been shown to explain a fraction of the inherited variation for many common diseases and quantitative traits, including height, a classic polygenic trait. The extent to which common variation determines the phenotype of highly heritable traits such as height is uncertain, as is the extent to which common variation is relevant to individuals with more extreme phenotypes. To address these questions, we studied 1,214 individuals from the top and bottom extremes of the height distribution (tallest and shortest ∼1.5%), drawn from ∼78,000 individuals from the HUNT and FINRISK cohorts. We found that common variants still influence height at the extremes of the distribution: common variants (49/141) were nominally associated with height in the expected direction more often than is expected by chance (p<5×10⁻²⁸), and the odds ratios in the extreme samples were consistent with the effects estimated previously in population-based data. To examine more closely whether the common variants have the expected effects, we calculated a weighted allele score (WAS), which is a weighted prediction of height for each individual based on the previously estimated effect sizes of the common variants in the overall population. The average WAS is consistent with expectation in the tall individuals, but was not as extreme as expected in the shortest individuals (p<0.006), indicating that some of the short stature is explained by factors other than common genetic variation. The discrepancy was more pronounced (p<10⁻⁶) in the most extreme individuals (height<0.25 percentile). The results at the extreme short tails are consistent with a large number of models incorporating either rare genetic non-additive or rare non-genetic factors that decrease height. We conclude that common genetic variants are associated with height at the extremes as well as across the population, but that additional factors become more prominent at the shorter extreme.

A role for coding functional variants in HNF4A in type 2 diabetes susceptibility

Abstract: Aims/hypothesis Rare mutations in the gene HNF4A, encoding the transcription factor hepatocyte nuclear factor 4α (HNF-4A), account for ~5% of cases of MODY and more frequent variants in this gene may be involved in multifactorial forms of diabetes. Two low-frequency, non-synonymous variants in HNF4A (V255M, minor allele frequency [MAF] ~0.1%; T130I, MAF ~3.0%)—known to influence downstream HNF-4A target gene expression—are of interest, but previous type 2 diabetes association reports were inconclusive. We aimed to evaluate the contribution of these variants to type 2 diabetes susceptibility through large-scale association analysis.

Quantifying and correcting for the winner's curse in quantitative‐trait association studies

Abstract: Quantitative traits (QT) are an important focus of human genetic studies both because of interest in the traits themselves, and because of their role as risk factors for many human diseases. For large-scale QT association studies including genome-wide association studies (GWAS), investigators usually focus on genetic loci showing significant evidence for SNP-QT association, and genetic effect size tends to be overestimated as a consequence of the winner’s curse. In this paper, we study the impact of the winner’s curse on QT association studies in which the genetic effect size is parameterized as the slope in a linear regression model. We demonstrate by analytical calculation that the overestimation in the regression slope estimate decreases as power increases. To reduce the ascertainment bias, we propose a three-parameter maximum likelihood method and then simplify this to a one-parameter method by excluding nuisance parameters. We show that both methods reduce the bias when power to detect association is low or moderate, and that the one-parameter model generally results in smaller variance in the estimate.