A general HIV incidence inference scheme based on likelihood of individual level data and a population renewal equation

Evidence of HIV incidence reduction in young women, but not in adolescent girls, in KwaZulu-Natal, South Africa

Objectives

We estimated changes in the HIV incidence from 2013-2018 in Eshowe/Mbongolwane, KwaZulu-Natal, South Africa where Médecins Sans Frontières is engaged in providing HIV testing and care since 2011.

Methods

Using data from two cross-sectional household-based surveys conducted in 2013 and 2018, with consenting participants aged 15-59 years, we applied the incidence estimation frameworks of Mahiane et al and Kassanjee et al.

Results

In total, 5599 (62.4% women) and 3276 (65.9% women) individuals were included in 2013 and 2018, respectively. We found a mean incidence in women aged 20-29 years of 2.71 cases per 100 person-years (95% confidence interval [CI]: 1.23;4.19) in 2013 and 0.4 cases per 100 person-years (95% CI: 0.0;1.5) in 2018. The incidence in men aged 20-29 years was 1.91 cases per 100 person-years (95% CI: 0.87; 2.93) in 2013 and 0.53 cases per 100 person-years (95% CI: 0.0; 1.4) in 2018. The incidence decline among women aged 15-19 was -0.34 cases per 100 person-years (95% CI: -1.31;0.64).

Conclusions

The lack of evidence of incidence decline among adolescent girls is noteworthy and disconcerting. Our findings suggest that large-scale surveys should seriously consider focusing their resources on the core group of women aged 15-19 years.

HIV incidence estimation among female sex workers in South Africa: a multiple methods analysis of cross-sectional survey data

BACKGROUND

Although numerous studies have investigated HIV risk factors and shown high HIV prevalence among female sex workers in South Africa, no national HIV incidence estimate exists for this potentially important group for HIV transmission. We aimed to estimate HIV incidence among female sex workers in South Africa who could be accessed through sex worker programmes, and to refine and describe the methods that enabled analysis.

METHODS

This study was embedded in a cross-sectional national survey of female sex workers who were linked to sex worker programmes. We aimed to enrol 3000 female sex workers aged at least 18 years who had sold or transacted in sex in the preceding 6 months in 12 randomly selected districts of the 22 districts with sex worker programmes, ensuring coverage of all provinces of South Africa. Women who self-reported as current victims of human trafficking were excluded from enrolment. We used a multistep process to sample districts and then hotspots, and a chain referral method to recruit participants. We collected cross-sectional data for self-reported HIV status, demographic characteristics, and exposure to violence. Two rapid tests were used to ascertain diagnostic markers, a viral load assay was used to ascertain clinical markers, and the Maxim Limiting Antigen Avidity EIA was used to ascertain infection-staging HIV markers. Given the challenges of estimating HIV incidence, especially cross-sectionally, multiple methods of estimation were adapted to our setting, leveraging the age structure of HIV prevalence, recency-of -infection biomarker results (ie, where recent infection is classified as ≤1·5 normalised optical density [ODn] on the avidity assay and viral load of ≥1000 copies per mL), and reported testing histories.

FINDINGS

Of 3005 female sex workers who were enrolled and interviewed between Feb 4 and June 26, 2019, 2999 who had HIV test results were included in this analysis. The median age of participants was 32 years (IQR 27-38). 1714 (57·2%) of 2999 participants self-reported as being HIV positive, and 1447 (48·3%) of 2993 participants reported client sexual violence in the past year. The measured HIV prevalence was 62·1% (95% CI 60·3-65·7) and peaked at approximately age 40 years. Using recency-of-infection biomarker results, we obtained a base case estimate of HIV incidence of 4·60 cases per 100 person-years (95% CI 1·53-8·45) for the population. Estimates were generally consistent by method, and outlying incidence estimates calculated by self-reported testing histories were considered unreliable. Various sensitivity analyses produced estimates up to 11 cases per 100 person-years, and we did not detect differences by age and region.

INTERPRETATION

We found that female sex workers have extraordinarily high HIV incidence of approximately 5 cases per 100 person-years, emphasising the need to sustain and strengthen efforts to mitigate risk and provide adequate care. The notable role that sex work has in HIV transmission demands substantial investment in ongoing epidemiological monitoring.

FUNDING

South African Medical Research Council, South African National Treasury, Global Fund, South African Department of Science and Innovation, Wellcome Trust.

Estimating HIV incidence in the Akwa Ibom AIDS indicator survey (AKAIS), Nigeria using the limiting antigen avidity recency assay

INTRODUCTION

HIV incidence estimates are important to characterize the status of an epidemic, identify locations and populations at high risk and to guide and evaluate HIV prevention interventions. We used the limiting antigen avidity assay (LAg) as part of a recent infection testing algorithm to estimate HIV incidence in the Akwa Ibom AIDS Indicator Survey (AKAIS), Nigeria.

METHODS

In 2017, AKAIS, a cross-sectional population-based study was conducted at the household (HH) level in 31 local government areas (LGAs) of Akwa Ibom state. Of the 8963 participants aged ≥15 years who were administered questionnaires for demographic and behavioural data, 8306 consented to HIV rapid testing. Whole-blood specimens were collected from 394 preliminary HIV-seropositive individuals for CD4+ cell count determination and plasma storage. Samples were shipped to a central quality laboratory for HIV confirmatory testing and viral load determination. A total of 370 HIV-positive specimens were tested for the recent HIV infection using the LAg assay.

RESULTS

Of the 8306 consenting adults, the HIV prevalence was 4.8%. Of the 370 HIV-positive samples tested for HIV recency, the median age was 35 years, 48.8% had CD4+ cell count >500/mm3 and 81.3% was not virally suppressed. Viral suppression was greater among females (21%) than for males (13%). A total of 11 specimens were classified as recent based on the LAg assay and HIV viral load ≥1000 copies/mL. The weighted, adjusted HIV-1 incidence was 0.41/100 person-years (95% CI 0.16 to 0.66); translating to 13,000 new cases of HIV infections annually in Akwa Ibom, a state with a population of 5.5 million. The HIV incidence rate was similar in females and males (0.41% and 0.42% respectively). The incidence rate was the highest among participants aged 15 to 49 years (0.44%, 95% CI 0.15 to 0.74) translating to 11,000 new infections annually, about 85% of all new infections in the state.

CONCLUSIONS

The finding of the high HIV incidence among the 15 to 49-year age group calls for renewed and innovative efforts to prevent HIV infection among young adults in Akwa Ibom state.

Performance comparison of the Maxim and Sedia Limiting Antigen Avidity assays for HIV incidence surveillance

BACKGROUND

Two manufacturers, Maxim Biomedical and Sedia Biosciences Corporation, supply CDC-approved versions of the HIV-1 Limiting Antigen Avidity EIA (LAg) for detecting 'recent' HIV infection in cross-sectional incidence estimation. This study assesses and compares the performance of the two assays for incidence surveillance.

METHODS

We ran both assays on a panel of 2,500 well-characterized HIV-1-infected specimens. We analysed concordance of assay results, assessed reproducibility using repeat testing and estimated mean durations of recent infection (MDRIs) and false-recent rates (FRRs) for a range of normalized optical density (ODn) thresholds, alone and in combination with viral load thresholds. We defined three hypothetical surveillance scenarios, similar to the Kenyan and South African epidemics, and a concentrated epidemic. These scenarios allowed us to evaluate the precision of incidence estimates obtained by means of various recent infection testing algorithms (RITAs) based on each of the two assays.

RESULTS

The Maxim assay produced lower ODn values than the Sedia assay on average, largely as a result of higher calibrator readings (mean OD of 0.749 vs. 0.643), with correlation of normalized readings lower (R2 = 0.908 vs. R2 = 0.938). Reproducibility on blinded control specimens was slightly better for Maxim. The MDRI of a Maxim-based algorithm at the 'standard' threshold (ODn ≤1.5 & VL >1,000) was 201 days (95% CI: 180,223) and for Sedia 171 (152,191). The difference Differences in MDRI were estimated at 32.7 (22.9,42.8) and 30.9 days (21.7,40.7) for the two algorithms, respectively. Commensurately, the Maxim algorithm had a higher FRR in treatment-naive subjects (1.7% vs. 1.1%). The two assays produced similar precision of incidence estimates in the three surveillance scenarios.

CONCLUSIONS

Differences between the assays can be primarily attributed to the calibrators supplied by the manufacturers. Performance for surveillance was extremely similar, although different thresholds were optimal (i.e. produced the lowest variance of incidence estimates) and at any given ODn threshold, different estimates of MDRI and FRR were obtained. The two assays cannot be treated as interchangeable: assay and algorithm-specific performance characteristic estimates must be used for survey planning and incidence estimation.

Second Order Segmented Polynomials for Syphilis and Gonorrhea Prevalence and Incidence Trends Estimation: Application to Spectrum's Guinea-Bissau and South Africa Data

This paper presents two approaches to smoothing time trends in prevalence and estimating the underlying incidence of remissible infections. In the first approach, we use second order segmented polynomials to smooth a curve in a bounded domain. In the second, incidence is modeled instead and the prevalence is reconstructed using the recovery rate which is assumed to be known. In both approaches, the number of knots and their positions are estimated, resulting in non-linear regressions. Akaike Information Criterion is used for model selection. The method is illustrated with Syphilis and Gonorrhea prevalence smoothing and incidence trend estimation in Guinea-Bissau and South Africa, respectively.

Illness-death model: statistical perspective and differential equations

The aim of this work is to relate the theory of stochastic processes with the differential equations associated with multistate (compartment) models. We show that the Kolmogorov Forward Differential Equations can be used to derive a relation between the prevalence and the transition rates in the illness-death model. Then, we prove mathematical well-definedness and epidemiological meaningfulness of the prevalence of the disease. As an application, we derive the incidence of diabetes from a series of cross-sections.

Population-level HIV incidence estimates using a combination of synthetic cohort and recency biomarker approaches in KwaZulu-Natal, South Africa

Introduction

There is a notable absence of consensus on how to generate estimates of population-level incidence. Incidence is a considerably more sensitive indicator of epidemiological trends than prevalence, but is harder to estimate. We used a novel hybrid method to estimate HIV incidence by age and sex in a rural district of KwaZulu-Natal, South Africa.

Methods

Our novel method uses an ‘optimal weighting’ of estimates based on an implementation of a particular ‘synthetic cohort’ approach (interpreting the age/time structure of prevalence, in conjunction with estimates of excess mortality) and biomarkers of ‘recent infection’ (combining Lag-Avidity, Bio-Rad Avidity and viral load results to define recent infection, and adapting the method for age-specific incidence estimation). Data were obtained from a population-based cross-sectional HIV survey conducted in Mbongolwane and Eshowe health service areas in 2013.

Results

Using the combined method, we find that age-specific HIV incidence in females rose rapidly during adolescence, from 1.33 cases/100 person-years (95% CI: 0.98,1.67) at age 15 to a peak of 5.01/100PY (4.14,5.87) at age 23. In males, incidence was lower, 0.34/100PY (0.00-0.74) at age 15, and rose later, peaking at 3.86/100PY (2.52-5.20) at age 30. Susceptible population-weighted average incidence in females aged 15-29 was estimated at 3.84/100PY (3.36-4.40), in males aged 15-29 at 1.28/100PY (0.68-1.50) and in all individuals aged 15-29 at 2.55/100PY (2.09-2.76). Using the conventional recency biomarker approach, we estimated HIV incidence among females aged 15-29 at 2.99/100PY (1.79-4.36), among males aged 15-29 at 0.87/100PY (0.22-1.60) and among all individuals aged 15-59 at 1.66/100PY (1.13-2.27).

Discussion

HIV incidence was very high in women aged 15-30, peaking in the early 20s. Men had lower incidence, which peaked at age 30. The estimates obtained from the hybrid method are more informative than those produced by conventional analysis of biomarker data, and represents a more optimal use of available data than either the age-continuous biomarker or synthetic cohort methods alone. The method is mainly useful at younger ages, where excess mortality is low and uncertainty in the synthetic cohort estimates is reasonably small.

Conclusion

Application of this method to large-scale population-based HIV prevalence surveys is likely to result in improved incidence surveillance over methods currently in wide use. Reasonably accurate and precise age-specific estimates of incidence are important to target better prevention, diagnosis and care strategies.

Infection Staging and Incidence Surveillance Applications of High Dynamic Range Diagnostic Immuno-Assay Platforms

BACKGROUND

Custom HIV staging assays, including the Sedia HIV-1 Limiting Antigen (LAg) Avidity EIA and avidity modifications of the Ortho VITROS anti-HIV-1+2 and Abbott ARCHITECT HIV Ag/Ab Combo assays, are used to identify "recent" infections in clinical settings and for cross-sectional HIV incidence estimation. However, the high dynamic range of chemiluminescent platforms allows differentiating recent and long-standing infection on signal intensity, and this raises the prospect of using unmodified diagnostic assays for infection timing and surveillance applications.

METHODS

We tested a panel of 2500 well-characterized specimens with estimable duration of HIV infection with the 3 assays and the unmodified ARCHITECT. Regression models were used to estimate mean durations of recent infection (MDRIs), context-specific false-recent rates (FRRs) and correlation between diagnostic signal intensity and LAg measurements. Hypothetical epidemiological scenarios were constructed to evaluate utility in surveillance applications.

RESULTS

Over a range of MDRIs (reflecting recency discrimination thresholds), a diluted ARCHITECT-based RITA produced lower FRRs than the VITROS platform (FRR ≈ 0.5% and 1.5%, respectively at MDRI ≈ 200 days), and the unmodified diagnostic ARCHITECT produces incidence estimates with comparable precision to LAg (relative SE ≈ 17.5% and 15%, respectively at MDRI ≈ 200 days). ARCHITECT S/CO measurements were highly correlated with LAg optical density measurements (r = 0.80), and values below 200 are strongly predictive of LAg recency and duration of infection less than 1 year.

CONCLUSIONS

Low quantitative measurements from the unmodified ARCHITECT obviate the need for additional recency testing, and its use is feasible in clinical staging and incidence surveillance applications.

Estimating HIV Incidence Using a Cross-Sectional Survey: Comparison of Three Approaches in a Hyperendemic Setting, Ndhiwa Subcounty, Kenya, 2012

OBJECTIVES

Estimating HIV incidence is critical for identifying groups at risk for HIV infection, planning and targeting interventions, and evaluating these interventions over time. The use of reliable estimation methods for HIV incidence is thus of high importance. The aim of this study was to compare methods for estimating HIV incidence in a population-based cross-sectional survey.

DESIGN/METHODS

The incidence estimation methods evaluated included assay-derived methods, a testing history-derived method, and a probability-based method applied to data from the Ndhiwa HIV Impact in Population Survey (NHIPS). Incidence rates by sex and age and cumulative incidence as a function of age were presented.

RESULTS

HIV incidence ranged from 1.38 [95% confidence interval (CI) 0.67-2.09] to 3.30 [95% CI 2.78-3.82] per 100 person-years overall; 0.59 [95% CI 0.00-1.34] to 2.89 [95% CI 0.86-6.45] in men; and 1.62 [95% CI 0.16-6.04] to 4.03 [95% CI 3.30-4.77] per 100 person-years in women. Women had higher incidence rates than men for all methods. Incidence rates were highest among women aged 15-24 and 25-34 years and highest among men aged 25-34 years.

CONCLUSION

Comparison of different methods showed variations in incidence estimates, but they were in agreement to identify most-at-risk groups. The use and comparison of several distinct approaches for estimating incidence are important to provide the best-supported estimate of HIV incidence in the population.

Segmented polynomials for incidence rate estimation from prevalence data

The study considers the problem of estimating incidence of a non remissible infection (or disease) with possibly differential mortality using data from a(several) cross-sectional prevalence survey(s). Fitting segmented polynomial models is proposed to estimate the incidence as a function of age, using the maximum likelihood method. The approach allows automatic search for optimal position of knots, and model selection is performed using the Akaike information criterion. The method is applied to simulated data and to estimate HIV incidence among men in Zimbabwe using data from both the NIMH Project Accept (HPTN 043) and Zimbabwe Demographic Health Surveys (2005-2006). Copyright © 2016 John Wiley & Sons, Ltd.

Bayesian inference for age-structured population model of infectious disease with application to varicella in Poland

The dynamics of the infectious disease transmission are often best understood by taking into account the structure of population with respect to specific features, for example age or immunity level. The practical utility of such models depends on the appropriate calibration with the observed data. Here, we discuss the Bayesian approach to data assimilation in the case of a two-state age-structured model. Such models are frequently used to explore the disease dynamics (i.e. force of infection) based on prevalence data collected at several time points. We demonstrate that, in the case when the explicit solution to the model equation is known, accounting for the data collection process in the Bayesian framework allows us to obtain an unbiased posterior distribution for the parameters determining the force of infection. We further show analytically and through numerical tests that the posterior distribution of these parameters is stable with respect to a cohort approximation (Escalator Boxcar Train) of the solution. Finally, we apply the technique to calibrate the model based on observed sero-prevalence of varicella in Poland.

The HIV epidemic in Southern Africa - Is an AIDS-free generation possible?

Southern Africa, home to about 20 % of the global burden of infection continues to experience high rates of new HIV infection despite substantial programmatic scale-up of treatment and prevention interventions. While several countries in the region have had substantial reductions in HIV infection, almost half a million new infections occurred in this region in 2012. Sexual transmission remains the dominant mode of transmission. A recent national household survey in Swaziland revealed an HIV prevalence of 14.3 % among 18-19 year old girls, compared to 0.8 % among their male peers. Expanded ART programmes in Southern Africa have resulted in dramatically decreased HIV incidence and HIV mortality rates. In South Africa alone, it is estimated that more than 2.1 million of the 6.1 million HIV-positive people were receiving ART by the end of 2012, and that this resulted in more than 2.7 million life-years saved, and hundreds of thousands of HIV infections averted. Biological, behavioural and structural factors all contribute to the ongoing high rates of new HIV infection; however, as the epidemic matures and mortality is reduced from increased ART coverage, epidemiological trends become hard to quantify. What is clear is that a key driver of the Southern African epidemic is the high incidence rate of infection in young women, a vulnerable population with limited prevention options. Moreover, whilst ongoing trials of combination prevention, microbicides and behavioural economics hold promise for further epidemic control, an AIDS-free generation will not be realised unless incident infections in key populations are reduced.

Effectiveness of antiviral prophylaxis coupled with contact tracing in reducing the transmission of the influenza A (H1N1-2009): a systematic review

BACKGROUND

During the very early stage of the 2009 pandemic, mass chemoprophylaxis was implemented as part of containment measure. The purposes of the present study were to systematically review the retrospective studies that investigated the effectiveness of antiviral prophylaxis during the 2009 pandemic, and to explicitly estimate the effectiveness by employing a mathematical model.

METHODS

A systematic review identified 17 articles that clearly defined the cases and identified exposed individuals based on contact tracing. Analysing a specific school-driven outbreak, we estimated the effectiveness of antiviral prophylaxis using a renewal equation model. Other parameters, including the reproduction number and the effectiveness of antiviral treatment and school closure, were jointly estimated.

RESULTS

Based on the systematic review, median secondary infection risks (SIRs) among exposed individuals with and without prophylaxis were estimated at 2.1% (quartile: 0, 12.2) and 16.6% (quartile: 8.4, 32.4), respectively. A very high heterogeneity in the SIR was identified with an estimated I2 statistic at 71.8%. From the outbreak data in Madagascar, the effectiveness of mass chemoprophylaxis in reducing secondary transmissions was estimated to range from 92.8% to 95.4% according to different model assumptions and likelihood functions, not varying substantially as compared to other parameters.

CONCLUSIONS

Only based on the meta-analysis of retrospective studies with different study designs and exposure settings, it was not feasible to estimate the effectiveness of antiviral prophylaxis in reducing transmission. However, modelling analysis of a single outbreak successfully yielded an estimate of the effectiveness that appeared to be robust to model assumptions. Future studies should fill the data gap that has existed in observational studies and allow mathematical models to be used for the analysis of meta-data.