Obtaining an accurate estimation of HIV incidence is critical for monitoring trends in the epidemic, assessing public health interventions, and developing prevention strategies. Reliable identification of individuals within the early phase of HIV infection opens an opportunity for behavioral and biomedical interventions directed toward reducing viral transmission, initiating early ART, and providing important information for designing new therapeutic and vaccine strategies through analysis of recently transmitted viruses and immune responses in primary HIV infection.
Longitudinal cohorts represent a classical way of studying HIV incidence by prospective follow up of HIV-negative individuals with periodic testing for HIV. However, methods have been developed enabling the analysis of HIV incidence through cross-sectional sampling, which is significantly less expensive, faster, and logistically more feasible.1-11 The widespread use of cross-sectional methods for analyzing HIV incidence has resulted in a substantial increase of data on trends in the HIV/AIDS epidemic in different parts of the world, including studies in HIV-1 subtype C settings.12
However, as noted in a recent report of the Institute of Medicine,13 an important concern with the use of cross-sectional testing is the reliability and generalizability of results. Long-term infections present a challenge for HIV incidence testing because even a small number of such cases might inflate the incidence estimate,12 particularly if prevalence is high relative to incidence.14 Suppression of viral replication with or without ART is associated with lower titers of anti-HIV antibodies, which may lead to misclassification of long-lasting infections as recent cases. Low levels of viral load have been found to account for misclassification15 and have been found in samples tested discordantly by alternative laboratory methods.16 Although it is known that patients treated with ART have decreased anti-HIV titers,15,17 the frequency and the time line for such decrease are still unclear.
Different methods for estimating HIV incidence-such as detuned Vironostika, BED, and avidity assay-can perform similarly well18-20 and might complement each other (eg, detuned Vironostika and avidity assay) in refining false-positive samples.16 However, a recent report from Cote d'Ivoire demonstrated a wide range of estimated HIV incidence (1.2%-11.2%) by different laboratory methods.21 It is likely that any particular diagnostic method might require adjustment to a local epidemic before it can be reliably used for accurate estimation of HIV incidence.
The rationale for this study was based on the assumption that with frequent longitudinal measurements, even a small number of cases within the early phase of HIV infection can provide a reliable estimate of the duration of the recency period and can serve as a prototype for seroconversion panels for calibration of the recency period, which in turn can provide a more accurate estimate of HIV incidence in the local epidemic. In this study, we characterized standardized optical density (SOD) evolution during the first year of infection in a cohort of subjects identified with primary HIV-1 subtype C infection, calibrated the recency period, and investigated potential associations between the recency period and the viral RNA load.
Study Subjects and Testing
A total of 70 subjects including 8 acute and 62 recent cases from a primary HIV-1 subtype C infection study cohort in Botswana22 who were enrolled from March 2004 to December 2007 and who completed at least 4 study visits were included in the current analysis. There were 17 male (2 acute) and 53 female (6 acute) participants. Ages ranged from 19 to 53 years. The age distributions of acutely and recently infected participants were similar, as were those of males and female. All subjects were Botswana nationals. All infections were caused by HIV-1 subtype C.23 The study was approved by Institutional Review Boards in Botswana and the United States. Written informed consent was obtained from each participant.
Acutely or recently infected subjects were identified in a prospective study cohort of postnatal HIV-negative women or through voluntary testing and counseling-based referrals.22 Participants in the prospective study were screened for HIV infection bimonthly. Referred subjects were screened for HIV infection cross sectionally. The RNA+Ab− status of acutely infected subjects was defined by a positive HIV-1 reverse transcriptase-polymerase chain reaction test combined with a negative HIV-1 serology in double enzyme immunoassay (EIA). Initial screening for HIV-1/2 antibodies was performed with rapid EIA (Determine HIV-1/2, Abbott and Uni-Gold HIV Kit; Trinity Biotech, Wicklow Bay, Ireland) and/or regular EIA (Murex HIV 1.2.0, Abbott and Ortho HIV-1/2 Ab Capture; Ortho Diagnostics, Rochester, NY). The “window period” of the third-generation immunoassays used in the study was estimated to be about 2 weeks.24 Recent HIV-1 infections were identified by applying a 2-step testing algorithm, the Serologic Testing Algorithm for Recent HIV Seroconversion2 using the Vironostika HIV-1 Plus O Microelisa System (bioMérieux, Durham, NC) according to the protocol described elsewhere.3,4 A total of 615 samples from 51 subjects with multiple measurements of SOD were analyzed.
Acutely infected subjects had weekly visits for the first 2 months, biweekly visits for the next 2 months, and monthly visits for the following 8 months. Recently infected subjects had monthly visits. After the first year, study visits for both acutely and recently infected subjects became quarterly. Quantification of viral RNA in plasma was performed by COBAS AmpliPrep/COBAS Amplicor HIV-1 Monitor Test, ver. 1.5, according to the manufacturer's instructions. The level of detection ranged from 50 copies per milliliter for the ultrasensitive method and 400 copies per milliliter for the standard method to 750,000 copies per milliliter. Specimens exceeding the upper level of detection were retested by 10-fold dilutions. Quantification of the total proviral DNA was performed by the method described elsewhere25 by targeting HIV-1 gag p24 and using the TaqMan Universal PCR master mix on the Applied Biosystems 7500 Real-time PCR System (Foster City, CA). Quantification of CD4+ cells was performed by flow cytometry using the 4-color FACSCalibur. Individuals whose CD4+ T-cell count dropped below 200 cells per cubic millimeter or who developed opportunistic infection had access to ART [Combivir (AZT/3TC) 300/150 mg twice a day plus nevirapine 200 mg twice a day if female or efavirenz 600 mg every day if male] free of charge, in accordance with Botswana National Program guidelines.
Adjustment of Time 0 by Fiebig Stage
Subjects with acute infection were RNA+Ab− and the estimated time of their seroconversion was computed as a midpoint between their last Ab− and Ab+ tests (average time difference = 19.5 days). To estimate the seroconversion time point for the recently infected subjects, the stage of HIV infection was categorized by applying the Fiebig algorithm.26 The earliest available plasma specimens were analyzed by HIV Blot 2.2 Western Blot Assay (Genelabs Diagnostics S.A., Geneva, Switzerland). The distribution of 62 recent cases by Fiebig staging26 included 10 cases (16.1%) in stage IV, 18 cases (29.0%) in stage V, 4 cases (6.5%) on the edge of stages V and VI (extremely faint p31 band evident for a transition from stage V to stage VI), and 30 cases (48%) in stage VI. According to the original article26 and application of Fiebig staging to HIV-1 subtype C samples,27 the assumption was made that the beginning of Fiebig stage III coincides with the time of detectable seroconversion (time 0), and the mean duration of Fiebig stage III is 3 days, 6 days for stage IV, and 70 days for stage V. The mean duration of stage VI is not reliably estimable. Based on these estimates, the time from seroconversion until detection was assumed to average 6 days for subjects in stage IV (3 days of phase III and 3 days to the midpoint of phase IV), 44 days for subjects in stage V (9 days of phases III and IV and 35 days to the midpoint of phase V), and 79 days for stage V/VI (9 days of phases III and IV and 70 days of phase V). Because the adjustment for stage VI was not considered reliable because of the uncertainty associated with the estimated recency period, subjects identified within Fiebig stage VI were excluded from analyses related to identification of the recency period.
Data were summarized with means ± SDs or means ± 95% confidence intervals (CIs). For subjects whose SOD values had not crossed 1, the time to SOD crossing 1 was assumed to be right censored at the time of the last follow-up. Kaplan-Meier curves were used to describe the proportion of subjects whose SOD values remained below 1 during the 400 days after seroconversion, separately for 2 groups defined by RNA values, with the Gehan-Wilcoxon test used for comparisons. Mean time to SOD crossing a specified level, in the presence of censoring observations, was obtained by calculating the area under the Kaplan-Meier curve, assuming an exponential tail in situations when the Kaplan-Meier estimator did not reach 0. All reported P values are 2 sided.
Profile of SOD Increase in Primary HIV-1 Subtype C Infection
As seen in Figure 1, antiviral antibody titers quantified by detuned EIA were uniformly close to 0 during the early post-seroconversion period (when available), but thereafter demonstrated substantial heterogeneity in rates of growth in acutely (Fig. 1A) and recently (Fiebig stage IV, Fig. 1B and Fiebig stages V and V/VI, Fig. 1C) infected subjects. As a result, subjects' values reached thresholds (such as 1.0) at substantially different time points. The profiles of SOD values demonstrated slow initial increase soon after seroconversion and were followed by a more rapid increase in most cases. However, in 2 of 29 (6.9%) subjects identified in Fiebig stage V/VI or earlier, SOD values did not reach the threshold of 1.0 during the follow-up period (695 days for subjects 3244 and 195 days for subjects 4872). The profiles of SOD evolution in 22 subjects identified within Fiebig stage VI (Fig. 1D) were similar to profiles in subjects with earlier Fiebig stages.
To characterize the profile of SOD dynamics in primary HIV-1 subtype C infection, the SOD values were analyzed within 50-day intervals post seroconversion separately in groups of acutely and recently infected subjects. Similar profiles of SOD kinetics were evident within groups up to 300 days post seroconversion followed by heterogeneity afterward (data not shown). This similarity supports the value of using Fiebig staging to adjust for average seroconversion time.
Recency Period Is a Time of Reaching SOD Threshold
To analyze the kinetics of SOD values, we computed the estimated time of crossing specified SOD levels from 0.1 to 1.0 by linearly interpolating between the SOD values observed just before and after reaching the specified SOD value. The heterogeneity in reaching specified SOD levels between acutely infected subjects is shown in Table 1. A generally gradual increase of the duration from the estimated time of seroconversion to the estimated time of reaching the specified SOD levels from 0.1 to 1.0 was accompanied by wide 95% CIs reflecting substantial variance between subjects. For the group of acutely infected subjects, the SOD value of 1.0 was reached in an average of 155 days (95% CI: 119 to 191 days) post seroconversion.
A similar analysis of the estimated time of reaching different SOD levels was performed for groups of recently infected subjects identified within Fiebig stages IV and V-V/VI (data not shown). The mean time to reach the SOD threshold of 1.0 in the group of recently infected subjects, n = 20, was similar to the acute infection group: 149.7 days (95% CI: 124 to 176 days). For the entire cohort, the mean time to reach the SOD threshold of 1.0 was 151 days (95% CI: 130 to 172 days). The SOD values for subjects 3244 increased at a very slow rate and remained low during an extended period; as a result, his/her time to crossing each specified SOD level was markedly higher than those of the remaining 28 subjects and heavily influenced the results. We therefore excluded his/her data in estimating mean time to crossing each specified SOD level, and note that this exclusion is likely to result in an underestimation of the recency periods.
Virologic and Immunologic Parameters in Subjects 3244 and 4872
To investigate why subjects may have continuously low SOD values, we analyzed viral RNA dynamics in plasma, cell-associated proviral DNA, and CD4+ T-cell counts in subjects 3244 and 4872 for 695 days and 195 days, respectively (Fig. 2). Both subjects had low levels of viral RNA and proviral DNA and relatively high levels of CD4+ T cells with remarkable fluctuations. A transient spike of viral RNA in subjects 4872 at day 108 did not result in changes in other studied parameters. Although no single underlying reason could be identified, it seems that efficient control of the viral replication and preservation of CD4+ T cells (or their fast recovery) might contribute to the observed delay of the anti-HIV antibody titers increase, and, in turn, present a potential problem for the incidence analysis based on cross-sectional sampling.
Association of Early Viral Replication With Recency Period
We assessed whether viral RNA load within the early post-seroconversion period of HIV-1C is associated with recency period. We computed each subject's average viral load in the first 79 days after seroconversion (which corresponds to the end of Fiebig stage V) and their estimated time post seroconversion of reaching an SOD threshold of 1.0. The resulting scatter plot (Fig. 3A) suggests that subjects with higher early viral load levels have shorter recency periods (R 2 = 0.419, P < 0.001).
Second, we classified each subject as having a high (above median) or low (below median) early viral load and used time-to-event methods to assess the association between viral load group (high versus low) and time until reaching an SOD threshold of 1.0. We performed 3 analyses based on different ways of defining “high” and “low” viral load and how the date of seroconversion was defined. In all analyses, viral load groups were compared using the generalized Wilcoxon test to account for censored data (subjects whose SODs never crossed 1). In the first analysis (Fig. 3B), days from seroconversion were calculated based on midpoint adjustment. That is, the time from seroconversion for a subject found by Western blot to be in a specific Fiebig stage was taken to be the midpoint of this Fiebig interval. The subject's early viral load was then defined as the average of her/his RNA values from estimated seroconversion to the end of Fiebig stage V and classified as high (or low)-RNA group if his/her average RNA value during this period was greater (or smaller) than the median value across subjects. By these definitions, recency period was significantly shorter for subjects in the high-RNA group than for those in the low-RNA group (P = 0.007). In the second analysis, a simulation was performed with randomly chosen time for each subject. As in the first analysis, RNA group was significantly associated with the recency period (the median P value obtained from 2000 such analyses was 0.01, data not shown). In the third analysis, viral load values for each subject were randomly chosen among all available RNA measurements from seroconversion to the end of Fiebig Stage V, as opposed to the mean viral load value in this period in the first and second analyses. As in the first 2 analyses, RNA group was significantly associated with time to SOD crossing 1 (P = 0.003, data not shown).
Overall, the observed patterns suggest that early levels of viral RNA load are associated with the recency period. It is likely that the level of viral replication reflected by the level of viral RNA load drives the recency period, which suggests that viral load might be usefully incorporated into methods for estimating HIV incidence rates.
Reduction of SOD Values in ARV-Treated subjects
Longitudinal SOD measurements were available in 8 of 13 subjects who initiated ART (Fig. 4, filled circles). The reduction of SOD values after starting ART was a dominant pattern, despite high heterogeneity among subjects in the time and level of the drop, or subsequent increase of SOD values (eg, subjects 2767 and 2990). It is important to determine whether initiation of ART is followed by reduction of SOD values below the threshold, and if so, what the fraction of those subjects is in subtype C infection. In our study, 4 of 8 studied subjects dropped their SOD values below the threshold of 1.0, which may be due, at least in part, to a short period of follow-up after initiation of ART in subjects 3312 and 3481. However, subjects 2767 and 2990 demonstrated a transient drop that did not reach the threshold, followed by an increase in SOD values. Surprisingly, the expected increase of viral load was not observed and SOD increase was accompanied by suppressed viral RNA load below detectable levels and by stable proviral DNA load in both subjects.
The goal of the study was to characterize SOD evolution during the first year of infection in a cohort of subjects identified with primary HIV-1 subtype C infection and to estimate the recency period, expressed as a duration from the estimated time of seroconversion to the estimated time of crossing the SOD threshold of 1.0 as measured by detuned assay with the Vironostika EIA kit. Western blot kinetics in sequential samples was found to be helpful in identifying recently infected subjects and estimating time of seroconversion. The SOD in primary HIV-1 infection showed gradual increase over time in most analyzed cases. It has been suggested that decline in any 2 consecutive SOD values indicates an old infection.28 Our data (Fig. 1) suggest that fluctuations in SOD values are relatively common within the early phase of HIV-1 subtype C infection. The duration of the recency period with SOD cutoff of 1.0 was estimated to be 151 days (95% CI: 130 to 172 days) post seroconversion.
We found that the recency period was significantly associated with early post-seroconversion viral load. Subjects with higher RNA load reached the SOD threshold earlier than subjects with low viral load. In most ARV-treated subjects, the drop of viral load was followed by reduction of SOD values. Although viral load is a major marker of viral replication, the recency period mirrors the increase of virus-specific antibody titers. It seems plausible that the level of viral replication after seroconversion in HIV-1 subtype C infection determines the rate of antiviral antibodies elicited. The diverse levels of viral RNA in different settings may contribute to the differential performance of detuned EIA tests in non-B subtype epidemics. However, the higher levels of early viral load reported in subtype E-infected individuals (as compared with subtype B infection) in Thailand29 are in contrast with the longer recency period of HIV-1 subtype E,28 suggesting a need for dedicated studies able to address relationships between viral load and recency period.
We used Fiebig staging to determine the average time of seroconversion for individuals detected with recent infection, enabling the estimation of the average recency period until an SOD threshold. The notable similarity between average recency times for these subjects and for those with acute infections, where seroconversion time was determined directly, supports this use of Fiebig staging in future studies. The finding that viral load is significantly associated with recency duration suggests that estimation of HIV incidence using cross-sectional samples of standard and detuned EIA, or the BED assay, might be improved if information about viral load and/or Western blot could be incorporated into the analysis.
We note that the average recency periods obtained in this study, when used in conjunction with cross-sectional (sensitive/less sensitive) antibody testing to estimate HIV incidence rates, depend critically on the test used. Thus, the estimate of 151 days, which was obtained based on Vironostika EIA, could not be used directly in incidence estimation if another test was used to determine whether subjects were in recency period. However, the key findings of this study regarding the use of Fiebig staging to estimate average recency period, and the significant association between early viral load and recency period, would be expected to hold regardless of the specific testing algorithm used.
In summary, our study determined the average recency period in the local HIV-1 subtype C epidemic in Botswana to be 151 days (95% CI: 130 to 172 days, using SOD cutoff of 1.0) post seroconversion and provided a detailed picture of the evolution of SOD values in the early stage of HIV-1 subtype C infection. The high intersubject variability of SOD values presents a significant challenge for accurate estimation of the recency period and HIV incidence in a heterogeneous population. We found that the level of viral replication during the early post-seroconversion period is significantly associated with the time that titers of antiviral antibodies quantified as SOD values reach the threshold.
We are grateful to the subjects from the “Tshedimoso” study in Botswana. We thank Gaseboloke Mothowaeng, Florence Modise, S'khatele Molefhabangwe, and Sarah Masole for their dedication and outstanding work in the clinic and outreach. We express thanks to Busisiwe Mlotshwa for excellent laboratory support. We greatly appreciate the enthusiasm and strong commitment of Erin McDonald, Melissa Ketunuti, Carl Davis, Kenneth Onyait, and Mary Fran McLane in achieving the overall study goals. We thank the Botswana Ministry of Health, Gaborone City Council clinics, and the Gaborone voluntary testing and counseling Tebelopele for their ongoing support and collaboration. We thank Lendsey Melton for excellent editorial assistance. Author Contributions: V.N. and M.E. designed and guided the study and wrote the article. E.W. contributed to patient recruitment and provided clinical support for the cohort. L.K., J.G., and R.R. performed the experiments. S.M. and R.M. provided laboratory support. R.W. and S.L. designed and performed statistical analysis and wrote the article.
1. Brookmeyer R, Quinn TC. Estimation of current human immunodeficiency virus incidence rates from a cross-sectional survey using early diagnostic tests. Am J Epidemiol
2. Janssen RS, Satten GA, Stramer SL, et al. New testing strategy to detect early HIV-1 infection for use in incidence estimates and for clinical and prevention purposes. JAMA
3. Rawal BD, Degula A, Lebedeva L, et al. Development of a new less-sensitive enzyme immunoassay for detection of early HIV-1 infection. J Acquir Immune Defic Syndr
4. Kothe D, Byers RH, Caudill SP, et al. Performance characteristics of a new less sensitive HIV-1 enzyme immunoassay for use in estimating HIV seroincidence. J Acquir Immune Defic Syndr
5. Parekh BS, Kennedy MS, Dobbs 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
6. Suligoi B, Galli C, Massi M, et al. Precision and accuracy of a procedure for detecting recent human immunodeficiency virus infections by calculating the antibody avidity index by an automated immunoassay-based method. J Clin Microbiol
7. Soroka SD, Granade TC, Candal D, et al. Modification of rapid human immunodeficiency virus (HIV) antibody assay protocols for detecting recent HIV seroconversion. Clin Diagn Lab Immunol
8. Wilson KM, Johnson EI, Croom HA, et al. Incidence immunoassay for distinguishing recent from established HIV-1 infection in therapy-naive populations. AIDS
9. Barin F, Meyer L, Lancar R, et al. Development and validation of an immunoassay for identification of recent human immunodeficiency virus type 1 infections and its use on dried serum spots. J Clin Microbiol
10. Li H, Ketema F, Sill AM, et al. A simple and inexpensive particle agglutination test to distinguish recent from established HIV-1 infection. Int J Infect Dis
11. Cleghorn FR, Jack N, Murphy JR, et al. Direct and indirect estimates of HIV-1 incidence in a high-prevalence population. Am J Epidemiol
12. Hargrove JW, Humphrey JH, Mutasa K, et al. Improved HIV-1 incidence estimates using the BED capture enzyme immunoassay. AIDS
13. Lagakos S, Gable A, eds. Methodological Challenges in Biomedical HIV Prevention Trials
. Washington, DC: National Academy Press; 2008;175-185.
14. McDougal JS, Parekh BS, Peterson ML, 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
15. Karita E, Price M, Hunter E, et al. Investigating the utility of the HIV-1 BED capture enzyme immunoassay using cross-sectional and longitudinal seroconverter specimens from Africa. AIDS
16. Laeyendecker O, Rothman RE, Henson C, et al. The effect of viral suppression on cross-sectional incidence testing in the Johns Hopkins hospital emergency department. J Acquir Immune Defic Syndr
17. Morris MK, Katzenstein DA, Israelski D, et al. Characterization of the HIV-1 specific humoral immune response during highly active antiretroviral therapy (HAART). J Acquir Immune Defic Syndr
18. Gupta SB, Murphy G, Koenig E, et al. Comparison of methods to detect recent HIV type 1 infection in cross-sectionally collected specimens from a cohort of female sex workers in the Dominican Republic. AIDS Res Hum Retroviruses
19. Buchacz K, Klausner JD, Kerndt PR, et al. HIV incidence among men diagnosed with early syphilis in Atlanta, San Francisco, and Los Angeles, 2004 to 2005. J Acquir Immune Defic Syndr
20. Priddy FH, Pilcher CD, Moore RH, et al. Detection of acute HIV infections in an urban HIV counseling and testing population in the United States. J Acquir Immune Defic Syndr
21. Sakarovitch C, Rouet F, Murphy G, et al. Do tests devised to detect recent HIV-1 infection provide reliable estimates of incidence in Africa? J Acquir Immune Defic Syndr
22. Novitsky V, Woldegabriel E, Wester C, et al. Identification of primary HIV-1C infection in Botswana. NIHMSID # 79283. AIDS Care
23. Novitsky V, Woldegabriel E, Kebaabetswe L, et al. Viral load and CD4+ T cell dynamics in primary HIV-1 subtype C infection. J Acquir Immune Defic Syndr
24. Speers D, Phillips P, Dyer J. Combination assay detecting both human immunodeficiency virus (HIV) p24 antigen and anti-HIV antibodies opens a second diagnostic window. J Clin Microbiol
25. Novitsky VA, Gilbert PB, Shea K, et al. Interactive association of proviral load and IFN-gamma-secreting T cell responses in HIV-1C infection. Virology
26. Fiebig EW, Wright DJ, Rawal BD, et al. Dynamics of HIV viremia and antibody seroconversion in plasma donors: implications for diagnosis and staging of primary HIV infection. AIDS
27. Salazar-Gonzalez JF, Bailes E, Pham KT, et al. Deciphering human immunodeficiency virus type 1 transmission and early envelope diversification by single-genome amplification and sequencing. J Virol
28. Young CL, Hu DJ, Byers R, et al. Evaluation of a sensitive/less sensitive testing algorithm using the bioMerieux Vironostika-LS assay for detecting recent HIV-1 subtype B or E infection in Thailand. AIDS Res Hum Retroviruses
29. Hu DJ, Vanichseni S, Mastro TD, et al. Viral load differences in early infection with two HIV-1 subtypes. AIDS