Momoko Horikoshi

SNPs in KCNQ1 are associated with susceptibility to type 2 diabetes in East Asian and European populations

We conducted a genome-wide association study using 207,097 SNP markers in Japanese individuals with type 2 diabetes and unrelated controls, and identified KCNQ1 (potassium voltage-gated channel, KQT-like subfamily, member 1) to be a strong candidate for conferring susceptibility to type 2 diabetes. We detected consistent association of a SNP in KCNQ1 (rs2283228) with the disease in several independent case-control studies (additive model P = 3.1 × 10−12; OR = 1.26, 95% CI = 1.18–1.34). Several other SNPs in the same linkage disequilibrium (LD) block were strongly associated with type 2 diabetes (additive model: rs2237895, P = 7.3 × 10−9; OR = 1.32, 95% CI = 1.20–1.45, rs2237897, P = 6.8 × 10−13; OR = 1.41, 95% CI = 1.29–1.55). The association of these SNPs with type 2 diabetes was replicated in samples from Singaporean (additive model: rs2237895, P = 8.5 × 10−3; OR = 1.14, rs2237897, P = 2.4 × 10−4; OR = 1.22) and Danish populations (additive model: rs2237895, P = 3.7 × 10−11; OR = 1.24, rs2237897, P = 1.2 × 10−4; OR = 1.36).

A genome-wide association study in the Japanese population identifies susceptibility loci for type 2 diabetes at UBE2E2 and C2CD4A-C2CD4B

We conducted a genome-wide association study of type 2 diabetes (T2D) using 459,359 SNPs in a Japanese population with a three-stage study design (stage 1, 4,470 cases and 3,071 controls; stage 2, 2,886 cases and 3,087 controls; stage 3, 3,622 cases and 2,356 controls). We identified new associations in UBE2E2 on chromosome 3 and in C2CD4A-C2CD4B on chromosome 15 at genome-wide significant levels (rs7612463 in UBE2E2, combined P = 2.27 × 10−9; rs7172432 in C2CD4A-C2CD4B, combined P = 3.66 × 10−9). The association of these two loci with T2D was replicated in other east Asian populations. In the European populations, the C2CD4A-C2CD4B locus was significantly associated with T2D, and a combined analysis of all populations gave P = 8.78 × 10−14, whereas the UBE2E2 locus did not show association to T2D. In conclusion, we identified two new loci at UBE2E2 and C2CD4A-C2CD4B associated with susceptibility to T2D.

A genetic variation of the transcription factor 7-like 2 gene is associated with risk of type 2 diabetes in the Japanese population

Aims/hypothesisIt has been suggested that transcription factor 7-like 2 protein (TCF7L2) plays an important role in glucose metabolism by regulating the production level of glucagon-like peptide-1, a hormone which modifies glucose-dependent insulin secretion. Recently, variants of TCF7L2 gene were reported to confer an increased risk of type 2 diabetes in three different samples from European and European-origin populations. We studied whether the single nucleotide polymorphisms (SNPs) in TCF7L2 were associated with type 2 diabetes in samples from a Japanese population.MethodsFive SNPs were genotyped in three different sample sets. Association with type 2 diabetes was investigated in each, as well as in combined sample sets.ResultsThe SNP rs7903146 was nominally associated with type 2 diabetes in the initial (p = 0.08) and two replication sample sets (p = 0.05 and 0.06). For the combined sample set, in which we successfully genotyped 1,174 type 2 diabetes patients and 823 control subjects, rs7903146 showed a significant association with type 2 diabetes (odds ratio = 1.69 [95% CI 1.21–2.36], p = 0.002) with the same direction as the previous reports in samples from European and European-origin populations. SNPs rs7903146 and rs7901695 were in complete linkage disequilibrium. The rest of the five SNPs (rs7895340, rs11196205 and rs12255372) did not show any significant associations with type 2 diabetes.Conclusions/interpretationThe consistent association between rs7903146 in TCF7L2 and type 2 diabetes in different ethnic groups, including the Japanese population, suggests that TCF7L2 is a common susceptibility gene for type 2 diabetes.

The impact of low-frequency and rare variants on lipid levels

Using a genome-wide screen of 9.6 million genetic variants achieved through 1000 Genomes Project imputation in 62,166 samples, we identify association to lipid traits in 93 loci, including 79 previously identified loci with new lead SNPs and 10 new loci, 15 loci with a low-frequency lead SNP and 10 loci with a missense lead SNP, and 2 loci with an accumulation of rare variants. In six loci, SNPs with established function in lipid genetics (CELSR2, GCKR, LIPC and APOE) or candidate missense mutations with predicted damaging function (CD300LG and TM6SF2) explained the locus associations. The low-frequency variants increased the proportion of variance explained, particularly for low-density lipoprotein cholesterol and total cholesterol. Altogether, our results highlight the impact of low-frequency variants in complex traits and show that imputation offers a cost-effective alternative to resequencing.

Genome-wide association study identifies three novel loci for type 2 diabetes

Although over 60 loci for type 2 diabetes (T2D) have been identified, there still remains a large genetic component to be clarified. To explore unidentified loci for T2D, we performed a genome-wide association study (GWAS) of 6 209 637 single-nucleotide polymorphisms (SNPs), which were directly genotyped or imputed using East Asian references from the 1000 Genomes Project (June 2011 release) in 5976 Japanese patients with T2D and 20 829 nondiabetic individuals. Nineteen unreported loci were selected and taken forward to follow-up analyses. Combined discovery and follow-up analyses (30 392 cases and 34 814 controls) identified three new loci with genome-wide significance, which were MIR129-LEP [rs791595; risk allele = A; risk allele frequency (RAF) = 0.080; P = 2.55 × 10(-13); odds ratio (OR) = 1.17], GPSM1 [rs11787792; risk allele = A; RAF = 0.874; P = 1.74 × 10(-10); OR = 1.15] and SLC16A13 (rs312457; risk allele = G; RAF = 0.078; P = 7.69 × 10(-13); OR = 1.20). This study demonstrates that GWASs based on the imputation of genotypes using modern reference haplotypes such as that from the 1000 Genomes Project data can assist in identification of new loci for common diseases.

Expression-based genome-wide association study links the receptor CD44 in adipose tissue with type 2 diabetes

Type 2 diabetes (T2D) is a complex, polygenic disease affecting nearly 300 million people worldwide. T2D is primarily characterized by insulin resistance, and growing evidence has indicated the causative link between adipose tissue inflammation and the development of insulin resistance. Genetic association studies have successfully revealed a number of important genes consistently associated with T2D to date. However, these robust T2D-associated genes do not fully elucidate the mechanisms underlying the development and progression of the disease. Here, we report an alternative approach, gene expression-based genome-wide association study (eGWAS): searching for genes repeatedly implicated in functional microarray experiments (often publicly available). We performed an eGWAS across 130 independent experiments (totally 1,175 T2D case-control microarrays) to find additional genes implicated in the molecular pathogenesis of T2D and identified the immune-cell receptor CD44 as our top candidate (P = 8.5 × 10−20). We found CD44 deficiency in a diabetic mouse model ameliorates insulin resistance and adipose tissue inflammation and also found that anti-CD44 antibody treatment decreases blood glucose levels and adipose tissue macrophage accumulation in a high-fat, diet-fed mouse model. Further, in humans, we observed CD44 is expressed in inflammatory cells in obese adipose tissue and discovered serum CD44 levels were positively correlated with insulin resistance and glycemic control. CD44 likely plays a causative role in the development of adipose tissue inflammation and insulin resistance in rodents and humans. Genes repeatedly implicated in publicly available experimental data may have unique functionally important roles in T2D and other complex diseases.

A single nucleotide polymorphism in ANK1 is associated with susceptibility to type 2 diabetes in Japanese populations

To identify a novel susceptibility locus for type 2 diabetes, we performed an imputation-based, genome-wide association study (GWAS) in a Japanese population using newly obtained imputed-genotype data for 2 229 890 single-nucleotide polymorphisms (SNPs) estimated from previously reported, directly genotyped GWAS data in the same samples (stage 1: 4470 type 2 diabetes versus 3071 controls). We directly genotyped 43 new SNPs with P-values of <10(-4) in a part of stage-1 samples (2692 type 2 diabetes versus 3071 controls), and the associations of validated SNPs were evaluated in another 11 139 Japanese individuals (stage 2: 7605 type 2 diabetes versus 3534 controls). Combined meta-analysis using directly genotyped data for stages 1 and 2 revealed that rs515071 in ANK1 and rs7656416 near MGC21675 were associated with type 2 diabetes in the Japanese population at the genome-wide significant level (P < 5 × 10(-8)). The association of rs515071 was also observed in European GWAS data (combined P for all populations = 6.14 × 10(-10)). Rs7656416 was in linkage disequilibrium to rs6815464, which had recently been identified as a top signal in a meta-analysis of East Asian GWAS for type 2 diabetes (r(2) = 0.76 in stage 2). The association of rs7656416 with type 2 diabetes disappeared after conditioning on rs6815464. These results indicate that the ANK1 locus is a new, common susceptibility locus for type 2 diabetes across different ethnic groups. The signal of association was weaker in the directly genotyped data, so the improvement in signal indicates the importance of imputation in this particular case.

Genome-wide association study of primary tooth eruption identifies pleiotropic loci associated with height and craniofacial distances

Twin and family studies indicate that the timing of primary tooth eruption is highly heritable, with estimates typically exceeding 80%. To identify variants involved in primary tooth eruption, we performed a population-based genome-wide association study of ‘age at first tooth’ and ‘number of teeth’ using 5998 and 6609 individuals, respectively, from the Avon Longitudinal Study of Parents and Children (ALSPAC) and 5403 individuals from the 1966 Northern Finland Birth Cohort (NFBC1966). We tested 2 446 724 SNPs imputed in both studies. Analyses were controlled for the effect of gestational age, sex and age of measurement. Results from the two studies were combined using fixed effects inverse variance meta-analysis. We identified a total of 15 independent loci, with 10 loci reaching genome-wide significance (P < 5 × 10−8) for ‘age at first tooth’ and 11 loci for ‘number of teeth’. Together, these associations explain 6.06% of the variation in ‘age of first tooth’ and 4.76% of the variation in ‘number of teeth’. The identified loci included eight previously unidentified loci, some containing genes known to play a role in tooth and other developmental pathways, including an SNP in the protein-coding region of BMP4 (rs17563, P = 9.080 × 10−17). Three of these loci, containing the genes HMGA2, AJUBA and ADK, also showed evidence of association with craniofacial distances, particularly those indexing facial width. Our results suggest that the genome-wide association approach is a powerful strategy for detecting variants involved in tooth eruption, and potentially craniofacial growth and more generally organ development.

Absence of an association between the polymorphisms in the genes encoding adiponectin receptors and type 2 diabetes

Aims/hypothesisSecreted by adipocytes, adiponectin is a hormone that acts as an antidiabetic and anti-atherogenic adipokine. We recently cloned the genes encoding two adiponectin receptors (ADIPOR1 and ADIPOR2). The aim of this study was to examine whether ADIPOR1 and/or ADIPOR2 play a major role in genetic susceptibility to insulin resistance or type 2 diabetes in the Japanese population.MethodsBy direct sequencing and a search of public databases, we identified single nucleotide polymorphisms (SNPs) in ADIPOR1 and ADIPOR2, and investigated whether these SNPs are associated with insulin resistance and type 2 diabetes in the Japanese population.ResultsThe linkage disequilibrium (LD) in the chromosomal region of ADIPOR1 was almost completely preserved, whereas the LD in ADIPOR2 was less well preserved. None of the SNPs in ADIPOR1 or ADIPOR2 were significantly associated with insulin resistance or type 2 diabetes. No differences in ADIPOR1 or ADIPOR2 haplotype frequencies were observed between type 2 diabetic and non-diabetic subjects.Conclusions/interpretationGenetic variations in ADIPOR1 or ADIPOR2 are unlikely to lead to a common genetic predisposition to insulin resistance or type 2 diabetes in the Japanese population.

Discovery and fine-mapping of glycaemic and obesity-related trait loci using high-density imputation

Reference panels from the 1000 Genomes (1000G) Project Consortium provide near complete coverage of common and low-frequency genetic variation with minor allele frequency ≥0.5% across European ancestry populations. Within the European Network for Genetic and Genomic Epidemiology (ENGAGE) Consortium, we have undertaken the first large-scale meta-analysis of genome-wide association studies (GWAS), supplemented by 1000G imputation, for four quantitative glycaemic and obesity-related traits, in up to 87,048 individuals of European ancestry. We identified two loci for body mass index (BMI) at genome-wide significance, and two for fasting glucose (FG), none of which has been previously reported in larger meta-analysis efforts to combine GWAS of European ancestry. Through conditional analysis, we also detected multiple distinct signals of association mapping to established loci for waist-hip ratio adjusted for BMI (RSPO3) and FG (GCK and G6PC2). The index variant for one association signal at the G6PC2 locus is a low-frequency coding allele, H177Y, which has recently been demonstrated to have a functional role in glucose regulation. Fine-mapping analyses revealed that the non-coding variants most likely to drive association signals at established and novel loci were enriched for overlap with enhancer elements, which for FG mapped to promoter and transcription factor binding sites in pancreatic islets, in particular. Our study demonstrates that 1000G imputation and genetic fine-mapping of common and low-frequency variant association signals at GWAS loci, integrated with genomic annotation in relevant tissues, can provide insight into the functional and regulatory mechanisms through which their effects on glycaemic and obesity-related traits are mediated.