Modeling the potential of wAu-Wolbachia strain invasion in mosquitoes to control Aedes-borne arboviral infections

Arboviral infections such as dengue, Zika and chikungunya are fast spreading diseases that pose significant health problems globally. In order to control these infections, an intracellular bacterium called Wolbachia has been introduced into wild-type mosquito populations in the hopes of replacing the vector transmitting agent, Aedes aegypti with one that is incapable of transmission. In this study, we developed a Wolbachia transmission model for the novel wAu strain which possesses several favourable traits (e.g., enhanced viral blockage and maintenance at higher temperature) but not cyctoplasmic incompatibility (CI)—when a Wolbachia-infected male mosquito mates with an uninfected female mosquito, producing no viable offspring. This model describes the competitive dynamics between wAu-Wolbachia-infected and uninfected mosquitoes and the role of imperfect maternal transmission. By analysing the system via computing the basic reproduction number(s) and stability properties, the potential of the wAu strain as a viable strategy to control arboviral infections is established. The results of this work show that enhanced maintenance of Wolbachia infection at higher temperatures can overcome the lack of CI induction to support wAu-Wolbachia infected mosquito invasion. This study will support future arboviral control programs, that rely on the introduction of new Wolbachia variants.

To control these infections, an intracellular bacterium called Wolbachia can be used to suppress transmission in arthropods such as mosquitoes and flies [22][23][24][25] . Wolbachia infection inhibits arboviral transmission in mosquitoes via four mechanisms: immune priming-preactivation of the mosquito immune system; induction of the phenoloxidase cascade-triggers immune response to viruses; competition of intracellular resourcesinducing authophagy; and induction of microRNA-dependent immune pathways-essential for gene regulation and stability, immune defense, ageing and organ differentiation 26 . This endosymbiotic bacterium which exists naturally in more than 50% of all insect species can be found within the cytoplasm of the cells of their hosts 25,27,28 . Whilst Wolbachia is not naturally present in Aedes aegypti, it can be introduced via stable transinfections using microinjections 29,30 .
The Wolbachia-based control strategy is carried out by infecting mosquitoes with a strain of Wolbachia and then releasing them into wild mosquito populations in the hopes of replacing the vector transmitting agent Aedes aegypti with one that is incapable of transmission [29][30][31] . Infecting an Aedes mosquito with Wolbachia can change some of the Aedes characteristic features. In practice, Wolbachia can reduce the life-span of mosquitoes by half producing a deleterous fitness effect 32 . Another feature is cytoplasmic incompatibility (CI) 22,[33][34][35] which occurs when a Wolbachia infected male mates with an incompatible female mosquito (usually Wolbachia uninfected) producing no offspring 36 . Other features of Wolbachia which serve as liabilities in mosquitoes include: imperfect maternal transmission (IMT) 30,37 and loss of Wolbachia infection (LWI). LWI impedes the establishment of Wolbachia-infected mosquitoes and is a result of mosquito vulnerability to high temperature 38,39 .
However, a novel strain of Wolbachia: wAu, has shown to produce high viral blockage whilst maintaining Wolbachia infection in Aedes mosquitoes at higher temperature 23 . Moreover, wAu allows superinfection to occur when wAu and other strains of Wolbachia co-exist in the vector host 23 . Despite these favourable features, wAu does not induce CI 23 . Although CI absence does not establish Wolbachia infected mosquitoes, the effect could be outweighed by LWI and IMT 37 .
The difference in the common Wolbachia strain features are described in Table 1 below. In general, the introduction of mathematical models to understand infection dynamics of diseases has long been helpful in the area of disease control 49 . A number mathematical models of Wolbachia dynamics in a mosquito population have been formulated 37,[50][51][52][53][54][55][56][57][58] . Some of these models introduced Wolbachia strain(s) into a mosquito population and classified them into age-sturctured Wolbachia-infected and -uninfected mosquito compartments 37,53,54,57 . Ndii et al. 53 , formulated a mathematical model for the Wolbachia interaction between the immature stages (aquatic stage), adult male and female mosquito populations to investigate the persistence of mosquitoes infected with Wolbachia when competing with the uninfected ones. They derived the steady state solutions and showed that parameters such as maternal transmission, reproductive, death and maturation rates drive the persistence of the Wolbachia-infected mosquito population. A similar model developed by Xue et al. considered the Wolbachia-induced fitness change and the CI effect 57 . They showed that if the basic reproduction number ( R 0 ) of the Wolbachia-infected mosquitoes is less than one, an endemic Wolbachia infection can still occur via backward bifurcation if a sufficient number of the mosquitoes are introduced into the population. A mathematical model of Wolbachia to control dengue fever transmission 52 was developed by Hughes et al. The model showed that the use of Wolbachia has high potential to control dengue where the R 0 due to Wolbachiainfected Aedes mosquitoes is not too large in endemic areas. Another study of a Wolbachia invasive model incorporated IMT and LWI and showed that CI does not guarantee the establishment of Wolbachia-infected mosquitoes as the disadvantages derived from IMT and LWI in the production of Wolbachia-infected mosquitoes could outweigh CI 37 .
Additionally, a study conducted by O'Reilly et al combining multiple modeling methods, was used to estimate the burden of dengue and map its distribution across Indonesia 59 . They predicted that there was a reduction in dengue transmission after a nationwide release of wMel-Wolbachia-infected mosquitoes. In addition, they predicted about 86% of the estimated 7.8 million annual cases of symptomatic dengue in Indonesia could be averted following a complete nationwide rollout of Wolbachia-infected mosquitoes. Recently, a modeling study presented a dengue transmission model in the presence of female wild-type and wMelPop Wolbachia-infected Aedes aegypti mosquitoes. They concluded that although the wMelPop strain reduces the lifespan of infected mosquitoes, which could be challenging to achieve replacement of wild-type mosquitoes, its optimal release ensured the replacement of wild-type mosquitoes and also reduced dengue burden in the human population 51  www.nature.com/scientificreports/ strains (wMel and wAlbB) to control arboviral spread 60 . They observed that each of the two different strains of Wolbachia can effectively decrease the rate of arboviral transmission.
Here, we develop a general Wolbachia model capable of faithfully replicating all of the strain features described in Table 1. The general transmission model is an extention of the Wolbachia transmission model introduced in Adekunle et al. 37 , which described the competitive dynamics between (wMel-like) Wolbachia-infected and uninfected mosquitoes. Despite the non-induction of CI in wAu-Wolbachia-infected mosquitoes, wAu infection is retained and able to block viral transmission efficiently compared to other strains even at high temperature. Therefore, we incorporated this feature to determine if the advantages (Wolbachia retainment) of the wAu strain outweigh the ineffectiveness of CI. This feature has not been considered in previous models. Furthermore, we incorporate imperfect maternal transmission into the model. By analysing the system via computing the basic reproduction number(s) and investigating the stability properties of the equilibrium points, the potential of the wAu strain as a viable strategy to control Aedes-borne infections can be established. The aim of this modeling approach is to support future Aedes-borne viral control programs, particularly with the introduction of new Wolbachia variants.

Methods
Model formation. Here, we investigate a modified Wolbachia transmission model studied in Adekunle et al. 37 , focusing on a novel Wolbachia strain, wAu, which has high retainment, high viral blockage and does not induce CI. The mosquito population is subdivided into two groups: the uninfected mosquitoes (.) u and the Wolbachia infected mosquitoes (.) w . The term (.) can be aquatic/immature (eggs, larvae and pupae) A, male M or female F mosquitoes. In addition, we denote the aquatic/immature stages, mature male and mature female uninfected mosquitoes as A u , M u , F u , and Wolbachia-infected mosquitoes as A w , M w , F w respectively. As in Adekunle et al. 37 the model also incorporates the IMT of wAu-Wolbachia.
There are four possible mosquitoes' mating pairs: F u M u , F u M w , F w M u and F w M w . As Wolbachia infection is maternally transmitted, F u M u and F u M w will produce uninfected offspring while F w M u and F w M w will typically produce infected offspring. However if there is imperfect maternal transmission, the two latter strategies could produce some proportions of uninfected offspring 23 .
To mathematically write the system of differential equations governing the Wolbachia transmission dynamics, we express the feasible mating strategies of uninfected and Wolbachia infected mosquito populations together with their per capita egg laying rates as Eqs. (1)-(6): Here, φ represents the CI effect which can be either 0 if there is no CI, or 1 if CI is present. σ is the effect of LWI, such that it can either be 0, if there is no Wolbachia loss or greater than zero otherwise. In Adekunle et al. 37 where CI is assumed and LWI is considered, these quantities are set to φ = 1 and σ ≥ 0 . In our modified model, considering different strains with the exception of wAu strain, φ = 1 and σ could vary from values greater than zero onwards. However, for the wAu-Wolbachia strain, CI is ineffective and high retainment of wAu-Wolbachia infection even at high temperatures 23 is established, therefore we set φ = 0 and σ = 0 . Our model also incorporates imperfect maternal transmission generating a proportion of infected and uninfected offspring from mating of both F w M u and F w M w mosquitoes. To simplify the system, we assume that M = F in accordance with the observed ratio of male to female mosquitoes of 1.02:1 62 . That is, we set ψ = 1/2 (Fig. 1). By this, it follows that the system of ordinary differential equations (ODEs) in Eqs. (1)- (6) can be reduced to (7)-(10) which is the governing Wolbachia infection dynamics.
To mathematically express the above schematics, we have that, the feasible mating strategies of uninfected and Wolbachia infected mosquito populations together with their per capita egg laying rates are given by the following differential system:  www.nature.com/scientificreports/ The closed set which is a feasible region for the above system dynamics is positively invariant 37 .
Also, letting ν = 1 , we assume a perfect maternal transmission for the reproduction outcome of F wMw mating. Therefore, the general Wolbachia model in terms of population proportions is given by Eqs. (12)- (14). Hereafter it is clear that we refer to the scaled values of each state variable and as such drop the overbar from our notation. The scaled model below now evolves in the feasible region ¯ , where The modeling of wAu-Wolbachia transmission dynamics has not been done as this a distinction from other Wolbachia transmission models. Unlike the modeling work in Adekunle et al. 37 , apart from the non-induction of CI, we considered the loss of Wolbachia infections due to seasonal fluctuation in temperature, a key dynamics that is absent in wAu strain.

Results
Analysis of the model. The above general model (11)- (14) is parametrically adjusted to simultaneously accommodate wAu and wMel Wolbachia strains. For the wAu-Wolbachia model, we set φ = σ = 0 and for the wMel-Wolbachia model, we set φ = 1 , σ > 0 . The wMel-Wolbachia model parameter adjustments correspond to the model studied in Adekunle et al. 37 .
Here, we want to analyse the general model (11) is the equilibrium point for the co-existence of both uninfected and wAu-Wolbachia-infected mosquito populations.
Non-existence mosquito population, e 1 . The equilibrium point e 1 is trivial and is not biologically realistic. However, we can gain some insights into the competitive model dynamics by examining the case where there is no interaction between the uninfected and Wolbachia-infected mosquitoes. In other words, we want to investigate how each population would behave in the absence of the other. In particular, we derive the reproduction number of the uninfected R 0u and Wolbachia-infected R 0w mosquito populations when they do not interact: where the factor of 1 2 in R 0u and R 0w stems from the choice to set M = F 62 , i.e. ψ = 1 2 . These reproductive numbers determine if the uninfected and Wolbachia-infected mosquito populations will die out or persist when there is no interaction. Specifically, if R 0u < 1 and R 0w < 1 , then the two populations will die out (Fig. 2a). We observed in the decoupled case, the expressions for R 0u and R 0w are independent of the effects of CI (φ) and LWI (σ ) and are therefore equivalent for both the wAu and wMel-Wolbachia strains (Fig. 2) 37 .
Uninfected mosquito population, e 2 . The uninfected-mosquito-only equilibrium point or Wolbachia-free equilibrium is For e 2 to exist, we require R 0u > 1 . In addition to the uncoupled reproduction numbers ( R 0u and R 0w ) we also define the invasive reproduction number R 0w|u which describes the average number of secondary offspring that , where we have substituted in the definition of R 0w from Eq. (16). The invasive reproduction number R 0w|u is the same for both wAu and wMel-Wolbachia strains as that derived in Adekunle et al. 37 . This is because, the expression (17) clearly shows that the invasive reproductive number R 0w|u is not dependent on the CI effect, φ or LWI, σ.
To check if the equilibrium point e 2 is stable, we compute the Jacobian of the system and evaluate it at e 2 . In particular, letting z 1 = (µ Au + τ u ) and z 2 = (µ Aw + τ w ) , yields To obtain the characteristic equation of J e 2 , we have which becomes where Therefore, e 2 is locally asymptotically stable if and only if R 0w|u < 1 and R 0u > 1 (Fig. 4). This is also consistent with the study in Adekunle et al. 37 (See Table 3).
Wolbachia-infected mosquito population, e 3 . The wAu-infected-only equilibrium point is This again is consistent with Adekunle et al. 37 . For e 3 to exist we require R 0w > 1 . By computation, the invasive reproductive number R 0u|w with respect to uninfected mosquitoes is given as, which is equivalent to that of Adekunle et al. 37 . However, for the wAu-Wolbachia strain, i.e. φ = 0 , |J e 2 − I| = 0, The red ball point indicates the point of stability, that is (F u , F w ) = (0, 0) representing mosquito extinction. We set ρ uu = 0.01 and ρ ww = 0.1 . Other parameters used for these model simulations are provided in Table 2. www.nature.com/scientificreports/ we have a modified expression of c = 1 + ρ ww ρ uu (1 − δ) in Eq. (18) because we do not assume CI. Therefore, c ≥ 1 for wAu-Wolbachia strain. Computing the Jacobian at e 3 , we have: The characteristic equation of J e 3 is then where Therefore, e 3 is locally asymptotically stable if and only if R 0u|w < 1 and R 0w > 1 (see Fig. 4). The condition is equivalent to that found in 37 with generalized expressions for R 0u|w used in place of the reduced version presented there (see Table 3).
The above conditions (i) and (ii) correspond to the cases for δ > 1 c and δ < 1 c respectively. Comparing these existence conditions with those found above for e 2 and e 3 , we see that condition (ii) for the existence of e 4 matches the combined existence and local asymptotic stability condition for e 2 and e 3 . In other words, e 2 , e 3 and e 4 can coexist, while e 1 always exists (see Fig. 4).
To establish whether e 4 is stable or not, we compute the Jacobian J e 4 evaluated at e 4 to obtain the following characteristic equation: then we have: |J e 4 − I| := 4 + s 1 3 + s 2 2 + s 3 + s 4 = 0. www.nature.com/scientificreports/ In order to establish the nature of the equilibrium point e 4 , we performed numerical testing using the Monte Carlo method in 50 to verify the conditions (i) and (ii) by computing the real part of the eigenvalues of the Jacobian matrix, evaluated at e 4 . Simulation results are illustrated in Fig. 3.
Although the conditions (i) and (ii) indicated the existence of e 4 , Fig. 3c showed that e 4 is locally stable for condition (i) as all the eigenvalues (real part) are negative ( 1 , 2 , 3 , 4 < 0 ). Whilst Fig. 3f showed that e 4 is unstable for condition (ii) as two of the eigenvalues (real part) are positive i.e. 3 , 4 > 0.
Numerically, we illustrated the existence and stability regions for e 4 in Fig. 4

for the two conditions (i) and (ii) relating to CI and maternal transmission (MT).
Following a modeling study of Aedes aegypti mosquitoes and normal Wolbachia (in the presence of CI only) interaction analyzed by Ferreira et al. 64 , three equilibrium points: trivial ( q 1 ); uninfected only ( q 2 ); and coexistence ( q 3 ), were obtained. However, the Wolbachia-only equilibrium point was not computed. The established local stability conditions for q 1 and q 2 correspond to that of the wMel-like Wolbachia conditions for e 1 and e 2 respectively. For coexistent populations to persist, the reproductive number for infected mosquitoes only, R i must be greater than 1 and R i > R u , where R u is the reproductive number for wild-type mosquitoes only. The model 64 also described the fitness parameter space between R u and R i , showing the change in extinction and persistence of the three equilibria when there is an increase in the initial population proportion of the Wolbachiainfected mosquitoes. Our model showed the changes in the no-mosquito, wild-type only, Wolbachia-only and coexistence population persistence and extinction in the presence and absence of CI with high and low maternal transmission (MT). Figure 4 illustrates the existence and local stability regions for the equilibrium points e 1 , e 2 , e 3 and e 4 with respect to the reproduction numbers R 0u and R 0w as well as the relative magnitude of δ and 1 c . For δ > 1 c (high MT), Fig. 4a,b describe the dynamics for φ = 0 (CI absent) and φ = 1 (CI present) respectively. Within the subset of the yellow region of these figures bounded by R 0u = 1, R 0w = 1 , and R 0w = ξ we find that only e 1 and e 3 exist. Since e 3 is unstable in this region, we expect the system trajectories to tend to the no-mosquito equilibrium e 1 . This was confirmed through numerical simulations shown in Fig. 5a. For the existence of e 4 we require R 0u|w > 1 , R 0w|u > 1 and R 0w > ξ for stability (within the blue region). But if R 0w < ξ , e 1 is stable (yellow).
For δ < 1 c (low MT), Fig. 4c,d portrayed the regions of stability for φ = 0 and φ = 1 respectively. The conditions R 0u < 1, R 0w > 1 , and R 0u|w > 1 project the trajectiory to tend to e 1 (see Fig. 5b). In the orange region, e 1 and e 3 exist and are simultaneously locally stable as R 0u|w > 1 and R 0w > ξ . In addition, we have that e 4 exists where R 0u|w < 1 and R 0w|u < 1 (condition (ii)). With these conditions, e 4 exists together with e 2 and e 3 (white region). In this white region, e 2 and e 3 are locally stable even as R 0u > 1, R 0w > 1 but e 4 is unstable. Also, e 1 exists when R 0w > 1 and R 0u < 1 because the local stability of other equilibrium points is violated with these conditions. When R 0u|w > 1 but R 0u < 1 and R 0w > 1 , the only stable outcome is the mosquito-free (no-mosquito) equilibrium e 1 . This occurs when R 0u is less than but still close to one. In this region, uninfected mosquitoes are capable of dominating initially when introduced into a Wolbachia saturated equilibrium because imperfect maternal transmission achieves R 0u|w > 1 . This competitive advantage drives out the Wolbachia infected mosquitoes leaving uninfected mosquitoes only, which then are unable to sustain their population because R 0u < 1 (Fig. 5).
With the rate of high maternal transmission (MT) in the absence of CI (like-wAu), the reproductive advantage favours the production of uninfected mosquito offspring as it tends to accommodate more coexistent mosquito populations with wild-type than wMel-like strain (presence of CI) due to the presence of CI (Fig. 4a,b). Whilst, with a low MT rate, the CI presence or absence would favour Wolbachia-infected mosquitoes or uninfected mosquitoes respectively. In other words, the coexistent equilibrium point is unstable for the two mosquito populations as these conditions are equivalent to the local stabilities of both Wolbachia-free and Wolbachia-only equilibrium points (Fig. 4c,d). If R 0w < ξ , the system trajectories tend to the no mosquito equilibrium e 1 .
The conditions for the local stability of all equilibrium points are shown in Table 3 below.
Sensitivity analysis of Wolbachia model. To carry out the sensitivity analysis we investigate the model robustness due to uncertainties associated with parameter value estimations. In other words, we examine how senitive the invasive reproductive numbers are with respect to these parameters. This in turn, gives insight on influential parameters and their impact in reducing (or increasing) mosquito-type populations. To carry out this, we compute the normalized sensitivity indices of the invasive reproduction numbers with respect to the parameters used in the model.
Definition. The normalized forward sensitivity index of a variable v with respect to parameter w is defined as: Using the above formular (20), we contruct the following plots in Fig. 6. From Fig. 6 and using the baseline parameter values for the wAu-Wolbachia strain in Table 2, it is clear that the reproductive and mortality rates for both wild-type ( ρ uu , µ u ) and wAu-Wolbachia-infected ( ρ ww , µ w ) mosquitoes and the proportion of wAu-Wolbachia-infected offspring ( δ ) have the most sensitivity in the invasive reproductive numbers R 0w|u . Whilst for R 0u|w , µ u and µ w are the most sensitive parameters. Hence for both invasive reproductive numbers, the most sensitive parameters are µ u and µ w . This demonstrates that an increase (or decrease) in the mortality rate of wAu-Wolbachia-infected mosquitoes by 10% will decrease (or increase) R 0w|u by 10%.
Does CI (φ) outweigh the LWi (σ )? For most Wolbachia strains except wAu, the mating between uninfected female and Wolbachia-infected male mosquito crosses generates no viable offspring. However, Wolbachiainfected mosquitoes tend to lose their Wolbachia infection and lower their maternal transmission rate at high  23 . With the effect of climate change gradually increasing the temperature by the day, Wolbachia strains with moderate or high temperature sensitivity such as wMel may not be able to fully maintain a sufficient frequency level to invade the mosquito population. In our general Wolbachia mathematical model, we describe a modified version of Adekunle et al. 37 . This modification accommodates parameter adjustments for novel wAu and wMel-Wolbachia strains. For wAu, our mathematical model showed that despite the production of mosquito offspring due to CI absence, the invasive reproduction number due to infected mosquitoes R 0w|u remains unchanged compared to the case where CI is present, as with the wMel-like strain 37 . This further strengthened the fact that CI (inclusion or exclusion) does  www.nature.com/scientificreports/ not guarantee Wolbachia mosquitoes' persistence. Also, the invasive reproduction number due to uninfected mosquitoes expression R 0u|w for wAu is similar to wMel, except that the expression depends on CI the effect. This is because, the mosquito gender crosses due to non-induction of CI for wAu, i.e. F u M w , generates uninfected offspring with perfect maternal transmission while wMel does not. The chances of establishing Wolbachia infected mosquitoes are lower when CI is ineffective compared to when it is induced. That is, for cytoplasmic inducing wMel-Wolbachia mosquitoes, the effect of LWI outweighs CI effect as mosquitoes still lose their infections (Fig. 7). However, wAu-Wolbachia infection retainment (no LWI) in mosquitoes has shown high level of maintaining the Wolbachia frequency in the absence of CI in mosquitoes (Fig. 7). This suggests that the LWI effect outweighs CI. The LWI rate σ (t) which is dependent on the seasons of the year can be modeled by a sinusoidal equation: The yellow shaded region indicates the local stability of e 1 equilibrium. The green shaded area illustrates the local stability for the Wolbachia-free equilibrium point ( e 2 ). e 3 is locally stable at the red shaded part. The blue region indicates the coexistence local stability e 4 . The white region shows the existence of e 2 , e 3 and e 4 and local stability of e 2 and e 3 equilibrium points. And the orange region describes the existence and local stability of e 1 and e 3 . For δ > 1 c ; (a) describes φ = 0 as the boundary R 0w|u = 1 sits above the boundary R 0u|w = 1 and the arc R 0w = ξ . The co-existent equilibrium e 4 (blue), always sits in the region between these three boundaries because R 0w|u > 1 , R 0u|w > 1 and R 0w > ξ . If R 0w < ξ , then e 1 becomes stable (yellow). (b) describes similar conditions as in (a) but for φ = 1 . We observed that the boundary R 0u|w = 1 shifts up while R 0w|u = 1 remained stationary to accommodate more e 3 . For δ < 1 c ; (c) describes φ = 0 as the relative position of e 4 boundaries in (a) flips so that boundary R 0u|w = 1 sits above boundary R 0w|u = 1 and the arc R 0w = ξ . Then, R 0w|u < 1 and R 0u|w < 1 and R 0w > ξ shows the co-existence of e 2 and e 3 (white). However, e 2 and e 3 are locally stable in the white region as R 0w > 1 and R 01 > 1 . For R 0u < 1 , e 2 and e 4 do not exist, only e 1 and e 3 do and if R 0w > ξ , e 1 and e 3 are locally stable (orange) and if R 0w < ξ , only e 3 becomes stable (red). (d) describes similar conditions as in (c) but for φ = 1 . It was observed that the boundary R 0u|w = 1 shifts up reducing the region of stability for e 2 .   Table 3. Expressions for the condition for stability associated with the equilibrium points.

Equilibrium points
Conditions for stability wMel 37 wAu www.nature.com/scientificreports/ The effects of CI ( φ ) and LWI ( σ (t) ) as features of wAu and wMel Wolbachia strains are shown in Fig. 7. For the total mosquito population, wAu-infected mosquitoes ( φ = 0, σ max = 0 ) reach the maximum frequency after approximately 250 days. To see the effect of CI induction and slight LWI i.e. φ = 1 , C = 0.25 , for σ max = 0.02 and σ max = 0.04 , the Wolbachia frequency level oscillates between (0.8 and 1) and (0.6 and 1) respectively. That is, there is a 20% and 40% drop in the frequency level of Wolbachia when σ (t) is at σ max = 0.02 and σ max = 0.04 respectively. This showed that, despite CI induction, LWI reduced the contribution of CI to the Wolbachia invasion (Fig. 7a). Therefore, the LWI gains highly outweigh the CI effect. By this, our analysis suggests that an increase in LWI in the presence of CI results in a drastic decrease in the Wolbachia frequency level (Fig. 7a). On the other hand, Fig. 7b showed the effect of LWI σ (t) and CI φ with respect to the competitiveness between F u and F w . We observed that the F w population dominates the F u when there was no CI induction and Wolbachia infection is retained, that is, φ = 0, σ max = 0 (Fig. 7b). However, if CI induction occurs with loss of Wolbachia infections, then the seasonal varying effect occurs as seen in Fig. 7c.

Discussion
In this work, we modelled and investigated a general Wolbachia model that contained the transmission dynamics of wAu and wMel Wolbachia strains in Aedes mosquitoes as special cases. These transmission dynamics described the competition between the novel wAu-Wolbachia infected Aedes mosquitoes and wild-type mosquitoes and compared the dynamics with the invasive properties of the popular wMel-Wolbachia infected mosquitoes. We first derived the Wolbachia infection-status reproduction numbers for our wAu-Wolbachia model and used them to establish the conditions for the local stability of the equilibrium points for the wAu-Wolbachia invasive model. The reproduction number associated with the uninfected mosquitoes shows the reproductive advantage that the wild type has over the wAu strain. The comparison of the wAu-Wolbachia model (CI and LWI absent) and wMel-Wolbachia model (CI and LWI present) showed that the wAu strain has the potential of compensating for the undesirable features of the wMel strain. Additionally, this study has reviewed the main features of different Wolbachia strains ( Table 1) and shown that the wAu Wolbachia strain is a promising candidate for efficient Aedes-borne arboviral transmission control. Moreover, we analyzed the system dynamics of a general Wolbachia invasion model and determined the regions of local stability for each of the identified equilibrium points, highlighting the regions in parameter for which Wolbachia-infected mosquito populations persist or go extinct. This work modelled the general Wolbachia dynamics which can accommodate various Wolbachia characteristics regarding the presence or absence of CI and seasonal changes, unlike Adekunle et al. 37 , which considers only the presence of CI. We also investigated the advantages gained from CI and LWI. This study has demonstrated that despite the absence of CI, the Wolbachia frequency level will drop as much as tenfold of the percentage of Wolbachia infection lost. We showed that the advantage of Wolbachia retainment in mosquitoes strongly outweighed the negative impact of CI indicating wAu Wolbachia strains may be suitable for arboviral control. Therefore, this modeling work contributes to the previous studies 37,54,57,64,65 and helps close the gap between ways of maintaining the Wolbachia frequency levels in the absence of LWI and CI.
One implementation question for using the wAu strain as a replacement of the wMel strain is whether the wAu strain is self-sustaining, given that it does not induce CI. In this work, the equilibrium points for the wAu-Wolbachia model are the same as that for the wMel-Wolbachia model except that stricter conditions are required to satisfy the wAu-Wolbachia model equilibrium points. These more stringent conditions translate to additional resources such as the continuous introduction of a larger scale of wAu-infected mosquitoes to ensure replacement 66 . Thus, the wAu strain is a promising alternative strain as it does not suffer from LWI due to high weather temperature and is highly effective in preventing the transmission of the arbovirus 23,39,67 . Otherwise, combining the two strains may also be a good strategy.
There are limitations associated with any mathematical modeling work, and this study is not exempted. We first assumed the same mosquito gender ratio and expected this proportion to be constant over time. This assumption may be true in a laboratory setting 62 , but not necessarily true in a natural mosquito habitat. However, similar conclusions are expected to be reached as the Wolbachia model reduction accurately reproduces the dynamics of the full system 68 . Secondly, we assumed that the absence of CI implies that cross mating resulted in offspring that are uninfected. This may not be true as a small proportion of the offspring may be Wolbachia infected 23 . If that is the case, then it means that lesser resources will be required to use the wAu strain as a Wolbachia-based control strategy. Lastly, we assumed the seasonality affects the associated parameters for the wMel dynamics. However, for the wAu strain, it is not affected by seasonality as wAu-Wolbachia infections are retained at high temperature.
Although several studies 22,38,42,57 have demonstrated that CI drives the persistence of Wolbachia-infected Aedes mosquitoes, these studies neglected the impact of Wolbachia loss in mosquitoes. The CI drive has been shown in four mating lines (see Fig. 1) involving a Wolbachia-transinfected Aedes mosqiutoes mating with wild-type mosquitoes. One of the mating lines for which Wolbachia-infected male and uninfected female mosquitoes produced no viable offspring (via CI) truncates the uninfected offspring from being produced as infection is maternally transmitted. With the exception of the mating between the uninfected male and female mosquito line, all other mating lines produce Wolbachia-infected offspring leading to persistence. In addition, high temperature affects these Wolbachia-infected mosquitoes as they lose their infection due to the unfavourable weather conditions. However, mosquitoes infected with the wAu-Wolbachia strain have been shown to not only block arboviral transmission efficiently, but also retain the Wolbachia infection at typically unfavourable high temperatures. This retainment of infection in mosquitoes strongly outweighed the absence of CI for the wAu strain in the establishment and dominance of wAu-Wolbachia infected mosquitoes.
While vaccine implementation may have been highly effective on dengue seropositive persons in high transmission areas 11,12 , the introduction of Wolbachia-infected mosquitoes in low and moderate arboviral endemic areas has also effectively shown successful reduction in dengue burden 43,51,59,69 . Given that these two strategies could reduce the transmission of Aedes-borne diseases, in particular, dengue depending on the transmission level, a modeling study by Ndii 70 proposed the use of these combined strategies and compared their effectiveness. The author showed that, Wolbachia performs better in the presence of low vaccine efficacy, but is outperformed otherwise 70 . Therefore combining the two strategies may be useful, however understanding both the temperature and seasonality effects on Wolbachia intervention programs, and serotypic differences relating to cross-protective immunity to investigate vaccine efficacy is necessary for the reduction and control of Aedes-borne arboviral disease transmission.
In conclusion, we have shown that the wAu-Wolbachia strain could be effective in controlling arbovirus transmission, as its advantages in terms of Wolbachia infection retention in mosquitoes may outweigh the absence of CI. This could prove even more promising, especially as the temperature increases due to climate