Whole Genome Sequence Association Analysis of Fasting Glucose and Fasting Insulin Levels in Diverse Cohorts from the NHLBI TOPMed Program

The genetic determinants of fasting glucose (FG) and fasting insulin (FI) have been studied mostly through genome and exome arrays, resulting in over 100 associated variants. We extended this work with a high-coverage whole genome sequencing (WGS) analysis from fifteen cohorts in the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. More than 23,000 non-diabetic individuals from five self-reported race/ethnicities (African, Asian, European, Hispanic and Samoan) were included for each trait. We analyzed 60M variants in race/ethnicity-specific and pooled single variant and rare variant aggregate tests. Twenty-two variants across sixteen gene regions were found significantly associated with FG or FI, eight of which were rare (Minor Allele Frequency, MAFG6PC2locus we identified a distinct FG signal at rare variant rs2232326 (MAF=0.01) after conditioning on known common variants. Functional annotations show rs2232326 to be disruptive and likely damaging while being weakly transcribed in islets. A pair of FG-associated variants were identified near theSLC30A8locus. These variants, one of which was rare (MAF=0.001) and Asian race/ethnicity-specific, were shown to be in islet-specific active enhancer regions. Other associated regions include rare variants nearROBO1andPTPRT, and common variants nearMTNR1B, GCK, GCKR, FOXA2, APOB, TCF7L2, andADCY5. We provide a catalog of nucleotide-resolution genomic variation spanning intergenic and intronic regions down to a minor allele count of 20, creating a foundation for future sequencing-based investigation of glycemic traits.