Background

Rare variants have gathered increasing attention as a possible alternative source of missing heritability. Since next generation sequencing technology is not yet cost-effective for large-scale genomic studies, a widely used alternative approach is imputation. However, the imputation approach may be limited by the low accuracy of the imputed rare variants. To improve imputation accuracy of rare variants, various approaches have been suggested, including increasing the sample size of the reference panel, using sequencing data from study-specific samples (i.e., specific populations), and using local reference panels by genotyping or sequencing a subset of study samples. While these approaches mainly utilize reference panels, imputation accuracy of rare variants can also be increased by using exome chips containing rare variants. The exome chip contains 250 K rare variants selected from the discovered variants of about 12,000 sequenced samples. If exome chip data are available for previously genotyped samples, the combined approach using a genotype panel of merged data, including exome chips and SNP chips, should increase the imputation accuracy of rare variants.

Results

In this study, we describe a combined imputation which uses both exome chip and SNP chip data simultaneously as a genotype panel. The effectiveness and performance of the combined approach was demonstrated using a reference panel of 848 samples constructed using exome sequencing data from the T2D-GENES consortium and 5,349 sample genotype panels consisting of an exome chip and SNP chip. As a result, the combined approach increased imputation quality up to 11 %, and genomic coverage for rare variants up to 117.7 % (MAF < 1 %), compared to imputation using the SNP chip alone. Also, we investigated the systematic effect of reference panels on imputation quality using five reference panels and three genotype panels. The best performing approach was the combination of the study specific reference panel and the genotype panel of combined data.

Conclusions

Our study demonstrates that combined datasets, including SNP chips and exome chips, enhances both the imputation quality and genomic coverage of rare variants.

Background

Genome-wide single marker and gene-based meta-analyses of long-term average (LTA) blood pressure (BP) phenotypes may reveal novel findings for BP.

Methods and results

We conducted genome-wide analysis among 18 422 East Asian participants (stage 1) followed by replication study of ≤46 629 participants of European ancestry (stage 2). Significant single-nucleotide polymorphisms and genes were determined by a P<5.0×10^-8 and 2.5×10^-6, respectively, in joint analyses of stage-1 and stage-2 data. We identified 1 novel ARL3 variant, rs4919669 at 10q24.32, influencing LTA systolic BP (stage-1 P=5.03×10^-8, stage-2 P=8.64×10^-3, joint P=2.63×10^-8) and mean arterial pressure (stage-1 P=3.59×10^-9, stage-2 P=2.35×10^-2, joint P=2.64×10^-8). Three previously reported BP loci (WBP1L, NT5C2, and ATP2B1) were also identified for all BP phenotypes. Gene-based analysis provided the first robust evidence for association of KCNJ11 with LTA systolic BP (stage-1 P=8.55×10^-6, stage-2 P=1.62×10^-5, joint P=3.28×10^-9) and mean arterial pressure (stage-1 P=9.19×10^-7, stage-2 P=9.69×10^-5, joint P=2.15×10^-9) phenotypes. Fourteen genes (TMEM180, ACTR1A, SUFU, ARL3, SFXN2, WBP1L, CYP17A1, C10orf32, C10orf32-ASMT, AS3MT, CNNM2, and NT5C2 at 10q24.32; ATP2B1 at 12q21.33; and NCR3LG1 at 11p15.1) implicated by previous genome-wide association study meta-analyses were also identified. Among the loci identified by the previous genome-wide association study meta-analysis of LTA BP, we transethnically replicated associations of the KCNK3 marker rs1275988 at 2p23.3 with LTA systolic BP and mean arterial pressure phenotypes (P=1.27×10^-4 and 3.30×10^-4, respectively).

Conclusions

We identified 1 novel variant and 1 novel gene and present the first direct evidence of relevance of the KCNK3 locus for LTA BP among East Asians.

Hypertension is a leading cause of global disease, mortality, and disability. While individuals of African descent suffer a disproportionate burden of hypertension and its complications, they have been underrepresented in genetic studies. To identify novel susceptibility loci for blood pressure and hypertension in people of African ancestry, we performed both single and multiple-trait genome-wide association analyses. We analyzed 21 genome-wide association studies comprised of 31,968 individuals of African ancestry, and validated our results with additional 54,395 individuals from multi-ethnic studies. These analyses identified nine loci with eleven independent variants which reached genome-wide significance (P < 1.25×10-8) for either systolic and diastolic blood pressure, hypertension, or for combined traits. Single-trait analyses identified two loci (TARID/TCF21 and LLPH/TMBIM4) and multiple-trait analyses identified one novel locus (FRMD3) for blood pressure. At these three loci, as well as at GRP20/CDH17, associated variants had alleles common only in African-ancestry populations. Functional annotation showed enrichment for genes expressed in immune and kidney cells, as well as in heart and vascular cells/tissues. Experiments driven by these findings and using angiotensin-II induced hypertension in mice showed altered kidney mRNA expression of six genes, suggesting their potential role in hypertension. Our study provides new evidence for genes related to hypertension susceptibility, and the need to study African-ancestry populations in order to identify biologic factors contributing to hypertension.

Human height is associated with risk of multiple diseases and is profoundly determined by an individual's genetic makeup and shows a high degree of ethnic heterogeneity. Large-scale genome-wide association (GWA) analyses of adult height in Europeans have identified nearly 180 genetic loci. A recent study showed high replicability of results from Europeans-based GWA studies in Asians; however, population-specific loci may exist due to distinct linkage disequilibrium patterns. We carried out a GWA meta-analysis in 93 926 individuals from East Asia. We identified 98 loci, including 17 novel and 81 previously reported loci, associated with height at P < 5 × 10(-8), together explaining 8.89% of phenotypic variance. Among the newly identified variants, 10 are commonly distributed (minor allele frequency, MAF > 5%) in Europeans, with comparable frequencies with in Asians, and 7 single-nucleotide polymorphisms are with low frequency (MAF < 5%) in Europeans. In addition, our data suggest that novel biological pathway such as the protein tyrosine phosphatase family is involved in regulation of height. The findings from this study considerably expand our knowledge of the genetic architecture of human height in Asians.

Sixty genetic loci associated with abdominal obesity, measured by waist circumference (WC) and waist-hip ratio (WHR), have been previously identified, primarily from studies conducted in European-ancestry populations. We conducted a meta-analysis of associations of abdominal obesity with approximately 2.5 million single nucleotide polymorphisms (SNPs) among 53,052 (for WC) and 48,312 (for WHR) individuals of Asian descent, and replicated 33 selected SNPs among 3,762 to 17,110 additional individuals. We identified four novel loci near the EFEMP1, ADAMTSL3 , CNPY2, and GNAS genes that were associated with WC after adjustment for body mass index (BMI); two loci near the NID2 and HLA-DRB5 genes associated with WHR after adjustment for BMI, and three loci near the CEP120, TSC22D2, and SLC22A2 genes associated with WC without adjustment for BMI. Functional enrichment analyses revealed enrichment of corticotropin-releasing hormone signaling, GNRH signaling, and/or CDK5 signaling pathways for those newly-identified loci. Our study provides additional insight on genetic contribution to abdominal obesity.

Meta-analyses of genome-wide association studies (GWAS) have identified more than 240 loci that are associated with type 2 diabetes (T2D)^1,2; however, most of these loci have been identified in analyses of individuals with European ancestry. Here, to examine T2D risk in East Asian individuals, we carried out a meta-analysis of GWAS data from 77,418 individuals with T2D and 356,122 healthy control individuals. In the main analysis, we identified 301 distinct association signals at 183 loci, and across T2D association models with and without consideration of body mass index and sex, we identified 61 loci that are newly implicated in predisposition to T2D. Common variants associated with T2D in both East Asian and European populations exhibited strongly correlated effect sizes. Previously undescribed associations include signals in or near GDAP1, PTF1A, SIX3, ALDH2, a microRNA cluster, and genes that affect the differentiation of muscle and adipose cells³. At another locus, expression quantitative trait loci at two overlapping T2D signals affect two genes-NKX6-3 and ANK1-in different tissues^4-6. Association studies in diverse populations identify additional loci and elucidate disease-associated genes, biology, and pathways.

Recent genetic association studies have identified 55 genetic loci associated with obesity or body mass index (BMI). The vast majority, 51 loci, however, were identified in European-ancestry populations. We conducted a meta-analysis of associations between BMI and ∼2.5 million genotyped or imputed single nucleotide polymorphisms among 86 757 individuals of Asian ancestry, followed by in silico and de novo replication among 7488-47 352 additional Asian-ancestry individuals. We identified four novel BMI-associated loci near the KCNQ1 (rs2237892, P = 9.29 × 10(-13)), ALDH2/MYL2 (rs671, P = 3.40 × 10(-11); rs12229654, P = 4.56 × 10(-9)), ITIH4 (rs2535633, P = 1.77 × 10(-10)) and NT5C2 (rs11191580, P = 3.83 × 10(-8)) genes. The association of BMI with rs2237892, rs671 and rs12229654 was significantly stronger among men than among women. Of the 51 BMI-associated loci initially identified in European-ancestry populations, we confirmed eight loci at the genome-wide significance level (P < 5.0 × 10(-8)) and an additional 14 at P < 1.0 × 10(-3) with the same direction of effect as reported previously. Findings from this analysis expand our knowledge of the genetic basis of obesity.

Recent genetic association studies have identified 55 genetic loci associated with obesity or body mass index (BMI). The vast majority, 51 loci, however, were identified in European-ancestry populations. We conducted a meta-analysis of associations between BMI and ∼2.5 million genotyped or imputed single nucleotide polymorphisms among 86 757 individuals of Asian ancestry, followed by in silico and de novo replication among 7488-47 352 additional Asian-ancestry individuals. We identified four novel BMI-associated loci near the KCNQ1 (rs2237892, P = 9.29 × 10(-13)), ALDH2/MYL2 (rs671, P = 3.40 × 10(-11); rs12229654, P = 4.56 × 10(-9)), ITIH4 (rs2535633, P = 1.77 × 10(-10)) and NT5C2 (rs11191580, P = 3.83 × 10(-8)) genes. The association of BMI with rs2237892, rs671 and rs12229654 was significantly stronger among men than among women. Of the 51 BMI-associated loci initially identified in European-ancestry populations, we confirmed eight loci at the genome-wide significance level (P < 5.0 × 10(-8)) and an additional 14 at P < 1.0 × 10(-3) with the same direction of effect as reported previously. Findings from this analysis expand our knowledge of the genetic basis of obesity.

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.