Type 2 Diabetes Variants Disrupt Function of SLC16A11 through Two Distinct Mechanisms.

Type 2 diabetes (T2D) affects Latinos at twice the rate seen in populations of European descent. We recently identified a risk haplotype spanning SLC16A11 that explains ∼20% of the increased T2D prevalence in Mexico. Here, through genetic fine-mapping, we define a set of tightly linked variants likely to contain the causal allele(s). We show that variants on the T2D-associated haplotype have two distinct effects: (1) decreasing SLC16A11 expression in liver and (2) disrupting a key interaction with basigin, thereby reducing cell-surface localization. Both independent mechanisms reduce SLC16A11 function and suggest SLC16A11 is the causal gene at this locus. To gain insight into how SLC16A11 disruption impacts T2D risk, we demonstrate that SLC16A11 is a proton-coupled monocarboxylate transporter and that genetic perturbation of SLC16A11 induces changes in fatty acid and lipid metabolism that are associated with increased T2D risk. Our findings suggest that increasing SLC16A11 function could be therapeutically beneficial for T2D. VIDEO ABSTRACT.

Identification of type 2 diabetes loci in 433,540 East Asian individuals.

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 cells3. At another locus, expression quantitative trait loci at two overlapping T2D signals affect two genes-NKX6-3 and ANK1-in different tissues4-6. Association studies in diverse populations identify additional loci and elucidate disease-associated genes, biology, and pathways.

Methods for estimating insulin resistance from untargeted metabolomics data.

CONTEXT: Insulin resistance is associated with multiple complex diseases; however, precise measures of insulin resistance are invasive, expensive, and time-consuming. OBJECTIVE: Develop estimation models for measures of insulin resistance, including insulin sensitivity index (SI) and homeostatic model assessment of insulin resistance (HOMA-IR) from metabolomics data. DESIGN: Insulin Resistance Atherosclerosis Family Study (IRASFS). SETTING: Community based. PARTICIPANTS: Mexican Americans (MA) and African Americans (AA). MAIN OUTCOME: Estimation models for measures of insulin resistance, i.e. SI and HOMA-IR. RESULTS: Least Absolute Shrinkage and Selection Operator (LASSO) and Elastic Net regression were used to build insulin resistance estimation models from 1274 metabolites combined with clinical data, e.g. age, sex, body mass index (BMI). Metabolite data were transformed using three approaches, i.e. inverse normal transformation, standardization, and Box Cox transformation. The analysis was performed in one MA recruitment site (San Luis Valley, Colorado (SLV); N = 450) and tested in another MA recruitment site (San Antonio, Texas (SA); N = 473). In addition, the two MA recruitment sites were combined and estimation models tested in the AA recruitment sample (Los Angeles, California; N = 495). Estimated and empiric SI were correlated in the SA (r2 = 0.77) and AA (r2 = 0.74) testing datasets. Further, estimated and empiric SI were consistently associated with BMI, low-density lipoprotein cholesterol (LDL), and triglycerides. We applied similar approaches to estimate HOMA-IR with similar results. CONCLUSIONS: We have developed a method for estimating insulin resistance with metabolomics data that has the potential for application to a wide range of biomedical studies and conditions.

Metabolomic profiling of glucose homeostasis in African Americans: the Insulin Resistance Atherosclerosis Family Study (IRAS-FS).

INTRODUCTION: African Americans are at increased risk for type 2 diabetes. OBJECTIVES: This work aimed to examine metabolomic signature of glucose homeostasis in African Americans. METHODS: We used an untargeted liquid chromatography-mass spectrometry metabolomic approach to comprehensively profile 727 plasma metabolites among 571 African Americans from the Insulin Resistance Atherosclerosis Family Study (IRAS-FS) and investigate the associations between these metabolites and both the dynamic (SI, insulin sensitivity; AIR, acute insulin response; DI, disposition index; and SG, glucose effectiveness) and basal (HOMA-IR and HOMA-B) measures of glucose homeostasis using univariate and regularized regression models. We also compared the results with our previous findings in the IRAS-FS Mexican Americans. RESULTS: We confirmed increased plasma metabolite levels of branched-chain amino acids and their metabolic derivatives, 2-aminoadipate, 2-hydroxybutyrate, glutamate, arginine and its metabolic derivatives, carbohydrate metabolites, and medium- and long-chain fatty acids were associated with insulin resistance, while increased plasma metabolite levels in the glycine, serine and threonine metabolic pathway were associated with insulin sensitivity. We also observed a differential ancestral effect of glutamate on glucose homeostasis with significantly stronger effects observed in African Americans than those previously observed in Mexican Americans. CONCLUSION: We extended the observations that metabolites are useful biomarkers in the identification of prediabetes in individuals at risk of type 2 diabetes in African Americans. We revealed, for the first time, differential ancestral effect of certain metabolites (i.e., glutamate) on glucose homeostasis traits. Our study highlights the need for additional comprehensive metabolomic studies in well-characterized multiethnic cohorts.

Rare and low-frequency coding variants alter human adult height.

Height is a highly heritable, classic polygenic trait with approximately 700 common associated variants identified through genome-wide association studies so far. Here, we report 83 height-associated coding variants with lower minor-allele frequencies (in the range of 0.1-4.8%) and effects of up to 2 centimetres per allele (such as those in IHH, STC2, AR and CRISPLD2), greater than ten times the average effect of common variants. In functional follow-up studies, rare height-increasing alleles of STC2 (giving an increase of 1-2 centimetres per allele) compromised proteolytic inhibition of PAPP-A and increased cleavage of IGFBP-4 in vitro, resulting in higher bioavailability of insulin-like growth factors. These 83 height-associated variants overlap genes that are mutated in monogenic growth disorders and highlight new biological candidates (such as ADAMTS3, IL11RA and NOX4) and pathways (such as proteoglycan and glycosaminoglycan synthesis) involved in growth. Our results demonstrate that sufficiently large sample sizes can uncover rare and low-frequency variants of moderate-to-large effect associated with polygenic human phenotypes, and that these variants implicate relevant genes and pathways.

Genome-wide association study identifies novel loci for type 2 diabetes-attributed end-stage kidney disease in African Americans.

BACKGROUND: End-stage kidney disease (ESKD) is a significant public health concern disproportionately affecting African Americans (AAs). Type 2 diabetes (T2D) is the leading cause of ESKD in the USA, and efforts to uncover genetic susceptibility to diabetic kidney disease (DKD) have had limited success. A prior genome-wide association study (GWAS) in AAs with T2D-ESKD was expanded with additional AA cases and controls and genotypes imputed to the higher density 1000 Genomes reference panel. The discovery analysis included 3432 T2D-ESKD cases and 6977 non-diabetic non-nephropathy controls (N = 10,409), followed by a discrimination analysis in 2756 T2D non-nephropathy controls to exclude T2D-associated variants. RESULTS: Six independent variants located in or near RND3/RBM43, SLITRK3, ENPP7, GNG7, and APOL1 achieved genome-wide significant association (P < 5 × 10-8) with T2D-ESKD. Following extension analyses in 1910 non-diabetic ESKD cases and 908 non-diabetic non-nephropathy controls, a meta-analysis of 5342 AA all-cause ESKD cases and 6977 AA non-diabetic non-nephropathy controls revealed an additional novel all-cause ESKD locus at EFNB2 (rs77113398; P = 9.84 × 10-9; OR = 1.94). Exclusion of APOL1 renal-risk genotype carriers identified two additional genome-wide significant T2D-ESKD-associated loci at GRAMD3 and MGAT4C. A second variant at GNG7 (rs373971520; P = 2.17 × 10-8, OR = 1.46) remained associated with all-cause ESKD in the APOL1-negative analysis. CONCLUSIONS: Findings provide further evidence for genetic factors associated with advanced kidney disease in AAs with T2D.

Discovery and fine-mapping of adiposity loci using high density imputation of genome-wide association studies in individuals of African ancestry: African Ancestry Anthropometry Genetics Consortium.

Genome-wide association studies (GWAS) have identified >300 loci associated with measures of adiposity including body mass index (BMI) and waist-to-hip ratio (adjusted for BMI, WHRadjBMI), but few have been identified through screening of the African ancestry genomes. We performed large scale meta-analyses and replications in up to 52,895 individuals for BMI and up to 23,095 individuals for WHRadjBMI from the African Ancestry Anthropometry Genetics Consortium (AAAGC) using 1000 Genomes phase 1 imputed GWAS to improve coverage of both common and low frequency variants in the low linkage disequilibrium African ancestry genomes. In the sex-combined analyses, we identified one novel locus (TCF7L2/HABP2) for WHRadjBMI and eight previously established loci at P < 5×10-8: seven for BMI, and one for WHRadjBMI in African ancestry individuals. An additional novel locus (SPRYD7/DLEU2) was identified for WHRadjBMI when combined with European GWAS. In the sex-stratified analyses, we identified three novel loci for BMI (INTS10/LPL and MLC1 in men, IRX4/IRX2 in women) and four for WHRadjBMI (SSX2IP, CASC8, PDE3B and ZDHHC1/HSD11B2 in women) in individuals of African ancestry or both African and European ancestry. For four of the novel variants, the minor allele frequency was low (<5%). In the trans-ethnic fine mapping of 47 BMI loci and 27 WHRadjBMI loci that were locus-wide significant (P < 0.05 adjusted for effective number of variants per locus) from the African ancestry sex-combined and sex-stratified analyses, 26 BMI loci and 17 WHRadjBMI loci contained ≤ 20 variants in the credible sets that jointly account for 99% posterior probability of driving the associations. The lead variants in 13 of these loci had a high probability of being causal. As compared to our previous HapMap imputed GWAS for BMI and WHRadjBMI including up to 71,412 and 27,350 African ancestry individuals, respectively, our results suggest that 1000 Genomes imputation showed modest improvement in identifying GWAS loci including low frequency variants. Trans-ethnic meta-analyses further improved fine mapping of putative causal variants in loci shared between the African and European ancestry populations.

Genome-wide interaction with the insulin secretion locus MTNR1B reveals CMIP as a novel type 2 diabetes susceptibility gene in African Americans.

Although type 2 diabetes (T2D) results from metabolic defects in insulin secretion and insulin sensitivity, most of the genetic risk loci identified to date relates to insulin secretion. We reported that T2D loci influencing insulin sensitivity may be identified through interactions with insulin secretion loci, thereby leading to T2D. Here, we hypothesize that joint testing of variant main effects and interaction effects with an insulin secretion locus increases power to identify genetic interactions leading to T2D. We tested this hypothesis with an intronic MTNR1B SNP, rs10830963, which is associated with acute insulin response to glucose, a dynamic measure of insulin secretion. rs10830963 was tested for interaction and joint (main + interaction) effects with genome-wide data in African Americans (2,452 cases and 3,772 controls) from five cohorts. Genome-wide genotype data (Affymetrix Human Genome 6.0 array) was imputed to a 1000 Genomes Project reference panel. T2D risk was modeled using logistic regression with rs10830963 dosage, age, sex, and principal component as predictors. Joint effects were captured using the Kraft two degrees of freedom test. Genome-wide significant (P < 5 × 10-8 ) interaction with MTNR1B and joint effects were detected for CMIP intronic SNP rs17197883 (Pinteraction  = 1.43 × 10-8 ; Pjoint  = 4.70 × 10-8 ). CMIP variants have been nominally associated with T2D, fasting glucose, and adiponectin in individuals of East Asian ancestry, with high-density lipoprotein, and with waist-to-hip ratio adjusted for body mass index in Europeans. These data support the hypothesis that additional genetic factors contributing to T2D risk, including insulin sensitivity loci, can be identified through interactions with insulin secretion loci.

Trans-ethnic Meta-analysis and Functional Annotation Illuminates the Genetic Architecture of Fasting Glucose and Insulin.

Knowledge of the genetic basis of the type 2 diabetes (T2D)-related quantitative traits fasting glucose (FG) and insulin (FI) in African ancestry (AA) individuals has been limited. In non-diabetic subjects of AA (n = 20,209) and European ancestry (EA; n = 57,292), we performed trans-ethnic (AA+EA) fine-mapping of 54 established EA FG or FI loci with detailed functional annotation, assessed their relevance in AA individuals, and sought previously undescribed loci through trans-ethnic (AA+EA) meta-analysis. We narrowed credible sets of variants driving association signals for 22/54 EA-associated loci; 18/22 credible sets overlapped with active islet-specific enhancers or transcription factor (TF) binding sites, and 21/22 contained at least one TF motif. Of the 54 EA-associated loci, 23 were shared between EA and AA. Replication with an additional 10,096 AA individuals identified two previously undescribed FI loci, chrX FAM133A (rs213676) and chr5 PELO (rs6450057). Trans-ethnic analyses with regulatory annotation illuminate the genetic architecture of glycemic traits and suggest gene regulation as a target to advance precision medicine for T2D. Our approach to utilize state-of-the-art functional annotation and implement trans-ethnic association analysis for discovery and fine-mapping offers a framework for further follow-up and characterization of GWAS signals of complex trait loci.

The African Descent and Glaucoma Evaluation Study (ADAGES) III: Contribution of Genotype to Glaucoma Phenotype in African Americans: Study Design and Baseline Data.

PURPOSE: To describe the study protocol and baseline characteristics of the African Descent and Glaucoma Evaluation Study (ADAGES) III. DESIGN: Cross-sectional, case-control study. PARTICIPANTS: Three thousand two hundred sixty-six glaucoma patients and control participants without glaucoma of African or European descent were recruited from 5 study centers in different regions of the United States. METHODS: Individuals of African descent (AD) and European descent (ED) with primary open-angle glaucoma (POAG) and control participants completed a detailed demographic and medical history interview. Standardized height, weight, and blood pressure measurements were obtained. Saliva and blood samples to provide serum, plasma, DNA, and RNA were collected for standardized processing. Visual fields, stereoscopic disc photographs, and details of the ophthalmic examination were obtained and transferred to the University of California, San Diego, Data Coordinating Center for standardized processing and quality review. MAIN OUTCOME MEASURES: Participant gender, age, race, body mass index, blood pressure, history of smoking and alcohol use in POAG patients and control participants were described. Ophthalmic measures included intraocular pressure, visual field mean deviation, central corneal thickness, glaucoma medication use, or past glaucoma surgery. Ocular conditions, including diabetic retinopathy, age-related macular degeneration, and past cataract surgery, were recorded. RESULTS: The 3266 ADAGES III study participants in this report include 2146 AD POAG patients, 695 ED POAG patients, 198 AD control participants, and 227 ED control participants. The AD POAG patients and control participants were significantly younger (both, 67.4 years) than ED POAG patients and control participants (73.4 and 70.2 years, respectively). After adjusting for age, AD POAG patients had different phenotypic characteristics compared with ED POAG patients, including higher intraocular pressure, worse visual acuity and visual field mean deviation, and thinner corneas (all P < 0.001). Family history of glaucoma did not differ between AD and ED POAG patients. CONCLUSIONS: With its large sample size, extensive specimen collection, and deep phenotyping of AD and ED glaucoma patients and control participants from different regions in the United States, the ADAGES III genomics study will address gaps in our knowledge of the genetics of POAG in this high-risk population.

Genetic Architecture of Primary Open-Angle Glaucoma in Individuals of African Descent: The African Descent and Glaucoma Evaluation Study III.

PURPOSE: To find genetic contributions to glaucoma in African Americans. DESIGN: Cross-sectional, case-control study. PARTICIPANTS: One thousand eight hundred seventy-five primary open-angle glaucoma (POAG) patients and 1709 controls, self-identified as being of African descent (AD), from the African Descent and Glaucoma Evaluation Study (ADAGES) III and Wake Forest School of Medicine. METHODS: MegaChip genotypes were imputed to Thousand Genomes data. Association of single nucleotide polymorphisms (SNPs) with POAG and advanced POAG was tested by linear mixed model correcting for relatedness and population stratification. Genetic risk scores were tested by receiver operator characteristic curves (ROC-AUCs). MAIN OUTCOME MEASURES: Primary open-angle glaucoma defined by visual field loss without other nonocular conditions (n = 1875). Advanced POAG was defined by age-based mean deviation of visual field (n = 946). RESULTS: Eighteen million two hundred eighty-one thousand nine hundred twenty SNPs met imputation quality of r2 > 0.7 and minor allele frequency > 0.005. Association of a novel locus, EN04, was observed for advanced POAG (rs185815146 ß, 0.36; standard error, 0.065; P < 3×10-8). For POAG, an AD signal was observed at the 9p21 European descent (ED) POAG signal (rs79721419; P < 6.5×10-5) independent of the previously observed 9p21 ED signal (rs2383204; P < 2.3×10-5) by conditional analyses. An association with POAG in FNDC3B (rs111698934; P < 3.9×10-5) was observed, not in linkage disequilibrium (LD) with the previously reported ED SNP. Additional previously identified loci associated with POAG in persons of AD were: 8q22, AFAP1, and TMC01. An AUC of 0.62 was observed with an unweighted genetic risk score comprising 11 SNPs in candidate genes. Two additional risk scores were studied by using a penalized matrix decomposition with cross-validation; risk scores of 50 and 400 SNPs were identified with ROC of AUC = 0.74 and AUC = 0.94, respectively. CONCLUSIONS: A novel association with advanced POAG in the EN04 locus was identified putatively in persons of AD. In addition to this finding, this genome-wide association study in POAG patients of AD contributes to POAG genetics by identification of novel signals in prior loci (9p21), as well as advancing the fine mapping of regions because of shorter average LD (FNDC3B). Although not useful without confirmation and clinical trials, the use of genetic risk scores demonstrated that considerable AD-specific genetic information remains in these data.

Genome-wide association studies suggest that APOL1-environment interactions more likely trigger kidney disease in African Americans with nondiabetic nephropathy than strong APOL1-second gene interactions.

African Americans carrying two apolipoprotein L1 gene (APOL1) renal risk variants have a high risk for nephropathy. However, only a minority develops end-stage renal disease (ESRD). Hence, modifying factors likely contribute to initiation of kidney disease such as endogenous (HIV infection) or exogenous (interferon treatment) environmental modifiers. In this report, genome-wide association studies and a meta-analysis were performed to identify novel loci for nondiabetic ESRD in African Americans and to detect genetic modifiers in APOL1-associated nephropathy. Two African American cohorts were analyzed, 1749 nondiabetic ESRD cases and 1136 controls from Wake Forest and 901 lupus nephritis (LN)-ESRD cases and 520 controls with systemic lupus erythematosus but lacking nephropathy from the LN-ESRD Consortium. Association analyses adjusting for APOL1 G1/G2 renal-risk variants were completed and stratified by APOL1 risk genotype status. Individual genome-wide association studies and meta-analysis results of all 2650 ESRD cases and 1656 controls did not detect significant genome-wide associations with ESRD beyond APOL1. Similarly, no single nucleotide polymorphism showed significant genome-wide evidence of an interaction with APOL1 risk variants. Thus, although variants with small individual effects cannot be ruled out and are likely to exist, our results suggest that APOL1-environment interactions may be of greater clinical importance in triggering nephropathy in African Americans than APOL1 interactions with other single nucleotide polymorphisms.

Genome-wide association study of primary open-angle glaucoma in continental and admixed African populations.

Primary open angle glaucoma (POAG) is a complex disease with a major genetic contribution. Its prevalence varies greatly among ethnic groups, and is up to five times more frequent in black African populations compared to Europeans. So far, worldwide efforts to elucidate the genetic complexity of POAG in African populations has been limited. We conducted a genome-wide association study in 1113 POAG cases and 1826 controls from Tanzanian, South African and African American study samples. Apart from confirming evidence of association at TXNRD2 (rs16984299; OR[T] 1.20; P = 0.003), we found that a genetic risk score combining the effects of the 15 previously reported POAG loci was significantly associated with POAG in our samples (OR 1.56; 95% CI 1.26-1.93; P = 4.79 × 10-5). By genome-wide association testing we identified a novel candidate locus, rs141186647, harboring EXOC4 (OR[A] 0.48; P = 3.75 × 10-8), a gene transcribing a component of the exocyst complex involved in vesicle transport. The low frequency and high degree of genetic heterogeneity at this region hampered validation of this finding in predominantly West-African replication sets. Our results suggest that established genetic risk factors play a role in African POAG, however, they do not explain the higher disease load. The high heterogeneity within Africans remains a challenge to identify the genetic commonalities for POAG in this ethnicity, and demands studies of extremely large size.

Analysis of Whole Exome Sequencing with Cardiometabolic Traits Using Family-Based Linkage and Association in the IRAS Family Study.

Family-based methods are a potentially powerful tool to identify trait-defining genetic variants in extended families, particularly when used to complement conventional association analysis. We utilized two-point linkage analysis and single variant association analysis to evaluate whole exome sequencing (WES) data from 1205 Hispanic Americans (78 families) from the Insulin Resistance Atherosclerosis Family Study. WES identified 211,612 variants above the minor allele frequency threshold of ≥0.005. These variants were tested for linkage and/or association with 50 cardiometabolic traits after quality control checks. Two-point linkage analysis yielded 10,580,600 logarithm of the odds (LOD) scores with 1148 LOD scores ≥3, 183 LOD scores ≥4, and 29 LOD scores ≥5. The maximal novel LOD score was 5.50 for rs2289043:T>C, in UNC5C with subcutaneous adipose tissue volume. Association analysis identified 13 variants attaining genome-wide significance (P < 5 × 10-08 ), with the strongest association between rs651821:C>T in APOA5 and triglyceride levels (P = 3.67 × 10-10 ). Overall, there was a 5.2-fold increase in the number of informative variants detected by WES compared to exome chip analysis in this population, nearly 30% of which were novel variants relative to the Database of Single Nucleotide Polymorphisms (dbSNP) build 138. Thus, integration of results from two-point linkage and single-variant association analysis from WES data enabled identification of novel signals potentially contributing to cardiometabolic traits.

Genome-wide linkage and association analysis of cardiometabolic phenotypes in Hispanic Americans.

Linkage studies of complex genetic diseases have been largely replaced by genome-wide association studies, due in part to limited success in complex trait discovery. However, recent interest in rare and low-frequency variants motivates re-examination of family-based methods. In this study, we investigated the performance of two-point linkage analysis for over 1.6 million single-nucleotide polymorphisms (SNPs) combined with single variant association analysis to identify high impact variants, which are both strongly linked and associated with cardiometabolic traits in up to 1414 Hispanics from the Insulin Resistance Atherosclerosis Family Study (IRASFS). Evaluation of all 50 phenotypes yielded 83 557 000 LOD (logarithm of the odds) scores, with 9214 LOD scores ⩾3.0, 845 ⩾4.0 and 89 ⩾5.0, with a maximal LOD score of 6.49 (rs12956744 in the LAMA1 gene for tumor necrosis factor-α (TNFα) receptor 2). Twenty-seven variants were associated with P<0.005 as well as having an LOD score >4, including variants in the NFIB gene under a linkage peak with TNFα receptor 2 levels on chromosome 9. Linkage regions of interest included a broad peak (31 Mb) on chromosome 1q with acute insulin response (max LOD=5.37). This region was previously documented with type 2 diabetes in family-based studies, providing support for the validity of these results. Overall, we have demonstrated the utility of two-point linkage and association in comprehensive genome-wide array-based SNP genotypes.

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(-94)

Variants in KCNQ1 are associated with susceptibility to type 2 diabetes mellitus.

We carried out a multistage genome-wide association study of type 2 diabetes mellitus in Japanese individuals, with a total of 1,612 cases and 1,424 controls and 100,000 SNPs. The most significant association was obtained with SNPs in KCNQ1, and dense mapping within the gene revealed that rs2237892 in intron 15 showed the lowest Pvalue (6.7 x 10(-13), odds ratio (OR) = 1.49). The association of KCNQ1 with type 2 diabetes was replicated in populations of Korean, Chinese and European ancestry as well as in two independent Japanese populations, and meta-analysis with a total of 19,930 individuals (9,569 cases and 10,361 controls) yielded a P value of 1.7 x 10(-42) (OR = 1.40; 95% CI = 1.34-1.47) for rs2237892. Among control subjects, the risk allele of this polymorphism was associated with impairment of insulin secretion according to the homeostasis model assessment of beta-cell function or the corrected insulin response. Our data thus implicate KCNQ1 as a diabetes susceptibility gene in groups of different ancestries.

Genetic and clinical variables identify predictors for chronic kidney disease in type 2 diabetes.

Type 2 diabetes and chronic kidney disease (CKD) may share common risk factors. Here we used a 3-stage procedure to discover novel predictors of CKD by repeatedly applying a stepwise selection based on the Akaike information criterion to subsamples of a prospective complete-case cohort of 2755 patients. This cohort encompassed 25 clinical variables and 36 genetic variants associated with type 2 diabetes, obesity, or fasting plasma glucose. We compared the performance of the clinical, genetic, and clinico-genomic models and used net reclassification improvement to evaluate the impact of top selected genetic variants to the clinico-genomic model. Associations of selected genetic variants with CKD were validated in 2 independent cohorts followed by meta-analyses. Among the top 6 single-nucleotide polymorphisms selected from clinico-genomic data, three (rs478333 of G6PC2, rs7754840 and rs7756992 of CDKAL1) contributed toward the improvement of prediction performance. The variant rs478333 was associated with rapid decline (over 4% per year) in estimated glomerular filtration rate. In a meta-analysis of 2 replication cohorts, the variants rs478333 and rs7754840 showed significant associations with CKD after adjustment for conventional risk factors. Thus, this novel 3-stage approach to a clinico-genomic data set identified 3 novel genetic predictors of CKD in type 2 diabetes. This method can be applied to similar data sets containing clinical and genetic variables to select predictors for clinical outcomes.

The ras responsive transcription factor RREB1 is a novel candidate gene for type 2 diabetes associated end-stage kidney disease.

Familial clustering and presumed genetic risk for type 2 diabetic (T2D) and non-diabetic end-stage kidney disease (ESKD) appear strong in African Americans. Examination of exome sequencing data in African American T2D-ESKD cases and non-diabetic non-nephropathy controls identified two low-frequency variants in the RREB1 gene, a repressor of the angiotensinogen (AGT) gene previously associated with kidney function, as being associated with T2D-ESKD: rs9379084 (P = 0.00087, OR = 0.26; D1171N) and rs41302867 (P = 0.00078, OR = 0.21; splice site variant). Rs41302867 replicated association in an independent sample of African Americans with T2D-ESKD [rs41302867 P = 0.033 (OR = 0.50)], and a trend towards rs9379084 association was observed (P = 0.070). In European Americans with T2D-ESKD compared with European American population based controls, both RREB1 variants replicated association [rs9379084 P = 1.67 × 10(-4) (OR = 0.54) and rs41302867 P = 0.013 (OR = 0.69)]. Rs9379084 was not associated with non-T2D-ESKD or T2D in African Americans (P = 0.55 and P = 0.37, respectively), but was associated with T2D in European Americans (P = 0.014, OR = 0.65). In African Americans, rs41302867 was associated with non-T2D-ESKD [P = 0.036 (OR = 0.54)] and hypertension attributed ESKD [H-ESKD, P = 0.029 (OR = 0.50)]. A meta-analysis combining African American and European American T2D-ESKD data revealed P = 3.52 × 10(-7) and 3.70 × 10(-5) for rs9379084 and rs41302867 association, respectfully. A locus-wide analysis evaluating putatively functional SNPs revealed several nominal associations with T2D-ESKD, non-T2D-ESKD and T2D in African and European Americans. RREB1 is a large, complex gene which codes a multidomain zinc finger binding protein and transcription factor. We posit that variants in RREB1 modulate seemingly disparate phenotypes (i.e. T2D, T2D-ESKD and non-T2D-ESKD) through altered activity resulting from splice site and missense variants.

Association of kidney structure-related gene variants with type 2 diabetes-attributed end-stage kidney disease in African Americans.

African Americans (AAs) are at higher risk for developing end-stage kidney disease (ESKD) compared to European Americans. Genome-wide association studies have identified variants associated with diabetic and non-diabetic kidney diseases. Nephropathy loci, including SLC7A9, UMOD, and SHROOM3, have been implicated in the maintenance of normal glomerular and renal tubular structure and function. Herein, 47 genes important in podocyte, glomerular basement membrane, mesangial cell, mesangial matrix, renal tubular cell, and renal interstitium structure were examined for association with type 2 diabetes (T2D)-attributed ESKD in AAs. Single-variant association analysis was performed in the discovery stage, including 2041 T2D-ESKD cases and 1140 controls (non-diabetic, non-nephropathy). Discrimination analyses in 667 T2D cases-lacking nephropathy excluded T2D-associated SNPs. Nominal associations were tested in an additional 483 T2D-ESKD cases and 554 controls in the replication stage. Meta-analysis of 4218 discovery and replication samples revealed three significant associations with T2D-ESKD at CD2AP and MMP2 (P corr < 0.05 corrected for effective number of SNPs in each locus). Removal of APOL1 renal-risk genotype carriers revealed additional association at five loci, TTC21B, COL4A3, NPHP3-ACAD11, CLDN8, and ARHGAP24 (P corr < 0.05). Genetic variants at COL4A3, CLDN8, and ARHGAP24 were potentially pathogenic. Gene-based associations revealed suggestive significant aggregate effects of coding variants at four genes. Our findings suggest that genetic variation in kidney structure-related genes may contribute to T2D-attributed ESKD in the AA population.