Search

Despite the successes of human genome-wide association studies, the causal genes underlying most metabolic traits remain unclear. We used outbred heterogeneous stock (HS) rats, coupled with expression data and mediation analysis, to identify quantitative trait loci (QTLs) and candidate gene mediators for adiposity, glucose tolerance, serum lipids, and other metabolic traits. Physiological traits were measured in 1,519 male HS rats, with liver and adipose transcriptomes measured in >410 rats. Genotypes were imputed from low-coverage whole-genome sequencing. Linear mixed models were used to detect physiological and expression QTLs (pQTLs and eQTLs, respectively), using both single nucleotide polymorphism (SNP)– and haplotype-based models for pQTL mapping. Genes with cis-eQTLs that overlapped pQTLs were assessed as causal candidates through mediation analysis. We identified 14 SNP-based pQTLs and 19 haplotype-based pQTLs, of which 10 were in common. Using mediation, we identified the following genes as candidate mediators of pQTLs: Grk5 for fat pad weight and serum triglyceride pQTLs on Chr1, Krtcap3 for fat pad weight and serum triglyceride pQTLs on Chr6, Ilrun for a fat pad weight pQTL on Chr20, and Rfx6 for a whole pancreatic insulin content pQTL on Chr20. Furthermore, we verified Grk5 and Ktrcap3 using gene knockdown/out models, thereby shedding light on novel regulators of obesity.

Decline in insulin sensitivity due to dysfunction of adipose tissue (AT) is one of the earliest pathogenic events in type 2 diabetes. We hypothesize that differential DNA methylation (DNAm) controls insulin sensitivity and obesity by modulating transcript expression in AT. Integrating AT DNAm profiles with transcript profile data measured in a cohort of 230 African Americans (AAs) from the African American Genetics of Metabolism and Expression cohort, we performed cis-expression quantitative trait methylation (cis-eQTM) analysis to identify epigenetic regulatory loci for glucometabolic trait–associated transcripts. We identified significantly associated cytosine-guanine dinucleotide regions for 82 transcripts (false discovery rate [FDR]-P < 0.05). The strongest eQTM locus was observed for the proopiomelanocortin (POMC; ρ = −0.632, P = 4.70 × 10−27) gene. Epigenome-wide association studies (EWAS) further identified 155, 46, and 168 cytosine-guanine dinucleotide regions associated (FDR-P < 0.05) with the Matsuda index, SI, and BMI, respectively. Intersection of EWAS, transcript level to trait association, and eQTM results, followed by causal inference test identified significant eQTM loci for 23 genes that were also associated with Matsuda index, SI, and/or BMI in EWAS. These associated genes include FERMT3, ITGAM, ITGAX, and POMC. In summary, applying an integrative multiomics approach, our study provides evidence for DNAm-mediated regulation of gene expression at both previously identified and novel loci for many key AT transcripts influencing insulin resistance and obesity.

Decline in insulin sensitivity due to dysfunction of adipose tissue (AT) is one of the earliest pathogenic events in type 2 diabetes. We hypothesize that differential DNA methylation (DNAm) in regulatory regions of a subset of genes controls insulin sensitivity and obesity by modulating transcript expression in AT. Integrating reduced representation bisulfite sequencing-derived AT DNAm profiles for 1,054,719 CpG sites and 200,800 CpG regions with AT transcript profile data in 230 African Americans, we performed cis-expression quantitative trait methylation (cis-eQTM) analysis focused on 6,774 glucometabolic trait-associated transcripts to identify local epigenetic regulatory loci for these transcripts. We identified significantly associated CpG regions for 140 transcripts (FDR-P In summary, applying an integrative multi-omics approach, our study provides evidence for DNAm-mediated regulation of gene expression at both previously identified and novel loci for many key AT transcripts influencing insulin resistance and obesity. Disclosure N.K. Sharma: None. M. Comeau: None. D.J. Montoya: None. M. Pellegrini: None. T. Howard: None. C.D. Langefeld: None. S.K. Das: None. Funding American Diabetes Association (1-18-ICTS-113 to S.K.D.); National Institutes of Health (R01DK090111)

Insulin resistance (IR) is a harbinger of type 2 diabetes (T2D) and partly determined by genetic factors. However, genetically regulated mechanisms of IR remain poorly understood. Using gene expression, genotype, and insulin sensitivity data from the African American Genetics of Metabolism and Expression (AAGMEx) cohort, we performed transcript-wide correlation and expression quantitative trait loci (eQTL) analyses to identify IR-correlated cis-regulated transcripts (cis-eGenes) in adipose tissue. These IR-correlated cis-eGenes were tested in the European ancestry individuals in the Metabolic Syndrome in Men (METSIM) cohort for trans-ethnic replication. Comparison of Matsuda index–correlated transcripts in AAGMEx with the METSIM study identified significant correlation of 3,849 transcripts, with concordant direction of effect for 97.5% of the transcripts. cis-eQTL for 587 Matsuda index–correlated genes were identified in both cohorts. Enoyl-CoA hydratase domain-containing 3 (ECHDC3) was the top-ranked Matsuda index–correlated cis-eGene. Expression levels of ECHDC3 were positively correlated with Matsuda index, and regulated by cis-eQTL, rs34844369 being the top cis-eSNP in AAGMEx. Silencing of ECHDC3 in adipocytes significantly reduced insulin-stimulated glucose uptake and Akt Ser473 phosphorylation. RNA sequencing analysis identified 691 differentially expressed genes in ECHDC3-knockdown adipocytes, which were enriched in γ-linolenate biosynthesis, and known IR genes. Thus, our studies elucidated genetic regulatory mechanisms of IR and identified genes and pathways in adipose tissue that are mechanistically involved in IR.