Integrating Multi-Omics Data to Identify Key Functional Variants Affecting Feed Efficiency in Large White Boars

Abstract

Optimizing feed efficiency through the feed conversion ratio (FCR) is paramount for economic viability and sustainability. In this study, we integrated RNA-seq, ATAC-seq, and genome-wide association study (GWAS) data to investigate key functional variants associated with feed efficiency in pigs. Identification of differentially expressed genes in the duodenal and muscle tissues of low- and high-FCR pigs revealed that pathways related to digestion of dietary carbohydrate are responsible for differences in feed efficiency between individuals. Differential open chromatin regions identified by ATAC-seq were linked to genes involved in glycolytic and fatty acid processes. GWAS identified 211 significant single-nucleotide polymorphisms associated with feed efficiency traits, with candidate genes PPP1R14C, TH, and CTSD. Integration of duodenal ATAC-seq data and GWAS data identified six key functional variants, particularly in the 1500985–1509676 region on chromosome 2. In those regions, CTSD was found to be highly expressed in the duodenal tissues of pigs with a high feed conversion ratio, suggesting its role as a potential target gene. Overall, the integration of multi-omics data provided insights into the genetic basis of feed efficiency, offering valuable information for breeding more efficient pig breeds.