Data from: Susceptibility to type 2 diabetes may be modulated by haplotypes in G6PC2, a target of positive selection
Data files
Jan 27, 2017 version files 2.68 MB
-
G6PC_align
-
G6PC_chimpanzee_gammaMap.fa
-
G6PC_gorilla_gammaMap.fa
-
G6PC_human_gammaMap.fa
-
G6PC_outgroup.fa
-
G6PC2_align
-
G6PC2_chimpanzee_gammaMap.fa
-
G6PC2_gorilla_gammaMap.fa
-
G6PC2_human_gammaMap.fa
-
G6PC2_outgroup.fa
-
G6PC3_align
-
G6PC3_chimpanzee_gammaMap.fa
-
G6PC3_gorilla_gammaMap.fa
-
G6PC3_human_gammaMap.fa
-
G6PC3_outgroup.fa
-
Metazoans_align
-
plink.map
-
plink.ped
-
plinkCOVAR.txt
-
skatCOVAR
-
skatGENOall
-
skatPHENO
-
trees
Abstract
Background: The endoplasmic reticulum enzyme glucose-6-phosphatase catalyzes the common terminal reaction in the gluconeogenic/glycogenolytic pathways and plays a central role in glucose homeostasis. In most mammals, different G6PC subunits are encoded by three paralogous genes (G6PC, G6PC2, and G6PC3). Mutations in G6PC and G6PC3 are responsible for human mendelian diseases, whereas variants in G6PC2 are associated with fasting glucose (FG) levels. Results: We analyzed the evolutionary history of G6Pase genes. Results indicated that the three paralogs originated during early vertebrate evolution and that negative selection was the major force shaping diversity at these genes in mammals. Nonetheless, site-wise estimation of evolutionary rates at corresponding sites revealed weak correlations, suggesting that mammalian G6Pases have evolved different structural features over time. We also detected pervasive positive selection at mammalian G6PC2. Most selected residues localize in the C-terminal protein region, where several human variants associated with FG levels also map. This region was re-sequenced in ~560 subjects from Saudi Arabia, 185 of whom suffering from type 2 diabetes (T2D). The frequency of rare missense and nonsense variants was not significantly different in T2D and controls. Association analysis with two common missense variants (V219L and S342C) revealed a weak but significant association for both SNPs when analyses were conditioned on rs560887, previously identified in a GWAS for FG. Two haplotypes were significantly associated with T2D with an opposite effect direction. Conclusions: We detected pervasive positive selection at mammalian G6PC2 genes and we suggest that distinct haplotypes at the G6PC2 locus modulate susceptibility to T2D.