Genome-Wide Association Studies Revealed Significant QTLs and Candidate Genes Associated with Backfat and Loin Muscle Area in Pigs Using Imputation-Based Whole Genome Sequencing Data

SIMPLE SUMMARY: Backfat thickness (BF) and loin muscle area (LMA) are important characteristics in pig breeding. The identification of the quantitative trait loci (QTLs) and genes closely associated with these carcass traits may be useful for pig molecular breeding. In this study, genome-wide association studies (GWAS) using imputation-based whole genome sequencing data of four selected phenotypic traits (adjusted 100 kg BF and LMA, adjusted 100 kg BF estimate breeding values and LMA estimate breeding values) were conducted in a total of 1131 pigs from three pure breeds (French Yorkshire, Landrace, and Duroc). As a result, we highlighted CCND2 and SHANK2 genes as strong candidates affecting BF traits. The results of this study may provide additional insight into the genetic contribution of these genes to carcass meat production capabilities in pig breeding and enhancing lean meat percentages in pigs. ABSTRACT: Improvement of carcass features is an essential goal in pig genetic breeding programs. Backfat (BF) and loin muscle area (LMA) are important carcass production metrics and useful indicators of pig production performance and lean meat rate. However, the genetic architecture of BF and LMA traits remains elusive. To identify quantitative trait loci (QTLs) and genes associated with these traits, we performed a genome-wide association study (GWAS) using imputation-based whole genome sequencing data for four phenotypes (adjusted 100 kg BF and LMA, adjusted 100 kg BF EBV and LMA EBV) in 1131 pigs from 3 breeds (French Yorkshire, Landrace, and Duroc). After genotype imputation and quality control, 14,163,315 single nucleotide polymorphisms (SNPs) were retained for further analysis. For the adjusted 100 kg BF trait, using the 2-LOD drop method, a QTL with a 13.4 Kb interval (2.91 to 2.93 Mb on SSC2) and containing a SHANK2 gene was defined. In addition, two QTLs with 135.40 Kb (from 66.10 to 66.23 Mb) and 3.12 Kb (from 66.886 to 66.889 Mb) intervals containing CCND2 and TSPAN11 genes, respectively, were found on SSC5. For the BF-EBV trait, two QTLs (128.77 Kb from 66.10 to 66.23 Mb on SSC5 and 42.10 Kb from 2.89 to 2.93 Mb on SSC2) were identified. Notably, CCND2 and SHANK2 were the only candidate genes in their respective QTL interval. Furthermore, we detected a 3.33 Kb (66.106 to 66.110 Mb on SSC2) haplotype block which was detected as affecting the BF_EBV trait, which only contained the CCND2 gene. Thus, we suggested CCND2 and SHANK2 as strong candidate genes for regulating the BF trait for pigs. The empirical confidence intervals of the QTLs were 1.14 Mb (165.65 to 166.79 Mb on SSC6) for adjusted 100 kg LMA and 1.49 Mb (165.26–166.74 Mb on SSC6) for LMA-EBV. These two confidence intervals contained 13 and 28 annotated genes, respectively. Our results provide a deeper understanding of the genetic basis of pig carcass traits. The identified molecular markers will be useful for selecting breeding lines for breeding pigs with superior carcass traits.