arxiv: v1 [q-bio.pe] 14 Sep 2016

Transcription

1 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in axiv:69.434v [q-io.pe] 4 Sep 6 A EPIDEIOLOGICAL ODEL OF VIRAL IFECTIOS I A VARROA-IFESTED EE COLOY: THE CASE OF A EE-DEPEDET ITE POPULATIO SIZE SARA ERARDI, EZIO VETURIO Dipatimento di atematica Giuseppe Peano, via Calo Aleto, Univesità di Toino, Italy s.enadi@unito.it, ezio.ventuino@unito.it In ecent yeas te spead of te ectopaasitic mite Vaoa destucto as ecome te most seious teat to woldwide apicultue. In te model pesented ee we extend te ee population dynamics wit mite vial epidemiology examined in an ealie pape y allowing a ee-dependent mite population size. Te esults of te analysis matc field osevations well and give a clea explanation of ow Vaoa affects te epidemiology of cetain natually occuing ee viuses, causing consideale damages to colonies. Te model allows only fou possile stale equiliia, using known field paametes. Te fist one contains only te tiving ealty ees. Hee te disease is eadicated and also te mites ae wiped out. Altenatively, we find te equiliium still wit no mite population, ut wit endemic disease among te tiving ee population. Tidly, infected ees coexist wit te mites in te Vaoa invasion scenaio; in tis situation te disease invades te ive, diving te ealty ees to extinction and teefoe affecting all te ees. Coexistence is also possile, wit ot populations of ees and mites tiving and wit te disease endemically affecting ot species. Te analysis is in line wit field osevations in natual oney ee colonies. amely, tese diseases ae endemic and if te mite population is pesent, necessaily te vial infection occus. Fute, in ageement wit te fact tat te pesence of Vaoa incements te vial tansmission, te wole ee population may ecome infected wen te disease vecto is pesent in te eeive. Also, a low oizontal tansmission ate of te vius among te oney ees will elp in potecting te ee colonies fom Vaoa infestation and vial epidemics.. Intoduction Cuent evidence sows tat Eat loses etween one and ten pecent of iodivesity pe decade, a majo event of extinction of iological divesity. Te causes ae mainly aitat loss, pest invasions, pollution, ove-avesting and diseases,, teey endangeing natual ecosystem sevices tat ae vital fo umanity.

2 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in All poducts of agicultue depend on pollination, tat is pefomed y wild, fee-living oganisms suc as ees, utteflies, mots and flies. To tis end, tee is also te availaility of commecially-managed ee species. ut ees epesent te most impotant goup of pollinatos in most geogapical egions, also fo economical easons,. An estimation of te United ations Food and Agicultue Oganisation (FAO) sets te ee-pollinated cops at 7 out of te species poviding 9% of food woldwide. In Euope, tee ae 64 cop species and 84% of tem ae animal-pollinated; fute, 4 vegetale vaieties tive due to te pollination activities of ees,. As te ee goup is te most impotant pollinato woldwide, tis pape focuses on te instaility of ee populations and, in paticula, on te te most seious teat to apicultue gloally: te extenal paasitic mite Vaoa destucto, discoveed in Souteast Asia in 94. Tis mite is of te size of a pinead; it feeds on te ees ciculatoy fluid. It is an invasive species, speading fom one ive to anote one. Today it is pesent nealy in te wole wold,. Te seious damage to te ee colonies does not deive fom te paasiticaction ofte mite ut, aoveall, fom its action as vecto of many vial diseases. It inceases te tansmission ate of diseases suc as acute paalysis (APV) and defomed wing viuses (DWV), tat ae consideed among te main causes of Colony Collapse Disode (CCD). Tis action as een einfoced since aout tity yeas ago wen te mite as sifted osts fom Apis ceana to Apis mellifea. Wit no contol, te infestations ae ound to cause te untimely deat of ee colonies witin tee yeas. atematical modelling epesents a poweful tool fo investigating te tiangula elationsip etween oney ees, Vaoa and vial disease. It allows te exploation of te eeive system witout te need of unfeasile o costly field studies. In paticula, matematical models of te epidemiology of vial diseases on Vaoa-infested colonies would allow us to exploe te ost population esponses to suc stessful situations. In te next Section, te model is pesented and some asic popeties ae investigated. In Section 3, equiliia ae analysed fo feasiility and staility. A sensitivity analysis on te model paametes is pefomed in Section 4.4. A final discussion concludes te study.

3 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in 3. Te model Let denote ealty ees, I te infected ones, te ealty mites and te infected ones and let te population vecto e X = (,I,,) T. In te model poposed ee we essentially extend te epidemiology of a eeive infested y Vaoa mites examined ealie in, to account fo a Leslie Gowe tem in te mite populations, to model te Vaoa population size as a ee population-dependent function, a step also undetaken in te ecent pape 9. In fact, te mites essentially cayout tei life on te osts ody. Tese canges ae teefoe eflected in te last two equations, fo te mites evolution. Te esulting model eads as follows: wee d dt = f = +I di dt = f I I = +I d dt X = f(x) = (f,f I,f,f ) T, f : D R 4 +, () λ γi m = f = ( +) n +λ +γi (m+)i p ( +) (βi δ e) ( +I) d dt = f p = n ( +)+βi +δ e. ( +I) To pevent in te igt-and side of () te vanising of some tems in te denominato, we define D as te domain of f, explicitly D = { X = (,I,,) R 4 + : +I }. () Te fist two equations descie espectively te evolution of ealty and infected ee populations. Tey ae on ealty o infected in popotion to te faction of ealty o infected ees in te colony, wit constant. Te infection pocess fo te lavae indeed occus mainly toug tei meals of contaminated oyal jelly. Healty ees can contact te vius y infected mites at ate λ, second tems in te equations. Te next tems model oizontal tansmission of te vius among adult ees, occuing at ate γ via small wounds of te exoskeleton, e.g. as a esult of ai loss, o ingestion of faeces. In te last tems we find te ees natual motality m and te disease-elated motality, taken into account fo te infected ee population only. Te last two equations contain te Vaoa dynamics, patitioned among susceptile and infected states, tat vecto te vial disease and once infected emain so fo tei lifetime. Te mites ae always on ealty at

4 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in 4 ate, as no viuses ae passed vetically fom te paents to offspings. In te extended model te fist tee tems descie Vaoa gowt, natual motality and intaspecific competition. ote tat instead in te peviously intoduced system,, te mites gow logistically. Hee if te mites epoduction ate is smalle tan tei natual motality n, i.e. < n, (3) it is easy to see tat te total mite population ecomes extinct. Te effect of intaspecific competition among mites is descied y te Leslie Gowe tem at ate p. It accounts fo te mite population dependence on te total nume of adult ees, te esouce fo wic tey compete in te colony, 8. Te paamete expesses te aveage nume of mites pe ee. ote tat ecause te vius does not cause any am to te infected Vaoa population, competition among ealty and infected mites occus at te same ate, in ote wods we find te tems ( + ) and ( + ) wit te same weigt in ot equations. If infected wee weakened y te disease, we would ate ave in te acket +c, wit some c <. Healty mites can ecome infected y infected ees at ate β, te fout tems in te last two equations, ut tey can also acquie te vius at ate δ also oizontally fom ote infected mites. Te last tem in te last two equations models te gooming eavio of ealty ees at ate e. We assume tat te infected ees do not goom ecause tey ae weake due to te disease effects. ote tat fo ealty and infected mites te damage due to gooming is te same, e, as tis activity depends only on te ees... Well-posedness and oundedness We now addess te issue of well-posedness, following asically te pat of. Te solutions tajectoies ae sown to e always at a finite distance fom te set Θ = { (,I,,) R 4 + : +I = }, wee te total ee population vanises. Teoem.. Well-posedness and oundedness Let X D. A solution of () defined on [,+ ) wit X() = X exists uniquely. Also, fo an aitay t >, it follows X(t) D, and, indicating y L a positive constant, m liminf ((t)+i(t)) limsup((t)+i(t)), m = m+; (4) t + t + m (t)+(t) L, t. ()

5 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in Poof. We follow te aguments of, fo te fist pat of te poof. Gloal Lipscitz continuity of te igt-and side of te system olds in D, teey implying existence and uniqueness of te solution of system () fo evey tajectoy at a finite distance of tis ounday. ow (4) and () old fo all te tajectoies stating at any point wit +I. Te oundedness of te vaiales entail consequently tat all tajectoies exist at all times in te futue and ae ounded away fom te set Θ. Adding te two equations fo te ees in (), at any point in D we ave +I = m (m+)i m( +I), m = m+. Integation of tis diffeential inequality etween X() = X and X(t), points elonging to a tajectoy suc tat X(τ) D fo all τ [,t], we otain te following lowe ound, wit positive igt and side at any t >, (t)+i(t) m ( e mt )+(()+I())e mt. (6) Similaly, ut ounding fom aove, we find + I m( + I), and teefoe (t)+i(t) m ( e mt )+(()+I())e mt. (7) Te inequalities in (4) fo any potion of te tajectoy elonging to D follow fom (6) and (7). We now tun to te Vaoa populations and, and tei te wole mite population V = +, summing tei coesponding equations of () and using (4), fo any point witin D we ave te uppe ound [ n+α p ( +I) V p ( +I) V ev (8) [ n+α pm ] V V = φ(v) φ(v), V +αv = ( n+α)v ] V wee φ(v) epesents te maximum value of te paaola φ(v), namely V = ( n+α)(pm). A consequence of (8) is te diffeential inequality V φ(v) αv and tus integating it we otain te following inequality, fo all t, poving (), V(t) φ(v) α ( e αt ) { +V()e αt max V(), φ(v) }. (9) α

6 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in 6 Fom (6), (7) and (9) te oundedness of all populations is tus ensued,. Tus all tajectoies oiginating in D emain in D fo all t >. Let D e te lagestsuset of D satisfying te inequalities ofteoem (.), { D. = (,I,,) R 4 + : m +I }, + L. m As a consequence of Teoem., it is a compact set, positively invaiant fo all te system s tajectoies. All te system s equiliia elong to D, as it is sown late. Tus te equiliia analysis is sufficient to explain te wole system s eavio. 3. Equiliia Te equiliia E k = ( k,i k, k, k ) of () ae ( ) ( ) E = m,,,, E =, m+,,, wic ae always feasile; te points ( ( E 3 = m,, e ) ) ( m m, γ +m, E 4 =, γ te fome feasile fo ) γ m,,, γ e m, () te latte feasile weneve te following condition is satisfied < γ m <. () Ten, tee is te coexistence equiliium E = (,I,, ). Howeve, it is not analytically tactale and will teefoe e investigated numeically, wit te elp of simulations. Finally, we find te equiliium E = (,I,, ) wit no ealty ees. It is otained as intesection of two conic sections tat lie in te fist quadant of te - pase plane. We discuss it in te next susection. 3.. Te ealty-ees-fee equiliium Fom te second equation we find I = (m+). Sustituting te value of I, te last two equations can e ewitten as

7 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in 7 and ψ : ( n) + η : n p(m+)( +) p(m+)( +) β δ =, () m+ + β +δ =. (3) m+ Tus, te equiliium follows fom detemining te intesection of ψ and η in te fist quadant of te pase plane. We egin y analyzing te cuve ψ. Solving () fo te vaiale, we find [ p(m+) +β ] +(n )(m+) ψ : = = (m+)( p(m+) δ) ψ() (4) Te conic ψ cosses te vetical axis at te oigin and te oizontal one at te asissa = [( n)(m+) β] p(m+) () Fute, tee is te vetical asymptote = >. (6) p(m+)+δ Fom (), it is easy to assess te signs of te numeato and te denominato, espectively > fo > and D > fo <. Depending on te sign of and D, ψ tus sows tee diffeent sapes, Figue. We now tun to te second conic section, η. Reaanging (3) fo, we get te explicit fom [n+p(m+)](m+) = β +[δ(m+) p(m+) ] = η() Poceeding as efoe, we detemine te intesections of η wit te axes. Te conic η goes toug te oigin and cosses te axis at Again, we get one vetical asymptote = = n p(m+) <. β (m+)[p(m+) δ]. ow fom (3), te sign of te numeato is positive, >, fo >, wile fo te denominato D two cases aise. amely, if δ < p(m+),

8 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in 8 Hypeola ψ Figue. Left: Case. < <, fo te paametes values =.6, e =., λ =.3, =, β =., m =.3, = 8, δ =.4, p =., =.9, n =.7. Cente: Case. < <, fo te paametes values =., e =., λ =.3, =, β =., m =.3, = 8, δ =., p =., =.9, n =.. Rigt: Case 3. < <, fo te paametes values =.6, e =., λ =.3, =, β =., m =.4, = 8, δ =.. i.e. >, we get D > fo <. Otewise, if δ > p(m+), i.e. <, to ave D > te opposite condition > must old. Figue sums up te tee possile sapes fo η. Hypeola η Figue. Left: Case A. < <, fo te paametes values =., e =., λ =.3, =, β =., m =.3, = 8, δ =., p =., =, n =.. Cente: Case. < <, fo te paametes values =., e =., λ =.3, =, β =., m =.3, = 8, δ =., p =., =, n =.. Rigt: Case C. < <, fo te paametes values =., e =., λ =.3, =, β =., m =.3, = 8, δ =., p =., =, n =.. Finally, y gapically plotting ot te ypeolae in te - plane, teymeetatteoiginandteyfuteintesectatanotepointlocatedin te fist quadant. Tis immediately povides te unconditional existence and feasiility of te equiliium point E.

9 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in Staility Te Jacoian J of () is te following matix J γ λ ( +I) J +I I +γ (m+) λ ( +I) p( +) J 3 ( +I) β J 33 p ( +I) δ p( +) J 4 ( +I) +β J 43 n p(+) (+I) +δ e wit J = +I J 3 = ( +I) λ γi m, J = I +λ +γi, ( +I) p( +) βi δ e, ( +I) p( +) ( +I) e, J 33 = n J 4 = p( +) ( +I) e, J 43 = p +βi +δ. ( +I) AtE twoeigenvaluesoftetejacoianevaluatedattispoint, J(E ), ae negative, m and n em. Te emaining two povidete staility conditions γ < m, m < e+nm. (7) At equiliium E te Jacoian J(E ) as two explicit eigenvalues, (m + ) < and γ(m + ) wile te Rout-Huwitz conditions on te emaining mino sow tat if (3) is not satisfied, te equiliium is unstale, ecause in suc case det(j(e )) <, t(j(e )) = n+ β m+, det(j(e )) = (n ) ( n+ β ). m+ Te Jacoian J(E 3 ) once again gives two explicit eigenvalues, m < and anote one poviding te fist staility condition < n + p 3 ( 3 ) +e 3, i.e. explicitly e < m( n), (8) and te Rout-Huwitz citeion on te emaining mino J 3 gives te fute conditions det( J 3 ) = t( J 3 ) = γ 3 +m++n+ p 3 δ 3 +e 3 >, (9) 3 ( )( 3 +γ 3 m δ 3 e 3 n p ) 3 βλ 3 3 >. () 3 3

10 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in ote tat if (3) olds, E 3 is unstale, ecause in tis case (8) cannot e satisfied. Atte pointe 4 te caacteisticequationfactoizesintote poductof two quadatic equations, fo wic te Rout-Huwitz conditions povide te following pais of inequalities, to e satisfied fo staility and ( + m++γi 4 > γ 4 +, () 4 +I 4 4 I 4 [γ ] ( ) ( 4 +I 4 ) 4 + m 4 +I 4 )[ ] I 4 m γ 4 4 +I 4 ( 4 +I 4 ) 4 ( 4 +I 4 ) > ( γi ( 4 +I 4 ) )[ γ 4 I 4 ( 4 +I 4 ) ], () n+e 4 +βi 4 >, (n+βi 4 +e 4 )(n+e 4 ) > [(n+e 4 )+βi 4 ]. (3) Two eigenvalues ae also explicitly found at E, m < and te ote one giving te fist staility condition I < γi +m+λ. (4) Te Rout-Huwitz conditions fo staility on te emaining quadatic ecome ten n+ 3p( + ) +βi +δ > +δ, () I ( n p( )( + ) βi δ δ n p( ) + ) I I ( + βi p )( ) p +δ +δ >. I I 4. Results Afte desciing te set of paamete values used, te cosen initial conditions and te field data availale, we pefom te sensitivity analysis. Te simulations ave een pefomed using te atla uilt-in odinay diffeential equations solve ode4.

11 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in 4.. odel paametes fom te liteatue Te model paametes tat ae known fom te liteatue, 6, 4, and 3, 8..3., ae fixed in te simulations, wile te emaining ones ae canged ove a suitale ange. In te sensitivity analysis oweve we will vay also te known ones, to simulate possile envionmental vaiations, due peaps to climatic canges o ote extenal disuptions. We set te time unit to e te day. Te woke oney ees it ate is =, tei natual motality ate instead is m =.3, wic implies a life expectancy of 43. days in te adult stage, 6. Tee ae no pecise values fo te gooming eavio in te liteatue. A possile ange of e is pesumed to e in te inteval [ 6, ], 4. Te Vaoa population epoduces exponentially fast, douling evey mont duing te sping and summe. We ten take 3 ln, i.e. =., 3, p.. In te same season, te mite natual motality ate in te poetic pase, i.e. wen attaced to an adult ee, is pesumed to ave te value n =.7, In Tale we list all te efeence values fo te numeical expeiments. Te ote model paametes ae feely cosen, wit ypotetical values. Tale. Fo tese model paametes te values ae otained fom te liteatue. ee and mite populations ae measued in pue numes. Paamete Intepetation Value Unit Souce ee daily it ate day 6 e Healty ee gooming ate 6 day 4 m ee natual motality ate.3 day 6 Vaoa gowt ate. day 3 n Vaoa natual motality.7 day ate in te poetic pase 4.. Setting of initial conditions and fee paametes Te colony conditions at te eginning of te sping come fom field data: all te ees ae ealty. Indeed te infected ones do not suvive te winte, ecause tey ave a lowe life expectancy. In addition, te colony teatments wit acaicide ae usually pefomed in te late autumn, to allow te mite eadication. We can safely assume ten te mite population is aound units at te stat of te sping.

12 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in 4.3. Use of field data efoe poceeding wit te sensitivity analysis, we ceck tat te feasiility and staility conditions of te equiliium points ae satisfied y te known paamete values, see Tale. Inpaticula,weemaktattedisease-feeequiliiumpointE 3 neve occus. amely, te feasiility condition () does not old in field conditions. Tis esult igligts te close connection etween Vaoa and viuses. Fute, it matces well wit eekeepe osevations: te ees vial infection is ound to occu weneve te mite population is pesent. Te same conclusion as een otained in ou pevious study,. Tale. Summay of te equiliia: feasiility and staility conditions Equiliium Feasiility Staility In field conditions E always (7) allowed E always unstale unstale E 3 () (8), (9), () infeasile E 4 () (), (), (3) allowed E always (4), () allowed E numeical numeical allowed simulations simulations 4.4. Sensitivity analysis In tis section we investigate te eavio of te system esponses wen te paametes cange tei values witin an appopiate ange. We compute te vaious populations equiliium values as function of a pai of paametes at a time, tus otaining sufaces in all te possile pais of paamete spaces. In ode to do tis, fo eac pai of paametes we comine te espective anges to uild a equispaced gid of values and ten we plot te fou sufaces esulting fom te values assumed y te populations at eac point of te gid. Since ou model contains paametes, tee ae 66 possile cases. Of tese, we pesent te most inteesting esults stating fom te est situation fo te ive, i.e. te cases in wic te mite-anddisease-fee equiliium E is staly attained. Tey aise fo te following paamete pais: γ, and m.

13 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in 3 In Figue 3 (left) te system sows diffeent tanscitical ifucations, as te paamete γ inceases. Specifically, fo smalle values of te oizontal tasmission ate γ te system settles to E, te ealty eeive scenaio, ten te infected ee population appeas in te system, i.e. we find te epidemics among te ees, E 4, and finally E, wee te ealty ees population disappeas and te mites invade te ive. An incease in instead positively affects te mite populations, as expected. Fom Figue 3 (igt), inceasing te ees disease-elated motality ate dives te system into a safe scenaio, all ees ae ealty and te colony ecomes mite-fee. A ige ees it ate seems to favo te infected mites and ees populations, wen te fome appea in te ecosystem and also te ealty ees, ut tis is tue only fo low values of. Instead te susceptile mites decease fo an inceasing. Hee tanscitical ifucations elate all te equiliium points. Te esults confim wat aleady emaked in. amely, te equiliium E is eaced fo vey small values of te tansmission ate γ comined wit a ig enoug ees disease-elated motality. Te eavio of te sensitivity suface in te m case is vey simila. In te paamete spaces λ and γ, not sown, a low γ o λ favos te ealty ees and ot ealty and infected mites, wile fostes all te system s populations. Fo β and δ te eavio is simila in tems of, ut deceasing te ote paamete depesses instead te infected mites. 4 I γ γ. I γ γ. 4-4 Figue 3. Sensitivity sufaces, in tems of te pai of paametes γ, (left), and (igt). Te ote paamete values ae coosen as =., =, m =.3, n =.7, e =., = 3, γ =., λ =.4, δ =.8, β =., =., p =.3. Initial conditions =, I =, = 6, = 4. In Figue 4 (left) as te Vaoa gowt ate inceases, te ealty ee population deceases wile te mite populations easonaly incease.

14 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in I 4 Te dynamics of te infected ee population, instead, is almost insensitive to te paamete, ut it gows linealy wit te daily ee it ate. Futemoe, tis kind of eavio is oseved weneve te paamete is consideed, egadless of te ote paamete eing taken into consideation. In fact, a lage means a geate nume of ees in te colony and tus popotionally also moe infected. As in te pevious Figue, infected mites enefit fom a ige, wile te ealty ones ae instead depessed. Anote tanscitical ifucation is sown in Figue 4 (igt) connecting temite-feeequiliiume 4 andcoexistence. TeiggeteVaoa gowt ate, te geate te gooming ate needed to wipe out te mites fom te eeive. If tis esistance mecanism is not ig enoug, te mites invade te colony. In te n paamete space, all te populations eave like in tis e case. 3 I e e e. - e. Figue 4. Sensitivity sufaces, in tems of te pai of paametes, (left), and e (igt). Te ote paamete values ae coosen as =., =, m =.3, n =.7, e =., = 3, γ =., λ =.4, δ =.8, β =., =., p =.3. Initial conditions =, I =, = 6, = 4. In Figue (left), te comined effect of te mites gowt and of te oizontal tansmission of te vius among ees ae epoted. Only fo eally small values of ot and γ te ealty ee population can suvive. An incease of ot γ and as a negative influence on te ealty ee population and a positive impact on te infected ee population. In tis egion of te paamete domain, te system appoaces te mite-fee attacto E 4. Fo sligtly lage values of anyone of tese paametes, te mites estalis temselves in te system. As a esult, we find coexistence followed, fo even lage values of ot suc paametes, y te ealtyees-fee equiliium E. Tis last tansition occus in paticula wen te ifucation paamete γ cosses te citical value γ.. ote tat

15 I Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in I γ γ γ γ Figue. Sensitivity sufaces, in tems of te pai of paametes γ, (left), and (igt). Te ote paamete values ae coosen as =., =, m =.3, n =.7, e =., = 3, γ =., λ =.4, δ =.8, β =., =., p =.3. Initial conditions =, I =, = 6, = 4. as inceases, te mites populations esonaly incease too, wile tey ae less sensitive to canges in γ. Te ees sufaces look also alike in te cases λ γ and δ β, not sown, ut te mites ave a peak at te oigin and ae instead depessed y an incease in eite one of te paametes. In te paamete spaces β, δ and λ te ee populations eave similaly, not sown. A simila pictue is also found in te δ paamete space. A low value of independently of te ote paamete leads to te ealty-ee-only point. Ten we find coexistence and ultimately tends towad te susceptile-ee-fee equiliium. Lage values of te paametes β and δ incease ot mites populations ut aove all te infected; instead fo te latte te opposite occus fo λ. Figue (igt) ette sows te influence of te disease-elated motality on te system dynamics. Stating fom low values of te ifucation paamete, we fist find te ealty-ees-fee equiliium E, ten coexistence followed y te ealty-ees-only equiliium E. Evidently, a ige value of te ees disease-elated motality as a positive impact on te ealty ees and a negative influence on te infected ee population. Te igest values of te paamete ( > ) depess ot mite populations. In paticula, we note tat as inceases, te populations I and damatically decease. Tus, a ige vial tite at te colony level is otained fo lowe values of. We tus find again te esult discussed in ou pevious investigation, : wen tansmitted y Vaoa mites, te least amful diseases fo te single ees ae te most viulent ones fo te wole colony. Te influence of a lage is positively felt just y te mite populations.

16 I Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in I e e. - e - e e. - e. - e - e Figue 6. Sensitivity sufaces, in tems of te pai of paametes e, (left), and e (igt). Te ote paamete values ae coosen as =., =, m =.3, n =.7, e =., = 3, γ =., λ =.4, δ =.8, β =., =., p =.3. Initial conditions =, I =, = 6, = 4. In Figue 6 (left) we exploe te effectiveness of te gooming eavio e as a esistance mecanism against to te Vaoa infestation. Te eavio of sensitivity sufaces is encouaging: egadless of te aveage nume of mites pe ee, if te gooming eavio is pefomed stongly enoug mites ae eventually wiped out of te system. Indeed, a tanscitical ifucation etween E 4 and te coexistence equiliium occus wen te gooming ate cosses te ifucation tesold e.7. Te ige te aveage nume of mites pe ee, te lage ot infected populations ecome, as well as te susceptile mites, wile ealty ees ae depessed. A simila eavio occus fo te case n. Also fo m n and m e we find te equiliium sufaces to look as fo te e case, wit te only diffeence tat in te last two cases te infected ees enefit fom a low m and at least in te ange exploed, tey do not vanis. Te ees in cases n and e eave similaly, see Figue 7 igt. ites and infected ees tend to disappea fo lage values of. Fo λ, γ, δ and β te infection in ot populations is indeed y lage values of te motality, wile ealty ees enefit fom it. Te emaining paamete plays a ole only on susceptile mites, wit low values fosteing tei gowt. Te sufaces in te paamete spaces n λ, n β, e λ, e δ, e β and n δ again look simila, see Figue 8 left. Te mite-fee equiliium is attained fo elatively lage values of n (o e), low values of te emaining paamete elp only te ealty mites. Cases λ δ, λ β ae again simila, wit lage λ favoing te infection in te ees and low values ofte ote paameteelping te ealty mites.

17 I Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in I I I γ e γ e - - e - e γ e -.. γ e - - e - e Figue 7. Sensitivity suface, in tems of te pai of paametes e γ, (left) and e (igt). Te ote paamete values ae coosen as =., =, m =.3, n =.7, e =., = 3, γ =., λ =.4, δ =.8, β =., =., p =.3. Initial conditions =, I =, = 6, = 4. Again in te β γ and δ γ cases te most influencial paamete is γ, lage values leading to te ealty-ee-fee point. Te ote paamete elps te susceptile populations if it is small. Cases p, p and p m ae simila, a lage p depesses all te populations ut te ealty ees, low values of te ote paamete elping tem too. Also in te paamete spaces p δ and p β te sufaces ave simila sapes. Fo p λ and p γ we find a simila eavio in tems of p, ut ee te infected ees ae muc less affected y its gowt. Fo p simila consideations old, ut again we find tat a decease in educes te infected ees, wile p on tem as scant effect δ n δ n β δ β δ δ n δ n. -4 β δ -4-4 β δ -4 Figue 8. Sensitivity suface, in tems of te pai of paametes n δ, (left) and δ β (igt). Te ote paamete values ae coosen as =., =, m =.3, n =.7, e =., = 3, γ =., λ =.4, δ =.8, β =., =., p =.3. Initial conditions =, I =, = 6, = 4.

18 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in I 8 I λ m λ m n p n p λ m.. λ m... n. p... n. p. Figue 9. Sensitivity suface, in tems of te pai of paametes m λ, (left) and p n (igt). Te ote paamete values ae coosen as =., =, m =.3, n =.7, e =., = 3, γ =., λ =.4, δ =.8, β =., =., p =.3. Initial conditions =, I =, = 6, = 4. InFigue8igt,wecanoseveteinfluenceoftetansmissionatesδ and β. As tei values incease, ot te ealty ees and mites populations decease wile te infected ones incease. As a esult, te level of infection in te ive damatically gows and almost all te Vaoa mites ecome vius caies. In Figue 9 left, anote tanscitical ifucation occus. An incease of te paamete λ as its geatest influence on te ealty ee population and makes it completely vanis fom te system. Hee te tansition etween te coexistence and te no-ealty-ees equiliium point is oseved. Convesely, te tansmission ate etween infected Vaoa and ealty ees as a positive effect on te infected ees population. Futemoe, also te mites populations decease ut evidently tis is due to te eduction of te ost population. Finally, fom Figue 9 igt, we note te comined effect of two paametes tat inde te infestation speading. If te Vaoa natual motality is ige tan te tesold value n.3 te mites disappea fom te ive and we find equiliium E 4. Instead, te influence of Vaoa intaspecific competition is sown fo smalle values of n: as p inceases, ot te mites populations and te infected ees decease wile te ealty ees incease. Tis eavio occus also in te p e paamete space. Altoug a sufficient ate of te gooming eavio can eally elp te colony to contol te Vaoa infestation, Figue 6 (igt) suggests tat te disease-elated motality as a moe significant influence on te sape of te sufaces and tus on te system dynamics tan te gooming ate e. Indeed te sape of te sensitivity sufaces is mostly detemined y te

19 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in 9 canging paamete. ote ee te following tanscitical ifucations: statingfolowvalues ofwefind E, te ealty-ee-feepoint, ten, as inceases, te coexistence equiliium and finally te mite-fee equiliium E 4. Tis cain of tansitions eempasizes once again te effect of te paamete. Futemoe, fom Figue 7 left, te effectiveness of te gooming eavio can e elevant only in te pesence of a eally low value of te paamete γ, i.e. in tis egion wee te tanscitical ifucation mentioned aove etween coexistence and E 4 is oseved, wit e as ifucation paamete. Convesely, fo ige values of te oizontal tansmission ate among ees, te colony is diven fom E 4, te mite-fee situation, toug coexistence, to E wit te extinction of ealty ees. Te knowledge of tese most elevant paametes could suggest te ones wic would maye e affected y uman extenal measues in ode to dive te system to settle in a possily safe position. To sum up, te paametes most affecting te system ae,, γ and e. Ou findings elucidate tei influence on te system. Tese ae oweve all ecosystem-elated paametes, tat peaps can adly e influenced y man-undetaken measues, altoug tey migt depend on ote extenal factos suc as fo instance climatic canges. Teoetically, te sensitivity sufaces sow tat te population of ealty ees would igly enefit fom a eduction of te Vaoa gowt ate, as well as a eduction of te oizontal tansmission ate among ees. Also a ige ees disease-elated motality and a lage gooming ate would elp in potecting te colonies.. Conclusion We ave intoduced a model fo ees and mites, modifying ou pevious appoac, allowing a link etween te cayingcapacity of te mites and te ees population, wile in te fome investigation tese wee not ound togete. Ou esults agee wit empiical evidence, in addition to desciing te fundamental ole played y te mite in tis pocess. Te ecosystem descied ee can settle only to te following possile outcomes: te diseaseand-mite-fee envionment, te ideal situation; te mite-fee situation, in wic oweve pat of te ees ae endemically infected; te coexistence equiliium in wic ot ees and mites ae pesent, all affected y te disease; te Vaoa invasion leading to extinction of te ealty ees, wile te emaining ones ae all infected. Among ou findings, we oseve tat te endemic disease cannot affect

20 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in all te ees in a Vaoa-fee colony. Indeed, te infected-ees-only situation, descied y equiliium E in ou analysis, tuns out to e unstale, if te mite epoduction ate exceeds tei motality ate, i.e. te opposite of condition (3) olds. Convesely, te Vaoa invasion scenaio can ecome possile wee te ees ae all infected and mites invade te ive, te point E aove. In ageement wit te fact tat te pesence of V. destucto incements te vial tansmission, tis esult indicates tat te wole ee population may ecome infected wen te disease vecto is pesent in te eeive. Also, te two ealty populations cannot suvive togete, in te asence of infection. Tis esult agees wit field osevations, in wic te ee colony is consideed infected wen V. destucto is pesent in te eeive. Tus, te discovey of Vaoa in te ive necessaily implies tat at least pat of te ees ae vius-affected. Te findings of tis study also indicate tat a low oizontal vius tansmission ate among oney ees in eeives will elp in potecting te ee colonies fom te Vaoa infestation and te vial epidemics. In fact, fom te fist condition in (7), a low γ is necessayfo te disease-and mite-fee equiliium E to e stale. Te sensitivity analysis allows us to identify te paametes most affecting te system:,, γ and e. Specifically, fo small canges of tese paametes te system expeiences tanscitical ifucations etween all te equiliia. Tese esults empasize te impotance of keeping te Vaoa gowt unde contol. In ode to do tis, te sensitivity sufaces suggest tat a decease of γ comined wit an incease of e will enefit te colonies. Futemoe, te analysis sustantiates te empiical emak tat te most amful diseases at te colony level ae tose tat least affect te single ees. In fact, we can oseve tat fo a vey low ee motality, te equiliium towad wic te system always settles is te ealty-ee-fee point, in wic only infected ees and mites tive. Te counteintuitive effect of te paamete on te system s dynamics agees wit te expeimental findings of 7. y atificially infecting te laval colonies of Apis mellifea and Apis ceana, te eseaces ave compaed te Vaoa evolution in ot colonies. Tey found tat in te euopean ees it does not cange muc in compaison wit te colonies tat ae not infected, wile te oiental ees wee indeed, wit a ig laval motality. Tey conclude tat te ige vulneaility of te latte to te mite could imply a ige esistance at te colony level, ecause te infected lavae ae moe easily spotted and killed y te woke ees. Tis also agees wit te

21 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in geneal consensus tat te Vaoa epesents te most impotant facto affecting te suvival of Apis mellifea colonies, wile te oiental ees tive elatively easily in its pesence. In a simila way, a ige ees motality induced y te diseases vectoed y Vaoa maintains te total vial load in te colonies at a elatively modeate level. Tis tus enances tei suvival cances, ecause te individual infected ees ave less time availale fo oizontally tansmitting te vius teey infecting ote individuals. ot 7 and tis investigation sow tat vulneale individuals may elp te supeoganism, against te common assumption tat it is just te pesence of stong individuals to ensue te colony suvival. Compaing te esults of tis investigation wit te outcomes of te model, it appeas tat te pesent model, moe matematically elaoated in tat it contains a tem of Leslie-Gowe type, namely te mite ee-dependent caying capacity, is not eally fundametally necessay, as its qualitative analysis is captued well y te fome system. Te ees ae te most impotant pollinatos woldwide. In ecent yeas managed oney eeives ave een suject to decline mainly due to infections caused y te invading paasite Vaoa destucto, leading to widspead colony collapse,. In tun, tis as caused alam among te scientists fo apicultue and even moe fo agicultue. Tis model epesents a fute step eyond te asic model pesented in fo te undestanding of tis polem. ut tee is te need of te model validation using eal data. At pesent, we ae cuently woking on tis aspect of te eseac, gateing te field data. Fom te collected field infomation we plan to otain an estimate of te model paametes of most inteest fo te eekeepes: fist of all te Vaoa gowt ate, a paamete most needed to assesste amount of acaicides to e used and te timings fo tei most effective application. In ode to measue eal field data we ave stated field expeiments, in collaoation wit eekeepes of te Coopeative Dalla Stessa Pate, Tuin, and fom Aspomiele, te Association of Piedmont Honey Poduces, Tuin. Tese expeiments epesent poaly te fist collaoation etween eekeepes and matematicians on tis topic in ou county and tey ae eally pecious fo te eseac, since eal field data ae not yet pesently availale. Refeences.. A. enavente, R. R. Deza,. Eguaas, Assessment of Stategies fo te Contol of te Vaoa destucto mite in Apis mellifea colonies, toug a

22 Feuay 7, 8 :4 Poceedings Tim Size: 9in x 6in simple model, ACI, II Congeso de matemática aplicada, computacional e industial, Editos E.. ancinelli, E. A. Santillán-acus, D. A. Tazia, Rosaio, Agentina, Deceme 4t-6t 9, -8.. S. enadi, E. Ventuino Vial epidemiology of te adult Apis ellifea infested y te Vaoa destucto mite, HELIYO (6). 3. E. Capana,. Lodesani, Editos. Patologia e avvesità dell alveae (Patologies and advesity of eeives), spinge, J. S. Figueió, F. C. Coelo, Te ole of esistance eavios in te population dynamics of Honey ees infested y Vaoa Destucto, Astact Collection, odels in Population Dynamics end Ecology, Ezio Ventuino (Edito), Intenational Confeence, Univesità di Toino, Italy, August t-9t 4, 3.. E. Genesc,. Auet, Emeging and e-emeging viuses of te oney ee (Apis mellifea L.), Veteinay Reseac 4, (), ilito, iologia della api, IZSLT ttp:// 7. P. Page, Z. Lin,. uawangpong, H. Zeng, F. Hu, P. eumann, P. Cantawannakul, V. Dietemann, Social apoptosis in oney ee supeoganisms, Scientific Repots 6, Aticle nume: 7 (6) doi:.38/sep7 8. V. Ratti, P. G. Kevan, H. J. Eel, A matematical model fo population dynamics in oneyee colonies infested wit Vaoa destucto and te Acute ee Paalysis Vius, Can. Appl. at. Q. (3), no., DOI.7/s V. Ratti, P. G. Kevan, H. J. Eel, A atematical odel of te Honeyee Vaoa destucto Acute ee Paalysis Vius System wit Seasonal Effects, ull at iol DOI.7/s J.F. Santos, F.C. Coelo, P.J. liman, eavioal modulation of te coexistence etween Apis mellifea and Vaoa destucto: A defense against colony collapse? PeeJ PePints 3:e739, (). doi.og/.787/peej.pepints.396v. United ations Envionment Pogamme, Gloal oney ee colony disodes and ote teats to insect pollinatos, UEP Emeging Issues, (UEP, aioi) ().