Banana Xanthomonas Wilt infection: the role of debudding and roguing as control options within a mixed cultivar plantation

. Anoptimalcontrolframeworkisdesignedinwhichtheuseofcleanplantingmaterials,debudding,disinfectionoftools,androguing areconsideredascontrolmeasuresofBanana Xanthomonas Wilt (BXW) within a plantation of multiple cultivars. A model for a specialcaseoftwocultivars(AAA-andABB-genomecultivars)wasanalyzed.ByPontryagin’sMaximumPrinciple,wecharacterized anddiscussedpossiblecontrolstrategiesthatsubstantiallyreducetheinfectionlevelsofBXWwithinaplantationofABB-and AAA-genomecultivars.Acombinationofbothpreventionandcontainmentcontrolsyieldedthegreatestdeclineintheinfection levelsinbothcultivars.Additionally,foreffectiveBXWmanagement,itisimportanttoassesstheendemicleveloftheplantation beforeapplicationofcontrols,andonceimplemented,thisshouldbemaintainedevenwhenthediseaseisundetectabletoeliminate possibleresurgence.


Introduction
Banana is a major crop grown in the East and Central African region.It provides food security and income to over 20 million people in the region [1].Since 2001, Banana Xanthomonas Wilt (BXW) has been reported as a major threat to the banana production and livelihoods of millions of people have been affected [2].The pathogen caused by Xanthomonas campestris pv.musacearum (Xcm) was first reported in Ethiopia on Ensete, a related crop to banana in 1964 [3].The disease is transmitted by insect vectors, birds, bats, contaminated farming tools, and infected suckers used for setting up new plantings.It attacks all the commonly grown cultivars although some are more susceptible than others.Disease symptoms include yellowing and wilting of leaves, male bud wilting, premature ripening and rotting of fruit, yellow ooze observed on the cross section cut of the pseudostem, and eventually death of the entire plant.
Since the reporting of BXW disease, scientists have disseminated information on identification based on symptoms, mode of spread, and how to implement cultural control practices that have been used to manage similar diseases such as Moko and Bugtok [4,5].Cultural control practices for management of BXW have been categorized as preventative or containment.Among the preventative practices, we have debudding which is the removal of the male bud by twisting with a forked stick as soon as the last cluster is formed [6].This approach avoids contact between the tool and the potentially infected tissues.This is done within 3 weeks since flowering to prevent insect vector from transmitting the disease.
Use of clean planting materials is another preventative approach being advocated for although there are no proper screening facilities for the seed system.Farmers basically obtain suckers for new plantings either from their plantations or from neighboring plantations.Such suckers could 2 International Journal of Mathematics and Mathematical Sciences be latently infected but since the available detecting tool (lateral flow device) is not yet accessible to the farmers, there is a risk of transferring latently infected plants to new plantations.Nevertheless, to reduce this risk, farmers are constantly cautioned on the use of these suckers and the need to monitor them as they develop.The third preventative approach is the disinfection of tools between plants using sodium hypochlorite solution (usually known by the trade name of the commonest brand in Uganda: Jik) or flames of fire [5,7].
Two approaches have been implemented as containment practices.These are single disease stem removal (SDSR) dealing with removal of the infected stem from the mat at stool level and roguing involving removal of the entire mat from which the infected stem arose.The latter is more effective as it eliminates the risk of vertical transmission from mother stem to attached suckers.However, it is expensive to implement as it requires more labor [7].Kubiriba et al. [8] evaluated the SDSR approach and noted that it was effective especially if the manifesting symptoms are via the inflorescence part of the plant.
Jogo et al. [9] studied the extent to which cultural control practices have been adopted by farmers and assessed the socioeconomic factors that influence adoption of these practices among smallholder farmers in Uganda.Their results indicate that among the key determinants were labor and perceived effectiveness of the technologies being implemented.This indicates knowledge gaps in identifying the combination of control options which is cheap enough to implement but causes a large decline in number of infected stems.
The best combination of cultural control practices also depends on the main production systems within the region.The two main groups of cultivars include the East African highland banana (EAHB) AAA-genome grown in western Uganda (1700 masl) and the "kayinja" (beer banana) ABBgenome cultivar grown in the central Uganda (1300 masl) [9].The intensity of management among these cultivars is different with the former being managed intensively than the later [9].Also ABB-genome cultivars are more susceptible to insect infection due to the greater production of sweet nectar that attracts the insects as compared to AAA-genome cultivars.These features influence BXW management and control adoption by farmers in these respective regions.
Even though application of all control approaches has been advocated, few farmers are in position to implement them.Therefore, it is necessary to examine the approaches and identify the most cost-effective combination for management of BXW disease.
The use of optimal control approaches/models allows the determination of the most cost-effective intervention.Optimal control modeling has been extensively used to understand biomedical problems and suggest optimal control combinations that will minimize the costs [10][11][12][13][14].Our model focuses on determining the optimal control combination within smallholder plantations with mixed cultivars given that in banana cropping system different cultivars are managed differently and the risk of infection is also different.
Thus, we propose a Healthy-Infected model for the spread of BXW via inflorescence infection, vertical transmission, and tool-mediated spread.We formulate an optimal control model with debudding, use of clean planting material (reducing the proportion of infected suckers), roguing/SSDR, and disinfection of tools used for cultivating between and within cultivars as our control options.The optimal control model is analyzed using the Pontryagin's Maximum Principle [15][16][17] and numerical simulations assessing different combinations of control options.The paper is organized as follows: in Section 2, the formulation of the optimal control model for a plantation with multiple cultivars () is presented.Mathematical analysis and numerical simulations of the optimal control model considering a special case of two cultivars are given in Section 3 and lastly conclusions drawn are given in Section 4.

Model Formulation
Consider a population of banana stems with  different cultivars each subdivided into healthy   and infected   classes for  = 1, 2, 3, . . ., .We assume negligible latent and incubation period for simplicity of the model, also because omission of this compartment does not entirely affect the dynamics of the disease.It is assumed that in a unified plantation several blocks of each cultivar are maintained.The effort going into maintenance of each block will depend on preference of farmers in a specific region.The healthy compartment contains all banana stems that are diseasefree but at a risk of being infected either by contaminated tools or by insect vectors during inflorescence formation.The populations per cultivar are maintained at equilibrium by a logistic function at a constant rate  and a constant carrying capacity .The equilibrium  is a result of various activities within the plantation including replanting, desuckering, sucker emergence, and harvesting.The infected class is increased by the following scenarios.
(1) Let   denote the proportion of healthy suckers from infected mats; then a proportion (1 −   ) indicates the suckers recruited into the infected compartment arising from vertical transmission.
(2) Disease transmission by contaminated tools is modeled using mass action incidence term with   as the coefficient of effective contact between healthy and infected stems by tools.Assuming random use of tools, spread is independent of plantation density but rather dependent on the proportion of infected stems within the plantation.Therefore, the incidence term due to contaminated tools is given as       for the interaction between cultivars  and .
(3) Transmission by vectors via the inflorescence infection is considered to be frequency dependent, that is, ( V   ,   )/ with  V as the coefficient of effective contact by vectors between cultivars  and .This incidence term is reasonable because spread by vectors depends on the number of flowering stems and vector activity resulting in infection at a particular time.
Once infected, infected stems showing symptoms are removed at a constant rate   for cultivar .This rate also depends on the available resources and farmers' desire to control that particular cultivar.In some regions within the central Uganda, farmers tend to concentrate on the cooking banana ignoring other cultivars like kayinja which tends to increase the risk of vector spread.We now incorporate the main cultural control practices that are being implemented in the management of BXW disease into the model formulation.
(I)  1 () as the control is associated with ensuring that only clean planting materials are used for replanting after mat removal or infilling of space.This involves the ability to detect and remove infected suckers from mats such that they are not used for replanting.Since there is no proper screening for the seed system, farmers rely on suckers from their plantations or neighboring plantations for replanting.It is necessary to determine whether a sucker/planting material is clean.The lateral flow device developed could be used to detect even latently infected material but currently it is not available for farmer and even then it would be expensive.Combining all the cases above, the following system can be used to describe the dynamics of BXW when control measures are being implemented.
The term (1 −  1 )    represents the proportion of healthy suckers that are recruited into healthy population arising from infected suckers given application of control.With  1 () = 1, only healthy plants are maintained in the plantation.Thus,  1 ,  2 ,  3 ,   = 1 indicated full effort placed in the implementing of controls while  1 ,  2 ,  3 ,   = 0 indicated failure to implement the controls.The objective function to be minimized is given as with   being the maximum time for which the control practices are implemented.
The constants   are the weight constants on the infected plants of cultivar  while  1  1 (),  2  2 (), (1/2) 3  2 3 (), and (1/2)   2  () are the corresponding costs for controls  1 ,  2 ,  3 , and   , respectively, in terms of labor or cash.The costs associated with use of clean planting materials ( 1 ) and debudding ( 2 ) are directly proportional to the rates at which these controls are implemented.Thus, it is reasonable to consider the objective function with linear controls  1 ,  2 strategy.Later these will be adjusted to include quadratic control terms to make the problem mathematically more tractable.The quadratic cost functions are used as the simplest form of describing nonlinear costs involved in implementation of the controls disinfection of tools (  ) and roguing  3 [10].

International Journal of Mathematics and Mathematical Sciences
Applying all the control options would be the most rapid method in stopping the spread of BXW but this is often too costly [9].Therefore, here we determine combinations of control measures that would effectively lead to disease eradication at minimum total cost of implementation.

Optimal Control Model: Special Case of Two Cultivars
To analyze the optimal control model, we consider a special case whereby, on a particular plantation, only two cultivars, that is, AAA-genome cultivars ( = 1) and ABB-genome cultivar ( = 2), are considered.This is justified by the common practice in the affected areas of growing these two types, one for food and trade and one for beer.First, we formulate our optimal control problem considering linear cost function for use of clean planting material and debudding controls with an intention of determining subject to From system (6), it is noted that, in the absence of disease, no control practice is implemented (  = 0,  1 = 0,  2 = 0,  3 = 0,   = 0).The population growth per cultivar is determined by Moreover, by solving for   (), we obtain Thus, given   (0) ≥ 0,   (0) ≥ 0, for  = 1, 2, the state solutions   (),   () are positively invariant for all  ≥ 0. That is, the set is positively invariant.Furthermore, the following holds: (), respectively.We show that the optimal control solution exists and then we characterize the system using Pontryagin's Maximum Principle [17,20] to obtain the optimal control solution.
Theorem 1.There exists an optimal control solution for objective function ( 5) subject to model system (6).
Proof.Let () denote the state variables and () denote the control variables.Let the integrand for the objective function be defined as To prove existence of the optimal control solution, we apply Theorem 4.1 in [21] by checking the following assumptions.
(A1) The set of control and corresponding variables is nonempty.
(A2) The control set Ω is closed and convex.
(A3) The right-hand side of the state system is bounded above by a sum of bounded control and the state and can be written as a linear function with coefficients dependent on time and state.
(A4) The integrand is a convex function on Ω.
(A5) There exists constant  1 ,  2 > 0 and  > 1 such that the integrand satisfies Since the equations of system (6) have bounded coefficients and the solution exists in finite time intervals, Theorem 9.2.1 in [22] guarantees the existence of the solution to the model system which gives assumption (A1).By definition of the set Ω, (A2) is also satisfied.From system (6), the righthand side of the equations is continuous and linear in the controls indicating that it can be written as (, , ) = (, ) + ℎ(, ), where  = ( 1 ,  For (A4), the integrand (, (), ()) is convex; that is, it can be easily shown given  ∈ (0, 1).Let () and ]() be two vectors in the integrand; then Lastly there exist constants  1 , Thus, assumption (A5) is satisfied.

Characterization of the Optimal Control Solution.
Pontryagin's Maximum Principle is applied to obtain the necessary conditions that an optimal control pair ( * (),  * ()) must satisfy.The principle converts the optimal control problem given by objective function with the transversality condition (  ) = 0.
We define the Hamiltonian equation ) The differential equations governing the adjoint variables are obtained by differentiating the Hamiltonian  with respect to the state variables.The adjoint system is obtained by determining written as where with the transversality condition  1 (  ) =  2 (  ) =  3 (  ) =  4 (  ) = 0.For the optimal controls, we solve for ) by differentiating  with respect to the respective controls, and considering the bounds, we obtain the optimal control solution as follows.
For the linear controls  1 ,  2 , bounded such that 0 ≤ () ≤ 1, the optimal control solution is given as Let  1 = / 1 and  2 = / 2 referred to as switching functions for the controls  1 and  2 , respectively.In the region where  1 and  2 are not equal to zero, we say that the controls  1 and  2 are bang-bang controls.In this case, applications of these controls are switched between the lower and the upper bound through the period of implementation.To address the case where  1 and  2 are zeros for some time interval { 1 ,  2 }, we note that all the derivatives with respect to the corresponding controls vanish in this time interval.The functions  1 and  2 are differentiated with respect to time and substitution is consequently done until the controls  1 and  2 , respectively, reappear.The optimal controls are then referred to as singular arcs.The singular arcs  1 and  2 are then determined by equating the time derivatives to zero and solving.If the th derivative of the switching functions  1 ,  2 is nonzero and the Generalized Legendre Clebsch Condition given by is satisfied, then the singular controls have order  [23].Thus, for the system (Section 3), we have the optimal control solution given as International Journal of Mathematics and Mathematical Sciences 7 For the quadratic controls, given the bound 0 ≤  3 ,   ≤ 1, where ,  = 1, 2, the optimal controls are given as These can be summarized as where  * 1 ,  * 2 ,  * 3 ,  * 4 are solutions to the adjoint system (18) with the transversality condition.The optimal controls for use of clean planting materials ( 1 ) and debudding ( 2 ) tend to suggest existence of singular arcs.In this case, the controls are implemented whenever the resources are available which better suits application than the bang-bang case where they are just switched between the lower and the upper limit.According to Ledzewicz and Schätler [23], singular controls (if they exist) tend to be either the best (minimizing) or the worst (maximizing) strategies and in either case they are essential in determining the structure of optimal controls.However due to the complexity of the optimal control model, this can not be verified.Thus we consider quadratic terms for all the controls in the objective function for the following reasons: (1) Minimizing linear control terms in the objective function is similar to minimizing quadratic terms.( 2) With varying seasons, the cost associated with implementation of these control may vary unproportionally to the rate of implementation.

Quadratic Control Terms for Debudding and Use of
Clean Planting Materials.Considering quadratic terms for the objective function to be minimized and by applying Pontryagin's Maximum Principle, we have the following Hamiltonian function: The existence of controls in this case can be proved using similar techniques as in Theorem 1.The adjoint variables are International Journal of Mathematics and Mathematical Sciences the same as system (18).Thus the optimal control solution can then be given as }} . ( The optimal control characterizations are inversely proportional to the associated weighted costs of implementation indicating that the effectiveness of the control strategy greatly depends on the necessary costs. Next we solve numerically the optimality system consisting of state system (6) with corresponding initial conditions, adjoint system (18) with the transversality condition, and optimal control characterization assuming the costs terms associated with the controls are quadratic (26).

Numerical Simulation. In this section, numerical results
for the optimal control model given in Section 3 are studied.The backward-forward sweep method with the fourth-order Runge-Kutta algorithm is applied to solve the optimality system.The algorithm has been implemented by various authors to obtain optimal control solutions of related systems [11,14,17].For easy reference we summarize it below.
(i) Make an initial guess of the control;  = 0 is always sufficient.(ii) Using the initial conditions  10 ,  10 ,  20 ,  20 , and the values of  solve for  forward in time according to the differential equations in the optimality system given.
(iii) Using the transversality condition (  ) = 0 and the values of  and  solve the adjoint system  backward in time according to the corresponding differential equations in the optimality system.(iv) Update the control set  by entering the new values of  and  into the characterization of the optimal control.
(v) Verify for convergence the following: if the values obtained are sufficiently close to the corresponding ones in the previous iterations, then the output of the current values is the optimal control solution; otherwise return to step (ii).
This algorithm was implemented in Matlab by modifying the optimal control code presented by [17].
The numerical values used for the numerical simulation as reflected in Table 1 are obtained from literature and some are selected from given ranges based on the banana growing system in Uganda, for which the feasible solution is guaranteed.
According to the guide by [24], the maintenance rate, , is considered to range between 3 between 5 stems of varying ages per mat.It is noted that the smaller the number of stems the bigger the bunch.Also, for a new plantation, the carrying capacity  of stems per hectare given 3 m × 3 m spacing is assumed to be 750-900 mats [24].This means, for a wellestablished plantation, the average carrying capacity would range from 2250 to 4500 stems.Roguing by removal of only the infected stem at stool level or by removal of the entire mat depending on the severity of the disease is considered to range between 0 and 1 per month per hectare per infected stem.The proportion of healthy suckers arising from infected mats is also assumed to range between 0 and 1. Transmission by tools is dependent on the activities being carried out within the plantation.For instance, harvesting of the leaf might require a single cut as compared to desuckering which requires several cuts and exposes the attached mat to infection.Accordingly, we assume the rate of tool-mediated transmission to range between 0 and 0.002.The given range guarantees existence of the optimal solution if other factors are kept constant.The rest of the parameters are as given in [18].The choice of the weight constants is intended to reflect the cash and labor involved in the control measure.According to the study by Sebikari [19], use of clean planting materials and roguing was more expensive than debudding and disinfection of tools.In Table 2, the costs for implementation of these controls are given.
The low cost for disinfection of tools was due to the fact that most farmers use fire frames made from collected firewood instead of JIK.Additionally, the required labor for activities like disinfection of tools and debudding is mainly from family members which makes these activities less expensive.Thus, we consider reasonable ratios as weights such that  1 = 50,  2 = 50,  1 = 1212,  2 = 6, 3 = 1346, and  11 = 12 =  21 =  22 = 1.Assume the control options are implemented for a period of two years such that   = 24 (months).
To illustrate the dynamics of disease with and without controls, we consider the following initial conditions:  10 = 300,  10 = 50,  20 = 300,  20 = 50 chosen arbitrary given that the model exhibits global stability behavior (not presented here).The controls are then introduced after one month.We investigate the following combinations of cultural control practices.Note that several combinations of control strategies have been considered.However, only those giving substantial reductions when applied as compared to no control are presented.In Figures 1 and 2, we present how the incidence of BXW varies from the same starting position with each control strategy.The objective function was optimized using only controls in the given strategy while setting others to zero. .Parameters are as given in Table 1.

International Journal of Mathematics and Mathematical Sciences
Strategies (A) (debudding only), (B1) (debudding and disinfection of tools between cultivars), and (B2) (debudding and disinfection of tools within cultivars) are all preventative control practices.In both Figures 1 and 2, the population levels of AAA and ABB infected stems obtained using these control are lower than when no control is applied.Furthermore, we observe that implementation of these controls only reduces the risk of further infection but does not lead to reduction in already infected stems.For example, the number of infected AAA and ABB cultivars was lower by 10% of the expected number in the absence of control after 20 months of application.It is also noted that disinfection of tools ( 11 ,  12 ,  21 ,  22 ) led to relatively better results in the AAAgenome cultivar as compared to ABB-genome.This indicates that, in the absence of containment controls, disinfection of tools within AAA-genome cultivars is crucial.From these results, we can conclude that application of only preventative controls delays further infection but, with these parameters, if no other measures are implemented, then uncontrolled endemic levels may be attained asymptotically.
When either of the strategies (B3) (debudding and roguing), (C1) (debudding, roguing, and disinfection of tools between cultivars), and (C2) (a combination of debudding ( 2 ), roguing ( 3 ), and disinfection of tools within cultivars) is implemented independently, the number of infected AAA and ABB stems reduced greatly as compared to the counterparts where only prevention controls are applied.In this case, both prevention and containment controls are applied and as expected the number of infected stems declines by more than 70% asymptotically to zero.Furthermore, we notice the importance of disinfection of tools among the ABB-genome cultivars.The findings here indicate that disinfection of tools within the highly susceptible cultivars (ABB) and between cultivars is important and reduces the risk of further infection that could arise during roguing and pruning exercises.
Lastly, when all controls are implemented, we observe that the number of infected stems in both cultivars declined rapidly (Figures 1 and 2).That is, AAA and ABB infected stems declined by 90% of the expected number of stems considered in the absence of the controls.
In Figure 3, we present the control profiles for each strategy defined above.In all cases, we observe that the control profiles for debudding ( 2 ) and roguing ( 3 ) are maintained at maximum intensity throughout the period of implementation.The control of use of clean planting materials ( 1 ) is maintained at the lower limit indicating that this control is neutral and should not be implemented alone but in the combination with other controls.As noted in Figures 1 and 2, application of all controls gave better results as compared to the strategies when  1 was set to zero.The control profile for disinfection of tools between and within cultivars ( 12 ) and ( 22 ) in strategies (B1) and (B2) is maintained at the upper bound throughout the period of implementation while those in strategies (C1) and (C2) decline and oscillate between the upper and the lower bound after a period of 15 months.This indicates that, in the absence of roguing and considering debudding and disinfection of tools as the only controls applied, maximum effort should be used to reduce new infection.On the contrary, with roguing disinfection of tools can be at maximum intensity when incidence is still high but relaxed and applied at intervals as the infection levels decline.
In Figure 4, numerical results compare control strategies applied at different infection levels ( = 0,  = 1 month,  = 2 months) with the case of no control.We notice that if control measures are applied at  = 0 when the prevalence is about 10%, then it would require about 20 months to bring infection below 1% and 30 months for less than 0.1%.On the contrary, if controls are introduced when the prevalence is about 80% ( = 1), then much more effort and time will be required to attain a prevalence of less than 1%.Also if controls were introduced when the disease was already endemic, then controlling the diseases would become much harder as it would require more time to reduce infection to undetectable levels.Clearly, we notice that disease management and containment are much dependent on the infection level.This implies that once the disease is endemic it may be more cost-effective in clearing the plantation and replanting than in continuously applying the controls which might take up to 5 years to bring disease to undetectable levels.
Additionally, we observe that even when infection is reduced to undetectable levels, that is, <1%, irrespective of the time of withdrawal of controls, the prevalence increases rapidly towards the endemic equilibrium.In Figure 5, a comparison of the cases where controls are completely withdrawn and when few are maintained is presented after infection is <1%.We observe that the number of infected stems continues to decline for the case where controls are just reduced compared with the case when they are completely withdrawn.This implies continuous monitoring and application of controls when necessary should be emphasized even when infection levels are undetectable to eliminate cases of resurgence within or outside plantations.B1), (B2), (B3), (C1), and (C2), respectively, when quadratic cost terms are considered.The graphs are maintained at maximum value in all strategies for debudding and roguing.Parameters are as given in Table 1.

Conclusions
In this paper, a deterministic model for optimal control of BXW within a mixed cultivar plantation is derived and analyzed.The possible modes of spread are incorporated including inflorescence infection and use of infected planting materials obtained via mother-to-sucker and tool-mediated spread were considered.The model focuses on the interplay of control options such as debudding, use of clean planting material, disinfection of tools, and roguing as control measures of BXW.
At first, linear cost terms for the debudding and use of clean planting material controls were considered.This is because the cost of implementation of these controls is directly proportional to the rate of implementation.For the other controls, quadratic cost terms were considered due to the nonlinear changes in their implementations.By applying Pontryagin's Maximum Principle, the necessary conditions for optimal control of BXW were obtained.These conditions suggest possibility of having debudding and use of clean planting materials controls as bang-bang controls within the upper and the lower bounds.They appear also to have singular arcs on a given interval for which controls  1 and  2 are reactivated whenever necessary although this could not be explicitly verified due to the complexity of the model.

(
A) Optimal debudding (B) Combination of two controls (B1) Optimal debudding and disinfection of tools between cultivars (B2) Optimal debudding and disinfection of tools within cultivars (B3) Optimal debudding and roguing

Table 1 :
Values used for model simulations.

Table 2 :
Table showing costs of implementations of control options [19].