DOI QR코드

DOI QR Code

Whole-genome selection signatures identified candidate genes associated with cashmere traits in Inner Mongolia cashmere goats

  • Youjun Rong (College of Animal Science, Inner Mongolia Agricultural University) ;
  • Xiaofang Ao (College of Animal Science, Inner Mongolia Agricultural University) ;
  • Mingxuan Han (College of Animal Science, Inner Mongolia Agricultural University) ;
  • Qincheng Xia (College of Animal Science, Inner Mongolia Agricultural University) ;
  • Fangzheng Shang (College of Animal Science, Inner Mongolia Agricultural University) ;
  • Qi Lv (College of Animal Science, Inner Mongolia Agricultural University) ;
  • Zhiying Wang (College of Animal Science, Inner Mongolia Agricultural University) ;
  • Rui Su (College of Animal Science, Inner Mongolia Agricultural University) ;
  • Yanhong Zhao (College of Animal Science, Inner Mongolia Agricultural University) ;
  • Yanjun Zhang (College of Animal Science, Inner Mongolia Agricultural University) ;
  • Ruijun Wang (College of Animal Science, Inner Mongolia Agricultural University)
  • 투고 : 2025.04.10
  • 심사 : 2025.07.01
  • 발행 : 2025.12.01

초록

Objective: Inner Mongolia cashmere goats are superior indigenous breeds developed through long-term natural selection and systematic artificial selection, which have experienced a certain intensity of selection pressure during the breeding process, leading to bipolar differentiation trends in cashmere traits. Therefore, identifying genomic selection signatures associated with cashmere traits in Inner Mongolia cashmere goats is crucial for breeding high-quality cashmere-producing goats. Methods: To unravel the genetic basis of cashmere traits, this study stratified 375 Inner Mongolia cashmere goats into eight subgroups based on breeding values for cashmere traits: high-yield vs low-yield cashmere types (HYCG vs LYCG), fine vs coarse cashmere types (FCG vs CCG), long vs short cashmere types (LCG vs SCG), and long vs short fleece types (LFCG vs SFCG). Whole-genome resequencing was performed for genotyping, followed by detection of selection signatures. Results: Results revealed 144, 158, 147, and 147 high-frequency run of homozygosity (ROH) regions in HYCG, FCG, LCG, and LFCG subgroups, respectively, annotating to 515, 565, 510, and 521 genes. Additionally, genomic regions under positive selection were identified using FST, θπ ratios, and XP-EHH methods, with overlapping regions detected by ≥2 methods defined as candidate regions. Gene annotation identified 777, 660, 712, and 726 candidate genes in HYCG vs LYCG, FCG vs CCG, LCG vs SCG, and LFCG vs SFCG comparisons, respectively. These genes were enriched in 3,051 GO terms and 318 KEGG pathways, including Hippo, MAPK, Wnt, PI3K-Akt, and mTOR signaling pathways associated with cashmere growth and development, involving genes such as LGR6, RUNX2, IGF1R, FGF9, and TCF7L1. Conclusion: In this study, we employed four complementary approaches, including ROHs, FST, θπ ratios, and XP-EHH, to identify genomic signatures of selection for cashmere traits in Inner Mongolia cashmere goats. These findings provide valuable insights for improving cashmere production performance and developing novel strains with high-quality cashmere in Inner Mongolia cashmere goats.

키워드

과제정보

Thanks to all the staff of Erlangshan Ranch of Inner Mongolia Beiping Textile Co., Ltd. for their contributions in the production performance measurement.

참고문헌

  1. Gibson J, Morton NE, Collins A. Extended tracts of homozygosity in outbred human populations. Hum Mol Genet 2006; 15:789-95. https://doi.org/10.1093/hmg/ddi493
  2. Yin C, Wang Y, Zhou P, et al. Genomic scan for runs of homozygosity and selective signature analysis to identify candidate genes in large white pigs. Int J Mol Sci 2023;24:12914. https://doi.org/10.3390/ijms241612914
  3. Lukic B, Curik I, Drzaic I, et al. Genomic signatures of selection, local adaptation and production type characterisation of East Adriatic sheep breeds. J Anim Sci Biotechnol 2023;14: 142. https://doi.org/10.1186/s40104-023-00936-y
  4. Ramljak J, Špehar M, Ceranac D, et al. Genomic characterization of local croatian sheep breeds-effective population size, inbreeding & signatures of selection. Animals 2024;14: 1928. https://doi.org/10.3390/ani14131928
  5. Bao J, Xiong J, Huang J, Yang P, Shang M, Zhang L. Genetic diversity, selection signatures, and genome-wide association study identify candidate genes related to litter size in Hu sheep. Int J Mol Sci 2024;25:9397. https://doi.org/10.3390/ijms25179397
  6. Weir BS, Cockerham CC. Estimating F-statistics for the analysis of population structure. Evolution 1984;38:1358-70. https://doi.org/10.1111/j.1558-5646.1984.tb05657.x
  7. Sabeti PC, Varilly P, Fry B, et al. Genome-wide detection and characterization of positive selection in human populations. Nature 2007;449:913-8. https://doi.org/10.1038/nature06250
  8. Wang W, Li Z, Xie G, et al. Convergent genomic signatures of cashmere traits: evidence for natural and artificial selection. Int J Mol Sci 2023;24:1165. https://doi.org/10.3390/ijms24021165
  9. Jin M, Lu J, Fei X, et al. Genetic signatures of selection for cashmere traits in Chinese goats. Animals 2020;10:1905. https://doi.org/10.3390/ani10101905
  10. Diao X, Yao L, Wang X, et al. Hair follicle development and cashmere traits in Albas goat kids. Animals 2023;13:617. https://doi.org/10.3390/ani13040617
  11. Zhao B, Suo L, Wu Y, et al. Stress adaptation in Tibetan cashmere goats is governed by inherent metabolic differences and manifested through variable cashmere phenotypes. Genomics 2024;116:110801. https://doi.org/10.1016/j.ygeno.2024.110801
  12. Rong Y, Wang X, Na Q, et al. Genome-wide association study for cashmere traits in Inner Mongolia cashmere goat population reveals new candidate genes and haplotypes. BMC Genomics 2024;25:658. https://doi.org/10.1186/s12864-024-10543-4
  13. Chen S, Zhou Y, Chen Y, Gu J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018;34:i884-90. https://doi.org/10.1093/bioinformatics/bty560
  14. 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
  15. 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
  16. 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
  17. 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
  18. Zhou J, Liu L, Lopdell TJ, Garrick DJ, Shi Y. HandyCNV: standardized summary, annotation, comparison, and visualization of copy number variant, copy number variation region, and runs of homozygosity. Front Genet 2021;12:731355. https://doi.org/10.3389/fgene.2021.731355
  19. 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
  20. Szpiech ZA, Hernandez RD. selscan: An efficient multithreaded program to perform EHH-based scans for positive selection. Mol Biol Evol 2014;31:2824-7. https://doi.org/10.1093/molbev/msu211
  21. Gondro C, Porto-Neto LR, Lee SH. R for genome-wide association studies. In: Gondro C, van der Werf J, Hayes B, editors. Genome-wide association studies and genomic prediction. Methods in molecular biology. Humana Press; 2013. pp. 1-17.
  22. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 2010;26:841-2. https://doi.org/10.1093/bioinformatics/btq033
  23. Wu T, Hu E, Xu S, et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innovation 2021;2: 100141. https://doi.org/10.1016/j.xinn.2021.100141
  24. Xu T, Guo X, Wang H, et al. Differential gene expression analysis between anagen and telogen of Capra hircus skin based on the de novo assembled transcriptome sequence. Gene 2013;520:30-8. https://doi.org/10.1016/j.gene.2013.01.068
  25. Diao X, Qin J, Dong C, He L, Zhang W. Integrative transcriptomic and metabolomic analyses reveal the role of melatonin in promoting secondary hair follicle development in cashmere goats. BMC Genomics 2025;26:200. https://doi.org/10.1186/s12864-025-11389-0
  26. Ge W, Wang SH, Sun B, et al. Melatonin promotes cashmere goat (Capra hircus) secondary hair follicle growth: a view from integrated analysis of long non-coding and coding RNAs. Cell Cycle 2018;17:1255-67. https://doi.org/10.1080/15384101.2018.1471318
  27. Ma S, Cao W, Ma X, et al. Metabolomics reveals metabolites associated with hair follicle cycle in cashmere goats. BMC Vet Res 2024;20:208. https://doi.org/10.1186/s12917-024-04057-0
  28. Vasu M, Ahlawat S, Arora R, Sharma R. Deciphering the molecular drivers for cashmere/pashmina fiber production in goats: a comprehensive review. Mamm Genome 2025;36: 162-82. https://doi.org/10.1007/s00335-025-10109-z
  29. Ceballos FC, Joshi PK, Clark DW, Ramsay M, Wilson JF. Runs of homozygosity: windows into population history and trait architecture. Nat Rev Genet 2018;19:220-34. https://doi.org/10.1038/nrg.2017.109
  30. Pemberton TJ, Absher D, Feldman MW, Myers RM, Rosenberg NA, Li JZ. Genomic patterns of homozygosity in worldwide human populations. Am J Hum Genet 2012;91:275-92. https://doi.org/10.1016/j.ajhg.2012.06.014
  31. Joost S, Jacob T, Sun X, et al. Single-cell transcriptomics of traced epidermal and hair follicle stem cells reveals rapid adaptations during wound healing. Cell Rep 2018;25:585-97. E7. https://doi.org/10.1016/j.celrep.2018.09.059
  32. Snippert HJ, Haegebarth A, Kasper M, et al. Lgr6 marks stem cells in the hair follicle that generate all cell lineages of the skin. Science 2010;327:1385-9. https://doi.org/10.1126/science.1184733
  33. Glotzer DJ, Zelzer E, Olsen BR. Impaired skin and hair follicle development in Runx2 deficient mice. Dev Biol 2008;315: 459-73. https://doi.org/10.1016/j.ydbio.2008.01.005
  34. Jia Q, Zhang S, Wang D, et al. Regulatory effects of FGF9 on dermal papilla cell proliferation in small-tailed Han sheep. Genes 2023;14:1106. https://doi.org/10.3390/genes14051106
  35. Yan H, Jin M, Li Y, et al. miR-1 regulates differentiation and proliferation of goat hair follicle stem cells by targeting IGF1R and LEF1 genes. DNA Cell Biol 2022;41:190-201. https://doi.org/10.1089/dna.2021.0288
  36. Ma T, Li J, Jiang Q, Wu S, Jiang H, Zhang Q. Differential expression of miR-let7a in hair follicle cycle of Liaoning cashmere goats and identification of its targets. Funct Integr Genomics 2018;18:701-7. https://doi.org/10.1007/s10142-018-0616-x
  37. Brunner MAT, Jagannathan V, Waluk DP, et al. Novel insights into the pathways regulating the canine hair cycle and their deregulation in alopecia X. PLOS ONE 2017;12:e0186469. https://doi.org/10.1371/journal.pone.0186469
  38. Hwang D, Lee H, Lee J, et al. Micro-current stimulation has potential effects of hair growth-promotion on human hair follicle-derived papilla cells and animal model. Int J Mol Sci 2021;22:4361. https://doi.org/10.3390/ijms22094361
  39. Liang W, Chen X, Ni N, et al. Corticotropin-releasing hormone inhibits autophagy by suppressing PTEN to promote apoptosis in dermal papilla cells. Ann Med 2025;57:2490823. https://doi.org/10.1080/07853890.2025.2490823
  40. Sohn KM, Jeong KH, Kim JE, Park YM, Kang H. Hair growth-promotion effects of different alternating current parameter settings are mediated by the activation of Wnt/β-catenin and MAPK pathway. Exp Dermatol 2015;24:958-63. https://doi.org/10.1111/exd.12827