Stability and Bifurcation Analysis of Discrete Dynamical Systems 2020
View this Special IssueResearch Article  Open Access
Tinashe B. Gashirai, Senelani D. HoveMusekwa, Steady Mushayabasa, "Lyapunov Stability Analysis of a Delayed FootandMouth Disease Model with Animal Vaccination", Discrete Dynamics in Nature and Society, vol. 2020, Article ID 3891057, 13 pages, 2020. https://doi.org/10.1155/2020/3891057
Lyapunov Stability Analysis of a Delayed FootandMouth Disease Model with Animal Vaccination
Abstract
Footandmouth disease virus remains one of the most important livestock diseases in subSaharan Africa and several Southeast Asian countries. Vaccination of livestock has been recognized as an important tool for the control of footandmouth disease virus. However, this intervention strategy has some limitations. Generally, vaccine production is a complex multistep process which involves development, manufacturing, and delivery processes, and through this extensive process, some challenges such as poor vaccine storage often arise. More often, these challenges alter the validity of the vaccination. Footandmouth disease virus epidemic dynamics have been extensively explored, but understanding the role of vaccination validity on virus endemicity is lacking. We present a timedelayed footandmouth disease model that incorporates relevant biological and ecological factors, vaccination effects, and disease carriers. We determined the basic reproduction number and demonstrated that it is an important metric for persistence and extinction of the disease in the community. Numerical illustrations were utilised to support some of the analytical results.
1. Introduction
Footandmouth disease (FMD) is one of the most important livestock diseases in subSaharan Africa and several Southeast Asian countries. It affects clovenhoofed animals and is caused by an RNA virus commonly known as footandmouth disease virus (FMDV). In many developed nations, the disease often occurs as epidemic (outbreak followed by extinction), whereas in many developing nations, such as Zimbabwe and Zambia and several Southeast Asian countries including Cambodia, Lao PDR, Malaysia, Myanmar, Thailand, and Vietnam, it is epidemic (persistent diseaseās presence) [1]. Epidemic or endemicity of the disease negatively affects agricultural economic world [2].
Among several other control strategies such as movement restrictions, public education, and culling, vaccination of animals has been recognized as the most important strategy to control the spread of the disease [3]. However, there are many limitations that affect vaccination validity period. For instance, cultivation of FMDV in the manufacturing units poses a risk of escape from production sites. Furthermore, vaccines may sometimes contain traces of FMD viral nonstructural proteins (NSPs), therefore, interfering with the NSPbased serological differentiation of infected from vaccinated animals [4]. In addition, vaccine production is a complex multistep process which involves development, manufacturing, and delivery processes, and through this extensive process, some challenges such as poor vaccine storage often arise. More often, these challenges alter the validity of the vaccination. For example, a study carried out in one Cambodian province recorded that more than half of cattle vaccinated with donated FMD vaccines are subsequently becoming infected with FMD virus and showing clinical signs of FMD indicating possible vaccine failure [1, 5]. Moreover, most vaccine production plants belong to companies, usually driven by profit rather than disease control or eradication; hence, most of the farmers in FMD endemic countries cannot afford the cost of these highquality vaccines leading them to rely on donated or cheap vaccines obtained from unregistered suppliers, thereby compromising the outcome of vaccination programs [6, 7].
Since the 2001 FMD outbreak in UK, several mathematical models have been proposed to study the transmission and control of FMD (see, e.g., [6, 8ā21]). Dexter [8] utilised stochastic epidemic models to explore the dynamics of FMD among feral pigs in the Australian semiarid rangelands. The author concluded that any control program designed to suppress FMD outbreak in feral pigs in the semiarid rangelands of Australia should take account of the prevailing environmental conditions. Using an SIR (susceptibleinfectiousrecovered) epidemic model, Bravo de Rueda et al. [9] demonstrated that both direct and indirect disease transmission pathways play a crucial role in FMD transmission. Utilising an nonautonomous mathematical model, Mushayabasa et al. [17] showed that seasonal variations may alter the dynamics of FMD. Recently, Zhang et al. [21] developed a dynamical switching model for FMD that incorporated three prevalent serotypes of the virus. The model was used to assess the spread of FMD in mainland China. Among several other important outcomes, their study showed that, under certain conditions, some of the FMD serotypes persist while others may disappear.
Despite all these efforts, however, mathematical models that aim to explore the role of vaccine waning in FMD endemic regions are still lacking. Therefore, in this study, we developed a timedelayed FMD model to explore the effects of vaccine waning on the dynamics of the disease. The proposed model incorporates all the relevant biological and ecological factors, vaccination, and disease carriers. We are aware that the role of FMD carriers in disease transmission remains a debatable issue [17, 22] and as such we have explored the outcomes of the model in two scenarios: (i) no transmission by carriers and (ii) when carriers are transmitting the disease. Time delay plays an important role in modeling the dynamics of infectious diseases, and the dynamical behaviour of an epidemic model with time delay(s) is known to be complex and difficult to comprehend [23]. In the proposed framework, the time delay represents the validity of the vaccination period.
We organize this paper as follows. In Section 2, we formulate and analyze the timedelayed FMD model. In particular, we computed the basic reproduction number and demonstrated that it is an important metric for the existence and global stability of the model steady states. In addition, we also explored the model numerically to support the analytical results. In the last section, we discuss the insightful implications of our findings on the role of vaccination validity on FMD transmission and control.
2. Methods and Results
In this section, we propose and analyze a timedelayed footandmouth disease model incorporating direct and indirect disease transmission pathways, vaccination effects, and virus excreting carriers. We will investigate the qualitative behaviours of the proposed framework through studying the stability of the model steady states.
2.1. Model Framework
Let the variables , , , , and represent the densities of susceptible, exposed (latently infected), clinical infectious, carrier, and recovered cattle at time , respectively, such that the total population of cattle at any given time is . In addition, let the variable denote the density of the FMDV excreted into the environment by clinically infected and carrier animals. The proposed model is as follows:subject to the initial functions as follows:whereare continuous on .
All model parameters are assumed to be positive and are defined as follows: denotes the constant recruitment rate of animals through birth and they are assumed to be susceptible to infection. All animals are assumed to have a life span of days, despite the epidemiological status of the animal. Susceptible animals are assumed to acquire infection either directly through contact with clinically infected animals or carriers at rate, , or indirectly when they inhale the virus excreted into the environment by the aforementioned infectious animals; at rate , accounts for differential infectivity between clinically infected animals and carriers. Thus, implies that clinically infected animals are assumed to be more infectious than carriers, and the reverse will be true, for . Model parameter is the proportion of vaccinated animals, the term represents the vaccinated animals, and is the validity period of the vaccination. Animals exposed to FMDV incubate the disease for a period ranging between 2 and 14 days [6, 16ā18]. Here, represents the average incubation period.
FMD is rarely fatal and infected cattle generally clear the systemic infection within 8ā15 days [22]. However, unlike pigs which usually clear FMDV within 3 weeks following infection and do not become carriers [22, 24ā26], up to 50% of FMDrecovered ruminants become persistently infected [22]. Hence, we have assumed that, after an average infectious period of days, a fraction of clinically infected animals becomes carriers and the remainder successfully clear the infection. Furthermore, compared to pigs, ruminants like cattle are highly susceptible to FMDV infection via the inhalation route [22, 27]. Clinically infected and carrier animals shed the virus in the environment at rates and , respectively, where is a dimensionless parameter that accounts for the difference in pathogen shedding rate between animals from the aforementioned infectious classes. We assume that it takes days for the pathogen in the environment to clear the infection. Furthermore, carriers are assumed to clear the infection after a duration of days. In this study, FMD carriers are animals from which the virus can be recovered more than 28 days after infection.
2.2. Properties of the Solutions
The following theorem shows that the model proposed in this study is biologically meaningful. Precisely, the theorem demonstrates that all the solutions of the proposed model are nonnegative and bounded for all . Since the five equations of model (1) do not depend on (2), it suffices to study the dynamical behaviour of the disease without this equation.
Theorem 1. There exists a unique solution for the FMD model (1). Furthermore, the solution is nonnegative for all and lies in the set
Proof. In order to demonstrate that the solutions of model (1) are nonnegative, we will investigate the direction of the vector field given by the righthand side of (1) on each space and note whether the vector field points to the interior of or is tangent to the coordinate space. Sinceit follows that the vector field given by the righthand side of (1) on each coordinate plane either is tangent to the coordinate plane or points to the interior of . Hence, the domain is a positively invariant region. Moreover, if the initial conditions are nonnegative, then it follows that the corresponding solutions of model (1) are nonnegative. Furthermore, let and it follows thatSince , it follows by applying the standard comparison [28] that . In particular, we if . Hence we can conclude that the animal population is bounded. From the last equation of system (1),Similarly, by the standard comparison theorem [28], we have . Therefore, , if and this implies that the pathogen population is also bounded. Hence, all solutions in eventually enter .
2.3. Stability and Bifurcation Analysis
In this section, we investigate the qualitative behaviours of the delay differential model (1). In particular, we will determine the basic reproduction number, existence, and global stability of the model steady states. Through direct calculations, one can easily note that model (1) admits two equilibrium points, namely, the diseasefree and the endemic defined bywithwhere , the basic reproduction number which is given by (1), iswith
The threshold quantity is defined as the average number of secondary footandmouth disease cases generated by a typical infected animal throughout its entire course of infection in a completely susceptible population in the presence of animal vaccination. The quantities and represent the basic reproduction number associated with direct and indirect disease transmission pathways, respectively. The basic reproduction number is an important metric for infectious disease models; it demonstrates the power of the disease to invade the population. From (17), we can observe that the disease persists if and only if and it dies out for .
Next, we investigate the global stability of the diseasefree equilibrium, using a Lyapunov functional. We define a continuous and differentiable function as
We note that for any and has the global minimum 0 at .
Theorem 2. Let be a closed and bounded positively invariant set with respect to system (1) [29]. Let be continuous and differentiable function such that and in for solutions of system (1). Let . Let be the largest invariant set in ; then every bounded solution starting in approaches as . Particularly, when , then , as .
Theorem 3. For system (1), when , the diseasefree equilibrium point is globally asymptotically stable.
Proof. To establish the global stability of the DFE, we consider the following Lyapunov functional:whereTaking the derivative of along the solutions of system (1) leads toAt diseasefree equilibrium, we have the identity . Using this identity and simplifying, one obtainsSince the arithmetic mean is greater than or equal to the geometric mean, we haveand it follows thatFurthermore, we define a continuous and differentiable function , for . It follows that , with the equality satisfied if and only if . Thus, we can note thatFrom the earlier discussion on initial conditions (3) as well as the properties of solutions of model system (1), we can have for and this implies that , . Hence, it follows that .
Hence, if , we have . Furthermore, let be the largest invariant subset of the set . We now claim that . In fact, when , it follows from (17) that , which leads to . If and from the first and second equations of system (1), we have . Again, we have . Noting that is invariant, by Lasalleās invariance principle [30], the diseasefree equilibrium is globally asymptotically stable whenever .
In what follows, we shall investigate the global stability of the endemic equilibrium (we have shown that it exists whenever , i.e., , , , , and if and only if ) by using the Lyapunov direct method. The construction of the Lyapunov functional is also based on the function , by which we have already observed that it satisfies for any and has the global minimum 0 at .
Theorem 4. Let be a closed and bounded positively invariant set with respect to system (1) [29]. Let be continuous and differentiable function such that and in for solutions of system (1). Let . Let be the largest invariant set in ; then every bounded solution starting in for approaches as . Particularly, when , then , as .
Theorem 5. If , then the endemic equilibrium of model system (1) exists and is globally asymptotically stable.
Proof. We consider a Lyapunov functional as follows:wherewithFrom Section 2.2, properties of solutions of model system (1), is well defined and , . The equality holds if and only if , , , , and .
Taking the derivative of along the solutions of system (1), one obtainsAt the endemic equilibrium of model system (1), we have the following identities:Utilising the identities leads toOnce again, since the arithmetic mean is greater than or equal to the geometric mean, we haveand it follows thatFurthermore, as already illustrated on demonstrating the global stability of the diseasefree equilibrium,Moreover, using the function , defined earlier, we haveFurthermore,Hence, from the illustrations (28)ā(31), if , we have . By Theorem 5.3.1 in [31], solutions are limited to , the largest invariant subset of . Therefore, for all provided that are nonnegative, it follows by Lasalleās invariance principle [30] that the endemic equilibrium point is globally asymptotically stable whenever .
2.4. Numerical Results
To explore the behaviour of system (1) and illustrate the stability of equilibria solutions, we utilised dde23 [32] based on RungeKutta methods. The baseline values for the model parameters are in Table 1, and the initial population levels were assumed as follows: , , , , and .

Simulation results in Figures 1ā3 illustrate the dynamical solutions of model (1) under the assumption that FMD carriers are noninfectious; that is, . In all the cases, we can note the existence of periodic oscillations for a certain period, before stability at endemic point is attained. We can also note that the intensity and amplitude of oscillations decrease as the vaccination validity period increases, implying that if the vaccine waning period is āreasonably longer,ā then the possibility of having periodic oscillations on the solutions of system (1) decreases.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
Graphical results in Figures 4ā6 demonstrate the potential impact of infecting FMD carriers on altering the dynamics of footandmouth disease. The contribution of FMD carriers on the disease dynamics remains a debatable issue. Some researchers argue that these animals are noninfectious, while others postulate that they are infectious [22]. However, a recent study [22] demonstrated that FMD carrier cattle are capable of shedding the virus into the environment. Thus, in Figures 4ā6, we fixed the vaccination validity period to 30 days and explored the role of FMD carrier cattle on the dynamics of the disease by varying model parameters associated with increasing or decreasing the infectivity of the FMD carriers. In Figure 4, we set , , and , leading to . We observed that if the proportion of animals which become carriers, , is large while the contribution of carriers to new infections (through transmission, , and pathogen shedding into the environment, ) is less than that of clinically infected animals, then intensity and amplitude of periodic oscillation in the solutions of model system (1) decrease. However, results in Figures 5 and 6 suggest that if the proportion of animals that become carriers decreases, then the intensity and amplitude of periodic solutions increase despite irrespective of whether carriers are less or more infectious compared to clinically infected animals. Overall, the results demonstrate that understanding the contribution of carriers on FMD dynamics is an important issue that requires multidimensional research.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
3. Concluding Remarks
Footandmouth disease, a highly contagious, acute viral disease of clovenhoofed domestic and wild ruminants including cattle, buffalo, swine, goats, and sheep, remains a major challenge in many developing nations in the subSahara Africa, such as Zimbabwe and Zambia. Vaccination is regarded as the most effective means of footandmouth disease mitigation. Nevertheless, this important intervention strategy has some limitations. Generally, vaccine production is a complex multistep process which involves development, manufacturing, and delivery processes, and through this extensive process, challenges such as poor vaccine storage often arise. Such challenges are known to affect the validity period of the vaccines. In this work, we proposed and analyzed a timedelayed footandmouth disease model that incorporates vaccine waning and footandmouth disease carriers. We are aware that the role of FMD carriers in disease transmission remains a debatable issue [17, 22] and as such we have explored the outcomes of the model in two scenarios: (i) no transmission by carriers and (ii) when carriers are transmitting the disease. We computed the basic reproduction number and demonstrated that it is an important threshold parameter for the existence and global stability of the model steady states. Utilising numerical illustrations, we demonstrated the impact of footandmouth disease carriers on long and shortterm dynamics of the disease.
This work demonstrated that understanding the role of vaccine waning and FMD carriers is an important topic that requires multidimensional research. Nevertheless, this work is not exhaustive; it can be extended by adding diffusion terms to the proposed model so as to represent the movement and dispersal of the animals and the pathogen and by adding convection terms to represent the migration of animals.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
References
 S. Sieng and J. Kerr, āInvestigation of vaccination effectiveness in two Cambodian villages facing an outbreak of footandmouth disease,ā in Cattle: Cattle Health, Production and Trade in Cambodia, J. R. Young, L. Rast, S. Suon, and P. A. Windsor, Eds., pp. 67ā71, ACIAR Proceedings No.138 Australian Centre for International Agricultural Research, Canberra, Australia, 2013. View at: Google Scholar
 P. M. Schnell, Y. Shao, L. W. Pomeroy, J. H. Tien, M. Moritz, and R. Garabed, āModeling the role of carrier and mobile hers on footandmouth disease virus endemicty in the far north region of Cameroon,ā Epidemics, vol. 29, Article ID 100355, 2019. View at: Publisher Site  Google Scholar
 C. Bravo de Rueda, A. Dekker, P. L. EblĆ©, and M. C. M. De Jong, āVaccination of cattle only is sufficient to stop FMDV transmission in mixed populations of sheep and cattle,ā Epidemiology and Infection, vol. 143, no. 11, pp. 2279ā2286, 2015. View at: Publisher Site  Google Scholar
 Y. Cao, Z. Lu, and Z. Liu, āFootandmouth disease: progress and problems,ā Expert Review of Vaccines, vol. 15, no. 6, 2016. View at: Publisher Site  Google Scholar
 S. Sieng, S. W. WalkdenBrown, and J. Kerr, āEffect of vaccine storage temperatures and dose rate on antibody responses to foot and mouth disease vaccination in Cambodia,ā Veterinary Medicine and Science, vol. 4, no. 1, pp. 35ā44, 2018. View at: Publisher Site  Google Scholar
 T. B. Gashirai, S. D. MusekwaHove, P. O. Lolika, and S. Mushayabasa, āGlobal stability and optimal control analysis of a footandmouth disease model with vaccine failure and environmental transmission,ā Chaos, Solitons & Fractals, vol. 132, Article ID 109568, 2020. View at: Publisher Site  Google Scholar
 M. Rweyemamu, F. Maree, C. Kasanga et al., āChallenges and prospects for the control of footandmouth disease: an African perspective,ā Veterinary Medicine: Research and Reports, vol. 5, pp. 119ā138, 2014. View at: Publisher Site  Google Scholar
 N. Dexter, āStochastic models of foot and mouth disease in feral pigs in the Australian semiarid rangelands,ā Journal of Applied Ecology, vol. 40, no. 2, pp. 293ā306, 2003. View at: Publisher Site  Google Scholar
 C. Bravo de Rueda, M. C. de Jong, P. L. EblĆ©, and A. Dekker, āQuantification of transmission of footandmouth disease virus caused by an environment contaminated with secretions and excretions from infected calves,ā Veterinary Research, vol. 46, no. 1, p. 43, 2015. View at: Publisher Site  Google Scholar
 M. E. ChaseTopping, I. Handel, B. M. Bankowski et al., āUnderstanding footandmouth disease virus transmission biology: identification of the indicators of infectiousness,ā Veterinary Research, vol. 44, no. 1, p. 46, 2013. View at: Publisher Site  Google Scholar
 D. M. Green, I. Z. Kiss, and R. R. Kao, āModelling the initial spread of footandmouth disease through animal movements,ā Proceedings of the Royal Society B: Biological Sciences, vol. 273, no. 1602, pp. 2729ā2735, 2006. View at: Publisher Site  Google Scholar
 R. R. Kao, L. Danon, D. M. Green, and I. Z. Kiss, āDemographic structure and pathogen dynamics on the network of livestock movements in Great Britain,ā Proceedings of the Royal Society B: Biological Sciences, vol. 273, no. 1597, pp. 1999ā2007, 2006. View at: Publisher Site  Google Scholar
 M. J. Keeling, M. E. J. Woolhouse, R. M. May, G. Davies, and B. T. Grenfell, āModelling vaccination strategies against footandmouth disease,ā Nature, vol. 421, no. 6919, pp. 136ā142, 2003. View at: Publisher Site  Google Scholar
 R. P. Kitching, āIdentification of foot and mouth disease virus carrier and subclinically infected animals and differentiation from vaccinated animals,ā Revue Scientifique et Technique de lāOIE, vol. 21, no. 3, pp. 531ā538, 2002. View at: Publisher Site  Google Scholar
 R. P. Kitching, A. M. Hutber, and M. V. Thrusfield, āA review of footandmouth disease with special consideration for the clinical and epidemiological factors relevant to predictive modelling of the disease,ā The Veterinary Journal, vol. 169, no. 2, pp. 197ā209, 2005. View at: Publisher Site  Google Scholar
 S. Mushayabasa, C. P. Bhunu, and M. Dhlamini, āImpact of vaccination and culling on controlling foot and mouth disease: a mathematical modelling approach,ā World Journal of Vaccines, vol. 1, no. 4, pp. 156ā161, 2011. View at: Publisher Site  Google Scholar
 S. Mushayabasa, D. Psosny, and J. Wang, āModeling the intrinsic dynamics of FMD,ā Molecular Biology and Evolution, vol. 13, pp. 425ā442, 2015. View at: Google Scholar
 S. Mushayabasa and G. Tapedzesa, āModeling the effects of multiple intervention strategies on controlling footandmouth disease,ā BioMed Research International, vol. 2015, Article ID 584234, 10 pages, 2015. View at: Publisher Site  Google Scholar
 N. Ringa and C. T. Bauch, āImpacts of constrained culling and vaccination on control of foot and mouth disease in nearendemic settings: a pair approximation model,ā Epidemics, vol. 9, pp. 18ā30, 2014. View at: Publisher Site  Google Scholar
 G. Rossi, R. Smith, S. Pongolini, and L. Bolzoni, āModelling farmtofarm disease transmission through personnel movements: from visits to contacts, and back,ā Scientific Reports, vol. 7, no. 1, p. 2375, 2017. View at: Publisher Site  Google Scholar
 J. Zhang, Z. Jin, and Y. Yuan, āAssessing the spread of foot and mouth disease in mainland China by dynamical switching model,ā Journal of Theoretical Biology, vol. 460, pp. 209ā219, 2019. View at: Publisher Site  Google Scholar
 A. B. R. Parthiban, M. Mahapatra, S. Gubbins, and S. Parida, āVirus excretion from footand mouth disease virus carrier cattle and their potential role in causing new outbreaks,ā PLoS One, vol. 10, no. 6, Article ID e0128815, 2015. View at: Publisher Site  Google Scholar
 J. K. Hale and W. Z. Huang, āGlobal geometry of the stable regions for two delay differential equations,ā Journal of Mathematical Analysis and Applications, vol. 178, no. 2, pp. 344ā362, 1993. View at: Publisher Site  Google Scholar
 S. Alexandersen, M. Quan, C. Murphy, J. Knight, and Z. Zhang, āStudies of quantitative parameters of virus excretion and transmission in pigs and cattle experimentally infected with footandmouth disease virus,ā Journal of Comparative Pathology, vol. 129, no. 4, pp. 268ā282, 2003. View at: Publisher Site  Google Scholar
 S. Parida, L. Fleming, Y. Oh et al., āReduction of footandmouth disease (FMD) virus load in nasal excretions, saliva and exhaled air of vaccinated pigs following direct contact challenge,ā Vaccine, vol. 25, no. 45, pp. 7806ā7817, 2007. View at: Publisher Site  Google Scholar
 C. Stenfeldt, J. M. Pacheco, L. L. Rodriguez, and J. Arzt, āEarly events in the pathogenesis of footandmouth disease in pigs; identification of oropharyngeal tonsils as sites of primary and sustained viral replication,ā PLoS One, vol. 9, no. 9, Article ID e106859, 2014. View at: Publisher Site  Google Scholar
 J. Slingluff, āRisk assessment for the transmission of footandmouth disease via movement of swine and cattle carcases from infected FMDinfected premises to a disposal site,ā 2014, https://conservancy.umn.edu/handle/11299/193839. View at: Google Scholar
 V. Lakshmikantham, S. Leela, and A. A. Martynyuk, Stability Analysis of Nonlinear Systems, Marcel Dekker, New York, NY, USA, 1989.
 J. Huo, H. Zhao, and L. Zhu, āThe effect of vaccines on backward bifurcation in a fractional order HIV model,ā Nonlinear Analysis: Real World Applications, vol. 26, pp. 289ā305, 2015. View at: Publisher Site  Google Scholar
 J. S. LaSalle, āThe stability of dynamical systems,ā in CBMSNSF Regional Conference Series in Applied Mathematics, vol. 25, SIAM, Philadelphia, PA, USA, 1976. View at: Google Scholar
 J. K. Hale and S. Verduyn Lunel, Introduction to Functional Differential Equations, Springer, New York, NY, USA, 1993.
 L. F. Shampine and S. Thompson, āSolving DDEs in MATLAB,ā Applied Numerical Mathematics, vol. 37, no. 4, pp. 441ā458, 2001. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Tinashe B. Gashirai et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.