So wa e M code o modelling i us dynamics in
mul ispecies ecosys ems
Familusi Oluwa osin Adekunlea,b,c, Go ka Bidegainb,c, Tal Ben-Ho ind
aIns i u des Sciences de la Me , L’Uni e si ´e du Qu´ebec `a Rimouski, All´ee des
U sulines, Rimouski, 3300, Qc, Canada
bDepa men o Applied Ma hema ics, Uni e si y o he Basque Coun y
(UPV/EHU), Plaza Eu opa 1, Donos ia, 20018, Gipuzkoa, Spain
cResea ch Cen e o Expe imen al Ma ine Biology and Bio echnology, Plen zia Ma ine
S a ion, Uni e si y o he Basque Coun y (PiE-UPV/EHU), A ea za
Pasealekua, Plen zia, 48620, Bizkaia, Spain
dNo h Ca olina S a e Uni e si y, College o Ve e ina y Medicine, 303, College
Ci cle, Mo ehead Ci y, NC 28557, No h Ca olina, USA
Abs ac
We p esen an M-code (MATLAB) modeling oolki o wa e bo ne i-
al ansmission in mixed-species communi ies. The amewo k couples hos
p ocesses—g ow h, backg ound mo ali y, and il a ion/clea ance—wi h i-
al p ocesses—shedding/ elease, en i onmen al anspo , and decay—and
suppo s bo h mono- and co-cul u e con igu a ions. In con ex s lacking co-
cul u e in ec ion da a, pa ame e s can be calib a ed o monocul u e obse -
a ions and hen explo ed ac oss species mix u es. We e i y ma hema ical
consis ency and ealis ic beha io ac oss wide composi ion g adien s (e.g.,
90/10 o 10/90). Using an oys e –he pes i us (OsHV-1) case s udy wi h M.
gigas and C. i ginica, simula ions show ha adding a second, mo e ole -
an hos does no hal epidemic p og ession in he suscep ible hos ; peak
mo ali y s ill occu s bu is delayed, consis en wi h sa u a ion o i al pa i-
cles in he wa e column ha limi s dilu ion bene i s. The oolki eplica es
obse ed monocul u e dynamics and gene alizes o mul ispecies scena ios, en-
abling hypo hesis es ing abou species in oduc ions, husband y s a egies,
and en i onmen al con ols on ansmission. O e all, i p o ides a anspa -
en , ex ensible pla o m o in es iga ing pa hogen dynamics in mul ispecies
ecosys ems and o e alua ing managemen op ions in aquacul u e and be-
yond.
Keywo ds: Modelling, Os eid He pes i us 1, co-cul u e, C. i ginica,M.
Oc obe 16, 2025
gigas
1. Publica ion1
This code was de eloped and applied o he s udy ”Cocul u e wi h Eas -2
e n oys e s is unlikely o educe OsHV-1 impac s o a med Paci ic oys e s: A3
modelling app oach” published in Aquacul u e Repo s (Else ie ) h ps://4
www.sciencedi ec .com/science/a icle/pii/S2352513424006550. This5
modelling wo k and pa ame e es ima es a e based on empi ical esea ch.6
2
