Generation of Isogenic hiPSCs with Targeted Edits at Multiple Intronic SNPs to Study the Effects of the Type 2 Diabetes Associated KCNQ1 Locus in American Indians.

The top genetic association signal for type 2 diabetes (T2D) in Southwestern American Indians maps to intron 15 of KCNQ1, an imprinted gene. We aim to understand the biology whereby variation at this locus affects T2D specifically in this genomic background. To do so, we obtained human induced pluripotent stem cells (hiPSC) derived from American Indians. Using these iPSCs, we show that imprinting of KCNQ1 and CDKN1C during pancreatic islet-like cell generation from iPSCs is consistent with known imprinting patterns in fetal pancreas and adult islets and therefore is an ideal model system to study this locus. In this report, we detail the use of allele-specific guide RNAs and CRISPR to generate isogenic hiPSCs that differ only at multiple T2D associated intronic SNPs at this locus which can be used to elucidate their functional effects. Characterization of these isogenic hiPSCs identified a few aberrant cell lines; namely cell lines with large hemizygous deletions in the putative functional region of KCNQ1 and cell lines hypomethylated at the KCNQ1OT1 promoter. Comparison of an isogenic cell line with a hemizygous deletion to the parental cell line identified CDKN1C and H19 as differentially expressed during the endocrine progenitor stage of pancreatic-islet development.

A missense variant Arg611Cys in LIPE which encodes hormone sensitive lipase decreases lipolysis and increases risk of type 2 diabetes in American Indians.

AIMS: Hormone sensitive lipase (HSL), encoded by the LIPE gene, is involved in lipolysis. Based on prior animal and human studies, LIPE was analysed as a candidate gene for the development of type 2 diabetes (T2D) in a community-based sample of American Indians. MATERIALS AND METHODS: Whole-exome sequence data from 6782 participants with longitudinal clinical measures were used to identify variation in LIPE. RESULTS: Amongst the 16 missense variants identified, an Arg611Cys variant (rs34052647; Cys-allele frequency = 0.087) significantly associated with T2D (OR [95% CI] = 1.38 [1.17-1.64], p = 0.0002, adjusted for age, sex, birth year, and the first five genetic principal components) and an earlier onset age of T2D (HR = 1.22 [1.09-1.36], p = 0.0005). This variant was further analysed for quantitative traits related to T2D. Amongst non-diabetic American Indians, those with the T2D risk Cys-allele had increased insulin levels during an oral glucose tolerance test (0.07 SD per Cys-allele, p = 0.04) and a mixed meal test (0.08 log10 µU/ml per Cys-allele, p = 0.003), and had increased lipid oxidation rates post-absorptively and during insulin infusion (0.07 mg [kg estimated metabolic body size {EMBS}]-1  min-1 per Cys-allele for both, p = 0.01 and 0.009, respectively), compared to individuals with the non-risk Arg-allele. In vitro functional studies showed that cells expressing the Cys-allele had a 17.2% decrease in lipolysis under isoproterenol stimulation (p = 0.03) and a 21.3% decrease in lipase enzyme activity measured by using p-nitrophenyl butyrate as a substrate (p = 0.04) compared to the Arg-allele. CONCLUSION: The Arg611Cys variant causes a modest impairment in lipolysis, thereby affecting glucose homoeostasis and risk of T2D.

Exome Sequencing Identifies A Nonsense Variant in DAO Associated With Reduced Energy Expenditure in American Indians.

BACKGROUND: Obesity and energy expenditure (EE) are heritable and genetic variants influencing EE may contribute to the development of obesity. We sought to identify genetic variants that affect EE in American Indians, an ethnic group with high prevalence of obesity. METHODS: Whole-exome sequencing was performed in 373 healthy Pima Indians informative for 24-hour EE during energy balance. Genetic association analyses of all high-quality exonic variants (≥5 carriers) was performed, and those predicted to be damaging were prioritized. RESULTS: Rs752074397 introduces a premature stop codon (Cys264Ter) in DAO and demonstrated the strongest association for 24-hour EE, where the Ter allele associated with substantially lower 24-hour EE (mean lower by 268 kcal/d) and sleeping EE (by 135 kcal/d). The Ter allele has a frequency = 0.5% in Pima Indians, whereas is extremely rare in most other ethnic groups (frequency < 0.01%). In vitro functional analysis showed reduced protein levels for the truncated form of DAO consistent with increased protein degradation. DAO encodes D-amino acid oxidase, which is involved in dopamine synthesis which might explain its role in modulating EE. CONCLUSION: Our results indicate that a nonsense mutation in DAO may influence EE in American Indians. Identification of variants that influence energy metabolism may lead to new pathways to treat human obesity. CLINICAL TRIAL REGISTRATION NUMBER: NCT00340132.

Using Luciferase Reporter Assays to Identify Functional Variants at Disease-Associated Loci.

The genomic era, highlighted by large scale, genome-wide association studies (GWAS) for both common and rare diseases, have identified hundreds of disease-associated variants. However, most of these variants are not disease causing, but instead only provide information about a potential proximal functional variant through linkage disequilibrium. It is critical that these functional variants be identified, so that their role in disease risk can be ascertained. Luciferase assays are an invaluable tool for identifying and characterizing functional variants, allowing investigations of gene expression, intracellular signaling, transcription factors, receptor activity, and protein folding. In this chapter, we provide an overview of the different ways that luciferase assays can be used to validate functionality of a variant.

Functional and association analysis of an Amerindian-derived population-specific p.(Thr280Met) variant in RBPJL, a component of the PTF1 complex.

PTF1 complex is critical for pancreatic development and maintenance of adult exocrine pancreas. As a part of our ongoing studies to identify genetic variation that contributes to type 2 diabetes (T2D) in American Indians, we analyzed variation in genes that form this complex, namely PTF1A, RBPJ, and its paralogue RBPJL. A c.839C>T (p.(Thr280Met)) variant (rs200998587:C>T, risk allele frequency = 0.03) in RBPJL, identified only in Amerindian-derived populations, associated with T2D (OR = 1.60[1.21-2.13] per Met allele, P = 0.001) and age of diabetes onset (HR = 1.40[1.14-1.72], P = 0.001). Knockdown of Rbpjl in mouse pancreatic acinar cells resulted in a significant decrease in the mRNA expression of genes encoding exocrine enzymes including Ctrb. CTRB1/2 is an established T2D locus where the protective allele associates with increased GLP-1-stimulated insulin secretion and higher expression of CTRB1/2. In vitro studies show that cells expressing the Met280 allele had lower RBPJL protein levels than cells expressing the Thr280 allele, despite having comparable levels of RNA, suggesting that the Met280 RBPJL is less stable. Additionally, luciferase assays in HEK293 cells which examined two different RBPJL responsive promoters, including the promoter for CTRB1, also identified reduced transactivation by the Met280 RBPJL. Similarly, overexpression of both Met280 and Thr280 RBPJL in mouse pancreatic acinar cells identified a significant impairment in the expression of Cel when transactivated by the Met280 RBPJL. In summary, we identified a functional, Amerindian-derived population-specific c.839C>T (p.(Thr280Met)) variant in the pancreas specific RBPJL that may modify T2D risk by regulating exocrine enzyme expression.

Identity-by-Descent Mapping Identifies Major Locus for Serum Triglycerides in Amerindians Largely Explained by an APOC3 Founder Mutation.

BACKGROUND: Identity-by-descent mapping using empirical estimates of identity-by-descent allele sharing may be useful for studies of complex traits in founder populations, where hidden relationships may augment the inherent genetic information that can be used for localization. METHODS AND RESULTS: Through identity-by-descent mapping, using ≈400 000 single-nucleotide polymorphisms (SNPs), of serum lipid profiles, we identified a major linkage signal for triglycerides in 1007 Pima Indians (LOD=9.23; P=3.5×10-11 on chromosome 11q). In subsequent fine-mapping and replication association studies in ≈7500 Amerindians, we determined that this signal reflects effects of a loss-of-function Ala43Thr substitution in APOC3 (rs147210663) and 3 established functional SNPs in APOA5. The association with rs147210663 was particularly strong; each copy of the Thr allele conferred 42% lower triglycerides (ß=-0.92±0.059 SD unit; P=9.6×10-55 in 4668 Pimas and 2793 Southwest Amerindians combined). The Thr allele is extremely rare in most global populations but has a frequency of 2.5% in Pimas. We further demonstrated that 3 APOA5 SNPs with established functional impact could explain the association with the most well-replicated SNP (rs964184) for triglycerides identified by genome-wide association studies. Collectively, these 4 SNPs account for 6.9% of variation in triglycerides in Pimas (and 4.1% in Southwest Amerindians), and their inclusion in the original linkage model reduced the linkage signal to virtually null. CONCLUSIONS: APOC3/APOA5 constitutes a major locus for serum triglycerides in Amerindians, especially the Pimas, and these results provide an empirical example for the concept that population-based linkage analysis is a useful strategy to identify complex trait variants.

Assessment of established HDL-C loci for association with HDL-C levels and type 2 diabetes in Pima Indians.

AIMS/HYPOTHESIS: Epidemiological studies in Pima Indians identified elevated levels of HDL-cholesterol (HDL-C) as a protective factor against type 2 diabetes risk in women. We assessed whether HDL-C-associated single-nucleotide polymorphisms (SNPs) also associate with type 2 diabetes in female Pima Indians. METHODS: Twenty-one SNPs in established HDL-C loci were initially analysed in 2,675 full-heritage Pima Indians. SNPs shown to associate with HDL-C (12 SNPs) were assessed for association with type 2 diabetes in 7,710 Pima Indians (55.6% female sex). The CETP locus provided the strongest evidence for association with HDL-C and was further interrogated by analysing tag SNPs. RESULTS: Twelve of the 21 SNPs analysed had a significant association with HDL-C in Pima Indians; five SNPs representing four loci (CETP, DOCK6, PPP1R3B and ABCA1) reached genome-wide significance. Three SNPs, at CETP, KLF14 and HNF4A, associated with type 2 diabetes only in female participants with the HDL-C-lowering allele increasing diabetes risk (p values: 3.2 × 10(-4) to 7.7 × 10(-5)); the association remained significant even after adjustment for HDL-C. Additional analysis across CETP identified rs6499863 as having the strongest association with type 2 diabetes in female participants (p = 5.0 × 10(-6)) and this association remained independent of the HDL-C association. CONCLUSIONS/INTERPRETATION: SNPs at the CETP, HNF4A and KLF14 locus are associated with HDL-C levels and type 2 diabetes (in female participants). However, since HNF4A and KLF14 are established loci for type 2 diabetes, it is unlikely that HDL-C solely mediates these associations.

Complex Genetics of Type 2 Diabetes and Effect Size: What have We Learned from Isolated Populations?

Genetic studies in large outbred populations have documented a complex, highly polygenic basis for type 2 diabetes (T2D). Most of the variants currently known to be associated with T2D risk have been identified in large studies that included tens of thousands of individuals who are representative of a single major ethnic group such as European, Asian, or African. However, most of these variants have only modest effects on the risk for T2D; identification of definitive 'causal variant' or 'causative loci' is typically lacking. Studies in isolated populations offer several advantages over outbred populations despite being, on average, much smaller in sample size. For example, reduced genetic variability, enrichment of rare variants, and a more uniform environment and lifestyle, which are hallmarks of isolated populations, can reduce the complexity of identifying disease-associated genes. To date, studies in isolated populations have provided valuable insight into the genetic basis of T2D by providing both a deeper understanding of previously identified T2D-associated variants (e.g. demonstrating that variants in KCNQ1 have a strong parent-of-origin effect) or providing novel variants (e.g. ABCC8 in Pima Indians, TBC1D4 in the Greenlandic population, HNF1A in Canadian Oji-Cree). This review summarizes advancements in genetic studies of T2D in outbred and isolated populations, and provides information on whether the difference in the prevalence of T2D in different populations (Pima Indians vs. non-Hispanic Whites and non-Hispanic Whites vs. non-Hispanic Blacks) can be explained by the difference in risk allele frequencies of established T2D variants.

Variants associated with type 2 diabetes identified by the transethnic meta-analysis study: assessment in American Indians and evidence for a new signal in LPP.

AIM/HYPOTHESIS: A recent genome-wide trans-ancestry meta-analysis identified seven new loci associated with type 2 diabetes. We assessed the replication of the seven lead single nucleotide polymorphisms (SNPs) and evaluated these loci for additional signals in American Indians. METHODS: Seven SNPs were genotyped in 7,710 individuals from a longitudinally studied American Indian population, and associations with type 2 diabetes, BMI and related phenotypes were assessed. Previous genome-wide association study (GWAS) data from these individuals were used to screen for additional type 2 diabetes signals at these loci. A variant independent of the trans-ancestry meta-analysis was identified within LPP, and its replication was assessed in an additional 3,106 urban American Indians. RESULTS: SNP rs6813195 near to TMEM154 was nominally associated with type 2 diabetes (p = 0.01, OR 1.12 [95% CI 1.03, 1.22]) and adiposity: the type 2 diabetes risk allele was associated with a lower percentage body fat (ß = -1.451%, p = 4.8 × 10(-4)). Another SNP, rs3130501 near to POU5F1-TCF19, was associated with BMI (ß = -0.012, p = 0.004), type 2 diabetes adjusted for BMI (p = 0.02, OR 1.11 [95% CI 1.02, 1.22]), 2 h glucose concentrations (ß = 0.080 mmol/l, p = 0.02) and insulin resistance estimated by homeostatic model (ß = 0.039, p = 0.009). The independent variant identified at the LPP locus in our American Indian GWAS for type 2 diabetes was replicated in the additional samples (all American Indian meta-analysis, p = 8.9 × 10(-6), OR 1.29 [95% CI 1.15, 1.45]). CONCLUSIONS/INTERPRETATION: For two of the seven newly identified variants, there was nominal evidence for association with type 2 diabetes and related traits in American Indians. Identification of an independent variant at the LPP locus suggests the existence of more than one type 2 diabetes signal at this locus.

Whole exome sequencing identifies variation in CYB5A and RNF10 associated with adiposity and type 2 diabetes.

OBJECTIVE: Few coding variants in genes associated with type 2 diabetes (T2D) have been identified, and the underlying physiologic mechanisms whereby susceptibility genes influence T2D risk are often unknown. The objective of this study was to identify coding variation that increases risk for T2D via an effect on a pre-diabetic trait. METHODS: Whole exome sequencing was done in 177 Pima Indians. Selected variants (N = 345) were genotyped in 555 subjects characterized for body fatness, glucose disposal rates during a clamp, acute insulin response to glucose, and 2-h plasma glucose concentrations during an OGTT, and were also genotyped in up to 5,880 subjects with longitudinal measures of BMI. Variants associated with quantitative traits were assessed for association with T2D in 7,667 subjects. RESULTS: rs7238987 in CYB5A associated with body fatness (P = 7.0 × 10(-6) ). This SNP and a novel SNP in RNF10 also associated with maximum recorded BMI (P = 6.2 × 10(-7) and P = 7.2 × 10(-4) ) and maximum childhood BMI z-score (P = 5.9 × 10(-4) and P = 8.5 × 10(-7) ). The BMI increasing alleles increased the risk for T2D (P = 0.01; OR = 1.13 [1.03-1.24] and 9.5 × 10(-3) ; OR = 1.49 [1.10-2.02]). CONCLUSIONS: CYB5A, which has a role in stearyl-CoA-desaturase activity, and RNF10, with an unknown role in weight regulating pathways, associated with adiposity and nominally increased the risk for T2D in American Indians.

A case-control analysis of common variants in GIP with type 2 diabetes and related biochemical parameters in a South Indian population.

BACKGROUND: Glucose-dependent insulinotropic polypeptide (GIP) is one of the incretins, which plays a crucial role in the secretion of insulin upon food stimulus and in the regulation of postprandial glucose level. It also exerts an effect on the synthesis and secretion of lipoprotein lipase, from adipocytes, important for lipid metabolism. The aim of our study was to do a case-control association analysis of common variants in GIP in association with type 2 diabetes and related biochemical parameters. METHOD: A total of 2000 subjects which includes 1000 (584M/416F) cases with type 2 diabetes and 1000 (470M/530F) normoglycemic control subjects belonging to Dravidian ethnicity from South India were recruited to assess the effect of single nucleotide polymorphisms (SNPs) in GIP (rs2291725, rs2291726, rs937301) on type 2 diabetes in a case-control manner. The SNPs were genotyped by using tetra primer amplification refractory mutation system-PCR (ARMS PCR). For statistical analysis, our study population was divided into sub-groups based on gender (male and female). Association analysis was carried out using chi-squared test and the comparison of biochemical parameters among the three genotypes were performed using analysis of covariance (ANCOVA). RESULT: Initial analysis revealed that, out of the total three SNPs selected for the present study, two SNPs namely rs2291726 and rs937301 were in complete linkage disequilibrium (LD) with each other. Therefore, only two SNPs, rs2291725 and rs2291726, were genotyped for the association studies. No significant difference in the allele frequency and genotype distribution of any of the SNPs in GIP were observed between cases and controls (P > 0.05). Analysis of biochemical parameters among the three genotypes showed a significant association of total cholesterol (P = 0.042) and low density lipoprotein (LDL) with the G allele of the SNP rs2291726 in GIP (P = 0.004), but this was observed only in the case of female subjects. However this association does not remain significant after correction for multiple testing by Bonferroni's inequality method. CONCLUSION: No statistically significant association was observed between any of the SNPs analysed and type 2 diabetes in our population. But the analysis of biochemical parameters indicates that the G allele in rs2291726 may be a putative risk allele for increased LDL cholesterol and further studies in other population needs to be carried out for ascertaining its role in cholesterol metabolism and subsequent cardiovascular risk.