Glucocorticoid signaling in pancreatic islets modulates gene regulatory programs and genetic risk of type 2 diabetes

Glucocorticoids are key regulators of glucose homeostasis and pancreatic islet function, but the gene regulatory programs driving responses to glucocorticoid signaling in islets and the contribution of these programs to diabetes risk are unknown. In this study we used ATAC-seq and RNA-seq to map chromatin accessibility and gene expression from eleven primary human islet samples cultured in vitro with the glucocorticoid dexamethasone at multiple doses and durations. We identified thousands of accessible chromatin sites and genes with significant changes in activity in response to glucocorticoids. Chromatin sites up-regulated in glucocorticoid signaling were prominently enriched for glucocorticoid receptor binding sites and up-regulated genes were enriched for ion transport and lipid metabolism, whereas down-regulated chromatin sites and genes were enriched for inflammatory, stress response and proliferative processes. Genetic variants associated with glucose levels and T2D risk were enriched in glucocorticoid-responsive chromatin sites, including fine-mapped variants at 51 known signals. Among fine-mapped variants in glucocorticoid-responsive chromatin, a likely casual variant at the 2p21 locus had glucocorticoid-dependent allelic effects on beta cell enhancer activity and affected SIX2 and SIX3 expression. Our results provide a comprehensive map of islet regulatory programs in response to glucocorticoids through which we uncover a role for islet glucocorticoid signaling in mediating genetic risk of T2D.     

Trans-Ethnic Fine-Mapping of Lipid Loci Identifies Population-Specific Signals and Allelic Heterogeneity That Increases the Trait Variance Explained

Genome-wide association studies (GWAS) have identified ∼100 loci associated with blood lipid levels, but much of the trait heritability remains unexplained, and at most loci the identities of the trait-influencing variants remain unknown. We conducted a trans-ethnic fine-mapping study at 18, 22, and 18 GWAS loci on the Metabochip for their association with triglycerides (TG), high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C), respectively, in individuals of African American (n = 6,832), East Asian (n = 9,449), and European (n = 10,829) ancestry. We aimed to identify the variants with strongest association at each locus, identify additional and population-specific signals, refine association signals, and assess the relative significance of previously described functional variants. Among the 58 loci, 33 exhibited evidence of association at P<1×10⁻⁴ in at least one ancestry group. Sequential conditional analyses revealed that ten, nine, and four loci in African Americans, Europeans, and East Asians, respectively, exhibited two or more signals. At these loci, accounting for all signals led to a 1.3- to 1.8-fold increase in the explained phenotypic variance compared to the strongest signals. Distinct signals across ancestry groups were identified at PCSK9 and APOA5. Trans-ethnic analyses narrowed the signals to smaller sets of variants at GCKR, PPP1R3B, ABO, LCAT, and ABCA1. Of 27 variants reported previously to have functional effects, 74% exhibited the strongest association at the respective signal. In conclusion, trans-ethnic high-density genotyping and analysis confirm the presence of allelic heterogeneity, allow the identification of population-specific variants, and limit the number of candidate SNPs for functional studies.     

Combined analysis of genome-wide association studies for Crohn disease and psoriasis identifies seven shared susceptibility loci

Psoriasis (PS) and Crohn disease (CD) have been shown to be epidemiologically, pathologically, and therapeutically connected, but little is known about their shared genetic causes. We performed meta-analyses of five published genome-wide association studies on PS (2,529 cases and 4,955 controls) and CD (2,142 cases and 5,505 controls), followed up 20 loci that showed strongest evidence for shared disease association and, furthermore, tested cross-disease associations for previously reported PS and CD risk alleles in additional 6,115 PS cases, 4,073 CD cases, and 10,100 controls. We identified seven susceptibility loci outside the human leukocyte antigen region (9p24 near JAK2, 10q22 at ZMIZ1, 11q13 near PRDX5, 16p13 near SOCS1, 17q21 at STAT3, 19p13 near FUT2, and 22q11 at YDJC) shared between PS and CD with genome-wide significance (p < 5 × 10⁻⁸) and confirmed four already established PS and CD risk loci (IL23R, IL12B, REL, and TYK2). Three of the shared loci are also genome-wide significantly associated with PS alone (10q22 at ZMIZ1, p_rs1250544 = 3.53 × 10⁻⁸, 11q13 near PRDX5, p_rs694739 = 3.71 × 10⁻⁰⁹, 22q11 at YDJC, p_rs181359 = 8.02 × 10⁻¹⁰). In addition, we identified one susceptibility locus for CD (16p13 near SOCS1, p_rs4780355 = 4.99 × 10⁻⁸). Refinement of association signals identified shared genome-wide significant associations for exonic SNPs at 10q22 (ZMIZ1) and in silico expression quantitative trait locus analyses revealed that the associations at ZMIZ1 and near SOCS1 have a potential functional effect on gene expression. Our results show the usefulness of joint analyses of clinically distinct immune-mediated diseases and enlarge the map of shared genetic risk loci.     

Replication Study in a Japanese Population to Evaluate the Association between 10 SNP Loci,Identified in European Genome-Wide Association Studies,and Type 2 Diabetes

AimWe performed a replication study in a Japanese population to evaluate the association between type 2 diabetes and 7 susceptibility loci originally identified by European genome-wide association study (GWAS) in 2012: ZMIZ1, KLHDC5, TLE1, ANKRD55, CILP2, MC4R, and BCAR1. We also examined the association of 3 additional loci: CCND2 and GIPR, identified in sex-differentiated analyses, and LAMA1, which was shown to be associated with non-obese European type 2 diabetes.MethodsWe genotyped 6,972 Japanese participants (4,280 type 2 diabetes patients and 2,692 controls) for each of the 10 single nucleotide polymorphisms (SNPs): rs12571751 in ZMIZ1, rs10842994 near KLHDC5, rs2796441 near TLE1, rs459193 near ANKRD55, rs10401969 in CILP2, rs12970134 near MC4R, rs7202877 near BCAR1, rs11063069 near CCND2, rs8108269 near GIPR, and rs8090011 in LAMA1 using a multiplex polymerase chain reaction invader assay. The association of each SNP locus with the disease was evaluated using a logistic regression analysis.ResultsAll SNPs examined in this study had the same direction of effect (odds ratio > 1.0, p = 9.77 × 10^-4, binomial test), as in the original reports. Among them, rs12571751 in ZMIZ1 was significantly associated with type 2 diabetes [p = 0.0041, odds ratio = 1.123, 95% confidence interval 1.037–1.215, adjusted for sex, age and body mass index (BMI)], but we did not observe significant association of the remaining 9 SNP loci with type 2 diabetes in the present Japanese population (p ≥ 0.005). A genetic risk score, constructed from the sum of risk alleles for the 7 SNP loci identified by un-stratified analyses in the European GWAS meta-analysis were associated with type 2 diabetes in the present Japanese population (p = 2.3 × 10^-4, adjusted for sex, age and BMI).ConclusionsZMIZ1 locus has a significant effect on conferring susceptibility to type 2 diabetes also in the Japanese population.     

Inferred expression regulator activities suggest genes mediating cardiometabolic genetic signals

Expression QTL (eQTL) analyses have suggested many genes mediating genome-wide association study (GWAS) signals but most GWAS signals still lack compelling explanatory genes. We have leveraged an adipose-specific gene regulatory network to infer expression regulator activities and phenotypic master regulators (MRs), which were used to detect activity QTLs (aQTLs) at cardiometabolic trait GWAS loci. Regulator activities were inferred with the VIPER algorithm that integrates enrichment of expected expression changes among a regulator’s target genes with confidence in their regulator-target network interactions and target overlap between different regulators (i.e., pleiotropy). Phenotypic MRs were identified as those regulators whose activities were most important in predicting their respective phenotypes using random forest modeling. While eQTLs were typically more significant than aQTLs in cis, the opposite was true among candidate MRs in trans. Several GWAS loci colocalized with MR trans-eQTLs/aQTLs in the absence of colocalized cis-QTLs. Intriguingly, at the 1p36.1 BMI GWAS locus the EPHB2 cis-aQTL was stronger than its cis-eQTL and colocalized with the GWAS signal and 35 BMI MR trans-aQTLs, suggesting the GWAS signal may be mediated by effects on EPHB2 activity and its downstream effects on a network of BMI MRs. These MR and aQTL analyses represent systems genetic methods that may be broadly applied to supplement standard eQTL analyses for suggesting molecular effects mediating GWAS signals.     

IL2RA Genetic Heterogeneity in Multiple Sclerosis and Type 1 Diabetes Susceptibility and Soluble Interleukin-2 Receptor Production

Multiple sclerosis (MS) and type 1 diabetes (T1D) are organ-specific autoimmune disorders with significant heritability, part of which is conferred by shared alleles. For decades, the Human Leukocyte Antigen (HLA) complex was the only known susceptibility locus for both T1D and MS, but loci outside the HLA complex harboring risk alleles have been discovered and fully replicated. A genome-wide association scan for MS risk genes and candidate gene association studies have previously described the IL2RA gene region as a shared autoimmune locus. In order to investigate whether autoimmunity risk at IL2RA was due to distinct or shared alleles, we performed a genetic association study of three IL2RA variants in a DNA collection of up to 9,407 healthy controls, 2,420 MS, and 6,425 T1D subjects as well as 1,303 MS parent/child trios. Here, we report “allelic heterogeneity” at the IL2RA region between MS and T1D. We observe an allele associated with susceptibility to one disease and risk to the other, an allele that confers susceptibility to both diseases, and an allele that may only confer susceptibility to T1D. In addition, we tested the levels of soluble interleukin-2 receptor (sIL-2RA) in the serum from up to 69 healthy control subjects, 285 MS, and 1,317 T1D subjects. We demonstrate that multiple variants independently correlate with sIL-2RA levels.     

Meta-Analysis of Genome-Wide Association Studies in African Americans Provides Insights into the Genetic Architecture of Type 2 Diabetes

Type 2 diabetes (T2D) is more prevalent in African Americans than in Europeans. However, little is known about the genetic risk in African Americans despite the recent identification of more than 70 T2D loci primarily by genome-wide association studies (GWAS) in individuals of European ancestry. In order to investigate the genetic architecture of T2D in African Americans, the MEta-analysis of type 2 DIabetes in African Americans (MEDIA) Consortium examined 17 GWAS on T2D comprising 8,284 cases and 15,543 controls in African Americans in stage 1 analysis. Single nucleotide polymorphisms (SNPs) association analysis was conducted in each study under the additive model after adjustment for age, sex, study site, and principal components. Meta-analysis of approximately 2.6 million genotyped and imputed SNPs in all studies was conducted using an inverse variance-weighted fixed effect model. Replications were performed to follow up 21 loci in up to 6,061 cases and 5,483 controls in African Americans, and 8,130 cases and 38,987 controls of European ancestry. We identified three known loci (TCF7L2, HMGA2 and KCNQ1) and two novel loci (HLA-B and INS-IGF2) at genome-wide significance (4.15×10⁻⁹⁴<P<5×10⁻⁸, odds ratio (OR) = 1.09 to 1.36). Fine-mapping revealed that 88 of 158 previously identified T2D or glucose homeostasis loci demonstrated nominal to highly significant association (2.2×10⁻²³ < locus-wide P<0.05). These novel and previously identified loci yielded a sibling relative risk of 1.19, explaining 17.5% of the phenotypic variance of T2D on the liability scale in African Americans. Overall, this study identified two novel susceptibility loci for T2D in African Americans. A substantial number of previously reported loci are transferable to African Americans after accounting for linkage disequilibrium, enabling fine mapping of causal variants in trans-ethnic meta-analysis studies.     

ARLTS1 and Prostate Cancer Risk - Analysis of Expression and Regulation

Prostate cancer (PCa) is a heterogeneous trait for which several susceptibility loci have been implicated by genome-wide linkage and association studies. The genomic region 13q14 is frequently deleted in tumour tissues of both sporadic and familial PCa patients and is consequently recognised as a possible locus of tumour suppressor gene(s). Deletions of this region have been found in many other cancers. Recently, we showed that homozygous carriers for the T442C variant of the ARLTS1 gene (ADP-ribosylation factor-like tumour suppressor protein 1 or ARL11, located at 13q14) are associated with an increased risk for both unselected and familial PCa. Furthermore, the variant T442C was observed in greater frequency among malignant tissue samples, PCa cell lines and xenografts, supporting its role in PCa tumourigenesis. In this study, 84 PCa cases and 15 controls were analysed for ARLTS1 expression status in blood-derived RNA. A statistically significant (p = 0.0037) decrease of ARLTS1 expression in PCa cases was detected. Regulation of ARLTS1 expression was analysed with eQTL (expression quantitative trait loci) methods. Altogether fourteen significant cis-eQTLs affecting the ARLTS1 expression level were found. In addition, epistatic interactions of ARLTS1 genomic variants with genes involved in immune system processes were predicted with the MDR program. In conclusion, this study further supports the role of ARLTS1 as a tumour suppressor gene and reveals that the expression is regulated through variants localised in regulatory regions.     

Epistasis between hyperglycemic QTLs revealed in a double congenic of the OLETF rat

Glucose homeostasis is believed to be regulated by multiple genetic components, in addition to numerous external factors. It is therefore crucial to dissect and understand what roles each causative gene plays in maintaining proper glucose metabolism. In OLETF (Otsuka Long-Evans Tokushima Fatty) rat, a model of polygenic type 2 diabetes, at least 14 quantitative trait loci (QTLs) influencing plasma glucose levels were identified. In congenic strains some of the OLETF allelic variants were shown to increase glucose levels. In this study the focus was on two of the hyperglycemic loci, Nidd1/of and Nidd2/of. Congenic rats possessing OLETF genome fragment at either locus showed similar levels of mild hyperglycemia. A newly established double congenic rat showed a further aggravation of hyperglycemia. The Nidd1/of locus was also shown to function in the reduction of plasma leptin levels and fat weights, while the Nidd2/of locus led to increased plasma insulin and fat weights. Interestingly, both plasma leptin and fat weights reverted to the control levels in the double congenic rat. These results indicate that there is an epistatic interaction between the two loci. However, it is unlikely that the abnormal level of enhanced glucose homeostasis is mediated, at least not directly, by leptin or fat mass.     

A Genome-Wide Assessment of the Role of Untagged Copy Number Variants in Type 1 Diabetes

Genome-wide association studies (GWAS) for type 1 diabetes (T1D) have successfully identified more than 40 independent T1D associated tagging single nucleotide polymorphisms (SNPs). However, owing to technical limitations of copy number variants (CNVs) genotyping assays, the assessment of the role of CNVs has been limited to the subset of these in high linkage disequilibrium with tag SNPs. The contribution of untagged CNVs, often multi-allelic and difficult to genotype using existing assays, to the heritability of T1D remains an open question. To investigate this issue, we designed a custom comparative genetic hybridization array (aCGH) specifically designed to assay untagged CNV loci identified from a variety of sources. To overcome the technical limitations of the case control design for this class of CNVs, we genotyped the Type 1 Diabetes Genetics Consortium (T1DGC) family resource (representing 3,903 transmissions from parents to affected offspring) and used an association testing strategy that does not necessitate obtaining discrete genotypes. Our design targeted 4,309 CNVs, of which 3,410 passed stringent quality control filters. As a positive control, the scan confirmed the known T1D association at the INS locus by direct typing of the 5′ variable number of tandem repeat (VNTR) locus. Our results clarify the fact that the disease association is indistinguishable from the two main polymorphic allele classes of the INS VNTR, class I-and class III. We also identified novel technical artifacts resulting into spurious associations at the somatically rearranging loci, T cell receptor, TCRA/TCRD and TCRB, and Immunoglobulin heavy chain, IGH, loci on chromosomes 14q11.2, 7q34 and 14q32.33, respectively. However, our data did not identify novel T1D loci. Our results do not support a major role of untagged CNVs in T1D heritability.     

The Multi-allelic Genetic Architecture of a Variance-Heterogeneity Locus for Molybdenum Concentration in Leaves Acts as a Source of Unexplained Additive Genetic Variance

Genome-wide association (GWA) analyses have generally been used to detect individual loci contributing to the phenotypic diversity in a population by the effects of these loci on the trait mean. More rarely, loci have also been detected based on variance differences between genotypes. Several hypotheses have been proposed to explain the possible genetic mechanisms leading to such variance signals. However, little is known about what causes these signals, or whether this genetic variance-heterogeneity reflects mechanisms of importance in natural populations. Previously, we identified a variance-heterogeneity GWA (vGWA) signal for leaf molybdenum concentrations in Arabidopsis thaliana. Here, fine-mapping of this association reveals that the vGWA emerges from the effects of three independent genetic polymorphisms that all are in strong LD with the markers displaying the genetic variance-heterogeneity. By revealing the genetic architecture underlying this vGWA signal, we uncovered the molecular source of a significant amount of hidden additive genetic variation or “missing heritability”. Two of the three polymorphisms underlying the genetic variance-heterogeneity are promoter variants for Molybdate transporter 1 (MOT1), and the third a variant located ~25 kb downstream of this gene. A fourth independent association was also detected ~600 kb upstream of MOT1. Use of a T-DNA knockout allele highlights Copper Transporter 6; COPT6 (AT2G26975) as a strong candidate gene for this association. Our results show that an extended LD across a complex locus including multiple functional alleles can lead to a variance-heterogeneity between genotypes in natural populations. Further, they provide novel insights into the genetic regulation of ion homeostasis in A. thaliana, and empirically confirm that variance-heterogeneity based GWA methods are a valuable tool to detect novel associations of biological importance in natural populations.     

Large-scale gene-centric meta-analysis across 39 studies identifies type 2 diabetes loci

To identify genetic factors contributing to type 2 diabetes (T2D), we performed large-scale meta-analyses by using a custom ～50,000 SNP genotyping array (the ITMAT-Broad-CARe array) with ～2000 candidate genes in 39 multiethnic population-based studies, case-control studies, and clinical trials totaling 17,418 cases and 70,298 controls. First, meta-analysis of 25 studies comprising 14,073 cases and 57,489 controls of European descent confirmed eight established T2D loci at genome-wide significance. In silico follow-up analysis of putative association signals found in independent genome-wide association studies (including 8,130 cases and 38,987 controls) performed by the DIAGRAM consortium identified a T2D locus at genome-wide significance (GATAD2A/CILP2/PBX4; p = 5.7 × 10(-9)) and two loci exceeding study-wide significance (SREBF1, and TH/INS; p < 2.4 × 10(-6)). Second, meta-analyses of 1,986 cases and 7,695 controls from eight African-American studies identified study-wide-significant (p = 2.4 × 10(-7)) variants in HMGA2 and replicated variants in TCF7L2 (p = 5.1 × 10(-15)). Third, conditional analysis revealed multiple known and novel independent signals within five T2D-associated genes in samples of European ancestry and within HMGA2 in African-American samples. Fourth, a multiethnic meta-analysis of all 39 studies identified T2D-associated variants in BCL2 (p = 2.1 × 10(-8)). Finally, a composite genetic score of SNPs from new and established T2D signals was significantly associated with increased risk of diabetes in African-American, Hispanic, and Asian populations. In summary, large-scale meta-analysis involving a dense gene-centric approach has uncovered additional loci and variants that contribute to T2D risk and suggests substantial overlap of T2D association signals across multiple ethnic groups.     

Effects of Genetic Variants Previously Associated with Fasting Glucose and Insulin in the Diabetes Prevention Program

Common genetic variants have been recently associated with fasting glucose and insulin levels in white populations. Whether these associations replicate in pre-diabetes is not known. We extended these findings to the Diabetes Prevention Program, a clinical trial in which participants at high risk for diabetes were randomized to placebo, lifestyle modification or metformin for diabetes prevention. We genotyped previously reported polymorphisms (or their proxies) in/near G6PC2, MTNR1B, GCK, DGKB, GCKR, ADCY5, MADD, CRY2, ADRA2A, FADS1, PROX1, SLC2A2, GLIS3, C2CD4B, IGF1, and IRS1 in 3,548 Diabetes Prevention Program participants. We analyzed variants for association with baseline glycemic traits, incident diabetes and their interaction with response to metformin or lifestyle intervention. We replicated associations with fasting glucose at MTNR1B (P<0.001), G6PC2 (P = 0.002) and GCKR (P = 0.001). We noted impaired β-cell function in carriers of glucose-raising alleles at MTNR1B (P<0.001), and an increase in the insulinogenic index for the glucose-raising allele at G6PC2 (P<0.001). The association of MTNR1B with fasting glucose and impaired β-cell function persisted at 1 year despite adjustment for the baseline trait, indicating a sustained deleterious effect at this locus. We also replicated the association of MADD with fasting proinsulin levels (P<0.001). We detected no significant impact of these variants on diabetes incidence or interaction with preventive interventions. The association of several polymorphisms with quantitative glycemic traits is replicated in a cohort of high-risk persons. These variants do not have a detectable impact on diabetes incidence or response to metformin or lifestyle modification in the Diabetes Prevention Program.     

Replication of Type 2 Diabetes Candidate Genes Variations in Three Geographically Unrelated Indian Population Groups

Type 2 diabetes (T2D) is a syndrome of multiple metabolic disorders and is genetically heterogeneous. India comprises one of the largest global populations with highest number of reported type 2 diabetes cases. However, limited information about T2D associated loci is available for Indian populations. It is, therefore, pertinent to evaluate the previously associated candidates as well as identify novel genetic variations in Indian populations to understand the extent of genetic heterogeneity. We chose to do a cost effective high-throughput mass-array genotyping and studied the candidate gene variations associated with T2D in literature. In this case-control candidate genes association study, 91 SNPs from 55 candidate genes have been analyzed in three geographically independent population groups from India. We report the genetic variants in five candidate genes: TCF7L2, HHEX, ENPP1, IDE and FTO, are significantly associated (after Bonferroni correction, p<5.5E−04) with T2D susceptibility in combined population. Interestingly, SNP rs7903146 of the TCF7L2 gene passed the genome wide significance threshold (combined P value = 2.05E−08) in the studied populations. We also observed the association of rs7903146 with blood glucose (fasting and postprandial) levels, supporting the role of TCF7L2 gene in blood glucose homeostasis. Further, we noted that the moderate risk provided by the independently associated loci in combined population with Odds Ratio (OR)<1.38 increased to OR = 2.44, (95%CI = 1.67–3.59) when the risk providing genotypes of TCF7L2, HHEX, ENPP1 and FTO genes were combined, suggesting the importance of gene-gene interactions evaluation in complex disorders like T2D.     

Genetic variability in sodium-glucose cotransporter 2 influences glycemic control and risk for diabetic retinopathy in type 2 diabetes patients

BackgroundGluconeogenesis and renal glucose excretion in kidneys both play an important role in glucose homeostasis. Sodium-glucose cotransporter (SGLT2), coded by the SLC5A2 gene is responsible for reabsorption up to 99% of the filtered glucose in proximal tubules. SLC5A2 genetic polymorphisms were suggested to influence glucose homeostasis. We investigated if common SLC5A2 rs9934336 polymorphism influences glycemic control and risk for macro or microvascular complications in Slovenian type 2 diabetes (T2D) patients.MethodsAll 181 clinically well characterized T2D patients were genotyped for SLC5A2 rs9934336 G>A polymorphism. Associations with glycemic control and T2D complications were assessed with nonparametric tests and logistic regression.ResultsSLC5A2 rs9934336 was significantly associated with increased fasting blood glucose levels (P<0.001) and HbA1c levels under the dominant genetic model (P=0.030). After adjustment for T2D duration, significantly higher risk for diabetic retinopathy was present in carriers of at least one polymorphic SLC5A2 rs9934336 A allele compared to non-carriers (OR=7.62; 95%CI=1.65-35.28; P=0.009).ConclusionsOur pilot study suggests an important role of SLC5A2 polymorphisms in the physiologic process of glucose reabsorption in kidneys in T2D patients. This is also the first report on the association between SLC5A2 polymorphism and diabetic retinopathy.     

An Exploratory Study of the Association between KCNB1 rs1051295 and Type 2 Diabetes and Its Related Traits in Chinese Han Population

Since the KCNB1 encoding Kv2.1 channel accounts for the majority of Kv currents modulating insulin secretion by pancreatic islet beta-cells, we postulated that KCNB1 is a plausible candidate gene for genetic variation contributing to the variable compensatory secretory function of beta-cells in type-2 diabetes (T2D). We conducted two studies, a case-control study and a cross-section study, to investigate the association of common single-nucleotide polymorphisms (SNPs) in KCNB1 with T2D and its linking traits. In the case-control study, we first examined the association of 20 tag SNPs of KCNB1 with T2D in a population with 226 T2D patients and non-diabetic subjects (screening study). We then identified the association in an enlarged population of 412 T2D patients and non-diabetic subjects (replication study). In the cross-sectional study, we investigated the linkage between the candidate SNP rs1051295 and T2D by comparing beta-cell function and insulin sensitivity among rs1051295 genotypes in a general population of 1051 subjects at fasting and after glucose loading (oral glucose tolerance tests, OGTT) in 84 fasting glucose impaired subjects, and several T2D-related traits. We found that among the 19 available tag SNPs, only the KCNB1 rs1051295 was associated with T2D (P = 0.027), with the rs1051295 TT genotype associated with an increased risk of T2D compared with genotypes CC (P = 0.009). At fasting, rs1051295 genotype TT was associated with a 9.8% reduction in insulin sensitivity compared to CC (P = 0.008); along with increased plasma triglycerides (TG) levels (TT/CC: P = 0.046) and increased waist/hip (W/H) ratio (TT/CC: P = 0.013; TT/TC: P = 0.002). OGTT confirmed that genotype TT exhibited reduced insulin sensitivity by 16.3% (P = 0.030) compared with genotype TC+CC in a fasting glucose impaired population. The KCNB1 rs1051295 genotype TT in the Chinese Han population is associated with decreased insulin sensitivity and increased plasma TG and W/H ratio, which together contribute to an increased risk for T2D.     

Cell Specific eQTL Analysis without Sorting Cells

The functional consequences of trait associated SNPs are often investigated using expression quantitative trait locus (eQTL) mapping. While trait-associated variants may operate in a cell-type specific manner, eQTL datasets for such cell-types may not always be available. We performed a genome-environment interaction (GxE) meta-analysis on data from 5,683 samples to infer the cell type specificity of whole blood cis-eQTLs. We demonstrate that this method is able to predict neutrophil and lymphocyte specific cis-eQTLs and replicate these predictions in independent cell-type specific datasets. Finally, we show that SNPs associated with Crohn’s disease preferentially affect gene expression within neutrophils, including the archetypal NOD2 locus.