FREE SHIPPING over $40
OVER 9000 FIVE-STAR REVIEWS
5% OFF ALL SUBSCRIPTIONS
100% HAPPINESS GUARANTEE
Analyze Diet
0
Home / Equine Research Bank / PLoS One / Finite-Element Analysis of Bone Stresses on Primary Impact i...
PloS one2016; 11(7); e0159541; doi: 10.1371/journal.pone.0159541

Finite-Element Analysis of Bone Stresses on Primary Impact in a Large-Animal Model: The Distal End of the Equine Third Metacarpal.

Abstract: To assess whether the transient stresses of foot impact with the ground are similar to those found during midstance loading and if the location of high stress correlate with the sites most commonly associated with mechanically induced osteoarthritis (OA). We compared impact stresses in subchondral bone between two subject-specific, three-dimensional, finite-element models of the equine metacarpophalangeal (MCP) joint-one with advanced OA and one healthy, and with similar published data on the stresses that occur at midstance. Methods: Two right MCP joints (third metacarpal and proximal phalanx) were scanned using micro-computed tomography (μCT). Images were segmented, and meshed using modified 10-node quadratic tetrahedral elements. Bone material properties were assigned based on the bone density. An impact velocity of 3.55 m/s was applied to each model and contact pressures and stress distribution were calculated for each. In a separate iteration, the third metacarpal was loaded statically. A sampling grid of 160 equidistant points was superimposed over selected slices, and average peak stresses were calculated for 6 anatomical regions. Within-region maximal peak and average von Mises stresses were compared between healthy and OA bones in both midstance and impact loading. Results: Average impact stresses across all regions, in both locations (palmar and dorsal) were greater in the OA model. Highest impact stresses were located in the dorsal medial condyle in the healthy (12.8 MPa) and OA (14.1MPa) models, and were lowest in the palmar medial and lateral parasagittal grooves in the healthy (5.94 MPa) and OA (7.07 MPa) models. The healthy static model had higher peak (up to 49.7% greater) and average (up to 38.6% greater) stresses in both locations and across all regions compared to the OA static model. Conclusions: Under simulated footfall a trot, loading on the dorsal aspect of the third metacarpal at impact created stresses similar to those found during midstance. The high accelerations that occur under impact loading are likely responsible for creating the high stresses, as opposed to midstance loading where the high stresses are the result of high mass loading. Although the stress magnitudes were found to be similar among the two loading conditions, the location of the high stress loading occurred in sites that are not typically associated with osteoarthritic changes.
Get Full Text
Publication Date: 2016-07-26 PubMed ID: 27459189PubMed Central: PMC4961423DOI: 10.1371/journal.pone.0159541Google 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.

The research focuses on investigating whether the temporary stresses experienced when a horse’s foot hits the ground are similar to those during midstance loading and if these high-stress locations relate to the common sites for developing osteoarthritis (OA). The study used computer-based modeling to compare two conditions: one with advanced OA and one healthy.

Methods Utilized

  • The researchers used two right Metacarpophalangeal (MCP) joints from equine specimens for the study. MCP joints used were scanned using micro-computed tomography (μCT), which generated detailed images of the structure.
  • These images were then segmented and meshed using modified 10-node quadratic tetrahedral elements. Here “meshing” relates to the process of subdividing a larger, complex area (like MCP joint) into smaller, simpler parts for easier analysis.
  • Material properties of the bone were assigned based on the bone density, which allow accurate simulations in the model.
  • Each model was subjected to an impact velocity of 3.55 m/s to simulate a horse’s foot hitting the ground, with contact pressures and stress distribution calculated for each model.
  • In a separate iteration, the third metacarpal was loaded statically, that is, without movement, to simulate midstance loading.
  • The stress levels were sampled at 160 equidistant points selected over the area of the bone and average peak stresses were calculated for 6 anatomical regions.
  • The highest peak stresses and average von Mises stresses were then compared between healthy and OA bones under midstance and impact loading.

Results Gathered

  • Results revealed that average impact stresses across all regions, in both locations (palmar and dorsal) were higher in the OA model.
  • The highest stress areas during impact were located in the dorsal medial condyle in both the healthy and OA models, and were lowest in the palmar medial and lateral parasagittal grooves.
  • The healthy static model (representing a horse during midstance) showed higher peak and average stresses across all regions compared to the OA static model.

Conclusion

  • Under simulated footfall at a trot, loading on the dorsal aspect of the third metacarpal bone at impact created stresses similar to those found during midstance.
  • The high accelerations that occur under impact loading could be responsible for creating high stresses, as opposed to midstance loading where the high stresses result from high mass loading.
  • Even though the stress magnitudes were similar among the two loading conditions, the locations of high stress loading occurred in sites that are not typically associated with osteoarthritic changes. This offers new insights into understanding the development of OA conditions in equine models.

Cite This Article

APA
McCarty CA, Thomason JJ, Gordon KD, Burkhart TA, Milner JS, Holdsworth DW. (2016). Finite-Element Analysis of Bone Stresses on Primary Impact in a Large-Animal Model: The Distal End of the Equine Third Metacarpal. PLoS One, 11(7), e0159541. https://doi.org/10.1371/journal.pone.0159541

Publication

PloS one
ISSN: 1932-6203
NlmUniqueID: 101285081
Country: United States
Language: English
Volume: 11
Issue: 7
Pages: e0159541
PII: e0159541

Researcher Affiliations

McCarty, Cristin A
  • Department of Biomedical Science, Ontario Veterinary College, University of Guelph, Guelph, Ontario, Canada.
Thomason, Jeffrey J
  • Department of Biomedical Science, Ontario Veterinary College, University of Guelph, Guelph, Ontario, Canada.
Gordon, Karen D
  • Department of Biomedical Engineering, School of Engineering, University of Guelph, Guelph, Ontario, Canada.
Burkhart, Timothy A
  • Department of Mechanical and Materials Engineering, School of Engineering, Western University, London, Ontario, Canada.
Milner, Jaques S
  • Department of Mechanical and Materials Engineering, School of Engineering, Western University, London, Ontario, Canada.
Holdsworth, David W
  • Robarts Research Institute, Schulich School of Medicine and Dentistry, Western University, London, Ontario, Canada.
  • Department of Surgery, Schulich School of Medicine and Dentistry, Western University, London, Ontario, Canada.

MeSH Terms

  • Animals
  • Biomechanical Phenomena
  • Bone Density
  • Disease Models, Animal
  • Finite Element Analysis
  • Horses
  • Metacarpal Bones
  • Metacarpophalangeal Joint
  • Models, Theoretical
  • Osteoarthritis / diagnostic imaging
  • Osteoarthritis / etiology
  • Osteoarthritis / pathology
  • Pressure
  • Stress, Mechanical
  • X-Ray Microtomography

Conflict of Interest Statement

The authors have declared that no competing interests exist.

References

This article includes 54 references
  1. Harrison SM, Whitton RC, Kawcak CE, Stover SM, Pandy MG. Evaluation of a subject-specific finite-element model of the equine metacarpophalangeal joint under physiological load.. J Biomech 2014 Jan 3;47(1):65-73.
    doi: 10.1016/j.jbiomech.2013.10.001pubmed: 24210848google scholar: lookup
  2. Easton KL. Effect of bone geometry on stress distribution patterns in the equine metacarpophalangeal joint. PhD Thesis, Colorado State University .
  3. Parkin TD, Clegg PD, French NP, Proudman CJ, Riggs CM, Singer ER, Webbon PM, Morgan KL. Catastrophic fracture of the lateral condyle of the third metacarpus/metatarsus in UK racehorses - fracture descriptions and pre-existing pathology.. Vet J 2006 Jan;171(1):157-65.
    pubmed: 16427592doi: 10.1016/j.tvjl.2004.10.009google scholar: lookup
  4. Neundorf RH, Lowerison MB, Cruz AM, Thomason JJ, McEwen BJ, Hurtig MB. Determination of the prevalence and severity of metacarpophalangeal joint osteoarthritis in Thoroughbred racehorses via quantitative macroscopic evaluation.. Am J Vet Res 2010 Nov;71(11):1284-93.
    doi: 10.2460/ajvr.71.11.1284pubmed: 21034319google scholar: lookup
  5. Muir P, Peterson AL, Sample SJ, Scollay MC, Markel MD, Kalscheur VL. Exercise-induced metacarpophalangeal joint adaptation in the Thoroughbred racehorse.. J Anat 2008 Dec;213(6):706-17.
    pmc: PMC2666139pubmed: 19094186doi: 10.1111/j.1469-7580.2008.00996.xgoogle scholar: lookup
  6. Drum MG, Kawcak CE, Norrdin RW, Park RD, McIlwraith CW, Les CM. Comparison of gross and histopathologic findings with quantitative computed tomographic bone density in the distal third metacarpal bone of racehorses.. Vet Radiol Ultrasound 2007 Nov-Dec;48(6):518-27.
    pubmed: 18018722doi: 10.1111/j.1740-8261.2007.00289.xgoogle scholar: lookup
  7. Radin EL, Paul IL. Response of joints to impact loading. I. In vitro wear.. Arthritis Rheum 1971 May-Jun;14(3):356-62.
    pubmed: 5562019doi: 10.1002/art.1780140306google scholar: lookup
  8. Young BD, Samii VF, Mattoon JS, Weisbrode SE, Bertone AL. Subchondral bone density and cartilage degeneration patterns in osteoarthritic metacarpal condyles of horses.. Am J Vet Res 2007 Aug;68(8):841-9.
    pubmed: 17669024doi: 10.2460/ajvr.68.8.841google scholar: lookup
  9. Norrdin RW, Kawcak CE, Capwell BA, McIlwraith CW. Calcified cartilage morphometry and its relation to subchondral bone remodeling in equine arthrosis.. Bone 1999 Feb;24(2):109-14.
    pubmed: 9951778doi: 10.1016/s8756-3282(98)00157-4google scholar: lookup
  10. Riggs CM. Aetiopathogenesis of parasagittal fractures of the distal condyles of the third metacarpal and third metatarsal bones--review of the literature.. Equine Vet J 1999 Mar;31(2):116-20.
    pubmed: 10213423doi: 10.1111/j.2042-3306.1999.tb03803.xgoogle scholar: lookup
  11. Radin EL, Martin RB, Burr DB, Caterson B, Boyd RD, Goodwin C. Effects of mechanical loading on the tissues of the rabbit knee.. J Orthop Res 1984;2(3):221-34.
    pubmed: 6436458doi: 10.1002/jor.1100020303google scholar: lookup
  12. Gustas P, Johnston C, Drevemo S. Ground reaction force and hoof deceleration patterns on two different surfaces at the trot. Equine and Comparative Exercise Physiology 3(4):209–16.
  13. Rubio-Martínez LM, Cruz AM, Gordon K, Hurtig MB. Mechanical properties of subchondral bone in the distal aspect of third metacarpal bones from Thoroughbred racehorses.. Am J Vet Res 2008 Nov;69(11):1423-33.
    doi: 10.2460/ajvr.69.11.1423pubmed: 18980424google scholar: lookup
  14. Cruz AM, Hurtig MB. Multiple pathways to osteoarthritis and articular fractures: is subchondral bone the culprit?. Vet Clin North Am Equine Pract 2008 Apr;24(1):101-16.
    doi: 10.1016/j.cveq.2007.12.001pubmed: 18314038google scholar: lookup
  15. Muir P, McCarthy J, Radtke CL, Markel MD, Santschi EM, Scollay MC, Kalscheur VL. Role of endochondral ossification of articular cartilage and functional adaptation of the subchondral plate in the development of fatigue microcracking of joints.. Bone 2006 Mar;38(3):342-9.
    pubmed: 16275175doi: 10.1016/j.bone.2005.08.020google scholar: lookup
  16. Hinterhofer C, Stanek C, Haider H. Finite element analysis (FEA) as a model to predict effects of farriery on the equine hoof.. Equine Vet J Suppl 2001 Apr;(33):58-62.
    pubmed: 11721570doi: 10.1111/j.2042-3306.2001.tb05360.xgoogle scholar: lookup
  17. McClinchey HL, Thomason JJ, Jofriet JC. Isolating the effects of equine hoof shape measurements on capsule strain with finite element analysis. Veterinary and Comparative Orthopaedics And Traumatology 16: 67–75.
  18. O'Hare LM, Cox PG, Jeffery N, Singer ER. Finite element analysis of stress in the equine proximal phalanx.. Equine Vet J 2013 May;45(3):273-7.
    doi: 10.1111/j.2042-3306.2012.00635.xpubmed: 22943561google scholar: lookup
  19. Les CM, Keyak JH, Stover SM, Taylor KT. Development and validation of a series of three-dimensional finite element models of the equine metacarpus.. J Biomech 1997 Jul;30(7):737-42.
    pubmed: 9239555doi: 10.1016/s0021-9290(97)00007-9google scholar: lookup
  20. Collins SN, Murray RC, Kneissl S, Stanek C, Hinterhofer C. Thirty-two component finite element models of a horse and donkey digit.. Equine Vet J 2009 Mar;41(3):219-24.
    pubmed: 19469224doi: 10.2746/042516409x395688google scholar: lookup
  21. Back W, Schamhardt HC, Hartman W, Barneveld A. Kinematic differences between the distal portions of the forelimbs and hind limbs of horses at the trot.. Am J Vet Res 1995 Nov;56(11):1522-8.
    pubmed: 8585667
  22. Clayton HM. The effect of an acute hoof wall angulation on the stride kinematics of trotting horses.. Equine Vet J Suppl 1990 Jun;(9):86-90.
    pubmed: 9259814doi: 10.1111/j.2042-3306.1990.tb04742.xgoogle scholar: lookup
  23. Linford RL. Camera speeds for optoelectronic assessment of stride-timing characteristics in horses at the trot.. Am J Vet Res 1994 Sep;55(9):1189-95.
    pubmed: 7802382
  24. Mansour JM. Biomechanics of cartilage. Kinesiology: the mechanics and pathomechanics of human movement 66.–.
  25. Laasanen MS, Töyräs J, Korhonen RK, Rieppo J, Saarakkala S, Nieminen MT, Hirvonen J, Jurvelin JS. Biomechanical properties of knee articular cartilage.. Biorheology 2003;40(1-3):133-40.
    pubmed: 12454397
  26. Oloyede A, Flachsmann R, Broom ND. The dramatic influence of loading velocity on the compressive response of articular cartilage.. Connect Tissue Res 1992;27(4):211-24.
    pubmed: 1576822doi: 10.3109/03008209209006997google scholar: lookup
  27. Shepherd DE, Seedhom BB. The 'instantaneous' compressive modulus of human articular cartilage in joints of the lower limb.. Rheumatology (Oxford) 1999 Feb;38(2):124-32.
    pubmed: 10342624doi: 10.1093/rheumatology/38.2.124google scholar: lookup
  28. Li LP, Buschmann MD, Shirazi-Adl A. Strain-rate dependent stiffness of articular cartilage in unconfined compression.. J Biomech Eng 2003 Apr;125(2):161-8.
    pubmed: 12751277doi: 10.1115/1.1560142google scholar: lookup
  29. Park S, Hung CT, Ateshian GA. Mechanical response of bovine articular cartilage under dynamic unconfined compression loading at physiological stress levels.. Osteoarthritis Cartilage 2004 Jan;12(1):65-73.
    pubmed: 14697684doi: 10.1016/j.joca.2003.08.005google scholar: lookup
  30. Langelier E, Buschmann MD. Increasing strain and strain rate strengthen transient stiffness but weaken the response to subsequent compression for articular cartilage in unconfined compression.. J Biomech 2003 Jun;36(6):853-9.
    pubmed: 12742453doi: 10.1016/s0021-9290(03)00006-xgoogle scholar: lookup
  31. Austman RL, Milner JS, Holdsworth DW, Dunning CE. The effect of the density-modulus relationship selected to apply material properties in a finite element model of long bone.. J Biomech 2008 Nov 14;41(15):3171-6.
    doi: 10.1016/j.jbiomech.2008.08.017pubmed: 18922532google scholar: lookup
  32. McErlain DD, Milner JS, Ivanov TG, Jencikova-Celerin L, Pollmann SI, Holdsworth DW. Subchondral cysts create increased intra-osseous stress in early knee OA: A finite element analysis using simulated lesions.. Bone 2011 Mar 1;48(3):639-46.
    doi: 10.1016/j.bone.2010.11.010pubmed: 21094285google scholar: lookup
  33. Les CM, Keyak JH, Stover SM, Taylor KT, Kaneps AJ. Estimation of material properties in the equine metacarpus with use of quantitative computed tomography.. J Orthop Res 1994 Nov;12(6):822-33.
    pubmed: 7983558doi: 10.1002/jor.1100120610google scholar: lookup
  34. Lanovaz JL, Clayton HM, Watson LG. In vitro attenuation of impact shock in equine digits.. Equine Vet J Suppl 1998 Sep;(26):96-102.
    pubmed: 9932099doi: 10.1111/j.2042-3306.1998.tb05127.xgoogle scholar: lookup
  35. Wilson AM, McGuigan MP, Su A, van Den Bogert AJ. Horses damp the spring in their step.. Nature 2001 Dec 20-27;414(6866):895-9.
    pubmed: 11780059doi: 10.1038/414895agoogle scholar: lookup
  36. Gustås P, Johnston C, Roepstorff L, Drevemo S. In vivo transmission of impact shock waves in the distal forelimb of the horse.. Equine Vet J Suppl 2001 Apr;(33):11-5.
    pubmed: 11721549doi: 10.1111/j.2042-3306.2001.tb05350.xgoogle scholar: lookup
  37. McCarty CA, Thomason JJ, Gordon K, Burkhart T, Bignell W. Effect of hoof orientation and ballast on acceleration and vibration in the hoof and distal forelimb following simulated impacts ex vivo.. Equine Vet J 2015 Mar;47(2):223-9.
    doi: 10.1111/evj.12252pubmed: 24580552google scholar: lookup
  38. Noble P, Lumay G, Coninx M, Collin B, Magnée A, Lecomte-Beckers J, Denoix JM, Serteyn D. A pendulum test as a tool to evaluate viscous friction parameters in the equine fetlock joint.. Vet J 2011 May;188(2):204-9.
    pubmed: 20413334doi: 10.1016/j.tvjl.2010.03.016google scholar: lookup
  39. Gustås P, Johnston C, Roepstorff L, Drevemo S, Lanshammar H. Relationships between fore- and hindlimb ground reaction force and hoof deceleration patterns in trotting horses.. Equine Vet J 2004 Dec;36(8):737-42.
    pubmed: 15656507doi: 10.2746/0425164044848136google scholar: lookup
  40. Brama PA, Karssenberg D, Barneveld A, van Weeren PR. Contact areas and pressure distribution on the proximal articular surface of the proximal phalanx under sagittal plane loading.. Equine Vet J 2001 Jan;33(1):26-32.
    pubmed: 11191606doi: 10.2746/042516401776767377google scholar: lookup
  41. Easton KL, Kawcak CE. Evaluation of increased subchondral bone density in areas of contact in the metacarpophalangeal joint during joint loading in horses.. Am J Vet Res 2007 Aug;68(8):816-21.
    pubmed: 17669020doi: 10.2460/ajvr.68.8.816google scholar: lookup
  42. Carter DR, Hayes WC. The compressive behavior of bone as a two-phase porous structure.. J Bone Joint Surg Am 1977 Oct;59(7):954-62.
    pubmed: 561786
  43. Linde F, Nørgaard P, Hvid I, Odgaard A, Søballe K. Mechanical properties of trabecular bone. Dependency on strain rate.. J Biomech 1991;24(9):803-9.
    pubmed: 1752864doi: 10.1016/0021-9290(91)90305-7google scholar: lookup
  44. McCarty CA, Thomason JJ, Gordon K, Hurtig M, Bignell W. Effect of hoof angle on joint contact area in the equine metacarpophalangeal joint following simulated impact loading ex vivo.. Equine Vet J 2015 Nov;47(6):715-20.
    pubmed: 25237779doi: 10.1111/evj.12354google scholar: lookup
  45. Schaffler MB, Burr DB. Stiffness of compact bone: effects of porosity and density.. J Biomech 1988;21(1):13-6.
    pubmed: 3339022doi: 10.1016/0021-9290(88)90186-8google scholar: lookup
  46. Harrison SM, Whitton RC, King M, Haussler KK, Kawcak CE, Stover SM, Pandy MG. Forelimb muscle activity during equine locomotion.. J Exp Biol 2012 Sep 1;215(Pt 17):2980-91.
    doi: 10.1242/jeb.065441pubmed: 22875767google scholar: lookup
  47. Drum MG, Les CM, Park RD, McIlwraith CW, Kawcak CE. Comparison of mean bone densities of three preparations of the distal portion of the equine third metacarpal bone measured by use of quantitative computed tomography.. Am J Vet Res 2008 Jul;69(7):891-3.
    doi: 10.2460/ajvr.69.7.891pubmed: 18593241google scholar: lookup
  48. Whitton RC, Trope GD, Ghasem-Zadeh A, Anderson GA, Parkin TD, Mackie EJ, Seeman E. Third metacarpal condylar fatigue fractures in equine athletes occur within previously modelled subchondral bone.. Bone 2010 Oct;47(4):826-31.
    doi: 10.1016/j.bone.2010.07.019pubmed: 20659599google scholar: lookup
  49. Ozcivici E, Garman R, Judex S. High-frequency oscillatory motions enhance the simulated mechanical properties of non-weight bearing trabecular bone.. J Biomech 2007;40(15):3404-11.
    pubmed: 17655852doi: 10.1016/j.jbiomech.2007.05.015google scholar: lookup
  50. Burr DB, Radin EL. Microfractures and microcracks in subchondral bone: are they relevant to osteoarthrosis?. Rheum Dis Clin North Am 2003 Nov;29(4):675-85.
    pubmed: 14603577doi: 10.1016/s0889-857x(03)00061-9google scholar: lookup
  51. Frost HM. Wolff's Law and bone's structural adaptations to mechanical usage: an overview for clinicians.. Angle Orthod 1994;64(3):175-88.
    pubmed: 8060014doi: 10.1043/0003-3219(1994)064<0175:wlabsa>2.0.co;2google scholar: lookup
  52. Kawcak CE, McIlwraith CW. Proximodorsal first phalanx osteochondral chip fragmentation in 336 horses.. Equine Vet J 1994 Sep;26(5):392-6.
    pubmed: 7988543doi: 10.1111/j.2042-3306.1994.tb04409.xgoogle scholar: lookup
  53. Colahan P, Turner TA, Poulos P, Piotrowski G. Mechanical functions and sources of injury in the fetlock and carpus. Proceedings of the annual convention of the American Association of Equine Practitioners 689–699.
  54. Thomason JJ, Peterson ML. Biomechanical and mechanical investigations of the hoof-track interface in racing horses.. Vet Clin North Am Equine Pract 2008 Apr;24(1):53-77.
    doi: 10.1016/j.cveq.2007.11.007pubmed: 18314036google scholar: lookup

Citations

This article has been cited 5 times.
  1. Lopas LA, Shen H, Zhang N, Jang Y, Tawfik VL, Goodman SB, Natoli RM. Clinical Assessments of Fracture Healing and Basic Science Correlates: Is There Room for Convergence?. Curr Osteoporos Rep 2023 Apr;21(2):216-227.
    doi: 10.1007/s11914-022-00770-7pubmed: 36534307google scholar: lookup
  2. Pramudita JA, Hiroki W, Yoda T, Tanabe Y. Variations in Strain Distribution at Distal Radius under Different Loading Conditions.. Life (Basel) 2022 May 16;12(5).
    doi: 10.3390/life12050740pubmed: 35629407google scholar: lookup
  3. Harrison SM, Whitton RC, Stover SM, Symons JE, Cleary PW. A Coupled Biomechanical-Smoothed Particle Hydrodynamics Model for Horse Racing Tracks.. Front Bioeng Biotechnol 2022;10:766748.
    doi: 10.3389/fbioe.2022.766748pubmed: 35265590google scholar: lookup
  4. Johnston GCA, Ahern BJ, Palmieri C, Young AC. Imaging and Gross Pathological Appearance of Changes in the Parasagittal Grooves of Thoroughbred Racehorses.. Animals (Basel) 2021 Nov 24;11(12).
    doi: 10.3390/ani11123366pubmed: 34944142google scholar: lookup
  5. Sakashita M, Yamasaki S, Yaji K, Kawamoto A, Kondo S. Three-dimensional topology optimization model to simulate the external shapes of bone.. PLoS Comput Biol 2021 Jun;17(6):e1009043.
    doi: 10.1371/journal.pcbi.1009043pubmed: 34133416google scholar: lookup
Subscribe to our newsletter
Join 99,344 horse owners like you who receive our email updates on horse health and nutrition.