Volume 4, Issue 6, November 2016, Pages: 56-70

Modelling of Infection Mildew of Taro (Phytophthora colocasiae)

Djouokep Léonel Gautier1, Asseng Charles Carnot1, *, Bowong Tsakou Samuel1, Ambang Zachée2, Monkam Tchamaha Fabrice1

1Faculty of Science, Laboratory of Plant Biology and Physiology, University of Douala (Cameroon), Douala, Cameroon

2Faculty of Science, Department of Plant Biology, Laboratory of Phytopathology and Microbiology, University of Yaounde, Yaoundé, Cameroon

Email address:

(A. C. Carnot)
(A. C. Carnot)

*Corresponding author

To cite this article:

Djouokep Léonel Gautier, Asseng Charles Carnot, Bowong Tsakou Samuel, Ambang Zachée, Monkam Tchamaha Fabrice. Modelling of Infection Mildew of Taro (Phytophthora colocasiae). Plant. Vol. 4, No. 6, 2016, pp. 56-70. doi: 10.11648/j.plant.20160406.13

Received: September 21, 2016; Accepted: October 7, 2016; Published:October 31, 2016

Abstract: Mildew taro caused by Phytophthora colocasiae affection is the most devastating of taro cultivation in Cameroon since 2010. It has been studied in leading the influence that can have a parameter considered favourable in the kinetics of the disease, and secondly, the interaction between plots through zoospores that can move from one field to another while estimating their dispersal throughout the plant. These models have allowed us to demonstrate that the duration of pathogen latency period, the number of sporangia produced on the surface of a lesion as well as the severity of the infection taken individually, are parameters to be taken into account in the development of a variety resistant to late blight taro. The dynamics of the fungus over time is represented by a matrix. The latter was used to establish a detailed estimate of the number of new infections caused by a sporangium placed in a landscape of healthy leaves. This number is known as the net rate of breeding base name (R0). The incidence and severity of disease are significantly reduced when the rate is less than or equal to one. So our approach can be used to guide research programs or evaluate the effectiveness of control strategies to design throughout the plant.

Keywords: Phytophthora colocasiae, Taro, Modelling, Simulation

1. Introduction

Taro (Colocasia esculenta (L.) Schott) is a monocotyledonous plant belonging to the Araceae family [1]. This plant about two meters high is from Asia. Originally from Malaysia, its culture developed to Polynesia where does it taro name [2] and is now cultivated in all tropical regions. Its spread has been almost universal in tropical Asia and even Japan before reaching Arabia and the Mediterranean basin at the beginning of our era and Tropical Africa later [3]. Hence, this plant was introduced to the Caribbean with African slaves.

This is an important staple food and subsistence for millions people in the countries of Africa, Asia and Central America [4]. Its worldwide production is estimated at 12 million tons on a cultivated area of 2 million hectares [5,6]. According to IITA [7], 77% of the world taro production comes from sub-Saharan Africa. This is the fourteenth most consumed vegetable in the world [7]. In 2012, global production was estimated at 10 million tons. Cameroon was the fourth largest producer of taro and the third largest producer in Africa after Nigeria and Ghana with a production of 1.6 million to 7 million tons for all of Africa [8,9,10,5]. Taro is grown in all regions of Cameroon for its leaves and tubers that have good nutritional qualities (easily digestible starches, vitamin C). It occupies an important place in traditional ceremonies in some African and Asian populations [12,13,14]. Tubers and taro leaves also have medicinal properties against tuberculosis, ulcers, pulmonary congestion and fungal infections [15,16]. The easily digestible starch tuber is an excellent food for diabetics [17]. A paradigm was observed in Japan and Egypt where taro has been used satisfactorily as a first harvest in the rehabilitation of saline soils [18]. This certainly opens the possibility of using the taro to subtract some ecological constraints.

Like humans and animals, plants must cope with new diseases. These diseases are currently one of the major threats to public health, sustainable agriculture, and food and ecosystem conservation [19,20]. The current environment is more conducive than ever to the appearance of new diseases: changes associated with globalization, especially increased trade and intercontinental economic climate modification, have important implications for the evolution of pathogens and dispersion of infectious agents [19,21].Thus in Cameroon since 2010, growing taro suffer of an epidemic disease; mildew caused by P. colocasiae [22]. This disease has reduced by over 70% taro production during the last three years in Cameroon (currently it is completely absent on the shelves in the markets) [23]. Although the impact of this disease on the population remains to be assessed, experience has shown that it can affect remarkably farmers' incomes and food security in Cameroon. Conventional modes of managements of this fungus that are far, widely used in agriculture require the use of phytosanitary products in large quantities. At the same time, measures are taken globally to encourage the farming community to limit its use. The avowed aim is to make all crops less dependent on fungicides to meet the ecological balance and preserve the environment. Thus, several measures have been taken as Ecophyto 2018 resulting from the Grenelle Environment Forum, the Framework Directive 2000 water or eco-conditionality of aid for the future CAP 2013 to comply with mandatory minimum standards. In addition, the massive and uncontrolled use of systemic fungicides has led to the birth P. colocasiae isolates resistant to these active ingredients. However, the prophylactic struggle and sanitation modes can represent popular among farmers because considered less effective. As for resistant cultivars, they are certainly as a durable solution in most taro production systems because the best control option is the spread and use of resistant cultivars [24]. Indeed, Cameroon, the fourth largest producer in 2010 with 1.6 million annual tons or about 13% of world production has since then losses ranging from 70 to 77%. This general observation demonstrates the need to provide solutions to farmers, firstly to better help they target their phytosanitary treatment in order to respect the environment and secondly to better explore the favourable traits for selection taro varieties; because according to Tansao [25], there is a high genetic variability of taro within and between countries. Any time dynamics of disease understanding and ability to evolve requires the study of the problem in its global dimension and must take into account the complex interactions between three key determinants: the characteristics of the pathogen, sensitivity of the host and the environment. It is in this light that we have created model taro blight in order to understand the evolution of this disease in one or more plots of taro to stem or stop this emerging epidemic in Central and West Africa.

In order to predict the attacks of P. colocasiae, we will propose a temporal model that will allow for a clarification on the level of infectious lesions may appear on taro leaves before the development of tubers depending estimated parameters. The efficiency of infection of the zoospore and sporangium, defined as the probability of penetration of a zoospore or sporangium lying on a healthy leaf. Knowing of course that the success of infection depends on the number of sporangia produced by infectious lesion also called infectivity of the pathogen, the latency period, the number of gendered spores or oospores produced by infectious lesion. The aim of our study is to develop a mathematical model to use as strategic as planned order. A computer program then will simulate the state of evolution of the number of infectious lesions that represents the intensity of the disease in a secluded plot. From the scientific point of view, the construction of such a mathematical model in plant pathology will generate the information necessary to learn about the operation of the host-parasite relationship, calculate the evolution of the intensity of the disease on taro leaves. This study will also contribute to help, for laboratory experiments on the development of P. colocasiae by varying the parameters related to this pathogen.

2. Material and Methods

2.1. Material

2.1.1. Biological Model

The term model has become very usual, both in research circles in crop protection, than in the agricultural council [26]. A model is to glue a mathematical theory on a piece of reality [27]. He tries to explain the functioning of a system through mathematical equations that relate the input and output components of the system [28]. This trend of modelling shows the determination of scientists trying to understand the dynamics of complex systems that involve the various processes. In plant pathology, models are built from key components of the development of the pathogen and the host-pathogen interaction [29]. Thus, for the parasitic plants fungi, the wetness is for example often paramount, because it affects spore germination and penetration of the pathogen into the host [30]. Accordingly, for epidemiological models in plant pathology, it is essential to integrate meteorological parameters in the data entering the system to predict the risk of infection [26].

Epidemiology is for two interacting populations: the population of the host and the parasite [29]. Thus, in our case, are used simultaneously to two levels of integration: the leaves, representing the host, and lesions, representing the parasite. The model uses a series of variables that at a given time can in principle be measured to describe the dynamics of the system. It is here the number of healthy leaves (H), the number of latent lesions (L), the number of infectious lesions (I), the number of sporangia (Sp), and the number of zoospores (Zp) present in the system and the proportion of oospores (O) products in the system. Particular attention was paid to the flexibility of this new model so that it can easily subject to change at the discretion of advanced research.

2.1.2. Biological and Data-Processing Material

The first mathematical model suggested in this work seek coupling the vegetative growth of the seedlings of taro and cycles it life of P. colocasiae, we have for this reason drawn using the software (illustrator) a cycle of probable life (Fig. 1) of this fungus, that of the taro being known the latter we will be used as biological material. Computer equipment consists of a portable computer of mark DELL of 2Go of RAM, with triple processor functioning with an operating system Windows 10; virtual machine JAVA version 7.0.450. Software Latex for the seizure, of the compiler winEdit and the software Scilab version 2008 for the digital simulation.

Fig. 1. Cycle life of P. colocasiae.

2.1.3. Inputs and Outputs of the Model

Build a model requires understanding the conditions under which it can be used. The methods of measurement and observation of model inputs should be generalizable, easy to implement and especially inexpensive. The model runs on a daily time step, but it is important to note that the IT point of view, it was designed to accommodate any changes in no time calculation. It uses data measured: the latency period, lesion size, the efficiency of infection of zoospores and sporangia.

As output, the model gives the incidence or severity of the disease on the leaves. This data will be faced with field observations to validate the calculations. This will also help to develop the favourable traits for plant breeding in order to decide that range set up to advise farmers.

2.2. Methods

2.2.1. Model Construction

The dynamics of P. colocasiae is governed by the following set of biological assumptions: (i) units of dissemination of P. colocasiae is the sporangia and the zoospore; (ii) it support of infection is a healthy sheet of C. esculenta; (iii) the sporangia can is to infect directly that is to say to produce n1 zoospore to increase the infections ; (iv) a lesion is latent during a certain period P then becomes infectious; (v) an infectious lesion produced n2 sporangia per day and can produce n3 oospores when the conditions are unfavourable; (vi) daily, d number of infectious lesions stops the sporulation and (vii) m healthy or latent sheets dies under the blow of the bad weather. We consider six distinct population, according to their disease status: likely sheets that one notes H (sheets which never was in contact with a sporangia or a zoospore including the young sheets), the latent ones that one notes by L (healthy sheets which were in contact with the units of contaminations of fungus but do not develop the disease immediately because the units of contamination (sporangia and zoospore) inside the healthy sheet carry out a latent life), infectious noted the I (infected sheets which already started to produce sporangia and which is able to transmit the disease), the sporangia or contaminant primary educations which we note Sp (represent the first units of dissemination of fungus through the culture and are carried per thousands with stroke of an infectious lesion), the noted zoospore Zp which are produced by the sporangia when the conditions are favourable, the oospores or spores sexual that we noted by O are those produced when the living conditions are unfavourable and which are close initiating the first infections when these conditions become favourable. The sheets of taro are produced with an intrinsic growth rate ε > 0 of which the number of sheets will not exceed the capacity limits K that a foot of taro can have, i.e., the healthy sheets are produced according to the logistic law:

The contamination being air mainly, the units of contamination are transported by the wind and the splashes of rains, we modelled the contact of susceptible (healthy sheet) with infectious (sporangia and zoospore) by the law of action of standard mass for a population with not constant dynamics what satisfies:

Where: β1 = e1g1 represents the force of infection of the zoospore; β2 = e2g2 represents the force of infection of the sporangia; D1 proportion of zoospore which gives a force of infection of 50% and D2 proportion of sporangia which gives a force of infection of 50%.

The natural mortalities occur at constant rate noted (m), but the sick sheets have an additional mortality which one notes (d) due to fungus.

A schematic model flowchart is depicted in (Fig. 2).

The dynamics of P. colocasiae, in a homogeneous population, is then described by the following system of nonlinear differential equations:

Fig. 2. Structure of the model.




Table 1. Variables and parameters with units for model system (1).

Symbols Description Units
g1 Proportion of zoospore able to germinate on a healthy sheet  
g2 Proportion of sporangia able to germinate on a healthy sheet  
e1 Effectiveness of infection of zoospore or severity 1/day
e2 Effectiveness of infection of the sporangia 1/day
n1 Quantity of zoospore released by a sporangia 50/day
n2 Quantity of sporangia produced by a lesion infectious 300000/day
n3 Proportion of oospores produced by an infectious lesion  
η Quantity of sporangia which releases from the zoospore  
mH Life expectancy (Longevity) of the healthy sheets 40 days approximately
mL Life expectancy (Longevity) of the sheets infected 10 to 20 days
mZP Rate of Zoospore washed on the sheets  
d Quantity of infectious lesions who stop the sporulation per day  
K Number of sheets which a seedling of taro can to have 6 to 7 sheets
α Duration of the latency period of fungus 24h to 18 days
D1 Proportion of zoospore which gives a force of infection of 50%  
D2 Proportion of sporangia which gives a force of infection of 50%.  


Birth rate of the sheets Day
δ Proportion of oospores washed on the sheets  

2.2.2. Basic Properties

Herein, we study the basic properties of the solutions of model system (1), which are essential

in the proof of the stability result. We have the following result.

Theorem 1: Model system (3) is a dynamical system on the biologically feasible compact domain:

Proof: The proof is provided in two steps.

Step 1: We show that the solution (H (t); L (t); I (t); Sp (t); Zp (t); O (t)) of model system (1)

corresponding to initial conditions such that H (0) > 0; L (0) > 0; I (0) > 0; Sp (0) > 0; Zp (0) > 0 and O (0) > 0 are not negative. Define

=sup{s > 0, /H(s) > 0, L(s) > 0, I(s) > 0, O(s) > 0, Sp(s) > 0, Zp(s) > 0} [0, t].

If τ = ∞, then, all the solutions of model system (1) are positive.

Let us apply the formula of Duhamel to prove the positivity of each equation of model system



In the same way, it is shown that


Third equation of the model system (2), it comes that

From there, one has

It comes whereas,

Fifth equation of the model system (2) one has

One can then conclude that

Lastly, of the last equation of the model system (2), one has

I.e. this completes the proof.

Step 2: It is proven that the model system (1) admits a single maximum solution.

Let us pose X (t) = (H (t); L (t); I (t); O (t); Sp (t); Zp (t)) and f a function of class  of  in  defined by:

The function f(X (t)) being infinitely differentiable on , it is consequently locally Lipschitzienne.

One considers the problem of following Cauchy:

Since f is of class C then; it admits a single maximum solution. Combining Step 1 and Step 2, Theorem 1 follows from the classical theory of dynamical systems. This concludes the proof.

2.2.3. The Disease-Free Equilibrium (DFE) and Its Stability

The disease-free equilibrium for model system (1) is given by



Then that H0 > 0, this point exist if and only if ε > mH.

In order to investigate the stability properties of the disease-free equilibrium X0, we need to compute the basic reproduction or threshold number R0 of model system (1). To this end, we define the following Jacobins matrix:

The matrix Jacobienne having a triangular structure. So that this commonplace point is stable, it is necessary that the various eigenvalues of this Jacobienne matrix are strictly negative real parts. Thus, the eigenvalues of the matrix  are and the eigenvalues of the new following matrix:

This new matrix  can cut out as it follows:

Thus, according to Kamgang and Sallet (2008), the new matrix  is stable if and only if matrix (is stable. A simple calculation gives

From there, it comes that

The matrix  being of order 2, it is stable if and only if its trace is negative and its determinant is positive, i.e.  and  From of the matrix, one has:


It is noticed that  involve


Indeed,  involve


This completes de proof.

Thus the matrix is stable if and only if  i.e.

Thus, the basic net reproduction rate of model system (1) is


It measures the number of new infectious lesions generated by an infectious starting lesion in a population entirely made up of healthy sheets. Consequently, the DFE of our model is locally asymptotically stable for R0 ≤ 1 and unstable for R0 > 1. Therefore, any fight planned to make disappear the mildew should seek to lower this rate below 1. We then showed the following result.




2.2.4. The Endemic Point of Balance and Its Stability

That is to say X = (H*; L*; I*; O*; Sp*; Zp*) a point of endemic balance unspecified of the model system (1). This point is obtained by posing the entire derivative equal to 0, which gives us.


By using the first, third, fourth and fifth equation of the model system (5); (I*), (O*), (Sp*), (Zp*) and (H*) express themselves according to (L*) as follows.




By integrating the expressions of (I*), (O*), (Sp*), (Zp*) and (H*) given to the equation (6) in the second equation of the system (5), one obtains the following equation of degree four in (L*):





The rule of the signs of Descartes will thus enable us to discuss on the existence the possible solutions this equation.

Theorem When R0 > 1 and R0 close to 1, the model of the system (1) admit locally a point of endemic balance asymptotically stable.

Table 2. Solutions of the equation (8).

A4 A3 A2 A1 A0 R0 Number of solutions
+ + + + +


0 solution
+ + + + -


1 solution
+ + + - +


2 solutions
+ + + - -


1 solution
+ + - + +


2 solutions
+ + - + -


3 solutions
+ + - - +


2 solutions
+ + - - -


1 solution
+ - + + +


2 solutions
+ - + + -


3 solutions
+ - + - +


4 solutions
+ - + - -


3 solutions
+ - - + +


2 solutions
+ - - + -


3 solutions
+ - - - +


2 solutions
+ - - - -


1 solution

2.3. Numerical Simulations

Fig. 3. Simulation of model system (1) when (ε = 0.3; K = 3000; β = 10-6; D1 =45000; D2 = 45000; mH = 0.01; α = 0.3; mL = 0.01; δ = 0.4; n1 = 15; d = 0.3; n3= 100; n2 = 300000; ƞ = 0.7; mZp = 0.5; ε = 0.6; (R0 = 0.153) < 1) (a) healthy sheets; (b) infected sheets; (c) infectious sheets; (d) oospore numbers; (e) sporangia numbers; (f) zoospore numbers.

Fig. 4. Simulation of model system (1) when (ε = 0.3; K = 3000; β = 10-1; D1 =450000; D2 = 450000; mH = 0.01; α = 0.2; mL = 0.01; δ = 0.5; n1 = 15; d = 0.7; n3= 100; n2 = 300000; ƞ = 0.5; mZp = 0.4; ε = 0.3; (R0 = 7.66) >1) (a) healthy sheets; (b) infected sheets; (c) infectious sheets; (d) oospores numbers; (e) sporangia numbers; (f) zoospore numbers.

2.4. Model Multi-patches

2.4.1. Biological Facts

It is supposed that at one moment t, a piece is in one of the following states: healthy sites (H), latent sites (L), infectious sites (I).

The structure of the model is given to the Fig. 5.

Fig. 5. Structure of the model.

If one takes v(i) variety of taro in the piece i. An infectious lesion produced rv(i) sporangia per day. The sporangia are dispersed on another piece j with a probability μij. Thus, at time t, the number of sporangia dispersed starting from the piece i on other pieces follows a distribution multinomial


The number of sporangia arriving in the piece j coming from the piece i at time t is.


 Infection of the healthy sites: a healthy sheet in an unspecified piece receives a sporangia or a zoospore then becomes infected with the probability βv(i). In the same way, a piece can undergo harvests with a force of v (i)

When a sheet is infected, the mushroom multiplies inside fabric during this time one speaks about latent lesion. The appearance of a new latent lesion follows the following binomial distribution:


There is the following recurring equation.


Sporangia and the zoospore which do not manage to infect are withdrawn from the system with the respective probabilities μ and mZp.

Transition from the sites of latent to the infectious one: After an infection, the healthy sites become latent lasting Pv(i) days (latency period) before producing new sporangia

Of this fact, the probability so that a new infectious lesion appears in a piece i between time

t-1 and time t is . This appearance follows a binomial distribution then such as:


One obtains the following recurring equation:


Passage of the infectious sites to the sites whose lesion ceases the sporulation, after tv(i) days of sporulation (period of infection) the site is regarded as withdrawn of the system studied

There is the following recurring equation:


Dispersion: We admit that the sporangia move at long distance by the wind. The density of spores emitted by a source z located at the point (x; y) and which are received at the point z, located in the zone (x,; y,) is given by (Klein et al., 2006):


Where the Euclidienne distance is enters z and , b > 0 are a scalar parameter and a > 2 determine the distance from dispersion.

The average distance is only defined when a > 3 then that

The probability for a spore of dispersing of a piece i a piece j is


Where, Ai and Aj are respective surfaces of the pieces i and j.

2.4.2. Model



Thus, the basic net reproduction rate of the model (20) is


2.4.3. Numerical Simulations

(a) healthy sheets; (b) infected sheets; (c) infectious sheets; (d) sporangia numbers; (e) zoospore numbers.

Fig. 6. Simulation of model system (20) when (i = 1; λ= 1000; β = 10-8;; α = 0.1; μ= 0.5; mL = 0.01; d = 0.5; D = 0.1; r = 300000; ƞ = 0.7; n1 = 15; mZp = 0.5; D1 = 45000; D2 = 45000; (R01 = 0.00017) < 1).

(a) healthy sheets; (b) infected sheets; (c) infectious sheets; (d) sporangia numbers; (e) zoospore numbers.

Fig. 7. Simulation of model system (20) when (i = 2; λ= 1000; β = 10-5;; α = 0.1; μ= 0.5; mL = 0.01; d = 0.5; D = 0.1; r = 300000; ƞ = 0.7; n1 = 15; mZp = 0.5; D1 = 45000; D2 = 45000; (R02 = 0.17) < 1).

(a) healthy sheets; (b) infected sheets; (c) infectious sheets; (d) sporangia numbers; (e) zoospore numbers.

Fig. 8. Simulation of model system (20) when (i = 1; λ= 1000; β = 10-2;; α = 0.1; μ= 0.5; mL = 0.01; d = 0.5; D = 0.1; r = 300000; ƞ = 0.7; n1 = 15; mZp = 0.5; D1 = 4500000; D2 = 4500000; (R01 = 1.33) > 1).

(a) healthy sheets; (b) infected sheets; (c) infectious sheets; (d) sporangia numbers; (e) zoospore numbers.

Fig. 9. Simulation of model system (20) when (i = 2; λ= 1000; β = 10-2;; α = 0.1; μ= 0.5; mL = 0.01; d = 0.5; D = 0.1; r = 300000; ƞ = 0.8; n1 = 15; mZp = 0.5; D1 = 4500000; D2 = 4500000; (R02 = 1.7) > 1).

3. Results and Discussion

The model is not necessarily based on a complex mathematical formulation. It is built from information that has been acquired on the system studied and allows to simulate the operation [33]. To explain the behaviour of the disease, we make numerical simulations that explore the epidemic by varying the conditions of development of the epidemic. Hen t is expressed in days, R0 the multiplying daily factor, that is to say the number of girls lesions initiated daily by a parent injury.

The model (Fig. 2) shows that the latency period, lesion size, the efficiency of infection of sporangia and zoospores have significant effects on both the speed of the epidemic on the maximum level.

3.1. Biological Interpretation of the Fig. 3

Fig. 3 presents six compartments: (a), (b), (c), (d), (e), (f).

1) The compartment (a) represents the evolution of the number of healthy leaves in a plot in the absence of disease because (R0 <1). To follow this development we set the maximum number of healthy leaves that may have a plot in 3000 (this is an example) that said, whatever the number of healthy leaves we at the beginning of the experiment, should expect that number to reach and does not exceed 3,000 sheets taking into account the constant renewals. However, the compartment (a) shows that 40 days after planting taro bulbs, the maximum number of sheets is reached, but it is not equal to 3,000 sheets as expected because healthy leaves if not infected die naturally. So we can have a complete defoliation of the field if weather is severe.

2) The compartment (b) is in direct communication with the compartment (a). This means that, if there is no infection in the compartment (a), and of this compartment there is no appearance even if we were to introduce into the plot of sick leaves in order to initiate infections, it will be observed that after 20 days the leaves die without even initiate infections.

3) (c) which represents the evolution of the number of sheets capable of transmitting the disease, there is no sheet capable of transmitting the disease because of (a) there is no infection.

4) (d) represents the number of oospores the oospores is fungal survival unit into the ground. So compartment (d) reports that if the previous crop due to leave millions of oospores in soil, these oospores transported and deposited on healthy leaves, due to being able to initiate damage, they will start die from the 20th day to the point of completely disappearing the 25th day.

5) Compartments (e) and (f) represent the number of sporangia and zoospores, describing the situation that these cannot be produced by the experimental culture, may come from the surrounding plots. Any time when placed on healthy leaves, due to not being able to initiate infection, they will start dying from the 20th day as oospores.

3.2. Biological Interpretation of (Fig. 4)

In Fig. 4, R0> 1, so the illness persists. But this disease cannot grow indefinitely. It will stabilize at a point that we call endemic disease.

But what is the biological information that we give these curves?

1) In the compartment (a), with the same lines as before; this time the number of healthy leaves is rather stable in 1000 still from the 40th day. This means that 2000 leaves are dead under the influence of weather and largely under the influence of the fungus. It appears here that there is for any culture of taro, a number of sheets for which although the variety of taro is sensitive, if this number is not reached the disease does not develop.

2) The compartment (b) is in communion with the compartment (a). So we should expect that 2000 healthy leaves infected (a) find themselves in (b) not. The block (b) presents the appearance rather than 1000 infected leaves only; therefore, infected leaves can also die naturally without any time to cause new infections.

It appears here that whatever the speed of disease progression, it should stabilize from the 45th day. In other words, whatever the number of infected leaves, the disease will cease to grow of the 45th day. So in all taro crops whatever the speed of evolution of the disease, there are a limited number of infected leaves that may appear during the disease. This number is in this case 1000.

3) The compartment (c) and (d) are similar since they are in direct communication. So what comes out of (b) is found in (c) without any time to forget that an infectious sheet can naturally die without having had time to infect another.

4) (d), (e), (f), there are also biological stabilities thresholds. So, whatever the level of disease, each culture can only produce a given number of oospores, sporangia and zoospores. Knowledge of the day when that number is reached is very useful for plant pathologists; in our case this value is reached the 45th day.

The values given in this simulation can be used as a scale for clarification on the field.

We say that as there is a very close relationship between the development of the disease and micro-meteorological variables such as temperature humidity and precipitation [32] as it exists between the evolution of the disease, the number of sheets, oospores, sporangia and zoospores. This relationship has helped to build predictive models of epidemics and thus identify the periods of maximum risk [33,34,35,36,37,38,39,40,41,42,43]. The ultimate goal of these models is to target periods or intervention is necessary and thus reduce the number of fungicide application.

Despite these models, it is still difficult to predict the contamination dates for two main reasons: the first is the heterogeneous distribution of primary sources of inoculum; the second is that these models do not incorporate conservation factors or germination and infectivity or oospores [44]. According Andrivon and Lebreton [44], the development of a reliable and realistic model requires first to fill the lack of information about the biology of these organs.

4. Conclusion and Perspectives

In this work, we studied in leading the influence that can have a favourable trial setting in the kinetics of the disease, the end of which we believe that the duration of the pathogen latency period, the number of sporangia products on the surface of a lesion as the efficiency of infection of the zoospore or sporangia individual parameters are taken into account in the development of a variety of taro appropriate level of resistance. But these parameters are exacerbated when taken together.

The model is then used to compare the impact of different resistance components on the evolution of the disease, and therefore to optimize varietal selection: choice of selection goals based on their ability to practice and assessment their epidemiological effects.

In the background, we addressed the question of the benefit of a diversified operation taro plant in the mushroom control P. colocasiae so the idea is to grow in the same plot several varieties of taro to minimize the severity of the disease. In this second model, the interaction between plots will be made only via sporangia, in perspective, it would be interesting for a complete model that strongly closer to reality also take into account the interactions through zoospores that can move from one plot to another; while considering the dispersion of these through culture, because this probability has only been a mere guess. Parameter estimation, analysis of the model parameters of sensitivity as well as field evaluation for model validation will be made. In developing countries, where climate data acquisitions are difficult or inputs and pesticides and fungicides are in certain circumstances prohibitively expensive, models that look for help in plant breeding are of great importance.


  1. Matthews, P.J. 1990. The origins, dispersal and domestication of taro. Ph.D. Thesis, Australian National University, Canberra, Australia.
  2. Purseglove, J.K. 1972. Tropical Crops. Monocotyledons I. Longman, London.
  3. Doku, E.V. 1980. Comment faire progresser la recherche sur les taros. In : Plantes racines tropicales: Stratégies de recherches pour les années 1980. Compte rendu du premier symposium triennal sur les plantes-racines de la société internationale pour les plantes-racines tropicales-Direction Afrique, du 08-12 Septembre 1980.(Eds), Ibadan, Nigéria, pp. 237-241.
  4. Erwin, D.C., Roberia, O.K. 1996. Phytophthora diseases worldwide. St Paul, Minnesota:APS Press. 562 p.
  5. FAOSTAT. 2011. FAO Economic and Social Department. The Statistics Division. MajorFood and Agricultural Commodities and Producers. http ://
  6. Scot, N., Brooks, F.E., Glenn, T. 2011. Taro Leaf Blight in Hawaii. University of Hawaiiat Mänoa, Plant Disease 71: 1-14.
  7. IITA. 2009. Root and Tuber systems. http ://
  8. Rao, V.R. 1996. Taro genetic resources. In: Proceedings of the Second Taro Symposium held 23-24 November 1994 at the faculty of Agriculture, Cenderawasih University, Manokwari, Indonesia.
  9. AGRISTAT. 2009. Annuaire des statistiques sur secteur agricole, campagnes 2006 à 2007. Ministère de l’agriculture et du développement rural. Yaoundé, Cameroun. 100 p.
  10. CTA (Centre Technique de Coopération Agricole et Rurale). 2010. Guide d’exportation pour les plantes à racines et tubercules en Afrique de l’Ouest et du Centre. Dakar, Sénégal. 32p.
  11. Food and Agriculture Organization of the United Nations.2014. Production de taro, base de donné
  12. Asseng, C.C., Mvoe, C.R., Ambang, Z., Monkam, T.F. 2016. Influence of the Number of Watering and Fungicide Treatments on the Development of Phytophthora colocasiae (Racib) on Cocoyam (Xanthosoma sagittifolium) and Taro (Colocasia esculenta) Greenhouse in Cameroon. Int.J.Curr.Microbiol.App.Sci. 5(8): 100-112.
  13. Onwueme, I.C. 1999. Taro cultivation in Asia and the Pacific. FAO Regional Office for Asia and the Pacific, Bangkok, Thailand.
  14. Caillon, S. 2005. Les taros du Vanuatu : que conserver et comment?Natures Sciences Sociétés 13: 306-310.
  15. Misra. R.S., Sriram, S. 2002. Medicinal value and export potential of tropical tuber crops. In Recent Progress in Medicinal Plants; Crop improvement, Production Technology and Commerce (vol. 5), Govil JN, Pandey J, Shivkumar BG, Singh VK (eds). SCI Tech Pub: USA; 376-386.
  16. Binoy, B., Hegde, V., Makeshkumar, T., Jeeva, M.L. 2010. Rapid detection and identification of potyvirus infecting Colocasia esculenta (L.) Schott by reverse transcription polymerase chain reaction. J. Root Crops 36, 88–94.
  17. Wang, J. 1983. Taro: a review of Colocasiae esculenta and its potentials. University ofHawaii Press. Honolulu, Hawaii 400p.
  18. Kay, D.E., 1973. Root Crops. In Crop and Product Digest. Ed. Tropical Products Institute. London. 245p.
  19. Anderson, S.P., Dunn, C., Laughter, A. 2004. Overlapping transcriptional programs regulated by the nuclear receptors peroxisome proliferator-activated receptor alpha, retinoid X receptor, and liver X receptor in mouse liver. Mol Pharmacol. 66: 1440–52.
  20. Kwon, Y.H., Epstein, L. 1993.A90-kDaglycoprofein associafed with adhesiono Nectria haematococca macroconidia to substrata. Mol. Plant-Micr. Interact. 6:481-487.
  21. Desprez-Loustau, M.L., Robin, C., Buee, M., Courtecuisse, R., Garbaye, J., Suffert, F., Sache, I., Rizzo, D. 2007.The fungal dimension of biological invasions. Trends. Ecol. Evol. 22: 472–480.
  22. Guarino, L. 2010. Taro leaf blight in Cameroon. 2010. Agricultural Biodiversity Weblog. http/agro. –leaf–blight - in-Cameroon/visited 15 May 2012.
  23. Fontem, D.A., Mbong, G. 2011.A novel epidemic of taro (Colocasiae esculenta) blight caused by Phytophthora colocasiae hits Cameroon. Third Life Science Days held at the University of Dschang, Cameroon. 26-28.
  24. Monkam, T.F., Asseng, C.C., Ngando, E.J., Noupadja, P. 2014. Evaluation of the Statistical Methods for the Classification of Plantain Hybrids According to their Sensibility to the Black Leaf Streak Disease (BLSD). Int. J. Res. Studies in Biosci., (IJRSB) Volume 2, Issue 11, PP 85-94.
  25. Tansao. 2000. Acessado 20 de janeiro de 2011.Annual Report. Taro network for South East Asia and Oceania: Evaluation and breeding for rain-fed cropping systems in South East Asia and Oceania, Port Vila.
  26. Savary, S. 1995. Modélisation en protection des cultures.Paris: ORSTOM, p. 141-156.
  27. Ruelle, D. 1991. Hasard et chaos.Edition Odile Jacob, Paris.
  28. Rouzet, J., Larguier, M. 1997. Lutte raisonnée: les différents apports des modèles de prévision selon les parasites et les cultures. ANPP, Cinquième conférence internationale sur les maladies des plantes, Tour, 3-5 Décembre 1997, pp.47-54.
  29. Vanderplank, J.E. 1963. Plant Diseases: Epidemics and Control. Academic Press New York.
  30. Van Maanen, A., Xu X.M. 2003.Modelling plant disease epidemics. European Journal of Plant Pathology.109: p 669-682.
  31. Lannou C. 1994. Utilisation de modèles de simulation pour la recherche en phytopathologie INRA, laboratoire de pathologie végétale, Thiverval-Grignon Actes d'un séminaire international, Montpellier, mai 1994.
  32. Harrison, J.G. 1992. Effects of the aerial environment on late blight of potato foliage. Plant Pathology. 41:384-416.
  33. Krause, R. A. Massie, L. B., and Hyre, R. A., 1975. Blight cast, a computerised forecastof potato late blight. Plant Disease 59: 95-98.
  34. Fry, W. E. 1977. Integrated control of potato late blight- effects of polygenic resistance and techniques of timing forecast. Phytopathology 74: 1145-1147.
  35. Levy, Y., Levi, R., and Cohen, Y. 1993. Build-up of a pathogen sub-population resistantto systemic fungicide under various control strategies: a flexible simulation model. Phytopathology 67: 675-677.
  36. Fohner, G.R., Fry, W.E., White, G.B. 1984. Computer simulation raises question about timing protect ant fungicide application frequency according to potato late blight forecast. Phytopathology 74: 1145-1147.
  37. Milgroom, M.G., Fry, W.E. 1988. A simulation analysis of the epidemiology principles for fungicide resistance management in pathogen populations. Phytopathology. 78: 565-570.
  38. Shtienberg, D., Doster, M.A., Pelletier, J.R., Fry, W.E. 1989. Use of simulation models to develop a low risk strategy to suppress Early and Late blight in potato foliage. Phytopathology 79: 590-595.
  39. Fry, W.E., Goodwin, S.D., Drenth, A., Tooley, P.W., Sujowski, L.S., Koh, Y.L., Cohen, B.A., Spielman, L.J., Deahl, K.L., Inglis, D.A. Sandlan, K.P. 1993. Historical and recent migrations of Phytophthora infestans: Chronology, pathways and implications. Plant Disease 77: 653-661.
  40. Shtienberg, D., Fry, W.E. 1990. Field and computer simulation evaluation of spray scheduling methods for control of early and late blight of potato. Phytopathology 80: 772-777.
  41. Duvauchelle, S., Andrivon, D. 1996. Maladies à distribution géographiques mondiale: Le mildiou et son agent Phyttophthora infestans (Mont.) de Bary. In La pomme de terre. Rousselle, P., Robert, Y. Croisnier J. C. ed.INRA et ITCF. France. 607p.
  42. Duvauchelle, S., Lherbier, V., Sarniguet, C., Lebreton, L., Andrivon, D. 1997.Evolution of phenylamide resistance in Frensh populations of Phytophthora infestans, 1981-1996.ANPP- Cinquième Conférence Internationale sur les Maladies des Plantes. Tours- 3-4-5 Décembre 1997, p: 471-478.
  43. JaeHoon, A., YoungIl, H., and Kwanyong. H., 1998. Modeling for prediction of potato late blight (Phytophthora infestans ) progress.Korean Journal of Plant Pathology 14:331-338.
  44. Andrivon, D., Lebreton, L. 1997. Mildiou de la pomme de terre, ou en somme-nous après 150 ans ?Phytoma, 494 (Mai 1997) p: 24-27.

Article Tools
Follow on us
Science Publishing Group
NEW YORK, NY 10018
Tel: (001)347-688-8931