Functional characterization of a novel p.Ser76Thr variant in IGFBP4 that associates with body mass index in American Indians.

Insulin-like growth factor binding protein 4 (IGFBP4) is involved in adipogenesis, and IGFBP4 null mice have decreased body fat through decreased PPAR-γ expression. In the current study, we assessed whether variation in the IGFBP4 coding region influences body mass index (BMI) in American Indians who are disproportionately affected by obesity. Whole exome sequence data from a population-based sample of 6779 American Indians with longitudinal measures of BMI were used to identify variation in IGFBP4 that associated with BMI. A novel variant that predicts a p.Ser76Thr in IGFBP4 (Thr-allele frequency = 0.02) was identified which associated with the maximum BMI measured during adulthood (BMI 39.8 kg/m2 for Thr-allele homozygotes combined with heterozygotes vs. 36.2 kg/m2 for Ser-allele homozygotes, ß = 6.7% per Thr-allele, p = 8.0 × 10-5, adjusted for age, sex, birth-year and the first five genetic principal components) and the maximum age- and sex-adjusted BMI z-score measured during childhood/adolescence (z-score 0.70 SD for Thr-allele heterozygotes vs. 0.32 SD for Ser-allele homozygotes, ß = 0.37 SD per Thr-allele, p = 8.8 × 10-6). In vitro functional studies showed that IGFBP4 with the Thr-allele (BMI-increasing) had a 55% decrease (p = 0.0007) in FOXO-induced transcriptional activity, reflecting increased activation of the PI3K/AKT pathway mediated through increased IGF signaling. Over-expression and knock-down of IGFBP4 in OP9 cells during differentiation showed that IGFBP4 upregulates adipogenesis through PPARγ, CEBPα, AGPAT2 and SREBP1 expression. We propose that this American Indian specific variant in IGFBP4 affects obesity via an increase of IGF signaling.

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 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.

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.

Weight tracking in childhood and adolescence and type 2 diabetes risk.

AIMS/HYPOTHESIS: The aim of this work was to examine the associations of average weight and weight velocity in three growth periods from birth through adolescence with type 2 diabetes incidence. METHODS: Child participants were selected from a 43 year longitudinal study of American Indians to represent three growth periods: pre-adolescence (birth to ~8 years); early adolescence (~8 to ~13 years); and late adolescence (~13 to ~18 years). Age-, sex- and height-standardised weight z score mean and weight z score velocity (change/year) were computed for each period. Participants were followed for up to 25 years from the end of each growth period until they developed diabetes. Associations of weight z score mean or weight z score velocity with diabetes incidence were determined with sex-, birth date- and maternal diabetes-adjusted Poisson regression models. RESULTS: Among 2100 participants representing the pre-adolescence growth period, 1558 representing the early adolescence period and 1418 representing the late adolescence period, there were 290, 315 and 380 incident diabetes cases, respectively. During the first 10 years of follow-up, the diabetes incidence rate ratio (95% CI) was 1.72 (1.40, 2.11)/SD of log10 weight z score mean in pre-adolescence, 2.09 (1.68, 2.60)/SD of log10 weight z score mean in early adolescence and 1.85 (1.58, 2.17)/SD of log10 weight z score mean in late adolescence. The diabetes incidence rate ratio (95% CI) was 1.79 (1.49, 2.17)/SD of log10 weight z score velocity in pre-adolescence, 1.13 (0.91, 1.41)/SD of log10 weight z score velocity in early adolescence and 1.29 (1.09, 1.51)/SD of log10 weight z score velocity in late adolescence. There were strong correlations in the weight z score means and weak correlations in the weight z score velocities between successive periods. CONCLUSIONS/INTERPRETATION: Higher weight and accelerated weight gain in all growth periods associate with increased type 2 diabetes risk. Importantly, higher weight and greater weight velocity during pre-adolescence jointly associate with the highest type 2 diabetes risk. Graphical abstract.

Low Serum Insulinlike Growth Factor II Levels Correlate with High BMI in American Indian Adults.

OBJECTIVE: Insulinlike growth factor II (IGF-II) regulates metabolism and growth. In humans, both positive and negative relationships have been reported between serum IGF-II levels and obesity. This study assessed the relationship between serum IGF-II levels and BMI and determined whether IGF-II levels predict weight gain. METHODS: Serum samples were available from 911 American Indians with a recorded BMI. IGF-II was measured using enzyme-linked immunosorbent assay. RESULTS: Serum IGF-II levels were negatively correlated with BMI (r = -0.17, P = 4.4 × 10-7 , adjusted for age, sex, and storage time). The strongest correlation was in participants aged ≥ 30 years (r = -0.28, P = 3.4 × 10-8 , N = 349), a modest correlation was in participants aged 20 to 29 years (r = -0.15, P = 7.6 × 10-3 , N = 322), and participants aged 15 to 19 years had no correlation (r = 0.05, P = 0.48, N = 240). IGF-II levels did not predict weight gain. However, among individuals who had genotypes for 64 established obesity variants (age ≥ 20 years, N = 671), a genetic risk score for high BMI was associated with lower IGF-II (ß = -0.08 SD of IGF-II per SD of the genetic risk score, P = 0.025). CONCLUSIONS: There is a negative relationship between IGF-II levels and BMI, in which the correlation is stronger at older ages. The association between genetic risk for BMI and IGF-II levels suggests that this correlation may be due to an effect of obesity on IGF-II.

Assessing the Role of 98 Established Loci for BMI in American Indians.

OBJECTIVE: Meta-analyses of genome-wide association studies in Europeans have identified > 98 loci for BMI. Transferability of these established associations in Pima Indians was analyzed. METHODS: Among 98 lead single nucleotide polymorphisms (SNPs), 82 had minor allele frequency ≥ 0.01 in Pima Indians and were analyzed for association with the maximum BMI in adulthood (n = 3,491) and BMI z score in childhood (n = 1,958). Common tag SNPs across 98 loci were also analyzed for additional signals. RESULTS: Among the lead SNPs, 13 (TMEM18, TCF7L2, MRPS33P4, PRKD1, ZFP64, FTO, TAL1, CALCR, GNPDA2, CREB1, LMX1B, ADCY9, NLRC3) were associated with BMI (P  ≤ 0.05) in Pima adults. A multi-allelic genetic risk score (GRS), which summed the risk alleles for 82 lead SNPs, showed a significant trend for a positive relationship between GRS and BMI in adulthood (beta = 0.48% per risk allele; P = 1.6 × 10-9 ) and BMI z score in childhood (beta = 0.024 SD; P = 1.7 × 10-7 ). GRS was significantly associated with BMI across all age groups ≥ 5 years, except for those ≥ 50 years. The strongest association was seen in adolescence (age 14-16 years; P = 1.84 × 10-9 ). CONCLUSIONS: In aggregate, European-derived lead SNPs had a notable effect on BMI in Pima Indians. Polygenic obesity in this population manifests strongly in childhood and adolescence and persists throughout much of adult life.

Identification and functional analysis of a novel G310D variant in the insulin-like growth factor 1 receptor (IGF1R) gene associated with type 2 diabetes in American Indians.

AIMS: Insulin-like growth factor 1 receptor (IGF1R) is involved in cell growth and glucose homeostasis. In the current study, the IGF1R locus was analysed as a candidate gene for type 2 diabetes (T2D) in American Indians. MATERIALS AND METHODS: Whole genome sequence data from 335 American Indians identified 3 novel missense variants in IGF1R. The associations of IGF1R variants with T2D, age of T2D onset and birth weight were analysed in a population-based sample of 7701 American Indians. RESULTS: A novel glycine-to-aspartic acid substitution (G310D) in IGF1R was identified, which associated with T2D in a sex-specific manner (Psex interaction = 0.02). In women, the aspartic acid (D) allele (frequency = 0.034) was associated with increased risk for T2D (n = 4292, P = 2.0 × 10-5 adjusted for age, birth year, and the first 5 genetic principal components; odds ratio [OR] = 2.23 [1.54-3.23] per risk allele) and an earlier age of T2D onset (n = 4292, P = 2 × 10-4 , hazard rate ratio = 1.45 [1.20-1.75], Psex interaction = 0.05). Female carriers of the D-allele also had lower birth weight (n = 1313, ß = -163 g, P = .006, Psex interaction = 0.008). Among 85 siblings discordant for G310D, carriers of the D-allele had shorter stature as compared with carriers of the G-allele (ß = -1.6 cm, P = .001, within family model). The G310D variant was functionally studied in vitro, where the D-allele had a 22% increase (P = .0005) in FOXO1-induced transcriptional activity, due to decreased activation of the PI3K/AKT pathway mediated through reduced IGF1R activity. CONCLUSION: A unique G310D variant in IGF1R, which occurs in 6% American Indians, may impair IGF1R signalling pathways, thereby increasing the risk of T2D.

A Genome-Wide Association Study Using a Custom Genotyping Array Identifies Variants in GPR158 Associated With Reduced Energy Expenditure in American Indians.

Pima Indians living in Arizona have a high prevalence of obesity, and we have previously shown that a relatively lower energy expenditure (EE) predicts weight and fat mass gain in this population. EE is a familial trait (heritability = 0.52); therefore, in the current study, we aimed to identify genetic variants that affect EE and thereby influence BMI and body fatness in Pima Indians. Genotypic data from 491,265 variants were analyzed for association with resting metabolic rate (RMR) and 24-h EE assessed in a whole-room calorimeter in 507 and 419 Pima Indians, respectively. Variants associated with both measures of EE were analyzed for association with maximum BMI and percent body fat (PFAT) in 5,870 and 912 Pima Indians, respectively. rs11014566 nominally associated with both measures of EE and both measures of adiposity in Pima Indians, where the G allele (frequency: Pima Indians = 0.60, Europeans <0.01) associated with lower 24-h EE (ß = -33 kcal/day per copy), lower RMR (ß = -31 kcal/day), higher BMI (ß = +0.6 kg/m2), and higher PFAT (ß = +0.9%). However, the association of rs11014566 with BMI did not directionally replicate when assessed in other ethnic groups. rs11014566 tags rs144895904, which affected promoter function in an in vitro luciferase assay. These variants map to GPR158, which is highly expressed in the brain and interacts with two other genes (RGS7 and CACNA1B) known to affect obesity in knockout mice. Our results suggest that common ethnic-specific variation in GPR158 may influence EE; however, its role in weight gain remains controversial, as it either had no association with BMI or associated with BMI but in the opposite direction in other ethnic groups.

Assessing variation across 8 established East Asian loci for type 2 diabetes mellitus in American Indians: Suggestive evidence for new sex-specific diabetes signals in GLIS3 and ZFAND3.

BACKGROUND: Eight new loci for type 2 diabetes mellitus (T2DM) were identified in an East Asian genome-wide association study meta-analysis. We assess tag SNPs across these loci for associations with T2DM in American Indians. METHODS: A total of 435 SNPs that tag (R2  ≥ .85) common variation across the 8 loci were analyzed for association with T2DM (n = 7710), early onset T2DM (n = 1060), body mass index (n = 6839), insulin sensitivity (n = 555), and insulin secretion (n = 298). RESULTS: Tag SNPs within FITM2-R3HDML-HNF4A, GLIS3, KCNK16, and ZFAND3 associated with T2DM after accounting for locus-wide multiple testing. The T2DM association in FITM2-R3HDML-HNF4A (rs3212183; P = .0002; OR = 1.19 [1.09-1.30]) was independent from the East Asian lead SNP (rs6017317), which did not associate with T2DM in American Indians. The top signals in GLIS3 (rs7875253; P = .0004; OR = 1.23 [1.10-1.38]) and KCNK16 (rs1544050; P = .002; OR = 1.16 [1.06-1.27]) were attenuated after adjustment for the East Asian lead SNPs (rs7041847 in GLIS3; rs1535500 in KCNK16), both of which also associated with T2DM in American Indians (P = .02; OR = 1.11 [1.01-1.21]; P = .007; OR = 1.19 [1.05-1.36] respectively). The top SNP in ZFAND3 (rs9470794; P = .002; OR = 1.43 [1.14-1.80]) was the identical East Asian lead SNP. Additional SNPs in GLIS3 (rs180867004) and ZFAND3 (rs4714120 and rs9470701) had significant genotype × sex interactions (P ≤ .008). The GLIS3 SNP (rs180867004) associated with T2DM only in men (P = .00006, OR = 1.94 [1.40-2.68]). The ZFAND3 SNPs (rs4714120 and rs9470701) associated with T2DM only in women (P = .0002, OR = 1.35 [1.16-1.59]; P = .0003, OR = 1.37 [1.16-1.63] respectively). CONCLUSIONS: Replication of lead T2DM SNPs in GLIS3, KCNK16, and ZFAND3 was observed in American Indians. Sex-specific T2DM signals in GLIS3 and ZFAND3, which are distinct from the East Asian GWAS signals, were also identified.

Assessing FOXO1A as a potential susceptibility locus for type 2 diabetes and obesity in American Indians.

OBJECTIVE: A prior genome-wide association study (GWAS) in Pima Indians identified variation within FOXO1A that modestly associated with early-onset (onset age < 25 years) type 2 diabetes (T2D). FOXO1A encodes the forkhead transcription factor involved in pancreatic ß-cell growth and hypothalamic energy balance; therefore, FOXO1A was analyzed as a candidate gene for T2D and obesity in a population-based sample of 7,710 American Indians. METHODS: Tag SNPs in/near FOXO1A (minor allele frequency ≥ 0.05) were analyzed for association with T2D at early onset (n = 1,060) and all ages (n = 7,710) and with insulin secretion (n = 298). SNPs were also analyzed for association with maximum body mass index (BMI) in adulthood (n = 5,918), maximum BMI z-score in childhood (n = 5,350), and % body fat (n = 555). RESULTS: An intronic SNP rs2297627 associated with early-onset T2D [OR = 1.34 (1.13-1.58), P = 8.7 × 10(-4)] and T2D onset at any age [OR = 1.19 (1.09-1.30), P = 1 × 10(-4) ]. The T2D risk allele also associated with lower acute insulin secretion (ß = 0.88, as a multiplier, P = 0.02). Another intronic SNP (rs1334241, D' = 0.99, r(2) = 0.49 with rs2297627) associated with maximum adulthood BMI (ß = 1.02, as a multiplier, P = 3 × 10(-5)), maximum childhood BMI z-score (ß = 0.08, P = 3 × 10(-4)), and % body fat (ß = 0.83%, P = 0.04). CONCLUSIONS: Common variation in FOXO1A may modestly affect risk for T2D and obesity in American Indians.

The effect of differing patterns of childhood body mass index gain on adult physiology in American Indians.

OBJECTIVE: Identifying groups of individuals with similar patterns of body mass index (BMI) change during childhood may increase understanding of the relationship between childhood BMI and adult health. METHODS: Discrete classes of BMI z-score change were determined in 1,920 American Indian children with at least four non diabetic health examinations between the ages of 2 and 18 years using latent class trajectory analysis. In subsets of subjects, data were available for melanocortin-4 receptor (MC4R) sequencing; in utero exposure to type 2 diabetes (T2D); or, as adults, oral glucose tolerance tests, onset of T2D, or body composition. RESULTS: Six separate groups were identified. Individuals with a more modern birth year, an MC4R mutation, or in utero exposure to T2D were clustered in the two groups with high increasing and chronic overweight z-scores (P < 0.0001). The z-score classes predicted adult percent fat (P < 0.0001, partial r(2) = 0.18 adjusted for covariates). There was a greater risk for T2D, independent from adult BMI, in three classes (lean increasing to overweight, high increasing, and chronic overweight z-scores) compared to the two leanest groups (respectively: HRR= 3.2, P = 0.01; 6.0, P = 0.0003; 11.6, P < 0.0001). CONCLUSIONS: Distinct patterns of childhood BMI z-score change associate with adult adiposity and may impact risk of T2D.

A cis-eQTL in PFKFB2 is associated with diabetic nephropathy, adiposity and insulin secretion in American Indians.

A prior genome-wide association study (GWAS) in Pima Indians identified a variant within PFKFB2 (rs17258746) associated with body mass index (BMI). PFKFB2 encodes 6-phosphofructo-2-kinase/fructose 2,6-bisphosphatase isoform 2, which plays a role in glucose metabolism. To follow-up on the GWAS, tag SNPs across PFKFB2 were genotyped in American Indians (AI) who had longitudinal data on BMI (n = 6839), type 2 diabetes (T2D; n = 7710), diabetic nephropathy (DN; n = 2452), % body fat (n = 555) and insulin secretion (n = 298). Two SNPs were further genotyped in urban AI to assess replication for DN (n = 864). PFKFB2 expression was measured in 201 adipose biopsies using real-time RT-PCR and 61 kidney biopsies using the Affymetrix U133 array. Two SNPs (rs17258746 and rs11120137), which capture the same signal, were associated with maximum BMI in adulthood (ß = 1.02 per risk allele, P = 7.3 × 10(-4)), maximum BMI z-score in childhood (ß = 0.079, P = 0.03) and % body fat in adulthood (ß = 3.4%, P = 3 × 10(-7)). The adiposity-increasing allele correlated with lower PFKFB2 adipose expression (ß = 0.81, P = 9.4 × 10(-4)). Lower expression of PFKFB2 further correlated with higher % body fat (r = -0.16, P = 0.02) and BMI (r = -0.17, P = 0.02). This allele was also associated with increased risk for DN in both cohorts of AI [odds ratio = 1.64 (1.32-2.02), P = 5.8 × 10(-6)], and similarly correlated with lower PFKFB2 expression in kidney glomeruli (ß = 0.87, P = 0.03). The same allele was also associated with lower insulin secretion assessed by acute insulin response (ß = 0.78, P = 0.03) and 30-min plasma insulin concentrations (ß = 0.78, P = 1.1 × 10(-4)). Variation in PFKFB2 appears to reduce PFKFB2 expression in adipose and kidney tissues, and thereby increase risk for adiposity and DN.

Common genetic variation in and near the melanocortin 4 receptor gene (MC4R) is associated with body mass index in American Indian adults and children.

Six rare functional coding mutations were previously identified in melanocortin 4 receptor (MC4R) in 6,760 American Indians. Individuals heterozygous for one of these mutations become obese while young. We now investigate whether common non-coding variation near MC4R also contributes to obesity. Fifty-six tag single-nucleotide polymorphisms (SNPs) were genotyped in 3,229 full-heritage Pima Indians, and nine of these SNPs which showed evidence for association were genotyped in additional 3,852 mixed-heritage American Indians. Associations of SNPs with maximum body mass index (BMI) in adulthood (n = 5,918), BMI z score in childhood (n = 5,350), percent body fat (n = 864), energy expenditure (n = 358) and ad libitum food intake (n = 178) were assessed. Conditional analyses demonstrated that SNPs, rs74861148 and rs483125, were independently associated with BMI in adulthood (ß = 0.68 kg/m(2) per risk allele, p = 5 × 10(-5); ß = 0.58 kg/m(2), p = 0.002, respectively) and BMI z score in childhood (ß = 0.05, p = 0.02; ß = 0.07, p = 0.01, respectively). One haplotype (frequency = 0.35) of the G allele at rs74861148 and the A allele at rs483125 provided the strongest evidence for association with adult BMI (ß = 0.89 kg/m(2), p = 5.5 × 10(-7)), and was also associated with childhood BMI z score (ß = 0.08, p = 0.001). In addition, a promoter SNP rs11872992 was nominally associated with adult BMI (ß = 0.61 kg/m(2), p = 0.05) and childhood BMI z score (ß = 0.11, p = 0.01), where the risk allele also modestly decreased transcription in vitro by 12 % (p = 0.005). This risk allele was further associated with increased percent body fat (ß = 2.2 %, p = 0.002), increased food intake (ß = 676 kcal/day, p = 0.007) and decreased energy expenditure (ß = -53.4 kcal/day, p = 0.054). Common and rare variation in MC4R contributes to obesity in American Indians.

Common genetic variation in the glucokinase gene (GCK) is associated with type 2 diabetes and rates of carbohydrate oxidation and energy expenditure.

AIMS/HYPOTHESIS: Glucokinase (GCK) plays a role in glucose metabolism and glucose-stimulated insulin secretion. Rare mutations in GCK cause MODY. We investigated whether common variation (minor allele frequency ≥0.01) in GCK is associated with metabolic traits and type 2 diabetes. METHODS: Four exonic single-nucleotide polymorphisms (SNPs) and three SNPs predicted to cause loss of promoter function were identified in whole-genome sequence data from 234 Pima Indians. These seven tag SNPs and rs4607517, a type 2 diabetes variant established in other studies, were analysed in 415 full-heritage non-diabetic Pima Indians characterised for metabolic traits, and 7,667 American Indians who had data on type 2 diabetes and BMI. RESULTS: A novel 3' untranslated region (3'UTR) SNP, chr7:44184184-G/A, was associated with the rate of carbohydrate oxidation post-absorptively (ß = 0.22 mg [kg estimated metabolic body size (EMBS)](-1) min(-1), p = 0.005) and during a hyperinsulinaemic-euglycaemic clamp (ß = 0.24 mg [kg EMBS](-1) min(-1), p = 0.0002), the rate of carbohydrate oxidation in a respiratory chamber (ß = 311 kJ/day, p = 0.03) and 24 h energy expenditure, which was attributable to the thermic effect of food (ß = 520 kJ/day, p = 3.39 × 10(-6)). This 3'UTR SNP was also associated with diabetes (OR 1.36, 95% CI 1.11, 1.65, p = 0.002), where the A allele (allele frequency 0.05) was associated with a lower rate of carbohydrate oxidation, lower 24 h energy expenditure and higher risk for diabetes. In a Cox proportional hazards model, a rate of insulin-stimulated carbohydrate oxidation lower than the mean rate at baseline predicted a higher risk for developing diabetes than for those above the mean (hazard rate ratio 2.2, 95% CI 1.3, 3.6, p = 0.002). CONCLUSIONS/INTERPRETATION: Common variation in GCK influences the rate of carbohydrate oxidation, 24 h energy expenditure and diabetes risk in Pima Indians.

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.

MAP2K3 is associated with body mass index in American Indians and Caucasians and may mediate hypothalamic inflammation.

To identify genes that affect body mass index (BMI) in American Indians who are predominately of Pima Indian heritage, we previously completed a genome-wide association study in 1120 American Indians. That study also included follow-up genotyping for 9 SNPs in 2133 additional subjects. A comprehensive follow-up study has subsequently been completed where 292 SNPs were genotyped in 3562 subjects, of which 128 SNPs were assessed for replication in 3238 additional subjects. In the combined subjects (n = 6800), BMI associations for two SNPs, rs12882548 and rs11652094, approached genome-wide significance (P = 6.7 × 10(-7) and 8.1 × 10(-7), respectively). Rs12882548 is located in a gene desert on chromosome 14 and rs11652094 maps near MAP2K3. Several SNPs in the MAP2K3 region including rs11652094 were also associated with BMI in Caucasians from the GIANT consortium (P = 10(-2)-10(-5)), and the combined P-values across both American Indians and Caucasian were P = 10(-4)-10(-9). Follow-up sequencing across MAP2K3 identified several paralogous sequence variants indicating that the region may have been duplicated. MAP2K3 expression levels in adipose tissue biopsies were positively correlated with BMI, although it is unclear if this correlation is a cause or effect. In vitro studies with cloned MAP2K3 promoters suggest that MAP2K3 expression may be up-regulated during adipogenesis. Microarray analyses of mouse hypothalamus cells expressing constitutively active MAP2K3 identified several up-regulated genes involved in immune/inflammatory pathways and a gene, Hap1, thought to play a role in appetite regulation. We conclude that MAP2K3 is a reproducible obesity locus that may affect body weight via complex mechanisms involving appetite regulation and hypothalamic inflammation.

Strong parent-of-origin effects in the association of KCNQ1 variants with type 2 diabetes in American Indians.

Parent-of-origin effects were observed in an Icelandic population for several genetic variants associated with type 2 diabetes, including those in KLF14 (rs4731702), MOB2 (rs2334499), and KCNQ1 (rs2237892, rs231362). We analyzed parent-of-origin effects for these variants, along with two others in KCNQ1 identified in previous genome-wide association studies (rs2237895, rs2299620), in 7,351 Pima Indians from 4,549 nuclear families; 34% of participants had diabetes. In a subset of 287 normoglycemic individuals, acute insulin secretion was measured by an intravenous glucose tolerance test. Statistically significant (P < 0.05) parent-of-origin effects were seen for association with type 2 diabetes for all variants. The strongest effect was seen at rs2299620 in KCNQ1; the C allele was associated with increased diabetes when maternally derived (odds ratio [OR], 1.92; P = 4.1 × 10(-12)), but not when paternally derived (OR, 0.93; P = 0.47; P = 9.9 × 10(-6) for difference in maternal and paternal effects). A maternally derived C allele also was associated with a 28% decrease in insulin secretion (P = 0.002). This study confirms parent-of-origin effects in the association with type 2 diabetes for variants in KLF14, MOB2, and KCNQ1. In Pima Indians, the effect of maternally derived KCNQ1 variants appears to be mediated through decreased insulin secretion and is particularly strong, accounting for 4% of the variance in liability to diabetes.

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.