Exploring genomic regions regulating the liver transcriptome and energy homeostasis in pigs

In pigs, energy homeostasis has an impact on meat quality and health. In a Duroc pig population, 30 quantitative trait locus (QTL) regions associated with fatty acid (FA) composition in adipose tissue, plasma, liver and muscle were previously identified. Mapping of expression quantitative trait locus (eQTL) regions will provide a molecular hypothesis for genotype-phenotype interactions and may allow the identification of shared causal variants, key to increasing our understanding of the genetic regulation of FA composition and energy homeostasis. However, gene expression is impacted by environmental factors, while individual-level allelic imbalance (AI) can be more reliable and can be surveyed via allelic-specific expression (ASE) analysis. Furthermore, treatment of ASE as a quantitative trait allows the identification of allele-specific expression quantitative trait loci (aseQTLs), which are variants whose heterozygosity is linked to the AI of a nearby single-nucleotide polymorphism (SNP), pointing to regulatory elements. In this study, liver was selected as a key metabolic hub with an important role in the regulation of energy homeostasis, and 310 liver RNA sequencing samples were analysed using a combination of (1) eQTL mapping, (2) ASE analysis, and (3) aseQTL mapping methods. A total of 2 188 eQTL regions were identified, mostly cis-eQTL regions (73.17%). ASE analysis reported 1 964 ASE SNPs, associated with 633 genes. Finally, aseQTL mapping reported 64 172 aseQTL, associated with the AI of 31 genes. Colocalisation analysis combined with ASE analysis showed that the expression of FADS1 and FADS2 genes is associated with the polyunsaturated FA composition in several tissues, where microRNA regulation may be present. Finally, in the DGAT2 gene, annotated in a QTL region associated with multiple FAs in adipose tissue, ASE revealed allelic imbalance in the 3' untranslated region (UTR) of this gene. Allelic differential expression can be caused by a 13-bp insertion affecting messenger RNA stability, previously described, exemplifying how allelic imbalance is caused by a post-transcriptional regulatory mechanism undetectable by eQTL mapping. Furthermore, aseQTLs were associated with this gene, linked to a previously identified copy number variant not yet associated with DGAT2 expression. These results demonstrate how ASE analysis and aseQTL mapping can complement eQTL mapping, as they resolved a complex region affected by allelic heterogeneity, a main confounding effect of QTL mapping. In conclusion, the combination of eQTL mapping, ASE analysis and aseQTL mapping allowed the characterisation of the regulation of liver gene expression, improving our understanding of the genetic determinism of energy homeostasis.