A co-infection model for Oncogenic HPV and TB with Optimal Control and Cost-Effectiveness Analysis ================================================================================================== * A. Omame * D. Okuonghae ## Abstract A co-infection model for oncogenic Human papillomavirus (HPV) and Tuberculosis (TB), with optimal control and cost-effectiveness analysis is studied and analyzed to assess the impact of controls against incident infection and against infection with HPV by TB infected individuals as well as optimal TB treatment in reducing the burden of the co-infection of the two diseases in a population. The co-infection model is shown to exhibit the dynamical property of backward bifurcation when the associated reproduction number is less than unity. Furthermore, it is shown that TB and HPV re-infection parameters (*ϕ*p = 0 and *σ*t = 0) as well as TB exogenous re-infection term (*ε*1 0) induced the phenomenon of backward bifurcation in the oncogenic HPV-TB co-infection model. The global asymptotic stability of the disease-free equilibrium of the co-infection model is also proven **not to exist**, when the associated reproduction number is below unity. The necessary conditions for the existence of optimal control and the optimality system for the co-infection model is established using the Pontryagin ‘s Maximum Principle. Uncertainty and global sensitivity analysis are also carried out to determine the top ranked parameters that drive the dynamics of the co-infection model, when the associated reproduction numbers as well as the infected populations are used as response functions. Numerical simulations of the optimal control model reveal that the intervention strategy which combines and implements control against HPV infection by TB infected individuals as well as TB treatment control for dually infected individuals is the most cost-effective of all the control strategies for the control and management of the burden of oncogenic HPV and TB co-infection. Keywords * Human Papillomavirus * tuberculosis * co-infection * optimal control * cost-effectiveness analysis ## 1 Introduction Tuberculosis (TB), caused by *Mycobacterium tuberculosis* (MTB), remains a global dominant cause of mortality and morbidity, although treatment and preventive measures are available [24, 38]. Globally, TB is one of the dominant cause of mortality from a single infectious disease agent, responsible for roughly 40% of deaths [36]. Worldwide, approximately, 10 million new cases of TB and about 1.57 million TB-induced deaths were reported in 2017, with incidence and mortality highest in South-Eastern Asia and Africa[36]. Estimates have equally shown that about one-fourth of the world ‘s population is infected with latent TB and is at the risk of progression to active TB [15]. Cancer is the major cause of mortality in developed countries in Europe and America and the second leading cause of death in developing countries [12]. Epidemiological evidences have proven that 80% of global cancer cases are attributed to the oncogenic Human papillomavirus (HPV) [2]. In a study carried out by Zetola *et al*. [37], it was revealed that prior TB infection is mostly common among patients with cases of cervical cancer. Likewise, Zhao *et al* [39], in another study, showed that TB infection is associated with high vulnerability to oncogenic HPV infection. According to Zetola *et al*. [37], oncogenic HPV and tuberculosis always co-exist, and the immune suppression caused by cancer (which is a consequence of oncogenic HPV) can result in latent TB re-activation, hence, leading to high mortality. In addition, persistent gynecologic TB infection, resulting in chronic inflammation, could be a high risk factor in the progression of oncogenic HPV infection to cervical cancer [37, 39]. Mathematical models have recently been studied to understand the transmission dynamics of human papillomavirus infection (See Omame et al. [23] and the references included therein). Malik *et al*. [17] investigated an optimal control model for HPV, incorporating optimal vaccination strategies for females in the population. In particular, they considered when three vaccines are used at the same time in comparison to the case where the bivalent *cervarix* and quadrivalent *Gardasil 4* vaccines were used at the initial stage and then, during the course of the vaccination program, one or two of them are interchanged with the nonavalent *Gardasil 9* vaccine. Saldana *et al*. [29] developed and analyzed a two-sex model of HPV with optimal control analysis, incorporating vaccination of adolescents, adults, and screening. They showed that to optimally curb the spread of HPV in a population, vaccination should be administered both before and after the commencement of active sexual life for both females and males in the population. Likewise, TB only models have been extensively studied in the literature (see for instance, [10, 21, 22, 33] and many others). Lately, mathematical models have been developed to consider the optimal control strategies for the dynamics of infectious diseases including their co-infections [1, 9, 11, 20, 26, 27, 31, 32, 34]. Agusto and Adekunle [1] studied the optimal control and cost-effectiveness analysis of the co-infection of drug-resistant tuberculosis and HIV/AIDS. They considered different control strategies for the control of the co-infections of both diseases. They showed that the strategy which combines efforts on minimizing the number of individuals with drug-sensitive and drug-resistant TB and case finding, prevention of treatment failure in the drug-sensitive TB infectious individuals and the treatment of individuals with drug-resistant TB is the most cost-effective in curbing the co-infections of drug-resistant TB and HIV/AIDS. Okosun and Makinde [20] considered a malaria-cholera co-infection model, in the presence of prevention and treatment controls for both diseases. Applying the Pontryagin ‘s Maximum principle, they derived the necessary conditions for optimality and equally showed using simulations, that to successfully curb malaria spread, control strategies must include cholera intervention strategies as well. Tilahun *et al*. [32] investigated a Pneumonia-Typhoid co-infection model with cost-effective optimal control analysis, incorporating preventive and treatment controls for both diseases. They showed that the strategy that combines Pneumonia treatment and Typhoid fever prevention is the most cost-effective in reducing the burden of the co-infections of Pnemonia and Typhoid fever. The researchers in [9] analyzed an optimal control model for HIV and TB co-infection, capturing resistance of HIV to anti-retroviral therapy. The authors showed that a combination of TB treatment and anti-retroviral therapy for HIV is the most effective in reducing the burden of the co-infection of the two diseases. In addition, Tanvi and Aggarwal [31] studied an HIV-TB co-infection, incorporating detection and treatment for both diseases. They observed that the strategy that implements detection for both diseases yielded the most economic and epidemic gains infighting the co-infections of the two diseases. To the best of the authors ‘knowledge, optimal control analysis has never been applied to the co-infection of Human papillomavirus and TB in the literature. Hence, it will be imperative to investigate the optimal control and cost-effectiveness analysis of the co-infections of oncogenic HPV and TB. In particular, this paper extends the work of Omame *et al*. [24] *by* assessing the impact of HPV prevention and TB treatment controls on the control and management of the co-infections of the two diseases. The rest of the paper is organized as follows. The model is formulated in Section 2, alongside the basic properties of the model. The co-infection model without controls is analyzed qualitatively in Section 3. The optimal control model is considered in Section 4. Simulations of the model are carried out in Section 5 while Section 6 gives the concluding remarks. ## 2 Model formulation and basic properties of the model The total sexually active population at time *t*, denoted by *N*H(*t*), is divided into eleven mutually-exclusive compartments: Susceptible individuals (*S*(*t*)), individuals infected with HPV (*I*HP (*t*)), individuals who have recovered from or cleared HPV infection (*R*HP (*t*)), individuals with persistent HPV infection (*P*HP (*t*)), individuals with latent TB (*E*T(*t*)), individuals with active TB infection (*A*T(*t*)), individuals treated of TB (*T*T(*t*)), individuals dually infected with HPV and latent ![Graphic][1], individuals dually infected with HPV and active ![Graphic][2], individuals dually infected with persistent HPV and latent TB, ![Graphic][3] individuals dually infected with persistent HPV and active ![Graphic][4]. Thus ![Formula][5] Based on the above formulations and assumptions, the HPV-TB Coinfection model is given by the following deterministic system of non-linear differential equations (flow diagram of the model is shown in Figure 1, the associated state variables and parameters are well described in Tables 1 and 2): View this table: [Table 1:](http://medrxiv.org/content/early/2020/09/27/2020.09.15.20195297/T1) Table 1: Description of variables in the model (1). View this table: [Table 2:](http://medrxiv.org/content/early/2020/09/27/2020.09.15.20195297/T2) Table 2: Description of parameters in the model (1). ![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/27/2020.09.15.20195297/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/09/27/2020.09.15.20195297/F1) Figure 1: Schematic diagram of the model (1) ![Formula][6] where: ![Formula][7] ![Formula][8] The parameters *κ*t(*κ*t ≥ 1) is a modification parameter accounting for the increased infectiousness of dually infected individuals due to tuberculosis, *η*t(*η*t ≤ 1) is a modification parameter accounting for reduced infectiousness of dually infected individuals due to latent TB. *ω**p*(*ω**p* ≥ 1) is a modification parameter accounting for the increased infectiousness of dually infected individuals due to HPV while *η**p*(*η**p* ≤ 1) is a modification parameter accounting for the reduced infectiousness of singly infected and dually infected individuals due to persistent HPV infection. This population is further reduced by natural death (at a rate *µ*H natural death occurs in all epidemiological compartments at this rate). In (2), *β*HP is the effective contact rate for transmission of HPV infection. ### 2.1 Basic properties of the co-infection model (1) without controls The basic qualitative properties of the oncogenic HPV-TB co-infection model (1) will now be considered. Specifically, we establish the following results. #### 2.1.1 Positivity and boundedness of solutions For the oncogenic HPV-TB co-infection model (1) to be meaningful in the concept of epidemiological sense, it is necessary to prove that all its state variables are non-negative over time. The approach adopted by Omame *et al* [23] *will be used to establish the results below*. Theorem 2.1 *Let the initial data be* ![Graphic][9] *Then the solutions* ![Graphic][10] *of the model* (1) *are non-negative for all time t >* 0. Lemma 2.1 *The region* ![Graphic][11] *is positively-invariant for the co-infection model* (1) *with initial conditions in* ![Graphic][12]. ### 3 Mathematical Analysis of the model without controls #### 3.1 Basic reproduction number of the co-infection model The co-infection model (1) has a DFE, obtained by setting the right-hand sides of the equations in the model (1) to zero, given by ![Formula][13] The basic reproduction number of the oncogenic HPV-TB co-infection model (1), using the approach in van den Driessche and Watmough [35], is given by ℛ = max{ℛ0H, ℛ0T} where ℛ0H and ℛ0T are, respectively, the oncogenic HPV and TB associated reproduction numbers, given by ![Formula][14] and ![Formula][15] where, ![Formula][16] Using Theorem 2 of [35], the following result is established. Lemma 3.1 *The DFE, ξ*, *of the oncogenic HPV-TB co-infection model* (1) *is locally asymptotically stable if* ℛ *<* 1, *and unstable if* ℛ *>* 1. ### 3.2 Global asymptotic stability(GAS) of the disease-free equilibrium(DFE) *ξ* of the co-infection model We shall apply the approach illustrated in [5] to investigate the global asymptotic stability of the disease free equilibrium of the co-infection model. In this section, we list two conditions that if met, also guarantee the global asymptotic stability of the disease-free state. First, system (1) must be written in the form: ![Formula][17] where *V* ∈ *R**m* denotes (its components) the number of uninfected individuals and *K* ∈ *R**n* denotes (its components) the number of infected individuals. *U* = (*V* *∗*, 0) denotes the disease-free equilibrium of this system. The conditions (*W* 1) and (*W* 2) following must be satisfied in order to guarantee local asymptotic stability: (*W* 1): For ![Graphic][18] *V* *∗*is globally asymptotically stable (GAS), (*W* 2): ![Graphic][19] for (*V, K*) ∈ Ω, where *B* = *D**K**Q*(*V* *∗*, 0) is an M-matrix (the off-diagonal elements of *B* are nonnegative) and Ω is the region where the model makes biological sense. If System (1) satisfies the above two conditions then the following theorem holds: Theorem 3.1 *The fixed point U* = (*V* *∗*, 0) *is a globally asymptotic stable (GAS) equilibrium of* (1) *provided that R* *<* 1 *(LAS) and that assumptions* (*W* 1) *and* (*W* 2) *are met* **Proof** ![Formula][20] ![Formula][21] where *V* denotes the number of non-infectious compartments and *K* denotes the number of infectious compartments ![Formula][22] It is clear from the above, that, ![Graphic][23]. Hence the DFE may not be globally asymptotically stable, suggesting the possibility of a backward bifurcation. This supports the backward bifurcation analysis in the proceeding section. ### 3.3 Backward bifurcation analysis of the model without controls We shall carry out analysis in this section to know the type of bifurcation the model (1) may undergo, using the Centre Manifold Theory as illustrated in [6]. The following result can be established. Theorem 3.2 *Suppose a backward bifurcation coefficient a >* 0, *(with a defined below), when* ℛ < 1 ![Formula][24] *then model* (1) *exhibits backward bifurcation at* ℛ = 1. *If a <* 0, *then the system* (1) *exhibits a forward bifurcation at* ℛ = 1. **Proof** Suppose ![Formula][25] represents any arbitrary endemic equilibrium of the model. The existence of backward bifurcation will be studied using the Centre Manifold Theory [6]. To apply this theory, it is appropriate to do the following change of variables. Let ![Formula][26] so that ![Formula][27] Further, using the vector notation ![Formula][28] the model (1) can be re-written in the form ![Formula][29] as follows: ![Formula][30] with: ![Formula][31] Consider the case when ℛ0H = 1. Assume, further, that *β*HP is chosen as a bifurcation parameter. Solving for ![Graphic][32] from ℛ0H = 1 gives ![Formula][33] Evaluating the Jacobian of the system (8) at the DFE, *J*(*ξ*), and using the approach in [6], we have that *J*(*ξ*) has a right eigenvector (associated with the simple zero eigenvalue) given by ![Formula][34] where, ![Formula][35] The components of the left eigenvector of ![Graphic][36], satisfying **v.w** = 1 are ![Formula][37] where, ![Formula][38] The associated bifurcation coefficients defined by *a* and *b*, given by: ![Formula][39] are computed to be ![Formula][40] and ![Formula][41] Since the bifurcation coefficient *b* is positive, it follows from Theorem 4.1 in [6] that the model (1), or the transformed model (8), will undergo a backward bifurcation if the backward bifurcation coefficient, *a*, given by (11) is positive. Setting the HPV and TB re-infection paramters *ϕ*p = 0, *σ*T = 0 and the TB exogenous re-infection term ![Graphic][42], the bifurcation coefficient, ![Graphic][43] and ![Graphic][44], based on the definition of the components in (9), and also noting that all other components of the left and right eigenvectors occuring in the coefficient, *a*, are positive, including the parameters of the model. Hence, backward bifurcation does not occur in the Oncogenic HPV-TB co-infection model, in the absence of HPV and TB re-infection and in the absence of exogenous re-infection. □ ## 4 Optimal control model In this section, we shall use the Pontryagin ‘s Maximum Principle to determine the necessary conditions for the optimal control of the oncogenic HPV-TB co-infection model. We incorporate time dependent controls into the model (1) to determine the optimal strategy for curbing the co-infections of the two diseases. Thus, we have, ![Formula][45] subject to the initial conditions ![Graphic][46] with: ![Formula][47] The control functions, *u*1(*t*), *u*2(*t*), and *u*3(*t*) are bounded, Lebesgue integrable functions. The control *u*1(*t*) represents the efforts at preventing incident HPV infections. The control *u*2(*t*) aims at preventing infection with HPV by TB infected individuals so as to reduce the co-infection cases. TB Treatment control for individuals dually infected with oncogenic HPV and TB is denoted by *u*3(*t*). The controls *u*1 and *u*2 satisfies 0 ≤ *u*1, *u*2 ≤ 1 while the control *u*3 satisfies 0 *< u*3 ≤ *h*, where *h* is the TB drug efficacy used for the treatment of co-infected individuals. Our optimal control problem involves a scenario where the number of HPV-infected, TB-infected, the co-infection cases and the cost of implementing preventive and treatment controls *u*1(*t*), *u*2(*t*) and *u*3(*t*) are minimized subject to the state system (12). For this, we consider the objective functional ![Formula][48] *T* is the final time. We seek to find an optimal control, ![Graphic][49], such that ![Formula][50] where ![Graphic][51] such that ![Graphic][52] are measurable with ![Graphic][53], for *t* ∈ [0, *T*] is the control set. We shall now establish the existence of such an optimal solution which minimizes the objective functional *J*. Theorem 4.1 *Given the objective functional J, defined on the control set U, and subject to the state system* (12) *with non-negative initial conditions at t* = 0, *then there exists an optimal control triple u*∗ = (*u*1, *u*2, *u*3) *such that J*(*u*∗) = *min* {*J*(*u*1, *u*2, *u*3)|*u*1, *u*2, *u*3 ∈ *U*}. We shall prove the existence of the given optimal control *u*∗, along with the corresponding solution trajectories by following the approach used in Mohammed-Awel and Numfor [18]. **Proof** The state functions are positive and the controls are Lebesgue measurable, therefore we have that *J*(*u**i*) ≥ 0 for all *u**i* ∈ *U*. As a result, ![Graphic][54] exists and is finite. Therefore, a minimizing sequnce of controls (*u**i*) ∈ *U* exists such that ![Formula][55] Let ![Graphic][56] be the associated state trajectories. Since the state sequences are uniformly bounded, we have that the derivatives are also uniformly bounded. As a result, the state sequences are Lipschitz continuous with the same constant. Applying the Arzela-Ascoli Theorem [13], there exist ![Graphic][57] such that on a sub-sequence ![Formula][58] where, ![Formula][59] Since ‖*u**i*‖*L**∞ < K* for some *K >* 0, it follows that *u**i* ∈ *L*2 ([0, *T*]), such that on a sub-sequence ![Formula][60] Applying the lower semi-continuity of *L*2 norm with respect to weak convergence, we have that ![Formula][61] Considering the convergence of the sequences ![Formula][62] and passing to the limit in the ordinary differential equation system (12), we have that ![Graphic][63] are the states corresponding to the control triple (*u**i*). Hence, (*u**i*) is an optimal control triple. The Pontryagin ‘s Maximum Principle [28] gives the necessary conditions which an optimal control pair must satisfy. This principle transforms (12), (13) and (14) into a problem of minimizing a Hamiltonian, ℋ, pointwisely with regards to the control functions, *u*1, *u*2, *u*3: ![Formula][64] Theorem 4.2 *For an optimal control set u*1, *u*2, *u*3 *that minimizes J over U, there are adjoint variables, λ*1, *λ*2, …, *λ*11 *Satisfying* ![Formula][65] *and with transversality conditions* ![Formula][66] *Furthermore*, ![Formula][67] **Proof of Theorem 4.2** Suppose ![Graphic][68] is an optimal control and ![Graphic][69] are the corresponding state solutions. Applying the Pontryagin ‘s Maximum Principle [28], there exist adjoint variables satisfying: ![Formula][70] with transversality conditions; ![Formula][71] We can determine the behaviour of the control by differentiating the Hamiltonian, ℋ with respect to the controls(*u*1, *u*2, *u*3) at *t*. On the interior of the control set, where 0 *< u**j* *<* 1 for all (*j* = 1, 2, 3), we obtain ![Formula][72] Therefore, we have that [14] ![Formula][73] ![Formula][74] ## 5 Simulations In this section, we shall explore uncertainty and sensitivity analyses on the parameters of the model due to lack of precision in the estimates of some of the co-infection model parameters and the uncertainty which may occur when gathering data for the simulations. We shall also carry out numerical simulations of the optimal control model (12), in order to assess the impact of different intervention strategies on the dynamics of oncogenic HPV-TB co-infections. ### 5.1 Uncertainty and sensitivity analyses The oncogenic HPV-TB co-infection model (1) has thirty nine (39) parameters, and uncertainties are expected to arise in the estimates of their values used in the numerical simulations. Adopting the approach used by Blower and Dowlatabadi [3], we shall carry out a Latin Hypercube Sampling (LHS) on the parameters of the model. For the sensitivity analysis, a Partial Rank Correlation Coefficient (PRCC) was calculated between values of the parameters in the response function and the values of the response function derived from the sensitivity analysis. A total of 1,000 simulations of the co-infection model (1) *per* LHS were run. Using the total number of individuals dually infected with HPV and latent ![Graphic][75] as the response function, the parameters that strongly influence the dynamics of the co-infection model (1) are the demographic parameter, *µ*H, the effective contact rate for HPV transmission, *β*HP, the parameter accounting for increased susceptibility to HPV by TB infected individuals, ![Graphic][76], and the recovery rate from HPV for individuals in ![Graphic][77] compartment. Also, using the population of individuals dually infected with HPV and active ![Graphic][78] as the input, the five top ranked parameters are the effective contact rate for HPV transmissibility, *β*HP, the effective contact rate neccesary for TB transmission, *β*T, the modification parameter accounting for increased susceptibility to HPV infection by active TB infected individuals, ![Graphic][79], the recovery rate from HPV for individuals in ![Graphic][80]compartment and the recovery rate from TB for individuals in ![Graphic][81] class. Taking the population of individuals dually infected with persistent HPV and latent ![Graphic][82] as the input, the three top ranked parameters are the effective contact rate for HPV transmissibility, *β*HP, the fraction of individuals who recover from HPV infection and do not progress to persistent HPV infection, ![Graphic][83] and the recovery rate from persistent HPV infection, ![Graphic][84], for individuals in ![Graphic][85] compartment. When the population of individuals dually infected with persistent HPV and active TB is used as the response function, the two top-ranked parameters are the effctive contact rate for TB transmission, *β*T and the TB treatment rate, ![Graphic][86], for individuals dually infected with persistent HPV and active TB. Considering the HPV associated basic reproduction number, ℛ0H, as the response function, it is observed in Table 3 that the top four PRCC-ranked parameters are the demographic parameter, *µ*H, effective contact rate for HPV transmission, *β*HP, the fraction of infected individuals who recover from HPV and do not develop persistent HPV, ![Graphic][87] and the recovery rate from HPV ![Graphic][88]. Finally, using the TB associated reproduction number, ℛ0T, as the response function the two key parameters that drive the dynamics of the model are the effective contact rate for TB transmission, *β*t and the TB treatment ratein ![Graphic][89]. View this table: [Table 3:](http://medrxiv.org/content/early/2020/09/27/2020.09.15.20195297/T3) Table 3: PRCC values for the co-infection model (1) parameters using the total number of individuals dually infected with: HPV and latent ![Graphic][90] HPV and active ![Graphic][91], Persistent HPV and latent TB![Graphic][92], Persistent HPV and active ![Graphic][93], respectively, as well as the associated reproduction numbers *ℛ*0H and *ℛ*0T, as response functions. Paramters which strongly inuence the dynamics of the co-infectionmodel with respect to each of the response functions are shown in bold fonts. ### 5.2 Numerical simulations Numerical simulations of the optimal control problem (12), adjoint equations (18) and characterizations of the control (21) are implemented by the Runge Kutta method using the forward backward sweep (carried out in MATLAB). Following the choice of the weight function for HPV vaccination control used in Malik *et al* [17], *the balancing factor χ*1 = 103. We also assume *χ*2 = 500 and *χ*3 = 400. The demographic data related to the Shanxi province in rural China is used [8]. We assume the initial conditions to be: ![Graphic][94],.We implement the following three different control strategies for numerical simulations of the co-infection model (4). **The unit of time is years, throughout the entire simulations**. 1. Strategy A: Control against incident HPV infection (*u*1 ≠0) and control against infection with HPV by TB infected individuals(*u*2 ≠ 0); 2. Strategy B: Control against incident HPV infection (*u*1 ≠ 0) and TB treatment control for dually infected individuals (*u*3 ≠ 0); 3. Strategy C: control against infection with HPV by TB infected individuals (*u*2 ≠ 0) and TB treatment control for dually infected individuals (*u*3 ≠ 0). #### 5.2.1 Strategy A: Control against incident HPV infection (*u*1 ≠ 0) and control against Cinfection with HPV by TB infected individuals (*u*2 ≠ 0) Applying this control strategy, we note that the total number of individuals dually infected with HPV and TB in latent and active stages, respectively, (Figures 2 (a) and 2(b)), and the total number of individuals dually infected with persistent HPV and TB in latent and active stages, respectively (Figures 2(c) and 2(d)) are less when this control is applied than when the control is not applied. To be specific, this control strategy averts almost 420,352 new co-infection cases, which is the greatest number of averted cases in comparison with other control strategies implemented. The control profile for this strategy given in Figure 5, reveals that control *u*1 is at its upper bound for the first 3.9 years before ultimately declining to zero. Similarly, the control *u*2 is at the maximum value of 90% for the first 8 months before declining to zero at time, *t*≠4.5 years. #### 5.2.2 Strategy B: Control against incident HPV infection (*u*1 ≠ 0) and TB treatment control for dually infected individuals (*u*3 ≠ 0) The simulations of the total number of dually infected individuals in the presence of TB treatment controls are depicted in Figures 3(a)-3(d). Applying this control, we observe that the total number of individuals dually infected with HPV and TB in latent and active stages of infection respectively, is less than the total population when no control is applied as expected (Figure 3(a) and 3(b)). Similarly, it is noticed from Figures (3(c) and 3(d), that high population level impact is observed in the total number of individuals dually infected with persistent HPV and TB in latent and active stages of infection, respectively, when this control strategy is applied. Specifically, when this control strategy is implemented, about 300,549 new co-infection cases were averted. The control profile for this strategy presented in Figure 7 shows that control *u*1 is at its upper bound for the first 4.0 years before steadily declining to zero at final time. Also, the control *u*3 is at the maximum value of 100% for the first 1.2 years and then steadily declines to zero at time, *t*≠4.5 years. This simulation results conforms with the epidemiological report in the introduction section that gynaecologic TB is a risk factor for oncogenic HPV infection (and subsequent cancer infection) [37]. Hence, if we focus on TB treatment controls, it can significantly reduce the burden of the co-infection of oncogenic HPV and TB in a population. The simulations are also in line with the point opined by [37], that prior TB infection was associated with persistent HPV and increased susceptibility to cervical cancer. As a result, treating TB infections in dually infected individuals will significantly curb the mixed infections of both diseases. ![Figure 2:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/27/2020.09.15.20195297/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/09/27/2020.09.15.20195297/F2) Figure 2: Plots of the total number of individuals dually infected with HPV and TB in latent and active stages, respectively (Figures 2a and 2b) as well as total number of individuals dually infected with persistent HPV and TB in latent and active stages (Figures 2c and 2d), respectively, the presence of optimal vaccination control against incident HPV infection (*u*1 ≠ 0) and control against infection with HPV by TB infected individuals (*u*2 ≠ 0). Here, *β*HP ≠ 2.0, *β*T ≠ 8.557. All other parameters as in Table 2 ![Figure 3:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/27/2020.09.15.20195297/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/09/27/2020.09.15.20195297/F3) Figure 3: Plots of the total number of individuals dually infected with HPV and TB in latent and active stages, respectively (Figures 3(a) and 3(b) as well as total number of individuals dually infected with persistent HPV and TB in latent and active stages (Figures 3(c) and 3(d), respectively, the presence of optimal vaccination control against incident HPV infection (*u*1 = 0) and TB treatment control dually infected individuals (*u*3 = 0). Here, *β*HP = 2.0, *β*T = 8.557. All other parameters as in Table 2 ![Figure 5:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/27/2020.09.15.20195297/F4.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2020/09/27/2020.09.15.20195297/F4) Figure 5: Combined effects of optimal controls *u*1 and *u*2 on the dynamics of the co-infection optimal control model (12) ![Figure 7:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/27/2020.09.15.20195297/F5.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2020/09/27/2020.09.15.20195297/F5) Figure 7: Combined effects of optimal controls *u*1 and *u*3 on the dynamics of the co-infection optimal control model (12) #### 5.2.3 Strategy C: control against infection with HPV by TB infected individuals (*u*2 ≠ 0) and TB treatment control for dually infected individuals (*u*3 ≠ 0) Using this control strategy, we observe in Figures 10(a) and 10 (b), that the number of individuals dually infected with HPV and TB in latent and active stages of infection, respectively, is less than the number when no control strategy is applied. Likewise, this strategy has positive population level impact on the populations of infected individuals dually infected with persistent HPV and TB in latent and active stages of infection, respectively, Figures 10(c) and 10 (d). When this intervention strategy is administered, approximately 362,573 new co-infection cases are averted. The control profile for this strategy, depicted by Figure 9 reveals that control *u*2 is at its upper bound for the first 4.5 years before gradually falling down to zero at final time. Moreso, the control *u*3 is at its minimum value for the first 2.1 years and then steadily rises to its peak value of 100% at time, *t*=2.7 years, before gradually declining to zero at final time. It is imperative to state categorically, that the results of the simulations are based on the parameter values given in table 2 and the initial conditions and weight constants given in Section 5.2 ![Figure 9:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/27/2020.09.15.20195297/F6.medium.gif) [Figure 9:](http://medrxiv.org/content/early/2020/09/27/2020.09.15.20195297/F6) Figure 9: Combined effects of optimal controls *u*2 and *u*3 on the dynamics of the co-infection optimal control model (12) ![Figure 10:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/27/2020.09.15.20195297/F7.medium.gif) [Figure 10:](http://medrxiv.org/content/early/2020/09/27/2020.09.15.20195297/F7) Figure 10: Plots of the total number of individuals dually infected with HPV and TB in latent and active stages, respectively (Figures 10a and 10b) as well as total number of individuals dually infected with persistent HPV and TB in latent and active stages (Figures 10c and 10d), respectively, in the presence of control against HPV infection by TB-infected individuals (*u*2 ≠ 0) and TB treatment control dually infected individuals (*u*3 ≠ 0). Here, *β*HP ≠ 2.0, *β*T ≠ 8.557. All other parameters as in Table 2 ![Figure 12:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/09/27/2020.09.15.20195297/F8.medium.gif) [Figure 12:](http://medrxiv.org/content/early/2020/09/27/2020.09.15.20195297/F8) Figure 12: Cost functions of the different control strategies ### 5.3 Cost-effectiveness analysis The cost-effectiveness analysis is used to evaluate the health interventions related benefits so as to justify the costs of the strategies [4]. This is obtained by comparing the differences among the health outcomes and costs of those interventions; achieved by computing the incremental cost-effectiveness ratio (ICER), which is defined as the cost per health outcome. It is given by: ![Formula][95] We calculated the total number of co-infection cases averted and the total cost of the strategies applied in Table 4. The total number of co-infection cases prevented is obtained by calculating the total number of individuals when controls are administered and the total number when control is not applied. Similarly, we apply the cost functions ![Graphic][96], over time, to compute the total cost for the various strategies that we implemented. We now compare the cost-effectiveness of strategy B (Optimal HPV vaccination control for sexually active susceptible individuals against incident HPV infection and TB treatment control for dually infected individuals) and strategy C (Control against HPV infection by TB infected individuals and TB treatment controls for dually infected individuals). View this table: [Table 4:](http://medrxiv.org/content/early/2020/09/27/2020.09.15.20195297/T4) Table 4: Increasing order of the total infection averted due to the control strategies ![Formula][97] From ICER (B) and ICER(C), we observe a cost saving of 0.01298 observed for strategy C over strategy B. This implies that strategy B strongly dominated strategy C, showing that strategy B is more costly and less effective compared to strategy C. Therefore, strategy B is removed from subsequent ICER computations, shown in Table 5. We shall now compare strategies C and A. Comparing strategy C (Control against HPV infection by TB infected individuals and TB treatment controls for dually infected individuals) and strategy A (Optimal HPV vaccination control for sexually active susceptible individuals against incident HPV infection and Control against HPV infection by TB infected individuals), we observe that ICER (A) is greater than ICER (C), showing that strategy A strongly dominated strategy C and is more expensive and less effective compared to strategy C. Consequently, strategy C (the strategy that combines and implements control against HPV infection by TB infected individuals as well as TB treatment controls for dually infected individuals) has the least ICER and is the most cost-effective of all the control strategies for the control and management of oncogenic HPV and TB co-infection. This clearly agrees with the results obtained from ACER method in Table 4 that strategy C is the most cost-effective strategy. View this table: [Table 5:](http://medrxiv.org/content/early/2020/09/27/2020.09.15.20195297/T5) Table 5: Increasing order of the total infection averted due to the control strategies ![Formula][98] ## 6 Conclusion In this work, we have developed and presented a co-infection model for Oncogenic HPV and TB with cost-effectiveness optimal control analysis. The full co-infection model was shown to undergo the phenomenon of backward bifurcation when the associated reproduction number is less than unity. It was further shown that TB and HPV re-infection parameters (*ϕ*p ≠ 0 and *σ*t ≠ 0) as well as the exogenous re-infection term (*ε*1 ≠ 0) induced the phenomenon of backward bifurcation 1he in the oncogenic HPV-TB co-infection model. The global asymptotic stability of the disease-free equilibrium of the co-infection model was also proven **not to exist**, when the associated reproduction number was below unity. The necessary conditions for the existence of optimal control and the optimality system for the co-infection model was established using the Pontryagin ‘s Maximum Principle. Uncertainty and global sensitivity analyses were also carried out to determine the top ranked parameters that influence the dynamics of the co-infection model, when the associated reproduction numbers as well as the infected populations were used as response functions. When the population of individuals dually infected with HPV and latent ![Graphic][99] was used as the response function, the parameters that strongly influence the dynamics of the co-infection model (1) are the demographic parameter, *µ*H, the effective contact rate for HPV transmission, *β*HP, the parameter accounting for increased susceptibility to HPV by TB infected individuals, ![Graphic][100], and the recovery rate from HPV hefor individuals in ![Graphic][101] compartment. In addition, using the population of individuals dually infected with HPV and active ![Graphic][102] as the input, the five top ranked parameters are the effective contact rate for HPV transmissibility, *β*HP, the effective contact rate neccesary for TB transmission, *β*T, the modification parameter accounting for increased susceptibility to HPV infection by active TB infected individuals, ![Graphic][103], the recovery rate from HPV for individuals in ![Graphic][104] compartment and the recovery rate from TB for individuals in ![Graphic][105] class. Taking the population of individuals dually infected with persistent HPV and active TB as the response function, the two top-ranked parameters are the effctive contact rate for TB transmission, *β**T*and the TB treatment rate, ![Graphic][106], for individuals dually infected with persistent HPV and active TB. Numerical simulations of the optimal control model showed that: 1. A combination of HPV prevention control and TB treatment control has a positive population level impact in reducing the burden of oncogenic HPV and TB co-infection cases in a population. 2. A strategy that implements control against HPV infection by TB infected individuals and TB treatment controls for dually infected individuals can significantly reduce the burden of oncogenic HPV and TB co-infections. 3. The strategy that combines and implements control against HPV infection by TB infected individuals as well as TB treatment control for dually infected individuals has the least ICER and is the most cost-effective of all the control strategies for the control and management of the burden of oncogenic HPV and TB co-infection. ## Data Availability Not applicable [https://www.indexmundi.com/china/demographics\_profile](https://www.indexmundi.com/china/demographics_profile) ## Footnotes * Added some references * Received September 15, 2020. * Revision received September 23, 2020. * Accepted September 27, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1]. FB Agusto, AI Adekunle, Optimal control of a two-strain tuberculosis-HIV/AIDS co-infection model, BioSystems, (2014) (119) 20–24. 2. [2]. HG Ahmed, SH Bensumaidea, IM Ashankyty, Frequency of human papillomavirus (HPV) subtypes 31, 33, 35, 39 and 45 among Yemeni women with cervical cancer, Infect. Ag. Cancer 10 (29) (2015) 1–6. 3. [3]. SM Blower, H Dowlatabadi, Sensitivity and uncertainty analysis of complex models of disease transmission: an HIV model, as an example, International Statistical Review 2 (1994) 229–243. 4. [4]. SB Cantor, TG Ganiats, Incremental cost-effectiveness analysis: the optimal strategy depends on the strategy set. Clin. Epidemiol. 52 (1999) (6), 517–522. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0895-4356(99)00021-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10408990&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F27%2F2020.09.15.20195297.atom) 5. [5].C Castillo-Chavez, Z Feng, W Huang, On the computation of R0 and its role on global stability, in Mathematical Approaches for Emerging and Reemerging Infectious Diseases: An Introduction (Minneapolis, MN, 1999), 229–250, IMA Vol. Math. Appl., 125 Springer New York. 6. [6]. C Castillo-Chavez, B Song, Dynamical models of tuberculosis and their applications, Math. Biosci. Eng. 2 (2004) 361–404. 7. [7]. T Cohen, C Colijin, B Finklea, M Murray, Exogenous re-infection and the dynamics of tuberculosis epidemics: local effects in a network model of transmission, J. R. Soc. Interface, 4 (2007) 523–531. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rsif.2006.0193&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17251134&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F27%2F2020.09.15.20195297.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000247311700008&link_type=ISI) 8. [8].China Demographics Profile, 2018, Available at [http://www.indexmundi.com/china/demographics](http://www.indexmundi.com/china/demographics) profile. Accessed: 2 May, 2018. 9. [9].Fatmawati, H Tasman, Optimal control of HIV resistance and tuberculosis co-infection using treatment intervention, Asian Pac. J. Trop. Dis. 7 (2017) 366–373. 10. [10]. FM Gilbert, I Ihejirika, SC Inyama, A Omame, HI Mbachu, JI Uwakwe, Deterministic Mathematical Model of Tuberculosis Disease with Treatment and Recovered Groups, International Journal of Mathematics and Statistics Invention (2019) 7 (1), 47–60 11. [11]. AO Isere, D Okuonghae, Optimal control model for the outbreak of cholera in Nigeria, African Journal of Mathematics and Computer Science Research, 7;24–30; 2014 12. [12]. A Jemal, F Bray, MM Center, J Ferlay, E Ward, D Forman, Global Cancer Statistics, CA Cancer J. Clin. 61 (2011) 69–90. 13. [13]. J. L. Kelley, General topology, Springer-Verlag, (1991) ISBN 978-0-387-90125-1 14. [14]. S Lenhart, JT Workman, 2007. Optimal Control Applied to Biological Models. Chapman & Hall, Boca Raton. 15. [15]. A MacNeil, P Glaziou, C Sismanidis, S Maloney, K Floyd, Global Epidemiology of Tuberculosis and Progress Toward Achieving Global Targets - 2017. MMWR Morb Mortal Wkly Rep (2019) 68 263–266 DOI: [http://dx.doi.org/10.15585/mmwr.mm6811a3](http://dx.doi.org/10.15585/mmwr.mm6811a3) 16. [16]. MT Malik, J Reimer, AB Gumel, EH Elbasha, SM Mahmud, The impact of an imperfect vaccine and pap cytology screening on the transmission of Human Papillomavirus and occurrence of associated cervical dysplasia and cancer, Math. Biosci. Eng. 10 (4) (2013). 17. [17]. T Malik, M Imran, R Jayaraman, Optimal control with multiple human papillomavirus vaccines, Journal of Theoretical Biology, (393) 179–193, 2016. 18. [18]. J Mohammed-Awel, E Numfor, Optimal Insecticide Treated Bednets Coverage and Malaria Treatment in a Malaria-HIV Co-infection Model. Journal of Biological Dynamics, 11 (2017) -Issue sup1 : 160–191. 19. [19]. Z Mukandavire, AB Gumel, W Garira, JM Tchuenche, Mathematical Analysis of a Model for HIV-Malaria Co-infection, Math. Biosci. Eng. 6 (2) (2009) 333–362. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3934/mbe.2009.6.333&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19364156&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F27%2F2020.09.15.20195297.atom) 20. [20]. KO Okosun, OD Makinde, A co-infection model of malaria and cholera diseases with optimal control, Math. Biosci. 258 (2014) 19–32. 21. [21]. D Okuonghae, A Korobeinikov, Dynamics of tuberculosis: the effect of direct observation therapy strategy (DOTS) in Nigeria, Math. Model. Nat. Phenomena, (2) (2007) 1 99–111. 22. [22]. D Okuonghae, V Aihie, Case detection and direct observation therapy strategy (DOTS) in Nigeria: its effect on TB dynamics, J. Biol. Syst. 16 (1) (2008) 1–31. 23. [23]. A Omame, RA Umana, D Okuonghae, SC Inyama, Mathematical analysis of a two-sex Human Papillomavirus (HPV) model, International Journal of Biomathematics, 11 (2018) (7) 24. [24]. A Omame, D Okuonghae, RA Umana, SC Inyama, Analysis of a co-infection model for HPV-TB, Applied Mathematical Modelling, 77 (2020) 881–901. 25. [25].1. F. Smith, 2. H. Dutta and 3. J.N. Mordeson (eds A Omame, D Okuonghae, SC Inyama, A mathematical study of a model for HPV with two high risk strains, in Mathematics Applied to Engineering, Modelling, and Social Issues Studies in Systems, Decision and Control (2020) 200, F. Smith, H. Dutta and J.N. Mordeson (eds.) 26. [26]. A Omame, N Sene, I Nometa, CI Nwakanma, EU Nwafor, NO Iheonu, D Okuonghae, Analysis of COVID-19 and comorbidity co-infection Model with Optimal Control, medRxiv preprint (2020) doi: [https://doi.org/10.1101/2020.08.04.20168013](https://doi.org/10.1101/2020.08.04.20168013). 27. [27]. A Omame, CU Nnanna, SC Inyama, Optimal Control and Cost-effectiveness Analysis of an HPV-Chlamydia Trachomatis co-infection model, medRxiv preprint (2020) doi:[https://doi.org/10.1101/2020.09.07.20190025](https://doi.org/10.1101/2020.09.07.20190025). 28. [28]. LS Pontryagin, VG Boltyanskii, RV Gamkrelidze, EF Mishchenko, The Mathematical Theory of Optimal Processes, Wiley, New York, 1962. 29. [29]. F Saldana, A Korobeinikov, I Barradas, Optimal Control against the Human Papillomavirus: Protection versus Eradication of the Infection, Abstract and Applied Analysis, [https://doi.org/10.1155/2019/4567825](https://doi.org/10.1155/2019/4567825) 30. [30].Shanxi Demographics Data, 2015, Available at [http://www.khoema.com/atlas/china/shanxi](http://www.khoema.com/atlas/china/shanxi) Accessed: 2 May, 2018. 31. [31].Tanvi, R Aggarwal, Dynamics of HIV-TB co-infection with detection as optimal intervention strategy, International Journal of Non-Linear Mechanics (2020), doi:[https://doi.org/10.1016/j.ijnonlinmec.2019.103388](https://doi.org/10.1016/j.ijnonlinmec.2019.103388). 32. [32]. GT Tilahun, OD Makinde, D Malonza, Co-dynamics of Pneumonia and Typhoid fever diseases with cost effective optimal control analysis, Appl. Math. Comput. 316 (2018) 438–459. 33. [33]. RA Umana, A Omame, SC Inyama, Deterministic and Stochastic Models of the Dynamics of Drug Resistant Tuberculosis; FUTO Journals Series, (2016) 2 (2) 173–194. 34. [34]. JI Uwakwe, SC Inyama, A Omame, Mathematical Model and Optimal Control of New-Castle Disease (ND). Applied and Computational Mathematics. (2020) 9 (3) 70–84. doi: 10.11648/j.acm.20200903.14 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.11648/j.acm.20200903.14&link_type=DOI) 35. [35]. P van den Driessche, J Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci. 180 (2002) 29–48. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0025-5564(02)00108-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12387915&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F09%2F27%2F2020.09.15.20195297.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000179220600004&link_type=ISI) 36. [36].World Health Organization. Global tuberculosis report 2018. Geneva, Switzerland: World Health Organization; 2018. [https://www.who.int/tb/publications/globalreport/en](https://www.who.int/tb/publications/globalreport/en) 37. [37]. NM Zetola, S Grover, C Modongo, SP Chiyapo, Collision of three pandemics: The coexistence of cervical cancer, HIV infection, and the prior tuberculosis in the sub-Saharan country of Botswana, J. Glob. Oncol. DOI:10.1200/JGO.2015.001701. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamdvIjtzOjU6InJlc2lkIjtzOjY6IjIvMS80NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA5LzI3LzIwMjAuMDkuMTUuMjAxOTUyOTcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 38. [38]. S Zhu, L Xia, S Yu, S Chen, J Zhang, The burden and challenges of tuberculosis in China: findings from the Global Burden of Disease Study 2015, Nature, 7 (2017) 14601, DOI:610.1038/s41598-017-15024-1. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-017-15024-1&link_type=DOI) 39. [39]. FH Zhao, AP Varanasi, CA Cunningham, BI Graubard, S-Y Hu, F Chen, CJ Barrett, Y-L Qiao, MR Forman, Tuberculosis and oncogenic HPV: Potential co-infections in women at high-risk of cervical cancer in rural China. Asian P. J. Canc. Prev. 12 (2011) 1409–1415. [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/inline-graphic-4.gif [5]: /embed/graphic-1.gif [6]: /embed/graphic-5.gif [7]: /embed/graphic-6.gif [8]: /embed/graphic-7.gif [9]: /embed/inline-graphic-5.gif [10]: /embed/inline-graphic-6.gif [11]: /embed/inline-graphic-7.gif [12]: /embed/inline-graphic-8.gif [13]: /embed/graphic-8.gif [14]: /embed/graphic-9.gif [15]: /embed/graphic-10.gif [16]: /embed/graphic-11.gif [17]: /embed/graphic-12.gif [18]: /embed/inline-graphic-9.gif [19]: /embed/inline-graphic-10.gif [20]: /embed/graphic-13.gif [21]: /embed/graphic-14.gif [22]: /embed/graphic-15.gif [23]: /embed/inline-graphic-11.gif [24]: /embed/graphic-16.gif [25]: /embed/graphic-17.gif [26]: /embed/graphic-18.gif [27]: /embed/graphic-19.gif [28]: /embed/graphic-20.gif [29]: /embed/graphic-21.gif [30]: /embed/graphic-22.gif [31]: /embed/graphic-23.gif [32]: /embed/inline-graphic-12.gif [33]: /embed/graphic-24.gif [34]: /embed/graphic-25.gif [35]: /embed/graphic-26.gif [36]: /embed/inline-graphic-13.gif [37]: /embed/graphic-27.gif [38]: /embed/graphic-28.gif [39]: /embed/graphic-29.gif [40]: /embed/graphic-30.gif [41]: /embed/graphic-31.gif [42]: /embed/inline-graphic-14.gif [43]: /embed/inline-graphic-15.gif [44]: /embed/inline-graphic-16.gif [45]: /embed/graphic-32.gif [46]: /embed/inline-graphic-17.gif [47]: /embed/graphic-33.gif [48]: /embed/graphic-34.gif [49]: /embed/inline-graphic-18.gif [50]: /embed/graphic-35.gif [51]: /embed/inline-graphic-19.gif [52]: /embed/inline-graphic-20.gif [53]: /embed/inline-graphic-21.gif [54]: /embed/inline-graphic-22.gif [55]: /embed/graphic-36.gif [56]: /embed/inline-graphic-23.gif [57]: /embed/inline-graphic-24.gif [58]: /embed/graphic-37.gif [59]: /embed/graphic-38.gif [60]: /embed/graphic-39.gif [61]: /embed/graphic-40.gif [62]: /embed/graphic-41.gif [63]: /embed/inline-graphic-25.gif [64]: /embed/graphic-42.gif [65]: /embed/graphic-43.gif [66]: /embed/graphic-44.gif [67]: /embed/graphic-45.gif [68]: /embed/inline-graphic-26.gif [69]: /embed/inline-graphic-27.gif [70]: /embed/graphic-46.gif [71]: /embed/graphic-47.gif [72]: /embed/graphic-48.gif [73]: /embed/graphic-49.gif [74]: /embed/graphic-50.gif [75]: /embed/inline-graphic-28.gif [76]: /embed/inline-graphic-29.gif [77]: /embed/inline-graphic-30.gif [78]: /embed/inline-graphic-31.gif [79]: /embed/inline-graphic-32.gif [80]: /embed/inline-graphic-33.gif [81]: /embed/inline-graphic-34.gif [82]: /embed/inline-graphic-35.gif [83]: /embed/inline-graphic-36.gif [84]: /embed/inline-graphic-37.gif [85]: /embed/inline-graphic-38.gif [86]: /embed/inline-graphic-39.gif [87]: /embed/inline-graphic-40.gif [88]: /embed/inline-graphic-41.gif [89]: /embed/inline-graphic-42.gif [90]: T3/embed/inline-graphic-43.gif [91]: T3/embed/inline-graphic-44.gif [92]: T3/embed/inline-graphic-45.gif [93]: T3/embed/inline-graphic-46.gif [94]: /embed/inline-graphic-47.gif [95]: /embed/graphic-59.gif [96]: /embed/inline-graphic-48.gif [97]: /embed/graphic-61.gif [98]: /embed/graphic-63.gif [99]: /embed/inline-graphic-49.gif [100]: /embed/inline-graphic-50.gif [101]: /embed/inline-graphic-51.gif [102]: /embed/inline-graphic-52.gif [103]: /embed/inline-graphic-53.gif [104]: /embed/inline-graphic-54.gif [105]: /embed/inline-graphic-55.gif [106]: /embed/inline-graphic-56.gif