Cross-sectional estimates revealed high HIV incidence in Botswana rural communities in the era of successful ART scale-up in 2013-2015

Abstract

Background

Botswana is close to reaching the UNAIDS “90-90-90” HIV testing, antiretroviral treatment (ART), and viral suppression goals. We sought to determine HIV incidence in this setting with both high HIV prevalence and high ART coverage.

Methods

We used a cross-sectional approach to assessing HIV incidence. A random, population-based sample of adults age 16–64 years was enrolled in 30 rural and peri-urban communities as part of the Botswana Combination Prevention Project (BCPP), from October 2013 –November 2015. Data and samples from the baseline household survey were used to estimate cross-sectional HIV incidence, following an algorithm that combined Limiting-Antigen Avidity Assay (LAg-Avidity EIA), ART status (documented or by testing ARV drugs in plasma) and HIV-1 RNA load. The LAg-Avidity EIA cut-off normalized optical density (ODn) was set at 1.5. The HIV-1 RNA cut-off was set at 400 copies/mL. For estimation purposes, the Mean Duration of Recent Infection was 130 days and the False Recent Rate (FRR) was evaluated at values of either 0 or 0.39%.

Results

Among 12,610 individuals participating in the baseline household survey, HIV status was available for 12,570 participants and 3,596 of them were HIV positive. LAg-Avidity EIA data was generated for 3,581 (99.6%) of HIV-positive participants. Of 326 participants with ODn ≤1.5, 278 individuals were receiving ART verified through documentation and were considered to represent longstanding HIV infections. Among the remaining 48 participants who reported no use of ART, 14 had an HIV-1 RNA load ≤400 copies/mL (including 3 participants with ARVs in plasma) and were excluded, as potential elite/viremic controllers or undisclosed ART. Thus, 34 LAg-Avidity-EIA-recent, ARV-naïve individuals with detectable HIV-1 RNA (>400 copies/mL) were classified as individuals with recent HIV infections. The annualized HIV incidence among 16–64 year old adults was estimated at 1.06% (95% CI 0.68–1.45%) with zero FRR, and at 0.64% (95% CI 0.24–1.04%) using a previously defined FRR of 0.39%. Within a subset of younger individuals 16–49 years old, the annualized HIV incidence was estimated at 1.29% (95% CI 0.82–1.77%) with zero FRR, and at 0.90% (95% CI 0.42–1.38%) with FRR set to 0.39%.

Conclusions

Using a cross-sectional estimate of HIV incidence from 2013–2015, we found that at the time of near achievement of the UNAIDS 90-90-90 targets, ~1% of adults (age 16–64 years) in Botswana’s rural and peri-urban communities became HIV infected annually.

About the Authors:

Sikhulile Moyo

Roles Conceptualization, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliations Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana, Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America

ORCID logo http://orcid.org/0000-0003-3821-4592

Simani Gaseitsiwe

Roles Formal analysis, Investigation, Project administration, Supervision, Writing – review & editing

Affiliations Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana, Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America

Terence Mohammed

Roles Investigation, Supervision, Writing – review & editing

Affiliation: Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana

Molly Pretorius Holme

Roles Conceptualization, Investigation, Supervision, Writing – review & editing

Affiliations Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana, Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America

Rui Wang

Roles Investigation, Methodology, Writing – review & editing

Affiliations Department of Population Medicine, Harvard Pilgrim Health Care Institute, Boston, Massachusetts, United States of America, Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston Massachusetts, United States of America, Harvard Medical School, Boston, Massachusetts, United States of America

Kenanao Peggy Kotokwe

Roles Investigation, Writing – review & editing

Affiliation: Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana

Corretah Boleo

Roles Investigation, Writing – review & editing

Affiliation: Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana

Lucy Mupfumi

Roles Investigation, Supervision, Writing – review & editing

Affiliation: Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana

Etienne Kadima Yankinda

Roles Investigation, Project administration, Supervision, Writing – review & editing

Affiliation: Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana

Unoda Chakalisa

Roles Investigation, Project administration, Supervision, Writing – review & editing

Affiliation: Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana

Erik van Widenfelt

Roles Investigation, Supervision, Writing – review & editing

Affiliation: Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana

Tendani Gaolathe

Roles Investigation, Project administration, Supervision, Writing – review & editing

Affiliation: Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana

Mompati O. Mmalane

Roles Investigation, Project administration, Supervision, Writing – review & editing

Affiliation: Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana

Scott Dryden-Peterson

Roles Investigation, Writing – review & editing

Affiliations Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana, Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America, Division of Infectious Diseases, Brigham and Women's Hospital, Boston, Massachusetts, United States of America

Madisa Mine

Roles Investigation, Writing – review & editing

Affiliation: Botswana Ministry of Health and Wellness, Gaborone, Botswana

Refeletswe Lebelonyane

Roles Investigation, Methodology, Project administration, Writing – review & editing

Affiliation: Botswana Ministry of Health and Wellness, Gaborone, Botswana

Kara Bennett

Roles Writing – review & editing

Affiliation: Bennett Statistical Consulting, Inc., Ballston Lake, New York, United States of America

ORCID logo http://orcid.org/0000-0002-9465-8983

Jean Leidner

Roles Writing – review & editing

Affiliation: Goodtables Data Consulting, Norman, OK, United States of America

Kathleen E. Wirth

Roles Investigation, Methodology, Supervision, Writing – review & editing

Affiliations Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America, Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America

Eric Tchetgen Tchetgen

Roles Conceptualization, Supervision, Writing – review & editing

Affiliations Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston Massachusetts, United States of America, Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America

Kathleen Powis

Roles Supervision, Writing – review & editing

Affiliations Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana, Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America, Harvard Medical School, Boston, Massachusetts, United States of America, Departments of Internal Medicine and Pediatrics, Massachusetts General Hospital, Boston, Massachusetts, United States of America

Janet Moore

Roles Investigation, Project administration, Supervision, Writing – review & editing

Affiliation: U.S. Centers for Disease Control, Atlanta, Georgia, United States of America

William A. Clarke

Roles Methodology

Affiliation: Johns Hopkins Hospital, Baltimore, MD, United States of America

Shahin Lockman

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

Affiliations Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana, Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America, Harvard Medical School, Boston, Massachusetts, United States of America, Division of Infectious Diseases, Brigham and Women's Hospital, Boston, Massachusetts, United States of America

Joseph M. Makhema

Roles Conceptualization, Investigation, Methodology, Project administration, Supervision, Writing – review & editing

Affiliations Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana, Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America

Max Essex

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

Affiliations Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana, Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America

Vlad Novitsky

Roles Conceptualization, Formal analysis, Investigation, Methodology, Supervision, Writing – original draft, Writing – review & editing

Affiliations Botswana Harvard AIDS Institute Partnership, Gaborone, Botswana, Department of Immunology and Infectious Diseases, Harvard T.H. Chan School of Public Health, Boston, Massachusetts, United States of America

Introduction

Botswana has been hard hit by the HIV-epidemic, with the third highest HIV prevalence worldwide among adults age 15–49, after Lesotho and Swaziland [1]. Botswana appears to be approaching the UNAIDS “90-90-90” HIV testing, treatment, and viral suppression targets [2]. These high levels of coverage have led to significant reductions in HIV-related mortality [1, 3–5]. In June 2016 Botswana adopted the World Health Organization (WHO) recommendation to provide Universal Test and Treat (UTT) [6]. The success of UTT could be measured by reduction in HIV incidence [7–25]. Monitoring of HIV incidence is a critical tool for assessment and evaluation the impact of HIV prevention and treatment programs.

Prospective longitudinal cohorts remain the gold standard for assessing HIV incidence. However, this approach is time consuming, costly and prone to selection and observational biases [8, 26]. Biomarkers of recent HIV infection that can be detected in cross-sectional samples represent a viable alternative to longitudinal cohort studies. Serological and molecular biomarkers could be combined in multi-assay algorithm (MAA). An optimized MAA with high sensitivity and specificity can discriminate between recent and established HIV infections in cross-sectional sample [27–29][8, 9, 25, 30–36]. Recent advances in design and development of MAA have facilitated estimating of HIV incidence in cross-sectional surveys with improved accuracy [7, 16, 22, 37].

In this study, we estimated HIV incidence using a baseline cross-sectional sample from the Botswana Combination Prevention Project (BCPP; the Ya Tsie study) [2, 38, 39]. BCPP is an ongoing pair-matched, cluster-randomized clinical trial in 30 rural and peri-urban communities across Botswana. The primary question of the study is reduction in the cumulative HIV incidence as a result of combination prevention interventions that included enhanced HIV testing and counseling (HTC) campaigns, linkage to care, antiretroviral treatment (ART), strengthened male circumcision (MC) and enhanced prevention of mother-to-child transmission of HIV. Between October 2013 and November 2015, we selected a random sample from 30 communities in three main geographic areas in Botswana: (1) south east, (2) north east, and (3) central eastern region. The communities were purposively selected and proposed based on (1) desired size, and (2) feasibility. Pairs of communities were matched by size, health services, population age structure, and geographic location. In each community, a complete list of all household-like structures (located within the prespecified community boundaries as defined by the 2011 Botswana Census) was obtained and geocoded using satellite imagery (Google Earth, Mountain View, CA, USA). Based on these lists, a simple random sample of approximately 20% of all households was drawn. At each selected household all household members were enumerated, assessed for eligibility, and approached for participation [2]. The vast majority of HIV-positive participants of the baseline household survey, 83%, knew their HIV status [2], 87% of them were receiving ART [2], and 96% of those on ART were virologically suppressed [2].

Materials and methods

Study participants

Blood specimens were collected during the BCPP baseline household survey. The HIV-positive status of participants was based on either written documentation provided (e.g., HIV test results, ART prescription) or HIV testing that was performed in the households according to the Botswana national guidelines by using double positive rapid HIV testing. Participants who self-reported not being on ART and classified as recently infected by the MAA were tested for presence of ARV drugs in their plasma. In addition to HTC, the survey staff provided point-of-care CD4 testing, collected blood from people living with HIV for viral load testing and viral genotyping (venous blood was collected by phlebotomy in households), evaluated uptake of HTC, and assessed ART and MC coverage.

Among 12,610 individuals participating in the baseline household survey, HIV status was available for 12,570 participants and 3,596 of them were HIV positive. The study was conducted in accordance with the Declaration of Helsinki. The study received institutional review board approval from the Botswana Health Research Development Committee and the U.S. Centers for Disease Control and Prevention. All participants provided written informed consent. Participants aged 16–18 years provided written assent (with parents or guardians providing written permission). The study is registered at ClinicalTrials.gov (NCT01965470).

Limiting Antigen Avidity assay and HIV recent infection algorithm

All plasma specimens from HIV-positive individuals who participated in the survey were tested using the Sedia HIV-1 Limiting Antigen (LAg)-Avidity EIA (Sedia Biosciences Corporation, Portland, OR, USA) according to manufacturer’s instructions [40]. The LAg-Avidity EIA differentiates between ‘recent’ and long-term HIV infection. A normalized optical density (ODn) of <1.5 was considered to represent recent infection [36]. ART status was verified through documentation provided by the participants or testing for presence of ARV drugs in plasma. The MAA [25] included the following steps: (1) plasma specimens from HIV-positive individuals were tested by LAg-Avidity EIA following manufacturer’s recommendations, and ODn was calculated; (2) cases with ODn ≤1.5 were checked for ART status and individuals on ART were excluded from HIV recency candidates; (3) levels of HIV-1 RNA were checked for the remaining candidates and individuals with undetectable viral load (≤400 copies/mL) were excluded from HIV recency candidates; and (4) HIV-positive individuals with ODn ≤1.5 in LAg-Avidity EIA, not taking ARV and having HIV-1 RNA >400 copies/mL were considered recently infected with HIV. Fig 1 shows the MAA applied for the cross-sectional HIV incidence estimation in BCPP. The MAA used the LAg-Avidity EIA in combination with testing for ARV and HIV-1 RNA load. The final sample of recently HIV infected individuals was identified based on HIV-positive status, ODn ≤1.5 in LAg-Avidity EIA, no use of ARVs and detectable viral load (HIV-1 RNA >400 copies/mL) at the time of testing.

[Figure omitted. See PDF.]

Fig 1. Multiassay algorithm applied for the cross-sectional HIV incidence estimation in Botswana Combination Prevention project.

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

HIV-1 RNA quantification

The HIV-1 RNA load in plasma was quantified by Abbott m2000sp/Abbott m2000rt (Wiesbaden, Germany). HIV-1 RNA >400 copies/mL was considered detectable viral load.

ARV drug testing

Plasma samples from participants who were classified as recently infected and had undetectable viral load (≤400 copies/mL) were screened for ARV drugs by high-throughput liquid chromatography coupled with Q-Exactive high-resolution mass spectrometry using data-dependent fragmentation and selected reaction monitoring at resolution of 17,500 [41]. To obtain qualitative results, each specimen was compared to positive and negative controls for each drug (abacavir, amprenavir, atazanavir, darunavir, efavirenz, emtricitabine, indinavir, lamivudine, lopinavir, maraviroc, nelfinavir, nevirapine, raltegravir, rilpivirine, ritonavir, saquinavir, stavudine, tenofovir, tipranavir, and zidovudine). The limit of identification ranged from 5 to 10 ng/ml for most drugs and is presented elsewhere [41].

Statistical analysis and estimation of HIV incidence

The annualized HIV incidence and 95% confidence intervals (CI) were estimated based on cross-sectional incidence assay-based methods that entail biomarkers of HIV disease pro-gression that can distinguish recent from long-term infections [7, 14, 42–44]. Specifically, the annual incidence was estimated as follows:

Where nR, n+, and ns represents the number of individuals who were classified as recent infection, who were HIV-positive (including both recent and non-recent infections), and who were uninfected, in the cross-sectional sample, respectively. is the mean duration of recent infection (MDRI), which is the population average time spent in the ‘recent’ state; and is thefalse-recent rate (FRR), representing the proportion of subjects who had been infected for longer than time T (set to be 730 days) but were misclassified as recent infections. MDRI was set to 130 days with a standard error of 5.98 days, corresponding to an ODn threshold of 1.5 [36, 45]. Using HIV-1 RNA load measurement and excluding 14 virologically suppressed individuals from recency candidates justified setting the false-recent rate (FRR) at zero [25, 46]. Missing LAg-Avidity EIA test results were considered missing completely at random. Confidence intervals were estimated taking into account of clustering of communities by applying a design effect to both HIV prevalence and proportion of recent infections among HIV positive individuals, implemented in the R package Inctools v. 1.0.10 [44]. The IncTools implements the HIV incidence calculations from cross-sectional surveys following the guidelines as proposed by the WHO Incidence assays technical working group [23–25, 47].

In addition to estimating HIV incidence among 16–64 years old participants, we assessed HIV incidence in a subset of younger participants in order to compare our results with other studies in Botswana. Specifically, we estimated HIV incidence in a subset of 16–49 years old participants (n = 10,164) including 2,798 HIV-positive and 7,366 HIV-negative individuals.

Results

A total of 3,596 (29%) individuals from 30 communities in Botswana were HIV positive among 12,570 adults 16–64 years old with definitive HIV status during the baseline household survey of the BCPP from 2013 to 2105 [2]. Table 1 presents basic socio-demographic and clinical characteristics of individuals participating in the baseline household survey. The median (IQR) age was 40 (33–48) years. The majority of participants were females (73%). Among HIV-positive participants, 3,581 (99.6%) were tested by the LAg-Avidity EIA.

[Figure omitted. See PDF.]

Table 1. Summary of baseline characteristics of HIV-infected participants enrolled in BCPP in 30 Botswana communities.

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

A subset of 326 participants were classified as LAg-Avidity EIA-recent HIV infections with ODn ≤1.5 (Fig 1). The documented ART status was considered as an indicator of long-term HIV infection, and 278 of 326 participants were excluded from recency candidates due to being on ART. A subset of 14 individuals who reported no prior use of ART (including 3 cases with detected ARV drugs in plasma) had undetectable HIV-1 RNA load (≤400 copies/mL), were excluded from the pool of individuals classified as recent HIV-infection [46]. Three of 14 individuals with undetectable viral load and reporting no prior ART use had ARVs in plasma that were the first-line treatment regimens most commonly prescribed in Botswana’s national ART program at the time of sampling: two cases of zidovudin/3TC/efavirenz and one case of zidovudin/3TC/nevirapine. Thus, 34 LAg-Avidity EIA-recent and ARV-naïve participants with detectable HIV-1 RNA load were classified as recent HIV infections (Fig 1). The estimate of annualized HIV incidence is 1.06% (95% CI 0.68–1.45%), assuming an FRR of zero. For a more conservative estimate, we used the adjusted FRR at 0.39% that was determined in our recent study in Botswana [48], and estimated the annualized HIV incidence at 0.64% (95% CI 0.24–1.04%). A higher proportion of recent infections were among young participants, less than 30 years of age (Table 1).

All recently infected individuals (n = 34) were younger than 49 years old. The annualized HIV incidence in the subset of 16–49 years old individuals was estimated at 1.29% (95% CI 0.82–1.77%) with FRR set to zero, and at 0.90% (95% CI 0.42–1.38%) with FRR set to 0.39%. For comparison, two alternative published estimates of HIV incidence in Botswana including UNAIDS [49] are presented in Table 2 along with results of this study.

[Figure omitted. See PDF.]

Table 2. Estimates of annualized HIV-1 incidence in Botswana.

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

Discussion

HIV incidence in a population-based sample of adults 16–64 years old residing in 30 communities across Botswana was estimated at about 1% from cross-sectional sampling that occurred in 2013–2015. Estimated HIV incidence was slightly higher (0.90–1.29%, depending on the FRR) in a subset of younger 16–49-year-old adults. Results of our study corroborate the recent UNAIDS estimates of HIV incidence in Botswana (0.93%) [49], and suggest a declining trend from previously estimated HIV incidence among 15–49 year old adults in Botswana (3.5% in 2000, 2.4% in 2007 and 1.7% in 2008 [51]). Our results support the observation that new HIV infections across sub-Saharan Africa continue to decline, although HIV incidence in Botswana remains unacceptably high [52].

The strength of the current estimate of HIV incidence includes population-based random sampling from 30 rural and peri-urban communities across the country, and application of MAA that includes LAg-Avidity EIA, ART status, and measurements of HIV-1 RNA in all HIV-positive participants. Limitations of our study include MAA that is reliant upon the ART status and the uncertainty arising from estimates of corresponding FRR and MDRI. In our previous studies, we have also found similar estimates of MDRI using an MAA with ART status and viral load [53]. Individuals on ART were excluded from recency candidates, because being on ART was interpreted as an indicator of longstanding HIV infection. This approach worked well in the era of CD4-driven initiation of ARV therapy, and was in line with the Botswana HIV treatment guidelines at the time of sampling, 2013 to 2015. In June 2016 Botswana introduced a new national policy “Treatment for All”, that is, initiating ART as soon as possible regardless of CD4+ T-cell counts [6]. The ongoing scaling up of this national policy means that ART status cannot be used as exclusion criteria for estimation of HIV recency in future studies. Novel cross-sectional assays and MAA independent of ART status are needed to address this issue.

Self-reported ART status could be considered one of the study limitations. Although ART status was verified through documentation among those who self-reported to be on ART, there is an uncertainty due to possible undisclosed ART use among those who self-reported not to be on ART. To address this limitation we performed ARV drug testing in plasma among those who had undetectable levels of HIV-1 RNA and reported no ART use. In fact, we found triple ARV drugs in 3 out of 14 cases. Using viral load threshold could minimize uncertainty of self-reported status. In fact, 3 cases with ARV drugs in plasma were excluded based on low levels of HIV-1 RNA.

While we used a Botswana-specific FRR from our previous study [48], it was estimated in a cohort sampled approximately a decade before sampling in the current study. We speculate that FRR in Botswana could be decreasing over time and could be lower than the 0.39% used for conservative estimates of HIV incidence in this study. The extent to which regional FRRs are changing over time remains unknown. The scale up of national ART programs could affect FRR estimates. Since ART guidelines have been changing, a greater proportion of individuals are initiating treatment sooner, more frequently soon after HIV diagnosis. This means the increasing ART coverage could reduce FRR. In this manuscript, the MAA algorithm included ARV drug testing leading to reduction of FRR to zero. However, uncertainly remains as to whether FRR can be eliminated [47]. The upper range of FRR then is FRR without drug tracing, as we determined in our previous study in Botswana [48, 54].

The relatively high estimated HIV incidence that we found (~1%) despite high levels of HIV testing, treatment, and viral suppression may reflect several factors. First, the impact of widespread ART on HIV incidence may take several years to be realized. In addition, the ~30% of HIV-infected individuals with detectable viremia [2] could yield such high HIV incidence in the setting of very high HIV prevalence (and these individuals could have different HIV risk behavior, compared with those with viral suppression on ART). Our findings highlight the importance of targeted interventions to reach individuals who have not yet sought HIV testing or treatment services.

Conclusion

In summary, using cross-sectional sampling and MAA based on LAg-Avidity EIA, ART status (either documented or by testing ARV drugs in plasma) and HIV-1 RNA measurements, we estimated the HIV incidence in 30 rural and peri-urban Botswana communities in 2013–2015 at about 1%. A higher proportion of recent infections were among participants less than 30 years of age. A reduction from this relatively high estimated HIV incidence may take several years to be realized despite the impact of widespread ART and other on-going interventions. Targeted interventions are required to reach individuals who have not yet sought HIV testing or treatment services.

Acknowledgments

We thank the study participants. We are grateful to the entire BCPP staff including the field study teams for making this study a success. We thank the BCPP team members for their contribution to this study: Ngozana Seonyatseng, Tumalano Sekoto Pharatlhatlhe, Vinoliah Simon, Rona Letlhogile, Atang Mbikiwa, Kutlo Manyake, Neo Mogowa, Moemedi Tshwenyana, Kagiso Watema, Chebukani Nkobodo, Phelimon Pong Sebogodi, Thabani Ncube, Dineo Mongwato, Coulsen Kgathi, Thuso Mokane, Kutlwano Mukokomani, Mompati Mmalane, Ria Madison, Chloe Auletta-Young, Botswana Harvard HIV Reference Laboratory Staff, Data Management Centre Staff, and field laboratory assistants. We thank Lendsey Melton for excellent editorial assistance. We thank Michelle Roland for helpful critique and discussion. We would like to thank Reshma Kassanjee, Alex Welte and Eduard Grebe for useful discussions in the application of cross-sectional methods and use of IncTools.

Citation: Moyo S, Gaseitsiwe S, Mohammed T, Pretorius Holme M, Wang R, Kotokwe KP, et al. (2018) Cross-sectional estimates revealed high HIV incidence in Botswana rural communities in the era of successful ART scale-up in 2013-2015. PLoS ONE 13(10): e0204840. https://doi.org/10.1371/journal.pone.0204840

References

1. UNAIDS. 90-90-90 On the right track towards the global target. URL http://www.cfenet.ubc.ca/sites/default/files/uploads/IAS2016/90_90_90_Progress_ReportFINAL.pdf. 2016.

2. Gaolathe T, Wirth KE, Holme MP, Makhema J, Moyo S, Chakalisa U, et al. Botswana's progress toward achieving the 2020 UNAIDS 90-90-90 antiretroviral therapy and virological suppression goals: a population-based survey. Lancet HIV. 2016;3(5):e221–30. pmid:27126489.

3. UNAIDS. Global AIDS Update 2016. http://www.unaids.org/sites/default/files/media_asset/global-AIDS-update-2016_en.pdf. 2016.

4. UNAIDS. UNAIDS Data 2017 http://www.unaids.org/sites/default/files/media_asset/20170720_Data_book_2017_en.pdf. 2017.

5. UNAIDS. GLOBAL AIDS RESPONSE PROGRESS REPORTING 2015 http://www.unaids.org/sites/default/files/media_asset/JC2702_GARPR2015guidelines_en.pdf. 2015.

6. World Health Organisation. Botswana Launches Treat All Strategy. 2016. Available from: http://www.afro.who.int/news/botswana-launches-treat-all-strategy.

7. Brookmeyer R, Quinn TC. Estimation of current human immunodeficiency virus incidence rates from a cross-sectional survey using early diagnostic tests. Am J Epidemiol. 1995;141(2):166–72. pmid:7817972

8. Brookmeyer R, Konikoff J, Laeyendecker O, Eshleman SH. Estimation of HIV incidence using multiple biomarkers. Am J Epidemiol. 2013;177(3):264–72. pmid:23302151; PubMed Central PMCID: PMC3626051.

9. Brookmeyer R, Laeyendecker O, Donnell D, Eshleman SH. Cross-sectional HIV incidence estimation in HIV prevention research. J Acquir Immune Defic Syndr. 2013;63 Suppl 2:S233–9. pmid:23764641; PubMed Central PMCID: PMC3737594.

10. Hall HI, Song R, Rhodes P, Prejean J, An Q, Lee LM, et al. Estimation of HIV incidence in the United States. JAMA. 2008;300(5):520–9. pmid:18677024; PubMed Central PMCID: PMC2919237.

11. Janssen RS, Satten GA, Stramer SL, Rawal BD, O'Brien TR, Weiblen BJ, et al. New testing strategy to detect early HIV-1 infection for use in incidence estimates and for clinical and prevention purposes. JAMA. 1998;280(1):42–8. pmid:9660362.

12. Hargrove JW, Humphrey JH, Mutasa K, Parekh BS, McDougal JS, Ntozini R, et al. Improved HIV-1 incidence estimates using the BED capture enzyme immunoassay. AIDS. 2008;22(4):511–8. Epub 2008/02/28. 00002030-200802190-00009 [pii]. pmid:18301064.

13. McDougal JS, Parekh BS, Peterson ML, Branson BM, Dobbs T, Ackers M, et al. Comparison of HIV type 1 incidence observed during longitudinal follow-up with incidence estimated by cross-sectional analysis using the BED capture enzyme immunoassay. AIDS Res Hum Retroviruses. 2006;22(10):945–52. Epub 2006/10/28. pmid:17067263.

14. McDougal JS, Pilcher CD, Parekh BS, Gershy-Damet G, Branson BM, Marsh K, et al. Surveillance for HIV-1 incidence using tests for recent infection in resource-constrained countries. AIDS. 2005;19 Suppl 2:S25–30. pmid:15930838.

15. Parekh BS, Hanson DL, Hargrove J, Branson B, Green T, Dobbs T, et al. Determination of mean recency period for estimation of HIV type 1 Incidence with the BED-capture EIA in persons infected with diverse subtypes. AIDS Res Hum Retroviruses. 2011;27(3):265–73. pmid:20954834.

16. Parekh BS, Kennedy MS, Dobbs T, Pau CP, Byers R, Green T, et al. Quantitative detection of increasing HIV type 1 antibodies after seroconversion: a simple assay for detecting recent HIV infection and estimating incidence. AIDS Res Hum Retroviruses. 2002;18(4):295–307. pmid:11860677.

17. Parekh BS, McDougal JS. Application of laboratory methods for estimation of HIV-1 incidence. The Indian journal of medical research. 2005;121(4):510–8. Epub 2005/04/09. pmid:15817960.

18. Parekh BS, Pau CP, Kennedy MS, Dobbs TL, McDougal JS. Assessment of Antibody Assays for Identifying and Distinguishing Recent from Long-Term HIV Type 1 Infection. AIDS Res Hum Retroviruses. 2001;17(2):137–46. pmid:11177393.

19. Guy R, Gold J, Calleja JM, Kim AA, Parekh B, Busch M, et al. Accuracy of serological assays for detection of recent infection with HIV and estimation of population incidence: a systematic review. Lancet Infect Dis. 2009;9(12):747–59. Epub 2009/11/21. doi: S1473-3099(09)70300-7 [pii] pmid:19926035.

20. Hargrove J, Eastwood H, Mahiane G, van Schalkwyk C. How should we best estimate the mean recency duration for the BED method? PLoS One. 2012;7(11):e49661. pmid:23166743; PubMed Central PMCID: PMC3500313.

21. Busch MP, Pilcher CD, Mastro TD, Kaldor J, Vercauteren G, Rodriguez W, et al. Beyond detuning: 10 years of progress and new challenges in the development and application of assays for HIV incidence estimation. AIDS. 2010;24(18):2763–71. pmid:20975514.

22. Mastro TD, Kim AA, Hallett T, Rehle T, Welte A, Laeyendecker O, et al. Estimating HIV Incidence in Populations Using Tests for Recent Infection: Issues, Challenges and the Way Forward. J HIV AIDS Surveill Epidemiol. 2010;2(1):1–14. pmid:21743821; PubMed Central PMCID: PMC3130510.

23. Kassanjee R, McWalter TA, Barnighausen T, Welte A. A new general biomarker-based incidence estimator. Epidemiology. 2012;23(5):721–8. pmid:22627902; PubMed Central PMCID: PMC3500970.

24. Kassanjee R, Pilcher CD, Keating SM, Facente SN, McKinney E, Price MA, et al. Independent assessment of candidate HIV incidence assays on specimens in the CEPHIA repository. AIDS. 2014;28(16):2439–49. pmid:25144218; PubMed Central PMCID: PMC4210690.

25. Kassanjee R, Pilcher CD, Busch MP, Murphy G, Facente SN, Keating SM, et al. Viral load criteria and threshold optimization to improve HIV incidence assay characteristics. AIDS. 2016;30(15):2361–71. pmid:27454561; PubMed Central PMCID: PMCPMC5024345.

26. Hallett TB. Estimating the HIV incidence rate: recent and future developments. Curr Opin HIV AIDS. 2011;6(2):102–7. Epub 2011/04/21. 01222929-201103000-00004 [pii]. pmid:21505383.

27. Moyo S, Wilkinson E, Novitsky V, Vandormael A, Gaseitsiwe S, Essex M, et al. Identifying Recent HIV Infections: From Serological Assays to Genomics. Viruses. 2015;7(10):5508–24. pmid:26512688; PubMed Central PMCID: PMC4632395.

28. Rosenberg NE, Pilcher CD, Busch MP, Cohen MS. How can we better identify early HIV infections? Curr Opin HIV AIDS. 2014. Epub 2014/11/13. pmid:25389806.

29. Mastro TD, Kim AA, Hallett T, Rehle T, Welte A, Laeyendecker O, et al. Estimating HIV incidence in populations using tests for recent infection: issues, challenges and the way forward. Journal of HIV AIDS surveillance & epidemiology. 2010;2(1):1.

30. Eshleman SH, Hughes JP, Laeyendecker O, Wang J, Brookmeyer R, Johnson-Lewis L, et al. Use of a multifaceted approach to analyze HIV incidence in a cohort study of women in the United States: HIV Prevention Trials Network 064 Study. J Infect Dis. 2013;207(2):223–31. pmid:23129758; PubMed Central PMCID: PMC3532822.

31. Konikoff J, Brookmeyer R, Longosz AF, Cousins MM, Celum C, Buchbinder SP, et al. Performance of a limiting-antigen avidity enzyme immunoassay for cross-sectional estimation of HIV incidence in the United States. PLoS One. 2013;8(12):e82772. pmid:24386116; PubMed Central PMCID: PMC3873916.

32. Laeyendecker O, Brookmeyer R, Cousins MM, Mullis CE, Konikoff J, Donnell D, et al. HIV incidence determination in the United States: a multiassay approach. J Infect Dis. 2013;207(2):232–9. pmid:23129760; PubMed Central PMCID: PMC3532826.

33. Laeyendecker O, Brookmeyer R, Mullis C, Donnell D, Lingappa J, Celum C, et al. Specificity of four laboratory approaches for cross-sectional HIV incidence determination: Analysis of samples from adults with known non-recent HIV infection from five African countries. AIDS Res Hum Retroviruses. 2012. Epub 2012/01/31. pmid:22283149.

34. Laeyendecker O, Kulich M, Donnell D, Komarek A, Omelka M, Mullis CE, et al. Development of Methods for Cross-Sectional HIV Incidence Estimation in a Large, Community Randomized Trial. PLoS One. 2013;8(11):e78818. pmid:24236054; PubMed Central PMCID: PMC3827276.

35. Laeyendecker O, Piwowar-Manning E, Fiamma A, Kulich M, Donnell D, Bassuk D, et al. Estimation of HIV incidence in a large, community-based, randomized clinical trial: NIMH project accept (HIV Prevention Trials Network 043). PLoS One. 2013;8(7):e68349. pmid:23874597; PubMed Central PMCID: PMC3708944.

36. Duong YT, Kassanjee R, Welte A, Morgan M, De A, Dobbs T, et al. Recalibration of the Limiting Antigen Avidity EIA to Determine Mean Duration of Recent Infection in Divergent HIV-1 Subtypes. PLoS One. 2015;10(2):e0114947. pmid:25710171; PubMed Central PMCID: PMC4339840.

37. Wei X, Liu X, Dobbs T, Kuehl D, Nkengasong JN, Hu DJ, et al. Development of two avidity-based assays to detect recent HIV type 1 seroconversion using a multisubtype gp41 recombinant protein. AIDS Res Hum Retroviruses. 2010;26(1):61–71. Epub 2010/01/13. pmid:20063992.

38. Novitsky V, Prague M, Moyo S, Gaolathe T, Mmalane M, Kadima Yankinda E, et al. High HIV-1 RNA among Newly Diagnosed People in Botswana. AIDS Res Hum Retroviruses. 2017. pmid:29214845.

39. Novitsky V, Gaolathe T, Mmalane M, Moyo S, Chakalisa U, Kadima Yankinda E, et al. Lack of Virological Suppression among Young HIV-Positive Adults in Botswana. JAIDS. 2018 [In Press].

40. Duong YT, Qiu M, De AK, Jackson K, Dobbs T, Kim AA, et al. Detection of Recent HIV-1 Infection Using a New Limiting-Antigen Avidity Assay: Potential for HIV-1 Incidence Estimates and Avidity Maturation Studies. PLoS One. 2012;7(3):e33328. Epub 2012/04/06. PONE-D-11-24174 [pii]. pmid:22479384.

41. Marzinke MA, Breaud A, Parsons TL, Cohen MS, Piwowar-Manning E, Eshleman SH, et al. The development and validation of a method using high-resolution mass spectrometry (HRMS) for the qualitative detection of antiretroviral agents in human blood. Clin Chim Acta. 2014;433:157–68. Epub 2014/03/26. pmid:24661980; PubMed Central PMCID: PMCPMC4039613.

42. Brookmeyer R, Quinn TC. Estimation of current human immunodeficiency virus incidence rates from a cross-sectional survey using early diagnostic tests. Am J Epidemiol. 1995;141(2):166–72. pmid:7817972

43. Wang R, Lagakos SW. On the use of adjusted cross-sectional estimators of HIV incidence. J Acquir Immune Defic Syndr. 2009;52(5):538–47. Epub 2009/11/03. pmid:19881357.

44. Welte A, Grebe E, McIntosh A, Bäumler P, Kassanjee R, Brand H, et al. inctools: Incidence Estimation Tools. http://www.incidence-estimation.org/page/inctools. 2017.

45. Moyo S, Mohammed T, Kotokwe K, Boleo C, Mupfumi L, Gaseitsiwe S, et al. Cross-sectional estimates of HIV incidence remain high in rural communities in Botswana in the era of successful scale-up of ART. Abstract TUPDC0103. IAS2016.

46. Rehle T, Johnson L, Hallett T, Mahy M, Kim A, Odido H, et al. A Comparison of South African National HIV Incidence Estimates: A Critical Appraisal of Different Methods. PLoS One. 2015;10(7):e0133255. pmid:26230949; PubMed Central PMCID: PMC4521952.

47. World Health Organization. Technical update on HIV incidence assays for surveillance and monitoring purposes. Geneva: World Health Organization Copyright (c) World Health Organization; 2015. Available from: http://www.who.int/hiv/pub/surveillance/technical-update-incidence/en/.

48. Moyo S, Kotokwe KP, Mohammed T, Boleo C, Mupfumi L, Chishala S, et al. Low False Recent Rate of Limiting Antigen-Avidity Assay Combined with HIV-1 RNA Data in Botswana. AIDS Res Hum Retroviruses. 2017;33(1):17–8. pmid:27481530.

49. UNAIDS. Country factsheets. BOTSWANA 2016. http://www.unaids.org/en/regionscountries/countries/botswana 2017.

50. Central Statistics Office. Botswana AIDS Impact Survey IV—2013. Gaborone: Republic of Botswana; 2014. Available from: http://www.statsbots.org.bw/sites/default/files/publications/Botswana AIDS Impact Survey Summary Apr 2014.pdf.

51. Stover J, Fidzani B, Molomo BC, Moeti T, Musuka G. Estimated HIV trends and program effects in Botswana. PLoS One. 2008;3(11):e3729. pmid:19008957; PubMed Central PMCID: PMCPMC2579326.

52. Kharsany AB, Karim QA. HIV Infection and AIDS in Sub-Saharan Africa: Current Status, Challenges and Opportunities. Open AIDS J. 2016;10:34–48. Epub 2016/06/28. pmid:27347270; PubMed Central PMCID: PMCPMC4893541.

53. Moyo S, Vandormael A, Wilkinson E, Engelbrecht S, Gaseitsiwe S, Kotokwe KP, et al. Analysis of Viral Diversity in Relation to the Recency of HIV-1C Infection in Botswana. PLoS One. 2016;11(8):e0160649. pmid:27552218.

54. World Health Organization. Foundation for Innovative New Diagnostics, WHO Working Group on HIV Incidence Assays. Meeting report, Boston, MA, USA, 20–26 February 2016. Geneva: World Health Organization; 2017. Available from: http://apps.who.int/iris/bitstream/handle/10665/254868/WHO-HIV-2017.04-eng.PDF?sequence=1.