Novel signals and polygenic score for height are associated with pubertal growth traits in Southwestern American Indians.

Most genetic variants associated with adult height have been identified through large genome-wide association studies (GWASs) in European-ancestry cohorts. However, it is unclear how these variants influence linear growth during adolescence. This study uses anthropometric and genotypic data from a longitudinal study conducted in an American Indian community in Arizona between 1965-2007. Growth parameters (i.e. height, velocity, and timing of growth spurt) were derived from the Preece-Baines growth model, a parametric growth curve fitted to longitudinal height data, in 787 participants with height measurements spanning the whole period of growth. Heritability estimates suggested that genetic factors could explain 25% to 71% of the variance of pubertal growth traits. We performed a GWAS of growth parameters, testing their associations with 5 077 595 imputed or directly genotyped variants. Six variants associated with height at peak velocity (P < 5 × 10-8, adjusted for sex, birth year and principal components). Implicated genes include NUDT3, previously associated with adult height, and PACSIN1. Two novel variants associated with duration of growth spurt (P < 5 × 10-8) in LOC105375344, an uncharacterized gene with unknown function. We finally examined the association of growth parameters with a polygenic score for height derived from 9557 single nucleotide polymorphisms (SNPs) identified in the GIANT meta-analysis for which genotypic data were available for the American Indian study population. Height polygenic score was correlated with the magnitude and velocity of height growth that occurred before and at the peak of the adolescent growth spurt, indicating overlapping genetic architecture, with no influence on the timing of adolescent growth.

Insulin resistance before type 2 diabetes onset is associated with increased risk of albuminuria after diabetes onset: A prospective cohort study.

AIM: Reduced renal insulin signalling is implicated in the pathogenesis of albuminuria. We sought to investigate whether insulin action and secretion, measured before diabetes onset, are associated with the development of albuminuria after diabetes onset. MATERIALS AND METHODS: Baseline body composition, insulin sensitivity by hyperinsulinaemic-euglycaemic clamp at submaximal and maximal insulin stimulation (240 and 2400 pmol/m2/min; M-low and M-high), and insulin secretion by intravenous glucose tolerance test [acute insulin response (AIR)] were measured in 170 Southwestern Indigenous American adults who subsequently developed diabetes. After diabetes onset and during the median follow-up of 13.6 years, 81 participants (48%) developed albuminuria (urine albumin-to-creatinine ratio ≥30 mg/g). Separate associations of M-low, M-high and AIR (per 1-SD change) with the risk of albuminuria were assessed by Cox regression models adjusted for age, sex and body fat (%). RESULTS: Participants who developed albuminuria were of similar age (26.4 ± 5.4 vs. 27.5 ± 6.1 years), sex (46% vs. 48% male), body fat (36.4 ± 7.5 vs. 35.7 ± 7.9%) and AIR [2.3 ± 0.3 vs. 2.3 ± 0.3, pmol/L (log)] as those who did not develop albuminuria but had lower insulin sensitivity [M-low: 0.33 ± 0.08 vs. 0.36 ± 0.12, p = .03; M-high: 0.87 ± 0.11 vs. 0.91 ± 0.12, p = .02; mg/kg-metabolic body size/min (log)]. In separate adjusted models, lower M-low and M-high were both associated with an increased risk for albuminuria [hazard ratio (HR) 1.51, 95% confidence interval (CI) 1.14, 2.00, p = .004; HR 1.31, 95% CI 1.06, 1.63, p = .01), whereas AIR was not (HR 1.15, 95% CI 0.87, 1.56, p = .3). CONCLUSIONS: Lower insulin sensitivity is associated with the development of albuminuria, suggesting a role for insulin signalling in the pathogenesis of proteinuria.

Pleiotropic effects of an eQTL in the CELSR2/PSRC1/SORT1 cluster that associates with LDL-C and resting metabolic rate.

CONTEXT: The locus CELSR2-PSRC1-SORT1, a primary genetic signal for lipids, has recently been implicated in different metabolic processes. Our investigation identified its association with energy metabolism. OBJECTIVE: To determine biological mechanisms that govern diverse functions of this locus. METHODS: Genotypes for 491,265 variants in 7,000 clinically characterized American Indians were previously determined using a custom-designed array specific for this longitudinally studied American Indian population. Among the genotyped individuals, 5,205 had measures of fasting lipid levels and 509 had measures of resting metabolic rate (RMR) and substrate oxidation rate assessed through indirect calorimetry. A genome-wide association study (GWAS) for LDL-C levels identified a variant in CELSR2 and the molecular impact of this variant on gene expression was assessed in skeletal muscle biopsies from 207 participants, followed by functional validation in mouse myoblasts using a luciferase assay. RESULTS: A GWAS in American Indians identified rs12740374 in CELSR2 as the top signal for LDL-C levels (P = 1 × 10-22); further analysis of this variant identified an unexpected correlation with reduced RMR (effect = -44.3 kcal/day/minor-allele) and carbohydrate oxidation rate (effect = -5.21 mg/hour/kg-EMBS). Tagged variants showed a distinct linkage disequilibrium architecture in American Indians, highlighting a potential functional variant, rs6670347 (minor-allele frequency = 0.20). Positioned in the glucocorticoid receptor's core binding motif, rs6670347 is part of a skeletal muscle-specific enhancer. Human skeletal muscle transcriptome analysis showed CELSR2 as the most differentially expressed gene (P = 1.9 × 10-7), with the RMR-lowering minor allele elevating gene expression. Experiments in mouse myoblasts confirmed enhancer-based regulation of CELSR2 expression, dependent on glucocorticoids. Rs6670347 also associated with increased oxidative phosphorylation gene expression; CELSR2 as a regulator of these genes, suggests potential influence on energy metabolism through muscle oxidative capacity. CONCLUSION: Variants in the CELSR2/PSRC1/SORT1 locus exhibit tissue-specific effects on metabolic traits, with an independent role in muscle metabolism through glucocorticoid signaling.

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.

An E115A Missense Variant in CERS2 Is Associated With Increased Sleeping Energy Expenditure and Hepatic Insulin Resistance in American Indians.

Genetic determinants of interindividual 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-h 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 postabsorptive sleeping metabolic rate.

Trans-ancestral genome-wide association study of longitudinal pubertal height growth and shared heritability with adult health outcomes.

BACKGROUND: Pubertal growth patterns correlate with future health outcomes. However, the genetic mechanisms mediating growth trajectories remain largely unknown. Here, we modeled longitudinal height growth with Super-Imposition by Translation And Rotation (SITAR) growth curve analysis on ~ 56,000 trans-ancestry samples with repeated height measurements from age 5 years to adulthood. We performed genetic analysis on six phenotypes representing the magnitude, timing, and intensity of the pubertal growth spurt. To investigate the lifelong impact of genetic variants associated with pubertal growth trajectories, we performed genetic correlation analyses and phenome-wide association studies in the Penn Medicine BioBank and the UK Biobank. RESULTS: Large-scale growth modeling enables an unprecedented view of adolescent growth across contemporary and 20th-century pediatric cohorts. We identify 26 genome-wide significant loci and leverage trans-ancestry data to perform fine-mapping. Our data reveals genetic relationships between pediatric height growth and health across the life course, with different growth trajectories correlated with different outcomes. For instance, a faster tempo of pubertal growth correlates with higher bone mineral density, HOMA-IR, fasting insulin, type 2 diabetes, and lung cancer, whereas being taller at early puberty, taller across puberty, and having quicker pubertal growth were associated with higher risk for atrial fibrillation. CONCLUSION: We report novel genetic associations with the tempo of pubertal growth and find that genetic determinants of growth are correlated with reproductive, glycemic, respiratory, and cardiac traits in adulthood. These results aid in identifying specific growth trajectories impacting lifelong health and show that there may not be a single "optimal" pubertal growth pattern.

Effects of the ABCC8 R1420H loss-of-function variant on beta-cell function, diabetes incidence, and retinopathy.

INTRODUCTION: The ABCC8 gene regulates insulin secretion and plays a critical role in glucose homeostasis. The effects of an ABCC8 R1420H loss-of-function variant on beta-cell function, incidence of type 2 diabetes, and age-at-onset, prevalence, and progression of diabetes complications were assessed in a longitudinal study in American Indians. RESEARCH DESIGN AND METHODS: We analyzed beta-cell function through the relationship between insulin secretion and insulin sensitivity in members of this population without diabetes aged ≥5 years using standard major axis regression. We used hierarchical logistic regression models to study cross-sectional associations with diabetes complications including increased albuminuria (albumin-to-creatinine ratio (ACR) ≥30 mg/g), severe albuminuria (ACR ≥300 mg/g), reduced estimated glomerular filtration rate (eGFR <60 mL/min/1.73 m2), and retinopathy. This study included 7675 individuals (254 variant carriers) previously genotyped for the R1420H with available phenotypic data and with a median follow-up time of 13.5 years (IQR 4.5-26.8). RESULTS: Variant carriers had worse beta-cell function than non-carriers (p=0.0004; on average estimated secretion was 22% lower, in carriers), in children and adults, with no difference in insulin sensitivity (p=0.50). At any body mass index and age before 35 years, carriers had higher type 2 diabetes incidence. This variant did not associate with prevalence of increased albuminuria (OR 0.87, 95% CI 0.66 to 1.16), severe albuminuria (OR 0.96, 95% CI 0.55 to 1.68), or reduced eGFR (OR 0.44, 95% CI 0.18 to 1.06). By contrast, the variant significantly associated with higher retinopathy prevalence (OR 1.74, 95% CI 1.19 to 2.53) and this association was only partially mediated (<11%) by glycemia, duration of diabetes, risk factors of retinopathy, or insulin use. Retinopathy prevalence in carriers was higher regardless of diabetes presence. CONCLUSIONS: The ABCC8 R1420H variant is associated with increased risks of diabetes and of retinopathy, which may be partially explained by higher glycemia levels and worse beta-cell function.