FREE SHIPPING over $40
OVER 9000 FIVE-STAR REVIEWS
5% OFF ALL SUBSCRIPTIONS
100% HAPPINESS GUARANTEE
Analyze Diet
0
Home / Equine Research Bank / Front Genet / Comparative analysis of blood whole transcriptome profiles i...
Frontiers in genetics2025; 16; 1651628; doi: 10.3389/fgene.2025.1651628

Comparative analysis of blood whole transcriptome profiles in Yili horses pre- and post-5000-meter racing.

Abstract: This study employed Yili horses participating in a 5000-meter race as a model to investigate exercise-induced gene expression alterations in peripheral blood using whole transcriptome sequencing. Jugular vein blood samples from the three leading horses were collected pre- and immediately post-race, yielding 2,171 differentially expressed mRNAs (2,080 upregulated, 91 downregulated), 4,375 differentially expressed LncRNAs (4,354 upregulated), and 68 differentially expressed circRNAs (64 upregulated). GO/KEGG analyses demonstrated significant enrichment of differential mRNAs in transmembrane transport function and pivotal signaling pathways (cAMP, MAPK, PI3K-Akt). Differential lncRNAs targeted neuro-signaling pathways (e.g., Neuroactive ligand-receptor interaction, Calcium signaling) and developmental regulators (stem cell pluripotency). Source genes of circRNAs were enriched in axon guidance and immune-related T cell receptor signaling. Molecular functions converged on transporter/receptor activity (mRNA/lncRNA) and nucleic acid/GTP binding (circRNA source genes). The protein-protein interaction analysis identified ten central genes within the heat shock protein family, such as and . Notably, significant upregulation of and indicated their potential roles in modulating cardiac rhythm, promoting tissue repair, and maintaining calcium-phosphorus homeostasis during exercise adaptation. This study provides comprehensive overview of transcriptomic regulatory mechanisms in the blood of Yili horses, offering a molecular framework for advancing understanding of physiological adaptation to exercise and optimizing equine exercise protocols.
Copyright © 2025 Su, Ren, Ma, Meng, Yao, Zeng, Li, Li, Wang and Wang.
Get Full Text
Publication Date: 2025-08-29 PubMed ID: 40949872PubMed Central: PMC12425706DOI: 10.3389/fgene.2025.1651628Google Scholar: Lookup
The Equine Research Bank provides access to a large database of publicly available scientific literature. Inclusion in the Research Bank does not imply endorsement of study methods or findings by Mad Barn.
LinkedInCopy
  • Journal Article

Summary

This research summary has been generated with artificial intelligence and may contain errors and omissions. Refer to the original study to confirm details provided. Submit correction.

Overview

  • This study examined how a 5000-meter race affects gene expression in the blood of Yili horses by comparing blood transcriptome profiles before and after racing.
  • The research aimed to identify key changes in RNA molecules and biological pathways involved in exercise-induced physiological adaptations.

Study Design and Methods

  • The subjects were Yili horses participating in a 5000-meter race, a model to study exercise effects on gene expression.
  • Jugular vein blood samples were collected from the top three performing horses immediately before and after the race.
  • Whole transcriptome sequencing was conducted to profile:
    • mRNAs (messenger RNAs)
    • lncRNAs (long non-coding RNAs)
    • circRNAs (circular RNAs)

Differential Expression Findings

  • mRNAs:
    • Total of 2,171 differentially expressed mRNAs detected post-race compared to pre-race.
    • Most were upregulated (2,080), indicating increased gene activity following exercise.
    • A smaller subset (91) were downregulated.
  • lncRNAs:
    • 4,375 differentially expressed, with 4,354 upregulated post-race.
    • These non-coding RNAs likely play regulatory roles in exercise response.
  • circRNAs:
    • 68 differential circRNAs identified, 64 upregulated post-exercise.
    • Circular RNAs often modulate gene expression and stability.

Functional Enrichment Analyses (GO and KEGG)

  • Analysis of differentially expressed mRNAs showed:
    • Significant enrichment in transmembrane transport functions, implying changes in molecule movement across cell membranes during exercise.
    • Activation of key signaling pathways:
      • cAMP signaling pathway – involved in energy balance and cellular responses.
      • MAPK pathway – regulates cell growth and stress responses.
      • PI3K-Akt pathway – important in metabolism and cell survival.
  • Differentially expressed lncRNAs targeted pathways:
    • Neuroactive ligand-receptor interaction – indicating modulation of nervous system signaling during exercise.
    • Calcium signaling pathway – crucial for muscle contraction and cellular signaling.
    • Stem cell pluripotency regulation – suggesting roles in cellular regeneration or repair mechanisms post-exercise.
  • Source genes of circRNAs were enriched in:
    • Axon guidance pathway – relating to neural development and function.
    • T cell receptor signaling pathway – indicating immune system involvement in exercise adaptation.
  • Molecular function convergences included:
    • Transporter and receptor activities for mRNAs and lncRNAs, suggesting modifications in cellular communication and nutrient transport.
    • Nucleic acid and GTP binding functions for circRNA source genes, implying regulation of genetic material and signaling molecules.

Protein-Protein Interaction and Key Genes

  • Ten central hub genes in the protein interaction network belonged to the heat shock protein (HSP) family, known for their role in protecting cells from stress.
  • Two genes, notably upregulated post-race and linked to:
    • Cardiac rhythm regulation – important for maintaining heart function during and after exertion.
    • Tissue repair processes – facilitating recovery from exercise-induced damage.
    • Calcium and phosphorus homeostasis – critical in muscle contraction and bone maintenance during physical stress.

Significance and Applications

  • This research provides a comprehensive molecular overview of how Yili horses adapt physiologically to exercise at the transcriptome level in blood cells.
  • The identification of key RNA molecules and pathways offers insights into:
    • The regulatory mechanisms that underpin exercise response and adaptation in horses.
    • Potential biomarkers for monitoring training effects and health status.
    • Guiding the development of optimized training regimens tailored to molecular responses.

Cite This Article

APA
Su Y, Ren W, Ma S, Meng J, Yao X, Zeng Y, Li Z, Li L, Wang R, Wang J. (2025). Comparative analysis of blood whole transcriptome profiles in Yili horses pre- and post-5000-meter racing. Front Genet, 16, 1651628. https://doi.org/10.3389/fgene.2025.1651628

Publication

Frontiers in genetics
ISSN: 1664-8021
NlmUniqueID: 101560621
Country: Switzerland
Language: English
Volume: 16
Pages: 1651628
PII: 1651628

Researcher Affiliations

Su, Yi
  • College of Animal Science, Xinjiang Agricultural University, Xinjiang, China.
Ren, Wanlu
  • College of Animal Science, Xinjiang Agricultural University, Xinjiang, China.
  • Xinjiang Key Laboratory of Equine Breeding and Exercise Physiology, Urumqi, China.
Ma, Shikun
  • College of Animal Science, Xinjiang Agricultural University, Xinjiang, China.
Meng, Jun
  • College of Animal Science, Xinjiang Agricultural University, Xinjiang, China.
  • Xinjiang Key Laboratory of Equine Breeding and Exercise Physiology, Urumqi, China.
Yao, Xinkui
  • College of Animal Science, Xinjiang Agricultural University, Xinjiang, China.
  • Xinjiang Key Laboratory of Equine Breeding and Exercise Physiology, Urumqi, China.
Zeng, Yaqi
  • College of Animal Science, Xinjiang Agricultural University, Xinjiang, China.
  • Xinjiang Key Laboratory of Equine Breeding and Exercise Physiology, Urumqi, China.
Li, Zexu
  • College of Animal Science, Xinjiang Agricultural University, Xinjiang, China.
Li, Luling
  • College of Animal Science, Xinjiang Agricultural University, Xinjiang, China.
Wang, Ran
  • College of Animal Science, Xinjiang Agricultural University, Xinjiang, China.
Wang, Jianwen
  • College of Animal Science, Xinjiang Agricultural University, Xinjiang, China.
  • Xinjiang Key Laboratory of Equine Breeding and Exercise Physiology, Urumqi, China.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

This article includes 49 references
  1. Akizuki K, Ono A, Xue H C, Kameshita I, Ishida A, Sueyoshi N. Biochemical characterization of four splice variants of mouse Ca2+/calmodulin-dependent protein kinase Iδ. J. Biochem. 169, 445–458.
    doi: 10.1093/jb/mvaa117pubmed: 33417706google scholar: lookup
  2. Ana C, Pedro M, Sonia T, Gomez-Cabrero D, Cervera A, McPherson A. A survey of best practices for RNA-Seq data analysis. Genome Biol. 17 (1), 13.
    doi: 10.1186/s13059-016-0881-8pmc: PMC4728800pubmed: 26813401google scholar: lookup
  3. Bouwman W, Verhaegh W, van Doorn A, Raymakers R, van der Poll T, van de Stolpe A. Quantitative characterization of immune cells by measuring cellular signal transduction pathway activity. Sci. Rep-UK 14, 24487.
    doi: 10.1038/s41598-024-75666-wpmc: PMC11489675pubmed: 39424625google scholar: lookup
  4. Bryan K, McGivney B A, Farries G, McGettigan P A, McGivney C L, Gough K F. Equine skeletal muscle adaptations to exercise and training: evidence of differential regulation of autophagosomal and mitochondrial components. BMC Genomics 18, 595.
    doi: 10.1186/s12864-017-4007-9pmc: PMC5551008pubmed: 28793853google scholar: lookup
  5. Cappelli K, Capomaccio S, Viglino A, Silvestrelli M, Beccati F, Moscati L. Circulating miRNAs as putative biomarkers of exercise adaptation in endurance horses. Front. Physiol. 9, 429.
    doi: 10.3389/fphys.2018.00429pmc: PMC5928201pubmed: 29740341google scholar: lookup
  6. Chan R J, Walker A, Vardy J, Chan A, Oppegaard K, Conley Y P. Perturbations in the neuroactive ligand-receptor interaction and renin angiotensin system pathways are associated with cancer-related cognitive impairment. Support Care Cancer 33, 254.
    doi: 10.1007/s00520-025-09317-9pmc: PMC11885406pubmed: 40047999google scholar: lookup
  7. Christopoulos P F, Msaouel P, Koutsilieris M. The role of the insulin-like growth factor-1 system in breast cancer. Mol. Cancer 14, 43.
    doi: 10.1186/s12943-015-0291-7pmc: PMC4335664pubmed: 25743390google scholar: lookup
  8. de Almeida R M C, Clendenon S G, Richards W G, Boedigheimer M, Damore M, Rossetti S. Transcriptome analysis reveals manifold mechanisms of cyst development in ADPKD. Hum. Genomics 10, 37.
    doi: 10.1186/s40246-016-0095-xpmc: PMC5117508pubmed: 27871310google scholar: lookup
  9. Farries G, Bryan K, McGivney C L, McGettigan P A, Gough K F, Browne J A. Expression quantitative trait loci in equine skeletal muscle reveals heritable variation in metabolism and the training responsive transcriptome. Front. Genet. 10, 1215.
    doi: 10.3389/fgene.2019.01215pmc: PMC6902038pubmed: 31850069google scholar: lookup
  10. Flick M, Vinther A M L, Jacobsen S, Berg L C, Gimeno M, Verwilghen D. Effect of exercise on serum neutrophil gelatinase-associated lipocalin concentration in racehorses. Vet. Clin. Path 50, 551–554.
    doi: 10.1111/vcp.13027pubmed: 34779025google scholar: lookup
  11. Fu X K, Zhou B, Yan Q N, Tao C, Qin L, Wu X H. Kindlin-2 regulates skeletal homeostasis by modulating PTH1R in mice. Signal Transduct. Tar 5, 297.
    doi: 10.1038/s41392-020-00328-ypmc: PMC7762753pubmed: 33361757google scholar: lookup
  12. Gesty-Palmer D, Chen M, Reiter E, Ahn S, Nelson C D, Wang S. Distinct beta-arrestin- and G protein-dependent pathways for parathyroid hormone receptor-stimulated ERK1/2 activation. J. Biol. Chem. 281, 10856–10864.
    doi: 10.1074/jbc.M513380200pubmed: 16492667google scholar: lookup
  13. Giers J, Bartel A, Kirsch K, Mueller S F, Horstmann S, Gehlen H. Blood-based assessment of oxidative stress, inflammation, endocrine and metabolic adaptations in eventing horses accounting for plasma volume shift after exercise. Vet. Med. Sci. 10, e1409.
    doi: 10.1002/vms3.1409pmc: PMC10958401pubmed: 38516822google scholar: lookup
  14. Glebov-McCloud A G P, Saide W S, Gaine M E, Strack S. Protein kinase A in neurological disorders. J. Neurodev. Disord. 16, 9.
    doi: 10.1186/s11689-024-09525-0pmc: PMC10936040pubmed: 38481146google scholar: lookup
  15. Goltzman D, Hendy G N. The calcium-sensing receptor in bone-mechanistic and therapeutic insights. Nat. Rev. Endocrinol. 11, 298–307.
    doi: 10.1038/nrendo.2015.30pubmed: 25752283google scholar: lookup
  16. Goyal A, Agrawal A, Verma A, Dubey N. The PI3K-AKT pathway: a plausible therapeutic target in parkinson's disease. Exp. Mol. Pathol. 129, 104846.
    doi: 10.1016/j.yexmp.2022.104846pubmed: 36436571google scholar: lookup
  17. Harari S, Deretz S, Saint Priest BD, Richard E, Ricard A. Comparison of blood parameters in two genetically different groups of horses for functional longevity in show jumping. Front. Genet. 15, 1455790.
    doi: 10.3389/fgene.2024.1455790pmc: PMC11554460pubmed: 39534078google scholar: lookup
  18. Heckman PRA, Blokland A, Bollen EPP, Prickaerts J. Phosphodiesterase inhibition and modulation of corticostriatal and hippocampal circuits: clinical overview and translational considerations. Neurosci. Biobehav R. 87, 233–254.
    doi: 10.1016/j.neubiorev.2018.02.007pubmed: 29454746google scholar: lookup
  19. Hondares E, Gallego-Escuredo JM, Flachs P, Frontini A, Cereijo R, Goday A. Fibroblast growth factor-21 is expressed in neonatal and pheochromocytoma-induced adult human brown adipose tissue. Metabolism 63, 312–317.
    doi: 10.1016/j.metabol.2013.11.014pubmed: 24369918google scholar: lookup
  20. Huangsaksri O, Sanigavatee K, Poochipakorn C, Wonghanchao T, Yalong M, Thongcham K. Physiological stress responses in horses participating in novice endurance rides. Heliyon 10, e31874.
    doi: 10.1016/j.heliyon.2024.e31874pmc: PMC11153225pubmed: 38845925google scholar: lookup
  21. Hutchinson JA. Quantification of mRNA expression by RT-qPCR. Transplantation 99 (10), 2009–2011.
    doi: 10.1097/TP.0000000000000948pubmed: 26426919google scholar: lookup
  22. Javier AJS, Kennedy FM, Yi X, Wehling-Henricks M, Tidball JG, White KE. Klotho is cardioprotective in the mdx mouse model of Duchenne muscular dystrophy. Am. J. Pathol. 29, 923–940.
    doi: 10.1016/j.ajpath.2024.12.017pmc: PMC12016860pubmed: 39889824google scholar: lookup
  23. Jia W, Dan H, Lina H, Qinzhen Z, Shaozhuang L, Jiapeng Hu. Protective effect of micheliolide against inflammation and oxidative stress in asthma through the MAPK/NF-κB signalling pathway. World Allergy Organ. J. 18 (8), 101091.
    doi: 10.1016/J.WAOJOU.2025.101091pmc: PMC12296505pubmed: 40718157google scholar: lookup
  24. Jin R, Pei HL, Yue F, Zhang XD, Zhang ZC, Xu Y. Network pharmacology combined with metabolomics reveals the mechanism of yangxuerongjin pill against type 2 diabetic peripheral neuropathy in rats. Drug Des. Dev. Ther. 19, 325–347.
    doi: 10.2147/DDDT.S473146pmc: PMC11745066pubmed: 39834645google scholar: lookup
  25. Lee J, Hong SW, Park SE, Rhee EJ, Park CY, Oh KW. Exendin-4 regulates lipid metabolism and fibroblast growth factor 21 in hepatic steatosis. Metabolism 63, 1041–1048.
    doi: 10.1016/j.metabol.2014.04.011pubmed: 24933399google scholar: lookup
  26. Li DJ, Fu H, Zhao T, Ni M, Shen FM. Exercise-stimulated FGF23 promotes exercise performance controlling the excess reactive oxygen species production and enhancing mitochondrial function in skeletal muscle. Metabolism 65, 747–756.
    doi: 10.1016/j.metabol.2016.02.009pubmed: 27085781google scholar: lookup
  27. Li PP, Wang JJ, Zou Y, Sun ZL, Zhang MH, Geng ZM. Interaction of Hsp90AA1 with phospholipids stabilizes membranes under stress conditions. BBA-Biomembranes 1861, 457–465.
    doi: 10.1016/j.bbamem.2018.11.009pubmed: 30517848google scholar: lookup
  28. Li JT, Gao H, Wang P, Sun C, Wei ZL, Yi XC. Plumbagin induces G2/M arrest and apoptosis and ferroptosis ROS/p38 MAPK pathway in human osteosarcoma cells. Alex Eng. J. 103, 222–236.
    doi: 10.1016/j.aej.2024.06.015google scholar: lookup
  29. Lin ZY, Lin BW, Hang CW, Lu RH, Xiong H, Liu JY. A new paradigm for generating high-quality cardiac pacemaker cells from mouse pluripotent stem cells. Signal Transduct. Tar 9, 230.
    doi: 10.1038/s41392-024-01942-wpmc: PMC11377569pubmed: 39237509google scholar: lookup
  30. Liu F, Fu P, Fan WX, Gou R, Huang YQ, Qiu HY. Involvement of parathyroid hormone-related protein in vascular calcification of chronic haemodialysis patients. Nephrology 17, 552–560.
    doi: 10.1111/j.1440-1797.2012.01601.xpubmed: 22448974google scholar: lookup
  31. Liu J, Xu Y, Hu YJ, Wang G. The role of fibroblast growth factor 21 in the pathogenesis of non-alcoholic fatty liver disease and implications for therapy. Metabolism 64, 380–390.
    doi: 10.1016/j.metabol.2014.11.009pubmed: 25516477google scholar: lookup
  32. Masko M, Domino M, Jasinski T, Witkowska-Pilaszewicz O. The physical activity-dependent hematological and biochemical changes in school horses in comparison to blood profiles in endurance and race horses. Animals 11, 1128.
    doi: 10.3390/ani11041128pmc: PMC8071065pubmed: 33920044google scholar: lookup
  33. Menetrey J, Kasemkijwattana C, Day CS, Bosch P, Vogt M, Fu FH. Growth factors improve muscle healing. J. Bone Jt. Surg. Br. 82, 131–137.
    doi: 10.1302/0301-620x.82b1.8954pubmed: 10697329google scholar: lookup
  34. Mir S, Cai WK, Carlson SW, Saatman KE, Andres DA. IGF-1 mediated neurogenesis involves a novel RIT1/Akt/Sox2 Cascade. Sci. Rep-UK 7, 3283.
    doi: 10.1038/s41598-017-03641-9pmc: PMC5468318pubmed: 28607354google scholar: lookup
  35. Neele AE, Gijbels MJJ, van der Velden S, Hoeksema MA, Boshuizen MCS, Prange KHM. Myeloid Kdm6b deficiency results in advanced atherosclerosis. Atherosclerosis 275, 156–165.
    doi: 10.1016/j.atherosclerosis.2018.05.052pubmed: 29908485google scholar: lookup
  36. Porro A, Thiel G, Moroni A, Saponaro A. Cyclic AMP regulation and its command in the pacemaker channel HCN4. Front. Physiol. 11, 771.
    doi: 10.3389/fphys.2020.00771pmc: PMC7358946pubmed: 32733276google scholar: lookup
  37. Reitzner SM, Emanuelsson EB, Arif M, Kaczkowski B, Kwon ATJ, Mardinoglu A. Molecular profiling of high-level athlete skeletal muscle after acute endurance or resistance exercise- A systems biology approach. Mol. Metab. 79, 101857.
    doi: 10.1016/j.molmet.2023.101857pmc: PMC10805945pubmed: 38141850google scholar: lookup
  38. Saito Y, Nakamura K, Yoshida M, Sugiyama H, Akagi S, Miyoshi T. Enhancement of pacing function by HCN4 overexpression in human pluripotent stem cell-derived cardiomyocytes. Stem Cell Res. Ther. 13, 141.
    doi: 10.1186/s13287-022-02818-ypmc: PMC8973792pubmed: 35365232google scholar: lookup
  39. Seiradake E, Jones EY, Klein R. Structural perspectives on axon guidance. Annu. Rev. Cell Dev. Bi 32, 577–608.
    doi: 10.1146/annurev-cellbio-111315-125008pubmed: 27576119google scholar: lookup
  40. Shacham EC, Maman N, Lazareva T, Masalha R, Mahagna L, Sela G. Normocalcemic primary hyperparathyroidism is an early stage of primary hyperparathyroidism according to fibroblast growth factor 23 level. Front. Endocrinol. 14, 1152464.
    doi: 10.3389/fendo.2023.1152464pmc: PMC10098304pubmed: 37065752google scholar: lookup
  41. Singh AT, Gilchrist A, Voyno-Yasenetskaya T, Radeff-Huang JM, Stern PH. G alpha12/G alpha13 subunits of heterotrimeric G proteins mediate parathyroid hormone activation of phospholipase D in UMR-106 osteoblastic cells. Endocrinology 146, 2171–2175.
    doi: 10.1210/en.2004-1283pubmed: 15705779google scholar: lookup
  42. Song XT, Zhang JN, Zhao DW, Zhai YF, Lu Q, Qi MY. Molecular cloning, expression, and functional features of IGF1 splice variants in sheep. Endocr. Connect. 10, 980–994.
    doi: 10.1530/EC-21-0181pmc: PMC8428077pubmed: 34319906google scholar: lookup
  43. Stefaniuk-Szmukier M, Ropka-Molik K, Piórkowska K, Zukowski K, Bugno-Poniewierska M. Transcriptomic hallmarks of bone remodelling revealed by RNA-seq profiling in blood of Arabian horses during racing training regime. Gene 676, 256–262.
    doi: 10.1016/j.gene.2018.07.040pubmed: 30021131google scholar: lookup
  44. Tuerxun Z, He YX, Niu YX, Bao Z, Liu XM, Yang YC. Analysis of differentially expressed murine miRNAs in acute myocardial infarction and target genes related to heart rate. Cell Biochem. Biophys. 83, 963–975.
    doi: 10.1007/s12013-024-01528-xpubmed: 39325365google scholar: lookup
  45. Wang SS, Xu Q, Zhang Y, Cao HX, Wang N, Hu YF. The FGF23-Klotho axis promotes microinflammation in chronic kidney disease. Cytokine 184, 156781.
    doi: 10.1016/j.cyto.2024.156781pubmed: 39454251google scholar: lookup
  46. Wilson J, De Donato M, Appelbaum B, Garcia CT, Peters S. Differential expression of innate and adaptive immune genes during acute physical exercise in American quarter horses. Animals 13, 308.
    doi: 10.3390/ani13020308pmc: PMC9854435pubmed: 36670847google scholar: lookup
  47. Yang QJ, Yang GJ, Wan LL, Li B, Lu J, Yu Q. Serum metabolic profiles reveal the effect of formoterol on cachexia in tumor-bearing mice. Mol. Biosyst. 9, 3015–3025.
    doi: 10.1039/c3mb70134dpubmed: 24056883google scholar: lookup
  48. Yang JH, Wu XX, You JH. Unveiling the potential of HSPA4: a comprehensive pan-cancer analysis of HSPA4 in diagnosis, prognosis, and immunotherapy. Aging 16, 2517–2541.
    doi: 10.18632/aging.205496pmc: PMC10911360pubmed: 38305786google scholar: lookup
  49. Zhang ZY, Xu Z, Wang S, Jia ZZ, Zhou Z, Wang C. Optimized new Shengmai Powder modulation of cAMP/Rap1A signaling pathway attenuates myocardial fibrosis in heart failure. Chin. Med-UK 19, 30.
    doi: 10.1186/s13020-024-00902-4pmc: PMC10894496pubmed: 38402401google scholar: lookup

Citations

This article has been cited 1 times.
  1. Huang Q, Ren W, Shan D, Su Y, Li Z, Li L, Wang R, Ma S, Wang J. Molecular Mechanisms Underlying Differences in Athletic Ability in Racehorses Based on Whole Transcriptome Sequencing. Biology (Basel) 2025 Oct 5;14(10).
    doi: 10.3390/biology14101364pubmed: 41154767google scholar: lookup
Subscribe to our newsletter
Join 99,030 horse owners like you who receive our email updates on horse health and nutrition.