Abstract
Backfat thickness (BFT) and feed conversion ratio (FCR) are important commercial traits in the pig industry. With the increasing demand for human health and meat production, identifying functional genomic regions and genes associated with these commercial traits is critical for enhancing production efficiency. In this research, we conducted a genome-wide association study (GWAS) on a Landrace population comprising 4,295 individuals with chip data for BFT and FCR. Our analysis revealed a total of 118 genome-wide significant signals located on chromosomes SSC1, SSC2, SSC7, SSC12, and SSC13, respectively. Furthermore, we identified 10 potential regions associated with the two traits and annotated the genes within these regions. In addition, enrichment analysis was also performed. Notably, candidate genes such as SHANK2, KCNQ1, and ABL1 were found to be associated with BFT, whereas NAP1L4, LSP1, and PPFIA1 genes were related to the FCR. Our findings provide valuable insights into the genetic architecture of these two traits and offer guidance for future pig breeding efforts.
Introduction
The increasing demand for human health and food nutrition has become a challenge due to the rapid growth of global population. Consequently, enhancing the production efficiency of livestock products has become vital for the livestock industry and sustainable development (). Swine is one of the most important economic livestock in the world, providing a diverse range of products to meet human needs. The rapid development of breeding methods, such as genomic selection, has effectively reduced the genetic interval between pig generations and significantly improved the performance of commercial pig breeds by increasing prediction accuracy (). As consumers’ demand for healthier meat products increases, pigs have been bred for lower fat content and higher lean meat. Previous research uncovered that daily energy intake is related to whole-body fat composition in male pigs (), and leaner pigs tend to exhibit higher feed efficiency. Therefore, understanding the genetic architecture of these commercial traits is essential.
Feed conversion ratio (FCR) and backfat thickness (BFT) are primary commercial phenotypes in the pig industry and have been extensively analyzed by numerous researchers. Candidate genes, such as phospholipase A2 group IB (PLA2G1B), have been reported to be associated with feed efficiency by influencing lipid catabolism (; ). Additionally, the members of the insulin-like growth factor family, such as IGF1 and IGF2, have been found to affect the growth rate and feed conversion efficiency (). Backfat thickness is another important trait in pig production as it impacts lean meat yield and the popularity of pork meat (). Many loci on SSC1, SSC5, SSC6, SSC7, and SSC12, as well as candidate genes such as MC4R, IGF2, and LEPR, were found to be related to backfat thickness (; ).
Over the past 15 years, genome-wide association studies (GWASs) have been employed to investigate the linkage between genomic markers and records of various traits (). This approach has facilitated the identification of numerous quantitative trait nucleotides (QTNs) and candidate genes associated with FCR and BFT (; ; ; ), which provided deep insights into these commercial traits and improved the quality of meat production. To date, about 55,688 quantitative trait loci (QTLs) have been released by pig QTLdb (). However, due to the complexity of these quantitative traits, many QTLs remain unknown.
Using a Landrace population with genomic chip data, a total of 4,295 individuals with two important commercial traits, including BFT at 100 kg and FCR, were analyzed in this study. Related variants with annotated candidate genes within candidate regions were detected using a mixed-effects linear model in a genome-wide association study. Furthermore, linkage disequilibrium (LD) block analysis with candidate regions and enrichment analysis of candidate genes were also performed. The main objectives of this research were to identify the associated genomic regions and candidate genes of BFT and FCR within our population. In addition, we also used a multi-omics swine database () to prioritize the candidate genes in order to provide an understanding of the majority of candidate genes.
Results
Summary of phenotype and genotype data
In this research, we used a Landrace population consisting of 4,295 individuals with chip-level genotype data derived from a functional SNP lipid chip. Summary statistics, including the sample size, mean of phenotype values, standard deviation of phenotypes, and coefficient of variation, are provided in Table 1. Additionally, phenotype distribution plots are shown in Supplementary Figure 1. The mean values for BFT (100 kg) and FCR in our population are 10.22 and 2.3 with standard deviations of 1.99 and 0.22 and the coefficients of variation are 0.2 and 0.09, respectively. According to Table 1 and Supplementary Figure 1, both traits can be used for further analysis. Quality control was performed on the genomic data. After imputation and filtering, 100,235 SNPs on autosomes remained for the association study, excluding those with a minor allele frequency (MAF) less than 0.01. The marker density on each chromosome is shown in Figure 1A. LD decay analysis was also conducted using PopLDdecay (), and the results are shown in Figure 1B.
TABLE 1
Summary statistics of phenotypes.
SD, standard deviation of each phenotype.
CV, coefficient of variation for each phenotype.
FIGURE 1
Genome-wide association studies
Genome-wide association studies were performed through a mixed-effects linear model in the rMVP package (). Sex, farms, and the first three principle components were included as fixed effects, whereas the additive genetic effects were considered the random effects variable. A Bonferroni cutoff of 0.05/N was used as a significant threshold, where N represents the number of SNPs. The details of quality control and genotype data are described in Methods. A total of 69 significant SNPs distributed on SSC1, SSC2, SSC7, and SSC12 were identified to be associated with backfat thickness. The details of these SNPs are shown in Supplementary Table 1. Manhattan and quantile–quantile (QQ) plots for BFT (100 kg) are shown in Figure 2, with a lambda value of 0.93 indicating minimal population inflation in GWAS. Based on the LD decay results, a distance of 300 kb was determined to define the candidate regions around each significant signal, and the candidate regions are shown in Table 2. Regions with overlapping areas were merged into one region, and a total of 244 genes were annotated within these candidate regions. After annotation, the genes were prioritized by a multi-omics database called ISwine (). Further details are provided in Supplementary Table 3.
FIGURE 2
TABLE 2
| Phenotype | Chr | Candidate region | Candidate gene |
|---|---|---|---|
| BFT100 kg | 1 | 160,473,437–161,073,437 | MC4R |
| 1 | 270,269,333–270,926,968 | ABL1 and FIBCD1 | |
| 2 | 0–2,457,364 | ENSSSCG00000035293, SLC22A18, PTDSS2, HRAS, PNPLA2, TALDO1, BRSK2, and DUSP8 | |
| 2 | 2,758,419–3,644,831 | ANO1, FGF19, FGF3, FGF4, SHANK2, KCNQ1, and INS | |
| 7 | 97,275,068–98,039,684 | AREL1 and VRTN | |
| 12 | 50,709,260–53,168,040 | ARRB2, GLTPD2, TNFSF12, SLC16A13, and TP53 | |
| FCR | 1 | 262,745,504–263,345,504 | ENSSSCG00000035556, ENSSSCG00000031416 |
| 2 | 133,461–733,461 | RNH1, TSPAN4, AP2A2, and RASSF7 | |
| 2 | 945,702–3,644,831 | KCNQ1, ANO1, SHANK2, LSP1, PPFIA1, NAP1L4, IFITM10, TNNI2, PRR33, TNNT3, MOB2, INS, SYT8, ASCL2, CTSD, TH, and DUSP8 | |
| 13 | 44,899,972–45,499,972 | SYNPR |
Summary of the identified candidate regions and genes of phenotypes.
A total of 49 significant markers located on SSC1, SSC2, and SSC13, respectively, were found to be related to the FCR. Manhattan and QQ plots of FCR GWAS results are presented in Figure 3, where the lambda value was 0.96 for FCR GWAS results. The candidate regions and genes for the FCR are provided in Table 2, and the details of gene annotation in the candidate regions are displayed in Supplementary Table 4.
FIGURE 3
Gene ontology annotation analysis and enrichment results
The candidate genes were annotated with Gene Ontology (GO) and the Kyoto Encyclopedia for Genes and Genomes (KEGG) database using IAnimal (). The details for the gene annotation results of BFT are shown in Supplementary Figure 2 and Supplementary Table 5. A total of 161 GO terms were significantly enriched, comprising 94 biological processes (BPs), 21 cellular components (CCs), and 46 molecular functions (MFs). Notably, several significant GO terms were associated with backfat thickness, such as positive regulation of insulin secretion (P = 0.019), intermembrane lipid transfer (P = 0.028), and lipid transfer activity (P = 0.041). Significant pathways, such as the MAPK signaling pathway (P = ), were also found in KEGG results. Genes such as anoctamin 1 (ANO1), glycolipid transfer protein domain containing 2 (GLTPD2), TNF superfamily member 12 (TNFSF12), fibroblast growth factor 19 (FGF19), fibroblast growth factor 3 (FGF3), fibroblast growth factor 4 (FGF4), and ENSSSCG00000035293 were involved in these significant terms and pathways.
Enrichment results for the FCR are shown in Supplementary Figure 3 and Supplementary Table 6. A total of 37 KEGG pathways were significantly enriched (P < 0.05). We also found 132 significant GO terms, including 94 BPs, 12 CCs, and 26 MFs. The positive regulation of insulin secretion, involved in the cellular response to glucose stimuli, is associated with nutrient absorption and energy metabolism, which were significantly enriched (P = 0.0057), and anoctamin 1 (ANO1) was involved in this process. Furthermore, some digestion-related pathways, such as protein digestion and absorption, were detected with suggestive P-values. Potassium voltage-gated channel subfamily Q member 1 (KCNQ1) was involved in this pathway and has been reported to be related to pig feed efficiency ().
Linkage disequilibrium block analysis
The LD blocks around the peak signals were analyzed and plotted using LDBlockShow (). These plots are shown in Figure 4, Supplementary Figure 4, and Supplementary Figure 5. Figure 4 highlights the most significant regions of GWAS results for the two traits. Multiple LD blocks were observed around the top signals, with an overlapping region between BFT and FCR results, spanning from 2.76 to 3.64 Mb on SSC2, indicating that this region may have an influence on both FCR and BFT.
FIGURE 4
Discussion
In our research, the GWAS approach was used to identify related genomic regions and candidate genes for backfat thickness and feed conversion ratio; a total of 4,295 Landraces were involved and analyzed. Our results showed 10 candidate regions on the genome for two commercial traits, and among the identified candidate genes, we found many reported genes that should be related to these two traits.
There were 46 overlapping signals between the GWAS results for BFT and FCR, suggesting that these variants may have influences on both traits. Some genes in these candidate regions were found to be related to both backfat thickness and feed conversion ratio. For example, as a member of the SHANK protein family, SHANK2 (SH3 and multiple ankyrin repeat domains 2) was reported to be associated with childhood obesity (). It was strongly highlighted as a candidate gene for backfat thickness in an association study based on imputed whole-genome data from a multi-breed population (). Interestingly, SHANK2 was found to be associated with average daily gain and the meat-to-fat ratio in pooled F2-designed pigs (). Since SHANK2 was also found as a candidate gene for the FCR, it may play an important role in both backfat thickness and feed conversion ratio.
The MAPK signaling pathway was significantly enriched in the KEGG analysis results for backfat thickness (P = ), which was reported to be essential for adipogenesis as it can regulate porcine fat deposition (; ; ). A total of nine genes, namely, FGF3, FGF4, FGF19, ARRB2, INS, HRAS, DUSP8, TP53, and ENSSSCG00000035293, were involved in this pathway. Many of these genes have been reported in previous studies. For example, INS was identified as a candidate gene for porcine backfat thickness (). It has been reported that the blood glucose level is regulated by insulin, which is encoded by the INS gene; it can also promote cell fat storage and affect lipid metabolism (). Additionally, DUSP8 was found to be related to ham weight loss, and HRAS was reported to be associated with backfat in pigs. Other genes related to fat metabolism, such as PTDSS2, TALDO1, and BRSK2, have been identified in recent studies (; ; ). ABL1 was detected as a candidate gene for backfat traits in Yorkshire and Duroc populations (; ). It was also found to be related to the meat-to-fat ratio in pigs (). A recent GWAS in the Yorkshire population also confirmed that ABL1 was associated with the average daily gain (). The QTL on SSC1 (270 Mb) was similar to the FIBCD1 gene, which was also found in a Swiss Large White pig population and associated with the body mass index in humans (; ). SLC16A13 and PNPLA2 were found in the candidate region on SSC2. SLC16A13 is a candidate gene for diabetes in mice, and its deletion will attenuate lipid accumulation and insulin resistance (). SLC22A18 is also related to lipids; knocking down SLC22A18 in mice will reduce hepatic lipid accumulation, revealing its positive effects on lipid accumulation (). PNPLA2 promoted lipid accumulation in an adipogenesis test in pigs ().
MC4R is a major gene influencing fatness in pigs. It is also involved in the regulation of feeding behavior and body weight in mice and humans. A missense mutation in this gene leads to increased fat accumulation in pigs (). MC4R has also been shown to affect growth, feed intake, and backfat thickness in pigs, according to previous studies by , , , and . We also found that AREL1 and VRTN on SSC7, which have previously been associated with body length, teat number, and intramuscular fat content (; ; ; ), are also associated with meat production traits, thus influencing backfat.
Pigs are known for their outstanding olfactory abilities, which are attributed to their abundant functional olfactory receptors (). Odors affect pig reproduction and also have an influence on early food preferences (). In the region of 262,745,504–263,345,504 on SSC1, we detected ENSSSCG00000035556 and ENSSSCG00000031416 to be associated with the FCR. These genes were enriched in an MF term of olfactory receptor activity (P = 0.015). In addition to these two genes, nine other genes (ENSSSCG00000036003, ENSSSCG00000032805, ENSSSCG00000032825, ENSSSCG00000027589, ENSSSCG00000037454, ENSSSCG00000027732, ENSSSCG00000031516, ENSSSCG00000026287, and ENSSSCG00000035439) were also enriched in this term and may be related. However, based on the prioritization results, ENSSSCG00000035556 and ENSSSCG00000031416 had higher scores, suggesting that these two genes may have a more significant impact on the FCR. The KEGG pathway analysis also identified significant enrichment in an olfactory-related pathway, namely, olfactory transduction (P = 0.046). ANO1 (anoctamin 1) was part of this pathway, along with 11 other genes enriched for olfactory receptor activity, as mentioned earlier. Olfactory receptors can perceive odor, transform biochemical signaling events into electrical impulses, and send them to the brain (). This pathway was also found to be significantly enriched and related to residual feed intake (). Therefore, ANO1 may affect the FCR through olfactory-related pathways. ANO1 also affects fat deposition in pigs ().
Feed efficiency can be largely influenced by energy metabolism and digestion. KCNQ1 (potassium voltage-gated channel subfamily Q member 1) was the most significantly related gene among these candidate genes. It plays an important role in fetal development in mice and humans, and it was reported to be associated with development. KCNQ1 is an imprinted gene expressed paternally in pigs () that influences nutrient absorption by regulating gastric acid secretion, as well as salt and glucose homeostasis (). It is also found to be involved in a related KEGG pathway that facilitates the digestion and absorption of proteins. In a recent multi-omics study, KCNQ1, along with SYT8, TNNI2, ASCL2, MOB2, DUSP8, TH, PRR33, TNNT3, IFITM10, CTSD, and INS, was identified as a candidate gene for the FCR in Large White boars (). These genes were also observed in our results. Interestingly, KCNQ1 was investigated in an association study as one of the candidate genes in the backfat thickness of pigs (). The LSP1 gene has been reported to be associated with weight loss after dry-curing of hams (). Additionally, it can stimulate myogenic factors and influence skeletal muscle development in pigs (). Since skeletal muscle plays a key role in energy storage and consumption and is closely linked to energy metabolism, LSP1 may impact muscle development, thereby affecting pig growth and FCR. TSPAN4 and AP2A2 were also reported to be associated with growth in pigs with LSP1 (). In addition to LSP1, PPF1A1 was found to be differently expressed in the top three canonical pathways (). In an epigenome-wide skeletal muscle study, NAP1L4 was detected to have CpG positions hypermethylated within its promoters (). RASSF7 is a candidate gene for the FCR and a member of the N-terminal Ras association domain family. Studies have shown that knocking down RASSF7 restricts cell growth (), and a deficiency in lysine negatively impacts the expression of RASSF7 (). A mitochondrial protein, RNH1, was found to be associated with angiogenesis in porcine corpus luteum (). Synaptoporin (SYNPR) is one of the tetratransmembrane transport vesicle proteins, which is distributed in the digestive system (). It has been found to be genome-wide associated with autoimmune hepatitis (AIH) in humans (). We identified SYNPR as one of the significant candidate genes on SSC13. Therefore, it may affect feed efficiency through digestion.
In conclusion, we performed a GWAS based on a Landrace population to investigate two commercial traits. Our analysis identified 118 significant signals, from which 10 candidate regions were selected. Candidate genes within these regions were annotated and further analyzed using GO and KEGG pathways. Among the identified candidate genes, MC4R, SLC22A13, and INS were associated with backfat thickness, whereas ENSSSCG00000035556, SHANK2, KCNQ1, and LSP1 were related to the feed conversion ratio. Overall, our research provides deeper insights into the genetic basis of these traits and could inform future pig breeding efforts.
Materials and methods
Collection of Landrace population and phenotypes
A total of 4,295 Landrace individuals were used in this research, collected from three different great-grandparent farms of COFCO Joycome Foods Co., Ltd. All the pigs were raised under uniform feeding and management standards during the measurement period. Original records, including daily feed intake and weight for each pig, were automatically collected using the Pig Performance Testing System (Nedap, Groenlo, Netherlands). Outliers in these records were removed. The start and end dates of the test, along with the initial and final weights, weight gain, and total feed consumed during the test, were recorded. The FCR was then calculated as the total feed intake divided by the weight gain. The backfat thickness was measured using living B-ultrasonography at the end of the test, and the measured traits were then adjusted to a body weight of 100 kg. Details and distribution of phenotypes are shown in Table 1 and Supplementary Figure 1.
Genotyping, imputation, and quality control
Genomic DNA was extracted from ear tissue samples using the Tecan Freedom EVO NGS Workstation and the MagPure Tissue DNA KF Kit (MD5112-02), with a concentration of ≥40 ng/µL and a quantity of ≥1 µg. After that, the samples were genotyped using the Porcine 80K functional SNP genotyping chip by Wuhan Yingzi Gene Co., Ltd., using target capture sequencing technology. A total of 187,255 variants were included in the original genotype data. Quality control was performed using PLINK 1.90 (). Markers with a call rate <90% and MAF < 0.01 and variants on the sex chromosomes were excluded. This resulted in the retention of 100,240 variants, with a total genotyping rate of 0.99. Beagle 5.4 () was applied to impute genotype data. After imputation, quality control was performed again to remove markers with MAF < 0.01. A total of 100,235 SNPs and 4,295 individuals were left for further analysis at last.
Genome-wide association study
GWAS analysis was performed by fitting a mixed-effects linear model using the following equation in the rMVP () package:where represents the phenotypic values for FCR and BFT (100 kg); represents the fixed effects, including the first three principal components, sex, and farm; represents the SNP vector being tested; represents the vector of random effects; and represents the vector of residual errors. In this study, the additive genetic relationship matrix is imputed by the genomic relationship matrix , derived from all SNP variants used in the association test (), and , , and are the incidence matrices for , , and , respectively. A significant threshold of 0.05/N was confirmed by a Bonferroni correction, where N represents the number of SNPs. Manhattan and QQ plots were constructed in the R environment using rMVP.
Identification of candidate regions, genes, and enrichment analysis
LD decay analysis was performed to detect the size of the candidate region using PopLDdecay (). A window size of 300 kb was determined based on the LD decay results. Significant SNPs were then annotated to nearby genes within a 300 kb upstream or downstream range and prioritized using the multi-omics swine knowledgebase ISwine (). The candidate regions were uploaded in the section “Search by Region” on ISwine to obtain a candidate gene list. Then, we used “Prioritize” in the tool section to prioritize the candidate genes for each commercial trait and downloaded the results. Enrichment analyses, including GO and KEGG pathway analyses, were performed using enrichment tools from the IAnimal database (). The LD blocks around the top QTN plot were generated by LDBlockShow ().
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://figshare.com/, https://figshare.com/articles/dataset/genotype_data_of_Landrace_GWAS/27861177.
Ethics statement
The animal studies were approved by Breeding Department, Wuhan COFCO Meat Co., Ltd. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author contributions
GM: conceptualization, writing–review and editing, formal analysis, and resources. XT: formal analysis, data curation, and writing–review and editing. YY: writing–review and editing. TZ: data curation, formal analysis, and writing–review and editing. JW: data curation, formal analysis, and writing–review and editing. XC: data curation, formal analysis, and writing–review and editing. JX: conceptualization, investigation, methodology, writing–original draft, and writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by STI2030-Major Projects [2023ZD0404707].
Acknowledgments
The authors would like to thank Shanxia Yang, Zijia Liu, and Jingya Zhang for the suggestions and help in this study.
Conflict of interest
Authors GM, XT, YY, TZ, JW, XC, and JX were employed by Wuhan COFCO Meat Co., Ltd.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2024.1505197/full#supplementary-material
SUPPLEMENTARY TABLE 1Details of GWAS results for BFT.
SUPPLEMENTARY TABLE 2Details of GWAS results for FCR.
SUPPLEMENTARY TABLE 3Details of annotated genes for BFT.
SUPPLEMENTARY TABLE 4Details of annotated genes for FCR.
SUPPLEMENTARY TABLE 5Details of GO and KEGG results for BFT.
SUPPLEMENTARY TABLE 6Details of GO and KEGG results for FCR.
References
1
AbdellaouiA.YengoL.VerweijK. J.VisscherP. M. (2023). 15 years of GWAS discovery: realizing the promise. Am. J. Hum. Genet.110, 179–194. 10.1016/j.ajhg.2022.12.011
2
AlbuquerqueA.ÓviloC.NúñezY.BenítezR.López-GarciaA.GarcíaF.et al (2021). Transcriptomic profiling of skeletal muscle reveals candidate genes influencing muscle growth and associated lipid composition in Portuguese local pig breeds. Animals11, 1423. 10.3390/ani11051423
3
AouadiM.LaurentK.ProtM.Le Marchand-BrustelY.BinétruyB.BostF. (2006). Inhibition of p38MAPK increases adipogenesis from embryonic to adult stages. Diabetes55, 281–289. 10.2337/diabetes.55.02.06.db05-0963
4
BlajI.TetensJ.PreußS.BennewitzJ.ThallerG. (2018). Genome-wide association studies and meta-analysis uncovers new candidate genes for growth and carcass traits in pigs. PloS one13, e0205576. 10.1371/journal.pone.0205576
5
BrowningB. L.ZhouY.BrowningS. R. (2018). A one-penny imputed genome from next-generation reference panels. Am. J. Hum. Genet.103, 338–348. 10.1016/j.ajhg.2018.07.015
6
BrunjesP. C.FeldmanS.OsterbergS. K. (2016). The pig olfactory brain: a primer. Chem. senses41, 415–425. 10.1093/chemse/bjw016
7
ComuzzieA. G.ColeS. A.LastonS. L.VorugantiV. S.HaackK.GibbsR. A.et al (2012). Novel genetic loci identified for the pathophysiology of childhood obesity in the Hispanic population. PloS one7, e51954. 10.1371/journal.pone.0051954
8
DelpuechE.AliakbariA.LabruneY.FèveK.BillonY.GilbertH.et al (2021). Identification of genomic regions affecting production traits in pigs divergently selected for feed efficiency. Genet. Sel. Evol.53, 49. 10.1186/s12711-021-00642-1
9
DingR.ZhuangZ.QiuY.RuanD.WuJ.YeJ.et al (2022). Identify known and novel candidate genes associated with backfat thickness in Duroc pigs by large-scale genome-wide association analysis. J. Animal Sci.100, skac012. 10.1093/jas/skac012
10
DoD. N.StratheA. B.OstersenT.PantS. D.KadarmideenH. N. (2014). Genome-wide association and pathway analysis of feed efficiency in pigs reveal candidate genes and pathways for residual feed intake. Front. Genet.5, 307. 10.3389/fgene.2014.00307
11
DongS.-S.HeW.-M.JiJ.-J.ZhangC.GuoY.YangT.-L. (2020). LDBlockShow: a fast and convenient tool for visualizing linkage disequilibrium and haplotype blocks based on variant call format files. Briefings Bioinforma.22, bbaa227. 10.1093/bib/bbaa227
12
FaggionS.BonfattiV.CarnierP. (2024). Genome-wide association study for weight loss at the end of dry-curing of hams produced from purebred heavy pigs. Animals14, 1983. 10.3390/ani14131983
13
Falker-GieskeC.BlajI.PreußS.BennewitzJ.ThallerG.TetensJ. (2019). GWAS for meat and carcass traits using imputed sequence level genotypes in pooled F2-designs in pigs. G3 Genes, Genomes, Genet.9, 2823–2834. 10.1534/g3.119.400452
14
FuL.JiangY.WangC.MeiM.ZhouZ.JiangY.et al (2020a). A genome-wide association study on feed efficiency related traits in Landrace pigs. Front. Genet.11, 692. 10.3389/fgene.2020.00692
15
FuM.ZhouX.LiuZ.WangT.LiuB. (2023a). Genome‐wide association study reveals a genomic region on SSC7 simultaneously associated with backfat thickness, skin thickness and carcass length in a Large White× Tongcheng advanced generation intercross resource population. Anim. Genet.54, 216–219. 10.1111/age.13285
16
FuY.LiuH.DouJ.WangY.LiaoY.HuangX.et al (2023b). IAnimal: a cross-species omics knowledgebase for animals. Nucleic Acids Res.51, D1312–D1324. 10.1093/nar/gkac936
17
FuY.XuJ.TangZ.WangL.YinD.FanY.et al (2020b). A gene prioritization method based on a swine multi-omics knowledgebase and a deep learning model. Commun. Biol.3, 502. 10.1038/s42003-020-01233-4
18
GalveA.BurgosC.SilióL.VaronaL.RodríguezC.OviloC.et al (2012). The effects of leptin receptor (LEPR) and melanocortin-4 receptor (MC4R) polymorphisms on fat content, fat distribution and fat composition in a Duroc× Landrace/Large White cross. Livest. Sci.145, 145–152. 10.1016/j.livsci.2012.01.010
19
Gozalo-MarcillaM.BuntjerJ.JohnssonM.BatistaL.DiezF.WernerC. R.et al (2021). Genetic architecture and major genes for backfat thickness in pig lines of diverse genetic backgrounds. Genet. Sel. Evol.53, 76. 10.1186/s12711-021-00671-w
20
GroenenM. A.ArchibaldA. L.UenishiH.TuggleC. K.TakeuchiY.RothschildM. F.et al (2012). Analyses of pig genomes provide insight into porcine demography and evolution. Nature491, 393–398. 10.1038/nature11622
21
HiroseK.MikawaS.OkumuraN.NoguchiG.FukawaK.KanayaN.et al (2013). Association of swine vertnin (VRTN) gene with production traits in D uroc pigs improved using a closed nucleus breeding system. Animal Sci. J.84, 213–221. 10.1111/j.1740-0929.2012.01066.x
22
HoaV. B.SeoH. W.SeongP. N.ChoS. H.KangS. M.KimY. S.et al (2021). Back‐fat thickness as a primary index reflecting the yield and overall acceptance of pork meat. Animal Sci. J.92, e13515. 10.1111/asj.13515
23
HollieN. I.HuiD. Y. (2011). Group 1B phospholipase A₂ deficiency protects against diet-induced hyperlipidemia in mice. J. lipid Res.52, 2005–2011. 10.1194/jlr.M019463
24
HongY.YeJ.DongL.LiY.YanL.CaiG.et al (2021). Genome-wide association study for body length, body height, and total teat number in large white pigs. Front. Genet.12, 650370. 10.3389/fgene.2021.650370
25
HuZ.-L.DrachevaS.JangW.MaglottD.BastiaansenJ.RothschildM. F.et al (2005). A QTL resource and comparison tool for pigs: PigQTLDB. Mamm. Genome16, 792–800. 10.1007/s00335-005-0060-9
26
KnolE. F.NielsenB.KnapP. W. (2016). Genomic selection in commercial pig breeding. Anim. Front.6, 15–22. 10.2527/af.2016-0003
27
LeeJ.KimY.ChoE.ChoK.SaS.KimY.et al (2020). Genomic analysis using Bayesian methods under different genotyping platforms in Korean Duroc pigs. Animals10, 752. 10.3390/ani10050752
28
LeeJ.-B.KangH.-C.KimE.-H.KimY.-J.YooC.-K.ChoiT.-J.et al (2018). Genome-wide association study identifies positional candidate genes affecting back fat thickness trait in pigs. Korean J. Agric. Sci.45, 707–713. 10.7744/kjoas.20180055
29
LiD.HuangM.ZhuangZ.DingR.GuT.HongL.et al (2021). Genomic analyses revealed the genetic difference and potential selection genes of growth traits in two duroc lines. Front. Veterinary Sci.8, 725367. 10.3389/fvets.2021.725367
30
LiJ.WuJ.JianY.ZhuangZ.QiuY.HuangR.et al (2022b). 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. Animals12, 2911. 10.3390/ani12212911
31
LiW.WangZ.LuoS.WuJ.ZhouL.LiuJ. (2022a). Genome-wide association analysis and genetic parameters for feed efficiency and related traits in Yorkshire and Duroc pigs. Animals12, 1902. 10.3390/ani12151902
32
LiY.SunY.LiuY.WangB.LiJ.WangH.et al (2022c). Genome‐wide meta‐analysis identifies susceptibility loci for autoimmune hepatitis type 1. Hepatology76, 564–575. 10.1002/hep.32417
33
LikszoP.SkarzynskiD. J.Moza JalaliB. (2021). Changes in porcine corpus luteum proteome associated with development, maintenance, regression, and rescue during estrous cycle and early pregnancy. Int. J. Mol. Sci.22, 11740. 10.3390/ijms222111740
34
LiuF.BrewsterC. J.GilmourS. L.HenmanD. J.SmitsR. J.LuxfordB. G.et al (2021). Relationship between energy intake and growth performance and body composition in pigs selected for low backfat thickness. J. Animal Sci.99, skab342. 10.1093/jas/skab342
35
LiuL.YaoX.WangY.HuR.FanC.GongH.et al (2022). Physins in digestive system neoplasms. Adv. Clin. Chem.111, 157–176. 10.1016/bs.acc.2022.08.002
36
LiuX.DuY.TrakooljulN.BrandB.MurániE.KrischekC.et al (2015). Muscle transcriptional profile based on muscle fiber, mitochondrial respiratory activity, and metabolic enzymes. Int. J. Biol. Sci.11, 1348–1362. 10.7150/ijbs.13132
37
MaH.ZhangS.ZhangK.ZhanH.PengX.XieS.et al (2019). Identifying selection signatures for backfat thickness in Yorkshire pigs highlights new regions affecting fat metabolism. Genes10, 254. 10.3390/genes10040254
38
MaM. (2007). Encoding olfactory signals via multiple chemosensory systems. Crit. Rev. Biochem. Mol. Biol.42, 463–480. 10.1080/10409230701693359
39
MehrabiZ.GillM.WijkM. v.HerreroM.RamankuttyN. (2020). Livestock policy for sustainable development. Nat. Food1, 160–165. 10.1038/s43016-020-0042-9
40
MiaoY.ZhaoY.WanS.MeiQ.WangH.FuC.et al (2023). Integrated analysis of genome-wide association studies and 3D epigenomic characteristics reveal the BMP2 gene regulating loin muscle depth in Yorkshire pigs. PLoS Genet.19, e1010820. 10.1371/journal.pgen.1010820
41
NoskováA.MehrotraA.KadriN.NeuenschwanderS.HoferA.PauschH. (2022). Detecting QTL for two lowly correlated traits using multi-trait and meta-analysis approaches in Swiss Large White pigs. Proc. 12th World Congr. Genet. Appl. Livest. Prod. (WCGALP) Tech. species orientated innovations animal Breed. contribution Genet. solving Soc. challenges, 1630–1633.
42
OviloC.FernándezA.RodríguezM.NietoM.SilióL. (2006). Association of MC4R gene variants with growth, fatness, carcass composition and meat and fat quality traits in heavy pigs. Meat Sci.73, 42–47. 10.1016/j.meatsci.2005.10.016
43
ParkJ. (2024). Genome-wide association study to reveal new candidate genes using single-step approaches for productive traits of Yorkshire pig in Korea. Anim. Biosci.37, 451–460. 10.5713/ab.23.0255
44
ParkJ.DoK. T.ParkK. D.LeeH. K. (2023). Genome‐wide association study using a single‐step approach for teat number in Duroc, Landrace and Yorkshire pigs in Korea. Anim. Genet.54, 743–751. 10.1111/age.13357
45
PiórkowskaK.TyraM.RogozM.Ropka-MolikK.OczkowiczM.RóżyckiM. (2010). Association of the melanocortin-4 receptor (MC4R) with feed intake, growth, fatness and carcass composition in pigs raised in Poland. Meat Sci.85, 297–301. 10.1016/j.meatsci.2010.01.017
46
PonsuksiliS.TrakooljulN.BasavarajS.HadlichF.MuraniE.WimmersK. (2019). Epigenome-wide skeletal muscle DNA methylation profiles at the background of distinct metabolic types and ryanodine receptor variation in pigs. BMC Genomics20, 492. 10.1186/s12864-019-5880-1
47
PulitS. L.StonemanC.MorrisA. P.WoodA. R.GlastonburyC. A.TyrrellJ.et al (2019). Meta-analysis of genome-wide association studies for body fat distribution in 694 649 individuals of European ancestry. Hum. Mol. Genet.28, 166–174. 10.1093/hmg/ddy327
48
PurcellS.NealeB.Todd-BrownK.ThomasL.FerreiraM. A.BenderD.et al (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet.81, 559–575. 10.1086/519795
49
RecinoA.SherwoodV.FlaxmanA.CooperW. N.LatifF.WardA.et al (2010). Human RASSF7 regulates the microtubule cytoskeleton and is required for spindle formation, Aurora B activation and chromosomal congression during mitosis. Biochem. J.430, 207–213. 10.1042/BJ20100883
50
SaltielA. R.KahnC. R. (2001). Insulin signalling and the regulation of glucose and lipid metabolism. Nature414, 799–806. 10.1038/414799a
51
SchumannT.KönigJ.von LoeffelholzC.VatnerD. F.ZhangD.PerryR. J.et al (2021). Deletion of the diabetes candidate gene Slc16a13 in mice attenuates diet-induced ectopic lipid accumulation and insulin resistance. Commun. Biol.4, 826. 10.1038/s42003-021-02279-8
52
ShiL.WangL.FangL.LiM.TianJ.WangL.et al (2022). Integrating genome-wide association studies and population genomics analysis reveals the genetic architecture of growth and backfat traits in pigs. Front. Genet.13, 1078696. 10.3389/fgene.2022.1078696
53
SunJ.MacKinnonR. (2020). Structural basis of human KCNQ1 modulation and gating. Cell180, 340–347. 10.1016/j.cell.2019.12.003
54
VanRadenP. M. (2008). Efficient methods to compute genomic predictions. J. Dairy Sci.91, 4414–4423. 10.3168/jds.2007-0980
55
WangH.WangX.YanD.SunH.ChenQ.LiM.et al (2022). Genome-wide association study identifying genetic variants associated with carcass backfat thickness, lean percentage and fat percentage in a four-way crossbred pig population using SLAF-seq technology. BMC Genomics23, 594. 10.1186/s12864-022-08827-8
56
WangT.FeugangJ. M.CrenshawM. A.RegmiN.BlantonJ. R.LiaoS. F. (2017). A systems biology approach using transcriptomic data reveals genes and pathways in porcine skeletal muscle affected by dietary lysine. Int. J. Mol. Sci.18, 885. 10.3390/ijms18040885
57
WangX.SunB.JiangQ.WuR.CaiM.YaoY.et al (2018). mRNA m6A plays opposite role in regulating UCP2 and PNPLA2 protein expression in adipocytes. Int. J. Obes.42, 1912–1924. 10.1038/s41366-018-0027-z
58
WuY.-Q.ZhaoH.LiY.-J.KhederzadehS.WeiH.-J.ZhouZ.-Y.et al (2020). Genome-wide identification of imprinted genes in pigs and their different imprinting status compared with other mammals. Zoological Res.41, 721–725. 10.24272/j.issn.2095-8137.2020.072
59
XiangY.SunJ.MaG.DaiX.MengY.FuC.et al (2024). Integrating multi-omics data to identify key functional variants affecting feed efficiency in large white boars. Genes15, 980. 10.3390/genes15080980
60
YamamotoT.IizukaY.Izumi-YamamotoK.ShirotaM.MoriN.TaharaY.et al (2024). Overexpression of Slc22a18 facilitates fat accumulation in mice. Biochem. Biophysical Res. Commun.712, 149922. 10.1016/j.bbrc.2024.149922
61
YangJ.HuangL.YangM.FanY.LiL.FangS.et al (2016). Possible introgression of the VRTN mutation increasing vertebral number, carcass length and teat number from Chinese pigs into European pigs. Sci. Rep.6, 19240. 10.1038/srep19240
62
YinL.ZhangH.TangZ.XuJ.YinD.ZhangZ.et al (2021). rMVP: a memory-efficient, visualization-enhanced, and parallel-accelerated tool for genome-wide association study. Genomics Proteomics Bioinforma.19, 619–628. 10.1016/j.gpb.2020.10.007
63
ZhangC.DongS. S.XuJ. Y.HeW. M.YangT. L. (2019). PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics35, 1786–1788. 10.1093/bioinformatics/bty875
64
ZhangD.WuW.HuangX.XuK.ZhengC.ZhangJ. (2021). Comparative analysis of gene expression profiles in differentiated subcutaneous adipocytes between Jiaxing Black and Large White pigs. BMC genomics22, 61–13. 10.1186/s12864-020-07361-9
65
ZhuL. H.XuJ. X.ZhuS. W.CaiX.YangS. F.ChenX. L.et al (2014). Gene expression profiling analysis reveals weaning-induced cell cycle arrest and apoptosis in the small intestine of pigs. J. Animal Sci.92, 996–1006. 10.2527/jas.2013-7551
Summary
Keywords
pigs, feed efficiency, backfat thickness, GWAS, candidate genes
Citation
Ma G, Tan X, Yan Y, Zhang T, Wang J, Chen X and Xu J (2025) A genome-wide association study identified candidate regions and genes for commercial traits in a Landrace population. Front. Genet. 15:1505197. doi: 10.3389/fgene.2024.1505197
Received
02 October 2024
Accepted
09 December 2024
Published
06 January 2025
Volume
15 - 2024
Edited by
Xiao-Lin Wu, Council on Dairy Cattle Breeding, United States
Updates
Copyright
© 2025 Ma, Tan, Yan, Zhang, Wang, Chen and Xu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jingya Xu, xujingya@cofco.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.