DOI QR코드

DOI QR Code

Genomic signatures of selection reveal genetic mechanisms underlying economic traits in Licha black pigs

  • Jiajia Liu (Key Laboratory of Animal Reproduction and Biotechnology in Universities of Shandong, College of Animal Science and Technology, Qingdao Agricultural University) ;
  • Zhe Tian (Key Laboratory of Animal Reproduction and Biotechnology in Universities of Shandong, College of Animal Science and Technology, Qingdao Agricultural University) ;
  • Mubin Yu (Key Laboratory of Animal Reproduction and Biotechnology in Universities of Shandong, College of Animal Science and Technology, Qingdao Agricultural University) ;
  • Wenwen Li (Key Laboratory of Animal Reproduction and Biotechnology in Universities of Shandong, College of Animal Science and Technology, Qingdao Agricultural University) ;
  • Pengcheng Lv (Key Laboratory of Animal Reproduction and Biotechnology in Universities of Shandong, College of Animal Science and Technology, Qingdao Agricultural University) ;
  • Tao Wang (Key Laboratory of Animal Reproduction and Biotechnology in Universities of Shandong, College of Animal Science and Technology, Qingdao Agricultural University) ;
  • Yu Tian (Key Laboratory of Animal Reproduction and Biotechnology in Universities of Shandong, College of Animal Science and Technology, Qingdao Agricultural University) ;
  • Shuer Zhang (General Station of Animal Husbandry of Shandong Province) ;
  • Junjie Wang (Key Laboratory of Animal Reproduction and Biotechnology in Universities of Shandong, College of Animal Science and Technology, Qingdao Agricultural University) ;
  • Wei Shen (Key Laboratory of Animal Reproduction and Biotechnology in Universities of Shandong, College of Animal Science and Technology, Qingdao Agricultural University)
  • Received : 2025.09.25
  • Accepted : 2025.11.26
  • Published : 2026.05.01

Abstract

Objective: The Licha black pig (LC) is a Chinese indigenous breed that is nationally protected, recognized for its superior meat quality and strong environmental adaptability. However, the population has experienced a rapid decline due to extensive crossbreeding with commercial lines. Understanding the genetic basis of economically important traits is crucial for conservation and genomic improvement. Methods: Whole-genome resequencing was performed on 120 LC pigs and combined with genomic data from 285 pigs representing 32 global populations, including wild boars, commercial breeds, and other Chinese indigenous pigs. Population structure was investigated using phylogenetic trees, principal component analysis, ADMIXTURE and TreeMix analysis. Selection signatures were identified through four complementary approaches (FST, θπ ratio, XP-CLR, and Tajima's D). Candidate genes were examined through functional enrichment analysis, protein structure prediction, and cross-referencing with trait association and tissue-specific expression databases. Phenotypic data on body size and teat number were also collected in LC pigs for targeted genotype-phenotype analysis. Results: Phylogenetic analyses showed clear stratification among global pig populations, with Chinese indigenous breeds significantly separated by the Qinling-Huaihe Line. LC pigs formed a distinct genetic cluster between northern Chinese and European breeds. Selective sweep analyses identified several candidate genes under positive selection, including SOCS6 and ATP2B4 (skeletal muscle development), RASAL2 (adipogenesis), and DOCK2 (male fertility). Trait-focused analyses identified ZNRF3 as a major locus for body size, with a missense mutation (g.46228935G>A; Gln→Arg) predicted to influence Wnt/β-catenin signaling. Signals of selection in ADGRB3, a gene potentially involved in teat patterning and mammary gland vascularization, were associated with variation in teat number. Conclusion: Our genomic analyses provide new insights into the genetic architecture of economically important traits and environmental adaptation in the LC. These findings provide a foundation for genomic-informed selective breeding and present valuable molecular tools for the genetic improvement and sustainable utilization of this indigenous genetic resource.

Keywords

Acknowledgement

This work was financially supported by the Key R&D Program of Shandong Province, China (2024LZGC002), the Key R&D Program of Shandong Province, China (2024LZGC017), Qingdao Science and Technology for the People Demonstration Project (24-1-8-xdny-5-nsh), Shandong Modern Agricultural Industry Technology System - Swine Innovation Team (Project ID: SDAIT-08-02).

References

  1. Wang T, Liu QC, Tian Z, et al. Selective sweep analysis reveals candidate genes related to immunity in Licha black pig. Gene 2025;963:149616. https://doi.org/10.1016/j.gene.2025.149616
  2. Xie Q, Zhang Z, Chen Z, et al. Integration of selection signatures and protein interactions reveals NR6A1, PAPPA2, and PIK3C2B as the promising candidate genes underlying the characteristics of Licha black pig. Biology 2023;12:500. https://doi.org/10.3390/biology12040500
  3. Wang Y, Dong R, Li X, Cui C, Yu G. Analysis of the genetic diversity and family structure of the Licha black pig population on Jiaodong Peninsula, Shandong province, China. Animals 2022;12:1045. https://doi.org/10.3390/ani12081045
  4. Zhang X, Li C, Li X, et al. Association analysis of polymorphism in the NR6A1 gene with the lumbar vertebrae number traits in sheep. Genes Genom 2019;41:1165-71. https://doi.org/10.1007/s13258-019-00843-5
  5. Wang Y, Wang M, He X, Dong R, Liu H, Yu G. Detection and analysis of genes affecting the number of thoracic vertebrae in Licha black pigs. Genes 2024;15:477. https://doi.org/10.3390/genes15040477
  6. Li X, Lu Y, Liu X, Xie X, Wang K, Yu D. Identification of chicken FSHR gene promoter and the correlations between polymorphisms and egg production in Chinese native hens. Reprod Domest Anim 2019;54:702-11. https://doi.org/10.1111/rda.13412
  7. Pecka-Kiełb E, Kowalewska-Łuczak I, Czerniawska-Piątkowska E, Króliczewska B. FASN, SCD1 and ANXA9 gene polymorphism as genetic predictors of the fatty acid profile of sheep milk. Sci Rep 2021;11:23761. https://doi.org/10.1038/s41598-021-03186-y
  8. Rohrer GA, Nonneman DJ, Miller RK, Zerby H, Moeller SJ. Association of single nucleotide polymorphism (SNP) markers in candidate genes and QTL regions with pork quality traits in commercial. Meat Sci 2012;92:511-8. https://doi.org/10.1016/j.meatsci.2012.05.020
  9. de Sena Brandine G, Smith AD. Falco: high-speed FastQC emulation for quality control of sequencing data. F1000Res 2019;8:1874. https://doi.org/10.12688/f1000research.21142.1
  10. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 2014;30:2114-20. https://doi.org/10.1093/bioinformatics/btu170
  11. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009;25:1754-60. https://doi.org/10.1093/bioinformatics/btp324
  12. Li H, Handsaker B, Wysoker A, et al. The sequence alignment/map format and SAMtools. Bioinformatics 2009;25:2078-9. https://doi.org/10.1093/bioinformatics/btp352
  13. McKenna A, Hanna M, Banks E, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res 2010;20:1297-303. https://doi.org/10.1101/gr.107524.110
  14. Danecek P, Auton A, Abecasis G, et al. The variant call format and VCFtools. Bioinformatics 2011;27:2156-8. https://doi.org/10.1093/bioinformatics/btr330
  15. Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet 2011;88:76-82. https://doi.org/10.1016/j.ajhg.2010.11.011
  16. Wickham H. ggplot2. Wiley Interdiscip Rev Comput Stat 2011;3:180-5. https://doi.org/10.1002/wics.147
  17. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol 2018;35:1547-9. https://doi.org/10.1093/molbev/msy096
  18. Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res 2009;19:1655-64. https://doi.org/10.1101/gr.094052.109
  19. Purcell S, Neale B, Todd-Brown K, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007;81:559-75. https://doi.org/10.1086/519795
  20. Zhang C, Dong SS, Xu JY, He WM, Yang TL. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics 2019;35:1786-8. https://doi.org/10.1093/bioinformatics/bty875
  21. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res 2010;38:e164. https://doi.org/10.1093/nar/gkq603
  22. Raudvere U, Kolberg L, Kuzmin I, et al. g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res 2019;47:W191-8. https://doi.org/10.1093/nar/gkz369
  23. Pertea G, Pertea M. GFF utilities: GffRead and GffCompare. F1000Res 2020;9:304. https://doi.org/10.12688/f1000research.23297.2
  24. Abramson J, Adler J, Dunger J, et al. Accurate structure prediction of biomolecular interactions with AlphaFold 3. Nature 2024;630:493-500. https://doi.org/10.1038/s41586-024-07487-w
  25. DeLano WL. The PyMOL molecular graphics system. DeLano Scientific LLC; 2002.
  26. Wang Z, Sun H, Chen Q, Zhang X, Wang Q, Pan Y. A genome scan for selection signatures in Taihu pig breeds using next-generation sequencing. Animal 2019;13:683-93. https://doi.org/10.1017/S1751731118001714
  27. Du L, Chang T, An B, et al. Transcriptome profiling analysis of muscle tissue reveals potential candidate genes affecting water holding capacity in Chinese Simmental beef cattle. Sci Rep 2021;11:11897. https://doi.org/10.1038/s41598-021-91373-2
  28. Schuh K, Cartwright EJ, Jankevics E, et al. Plasma membrane Ca2+ ATPase 4 is required for sperm motility and male fertility. J Biol Chem 2004;279:28220-6. https://doi.org/10.1074/jbc.M312599200
  29. Singh R, Junghare V, Hazra S, et al. Database on spermatozoa transcriptogram of catagorised Frieswal crossbred (Holstein Friesian X Sahiwal) bulls. Theriogenology 2019;129:130-45. https://doi.org/10.1016/j.theriogenology.2019.01.025
  30. Zhu JJ, Ramanathan P, Bishop EA, O'Donnell V, Gladue DP, Borca MV. Mechanisms of African swine fever virus pathogenesis and immune evasion inferred from gene expression changes in infected swine macrophages. PLOS ONE 2019;14:e0223955. https://doi.org/10.1371/journal.pone.0223955
  31. Pan Y, Tong JHM, Lung RWM, et al. RASAL2 promotes tumor progression through LATS2/YAP1 axis of hippo signaling pathway in colorectal cancer. Mol Cancer 2018;17:102. https://doi.org/10.1186/s12943-018-0853-6
  32. Zhu X, Xie S, Xu T, Wu X, Han M. Rasal2 deficiency reduces adipogenesis and occurrence of obesity-related disorders. Mol Metab 2017;6:494-502. https://doi.org/10.1016/j.molmet.2017.03.003
  33. Deng S, Qiu Y, Zhuang Z, et al. Genome-wide association study of body conformation traits in a three-way crossbred commercial pig population. Animals 2023;13:2414. https://doi.org/10.3390/ani13152414
  34. Liu S, Chen S, Cai W, et al. Divergence analyses of sperm DNA methylomes between monozygotic twin AI bulls. Epigenomes 2019;3:21. https://doi.org/10.3390/epigenomes3040021
  35. Serrano M, Ramón M, Calvo JH, et al. Genome-wide association studies for sperm traits in Assaf sheep breed. Animal 2021;15:100065. https://doi.org/10.1016/j.animal.2020.100065
  36. Szenker-Ravi E, Altunoglu U, Leushacke M, et al. RSPO2 inhibition of RNF43 and ZNRF3 governs limb development independently of LGR4/5/6. Nature 2018;557:564-9. https://doi.org/10.1038/s41586-018-0118-y
  37. Nazar M, Abdalla IM, Chen Z, et al. Genome-wide association study for udder conformation traits in Chinese Holstein cattle. Animals 2022;12:2542. https://doi.org/10.3390/ani12192542
  38. Qi L, Xiao L, Fu R, Nie Q, Zhang X, Luo W. Genetic characteristics and selection signatures between Southern Chinese local and commercial chickens. Poult Sci 2024;103:103863. https://doi.org/10.1016/j.psj.2024.103863
  39. Taheri S, Javadmanesh A, Zerehdaran S. Identification of selective sweep and associated QTL traits in Iranian Ovis aries and Ovis orientalis populations. Front Genet 2024;15:1414717. https://doi.org/10.3389/fgene.2024.1414717
  40. Gangwar M, Ahmad SF, Ali AB, et al. Identifying low-density, ancestry-informative SNP markers through whole genome resequencing in Indian, Chinese, and wild yak. BMC Genomics 2024;25:1043. https://doi.org/10.1186/s12864-024-10924-9