The past, present and future impact of HIV prevention and control on HPV and cervical disease in Tanzania: A modelling study

Abstract

Background

Women with HIV have an elevated risk of HPV infection, and eventually, cervical cancer. Tanzania has a high burden of both HIV and cervical cancer, with an HIV prevalence of 5.5% in women in 2018, and a cervical cancer incidence rate among the highest globally, at 59.1 per 100,000 per year, and an estimated 9,772 cervical cancers diagnosed in 2018. We aimed to quantify the impact that interventions intended to control HIV have had and will have on cervical cancer in Tanzania over a period from 1995 to 2070.

Methods

A deterministic transmission-dynamic compartment model of HIV and HPV infection and natural history was used to simulate the impact of voluntary medical male circumcision (VMMC), anti-retroviral therapy (ART), and targeted pre-exposure prophylaxis (PrEP) on cervical cancer incidence and mortality from 1995–2070.

Findings

We estimate that VMMC has prevented 2,843 cervical cancer cases and 1,039 cervical cancer deaths from 1995–2020; by 2070 we predict that VMMC will have lowered cervical cancer incidence and mortality rates by 28% (55.11 cases per 100,000 women in 2070 without VMMC, compared to 39.93 with VMMC only) and 26% (37.31 deaths per 100,000 women in 2070 without VMMC compared to 27.72 with VMMC), respectively. We predict that ART will temporarily increase cervical cancer diagnoses and deaths, due to the removal of HIV death as a competing risk, but will ultimately further lower cervical cancer incidence and mortality rates by 7% (to 37.31 cases per 100,000 women in 2070) and 5% (to 26.44 deaths per 100,000 women in 2070), respectively, relative to a scenario with VMMC but no ART. A combination of ART and targeted PrEP use is anticipated to lower cervical cancer incidence and mortality rates to 35.82 and 25.35 cases and deaths, respectively, per 100,000 women in 2070.

Conclusions

HIV treatment and control measures in Tanzania will result in long-term reductions in cervical cancer incidence and mortality. Although, in the near term, the life-extending capability of ART will result in a temporary increase in cervical cancer rates, continued efforts towards HIV prevention will reduce cervical cancer incidence and mortality over the longer term. These findings are critical background to understanding the longer-term impact of achieving cervical cancer elimination targets in Tanzania.

About the Authors:

Michaela T. Hall

Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliations School of Mathematics and Statistics, UNSW Sydney, Kensington, Australia, Cancer Research Division, Cancer Council NSW, Woolloomooloo, Australia

ORCID logo http://orcid.org/0000-0002-8989-1449

Megan A. Smith

Roles Investigation, Methodology, Writing – review & editing

Affiliations Cancer Research Division, Cancer Council NSW, Woolloomooloo, Australia, School of Public Health, University of Sydney, Sydney, Australia

Kate T. Simms

Roles Investigation, Methodology, Writing – review & editing

Affiliations Cancer Research Division, Cancer Council NSW, Woolloomooloo, Australia, School of Public Health, University of Sydney, Sydney, Australia

Ruanne V. Barnabas

Roles Data curation, Methodology, Writing – review & editing

Affiliation: University of Washington, Seattle, WA, United States of America

Karen Canfell

Roles Conceptualization, Supervision, Writing – review & editing

Affiliations Cancer Research Division, Cancer Council NSW, Woolloomooloo, Australia, School of Public Health, University of Sydney, Sydney, Australia

ORCID logo http://orcid.org/0000-0002-6443-6618

John M. Murray

Roles Conceptualization, Funding acquisition, Investigation, Methodology, Supervision, Writing – review & editing

Affiliation: School of Mathematics and Statistics, UNSW Sydney, Kensington, Australia

Abstract

Background

Women with HIV have an elevated risk of HPV infection, and eventually, cervical cancer. Tanzania has a high burden of both HIV and cervical cancer, with an HIV prevalence of 5.5% in women in 2018, and a cervical cancer incidence rate among the highest globally, at 59.1 per 100,000 per year, and an estimated 9,772 cervical cancers diagnosed in 2018. We aimed to quantify the impact that interventions intended to control HIV have had and will have on cervical cancer in Tanzania over a period from 1995 to 2070.

Methods

A deterministic transmission-dynamic compartment model of HIV and HPV infection and natural history was used to simulate the impact of voluntary medical male circumcision (VMMC), anti-retroviral therapy (ART), and targeted pre-exposure prophylaxis (PrEP) on cervical cancer incidence and mortality from 1995–2070.

Findings

We estimate that VMMC has prevented 2,843 cervical cancer cases and 1,039 cervical cancer deaths from 1995–2020; by 2070 we predict that VMMC will have lowered cervical cancer incidence and mortality rates by 28% (55.11 cases per 100,000 women in 2070 without VMMC, compared to 39.93 with VMMC only) and 26% (37.31 deaths per 100,000 women in 2070 without VMMC compared to 27.72 with VMMC), respectively. We predict that ART will temporarily increase cervical cancer diagnoses and deaths, due to the removal of HIV death as a competing risk, but will ultimately further lower cervical cancer incidence and mortality rates by 7% (to 37.31 cases per 100,000 women in 2070) and 5% (to 26.44 deaths per 100,000 women in 2070), respectively, relative to a scenario with VMMC but no ART. A combination of ART and targeted PrEP use is anticipated to lower cervical cancer incidence and mortality rates to 35.82 and 25.35 cases and deaths, respectively, per 100,000 women in 2070.

Conclusions

HIV treatment and control measures in Tanzania will result in long-term reductions in cervical cancer incidence and mortality. Although, in the near term, the life-extending capability of ART will result in a temporary increase in cervical cancer rates, continued efforts towards HIV prevention will reduce cervical cancer incidence and mortality over the longer term. These findings are critical background to understanding the longer-term impact of achieving cervical cancer elimination targets in Tanzania.

Figures

Fig 11

Fig 12

Table 1

Fig 1

Table 2

Fig 2

Table 3

Table 4

Table 5

Fig 3

Fig 4

Fig 5

Fig 6

Fig 7

Fig 8

Table 6

Fig 9

Table 7

Fig 10

Table 8

Fig 11

Fig 12

Table 1

Fig 1

Table 2

Citation: Hall MT, Smith MA, Simms KT, Barnabas RV, Canfell K, Murray JM (2020) The past, present and future impact of HIV prevention and control on HPV and cervical disease in Tanzania: A modelling study. PLoS ONE 15(5): e0231388. https://doi.org/10.1371/journal.pone.0231388

Editor: Alberto d’Onofrio, International Prevention Research Institute, FRANCE

Received: September 4, 2019; Accepted: March 3, 2020; Published: May 6, 2020

Copyright: © 2020 Hall et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the manuscript and its Supporting Information files.

Funding: This work was funded by an Australian Government Research Training Program (RTP) Scholarship (Ms Michaela Hall 5045590). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: KC is a co-PI of an investigator-initiated trial of primary HPV screening in Australia (‘Compass’) which is conducted and funded by the VCS Foundation, a government-funded health promotion charity, which has received a funding contribution from Roche Molecular Systems and Roche Tissue Diagnostics, AZ, USA. Neither KC nor her institution on her behalf have received funding from industry for this or any other project. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

Background

For many years human immunodeficiency virus (HIV) has been one of the most heavily researched infectious diseases, and now, controlling HIV is beginning to look achievable [1]. Improved methods of HIV prevention and control such as pre-exposure prophylaxis (PrEP), anti-retroviral therapy (ART) and even voluntary medical male circumcision (VMMC) are at the forefront of health policy recommendations [1–4]. If these interventions are effectively implemented at the population level, they may substantially reduce HIV transmission and eventually end the HIV epidemic. Many modelling studies have attempted to quantify the impact of these interventions on HIV prevalence and related mortality in a range of settings [5–10].

HIV positivity has been linked to higher rates of human papillomavirus (HPV) acquisition, and, among those infected with HPV, the presence of an HIV co-infection is known to reduce the likelihood of HPV clearance and regression of pre-cancerous lesions, and, increase the risk of progression [11]. For this reason, modelling studies evaluating cervical cancer prevention policies are increasingly considering, either directly or indirectly, the impacts of endemic HIV in relevent settings [7, 12, 13].

Methods of HIV control may have a substantial impact not only on prevalence and deaths due to HIV, but also on HPV prevalence and subsequently cervical cancer incidence and mortality rates [11, 14]. In particular, male circumcision has been shown to reduce the risk of HIV-1 acquisition in heterosexual men over a time-period of 18–24 months by at least 60%, and, reduce HPV prevalence among heterosexual men by 63% [15–19]. Reductions in male HIV and HPV prevalence then results in women also experiencing less HIV and oncogenic HPV infection, and subsequently, less cervical cancer [16, 20]. A global ecological analysis classifying VMMC into high (>80%), intermediate (20–80%) and low (<20%) prevalence has reported that for each categorical shift in VMMC prevalence, cervical cancer incidence was reduced by 3.65 (0.54–6.76) cases per 100,000 women per year [21].

The United Republic of Tanzania has a high burden of both HIV and cervical cancer. It was estimated that in 2018, 5.5% of Tanzanian women aged 15–49 years were living with HIV [22], while the incidence of cervical cancer was among the highest globally, at 59.1 cases diagnosed per 100,000 women (9,771 cervical cancers detected) in 2018 [23]. The 2018 incidence rates of cervical cancer in Southern Africa and Eastern Africa were 43.1 cases per 100,000 women per year, and 40.1 cases per 100,000 women per year, respectively [23]. Tanzania is within the sub-Saharan African region, which in 2018 contained 53% of all people living with HIV globally and had an estimated HIV prevalence among adults aged 15–49 years of 7% [24, 25]. The two main interventions against HIV currently in place in Tanzania are ART and VMMC, which are both being actively scaled up [26]; while the Tanzanian Ministry of Health recommends PrEP use for those at significant risk of HIV acquisition, scale up of access to PrEP has been limited [27]. In light of Tanzania’s high burden of cervical cancer and the known impact of HIV, it is important to assess the impact of HIV control interventions that are currently being scaled up (ART and VMMC) or considered (PrEP) on not only HIV incidence and prevalence, but also rates of cervical cancer incidence and mortality. This exercise will provide important context to understanding the impact of scaling up prevention and treatment strategies to achieve cervical cancer elimination targets. Furthermore, while there exists significant variation in national laws pertaining to sexual identity and orientation, sex-work, and access to contraception across the African continent which affect local rates of sexually transmitted diseases [24], the relative impact of HIV interventions on cervical cancer incidence rates in Tanzania is likely to be broadly representative of the region.

The aim of this analysis, therefore, was to quantify the effect of HIV control actions to date on cervical cancer incidence and mortality in terms of rates, cancer diagnoses and lives saved, and to predict future cervical cancer incidence and mortality rates in Tanzania, in the context of scaled-up HIV control interventions.

Methods

Model overview and parameterisation

A detailed deterministic transmission-dynamic compartment model was developed to concurrently simulate the transmission and natural history of HIV, HPV 16/18, HPV 31/33/45/52/58 (referred to as HPV H5) and other oncogenic high-risk HPV types (referred to as HPV OHR) in the United Republic of Tanzania. While there are a range of transmission modalities for HIV and HPV, this platform simulates heterosexual transmission only (a simplifying assumption), as this is the dominant mode of transmission for both HIV and HPV in sub-Saharan Africa [28, 29]. The platform can simulate dual HIV and HPV infections, as well as infections with multiple HPV types, with and without ART. The model incorporates comprehensive demographic, sexual behaviour and natural history assumptions, and accounts for VMMC, ART and PrEP. The simulated population includes males, females and a separate subgroup of female commercial sex-workers (which females may be hired into or retire from), from ages 5 to 79 years, stratified by sex, five-year age group, sexual activity level, HIV and HPV infection, and treatment status. This model is comprised of 11,022,480 compartments, where simulated populations move between the states described in Table 1. Note that all persons in the simulated population are categorised by some combination of attributes: sex/career, age, sexual activity level, HIV infection status and natural history, HIV treatment status (note that no HIV negative individuals are treated with the exception of PReP which may be provided prophylactically for HIV prevention in HIV negative individuals), HPV 16/18 infection status and natural history, HPV H5 infection status and natural history, HPV OHR infection status and natural history and cervical cancer detection status and treatment. Note that only women can progress from HPV infections to cervical pre-cancer, and only women with cervical cancer may have cancer detected. The model implementation utilised for this analysis runs on a quarterly timestep (13 weeks).

[Figure omitted. See PDF.]

Table 1. Model compartments exist for the cartesian product of sets (A) to (I).

https://doi.org/10.1371/journal.pone.0231388.t001

The model’s input parameters were specified primarily using empirical data; however, some parameters, particularly those unobservable or informed by survey data, were found through calibration using a trust region reflective algorithm. The parameter inputs found through calibration were the per-timestep volumes of sex-, age- and activity-group specific high-risk sexual contacts, the degree of age-assortative sexual mixing, annual fluctuations in population-level risk aversive behaviour, and the relative per-sex-act probability of HIV acquisition for females compared to males. These inputs were calibrated to estimated sex-specific HIV prevalence over time, as well as annual rates of new HIV infections obtained from UNAIDS [22, 30, 31]. Note that different groupings in the presented age-range of calibration/validation results are due to variation in reported age-ranges in the observed data.

Demography.

Population demography encompasses compartments for sex/career (male, female, commercial sex-worker), age (five-year age-groups from 5–9 to 75–79 years) and sexual activity level (age- and sex-specific rates of propensity towards sexual risk-taking). The demography module accounts for population ageing, recruitment, natural mortality and assigns risk groups. The youngest simulated age group is 5–9 years, therefore recruitment represents the number of children born who survive to age five, and accounts for the age- and year-specific fertility rates of the simulated female population, as well as infant mortality. The per-timestep probability of any individual ageing to the next five-year age-group is calculated using the number of single-year ages in the age group, and the number of model iterations per year. For example, of individuals in the 10–14 year age-group turn 15 in any given year, and since there are four timesteps simulated per year, the probability of ageing from the 10–14 year group to the 15–19 year group is . In calculating recruitment, annual fertility rates were sourced from the World Bank using the median fertility variant [32], whereas data on maternal age at birth was sourced from the United Republic of Tanzania Ministry of Finance and is based on the 2012 census [33]. The simulated population is subject to an age-specific probability of death resulting from any cause other than HIV or cervical cancer (other cause mortality). Age-and-year-specific mortality rates were derived using the projected year-on-year life tables reported by the United Nations Population Division, adjusted for HIV and cervical cancer mortality [34]. Finally, the demography module re-distributes the simulated population into two sex- and age-specific sexual activity groups (high-activity and general-activity) and simulates the recruitment of women into a career of commercial sex work, and, their eventual retirement. The initial age-distribution was based on the 1960 Tanzanian population [35], with a sex ratio of 1 male to 1.03 females, based on data from the World Bank [36].

Force of infection.

The model simulates HIV and HPV transmission between sexual partners, including interactions between commercial sex workers (CSWs) and their male clients. The population is compartmentalised into ‘high activity’ and ‘general activity’ sexual activity groups, which differ in their assumed number of sexual contacts per timestep. The number of sexual interactions per time-step implicitly accounts for new partners, the per-partnership frequency of sex, and relationship type (casual or monogamous). An age-dependent proportion of the female population were assumed to be CSWs, with an age-specific probability of seeking commercial sex defined for males. Furthermore, CSWs engage in both personal and commercial sexual interactions, with a pre-defined age-specific client volume per timestep.

The model platform calculates the sex- and age-specific per-timestep force of infection using age-specific partnership preferences, sexual activity group, HIV/HPV prevalence among sex partners, the per-sex-act probability of pathogen transmission (stratified by disease stage where applicable) and uptake of preventative interventions such as condom use (specified separately for commercial sex and general partnerships), VMMC prevalence and ART use. For example, the force of infection for HIV for a male aged a in sexual activity group r at time t is calculated using Eq 1.(1)cM(a, r, t) denotes the average number of sexual contacts for a male aged a in sexual-activity group r at time t; kHIV(t) denotes the per-sex-act probability that a condom is worn and prevents HIV acquisition; υHIV(t) denotes the probability that the male has undergone VMMC and the per-sex act-probability that this prevents HIV acquisition; i is the stage of HIV disease among female sex-partners; ρM(a) is a vector of the distribution of preferences for female partners of each age-group for males aged a (vector over all ages summing to unity); is the HIV-stage-specific per sex-act female-to-male transmission probability for females with treatment status Tx; refers to the dot product between vectors ρM(a) and ; is a vector specifying the age-specific probability of a female being HIV positive (simulated), stage i, and with treatment status Tx and time t; and finally, is the probability of acquiring an HIV infection from a commercial sex worker. Note that(2)Where ςM(a, r, t) denotes the average number of commercial sexual contacts assumed for a male aged a in sexual risk group r at time t; and, is a vector specifying the age-specific probability of a commercial sex worker being HIV positive, stage i and with treatment status Tx at time t.

The HPV 16/18 force of infection for a high activity male aged a at time t is calculated in a similar way but simplified by the assumption that the probability of HPV transmission is fixed irrespective of HPV disease stage. That is,(3)kHPV(t) denotes the per-sex-act probability that a condom is worn and prevents HPV acquisition; υHPV(t) denotes the probability that the male has undergone VMMC and the per-sex act-probability that this prevents HPV acquisition; is the per sex-act female-to-male HPV16/18 transmission probability; and, is a vector specifying the age-specific probability of a female being HPV 16/18 positive. Additionally,(4)

Equations specifying the force of infection for HPV H5 and HPV OHR in males and females are similar to the above and are provided in S1 File. Detailed input parameter assumptions relevant to the calculation of force of infection are also described in S1 File (see equations s1-s12).

Disease natural history.

Disease progression for HIV infection is governed by the following state diagram (Fig 1), where specific progression rates are dependent on the current stage of disease and treatment status.

[Figure omitted. See PDF.]

Fig 1. State-space diagram for HIV disease progression.

Note that viral suppression was assumed to halt disease progression and that all states are subject to other cause mortality.

https://doi.org/10.1371/journal.pone.0231388.g001

The natural history of HIV infection progresses from acute HIV infection though four clinical disease stages. These stages are aligned with the World Health Organisation (WHO) Clinical Staging of HIV/AIDS for Adults and Adolescents [37], and are defined in terms of patient symptoms. Input parameters specifying HIV progression rates in the model are described in Table 2.

[Figure omitted. See PDF.]

Table 2. Average length of time spent in each disease stage, and the probability of HIV-death for each HIV disease stage.

https://doi.org/10.1371/journal.pone.0231388.t002

The model platform accounts for the detailed and well-understood natural history of HPV. Disease progression and regression for HPV infection are governed as per Fig 2, where specific progression rates are dependent on HPV type, age, disease stage, HIV positivity and ART treatment status. HPV types 16/18 are known to be more aggressive than other oncogenic HPV types, with elevated disease progression rates and reduced regression rates. Women with an HIV co-infection also experience more aggressive HPV infections; however, viral suppression through ART can help to mitigate this [11]. The model contains interacting compartments for all HPV susceptibility/infection and natural history states for HPV types 16/18, HPV H5 and HPV OHR. These stages are described in Table 3, and their interactions are summarised in Fig 2, which describes the state transitions possible from each state at the start of a new timestep, including the case where no state transition is made.

[Figure omitted. See PDF.]

Fig 2. State space diagram for the natural history of HPV and cervical cancer carcinogenesis; note that all compartments are subject to natural mortality, and detected cancer (grey) compartments are subject to stage-specific cervical cancer mortality and survival rates; CIN = cervical intraepithelial neoplasia; CC = cervical cancer.

https://doi.org/10.1371/journal.pone.0231388.g002

[Figure omitted. See PDF.]

Table 3. List and description of HPV transmission and natural history compartments.

https://doi.org/10.1371/journal.pone.0231388.t003

The model explicitly simulates the natural history of human papillomavirus infection for HPV types 16/18, HPV H5 and HPV OHR. The stage-, age- and HPV type-specific progression and regression rates are as published in previous analyses [39], and have been reproduced in S1 File. HIV positivity status and viral suppression through ART both impact HPV acquisition and natural history; assumptions regarding the impact of HIV positivity on HPV natural history are summarised in S1 File.

Interventions.

The model accounts for the impact of a range of HIV control interventions, including uptake of ART, PrEP, VMMC, and behavioural factors including the use of condoms. Effective use of ART in an individual infected with HIV not only acts to reduce disease progression and HIV-death but also significantly reduces the infectivity of virally suppressed patients. Further, use of PrEP, VMMC and condoms all lower the probability of disease acquisition to varying degrees. In the model, VMMC is specified by year, and we assumed rates were as reported in the literature and Tanzanian DHS reports [40–42]. The modelled VMMC rate applies to males of all ages, and reduces female to male HIV transmission by 60%, as consistent with the available evidence [43].

ART is also considered in two categories: those who are receiving ART, and those who are receiving ART and are ‘virally suppressed’. The model assumes some mortality benefit for all individuals receiving ART, with viral suppression completely halting disease progression and/or HIV death, and reducing infectiousness by 96% [13]. The percentage of virally suppressed people living with HIV (PLHIV) in Tanzania was assumed to match figures published by UNAIDS [44].

Scenarios and outcomes

A range of counterfactual and potential future HIV epidemic control scenarios were simulated, as described in Table 4. For each scenario, we estimated cervical cancer incidence and cervical cancer mortality (stratified by HIV positivity) from 1995–2020, projecting these outcomes based on model hypotheses from 2020–2070. The absolute numbers of cervical cancer cases and deaths prevented by interventions to date (2020) are presented in this analysis, in addition to an age-standardised rate (ASR). The age-standardised rate is a weighted mean of the age-specific rates where weights (summing to unity) are derived from the 2015 estimated world female population [35], presented per 100,000 women.

[Figure omitted. See PDF.]

Table 4. Modelled scenarios of HIV control.

https://doi.org/10.1371/journal.pone.0231388.t004

Sensitivity analysis

A multivariate sensitivity analysis was carried out to assess the robustness of model outcomes to variation in a range of parameters. Parameters were selected for sensitivity analysis if they were either difficult to observe/report on, suspected to be highly influential, or directly affect the interventions assessed in the scenario analysis. The modelled effect of HIV control interventions is dependent on assumptions about the magnitude of their effectiveness, and, the literature indicates uncertainty surrounding the effect of VMMC, ART and PrEP [11, 16, 46–48]. Furthermore, any population or individual level behavioural change driven by implementation of these interventions is difficult to quantify, as perceptions of risk are constantly changing. However, evidence suggests that the availability and uptake of HIV control interventions may facilitate an increase in risky sexual practices to the order of up to 21% [49–54]. A Latin Hypercube Sampling (LHS) analysis over 6,000 possible parameter sets was utilised, the values of which are described in Table 5.

[Figure omitted. See PDF.]

Table 5. Parameter variation considered in LHS analysis, and the rationale for selecting these parameters/ranges.

Parameters are described in S1 Table in S1 File.

https://doi.org/10.1371/journal.pone.0231388.t005

Results

Calibration and validation

Simulations from the calibrated model were consistent with observed HIV-specific outcomes including male and female HIV prevalence, total HIV incidence and number of HIV deaths (Fig 3), in addition to age-specific 2018 cervical cancer incidence and mortality rates (Fig 4).

[Figure omitted. See PDF.]

Fig 3. Calibrated HIV outcomes.

(A) and (B) male and female HIV prevalence from 1995 to 2015; (C) HIV incidence from 1995 to 2015; (D) number of HIV deaths from 1995 to 2015. Error bars are 95% CI of observed data. Training data sourced from UNAIDS [22, 30,31,57].

https://doi.org/10.1371/journal.pone.0231388.g003

[Figure omitted. See PDF.]

Fig 4. Calibrated (A) age-specific cervical cancer incidence and (B) mortality for the year 2018 compared to estimated data sourced from the International Agency for Research on Cancer IARC [23].

https://doi.org/10.1371/journal.pone.0231388.g004

Following the model calibration procedure, where parameters were chosen such that the model was a good fit to UNAIDS and Globocan (IARC) data [22, 23, 30, 31], the model was validated against independent datasets. These included sex- and age-specific HIV prevalence (Fig 5A and 5B) and the sex-specific age distribution of AIDS diagnoses (Fig 5C and 5D), age-specific HPV prevalence (Fig 6), and the prevalence of HSIL (high-grade squamous intra-epithelial lesion; considered equivalent to a diagnosed CIN 2/3) among HIV negative versus positive women from the PROTECT study [58] (Fig 7).

[Figure omitted. See PDF.]

Fig 5. (A) and (B) Age-specific HIV prevalence for males and females in 2016 compared to observed data; (C) and (D: age-distribution of AIDS diagnoses for males and females in 2011 compared to observed data.

Observed data from the Tanzanian Ministry of Finance (no confidence intervals available) [26, 59].

https://doi.org/10.1371/journal.pone.0231388.g005

[Figure omitted. See PDF.]

Fig 6. (A), (B) and (C) Age-specific HPV prevalence in cervical cytology (all cytological results) of HPV 16/18, HPV H5 and HPV OHR compared to observed data.

Observed data from Dartell et al 2012 (no confidence intervals available) [58].

https://doi.org/10.1371/journal.pone.0231388.g006

[Figure omitted. See PDF.]

Fig 7. Age specific rates of HSIL (detected high-grade squamous intraepithelial lesion consistent with CIN2/3) prevalence for (A) HIV positive and (B) HIV negative women.

Observed data from Dartell et al 2012 (no confidence intervals available) [58].

https://doi.org/10.1371/journal.pone.0231388.g007

Scenario analysis

The model estimates that there were 2,843 (and 1,039) fewer cervical cancer cases (and deaths), and 33,648 fewer total deaths (HIV and cervical cancer combined) from 1995 to 2020 as a result of the introduction and scale-up of VMMC. Assuming VMMC is maintained at 2018 levels, by 2070, VMMC is expected to have averted 330,400 (and 186,260) cervical cancer cases (and deaths) and save 3.47 million lives (HIV and cervical cancer combined) (Table 6; Fig 8B and 8D). ART use to 2020 is estimated to have led to an additional 1,573 (and 1,222) cervical cancer cases (and deaths) from 1995–2020; this is due to ART reducing the competing risk of mortality due to HIV. The cumulative number of deaths averted by ART is predicted to reach 2,253,700 by 2070 and will have prevented 52,430 (and 16,390) additional cervical cancer cases (and deaths) by 2070 (Table 6; Fig 8B and 8D).

[Figure omitted. See PDF.]

Fig 8. (A) Annual cervical cancer cases averted due to VMMC and ART (negative values under ART denote additional cases rather than cases averted); (B) cumulative cervical cancer cases averted due to VMMC and ART; (C) annual cervical cancer deaths averted due to VMMC and ART; (D) cumulative cervical cancer deaths averted due to VMMC and ART.

https://doi.org/10.1371/journal.pone.0231388.g008

[Figure omitted. See PDF.]

Table 6. Number of prevented cervical cancer cases (and deaths) due to VMMC and ART.

Note that values presented reflect the incremental benefit of each intervention, and that negative values denote additional cases/deaths.

https://doi.org/10.1371/journal.pone.0231388.t006

VMMC is predicted to substantially reduce both HPV and HIV prevalence in Tanzanian men and women (Table 7, Fig 9). By 2070, VMMC is predicted to reduce HPV prevalence in men and women by 28% (44.71% under ‘No interventions’ to 32.13% under ‘VMMC only’ in 2070) and 17% (46.39% under ‘No interventions’ to 38.56% under ‘VMMC only’), respectively. Similarly, the estimated reduction in HIV prevalence due to VMMC is 75% (7.59% under ‘No interventions’ to 1.86% under ‘VMMC only’) and 71% (8.47% under ‘No interventions’ to 2.46% under ‘VMMC only’), respectively, for men and women. Compared to the ‘No intervention’ scenario, the introduction and scale-up of ART and PrEP is estimated to reduce HIV prevalence in men and women by approximately 99% (to 0.05% in males and 0.11% in females) in 2070.

[Figure omitted. See PDF.]

Fig 9. (A) and (B) Male and female HPV prevalence; (C) and (D) male and female HIV prevalence from 1995 to 2070 under five intervention scenarios.

https://doi.org/10.1371/journal.pone.0231388.g009

[Figure omitted. See PDF.]

Table 7. Simulated HPV and HIV prevalence for males and females in 2070 (and 2020) under five scenarios.

Note that intervention start-year occurs pre-2020 for simulated interventions.

https://doi.org/10.1371/journal.pone.0231388.t007

The introduction and scale-up of HIV preventions is expected to reduce the age-standardised rates of cervical cancer incidence and mortality over time. VMMC is expected to reduce cervical cancer incidence and mortality rates by 36–40% in 2070 compared to 2020 rates (under ‘VMMC and ART (baseline)’ scenario), whereas the provision and scale-up of ART to meet World Health Organisation 90-90-90 HIV control targets reduces cervical cancer incidence and mortality rates by 41–45% in 2070 compared to the current rates in 2020. Absolute rates of cervical cancer incidence and mortality are presented in Table 8 and visualised in Fig 10. Here, we note that the reductions in cervical cancer incidence and mortality due to ART among all women are driven by reductions among HIV positive women.

[Figure omitted. See PDF.]

Fig 10. (A) and (B) Age-standardised cervical cancer incidence and mortality rates among all women aged 0–99 years; (C) and (D) cervical cancer incidence and mortality rates among HIV negative women aged 0–99 years; (E) and (F) cervical cancer incidence and mortality rates among HIV positive women aged 0–99 years.

Age-standardised rates are calculated using the 2015 World Female Population [35].

https://doi.org/10.1371/journal.pone.0231388.g010

[Figure omitted. See PDF.]

Table 8. Simulated age-standardised cervical cancer incidence and mortality rates per 100,000 women, for all women (and stratified by HIV positivity) in 2070 (and 2020) under five scenarios.

Rates are standardised to the WFP2015 population. [35]

https://doi.org/10.1371/journal.pone.0231388.t008

Sensitivity analysis

Findings from the multivariate sensitivity analysis indicate that the simulation outcomes are highly sensitive to variation in parameters specifying sexual behaviour, disease transmission and natural history, and, intervention effectiveness. Fig 11 summarises the baseline and total variation in 2070 endpoint predictions for all simulated outcomes over the five scenarios.

[Figure omitted. See PDF.]

Fig 11. (A) and (B) Male and female HPV prevalence; (C) and (D) male and female HIV prevalence; (E) and (F) cervical cancer incidence and mortality among all women; (G) and (H) cervical cancer incidence and mortality among HIV negative women; (I) and (J) cervical cancer incidence and mortality among HIV positive women, simulated in the year 2070 (error bars correspond to the total variation generated by the sensitivity analysis).

https://doi.org/10.1371/journal.pone.0231388.g011

An analysis of partial rank correlation coefficients indicate that cervical cancer incidence is most strongly correlated with VMMC efficacy for HPV prevention (correlation coefficient of -0.25; Fig 12), followed by VMMC efficacy for HIV prevention (correlation coefficient of -0.13).

[Figure omitted. See PDF.]

Fig 12. Correlation strength of selected outputs (HIV and HPV prevalence for males and females, and cervical cancer incidence and mortality) against intervention and behavioural parameters varied in multivariate sensitivity analysis.

Partial rank correlation analysis was performed on the ‘VMMC, target ART and PrEP’ scenario, as it is the only scenario considering all modelled interventions.

https://doi.org/10.1371/journal.pone.0231388.g012

Discussion

To our knowledge, this analysis is the first to directly estimate the impact over time of changing VMMC prevalence, ART utilisation and PrEP uptake on cervical cancer in any setting. Findings from this analysis are likely to be broadly applicable to other low-income settings with high HIV prevalence and cervical cancer incidence rates, particularly in sub-Saharan Africa. These HIV control interventions were found to have a substantial impact on cervical cancer incidence and mortality in Tanzania. Using a simulation model of HIV and HPV transmission to estimate cervical cancer cases, cervical cancer deaths and total deaths (including HIV deaths) in Tanzania from 1995 to 2070 in the context of currently implemented HIV control measures, we estimated that VMMC has prevented 2,843 cervical cancer cases and 1,039 cervical cancer deaths from 1995 to 2020. Perhaps a less intuitive finding is that, while the addition and scale-up of ART in HIV-positive women reduces both overall HIV and HPV prevalence (women effectively treated with ART are less likely to acquire HPV and more likely to clear an HPV infection than their untreated counterparts), ART is estimated to have resulted in some 1,573 additional cervical cancer diagnoses and 1,221 additional cervical cancer deaths, cumulatively, from 1995 to 2020. These additional cases and deaths are among HIV infected women and are caused by the removal of HIV-related death as a competing risk, as some women who would have otherwise died from HIV-related causes will develop cervical cancer and subsequently die from it, in the absence of scaled-up cervical cancer prevention. In the longer-term, the protective effect of ART prevails, as scale-up to meet World Health Organisation 90-90-90 HIV control targets would result in cervical cancer incidence and mortality rates that are 43% (37.31 c.f. 65.42 cases per 100,000 women per year) and 39% (26.44 c.f. 43.27 deaths per 100,000 women per year) lower, respectively, in 2070 compared to the current rates in 2020.

The model prediction for HIV prevalence over time is consistent with empirical data for Tanzania [22, 30], and, future predictions are broadly consistent with the findings from HIV modelling studies specific to sub-saharan Africa [60]. A comparative modelling study utilising predictions from four independent models predicts that under a scenario assuming universal HIV testing and treatment (up to 90% coverage), HIV prevalance will be reduced to 0–3% in sub-Saharan Africa in 2050. The current analysis predicts that HIV prevalence in Tanzania in 2050 will be 0.12% in males, and 0.19% in females aged 15–49 years in the ‘VMMC and target ART’ scenario. In addition to this, predicted HPV prevalence in males (35.4% in 2017) under the ‘VMMC and ART (baseline)’ scenario showed relatively consistent agreement with observed HPV prevalence in South Africa in 2017 (40%) [61], and, predictions published by Tan et al [13]. Similarly, Tan et al predict that the cervical cancer incidence rates in KwaZulu Natal women will be approximately 31 cases per 100,000 among HIV negative women in 2070, and approximately 145 cases per 100,000 among HIV positive women (S1 File) [13]. For Tanzania, our analysis predicts 35.17 cases per 100,000 women among HIV negative women, and 220.23 cases per 100,000 among HIV positive women under the ‘VMMC and ART (baseline)’ scenario in 2070. In both analyses, the cervical cancer incidence rate in HIV positive women is close to five times higher than the rate among HIV negative women.

The dynamic and highly detailed nature of the model of HIV and HPV co-infection is a strength of this study. The model is stratified by sex, age, sexual activity level (including women engaging in transactional and commercial sex), HIV positivity (including disease stage and treatment status) and HPV positivity (including disease stage/detection status) for multiple HPV types; this allows exploration of disease transmission and progression dynamics in detail, accounting for herd protective effects and the protective effects of ART.

This analysis was limited by the inherent uncertainty surrounding input parameter assumptions, in particular, sexual behaviour assumptions including condom usage. In many settings which have implemented such HIV controls, a reduction in safe sex practices and an increase in other sexually transmitted infections is often observed [49–54]. A recent study into the sexual behaviour of PrEP users in Amsterdam found that daily PrEP use among HIV negative men who have sex with men (MSM) was associated with a 2–9% increase in condomless sex acts [50]; whereas another study has reported a 21% increase in risky sexual practice, and an increase in HIV incidence, among the San Franciscan MSM population since the advent of ART [53]. Sensitivity analysis findings indicate that HPV and HIV prevalence, and cervical cancer incidence and mortality are highly sensitive to variations in condom usage, therefore if condom usage trends over time vary, model predictions could substantially under-estimate or over-estimate disease burden. A number of simplifying assumptions were also made regarding the effect of ART on HPV natural history in HIV positive women; namely, we assume that ART affects all HPV types to the same degree, and, that commencing ART has the same effect on HPV natural history regardless of CD4+ count or HIV disease stage. Furthermore, the HIV transmission component of this model accounts for heterosexual transmission only, which is based on the assumption that the impact of HIV transmitted via sexual contact between men, and injection drug use, will have negligible impacts on the cervical cancer in women. Cervical cancer is an AIDS defining disease; therefore, there is a degree of uncertainty surrounding the true cervical cancer mortality rates in Tanzania, that is, whether a cervical cancer death in a HIV positive woman was attributed to cervical cancer or HIV [62]. The grouping of many individual HPV genotypes in the model, for example HPV 16/19, HPV 31/33/45/52/58 (HPV H5) and a category for “other high-risk HPV types”, may impact the overall simulated HPV prevalence in addition to the overall transmission dynamics of the model; for example, in this model it is impossible to discern whether individuals are infected with only one or any combination of the HPV genotypes in each simulated HPV subgroup. This may result in an overestimation of effectiveness of interventions targeted at HPV reduction. Finally, the findings of this study must be interpreted in the context of the lengthening life expectency in Tanzania, which is explicitely accounted for in this analysis. Due to reductions in HIV mortality in addition to other cause mortality (e.g. driven by improvements in sanitation and health care), the life expectancy at birth is expected to rise to 75 years in 2065–2070 (compared to 54 years in 1995–2000) [34]. This will necessarily result in an increased opportunity for the development of cervical cancer (and other diseases), irrespective of additional effects due to HIV treatment.

In 2019, a draft global strategy for the elimination of cervical cancer as a public health problem was released by the World Health Organisation [63]. This strategy, to be considered by the World Health Assembly in May 2020, defines that cervical cancer is eliminated as a public health problem when all countries achieve an incidence rate of less than four cases per 100 000 women per year. To achieve this target, the WHO recommends that each country implement HPV vaccination programmes whereby 90% of girls are vaccinated by the age of 15, organised cervical screening programmes whereby 70% of women are screened at least twice per lifetime, and effective management of 90% of women diagnosed with cervical pre-cancer or invasive cervical cancer [63]. While VMMC and ART can reduce the burden of cervical cancer in Tanzania in the long term, they are not sufficient to bring cervical cancer incidence beneath the threshold proposed by the WHO for cervical cancer elimination. Our finding that even under the best-case scenario the rate of cervical cancer incidence in all Tanzanian women is not reduced below 35 cases per 100,000 women per year (more than eight fold higher than the elimination threshold) demonstrates the importance and urgency of scaling up cervical cancer prevention programs, such as HPV vaccination and cervical screening, as well as HIV control, in order to avoid the situation that lives saved from HIV-related death are instead lost to cervical cancer. The WHO call for global action to eliminate cervical cancer as a public health problem is an important opportunity to galvanise and unite efforts to prevent cervical cancer in Tanzania and globally [64].

Supporting information

[Figure omitted. See PDF.]

S1 File.

https://doi.org/10.1371/journal.pone.0231388.s001

(PDF)

Acknowledgments

This research includes computations using the computational cluster Katana supported by Research Technology Services at UNSW Sydney. We also acknowledge the National Cancer Institute–funded Cancer Intervention and Surveillance Modeling Network (CISNET) cervical cancer working group for intellectual support and feedback throughout the project.

Citation: Hall MT, Smith MA, Simms KT, Barnabas RV, Canfell K, Murray JM (2020) The past, present and future impact of HIV prevention and control on HPV and cervical disease in Tanzania: A modelling study. PLoS ONE 15(5): e0231388. https://doi.org/10.1371/journal.pone.0231388

References

1. endinghiv.org. Ending HIV 2020: ACON; 2019 [https://endinghiv.org.au/#.

2. Joint United Nations Programme on HIV/AIDS (UNAIDS). 90-90-90 An ambition treatment to help end the AIDS epidemic2014 27/06/2019. https://www.unaids.org/sites/default/files/media_asset/90-90-90_en.pdf.

3. Joint United Nations Programme on HIV/AIDS (UNAIDS). Making decisions on male circumcision for HIV risk reduction: modelling the impact and costs. 2009 27/06/2019.

4. PrEP Watch. National Policies and Guidelines for PrEP 2019 [27/06/2019]. https://www.prepwatch.org/prep-planning/national-policies-guidelines/.

5. Mishra S, Steen R, Gerbase A, Lo YR, Boily MC. Impact of high-risk sex and focused interventions in heterosexual HIV epidemics: a systematic review of mathematical models. PLoS One. 2012;7(11):e50691. pmid:23226357

6. Blower S, Schwartz EJ, Mills J. Forecasting the future of HIV epidemics: the impact of antiretroviral therapies & imperfect vaccines. AIDS Rev. 2003;5(2):113–25. pmid:12876900

7. Van der Ploeg CPB, Van Vliet C, De Vlas SJ, Ndinya-Achola JO, Fransen L, van Oortmarssen GJ, et al. STDSIM: a microsimulation model for decision support in STD control. Interfaces. 1998;28(3):84–100.

8. Vissers DC, DEV SJ, Bakker R, Urassa M, Voeten HA, Habbema JD. The impact of mobility on HIV control: a modelling study. Epidemiol Infect. 2011;139(12):1845–53. pmid:21299914

9. Kaldor JM, Wilson DP. How low can you go: the impact of a modestly effective HIV vaccine compared with male circumcision. AIDS (London, England). 2010;24(16):2573–8.

10. Herbeck JT, Mittler JE, Gottlieb GS, Goodreau SM, Murphy JT, Cori A, et al. Evolution of HIV virulence in response to widespread scale up of antiretroviral therapy: a modeling study. Virus Evol. 2016;2(2):vew028. pmid:29492277

11. Liu G, Sharma M, Tan N, Barnabas RV. HIV-positive women have higher risk of human papilloma virus infection, precancerous lesions, and cervical cancer. AIDS (London, England). 2018;32(6):795–808.

12. Campos NG, Lince-Deroche N, Chibwesha CJ, Firnhaber C, Smith JS, Michelow P, et al. Cost-Effectiveness of Cervical Cancer Screening in Women Living With HIV in South Africa: A Mathematical Modeling Study. J Acquir Immune Defic Syndr. 2018;79(2):195–205. pmid:29916959

13. Tan N, Sharma M, Winer R, Galloway D, Rees H, Barnabas RV. Model-estimated effectiveness of single dose 9-valent HPV vaccination for HIV-positive and HIV-negative females in South Africa. Vaccine. 2018;36(32 Pt A):4830–6. pmid:29891348

14. Yuan T, Fitzpatrick T, Ko NY, Cai Y, Chen Y, Zhao J, et al. Circumcision to prevent HIV and other sexually transmitted infections in men who have sex with men: a systematic review and meta-analysis of global data. Lancet Glob Health. 2019;7(4):e436–e47. pmid:30879508

15. Prodger JL, Kaul R. The biology of how circumcision reduces HIV susceptibility: broader implications for the prevention field. AIDS Research and Therapy. 2017;14(1):49. pmid:28893286

16. Castellsague X, Bosch FX, Munoz N, Meijer CJ, Shah KV, de Sanjose S, et al. Male circumcision, penile human papillomavirus infection, and cervical cancer in female partners. The New England journal of medicine. 2002;346(15):1105–12. pmid:11948269

17. Gray RH, Kigozi G, Serwadda D, Makumbi F, Watya S, Nalugoda F, et al. Male circumcision for HIV prevention in men in Rakai, Uganda: a randomised trial. Lancet. 2007;369(9562):657–66. pmid:17321311

18. Bailey RC, Moses S, Parker CB, Agot K, Maclean I, Krieger JN, et al. Male circumcision for HIV prevention in young men in Kisumu, Kenya: a randomised controlled trial. Lancet. 2007;369(9562):643–56. pmid:17321310

19. Auvert B, Taljaard D, Lagarde E, Sobngwi-Tambekou J, Sitta R, Puren A. Randomized, Controlled Intervention Trial of Male Circumcision for Reduction of HIV Infection Risk: The ANRS 1265 Trial. PLOS Medicine. 2005;2(11):e298. pmid:16231970

20. Grund JM, Bryant TS, Jackson I, Curran K, Bock N, Toledo C, et al. Association between male circumcision and women’s biomedical health outcomes: a systematic review. The Lancet Global Health. 2017;5(11):e1113–e22. pmid:29025633

21. Drain PK, Halperin DT, Hughes JP, Klausner JD, Bailey RC. Male circumcision, religion, and infectious diseases: an ecologic analysis of 118 developing countries. BMC Infectious Diseases. 2006;6(1):172.

22. UNAIDS. Data Sheet—HIV Prevalence Population: Female adults (15–49) Geneva, Switzerland2018 [http://aidsinfo.unaids.org?did=554ca08005b5b281750af20a&r=world&t=2016&tb=d&bt=dnli&ts=0,0&tr=world&aid=577bfcca880284565dfd2ebf&sav=Population:Femaleadults(15–49)&tl=2.

23. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424. pmid:30207593

24. UNAIDS Joint United Nations Programme on HIV/AIDS. Miles to go: the response to HIV in Eastern and Southern Africa [https://www.unaids.org/sites/default/files/media_asset/miles-to-go_eastern-and-southern-africa_en.pdf.

25. AVERT. HIV AND AIDS IN EAST AND SOUTHERN AFRICA REGIONAL OVERVIEW United Kingdom: AVERT; 2019 [https://www.avert.org/professionals/hiv-around-world/sub-saharan-africa/overview.

26. Ministry of Health CD, Gender, Elderly and Children (MoHCDGEC), Government of Tanzania, Ministry of Health Zanzibar (MoH), National Bureau of Statistics (NBS), Office of Chief Government Statistician (OCGS). Tanzania HIV Impact Survey (THIS) 2016–2017. Dar es Salaam, Tanzania; 2017.

27. Ministry of Health CD, Gender, Elderly and Children (MoHCDGEC),. National guidelines for the management of HIV and AIDS. 2019.

28. Kharsany ABM, Karim QA. HIV Infection and AIDS in Sub-Saharan Africa: Current Status, Challenges and Opportunities. Open AIDS J. 2016;10:34–48. pmid:27347270

29. Mbulawa ZZA, van Schalkwyk C, Hu N-C, Meiring TL, Barnabas S, Dabee S, et al. High human papillomavirus (HPV) prevalence in South African adolescents and young women encourages expanded HPV vaccination campaigns. PloS one. 2018;13(1):e0190166-e.

30. UNAIDS. Data Sheet—HIV Prevalence Population: Male adults (15–49) Geneva, Switzerland: UNAIDS; 2018 [http://aidsinfo.unaids.org?did=554ca08005b5b281750af20a&r=world&t=2016&tb=d&bt=dnli&ts=0,0&tr=world&aid=577bfcca880284565dfd2ebf&sav=Population: Male adults (15–49)&tl = 2.

31. UNAIDS. Data Sheet—HIV Incidence Per 1000 Population Geneva, Switzerland2018 [http://aidsinfo.unaids.org?did=55da49cd64e925b94e70b0ce&r=world&t=2016&tb=d&bt=dnli&ts=0,0&tr=world&aid=5970eccef7341ed11f26de5d&sav=Population:Adults(15–49)&tl=2.

32. The World Bank. Fertility rate, total (births per woman): WORLD BANK GROUP; 2017 [https://data.worldbank.org/indicator/SP.DYN.TFRT.IN?locations=TZ.

33. National Bureau of Statistics. Mortality and Health. Dar es Salaam, Tanzania; 2015 July 2015.

34. United Nations Department of Economic and Social Affairs Population Division. World Population Prospects 2019, Online Edition. 2019 [https://population.un.org/wpp/Download/Standard/Population/.

35. United Nations DESA / Population Division. World Population Prospects: The 2015 Revision, DVD Edition. 2017 [https://esa.un.org/unpd/wpp/Download/Standard/Population/.

36. The World Bank. Sex ratio at birth (male births per female births) 2017 [https://data.worldbank.org/indicator/SP.POP.BRTH.MF?end=2017&locations=TZ&start=2017&view=bar.

37. Weinberg JL, Kovarik CL. The WHO Clinical Staging System for HIV/AIDS. American Medical Association Journal of Ethics. 2010;12(3):202–6. pmid:23140869

38. Palk L, Gerstoft J, Obel N, Blower S. A modeling study of the Danish HIV epidemic in men who have sex with men: travel, pre-exposure prophylaxis and elimination. Scientific reports. 2018;8(1):16003. pmid:30375426

39. Lew J-B, Simms KT, Smith MA, Hall M, Kang Y-J, Xu XM, et al. Primary HPV testing versus cytology-based cervical screening in women in Australia vaccinated for HPV and unvaccinated: effectiveness and economic assessment for the National Cervical Screening Program. The Lancet Public Health. 2017;2(2):e96–e107. pmid:29253402

40. The United Republic of Tanzania Ministry of Health and Social Welfare. Report Number 23. Dar es Salaam, Tanzania: National AIDS Control Programme; 2013.

41. Forbes HJ, Doyle AM, Maganja K, Changalucha J, Weiss HA, Ross DA, et al. Rapid Increase in Prevalence of Male Circumcision in Rural Tanzania in the Absence of a Promotional Campaign. PLOS ONE. 2012;7(7):e40507. pmid:22792359

42. Ministry of Health CD, Gender, Elderly and Children (MoHCDGEC), Ministry of Health (MoH), National Bureau of Statistics (NBS), Office of the Chief Government Statistician (OCGS), ICF. Tanzania Demographic and Health Survey and Malaria Indicator Survey 2015–16. 2016.

43. World Health Organization. Male circumcision for HIV prevention Geneva, Switzerland [20/01/2019]. https://www.who.int/hiv/topics/malecircumcision/en/.

44. UNAIDS. Data Sheet—PEOPLE LIVING WITH HIV RECIEVING ART (%) Geneva, Switzerland: UNAIDS; 2018 [http://aidsinfo.unaids.org?did=554c8ceac904ab29709ecd3b&r=world&t=null&tb=d&bt=dnli&ts=0,0&tr=world&aid=5ae071a462abc329969a8de1&sav=Population:Allages&tl=2.

45. UNAIDS. 90–90–90—An ambitious treatment target to help end the AIDS epidemic 2017 [https://www.unaids.org/en/resources/documents/2017/90-90-90.

46. Tobian AAR, Kacker S, Quinn TC. Male circumcision: a globally relevant but under-utilized method for the prevention of HIV and other sexually transmitted infections. Annu Rev Med. 2014;65:293–306. pmid:24111891

47. Minkoff H, Zhong Y, Burk RD, Palefsky JM, Xue X, Watts DH, et al. Influence of adherent and effective antiretroviral therapy use on human papillomavirus infection and squamous intraepithelial lesions in human immunodeficiency virus-positive women. J Infect Dis. 2010;201(5):681–90. pmid:20105077

48. Anderson PL, Glidden DV, Liu A, Buchbinder S, Lama JR, Guanira JV, et al. Emtricitabine-Tenofovir Concentrations and Pre-Exposure Prophylaxis Efficacy in Men Who Have Sex with Men. Science Translational Medicine. 2012;4(151):151ra25.

49. Ramchandani MS, Golden MR. Confronting Rising STIs in the Era of PrEP and Treatment as Prevention. Curr HIV/AIDS Rep. 2019;16(3):244–56. pmid:31183609

50. Hoornenborg E, Coyer L, Achterbergh RCA, Matser A, Schim van der Loeff MF, Boyd A, et al. Sexual behaviour and incidence of HIV and sexually transmitted infections among men who have sex with men using daily and event-driven pre-exposure prophylaxis in AMPrEP: 2 year results from a demonstration study. Lancet HIV. 2019.

51. Chen SY, Gibson S, Katz MH, Klausner JD, Dilley JW, Schwarcz SK, et al. Continuing increases in sexual risk behavior and sexually transmitted diseases among men who have sex with men: San Francisco, Calif, 1999–2001, USA. Am J Public Health. 2002;92(9):1387–8. pmid:12197957

52. Dukers NH, Goudsmit J, de Wit JB, Prins M, Weverling GJ, Coutinho RA. Sexual risk behaviour relates to the virological and immunological improvements during highly active antiretroviral therapy in HIV-1 infection. AIDS (London, England). 2001;15(3):369–78.

53. Katz MH, Schwarcz SK, Kellogg TA, Klausner JD, Dilley JW, Gibson S, et al. Impact of highly active antiretroviral treatment on HIV seroincidence among men who have sex with men: San Francisco. Am J Public Health. 2002;92(3):388–94. pmid:11867317

54. Scheer S, Chu PL, Klausner JD, Katz MH, Schwarcz SK. Effect of highly active antiretroviral therapy on diagnoses of sexually transmitted diseases in people with AIDS. Lancet. 2001;357(9254):432–5. pmid:11273063

55. Munguti K, Grosskurth H, Newell J, Senkoro K, Mosha F, Todd J, et al. Patterns of sexual behaviour in a rural population in north-western Tanzania. Soc Sci Med. 1997;44(10):1553–61. pmid:9160444

56. Quinn TC, Wawer MJ, Sewankambo N, Serwadda D, Li C, Wabwire-Mangen F, et al. Viral load and heterosexual transmission of human immunodeficiency virus type 1. Rakai Project Study Group. The New England journal of medicine. 2000;342(13):921–9. pmid:10738050

57. UNAIDS. UNITED REPUBLIC OF TANZANIA 2018 [http://aidsinfo.unaids.org/.

58. Dartell M, Rasch V, Kahesa C, Mwaiselage J, Ngoma T, Junge J, et al. Human papillomavirus prevalence and type distribution in 3603 HIV-positive and HIV-negative women in the general population of Tanzania: the PROTECT study. Sexually transmitted diseases. 2012;39(3):201–8. pmid:22337107

59. Tanzania Commission for AIDS (TACAIDS), Zanzibar AIDS Commission (ZAC), National Bureau of Statistics (NBS), Office of Chief Government Statistician (OCGS), ICF International. Tanzania HIV/AIDS and Malaria Indicator Survey 2011–2012. Dar es Salaam, Tanzania; 2013.

60. Hontelez JAC, Lurie MN, Bärnighausen T, Bakker R, Baltussen R, Tanser F, et al. Elimination of HIV in South Africa through Expanded Access to Antiretroviral Therapy: A Model Comparison Study. PLOS Medicine. 2013;10(10):e1001534. pmid:24167449

61. Mbulawa ZZA, Wilkin T, Goeieman BJ, Jong E, Michelow P, Swarts A, et al. Prevalence of Anal Human Papillomavirus (HPV) and Performance of Cepheid Xpert and Hybrid Capture 2 (hc2) HPV Assays in South African HIV-Infected Women. Am J Clin Pathol. 2017;148(2):148–53. pmid:28898982

62. World Health Organization. Annex 1. WHO clinical staging of HIV disease in adults, adolescents and children. Geneva, Switzerland: WHO; 2013 [http://www.who.int/hiv/pub/guidelines/arv2013/annexes/WHO_CG_annex_1.pdf.

63. World Health Organisation (WHO). Draft: global strategy towards the elimination of cervical cancer as a public health problem Geneva, Switzerland2019 [https://www.who.int/docs/default-source/documents/cervical-cancer-elimination-draft-strategy.pdf.

64. World Health Organisation (WHO). WHO leads the way towards the elimination of cervical cancer as a public health concern 2018 [https://www.who.int/reproductivehealth/cervical-cancer-public-health-concern/en/.