An E115A missense variant in CERS2 is associated with increased sleeping energy expenditure and hepatic insulin resistance in American Indians.

Genetic determinants of inter-individual differences in energy expenditure (EE) are largely unknown. Sphingolipids, such as ceramides, have been implicated in the regulation of human EE via mitochondrial uncoupling. In this study, we investigated whether genetic variants within enzymes involved in sphingolipid synthesis and degradation affect EE and insulin-related traits in a cohort of American Indians informative for 24-hour EE and glucose disposal rates during a hyperinsulinemic-euglycemic clamp. Association analysis of 10,084 genetic variants within 28 genes involved in sphingolipid pathways identified a missense variant (rs267738, A>C, E115A) in exon 4 of CERS2 that was associated with higher sleeping EE (+116 kcal/day) and increased rates of endogenous glucose production during basal (+5%) and insulin-stimulated (+43%) conditions, both indicators of hepatic insulin resistance. The rs267738 variant did not affect ceramide synthesis in HepG2 cells but resulted in a 30% decrease in basal mitochondrial respiration. In conclusion, we provide evidence that the CERS2 rs267738 missense variant may influence hepatic glucose production and post-absorptive sleeping metabolic rate.

Epistasis Between HLA-DRB1*16:02:01 and SLC16A11 T-C-G-T-T Reduces Odds for Type 2 Diabetes in Southwest American Indians.

We sought to identify genetic/immunologic contributors of type 2 diabetes (T2D) in an indigenous American community by genotyping all study participants for both high-resolution HLA-DRB1 alleles and SLC16A11 to test their risk and/or protection for T2D. These genes were selected based on independent reports that HLA-DRB1*16:02:01 is protective for T2D and that SLC16A11 associates with T2D in individuals with BMI <35 kg/m2. Here, we test the interaction of the two loci with a more complete data set and perform a BMI sensitivity test. We defined the risk protection haplotype of SLC16A11, T-C-G-T-T, as allele 2 of a diallelic genetic model with three genotypes, SLC16A11*11, *12, and *22, where allele 1 is the wild type. Both earlier findings were confirmed. Together in the same logistic model with BMI ≥35 kg/m2, DRB1*16:02:01 remains protective (odds ratio [OR] 0.73), while SLC16A11 switches from risk to protection (OR 0.57 [*22] and 0.78 [*12]); an added interaction term was statistically significant (OR 0.49 [*12]). Bootstrapped (b = 10,000) statistical power of interaction, 0.4801, yielded a mean OR of 0.43. Sensitivity analysis demonstrated that the interaction is significant in the BMI range of 30-41 kg/m2. To investigate the epistasis, we used the primary function of the HLA-DRB1 molecule, peptide binding and presentation, to search the entire array of 15-mer peptides for both the wild-type and ancient human SLC16A11 molecules for a pattern of strong binding that was associated with risk and protection for T2D. Applying computer binding algorithms suggested that the core peptide at SLC16A11 D127G, FSAFASGLL, might be key for moderating risk for T2D with potential implications for type 1 diabetes.

One-hour glucose is an earlier marker of dysglycemia than two-hour glucose.

AIMS: The timing of increase in 1-hour PG and its utility as an earlier predictor of both prediabetes (PreDM) and type 2 diabetes (T2D) compared to 2-hour PG (2 h-PG) are unknown. To evaluate the timing of crossing of the 1 h-PG ≥ 155 mg/dl (8.6 mmol/L) for PreDM and 209 mg/dl (11.6 mmol/L) for T2D and respective current 2 h-PG thresholds of 140 mg/dl (7.8 mmol/L) and 200 mg/dl (11.1 mmol/L). METHODS: Secondary analysis of 201 Southwest Native Americans who were followed longitudinally for 6-10 years and had at least 3 OGTTs. RESULTS: We identified a subset of 43 individuals who first developed PreDM by both 1 h-PG and 2 h-PG criteria during the study. For most (32/43,74%), 1 h-PG ≥ 155 mg/dl was observed before 2 h-PG reached 140 mg/dl (median [IQR]: 1.7 [-0.25, 4.59] y; mean ± SEM: 5.3 ± 1.9 y). We also identified a subset of 33 individuals who first developed T2D during the study. For most (25/33, 75%), 1 h-PG reached 209 mg/dl earlier (median 1.0 [-0.56, 2.02] y; mean ± SEM: 1.6 ± 0.8 y) than 2 h-PG reached 200 mg/dl, diagnostic of T2D. CONCLUSIONS: 1 h-PG ≥ 155 mg/dl is an earlier marker of elevated risk for PreDM and T2D than 2 h-PG ≥ 140 mg/dl.

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.

Functional variants in cytochrome b5 type A (CYB5A) are enriched in Southwest American Indian individuals and associate with obesity.

OBJECTIVE: This study aimed to identify genetic variants enriched in Southwest American Indian (SWAI) individuals that associate with BMI. METHODS: Whole genome sequencing data (n = 296) were used to identify potentially functional variants that are common in SWAI individuals (minor allele frequency ≥10%) but rare in other ethnic groups (minor allele frequency < 0.1%). Enriched variants were tested for association with BMI in 5,870 SWAI individuals. One variant was studied using a luciferase reporter, and haplotypes that included this variant were analyzed for association with various measures of obesity (n = 917-5,870), 24-hour energy expenditure (24-h EE; n = 419), and skeletal muscle biopsy expression data (n = 207). RESULTS: A 5' untranslated region variant in cytochrome b5 type A (CYB5A), rs548402150, met the enrichment criteria and associated with increased BMI (ß = 2%, p = 0.004). Functionally, rs548402150 decreased luciferase expression by 30% (p = 0.003) and correlated with decreased skeletal muscle CYB5A expression (ß = -0.5 SD, p = 0.0008). Combining rs548402150 with two splicing quantitative trait loci in CYB5A identified a haplotype carried almost exclusively in SWAI individuals that associated with increased BMI (ß = 3%, p = 0.0003) and decreased CYB5A expression, whereas the most common haplotype in all ethnic groups associated with lower BMI and percentage of body fatness, increased 24-h EE, and increased CYB5A expression. CONCLUSIONS: Further studies on the effects of CYB5A on 24-h EE and BMI may provide insights into obesity-related physiology.

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.

Relationship Between Insulin Secretion and Insulin Sensitivity and its Role in Development of Type 2 Diabetes Mellitus: Beyond the Disposition Index.

We assessed whether the relationship between insulin secretion and sensitivity predicted development of type 2 diabetes in American Indians participating in a longitudinal epidemiologic study. At baseline, when all subjects did not have diabetes, 1566 participants underwent oral tests and 420 had intravenous measures of glucose regulation with estimates of insulin secretion and sensitivity. Standardized major axis regression was used to study the relationship of secretion and sensitivity. Distances away from and along the regression line estimated compensatory insulin secretion and secretory demand, respectively. This relationship differed according to glucose tolerance and BMI categories. The distance away from the line is similar to the disposition index (DI) defined as the product of estimated secretion and sensitivity, but the regression line may differ from a line with constant DI (i.e., it is not necessarily hyperbolic). Subjects with the same DI but different levels of insulin secretion and sensitivity had different incidence rates of diabetes; lower sensitivity with higher secretory demand was associated with greater diabetes risk. Insulin secretion and insulin sensitivity, analyzed together, predict diabetes better than DI alone. Physiologically, this may reflect long-term risk associated with increased allostatic load resulting from the stimulation of insulin hypersecretion by increased glycemia.

Relationship Between Insulin Secretion and Insulin Sensitivity and its Role in Development of Type 2 Diabetes Mellitus: Beyond the Disposition Index.

We assessed whether the relationship between insulin secretion and sensitivity predicted development of type 2 diabetes in American Indians participating in a longitudinal epidemiologic study. At baseline, when all subjects did not have diabetes, 1566 participants underwent oral tests and 420 had intravenous measures of glucose regulation with estimates of insulin secretion and sensitivity. Standardized major axis regression was used to study the relationship of secretion and sensitivity. Distances away from and along the regression line estimated compensatory insulin secretion and secretory demand, respectively. This relationship differed according to glucose tolerance and BMI categories. The distance away from the line is similar to the disposition index (DI) defined as the product of estimated secretion and sensitivity, but the regression line may differ from a line with constant DI (i.e., it is not necessarily hyperbolic). Subjects with the same DI but different levels of insulin secretion and sensitivity had different incidence rates of diabetes; lower sensitivity with higher secretory demand was associated with greater diabetes risk. Insulin secretion and insulin sensitivity, analyzed together, predict diabetes better than DI alone. Physiologically, this may reflect long-term risk associated with increased allostatic load resulting from the stimulation of insulin hypersecretion by increased glycemia.

Next generation sequencing for HLA loci in full heritage Pima Indians of Arizona, Part II: HLA-A, -B, and -C with selected non-classical loci at 4-field resolution from whole genome sequences.

While the samples and data from the Pima Indians of the Gila River Indian Community have been included in many international HLA workshops and conferences and have been the focus of numerous population reports and the source of novel alleles at the classical HLA loci, they have not been studied for the non-classical loci. In order to expand our HLA-disease association studies, we typed over 300 whole genome sequences from full Pima heritage members, controlled for first degree relationship, and employed recently developed computer algorithms to resolve HLA alleles. Both classical-HLA-A, -B, and -C- and non-classical- HLA-E, -F, -G, -J, -L, -W, -Y, -DPA2, -DPB2, -DMA, -DMB, -DOA, -DRB2, -DRB9, TAP1- loci were typed at the 4-field level of resolution. We present allele and selected haplotype frequencies, test the genotype distributions for population structure, discuss the issues that are created for tests of Hardy-Weinberg equilibrium over the four sample spaces of high resolution HLA typing, and address the implications for the evolution of non-classical pseudogenes that are no longer expressed in a phenotype subject to natural selection.

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.

Characterization of Exome Variants and Their Metabolic Impact in 6,716 American Indians from the Southwest US.

Applying exome sequencing to populations with unique genetic architecture has the potential to reveal novel genes and variants associated with traits and diseases. We sequenced and analyzed the exomes of 6,716 individuals from a Southwestern American Indian (SWAI) population with well-characterized metabolic traits. We found that the SWAI population has distinct allelic architecture compared to populations of European and East Asian ancestry, and there were many predicted loss-of-function (pLOF) and nonsynonymous variants that were highly enriched or private in the SWAI population. We used pLOF and nonsynonymous variants in the SWAI population to evaluate gene-burden associations of candidate genes from European genome-wide association studies (GWASs) for type 2 diabetes, body mass index, and four major plasma lipids. We found 19 significant gene-burden associations for 11 genes, providing additional evidence for prioritizing candidate effector genes of GWAS signals. Interestingly, these associations were mainly driven by pLOF and nonsynonymous variants that are unique or highly enriched in the SWAI population. Particularly, we found four pLOF or nonsynonymous variants in APOB, APOE, PCSK9, and TM6SF2 that are private or enriched in the SWAI population and associated with low-density lipoprotein (LDL) cholesterol levels. Their large estimated effects on LDL cholesterol levels suggest strong impacts on protein function and potential clinical implications of these variants in cardiovascular health. In summary, our study illustrates the utility and potential of exome sequencing in genetically unique populations, such as the SWAI population, to prioritize candidate effector genes within GWAS loci and to find additional variants in known disease genes with potential clinical impact.

Assessing established BMI variants for a role in nighttime eating behavior in robustly phenotyped Southwestern American Indians.

BACKGROUND/OBJECTIVES: Nighttime eating (NE) behavior has a genetic component and predicts weight gain. We hypothesized that some genetic variants, which affect NE would also show an effect on body mass index (BMI). We aimed to determine which known BMI variants associate with NE in Southwestern American Indians (SWAIs), who are at elevated risk for obesity. METHODS: Known BMI variants from the GIANT-UK Biobank meta-analysis (N = 700,000) were analysed in SWAIs characterized for NE during an inpatient 3-day protocol. Variants were analysed for association with NE using whole-genome sequence data from 50 SWAIs (23 cases and 27 controls) and selected variants were genotyped in an additional 32 SWAIs (13 NE cases and 19 controls). Variants associated with NE in a meta-analysis of the two SWAI samples were further analysed for association with nightly caloric intake and functionality in hypothalamus, pituitary, and adrenal tissues. RESULTS: Variants were identified where the allele that associated with increased BMI in the GIANT-UK Biobank meta-analysis (P ≤ 1 × 10-8) also had a P < 0.01 for increased NE in the SWAI meta-analysis. These variants were captured by six tagSNPs. Comparison of the nightly calorie intake by genotype and eQTL data from relevant tissues highlighted rs3753612 upstream of HCRTR1. CONCLUSIONS: Our strategy led to the HCRTR1 locus, which has previously been linked to sleep regulation and feeding. Although this is an intriguing candidate gene for NE, further studies in larger samples and different populations are required to validate the role of HCRTR1 in NE.

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.

Impaired Metabolic Flexibility to High-Fat Overfeeding Predicts Future Weight Gain in Healthy Adults.

The ability to switch fuels for oxidation in response to changes in macronutrient composition of diet (metabolic flexibility) may be informative of individuals' susceptibility to weight gain. Seventy-nine healthy, weight-stable participants underwent 24-h assessments of energy expenditure and respiratory quotient (RQ) in a whole-room calorimeter during energy balance (EBL) (50% carbohydrate, 30% fat) and then during 24-h fasting and three 200% overfeeding diets in a crossover design. Metabolic flexibility was defined as the change in 24-h RQ from EBL during fasting and standard overfeeding (STOF) (50% carbohydrate, 30% fat), high-fat overfeeding (HFOF) (60% fat, 20% carbohydrate), and high-carbohydrate overfeeding (HCOF) (75% carbohydrate, 5% fat) diets. Free-living weight change was assessed after 6 and 12 months. Compared with EBL, RQ decreased on average by 9% during fasting and by 4% during HFOF but increased by 4% during STOF and by 8% during HCOF. A smaller decrease in RQ, reflecting a smaller increase in lipid oxidation rate, during HFOF but not during the other diets predicted greater weight gain at both 6 and 12 months. An impaired metabolic flexibility to acute HFOF can identify individuals prone to weight gain, indicating that an individual's capacity to oxidize dietary fat is a metabolic determinant of weight change.

Next generation sequencing and the classical HLA loci in full heritage Pima Indians of Arizona: Defining the core HLA variation for North American Paleo-Indians.

The Pima Indians of the Gila River Indian Community in Arizona have participated in a long-range study of type 2 diabetes mellitus since 1965 and have been the subject of HLA typing and population studies since the early days of serological assays. These data have been in numerous HLA workshops and conferences and have been the source of at least five novel alleles at the classical HLA loci. In recent time nearly the entire study group was subject to next generation sequencing by whole genome or exome technologies, which has allowed us to HLA type over 3000 full heritage persons with recently developed computer algorithms. We present here the results for the classical HLA Loci: HLA-A, B, C, DRA, DRB1, DRB3, DRB4, DRB5, DPA1, DPB1, DQA1, and DQB1 to the third field of resolution for synonymous alleles and type the likely four field resolution alleles from the subset of whole genome sequences. Allele frequencies, and haplotype frequencies at up to five loci, are presented as well as measures of population structure and heterozygosity. We define a core set of HLA variation that approximates the distribution for the Paleo-Indians and impute nine-locus, 4-field haplotypes that are expected to be common in full heritage peoples.

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.

FGF21 Is a Hormonal Mediator of the Human "Thrifty" Metabolic Phenotype.

Fibroblast growth factor 21 (FGF21) regulates energy expenditure (EE) and influences weight change during low-protein overfeeding in rodent models. The change in EE after a low-protein overfeeding diet is a predictor of weight change in humans and a feature of the "thrifty" metabolic phenotype. However, there are no studies showing an association between circulating FGF21 and diet-related EE in humans. We assessed the changes in plasma FGF21 concentrations after 24 h of seven dietary interventions with different macronutrient content while in a whole-room indirect calorimeter in 64 healthy subjects with normal glucose regulation. Plasma FGF21 concentration consistently increased by threefold only after the two low-protein (3%) overfeeding diets, one high in carbohydrate (75%) and the other high in fat (46%), with larger increases in FGF21 being associated with greater increases in 24-h EE. Subjects with smaller increases in FGF21 after the low-protein high-fat diet gained more weight after 6 months in free-living conditions. Therefore, the individual predisposition to weight gain over time can be assessed by 24-h overfeeding a low-protein diet and measurements of plasma FGF21 concentrations. Individuals with a blunted FGF21 response to a low-protein diet have a thrifty metabolism and are at risk for future weight gain.