Phenome- and genome-wide analyses of retinal optical coherence tomography images identify links between ocular and systemic health

The human retina is a multilayered tissue that offers a unique window into systemic health. Optical coherence tomography (OCT) is widely used in eye care and allows the noninvasive, rapid capture of retinal anatomy in exquisite detail. We conducted genotypic and phenotypic analyses of retinal layer thicknesses using macular OCT images from 44,823 UK Biobank participants. We performed OCT layer cross-phenotype association analyses (OCT-XWAS), associating retinal thicknesses with 1866 incident conditions (median 10-year follow-up) and 88 quantitative traits and blood biomarkers. We performed genome-wide association studies (GWASs), identifying inherited genetic markers that influence retinal layer thicknesses and replicated our associations among the LIFE-Adult Study (N = 6313). Last, we performed a comparative analysis of phenome- and genome-wide associations to identify putative causal links between retinal layer thicknesses and both ocular and systemic conditions. Independent associations with incident mortality were detected for thinner photoreceptor segments (PSs) and, separately, ganglion cell complex layers. Phenotypic associations were detected between thinner retinal layers and ocular, neuropsychiatric, cardiometabolic, and pulmonary conditions. A GWAS of retinal layer thicknesses yielded 259 unique loci. Consistency between epidemiologic and genetic associations suggested links between a thinner retinal nerve fiber layer with glaucoma, thinner PS with age-related macular degeneration, and poor cardiometabolic and pulmonary function with a thinner PS. In conclusion, we identified multiple inherited genetic loci and acquired systemic cardio-metabolic-pulmonary conditions associated with thinner retinal layers and identify retinal layers wherein thinning is predictive of future ocular and systemic conditions.

Leveraging type 1 diabetes human genetic and genomic data in the T1D knowledge portal

To address the challenge of translating genetic discoveries for type 1 diabetes (T1D) into mechanistic insight, we have developed the T1D Knowledge Portal (T1DKP), an open-access resource for hypothesis development and target discovery in T1D.

Mapping genomic regulation of kidney disease and traits through high-resolution and interpretable eQTLs

Expression quantitative trait locus (eQTL) studies illuminate genomic variants that regulate specific genes and contribute to fine-mapped loci discovered via genome-wide association studies (GWAS). Efforts to maximize their accuracy are ongoing. Using 240 glomerular (GLOM) and 311 tubulointerstitial (TUBE) micro-dissected samples from human kidney biopsies, we discovered 5371 GLOM and 9787 TUBE genes with at least one variant significantly associated with expression (eGene) by incorporating kidney single-nucleus open chromatin data and transcription start site distance as an "integrative prior" for Bayesian statistical fine-mapping. The use of an integrative prior resulted in higher resolution eQTLs illustrated by (1) smaller numbers of variants in credible sets with greater confidence, (2) increased enrichment of partitioned heritability for GWAS of two kidney traits, (3) an increased number of variants colocalized with the GWAS loci, and (4) enrichment of computationally predicted functional regulatory variants. A subset of variants and genes were validated experimentally in vitro and using a Drosophila nephrocyte model. More broadly, this study demonstrates that tissue-specific eQTL maps informed by single-nucleus open chromatin data have enhanced utility for diverse downstream analyses.

The Type 2 Diabetes Knowledge Portal: An open access genetic resource dedicated to type 2 diabetes and related traits

Associations between human genetic variation and clinical phenotypes have become a foundation of biomedical research. Most repositories of these data seek to be disease-agnostic and therefore lack disease-focused views. The Type 2 Diabetes Knowledge Portal (T2DKP) is a public resource of genetic datasets and genomic annotations dedicated to type 2 diabetes (T2D) and related traits. Here, we seek to make the T2DKP more accessible to prospective users and more useful to existing users. First, we evaluate the T2DKP's comprehensiveness by comparing its datasets with those of other repositories. Second, we describe how researchers unfamiliar with human genetic data can begin using and correctly interpreting them via the T2DKP. Third, we describe how existing users can extend their current workflows to use the full suite of tools offered by the T2DKP. We finally discuss the lessons offered by the T2DKP toward the goal of democratizing access to complex disease genetic results.

Leveraging type 1 diabetes human genetic and genomic data in the T1D Knowledge Portal

Translating genetic discoveries for type 1 diabetes (T1D) into mechanistic insight can reveal novel biology and therapeutic targets but remains a major challenge. We developed the T1D Knowledge Portal (T1DKP), a disease-specific resource of genetic and functional annotation data that enables users to develop hypotheses for T1D-based research and target discovery. The T1DKP can be used to query genes and genomic regions for genetic associations, identify epigenomic features, access results of bioinformatic analyses, and obtain expert-curated resources. The T1DKP is available at http://t1d.hugeamp.org .

Implicating genes, pleiotropy, and sexual dimorphism at blood lipid loci through multi-ancestry meta-analysis

BACKGROUND

Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery.

RESULTS

To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism.

CONCLUSIONS

Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

Discovery and systematic characterization of risk variants and genes for coronary artery disease in over a million participants

The discovery of genetic loci associated with complex diseases has outpaced the elucidation of mechanisms of disease pathogenesis. Here we conducted a genome-wide association study (GWAS) for coronary artery disease (CAD) comprising 181,522 cases among 1,165,690 participants of predominantly European ancestry. We detected 241 associations, including 30 new loci. Cross-ancestry meta-analysis with a Japanese GWAS yielded 38 additional new loci. We prioritized likely causal variants using functionally informed fine-mapping, yielding 42 associations with less than five variants in the 95% credible set. Similarity-based clustering suggested roles for early developmental processes, cell cycle signaling and vascular cell migration and proliferation in the pathogenesis of CAD. We prioritized 220 candidate causal genes, combining eight complementary approaches, including 123 supported by three or more approaches. Using CRISPR-Cas9, we experimentally validated the effect of an enhancer in MYO9B, which appears to mediate CAD risk by regulating vascular cell motility. Our analysis identifies and systematically characterizes >250 risk loci for CAD to inform experimental interrogation of putative causal mechanisms for CAD.

Erratum: Ancestral diversity improves discovery and fine-mapping of genetic loci for anthropometric traits-The Hispanic/Latino Anthropometry Consortium

[This corrects the article DOI: 10.1016/j.xhgg.2022.100099.].

A glomerular transcriptomic landscape of apolipoprotein L1 in Black patients with focal segmental glomerulosclerosis

Apolipoprotein L1 (APOL1)-associated focal segmental glomerulosclerosis (FSGS) is the dominant form of FSGS in Black individuals. There are no targeted therapies for this condition, in part because the molecular mechanisms underlying APOL1's pathogenic contribution to FSGS are incompletely understood. Studying the transcriptomic landscape of APOL1 FSGS in patient kidneys is an important way to discover genes and molecular behaviors that are unique or most relevant to the human disease. With the hypothesis that the pathology driven by the high-risk APOL1 genotype is reflected in alteration of gene expression across the glomerular transcriptome, we compared expression and co-expression profiles of 15,703 genes in 16 Black patients with FSGS at high-risk vs 14 Black patients with a low-risk APOL1 genotype. Expression data from APOL1-inducible HEK293 cells and normal human glomeruli were used to pursue genes and molecular pathways uncovered in these studies. We discovered increased expression of APOL1 and nine other significant differentially expressed genes in high-risk patients. This included stanniocalcin, which has a role in mitochondrial and calcium-related processes along with differential correlations between high- and low-risk APOL1 and metabolism pathway genes. There were similar correlations with extracellular matrix- and immune-related genes, but significant loss of co-expression of mitochondrial genes in high-risk FSGS, and an NF-κB-down regulating gene, NKIRAS1, as the most significant hub gene with strong differential correlations with NDUF family (mitochondrial respiratory genes) and immune-related (JAK-STAT) genes. Thus, differences in mitochondrial gene regulation appear to underlie many differences observed between high- and low-risk Black patients with FSGS.

Ancestral diversity improves discovery and fine-mapping of genetic loci for anthropometric traits-The Hispanic/Latino Anthropometry Consortium

Hispanic/Latinos have been underrepresented in genome-wide association studies (GWAS) for anthropometric traits despite their notable anthropometric variability, ancestry proportions, and high burden of growth stunting and overweight/obesity. To address this knowledge gap, we analyzed densely imputed genetic data in a sample of Hispanic/Latino adults to identify and fine-map genetic variants associated with body mass index (BMI), height, and BMI-adjusted waist-to-hip ratio (WHRadjBMI). We conducted a GWAS of 18 studies/consortia as part of the Hispanic/Latino Anthropometry (HISLA) Consortium (stage 1, n = 59,771) and generalized our findings in 9 additional studies (stage 2, n = 10,538). We conducted a trans-ancestral GWAS with summary statistics from HISLA stage 1 and existing consortia of European and African ancestries. In our HISLA stage 1 + 2 analyses, we discovered one BMI locus, as well as two BMI signals and another height signal each within established anthropometric loci. In our trans-ancestral meta-analysis, we discovered three BMI loci, one height locus, and one WHRadjBMI locus. We also identified 3 secondary signals for BMI, 28 for height, and 2 for WHRadjBMI in established loci. We show that 336 known BMI, 1,177 known height, and 143 known WHRadjBMI (combined) SNPs demonstrated suggestive transferability (nominal significance and effect estimate directional consistency) in Hispanic/Latino adults. Of these, 36 BMI, 124 height, and 11 WHRadjBMI SNPs were significant after trait-specific Bonferroni correction. Trans-ancestral meta-analysis of the three ancestries showed a small-to-moderate impact of uncorrected population stratification on the resulting effect size estimates. Our findings demonstrate that future studies may also benefit from leveraging diverse ancestries and differences in linkage disequilibrium patterns to discover novel loci and additional signals with less residual population stratification.

The Lipid Droplet Knowledge Portal: A resource for systematic analyses of lipid droplet biology

Lipid droplets (LDs) are organelles of cellular lipid storage with fundamental roles in energy metabolism and cell membrane homeostasis. There has been an explosion of research into the biology of LDs, in part due to their relevance in diseases of lipid storage, such as atherosclerosis, obesity, type 2 diabetes, and hepatic steatosis. Consequently, there is an increasing need for a resource that combines datasets from systematic analyses of LD biology. Here, we integrate high-confidence, systematically generated human, mouse, and fly data from studies on LDs in the framework of an online platform named the "Lipid Droplet Knowledge Portal" (https://lipiddroplet.org/). This scalable and interactive portal includes comprehensive datasets, across a variety of cell types, for LD biology, including transcriptional profiles of induced lipid storage, organellar proteomics, genome-wide screen phenotypes, and ties to human genetics. This resource is a powerful platform that can be utilized to identify determinants of lipid storage.

Rare coding variants in 35 genes associate with circulating lipid levels-A multi-ancestry analysis of 170,000 exomes

Large-scale gene sequencing studies for complex traits have the potential to identify causal genes with therapeutic implications. We performed gene-based association testing of blood lipid levels with rare (minor allele frequency < 1%) predicted damaging coding variation by using sequence data from >170,000 individuals from multiple ancestries: 97,493 European, 30,025 South Asian, 16,507 African, 16,440 Hispanic/Latino, 10,420 East Asian, and 1,182 Samoan. We identified 35 genes associated with circulating lipid levels; some of these genes have not been previously associated with lipid levels when using rare coding variation from population-based samples. We prioritize 32 genes in array-based genome-wide association study (GWAS) loci based on aggregations of rare coding variants; three (EVI5, SH2B3, and PLIN1) had no prior association of rare coding variants with lipid levels. Most of our associated genes showed evidence of association among multiple ancestries. Finally, we observed an enrichment of gene-based associations for low-density lipoprotein cholesterol drug target genes and for genes closest to GWAS index single-nucleotide polymorphisms (SNPs). Our results demonstrate that gene-based associations can be beneficial for drug target development and provide evidence that the gene closest to the array-based GWAS index SNP is often the functional gene for blood lipid levels.

Determinants of penetrance and variable expressivity in monogenic metabolic conditions across 77,184 exomes

Hundreds of thousands of genetic variants have been reported to cause severe monogenic diseases, but the probability that a variant carrier develops the disease (termed penetrance) is unknown for virtually all of them. Additionally, the clinical utility of common polygenetic variation remains uncertain. Using exome sequencing from 77,184 adult individuals (38,618 multi-ancestral individuals from a type 2 diabetes case-control study and 38,566 participants from the UK Biobank, for whom genotype array data were also available), we apply clinical standard-of-care gene variant curation for eight monogenic metabolic conditions. Rare variants causing monogenic diabetes and dyslipidemias display effect sizes significantly larger than the top 1% of the corresponding polygenic scores. Nevertheless, penetrance estimates for monogenic variant carriers average 60% or lower for most conditions. We assess epidemiologic and genetic factors contributing to risk prediction in monogenic variant carriers, demonstrating that inclusion of polygenic variation significantly improves biomarker estimation for two monogenic dyslipidemias.

The trans-ancestral genomic architecture of glycemic traits

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Exome sequencing of 20,791 cases of type 2 diabetes and 24,440 controls

Protein-coding genetic variants that strongly affect disease risk can yield relevant clues to disease pathogenesis. Here we report exome-sequencing analyses of 20,791 individuals with type 2 diabetes (T2D) and 24,440 non-diabetic control participants from 5 ancestries. We identify gene-level associations of rare variants (with minor allele frequencies of less than 0.5%) in 4 genes at exome-wide significance, including a series of more than 30 SLC30A8 alleles that conveys protection against T2D, and in 12 gene sets, including those corresponding to T2D drug targets (P = 6.1 × 10-3) and candidate genes from knockout mice (P = 5.2 × 10-3). Within our study, the strongest T2D gene-level signals for rare variants explain at most 25% of the heritability of the strongest common single-variant signals, and the gene-level effect sizes of the rare variants that we observed in established T2D drug targets will require 75,000-185,000 sequenced cases to achieve exome-wide significance. We propose a method to interpret these modest rare-variant associations and to incorporate these associations into future target or gene prioritization efforts.

Mexican Carriers of the HNF1A p.E508K Variant Do Not Experience an Enhanced Response to Sulfonylureas

OBJECTIVE

To assess whether an ethnic-specific variant (p.E508K) in the maturity-onset diabetes of the young (MODY) gene hepatocyte nuclear factor-1α (HNF1A) found in Mexicans is associated with higher sensitivity to sulfonylureas, as documented in patients with MODY3.

RESEARCH DESIGN AND METHODS

We recruited 96 participants (46 variant carriers and 50 age- and sex-matched noncarriers). Response to glipizide (one 2.5-5.0-mg dose), metformin (four 500-mg doses), and an oral glucose challenge was evaluated using a previously validated protocol. Glucose and insulin levels and their areas under the curve (AUCs) were compared between groups.

RESULTS

Carriers of the p.E508K variant had a lower maximum insulin peak during the glipizide challenge as compared with noncarriers with diabetes (P < 0.05). Also, carriers had a lower insulin response after the oral glucose challenge. Following an oral glucose tolerance test in the presence of metformin, carriers of the p.E508K variant with diabetes had a lower maximum insulin peak and total and incremental insulin AUC value as compared with noncarriers with diabetes (P < 0.05). A similar but nonsignificant trend was seen in participants without type 2 diabetes.

CONCLUSIONS

Carriers of variant p.E508K in HNF1A have a reduced insulin response rather than the increased sensitivity to sulfonylureas seen in patients with MODY3.

Refining the accuracy of validated target identification through coding variant fine-mapping in type 2 diabetes

We aggregated coding variant data for 81,412 type 2 diabetes cases and 370,832 controls of diverse ancestry, identifying 40 coding variant association signals (P < 2.2 × 10-7); of these, 16 map outside known risk-associated loci. We make two important observations. First, only five of these signals are driven by low-frequency variants: even for these, effect sizes are modest (odds ratio ≤1.29). Second, when we used large-scale genome-wide association data to fine-map the associated variants in their regional context, accounting for the global enrichment of complex trait associations in coding sequence, compelling evidence for coding variant causality was obtained for only 16 signals. At 13 others, the associated coding variants clearly represent 'false leads' with potential to generate erroneous mechanistic inference. Coding variant associations offer a direct route to biological insight for complex diseases and identification of validated therapeutic targets; however, appropriate mechanistic inference requires careful specification of their causal contribution to disease predisposition.

Erratum: Sequence data and association statistics from 12,940 type 2 diabetes cases and controls

This corrects the article DOI: 10.1038/sdata.2017.179.

Sequence data and association statistics from 12,940 type 2 diabetes cases and controls

To investigate the genetic basis of type 2 diabetes (T2D) to high resolution, the GoT2D and T2D-GENES consortia catalogued variation from whole-genome sequencing of 2,657 European individuals and exome sequencing of 12,940 individuals of multiple ancestries. Over 27M SNPs, indels, and structural variants were identified, including 99% of low-frequency (minor allele frequency [MAF] 0.1-5%) non-coding variants in the whole-genome sequenced individuals and 99.7% of low-frequency coding variants in the whole-exome sequenced individuals. Each variant was tested for association with T2D in the sequenced individuals, and, to increase power, most were tested in larger numbers of individuals (>80% of low-frequency coding variants in ~82 K Europeans via the exome chip, and ~90% of low-frequency non-coding variants in ~44 K Europeans via genotype imputation). The variants, genotypes, and association statistics from these analyses provide the largest reference to date of human genetic information relevant to T2D, for use in activities such as T2D-focused genotype imputation, functional characterization of variants or genes, and other novel analyses to detect associations between sequence variation and T2D.

A Low-Frequency Inactivating AKT2 Variant Enriched in the Finnish Population Is Associated With Fasting Insulin Levels and Type 2 Diabetes Risk

To identify novel coding association signals and facilitate characterization of mechanisms influencing glycemic traits and type 2 diabetes risk, we analyzed 109,215 variants derived from exome array genotyping together with an additional 390,225 variants from exome sequence in up to 39,339 normoglycemic individuals from five ancestry groups. We identified a novel association between the coding variant (p.Pro50Thr) in AKT2 and fasting plasma insulin (FI), a gene in which rare fully penetrant mutations are causal for monogenic glycemic disorders. The low-frequency allele is associated with a 12% increase in FI levels. This variant is present at 1.1% frequency in Finns but virtually absent in individuals from other ancestries. Carriers of the FI-increasing allele had increased 2-h insulin values, decreased insulin sensitivity, and increased risk of type 2 diabetes (odds ratio 1.05). In cellular studies, the AKT2-Thr50 protein exhibited a partial loss of function. We extend the allelic spectrum for coding variants in AKT2 associated with disorders of glucose homeostasis and demonstrate bidirectional effects of variants within the pleckstrin homology domain of AKT2.

Genome-wide association studies in the Japanese population identify seven novel loci for type 2 diabetes

Genome-wide association studies (GWAS) have identified more than 80 susceptibility loci for type 2 diabetes (T2D), but most of its heritability still remains to be elucidated. In this study, we conducted a meta-analysis of GWAS for T2D in the Japanese population. Combined data from discovery and subsequent validation analyses (23,399 T2D cases and 31,722 controls) identify 7 new loci with genome-wide significance (P<5 × 10(-8)), rs1116357 near CCDC85A, rs147538848 in FAM60A, rs1575972 near DMRTA1, rs9309245 near ASB3, rs67156297 near ATP8B2, rs7107784 near MIR4686 and rs67839313 near INAFM2. Of these, the association of 4 loci with T2D is replicated in multi-ethnic populations other than Japanese (up to 65,936 T2Ds and 158,030 controls, P<0.007). These results indicate that expansion of single ethnic GWAS is still useful to identify novel susceptibility loci to complex traits not only for ethnicity-specific loci but also for common loci across different ethnicities.

Power in the phenotypic extremes: a simulation study of power in discovery and replication of rare variants

Next-generation sequencing technologies are making it possible to study the role of rare variants in human disease. Many studies balance statistical power with cost-effectiveness by (a) sampling from phenotypic extremes and (b) utilizing a two-stage design. Two-stage designs include a broad-based discovery phase and selection of a subset of potential causal genes/variants to be further examined in independent samples. We evaluate three parameters: first, the gain in statistical power due to extreme sampling to discover causal variants; second, the informativeness of initial (Phase I) association statistics to select genes/variants for follow-up; third, the impact of extreme and random sampling in (Phase 2) replication. We present a quantitative method to select individuals from the phenotypic extremes of a binary trait, and simulate disease association studies under a variety of sample sizes and sampling schemes. First, we find that while studies sampling from extremes have excellent power to discover rare variants, they have limited power to associate them to phenotype—suggesting high false-negative rates for upcoming studies. Second, consistent with previous studies, we find that the effect sizes estimated in these studies are expected to be systematically larger compared with the overall population effect size; in a well-cited lipids study, we estimate the reported effect to be twofold larger. Third, replication studies require large samples from the general population to have sufficient power; extreme sampling could reduce the required sample size as much as fourfold. Our observations offer practical guidance for the design and interpretation of studies that utilize extreme sampling.

Identification and functional characterization of G6PC2 coding variants influencing glycemic traits define an effector transcript at the G6PC2-ABCB11 locus

Genome wide association studies (GWAS) for fasting glucose (FG) and insulin (FI) have identified common variant signals which explain 4.8% and 1.2% of trait variance, respectively. It is hypothesized that low-frequency and rare variants could contribute substantially to unexplained genetic variance. To test this, we analyzed exome-array data from up to 33,231 non-diabetic individuals of European ancestry. We found exome-wide significant (P<5×10-7) evidence for two loci not previously highlighted by common variant GWAS: GLP1R (p.Ala316Thr, minor allele frequency (MAF)=1.5%) influencing FG levels, and URB2 (p.Glu594Val, MAF = 0.1%) influencing FI levels. Coding variant associations can highlight potential effector genes at (non-coding) GWAS signals. At the G6PC2/ABCB11 locus, we identified multiple coding variants in G6PC2 (p.Val219Leu, p.His177Tyr, and p.Tyr207Ser) influencing FG levels, conditionally independent of each other and the non-coding GWAS signal. In vitro assays demonstrate that these associated coding alleles result in reduced protein abundance via proteasomal degradation, establishing G6PC2 as an effector gene at this locus. Reconciliation of single-variant associations and functional effects was only possible when haplotype phase was considered. In contrast to earlier reports suggesting that, paradoxically, glucose-raising alleles at this locus are protective against type 2 diabetes (T2D), the p.Val219Leu G6PC2 variant displayed a modest but directionally consistent association with T2D risk. Coding variant associations for glycemic traits in GWAS signals highlight PCSK1, RREB1, and ZHX3 as likely effector transcripts. These coding variant association signals do not have a major impact on the trait variance explained, but they do provide valuable biological insights.

Understanding how FI and FG levels are regulated is important because their derangement is a feature of T2D. Despite recent success from GWAS in identifying regions of the genome influencing glycemic traits, collectively these loci explain only a small proportion of trait variance. Unlocking the biological mechanisms driving these associations has been challenging because the vast majority of variants map to non-coding sequence, and the genes through which they exert their impact are largely unknown. In the current study, we sought to increase our understanding of the physiological pathways influencing both traits using exome-array genotyping in up to 33,231 non-diabetic individuals to identify coding variants and consequently genes associated with either FG or FI levels. We identified novel association signals for both traits including the receptor for GLP-1 agonists which are a widely used therapy for T2D. Furthermore, we identified coding variants at several GWAS loci which point to the genes underlying these association signals. Importantly, we found that multiple coding variants in G6PC2 result in a loss of protein function and lower fasting glucose levels.

Association of a low-frequency variant in HNF1A with type 2 diabetes in a Latino population

IMPORTANCE

Latino populations have one of the highest prevalences of type 2 diabetes worldwide.

OBJECTIVES

To investigate the association between rare protein-coding genetic variants and prevalence of type 2 diabetes in a large Latino population and to explore potential molecular and physiological mechanisms for the observed relationships.

DESIGN, SETTING, AND PARTICIPANTS

Whole-exome sequencing was performed on DNA samples from 3756 Mexican and US Latino individuals (1794 with type 2 diabetes and 1962 without diabetes) recruited from 1993 to 2013. One variant was further tested for allele frequency and association with type 2 diabetes in large multiethnic data sets of 14,276 participants and characterized in experimental assays.

MAIN OUTCOME AND MEASURES

Prevalence of type 2 diabetes. Secondary outcomes included age of onset, body mass index, and effect on protein function.

RESULTS

A single rare missense variant (c.1522G>A [p.E508K]) was associated with type 2 diabetes prevalence (odds ratio [OR], 5.48; 95% CI, 2.83-10.61; P = 4.4 × 10(-7)) in hepatocyte nuclear factor 1-α (HNF1A), the gene responsible for maturity onset diabetes of the young type 3 (MODY3). This variant was observed in 0.36% of participants without type 2 diabetes and 2.1% of participants with it. In multiethnic replication data sets, the p.E508K variant was seen only in Latino patients (n = 1443 with type 2 diabetes and 1673 without it) and was associated with type 2 diabetes (OR, 4.16; 95% CI, 1.75-9.92; P = .0013). In experimental assays, HNF-1A protein encoding the p.E508K mutant demonstrated reduced transactivation activity of its target promoter compared with a wild-type protein. In our data, carriers and noncarriers of the p.E508K mutation with type 2 diabetes had no significant differences in compared clinical characteristics, including age at onset. The mean (SD) age for carriers was 45.3 years (11.2) vs 47.5 years (11.5) for noncarriers (P = .49) and the mean (SD) BMI for carriers was 28.2 (5.5) vs 29.3 (5.3) for noncarriers (P = .19).

CONCLUSIONS AND RELEVANCE

Using whole-exome sequencing, we identified a single low-frequency variant in the MODY3-causing gene HNF1A that is associated with type 2 diabetes in Latino populations and may affect protein function. This finding may have implications for screening and therapeutic modification in this population, but additional studies are required.

Loss-of-function mutations in SLC30A8 protect against type 2 diabetes

Loss-of-function mutations protective against human disease provide in vivo validation of therapeutic targets, but none have yet been described for type 2 diabetes (T2D). Through sequencing or genotyping of ~150,000 individuals across 5 ancestry groups, we identified 12 rare protein-truncating variants in SLC30A8, which encodes an islet zinc transporter (ZnT8) and harbors a common variant (p.Trp325Arg) associated with T2D risk and glucose and proinsulin levels. Collectively, carriers of protein-truncating variants had 65% reduced T2D risk (P = 1.7 × 10(-6)), and non-diabetic Icelandic carriers of a frameshift variant (p.Lys34Serfs*50) demonstrated reduced glucose levels (-0.17 s.d., P = 4.6 × 10(-4)). The two most common protein-truncating variants (p.Arg138* and p.Lys34Serfs*50) individually associate with T2D protection and encode unstable ZnT8 proteins. Previous functional study of SLC30A8 suggested that reduced zinc transport increases T2D risk, and phenotypic heterogeneity was observed in mouse Slc30a8 knockouts. In contrast, loss-of-function mutations in humans provide strong evidence that SLC30A8 haploinsufficiency protects against T2D, suggesting ZnT8 inhibition as a therapeutic strategy in T2D prevention.

Plasma HDL cholesterol and risk of myocardial infarction: a mendelian randomisation study

BACKGROUND

High plasma HDL cholesterol is associated with reduced risk of myocardial infarction, but whether this association is causal is unclear. Exploiting the fact that genotypes are randomly assigned at meiosis, are independent of non-genetic confounding, and are unmodified by disease processes, mendelian randomisation can be used to test the hypothesis that the association of a plasma biomarker with disease is causal.

METHODS

We performed two mendelian randomisation analyses. First, we used as an instrument a single nucleotide polymorphism (SNP) in the endothelial lipase gene (LIPG Asn396Ser) and tested this SNP in 20 studies (20,913 myocardial infarction cases, 95,407 controls). Second, we used as an instrument a genetic score consisting of 14 common SNPs that exclusively associate with HDL cholesterol and tested this score in up to 12,482 cases of myocardial infarction and 41,331 controls. As a positive control, we also tested a genetic score of 13 common SNPs exclusively associated with LDL cholesterol.

FINDINGS

Carriers of the LIPG 396Ser allele (2·6% frequency) had higher HDL cholesterol (0·14 mmol/L higher, p=8×10(-13)) but similar levels of other lipid and non-lipid risk factors for myocardial infarction compared with non-carriers. This difference in HDL cholesterol is expected to decrease risk of myocardial infarction by 13% (odds ratio [OR] 0·87, 95% CI 0·84-0·91). However, we noted that the 396Ser allele was not associated with risk of myocardial infarction (OR 0·99, 95% CI 0·88-1·11, p=0·85). From observational epidemiology, an increase of 1 SD in HDL cholesterol was associated with reduced risk of myocardial infarction (OR 0·62, 95% CI 0·58-0·66). However, a 1 SD increase in HDL cholesterol due to genetic score was not associated with risk of myocardial infarction (OR 0·93, 95% CI 0·68-1·26, p=0·63). For LDL cholesterol, the estimate from observational epidemiology (a 1 SD increase in LDL cholesterol associated with OR 1·54, 95% CI 1·45-1·63) was concordant with that from genetic score (OR 2·13, 95% CI 1·69-2·69, p=2×10(-10)).

INTERPRETATION

Some genetic mechanisms that raise plasma HDL cholesterol do not seem to lower risk of myocardial infarction. These data challenge the concept that raising of plasma HDL cholesterol will uniformly translate into reductions in risk of myocardial infarction.

FUNDING

US National Institutes of Health, The Wellcome Trust, European Union, British Heart Foundation, and the German Federal Ministry of Education and Research.

The metabochip, a custom genotyping array for genetic studies of metabolic, cardiovascular, and anthropometric traits

Genome-wide association studies have identified hundreds of loci for type 2 diabetes, coronary artery disease and myocardial infarction, as well as for related traits such as body mass index, glucose and insulin levels, lipid levels, and blood pressure. These studies also have pointed to thousands of loci with promising but not yet compelling association evidence. To establish association at additional loci and to characterize the genome-wide significant loci by fine-mapping, we designed the "Metabochip," a custom genotyping array that assays nearly 200,000 SNP markers. Here, we describe the Metabochip and its component SNP sets, evaluate its performance in capturing variation across the allele-frequency spectrum, describe solutions to methodological challenges commonly encountered in its analysis, and evaluate its performance as a platform for genotype imputation. The metabochip achieves dramatic cost efficiencies compared to designing single-trait follow-up reagents, and provides the opportunity to compare results across a range of related traits. The metabochip and similar custom genotyping arrays offer a powerful and cost-effective approach to follow-up large-scale genotyping and sequencing studies and advance our understanding of the genetic basis of complex human diseases and traits.

Stratifying type 2 diabetes cases by BMI identifies genetic risk variants in LAMA1 and enrichment for risk variants in lean compared to obese cases

Common diseases such as type 2 diabetes are phenotypically heterogeneous. Obesity is a major risk factor for type 2 diabetes, but patients vary appreciably in body mass index. We hypothesized that the genetic predisposition to the disease may be different in lean (BMI<25 Kg/m²) compared to obese cases (BMI≥30 Kg/m²). We performed two case-control genome-wide studies using two accepted cut-offs for defining individuals as overweight or obese. We used 2,112 lean type 2 diabetes cases (BMI<25 kg/m²) or 4,123 obese cases (BMI≥30 kg/m²), and 54,412 un-stratified controls. Replication was performed in 2,881 lean cases or 8,702 obese cases, and 18,957 un-stratified controls. To assess the effects of known signals, we tested the individual and combined effects of SNPs representing 36 type 2 diabetes loci. After combining data from discovery and replication datasets, we identified two signals not previously reported in Europeans. A variant (rs8090011) in the LAMA1 gene was associated with type 2 diabetes in lean cases (P = 8.4×10⁻⁹, OR = 1.13 [95% CI 1.09-1.18]), and this association was stronger than that in obese cases (P = 0.04, OR = 1.03 [95% CI 1.00-1.06]). A variant in HMG20A--previously identified in South Asians but not Europeans--was associated with type 2 diabetes in obese cases (P = 1.3×10⁻⁸, OR = 1.11 [95% CI 1.07-1.15]), although this association was not significantly stronger than that in lean cases (P = 0.02, OR = 1.09 [95% CI 1.02-1.17]). For 36 known type 2 diabetes loci, 29 had a larger odds ratio in the lean compared to obese (binomial P = 0.0002). In the lean analysis, we observed a weighted per-risk allele OR = 1.13 [95% CI 1.10-1.17], P = 3.2×10⁻¹⁴. This was larger than the same model fitted in the obese analysis where the OR = 1.06 [95% CI 1.05-1.08], P = 2.2×10⁻¹⁶. This study provides evidence that stratification of type 2 diabetes cases by BMI may help identify additional risk variants and that lean cases may have a stronger genetic predisposition to type 2 diabetes.

A genome-wide association search for type 2 diabetes genes in African Americans

African Americans are disproportionately affected by type 2 diabetes (T2DM) yet few studies have examined T2DM using genome-wide association approaches in this ethnicity. The aim of this study was to identify genes associated with T2DM in the African American population. We performed a Genome Wide Association Study (GWAS) using the Affymetrix 6.0 array in 965 African-American cases with T2DM and end-stage renal disease (T2DM-ESRD) and 1029 population-based controls. The most significant SNPs (n = 550 independent loci) were genotyped in a replication cohort and 122 SNPs (n = 98 independent loci) were further tested through genotyping three additional validation cohorts followed by meta-analysis in all five cohorts totaling 3,132 cases and 3,317 controls. Twelve SNPs had evidence of association in the GWAS (P<0.0071), were directionally consistent in the Replication cohort and were associated with T2DM in subjects without nephropathy (P<0.05). Meta-analysis in all cases and controls revealed a single SNP reaching genome-wide significance (P<2.5×10(-8)). SNP rs7560163 (P = 7.0×10(-9), OR (95% CI) = 0.75 (0.67-0.84)) is located intergenically between RND3 and RBM43. Four additional loci (rs7542900, rs4659485, rs2722769 and rs7107217) were associated with T2DM (P<0.05) and reached more nominal levels of significance (P<2.5×10(-5)) in the overall analysis and may represent novel loci that contribute to T2DM. We have identified novel T2DM-susceptibility variants in the African-American population. Notably, T2DM risk was associated with the major allele and implies an interesting genetic architecture in this population. These results suggest that multiple loci underlie T2DM susceptibility in the African-American population and that these loci are distinct from those identified in other ethnic populations.

Comparing strategies to fine-map the association of common SNPs at chromosome 9p21 with type 2 diabetes and myocardial infarction

Noncoding variants at human chromosome 9p21 near CDKN2A and CDKN2B are associated with type 2 diabetes, myocardial infarction, aneurysm, vertical cup disc ratio and at least five cancers. Here we compare approaches to more comprehensively assess genetic variation in the region. We carried out targeted sequencing at high coverage in 47 individuals and compared the results to pilot data from the 1000 Genomes Project. We imputed variants into type 2 diabetes and myocardial infarction cohorts directly from targeted sequencing, from a genotyped reference panel derived from sequencing and from 1000 Genomes Project low-coverage data. Polymorphisms with frequency >5% were captured well by all strategies. Imputation of intermediate-frequency polymorphisms required a higher density of tag SNPs in disease samples than is available on first-generation genome-wide association study (GWAS) arrays. Our association analyses identified more comprehensive sets of variants showing equivalent statistical association with type 2 diabetes or myocardial infarction, but did not identify stronger associations than the original GWAS signals.

Whole-genome amplification using Φ29 DNA polymerase

The cornerstones of any genetic analysis study are the quality and quantity of the DNA samples. DNA is a precious limited resource, and in human disease studies the accessibility of sample DNA is often governed by the isolation method and the human source. Additionally, forensic analysis and archaeological research are generally infeasible without intact sample DNA. Therefore, mechanisms to preserve or enhance the quantity of the DNA stock are crucial to the success of these studies. Historically, to preserve and maintain DNA stocks, costly and labor-intensive Epstein-Barr-virus-transformed cell lines were produced. The creation of cell lines can be valuable for a number of reasons in addition to creating a renewable resource of DNA, but the cost and effort to create them, as well as the requirement of intact cells to begin with, limit the utility of this approach. More recently, whole-genome amplification (WGA), utilizing the unique property of the enzyme Φ29 DNA polymerase, has been used to generate robust high-fidelity copies of the genome. As described in this protocol, WGA using Φ29 DNA polymerase allows unbiased representation of the genome via multiple-strand displacement, followed by rolling-circle amplification on random primers.

The major genetic determinants of HIV-1 control affect HLA class I peptide presentation

Infectious and inflammatory diseases have repeatedly shown strong genetic associations within the major histocompatibility complex (MHC); however, the basis for these associations remains elusive. To define host genetic effects on the outcome of a chronic viral infection, we performed genome-wide association analysis in a multiethnic cohort of HIV-1 controllers and progressors, and we analyzed the effects of individual amino acids within the classical human leukocyte antigen (HLA) proteins. We identified >300 genome-wide significant single-nucleotide polymorphisms (SNPs) within the MHC and none elsewhere. Specific amino acids in the HLA-B peptide binding groove, as well as an independent HLA-C effect, explain the SNP associations and reconcile both protective and risk HLA alleles. These results implicate the nature of the HLA-viral peptide interaction as the major factor modulating durable control of HIV infection.

IGF2R missense single-nucleotide polymorphisms and breast cancer risk: the multiethnic cohort study

IGF2R has been proposed to be a tumor suppressor gene given its antagonist role on cellular growth and evidence of loss of heterozygosity in several cancers, including breast cancer. To investigate whether inherited differences in potentially functional IGF2R variants influence the risk of breast cancer, we sequenced 46 exons of IGF2R to identify novel missense single-nucleotide polymorphisms (SNP) and tested 12 missense SNPs for their associations with breast cancer risk among 1,614 breast cancer cases and 1,960 controls from the Multiethnic Cohort. None of these missense SNPs were significantly associated with breast cancer risk. Our findings provide no evidence that missense SNPs in IGF2R influence breast cancer susceptibility

Common variants at CD40 and other loci confer risk of rheumatoid arthritis

To identify rheumatoid arthritis risk loci in European populations, we conducted a meta-analysis of two published genome-wide association (GWA) studies totaling 3,393 cases and 12,462 controls. We genotyped 31 top-ranked SNPs not previously associated with rheumatoid arthritis in an independent replication of 3,929 autoantibody-positive rheumatoid arthritis cases and 5,807 matched controls from eight separate collections. We identified a common variant at the CD40 gene locus (rs4810485, P = 0.0032 replication, P = 8.2 x 10(-9) overall, OR = 0.87). Along with other associations near TRAF1 (refs. 2,3) and TNFAIP3 (refs. 4,5), this implies a central role for the CD40 signaling pathway in rheumatoid arthritis pathogenesis. We also identified association at the CCL21 gene locus (rs2812378, P = 0.00097 replication, P = 2.8 x 10(-7) overall), a gene involved in lymphocyte trafficking. Finally, we identified evidence of association at four additional gene loci: MMEL1-TNFRSF14 (rs3890745, P = 0.0035 replication, P = 1.1 x 10(-7) overall), CDK6 (rs42041, P = 0.010 replication, P = 4.0 x 10(-6) overall), PRKCQ (rs4750316, P = 0.0078 replication, P = 4.4 x 10(-6) overall), and KIF5A-PIP4K2C (rs1678542, P = 0.0026 replication, P = 8.8 x 10(-8) overall).

Polymorphisms associated with cholesterol and risk of cardiovascular events

BACKGROUND

Common single-nucleotide polymorphisms (SNPs) that are associated with blood low-density lipoprotein (LDL) or high-density lipoprotein (HDL) cholesterol modestly affect lipid levels. We tested the hypothesis that a combination of such SNPs contributes to the risk of cardiovascular disease.

METHODS

We studied SNPs at nine loci in 5414 subjects from the cardiovascular cohort of the Malmö Diet and Cancer Study. We first validated the association between SNPs and either LDL or HDL cholesterol and subsequently created a genotype score on the basis of the number of unfavorable alleles. We used Cox proportional-hazards models to determine the time to the first cardiovascular event in relation to the genotype score.

RESULTS

All nine SNPs showed replication of an association with levels of either LDL or HDL cholesterol. With increasing genotype scores, the level of LDL cholesterol increased from 152 mg to 171 mg per deciliter (3.9 to 4.4 mmol per liter), whereas HDL cholesterol decreased from 60 mg to 51 mg per deciliter (1.6 to 1.3 mmol per liter). During follow-up (median, 10.6 years), 238 subjects had a first cardiovascular event. The genotype score was associated with incident cardiovascular disease in models adjusted for covariates including baseline lipid levels (P<0.001). The use of the genotype score did not improve the clinical risk prediction, as assessed by the C statistic. However, there was a significant improvement in risk classification with the use of models that included the genotype score, as compared with those that did not include the genotype score.

CONCLUSIONS

A genotype score of nine validated SNPs that are associated with modulation in levels of LDL or HDL cholesterol was an independent risk factor for incident cardiovascular disease. The score did not improve risk discrimination but did modestly improve clinical risk reclassification for individual subjects beyond standard clinical factors.

Association testing of common variants in the insulin receptor substrate-1 gene (IRS1) with type 2 diabetes

AIMS/HYPOTHESIS

Activation of the insulin receptor substrate-1 (IRS1) is a key initial step in the insulin signalling pathway. Despite several reports of association of the G972R polymorphism in its gene IRS1 with type 2 diabetes, we and others have not observed this association in well-powered samples. However, other nearby variants might account for the putative association signal.

SUBJECTS AND METHODS

We characterised the haplotype map of IRS1 and selected 20 markers designed to capture common variations in the region. We genotyped this comprehensive set of markers in several family-based and case-control samples of European descent totalling 12,129 subjects.

RESULTS

In an initial sample of 2,235 North American and Polish case-control pairs, the minor allele of the rs934167 polymorphism showed nominal evidence of association with type 2 diabetes (odds ratio [OR] 1.25, 95% CI 1.03-1.51, p = 0.03). This association showed a trend in the same direction in 7,659 Scandinavian samples (OR 1.16, 95% CI 0.96-1.39, p = 0.059). The combined OR was 1.20 (p = 0.008), but statistical correction for the number of variants examined yielded a p value of 0.086. We detected no differences across rs934167 genotypes in insulin-related quantitative traits.

CONCLUSIONS/INTERPRETATION

Our data do not support an association of common variants in IRS1 with type 2 diabetes in populations of European descent.

Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels

New strategies for prevention and treatment of type 2 diabetes (T2D) require improved insight into disease etiology. We analyzed 386,731 common single-nucleotide polymorphisms (SNPs) in 1464 patients with T2D and 1467 matched controls, each characterized for measures of glucose metabolism, lipids, obesity, and blood pressure. With collaborators (FUSION and WTCCC/UKT2D), we identified and confirmed three loci associated with T2D-in a noncoding region near CDKN2A and CDKN2B, in an intron of IGF2BP2, and an intron of CDKAL1-and replicated associations near HHEX and in SLC30A8 found by a recent whole-genome association study. We identified and confirmed association of a SNP in an intron of glucokinase regulatory protein (GCKR) with serum triglycerides. The discovery of associated variants in unsuspected genes and outside coding regions illustrates the ability of genome-wide association studies to provide potentially important clues to the pathogenesis of common diseases.

The Krüppel-like factor 11 (KLF11) Q62R polymorphism is not associated with type 2 diabetes in 8,676 people

Krüppel-like factor 11 is a pancreatic transcription factor whose activity induces the insulin gene. A common glutamine-to-arginine change at codon 62 (Q62R) in its gene KLF11 has been recently associated with type 2 diabetes in two independent samples. Q62R and two other rare missense variants (A347S and T220M) were also shown to affect the function of KLF11 in vitro, and insulin levels were lower in carriers of the minor allele at Q62R. We therefore examined their impact on common type 2 diabetes in several family-based and case-control samples of northern-European ancestry, totaling 8,676 individuals. We did not detect the rare A347S and T220M variants in our samples. With respect to Q62R, despite >99% power to detect an association of the previously published magnitude, Q62R was not associated with type 2 diabetes (pooled odds ratio 0.97 [95% CI 0.88-1.08], P = 0.63). In a subset of normoglycemic individuals, we did not observe significant differences in various insulin traits according to genotype at KLF11 Q62R. We conclude that the KLF11 A347S and T220M mutations do not contribute to increased risk of diabetes in European-derived populations and that the Q62R polymorphism has, at best, a minor effect on diabetes risk.

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.

Common single nucleotide polymorphisms in TCF7L2 are reproducibly associated with type 2 diabetes and reduce the insulin response to glucose in nondiabetic individuals

Recently, common noncoding variants in the TCF7L2 gene were strongly associated with increased risk of type 2 diabetes in samples from Iceland, Denmark, and the U.S. We genotyped 13 single nucleotide polymorphisms (SNPs) across TCF7L2 in 8,310 individuals in family-based and case-control designs from Scandinavia, Poland, and the U.S. We convincingly confirmed the previous association of TCF7L2 SNPs with the risk of type 2 diabetes (rs7903146T odds ratio 1.40 [95% CI 1.30-1.50], P = 6.74 x 10(-20)). In nondiabetic individuals, the risk genotypes were associated with a substantial reduction in the insulinogenic index derived from an oral glucose tolerance test (risk allele homozygotes have half the insulin response to glucose of noncarriers, P = 0.003) but not with increased insulin resistance. These results suggest that TCF7L2 variants may act through insulin secretion to increase the risk of type 2 diabetes.

Haplotype structures and large-scale association testing of the 5' AMP-activated protein kinase genes PRKAA2, PRKAB1, and PRKAB2 [corrected] with type 2 diabetes

AMP-activated protein kinase (AMPK) is a key molecular regulator of cellular metabolism, and its activity is induced by both metformin and thiazolidinedione antidiabetic medications. It has therefore been proposed both as a putative agent in the pathophysiology of type 2 diabetes and as a valid target for therapeutic intervention. Thus, the genes that encode the various AMPK subunits are intriguing candidates for the inherited basis of type 2 diabetes. We therefore set out to test for the association of common variants in the genes that encode three selected AMPK subunits with type 2 diabetes and related phenotypes. Of the seven genes that encode AMPK isoforms, we initially chose PRKAA2, PRKAB1, and PRKAB2 because of their higher prior probability of association with type 2 diabetes, based on previous reports of genetic linkage, functional molecular studies, expression patterns, and pharmacological evidence. We determined their haplotype structure, selected a subset of tag single nucleotide polymorphisms that comprehensively capture the extent of common genetic variation in these genes, and genotyped them in family-based and case/control samples comprising 4,206 individuals. Analysis of single-marker and multi-marker tests revealed no association with type 2 diabetes, fasting plasma glucose, or insulin sensitivity. Several nominal associations of variants in PRKAA2 and PRKAB1 with BMI appear to be consistent with statistical noise.

High-density haplotype structure and association testing of the insulin-degrading enzyme (IDE) gene with type 2 diabetes in 4,206 people

The insulin-degrading enzyme is responsible for the intracellular proteolysis of insulin. Its gene IDE is located on chromosome 10, in an area with suggestive linkage to type 2 diabetes and related phenotypes. Due to the impact of genetic variants of this gene in rodents and the function of its protein product, it has been proposed as a candidate gene for type 2 diabetes. Various groups have explored the role of the common genetic variation of IDE on insulin resistance and reported associations of various single nucleotide polymorphisms (SNPs) and haplotypes on both type 2 diabetes and glycemic traits. We sought to characterize the haplotype structure of IDE in detail and replicate the association of common variants with type 2 diabetes, fasting insulin, fasting glucose, and insulin resistance. We assessed linkage disequilibrium, selected single-marker and multimarker tags, and genotyped these markers in several case-control and family-based samples totalling 4,206 Caucasian individuals. We observed no statistically significant evidence of association between single-marker or multimarker tests in IDE and type 2 diabetes. Nominally significant differences in quantitative traits are consistent with statistical noise. We conclude that common genetic variation at IDE is unlikely to confer clinically significant risk of type 2 diabetes in Caucasians.

A candidate gene approach to searching for low-penetrance breast and prostate cancer genes

Most cases of breast and prostate cancer are not associated with mutations in known high-penetrance genes, indicating the involvement of multiple low-penetrance risk alleles. Studies that have attempted to identify these genes have met with limited success. The National Cancer Institute Breast and Prostate Cancer Cohort Consortium--a pooled analysis of multiple large cohort studies with a total of more than 5,000 cases of breast cancer and 8,000 cases of prostate cancer--was therefore initiated. The goal of this consortium is to characterize variations in approximately 50 genes that mediate two pathways that are associated with these cancers--the steroid-hormone metabolism pathway and the insulin-like growth factor signalling pathway--and to associate these variations with cancer risk.

Association of common variation in the HNF1alpha gene region with risk of type 2 diabetes

It is currently unclear how often genes that are mutated to cause rare, early-onset monogenic forms of disease also harbor common variants that contribute to the more typical polygenic form of each disease. The gene for MODY3 diabetes, HNF1alpha, lies in a region that has shown linkage to late-onset type 2 diabetes (12q24, NIDDM2), and previous association studies have suggested a weak trend toward association for common missense variants in HNF1alpha with glucose-related traits. Based on genotyping of 79 common SNPs in the 118 kb spanning HNF1alpha, we selected 21 haplotype tag single nucleotide polymorphisms (SNPs) and genotyped them in >4,000 diabetic patients and control subjects from Sweden, Finland, and Canada. Several SNPs from the coding region and 5' of the gene demonstrated nominal association with type 2 diabetes, with the most significant marker (rs1920792) having an odds ratio of 1.17 and a P value of 0.002. We then genotyped three SNPs with the strongest evidence for association to type 2 diabetes (rs1920792, I27L, and A98V) in an additional 4,400 type 2 diabetic and control subjects from North America and Poland and compared our results with those of the original sample and of Weedon et al. None of the results were consistently observed across all samples, with the possible exception of a modest association of the rare (3-5%) A98V variant. These results indicate that common variants in HNF1alpha either play no role in type 2 diabetes, a very small role, or a role that cannot be consistently observed without consideration of as yet unmeasured genetic or environmental modifiers.

Association testing of the protein tyrosine phosphatase 1B gene (PTPN1) with type 2 diabetes in 7,883 people

Protein tyrosine phosphatase (PTP)-1B, encoded by the PTPN1 gene, inactivates the insulin signal transduction cascade by dephosphorylating phosphotyrosine residues in insulin signaling molecules. Due to its chromosomal location under a chromosome 20 linkage peak and the metabolic effects of its absence in knockout mice, it is a candidate gene for type 2 diabetes. Recent studies have associated common sequence variants in PTPN1 with type 2 diabetes and diabetes-related phenotypes. We sought to replicate the association of common single nucleotide polymorphisms (SNPs) and haplotypes in PTPN1 with type 2 diabetes, fasting plasma glucose, and insulin sensitivity in a large collection of subjects. We assessed linkage disequilibrium, selected tag SNPs, and typed these markers in 3,347 cases of type 2 diabetes and 3,347 control subjects as well as 1,189 siblings discordant for type 2 diabetes. Despite power estimated at >95% to replicate the previously reported associations, no statistically significant evidence of association was observed between PTPN1 SNPs or common haplotypes with type 2 diabetes or with diabetic phenotypes.

Clarifying the PROGINS Allele Association in Ovarian and Breast Cancer Risk: A Haplotype-Based Analysis

BACKGROUND

The PROGINS allele of the progesterone receptor (PGR) gene has been associated with an increased risk of ovarian cancer and a decreased risk of breast cancer. We set out to refine the association between common variation at the PGR gene locus and these two diseases.

METHODS

We characterized the haplotype structure of the PGR gene by genotyping 54 single-nucleotide polymorphisms (SNPs) in 349 women. We then selected a subset of 17 haplotype-tagging SNPs that captured variation across the locus and typed them in 267 ovarian cancer case patients and 397 control subjects from two case-control studies and in 1715 breast cancer case patients and 2505 control subjects from a cohort study.

RESULTS

The PGR locus was characterized by four blocks of strong linkage disequilibrium. Two SNPs in block 4 were associated with an increased risk of ovarian cancer among homozygous carriers as compared with noncarriers: rs1042838 (PROGINS allele; odds ratio [OR] = 3.23, 95% confidence interval [CI] = 1.19 to 8.75, P = .022) and rs608995 (minor allele; OR = 3.10, 95% CI = 1.63 to 5.89, P<.001). The PROGINS allele was observed on a subset of chromosomes carrying the minor allele at rs608995, and its association with ovarian cancer was fully explained by its association with rs608995. In addition, rs608995 fell on two common haplotypes (4-D and 4-E), whose association with ovarian cancer was the same as that of rs608995. These same two haplotypes were associated with a non-statistically significantly reduced risk of breast cancer.

CONCLUSIONS

Variation in PGR was associated with ovarian cancer risk, although the strongest result was not with the PROGINS allele. Instead, any causal allele(s) are likely in or downstream of block 4 and carried on haplotypes 4-D and 4-E. There was some evidence that the same variation was associated with a reduced risk of breast cancer, but the association was not statistically significant.

5' flanking variants of resistin are associated with obesity

Diabetes and obesity have long been known to be related. The recently characterized adipocyte hormone resistin (also called FIZZ3/ADSF) has been implicated as a molecular link between impaired glucose tolerance (IGT) and obesity in mice. A search for sequence variants at the human resistin locus identified nine single-nucleotide polymorphisms (SNPs) but no coding variants. An investigation into the association of these SNPs with diabetes and obesity revealed two 5' flanking variants (g.-537 and g.-420), in strong linkage disequilibrium, that are associated with BMI. In nondiabetic individuals from the Quebec City area and the Saguenay-Lac-St-Jean region of Quebec, the g.-537 mutation (allelic frequency = 0.04) was significantly associated with an increase in BMI (P = 0.03 and P = 0.01, respectively). When the data from these two populations were combined and adjusted for age and sex, both the g.-537 (odds ratio [OR] 2.72, 95% CI 1.28-5.81) and the g.-420 variants (1.58, 1.06-2.35) were associated with an increased risk for a BMI > or =30 kg/m(2). In contrast, in case/control and family-based study populations from Scandinavia, we saw no effect on BMI with either of these promoter variants. No association was seen with diabetes in any of the population samples.

Genomewide linkage analysis of stature in multiple populations reveals several regions with evidence of linkage to adult height

Genomewide linkage analysis has been extremely successful at identification of the genetic variation underlying single-gene disorders. However, linkage analysis has been less successful for common human diseases and other complex traits in which multiple genetic and environmental factors interact to influence disease risk. We hypothesized that a highly heritable complex trait, in which the contribution of environmental factors was relatively limited, might be more amenable to linkage analysis. We therefore chose to study stature (adult height), for which heritability is approximately 75%-90% (Phillips and Matheny 1990; Carmichael and McGue 1995; Preece 1996; Silventoinen et al. 2000). We reanalyzed genomewide scans from four populations for which genotype and height data were available, using a variance-components method implemented in GENEHUNTER 2.0 (Pratt et al. 2000). The populations consisted of 408 individuals in 58 families from the Botnia region of Finland, 753 individuals in 183 families from other parts of Finland, 746 individuals in 179 families from Southern Sweden, and 420 individuals in 63 families from the Saguenay-Lac-St.-Jean region of Quebec. Four regions showed evidence of linkage to stature: 6q24-25, multipoint LOD score 3.85 at marker D6S1007 in Botnia (genomewide P<.06), 7q31.3-36 (LOD 3.40 at marker D7S2195 in Sweden, P<.02), 12p11.2-q14 (LOD 3.35 at markers D12S10990-D12S398 in Finland, P<.05) and 13q32-33 (LOD 3.56 at markers D13S779-D13S797 in Finland, P<.05). In a companion article (Perola et al. 2001 [in this issue]), strong supporting evidence is obtained for linkage to the region on chromosome 7. These studies suggest that highly heritable complex traits such as stature may be genetically tractable and provide insight into the genetic architecture of complex traits.