Exome Sequencing of 21 Bardet-Biedl Syndrome (BBS) Genes to Identify Obesity Variants in 6,851 American Indians.

OBJECTIVE: In an ongoing effort to identify the genetic variation that contributes to obesity in American Indians, known Bardet-Biedl syndrome (BBS) genes were analyzed for an effect on BMI and leptin signaling. METHODS: Potentially deleterious variants (Combined Annotation Dependent Depletion score > 20) in BBS genes were identified in whole-exome sequence data from 6,851 American Indians informative for BMI. Common variants (detected in ≥ 10 individuals) were analyzed for association with BMI; rare variants (detected in < 10 individuals) were analyzed for mean BMI of carriers. Functional assessment of variants' effect on signal transducer and activator of transcription 3 (STAT3) activity was performed in vitro. RESULTS: One common variant, rs59252892 (Thr549Ile) in BBS9, was associated with BMI (P = 0.0008, ß = 25% increase per risk allele). Among rare variants for which carriers had severe obesity (mean BMI > 40 kg/m2 ), four were in BBS9. In vitro analysis of BBS9 found the Ile allele at Thr549Ile had a 20% increase in STAT3 activity compared with the Thr allele (P = 0.01). Western blot analysis showed the Ile allele had a 15% increase in STAT3 phosphorylation (P = 0.006). Comparable functional results were observed with Ser545Gly and Val209Leu but not Leu665Phe and Lys810Glu. CONCLUSIONS: Potentially functional variants in BBS genes in American Indians are reported. However, functional evidence supporting a causal role for BBS9 in obesity is inconclusive.

A genome-wide association study in American Indians implicates DNER as a susceptibility locus for type 2 diabetes.

Most genetic variants associated with type 2 diabetes mellitus (T2DM) have been identified through genome-wide association studies (GWASs) in Europeans. The current study reports a GWAS for young-onset T2DM in American Indians. Participants were selected from a longitudinal study conducted in Pima Indians and included 278 cases with diabetes with onset before 25 years of age, 295 nondiabetic controls ≥45 years of age, and 267 siblings of cases or controls. Individuals were genotyped on a ∼1M single nucleotide polymorphism (SNP) array, resulting in 453,654 SNPs with minor allele frequency >0.05. SNPs were analyzed for association in cases and controls, and a family-based association test was conducted. Tag SNPs (n = 311) were selected for 499 SNPs associated with diabetes (P < 0.0005 in case-control analyses or P < 0.0003 in family-based analyses), and these SNPs were genotyped in up to 6,834 additional Pima Indians to assess replication. Rs1861612 in DNER was associated with T2DM (odds ratio = 1.29 per copy of the T allele; P = 6.6 × 10(-8), which represents genome-wide significance accounting for the number of effectively independent SNPs analyzed). Transfection studies in murine pancreatic ß-cells suggested that DNER regulates expression of notch signaling pathway genes. These studies implicate DNER as a susceptibility gene for T2DM in American Indians.

Identification of genetic variation that determines human trehalase activity and its association with type 2 diabetes.

A prior linkage scan in Pima Indians identified a putative locus for type two diabetes (T2D) and body mass index (BMI) on chromosome 11q23-25. Association mapping across this region identified single nucleotide polymorphisms (SNPs) in the trehalase gene (TREH) that were associated with T2D. To assess the putative connection between trehalase activity and T2D, we performed a linkage study for trehalase activity in 570 Pima Indians who had measures of trehalase activity. Strong evidence of linkage of plasma trehalase activity (LOD = 7.0) was observed in the TREH locus. Four tag SNPs in TREH were genotyped in these subjects and plasma trehalase activity was highly associated with three SNPs: rs2276064, rs117619140 and rs558907 (p = 2.2 × 10(-11)-1.4 × 10(-23)), and the fourth SNP, rs10790256, was associated conditionally on these three (p = 2.9 × 10(-7)). Together, the four tag SNPs explained 51 % of the variance in plasma trehalase activity and 79 % of the variance attributed to the linked locus. These four tag SNPs were further genotyped in 828 subjects used for association mapping of T2D, and rs558907 was associated with T2D (odds ratio (OR) 1.94, p = 0.002). To assess replication of the T2D association, all four tag SNPs were additionally genotyped in two non-overlapping samples of Native Americans. Rs558907 was reproducibly associated with T2D in 2,942 full-heritage Pima Indians (OR 1.27 p = 0.03) and 3,897 "mixed" heritage Native Americans (OR 1.21, p = 0.03), and the strongest evidence for association came from combining all samples (OR 1.27 p = 1.6 × 10(-4), n = 7,667). However, among 320 longitudinally studied subjects, measures of trehalase activity from a non-diabetic exam did not predict those who would eventually develop diabetes versus those who would remain non-diabetic (hazard ratio 0.94 per SD of trehalase activity, p = 0.29). We conclude that variants in TREH control trehalase activity, and although one of these variants is also reproducibly associated with T2D, it is likely that the effect of the SNP on risk of T2D occurs by a mechanism different than affecting trehalase activity. Alternatively, TREH variants may be tagging a nearby T2D locus.

Variants in ASK1 are associated with skeletal muscle ASK1 expression, in vivo insulin resistance, and type 2 diabetes in Pima Indians.

OBJECTIVE: Prior genome-wide association and exon array expression studies both provided suggestive evidence that apoptosis signal regulating kinase 1 (ASK1) may influence in vivo insulin action in Pima Indians. Genetic variants in or near ASK1 were analyzed to assess the role of this gene in insulin action and type 2 diabetes. RESEARCH DESIGN AND METHODS: Genotypic data from 31 variants were used to determine the linkage disequilibrium pattern across ASK1 in Pima Indians. Eight tag SNPs were initially genotyped in 3,501 full-heritage Pima Indians. Replication for association with diabetes was assessed in a second population-based sample of 3,723 Native Americans and the published DIAGRAM study. Quantitative traits were analyzed in 536 nondiabetic Native Americans, and ASK1 expression was examined in skeletal muscle of 153 nondiabetic Native Americans. RESULTS: Three tag SNPs were associated with type 2 diabetes (rs35898099, P = 0.003, odds ratio [95% CI] 1.27 [1.08-1.47]; rs1570056, P = 0.007, 1.19 [1.05-1.36]; rs7775356, P = 0.04, 1.14 [1.01-1.28]) in the full-heritage Pima Indians. The association with rs35898099 was replicated in a second sample of Native Americans (P = 0.04, 1.22 [1.01-1.47]), while that for rs1570056 was replicated in the DIAGRAM study of Caucasians (Z statistic based P = 0.026; fixed-effect model, 1.06 [1.00-1.12]). The diabetes risk allele for rs1570056 was associated with reduced insulin action as assessed by either HOMA-IR in 2,549 nondiabetic full-heritage Pima Indians (P = 0.027) or a hyperinsulinemic-euglycemic clamp among 536 nondiabetic Native Americans (P = 0.02). Real-time PCR identified a positive correlation between ASK1 expression in skeletal muscle biopsies and in vivo insulin action (P = 0.02, r = 0.23), and the risk allele for rs1570056 was associated with lower ASK1 expression (P = 0.003, r = -0.22). CONCLUSIONS: ASK1 variants may increase susceptibility to type 2 diabetes by decreasing insulin sensitivity via reduced ASK1 expression.