Comprehensive genome and transcriptome analyses reveal genetic relationship, selection signature, and transcriptome landscape of small-sized Korean native Jeju horse.
- Journal Article
- Research Support
- Non-U.S. Gov't
Summary
This research article primarily focuses on the analysis and study of the unique genetic properties of the Jeju horse breed, which is native to Jeju Island in Korea. The Jeju horses are believed to be originally descended from Mongolian breeds; however, they’ve adapted to their local environment and have developed distinct characteristics.
Summary of Research Work
The researchers sequenced and compared the genome of 41 horses from six different breeds. The aim was to understand and identify the distinct genetic components of the Jeju horse breed through in-depth genomic and transcriptome analysis. This was conducted for the first time considering the specific breed.
- The team discovered numerous breed-specific non-synonymous SNPs (single nucleotide polymorphisms) and loss-of-function mutations specific to the Jeju horse breed.
- Despite its close genetic proximity to Mongolian breeds, the demographic and admixture analyses suggest that the Jeju horse’s genetic ancestry is independent of Mongolian breeds. It implies that over time, the Jeju horse has developed a unique genetic lineage.
- The research also looked at genome-wide selection signature analysis, which revealed certain genes (LCORL, MSTN, HMGA2, ZFAT, LASP1, PDK4, and ACTN2) to be positively selected in the Jeju horse. This points towards a genetic bias towards these genes that might be assisting the Jeju horse with its survival and adaptation in its habitat.
RNAseq and Comparative Muscle Fiber Analysis
In addition to the genomic study, the researchers also conducted a RNAseq analysis which further demonstrated the distinct genetic makeup of the Jeju horse breed.
- This analysis illustrated that several genes which were positively selected (as identified in the genomic analysis) were also differentially expressed when compared to a Thoroughbred horse breed. This suggests that these genes could be responsible for the specific characteristics of the Jeju horse breed.
- Furthermore, a comparative muscle fiber analysis between the Jeju horse and the Thoroughbred horse revealed a significantly higher type I muscle fiber content in the former. Type I muscle fibers can influence several characteristics related to stamina and endurance, which might be another unique adaptation trait for the Jeju horse.
This comprehensive genetic analysis has offered insights into how complex phenotypic traits like small body size are selected and inherited in the Jeju horse breed. The novel SNPs identified in the research could be pivotal in designing a high-density SNP chip, which could further assist in studying other native horse breeds.
Cite This Article
Publication
Researcher Affiliations
- Animal Genomics and Bioinformatics Division, National Institute of Animal Science, Rural Development Administration, Wanju, 55365, Republic of Korea.
- Subtropical Livestock Research Institute, National Institute of Animal Science, Rural Development Administration, Jeju-do, 63242, Republic of Korea.
- Animal Genomics and Bioinformatics Division, National Institute of Animal Science, Rural Development Administration, Wanju, 55365, Republic of Korea.
- Cancer Genetics and Comparative Genomics Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, 20892, USA.
- C&K genomics, Seoul, Republic of Korea.
- Animal Breeding and Genetics Division, National Institute of Animal Science, Rural Development Administration, Wanju, 55365, Republic of Korea.
- Animal Genomics and Bioinformatics Division, National Institute of Animal Science, Rural Development Administration, Wanju, 55365, Republic of Korea.
- Animal Genomics and Bioinformatics Division, National Institute of Animal Science, Rural Development Administration, Wanju, 55365, Republic of Korea.
- College of Animal Life Science, Kangwon National University, Chuncheon, 24341, Republic of Korea.
- Animal Genomics and Bioinformatics Division, National Institute of Animal Science, Rural Development Administration, Wanju, 55365, Republic of Korea.
- C&K genomics, Seoul, Republic of Korea.
- Division of Biotechnology, Jeju National University, Jeju, 63243, Republic of Korea.
- Department of Life Science, Hanyang University, Seoul, 133-791, Republic of Korea.
- Animal Genomics and Bioinformatics Division, National Institute of Animal Science, Rural Development Administration, Wanju, 55365, Republic of Korea.
- Department of Animal Science and Technology, College of Biotechnology and Natural Resources, Chung-Ang University, Ansung-si, 17546, Republic of Korea. junmokim@cau.ac.kr.
- Animal Genomics and Bioinformatics Division, National Institute of Animal Science, Rural Development Administration, Wanju, 55365, Republic of Korea. lim.dj@korea.kr.
MeSH Terms
- Animals
- Breeding
- Gene Expression Profiling
- Genetics, Population
- Genome
- Horses / genetics
- Phenotype
- Polymorphism, Single Nucleotide
- Republic of Korea
- Selection, Genetic
- Transcriptome
Conflict of Interest Statement
References
- Jo Y-S, Kim T-W, Choi B-J, Oh H-S. Current status of terrestrial mammals on Jeju Island.. Journal of Species Research 2012;1:249–256.
- Kim B-W. Characterization of the European type of maternal lineage evident in extant Jeju native pigs.. Genes & Genomics 2011;33:111.
- Yoon SH. Complete mitochondrial genome sequences of Korean native horse from Jeju Island: uncovering the spatio-temporal dynamics.. Molecular biology reports 2017;44:233–242.
- Levine M. Investigating the origins of horse domestication.. Equine Veterinary Journal 1999;31:6–14.
- Vilà C. Widespread origins of domestic horse lineages.. Science 2001;291:474–477.
- Kim N Y. Genome-wide analyses of the Jeju, Thoroughbred, and Jeju crossbred horse populations using the high density SNP array.. Genes & genomics 2018;1–10.
- Do K-T, Lee J-H, Lee H-K, Kim J, Park K-D. Estimation of effective population size using single-nucleotide polymorphism (SNP) data in Jeju horse.. Journal of animal science and technology 2014;56:28.
- Lee J-H, Song K-D, Kim J-M, Leem H-K, Park K-D. Identification of genes with nonsynonymous SNP in Jeju horse by whole-genome resequencing reveals a functional role for immune response.. Journal of animal science 2016;94:895–901.
- Nam D. Horse production in Cheju during Lee dynasty.. Studies on Korean History 1969;4:131–131.
- Shin T. An anatomy study of animal bones excavated in the Kwakji archaeological site in Cheju Island.. Go-Moon-Wha 1992;40:31–42.
- Kim KI. Phylogenetic relationships of Cheju horses to other horse breeds as determined by mtDNA D‐loop sequence polymorphism.. Animal Genetics 1999;30:102–108.
- Oh M. Phylogenetic relationship of Cheju native horses by mitochondrial DNA analysis.. Molecules and Cells(Korea Republic) 1994.
- Xu S. High altitude adaptation and phylogenetic analysis of Tibetan horse based on the mitochondrial genome.. Journal of Genetics and Genomics 2007;34:720–729.
- Jung Y-H, Han S-H, Shin T, Oh M-Y. Genetic characterization of horse bone excavated from the Kwakji archaeological site, Jeju, Korea.. Molecules and cells 2002;14:224–230.
- Yang Y, Kim K, Cothran E, Flannery A. Genetic diversity of Cheju horses (Equus caballus) determined by using mitochondrial DNA D-loop polymorphism.. Biochemical genetics 2002;40:175–186.
- Kim NY. Estimation of genetic parameters for temperament in Jeju crossbred horses.. Asian-Australasian journal of animal sciences 2018;31:1098.
- Brown-Douglas C G, Pagan J D. Body weight, wither height and growth rates in Thoroughbreds raised in America, England, Australia, New Zealand and India.. Advances in EquineNutrition IV 2009;213.
- Metzger J, Schrimpf R, Philipp U, Distl O. Expression levels of LCORL are associated with body size in horses.. PLoS One 2013;8:e56497.
- Lee J-R. Genome-wide analysis of DNA methylation patterns in horse.. BMC genomics 2014;15:598.
- Lee W. Analysis of cross-population differentiation between Thoroughbred and Jeju horses.. Asian-Australasian journal of animal sciences 2018;31:1110.
- Choi S-K, Cho C-Y, Yeon S-H, Cho B-W, Cho G-J. Genetic characterization and polymorphisms for parentage testing of the Jeju horse using 20 microsatellite loci.. Journal of Veterinary Medical Science 2008;70:1111–1115.
- Gim J-A. HEpD: A database describing epigenetic differences between Thoroughbred and Jeju horses.. Gene 2015;560:83–88.
- Schubert M. Prehistoric genomes reveal the genetic foundation and cost of horse domestication.. Proceedings of the National Academy of Sciences 2014;111:E5661–E5669.
- Do K-T. Genomic characterization of the Przewalski׳ s horse inhabiting Mongolian steppe by whole genome re-sequencing.. Livestock Science 2014;167:86–91.
- Goto H. A massively parallel sequencing approach uncovers ancient origins and high genetic variability of endangered Przewalski’s horses.. Genome biology and evolution 2011;3:1096–1106.
- Huang J. Analysis of horse genomes provides insight into the diversification and adaptive evolution of karyotype.. Scientific reports 2014;4:4958.
- Volf J, Kus E, Prokopova L. General studbook of the Przewalski horse.. Zoological Garden Prague, Prague (1991).
- Minjigrorj N, Austbo D. Production of mare’s milk in Mongolia.. Annals of Microbiology 2009;59:493–498.
- Zhang C. Detecting the Population Structure and Scanning for Signatures of Selection in Horses (Equus caballus) From Whole-Genome Sequencing Data.. Evolutionary Bioinformatics 2018;14:1176934318775106.
- Gouveia JJdS, Silva MVGBd, Paiva SR, Oliveira SMPd. Identification of selection signatures in livestock species.. Genetics and molecular biology 2014;37:330–342.
- Avise JC, Ayala FJ. In the light of evolution III: Two centuries of Darwin.. Proceedings of the National Academy of Sciences 2009;106:9933–9938.
- Driscoll CA, Macdonald DW, O’Brien SJ. From wild animals to domestic pets, an evolutionary view of domestication.. Proceedings of the National Academy of Sciences 2009;106:9971–9978.
- Gregory TR. Artificial selection and domestication: modern lessons from Darwin’s enduring analogy.. Evolution: Education and Outreach 2009;2:5.
- Wade C. Genome sequence, comparative analysis, and population genetics of the domestic horse.. Science 2009;326:865–867.
- Kim J. The genome landscape of indigenous African cattle.. Genome biology 2017;18:34.
- Moon S. A genome-wide scan for selective sweeps in racing horses.. Asian-Australasian journal of animal sciences 2015;28:1525.
- Choi J-W. Whole-genome analyses of Korean native and Holstein cattle breeds by massively parallel sequencing.. PloS one 2014;9:e101127.
- Consortium GP. An integrated map of genetic variation from 1,092 human genomes.. Nature 2012;491:56.
- Choi J-W. Whole-genome resequencing analyses of five pig breeds, including Korean wild and native, and three European origin breeds.. DNA Research 2015;22:259–267.
- Wang J, Raskin L, Samuels DC, Shyr Y, Guo Y. Genome measures used for quality control are dependent on gene function and ancestry.. Bioinformatics 2014;31:318–323.
- Cingolani P. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3.. Fly 2012;6:80–92.
- Kaiser VB. Homozygous loss-of-function variants in European cosmopolitan and isolate populations.. Human molecular genetics 2015;24:5464–5474.
- Oh HJ, Choi D, Goh CJ, Hahn Y. Loss of gene function and evolution of human phenotypes.. BMB reports 2015;48:373.
- Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals.. Genome research 2009;19:1655–1664.
- Boitard S, Rodriguez W, Jay F, Mona S, Austerlitz F. Inferring population size history from large samples of genome-wide molecular data-an approximate Bayesian computation approach.. PLoS genetics 2016;12:e1005877.
- Pickrell JK, Pritchard JK. Inference of population splits and mixtures from genome-wide allele frequency data.. PLoS genetics 2012;8:e1002967.
- Nozawa K, Kondo K. Gene constitution of Cheju native horse and its phylogenetic relationships with Japanese native horses.. SABRAO Newsletter 1970;2:7–18.
- Hill E W. Correction: A sequence polymorphism in MSTN predicts sprinting ability and racing stamina in thoroughbred horses.. PloS one 2010;5.
- Gu J. Association of sequence variants in CKM (creatine kinase, muscle) and COX4I2 (cytochrome c oxidase, subunit 4, isoform 2) genes with racing performance in Thoroughbred horses.. Equine Veterinary Journal 2010;42:569–575.
- Hill E, Gu J, McGivney B, MacHugh D. Targets of selection in the Thoroughbred genome contain exercise‐relevant gene SNPs associated with elite racecourse performance.. Animal genetics 2010;41:56–63.
- Sabeti PC. Genome-wide detection and characterization of positive selection in human populations.. Nature 2007;449:913.
- Rubin C-J. Strong signatures of selection in the domestic pig genome.. Proceedings of the National Academy of Sciences 2012;109:19529–19536.
- Qanbari S. A high resolution genome-wide scan for significant selective sweeps: an application to pooled sequence data in laying chickens.. PloS one 2012;7:e49525.
- Qanbari S, Simianer H. Mapping signatures of positive selection in the genome of livestock.. Livestock Science 2014;166:133–143.
- Signer-Hasler H. A genome-wide association study reveals loci influencing height and other conformation traits in horses.. PloS one 2012;7:e37282.
- Makvandi-Nejad S. Four loci explain 83% of size variation in the horse.. PLoS One 2012;7:e39929.
- Schröder W. Athletic performance and conformation in Hanoverian warmblood horses-population genetic and genome-wide association analyses. cumulative thesis.. Hannover: University of Veterinary Medicine (2010).
- Junior AB. Polymorphisms in the LASP1 gene allow selection for smaller stature in ponies.. Livestock Science 2018;216:160–164.
- Kang M-s. Grazing behaviour of Jeju Native Horses.. Journal of integrated field science 2004;1:33–34.
- Tsunoda T. Immune-related zinc finger gene ZFAT is an essential transcriptional regulator for hematopoietic differentiation in blood islands.. Proceedings of the National Academy of Sciences 2010;107:14199–14204.
- Allen HL. Hundreds of variants clustered in genomic loci and biological pathways affect human height.. Nature 2010;467:832.
- Takeuchi F. Evaluation of genetic loci influencing adult height in the Japanese population.. Journal of human genetics 2009;54:749.
- N’Diaye A. Identification, replication, and fine-mapping of Loci associated with adult height in individuals of african ancestry.. PLoS genetics 2011;7:e1002298.
- Cleynen I, Van de Ven WJ. The HMGA proteins: a myriad of functions.. International journal of oncology 2008;32:289–305.
- Weedon MN. A common variant of HMGA2 is associated with adult and childhood height in the general population.. Nature genetics 2007;39:1245.
- Boyko AR. A simple genetic architecture underlies morphological variation in dogs.. PLoS biology 2010;8:e1000451.
- Jones P. Single-nucleotide-polymorphism-based association mapping of dog stereotypes.. Genetics 2008;179:1033–1044.
- Davis S. Molecular mechanisms of pituitary organogenesis: in search of novel regulatory genes.. Molecular and cellular endocrinology 2010;323:4–19.
- Deladoëy J. “Hot spot” in the PROP1 gene responsible for combined pituitary hormone deficiency.. The Journal of Clinical Endocrinology & Metabolism 1999;84:1645-–1650.
- Sornson MW. Pituitary lineage determination by the Prophet of Pit-1 homeodomain factor defective in Ames dwarfism.. Nature 1996;384:327.
- Orr N. Genome‐wide SNP association–based localization of a dwarfism gene in Friesian dwarf horses.. Animal Genetics 2010;41:2–7.
- Rivero J, Serrano AL, Henckel P, Aguera E. Muscle fiber type composition and fiber size in successfully and unsuccessfully endurance-raced horses.. Journal of Applied Physiology 1993;75:1758–1766.
- Dall’Olio S. Analysis of horse myostatin gene and identification of single nucleotide polymorphisms in breeds of different morphological types.. BioMed Research International 2010.
- Thomas K, Hamilton N, North K, Houweling P. Sequence analysis of the equine ACTN3 gene in Australian horse breeds.. Gene 2014;538:88–93.
- Wang J. Analysis of Equine ACTN3 Gene Polymorphisms in Yili Horses.. Journal of Equine Veterinary Science 2018;70:101–106.
- McPherron AC, Lawler AM, Lee S-J. Regulation of skeletal muscle mass in mice by a new TGF-p superfamily member.. Nature 1997;387:83.
- McPherron AC, Lee S-J. Double muscling in cattle due to mutations in the myostatin gene.. Proceedings of the National Academy of Sciences 1997;94:12457–12461.
- Miyata H. Effect of Myostatin SNP on muscle fiber properties in male Thoroughbred horses during training period.. The Journal of Physiological Sciences 2018;68:639–646.
- Rooney MF, Porter RK, Katz LM, Hill EW. Skeletal muscle mitochondrial bioenergetics and associations with myostatin genotypes in the Thoroughbred horse.. PloS one 2017;12:e0186247.
- Petersen JL, Valberg SJ, Mickelson JR, McCue ME. Haplotype diversity in the equine myostatin gene with focus on variants associated with race distance propensity and muscle fiber type proportions.. Animal genetics 2014;45:827–835.
- Petersen JL. Genome-wide analysis reveals selection for important traits in domestic horse breeds.. PLoS genetics 2013;9:e1003211.
- McGivney BA. MSTN genotypes in T horoughbred horses influence skeletal muscle gene expression and racetrack performance.. Animal genetics 2012;43:810–812.
- Constantinopol M. Oxygen transport during exercise in large mammals. II. Oxygen uptake by the pulmonary gas exchanger.. Journal of Applied Physiology 1989;67:871–878.
- Erickson B. Mechanism of reduction in alveolar-arterial PO2 difference by helium breathing in the exercising horse.. Journal of Applied Physiology 1994;76:2794–2801.
- Flück M. Functional, structural and molecular plasticity of mammalian skeletal muscle in response to exercise stimuli.. Journal of Experimental Biology 2006;209:2239–2248.
- Eivers SS. Alterations in oxidative gene expression in equine skeletal muscle following exercise and training.. Physiological genomics 2009;40:83–93.
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2.. Nature methods 2012;9:357.
- Li H. The Sequence Alignment/Map format and SAMtools.. Bioinformatics 2009;25:2078–2079.
- McKenna A. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data.. Genome research 2010;20:1297–1303.
- Aslam ML. Whole genome SNP discovery and analysis of genetic diversity in Turkey (Meleagris gallopavo). BMC genomics 2012;13:391.
- Kalbfleisch TS. Eq쪳, an Updated Reference Genome for the Domestic Horse.. BioRxiv 2018;306928.
- Danecek P. The variant call format and VCFtools.. Bioinformatics 2011;27:2156–2158.
- Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: a tool for genome-wide complex trait analysis.. The American Journal of Human Genetics 2011;88:76–82.
- Rodriguez W, Jay F, Mona S, Austerlitz F. Inferring population size history from large samples of genome-wide molecular data-an approximate bayesian computation approach.. Plos Genetics 2016;3(12):1–36.
- Purcell S. PLINK: a tool set for whole-genome association and population-based linkage analyses.. The American Journal of Human Genetics 2007;81:559–575.
- Rubin C-J. Whole-genome resequencing reveals loci under selection during chicken domestication.. Nature 2010;464:587.
- Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data.. Bioinformatics 2014;30:2114–2120.
- Andrews S. FastQC: a quality control tool for high throughput sequence data.. 2010.
- Kim, D., Langmead, B. & Salzberg, S. (2016).
- Liao Y, Smyth GK, Shi W. The Subread aligner: fast, accurate and scalable read mapping by seed-and-vote.. Nucleic acids research 2013;41:e108–e108.
- Robinson MD, McCarthy DJ, Smyth G K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data.. Bioinformatics 2010;26:139–140.
- Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources.. Nature protocols 2008;4:44.
- Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms.. PloS one 2011;6:e21800.
- Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters.. Omics: a journal of integrative biology 2012;16:284–287.
- Lee SH, Kim J-M, Ryu YC, Ko KS. Effects of Morphological Characteristics of Muscle Fibers on Porcine Growth Performance and Pork Quality.. Korean journal for food science of animal resources 2016;36:583.
- Brooke MH, Kaiser KK. Muscle fiber types: how many and what kind?. Archives of neurology 1970;23:369–379.
- Hansen M, Knorr C, Hall A, Broad T, Brenig B. Sequence analysis of the equine SLC26A2 gene locus on chromosome 14q15→ q21.. Cytogenetic and genome research 2007;118:55–62.
Citations
This article has been cited 13 times.- Jafari H, Abebe BK, Cong L, Ahmed Z, Zhaofei W, Sun M, Muhatai G, Chuzhao L, Dang R. Review: Genomic insights into the adaptive traits and stress resistance in modern horses. Stress Biol 2026 Jan 12;6(1):5.
- Burke D, Raychaudhuri J, Chuong E, Taylor W, Layer R. TEPEAK: A novel method for identifying and characterizing polymorphic transposable elements in non-model species populations. PLoS Comput Biol 2026 Jan;22(1):e1013122.
- Austin MMP, Ivey JLZ, Shepherd EA, Myer PR. Methodologies to Identify Metabolic Pathway Differences Between Emaciated and Moderately Conditioned Horses: A Review of Multiple Gene Expression Techniques. Animals (Basel) 2025 Oct 10;15(20).
- Šimon M, Kaić A, Potočnik K. Unveiling Genetic Potential for Equine Meat Production: A Bioinformatics Approach. Animals (Basel) 2024 Aug 22;14(16).
- Reißmann M, Rajavel A, Kokov ZA, Schmitt AO. Identification of Differentially Expressed Genes after Endurance Runs in Karbadian Horses to Determine Candidates for Stress Indicators and Performance Capability. Genes (Basel) 2023 Oct 24;14(11).
- Batcher K, Varney S, Raudsepp T, Jevit M, Dickinson P, Jagannathan V, Leeb T, Bannasch D. Ancient segmentally duplicated LCORL retrocopies in equids. PLoS One 2023;18(6):e0286861.
- Srikanth K, von Pfeil DJF, Stanley BJ, Griffitts C, Huson HJ. Genome Wide Association Study with Imputed Whole Genome Sequence Data Identifies a 431 kb Risk Haplotype on CFA18 for Congenital Laryngeal Paralysis in Alaskan Sled Dogs. Genes (Basel) 2022 Oct 6;13(10).
- Wang X, Ran X, Niu X, Huang S, Li S, Wang J. Whole-genome sequence analysis reveals selection signatures for important economic traits in Xiang pigs. Sci Rep 2022 Jul 12;12(1):11823.
- Hou H, Wang X, Ding W, Xiao C, Cai X, Lv W, Tu Y, Zhao W, Yao J, Yang C. Whole-genome sequencing reveals the artificial selection and local environmental adaptability of pigeons (Columba livia). Evol Appl 2022 Apr;15(4):603-617.
- de Faria DA, do Prado Paim T, Dos Santos CA, Paiva SR, Nogueira MB, McManus C. Selection signatures for heat tolerance in Brazilian horse breeds. Mol Genet Genomics 2022 Mar;297(2):449-462.
- Hou H, Wang X, Yang C, Cai X, Lv W, Tu Y, Bao A, Wu Q, Zhao W, Yao J, Ding W. Comparative Genome and Transcriptome Integration Studies Reveal the Mechanism of Pectoral Muscle Development and Function in Pigeons. Front Genet 2021;12:735795.
- Santos WB, Schettini GP, Maiorano AM, Bussiman FO, Balieiro JCC, Ferraz GC, Pereira GL, Baldassini WA, Neto ORM, Oliveira HN, Curi RA. Genome-wide scans for signatures of selection in Mangalarga Marchador horses using high-throughput SNP genotyping. BMC Genomics 2021 Oct 14;22(1):737.
- Pu Y, Zhang Y, Zhang T, Han J, Ma Y, Liu X. Identification of Novel lncRNAs Differentially Expressed in Placentas of Chinese Ningqiang Pony and Yili Horse Breeds. Animals (Basel) 2020 Jan 11;10(1).