Identification and functional characterization of G6PC2 coding variants influencing glycemic traits define an effector transcript at the G6PC2-ABCB11 locus.
Mahajan A., Sim X., Ng HJ., Manning A., Rivas MA., Highland HM., Locke AE., Grarup N., Im HK., Cingolani P., Flannick J., Fontanillas P., Fuchsberger C., Gaulton KJ., Teslovich TM., Rayner NW., Robertson NR., Beer NL., Rundle JK., Bork-Jensen J., Ladenvall C., Blancher C., Buck D., Buck G., Burtt NP., Gabriel S., Gjesing AP., Groves CJ., Hollensted M., Huyghe JR., Jackson AU., Jun G., Justesen JM., Mangino M., Murphy J., Neville M., Onofrio R., Small KS., Stringham HM., Syvänen A-C., Trakalo J., Abecasis G., Bell GI., Blangero J., Cox NJ., Duggirala R., Hanis CL., Seielstad M., Wilson JG., Christensen C., Brandslund I., Rauramaa R., Surdulescu GL., Doney ASF., Lannfelt L., Linneberg A., Isomaa B., Tuomi T., Jørgensen ME., Jørgensen T., Kuusisto J., Uusitupa M., Salomaa V., Spector TD., Morris AD., Palmer CNA., Collins FS., Mohlke KL., Bergman RN., Ingelsson E., Lind L., Tuomilehto J., Hansen T., Watanabe RM., Prokopenko I., Dupuis J., Karpe F., Groop L., Laakso M., Pedersen O., Florez JC., Morris AP., Altshuler D., Meigs JB., Boehnke M., McCarthy MI., Lindgren CM., Gloyn AL., T2D-GENES consortium and GoT2D consortium None.
Genome wide association studies (GWAS) for fasting glucose (FG) and insulin (FI) have identified common variant signals which explain 4.8% and 1.2% of trait variance, respectively. It is hypothesized that low-frequency and rare variants could contribute substantially to unexplained genetic variance. To test this, we analyzed exome-array data from up to 33,231 non-diabetic individuals of European ancestry. We found exome-wide significant (P<5×10-7) evidence for two loci not previously highlighted by common variant GWAS: GLP1R (p.Ala316Thr, minor allele frequency (MAF)=1.5%) influencing FG levels, and URB2 (p.Glu594Val, MAF = 0.1%) influencing FI levels. Coding variant associations can highlight potential effector genes at (non-coding) GWAS signals. At the G6PC2/ABCB11 locus, we identified multiple coding variants in G6PC2 (p.Val219Leu, p.His177Tyr, and p.Tyr207Ser) influencing FG levels, conditionally independent of each other and the non-coding GWAS signal. In vitro assays demonstrate that these associated coding alleles result in reduced protein abundance via proteasomal degradation, establishing G6PC2 as an effector gene at this locus. Reconciliation of single-variant associations and functional effects was only possible when haplotype phase was considered. In contrast to earlier reports suggesting that, paradoxically, glucose-raising alleles at this locus are protective against type 2 diabetes (T2D), the p.Val219Leu G6PC2 variant displayed a modest but directionally consistent association with T2D risk. Coding variant associations for glycemic traits in GWAS signals highlight PCSK1, RREB1, and ZHX3 as likely effector transcripts. These coding variant association signals do not have a major impact on the trait variance explained, but they do provide valuable biological insights.