A Genomewide Integrative Analysis of GWAS and eQTLs Data Identifies Multiple Genes and Gene Sets Associated with Obesity

To identify novel susceptibility genes and gene sets for obesity, we conducted a genomewide expression association analysis of obesity via integrating genomewide association study (GWAS) and expression quantitative trait loci (eQTLs) data. GWAS summary data of body mass index (BMI) and waist-to-hip ratio (WHR) was driven from a published study, totally involving 339,224 individuals. The eQTLs dataset (containing 927,753 eQTLs) was obtained from eQTLs meta-analysis of 5,311 subjects. Integrative analysis of GWAS and eQTLs data was conducted by SMR software. The SMR single gene analysis results were further subjected to gene set enrichment analysis (GSEA) for identifying obesity associated gene sets. A total of 13,311 annotated gene sets were analyzed in this study. SMR single gene analysis identified 20 BMI associated genes (TUFM, SPI1, APOB48R, etc.). Also 3 WHR associated genes were detected (CPEB4, WARS2, and L3MBTL3). The significant association between Chr16p11 and BMI was observed by GSEA (FDR adjusted p value = 0.040). The TGCTGCT, MIR-15A, MIR-16, MIR-15B, MIR-195, MIR-424, and MIR-497 (FDR adjusted p value = 0.049) gene set appeared to be linked with WHR. Our results provide novel clues for the genetic mechanism studies of obesity. This study also illustrated the good performance of SMR for susceptibility gene mapping.