Arroyo, Miguel A PhD*†; Sateren, Warren B MPH*; Serwadda, David MD, MMed, MPH‡; Gray, Ronald H MD, MSc, DTM&H§; Wawer, Maria J MD, MHSc∥; Sewankambo, Nelson K MD, MMed, MSc‡; Kiwanuka, Noah MBChB‡; Kigozi, Godfrey MBChB, MPH‡; Wabwire-Mangen, Fred MD, PhD‡; Eller, Michael MSc¶; Eller, Leigh Anne MSc¶; Birx, Deborah L MD*; Robb, Merlin L MD#; McCutchan, Francine E PhD#
From the *US Military HIV Research Program (USMHRP), Walter Reed Army Institute of Research, Division of Retrovirology, Silver Spring, MD; †Department of Retrovirology, Armed Forces Research Institute of Medical Sciences, Bangkok, Thailand; ‡Makarere University, Kampala, Uganda; §Johns Hopkins University, Baltimore, MD; ∥Columbia University, New York, NY; ¶Makarere University Walter Reed Research Project, Henry M. Jackson Foundation, Kampala, Uganda; and #USMHRP, Henry M. Jackson Foundation, Rockville, MD.
Received for publication January 10, 2006; accepted July 5, 2006.
Supported by a cooperative agreement between the Henry M. Jackson Foundation for the Advancement of Military Medicine and the United States Department of Defense. The US Military HIV Research Program is jointly planned and funded by the Department of Defense and National Institute of Allergy and Infectious Diseases, National Institutes of Health, by an interagency agreement.
This manuscript was approved for publication by the Walter Reed Army Institute of Research. The opinions and assertions made by the authors do not reflect the official position or opinion of the US Department of Defense, the US Department of Army, the US Military HIV Research Program, the Walter Reed Army Institute of Research, or the Henry M. Jackson Foundation for the Advancement of Military Medicine, Inc.
Reprints: Francine E. McCutchan, PhD, USMHRP, Henry M. Jackson Foundation, 1600 East Gude Drive, Rockville, MD 20850 (e-mail: firstname.lastname@example.org).
The molecular epidemiology of HIV-1 reveals a complex evolving pandemic, composed of 9 genetic subtypes, 20 or more circulating recombinant forms (CRFs), and myriad unique recombinant forms (URFs).1 Viral variation must be accurately documented to assess its impact on treatment prevention,2 and the social and geographic dimensions of the epidemic are particularly important for the planning of HIV-1 vaccine trials.3
In many potential vaccine trial sites in Africa, intersubtype recombination is an important factor influencing the complexity of HIV-1 strains. The origin and genesis of recombinant strains are being actively investigated to gain insight into associated population demographic variables. In urban and rural areas of the Mbeya region in Tanzania, an association between HIV-1 prevalence and viral diversity has recently been reported.4 Individuals living in high-prevalence urban areas were more likely to be infected with a recombinant strain and to be dually infected compared with individuals from a lower prevalence rural village. The urban community had greater access to the TransAfrican highway, where residents often engage in social interactions, with increased opportunity for casual sexual encounters between travelers and local inhabitants.5,6 In the Rakai District of Uganda, previous studies have also reported differences in HIV-1 prevalence and risk factors associated with the sexual transmission of HIV-1 within geographic strata of differing accessibility to main roads.7,8
A previous study in Rakai documented that subtype D accounts for more than 50% of infections and subtype A for <20%, whereas almost 30% of HIV-1 infections in Rakai are intersubtype recombinant strains.9 We hypothesized that communities along main roads in Rakai District, compared with their more rural and less accessible counterparts, would provide an environment for increased social and sexual interaction, higher HIV-1 incidence, increased importation and transmission of novel HIV-1 strains, and increased risk for multiple exposures to HIV-1, leading to dual infection and generation of intersubtype recombinant strains. Here, we report an analysis of HIV-1 incidence and demographic variables in 10 communities in the Rakai District of Uganda and the genotype of HIV-1 strains in 329 recently infected individuals from these communities.
The population in Rakai District has been extensively studied for HIV incidence and risk factors since 1988 through a series of community cohort studies. These cohort studies are annual surveys of selected trading center communities and rural village communities representing a cross section of the Rakai population. An average of 19,228 volunteers participated annually. Individuals who had a negative HIV test result followed by a positive HIV test result within any of 3 successive cohort studies were recruited to the Molecular Epidemiology (MER) study within 1 month of their first positive HIV test result. The MER study recruited 80 of these incident cases from the 1997 through 1998 study interval, 188 from 1999 through 2002, and 140 from 2002 through 2003. Among this total population of recent seroconverters (n = 408), 384 had a cryopreserved plasma sample from entry into the MER study available for analysis.
MER study participants were offered posttest counseling by the cohort staff and were referred to MER study staff for study-specific informed consent, which was obtained in Lugandan. This research project and the consent documents were approved by the Uganda Viral Research Institute Institutional Review Board (IRB), the IRBs of Columbia University and Johns Hopkins School of Public Health, and the US Army Surgeon General's Human Subjects Research Review Board.
MER study participants resided in 1 of 10 different communities within the Rakai District as previously described by Wawer et al7 (Fig. 1). The main road trading centers (Katana, Kalisizo, Kyotera, Kasasa-Sanje, and Kakuuto) included shops, bars, and hotels, all serving domestic and international traffic. In contrast, the villages (Lwamaggwa, Kabira, Buyamba, Lwanda, and Kibale-Rakai) are 2 to 10 km from main roads and accessible only via gravel or dirt roads.
For determination of HIV-1 serostatus within the cohort studies, sera were tested with dual enzyme-linked immunosorbent assay (ELISA) tests (Vironostika, HIV Uni-Form II Test, Durham, NC; Abbott Wellcozyme HIV Recombinant Test [Abbott Laboratories, North Chicago, IL]) and were confirmed by HIV-1 Western blot (Calypte/Cambridge Biotech HIV-1, Berkeley, CA).
HIV-1 incidence rates were calculated only for the cohort study from 1999 through 2002, because the remaining 2 studies are ongoing or part of separate reports. This was part of the MER study; it was called the Community HIV Epidemiology Research (CHER) study and was similarly conducted, with all participants providing informed consent under protocols approved by local and US IRBs as described previously. In this study, 188 incident cases were identified during 3 annual surveys in the 10 study communities. The total number of person-years (PYs) was calculated separately for 2 groups: based on interval of years between the first and last visits for HIV-1-negative individuals and based on seroconverters contributing person-time from their last negative HIV test result to their first positive HIV test result. We defined this as the midpoint between the 2 dates, and seroconverters were presumed to contribute person-time to the cohort follow-up from their first negative HIV test result to the midpoint date of their seroconversion interval. The incidence rates presented are unadjusted, because the calendar period was too short to determine any meaningful trends (3 years), with numerous individuals participating in only the first 2 years or the last 2 years. Incidence rates were calculated as the number of seroconversions per 100 PYs of follow-up, along with the 95% confidence interval (CI).
Multiregion Hybridization Assay version 2
Viral RNA was extracted from 200 μL of plasma using the MagNA Pure Total Nucleic Acid robotic extraction procedure (Roche Diagnostics Corporation, Indianapolis, IN). As negative controls, phosphate-buffered saline (PBS) and HIV-1-seronegative samples were included in each extraction procedure and in the multiregion hybridization assay (MHAacd) analysis. The MHAacd was performed as described by Arroyo et al.4,10,11
The 2 geographic strata were compared for population demographic characteristics and HIV subtype proportions using the χ2 test or Fisher exact test. Only significant P values <0.05 are shown. Adjusted odds ratios (AORs) were calculated, adjusting for those variables that were statistically significantly predictive of HIV risk on bivariate analysis, gender, age, strata, and seroconversion time period so as to determine any associations with HIV-1 subtype. Cox proportional hazards regression analysis was used to calculate HIV-1 subtype incidence hazard ratios (HRs) along with their 95% CIs. All statistical analyses were performed using JMP 5.1 software (SAS Institute, Cary, NC).
Characteristics of Molecular Epidemiology Study Participants
A total of 408 individuals with incident HIV infection observed within ongoing community cohort studies from 10 communities in Rakai District were enrolled into the MER study over the period 1997 through 2003 (Fig. 1). Two hundred thirty-five (58%) resided in 1 of 5 trading centers along the TransAfrican highway leading to Tanzania, whereas 173 (42%) came from trading villages along secondary roads. Of all the MER study volunteers, 246 (60%) were women (56% lived in trading centers and 44% in trading villages) and 162 (40%) were men (60% lived in trading centers and 40% in trading villages). Each community enrolled male and female participants in similar numbers, except Kakuuto, which enrolled fewer female participants (P = 0.05). There was no significant difference in gender ratio between strata (trading center vs. trading village). The age of enrollees at first enrollment over the 3-year period ranged from 14 to 90 years, with a mean age of 28.5 years and a median age of 26 years; the interquartile range was 21 to 34 years. HIV-positive volunteers were significantly (P < 0.05) older than HIV-negative volunteers by 2.5 years (31.1 vs. 28.6 years, respectively). Men and women had a similar age distribution (mean age of 28.5 years, SD = 8.0 for women, mean age of 31.44 years, SD = 6.9 for men) in the 2 strata.
HIV-1 Subtype Distribution
Plasma samples collected within 1 month of the first HIV-1-seropositive sample were available from 384 (94%) of 408 participants, and of these, 329 (86%) were successfully genotyped using the MHAacd (Table 1). Overall, the cohort had 52 (16%) subtype A, 189 (57%) subtype D, and 87 (26%) recombinant HIV-1 strains. The only subtype C strain detected was from Kyotera main road trading center, representing <1% of the total subtype distribution. The distribution among the recombinant strains was 3 (3%) AC, 68 (78%) AD, 14 (16%) CD, and 2 (2%) ACD. Among all samples genotyped, 273 (83%) were subtype D or a subtype-D-containing recombinant strain. The proportions of subtype D differ by 8% between the 2 geographic strata (subtype D, 54% vs. 62%). Kakuuto had more subtype D (71%) and Kibale-Rakai had more subtype A (29%) than other communities. Main road trading centers had more recombinant strains than trading villages (32% vs. 20%; P < 0.05). Recombinants constituted 32% to 35% of strains in 4 trading centers, but Kakuuto was the exception, with only 14% recombinants. In trading villages, the proportion of recombinant strains averaged 20%, ranging from 15% to 23%.
Dual infections, ascertained by 2 or more subtype-specific probes hybridizing in a given genome region, were found in 23 (7%) of participants (Table 1) at rates of 7.5% in trading centers and 6.3% in trading villages. Most dual infection cases detected occurred between AD (87%) subtypes, whereas only a few occurred between CD (9%) and ACD (4%) subtypes. These results showed an increase in the number of HIV recombinant strains and dual infections in trading centers, which was also associated with a decrease in the number of subtype D strains. These differences represent an increase in the genetic complexity of HIV strains circulating in this stratum. To determine the significance of these differences between the 2 strata, we categorized HIV-1 genotypes in subsequent analyses as subtype D versus non-subtype D, pure subtype versus recombinant, or dual versus single infection, respectively.
Multivariate Regression Analysis
Multivariate regression analysis was performed to calculate the odds ratio for a particular HIV-1 subtype among those MER study participants who were successfully genotyped (n = 329). The unadjusted odds ratios and AORs for recombinant, nonsubtype D, and dual infections, by gender and geographic stratum, are shown in Table 2. The analysis confirms increased odds ratios for recombinant HIV-1 strain infection among men and women in trading centers compared with trading villages (AOR among women = 1.7, 95% CI: 0.9 to 3.5; AOR among men = 2.0, 95% CI: 0.8 to 4.8). The risk of non-subtype D infections was increased among men in trading centers (AOR = 1.7, 95% CI: 0.8 to 3.6), whereas among women in trading centers, there was also an increase in the risk of a dual HIV-1 infection (AOR = 1.8, 95% CI: 0.5 to 7.8). These results are noteworthy, but the lack of statistical significance might be accounted for by the small number of samples analyzed.
The CHER study conducted between 1999 and 2002 identified 188 incident cases (Fig. 1). HIV-1 incidence rates stratified by gender and geography are shown in Table 3. The overall incidence was 1.3 per 100 PYs. Women had a higher HIV-1 incidence rate than men (1.3 vs. 1.2 per 100 PYs). The HIV-1 incidence rate was higher in the main road trading centers compared with trading villages in men and women (men: 1.3 vs. 1.2 per 100 PYs; women: 1.3 vs. 1.1 per 100 PYs). The trading center incidence was also higher than that of trading villages overall (1.3 vs. 1.1 per 100 PYs). The HIV-1 incidence was significantly higher among individuals who moved between strata (6.1 per 100 PYs) versus those who remained in trading villages (1.1 per 100 PYs; P < 0.01) or trading centers (1.3 per 100 PY; P < 0.01) at any time during the CHER study. This migration-associated increased incidence was higher among women (9.2 per 100 PYs; P < 0.01) than men (1.8 per 100 PYs).
Hazard Ratios of HIV-1 Infection by Genotype and Geographic Stratum
The HRs for HIV-infection with a particular HIV-1 strain (subtype D, pure subtype [A, C, or D], non-D subtype, and recombinant) were calculated among those individuals who seroconverted (n = 188) during 1999 through 2002 (Table 4). The risk for subtype D infection varied little by gender or community type. Considering pure subtype infections (A, C, or D strains) collectively, however, women who moved between stratums had a higher risk (HR = 5.7, 95% CI: 1.4 to 16.2). Similarly, grouping the non-D strains together (A, C, and recombinant) showed elevated risk among women in trading centers (HR = 2.3, 95% CI: 1.0 to 5.9), and especially among those who were mobile between strata (HR = 20.6, 95% CI: 4.4 to 74.7).
The risk of a recombinant infection was also elevated in women living in trading centers (HR = 5.3, 95% CI: 1.5 to 34.2) or in those who moved between strata (HR = 22.9, 95% CI: 1.1 to 241.0) compared with women who lived in rural trading villages. Of interest, this disproportionate risk for acquisition of certain HIV-1 strains, namely, those other than the majority subtype D strain, was not observed among men.
Through the development and application of high-throughput genotyping assays like the MHAacd, the social and geographic dimensions of HIV-1 epidemics can be described with new clarity and the molecular epidemiology of HIV-1 in cohorts with different levels of risk can be assessed. This study adds to a growing body of evidence that the risk for HIV-1 infection, and the social milieu in which infection takes place, can have an impact not only on HIV-1 incidence but on the genetic complexity of HIV-1 strains in the population. In a larger study conducted previously throughout Rakai District, the association of higher HIV-1 prevalence with proximity to paved roads was demonstrated, suggesting that HIV-1 transmission tends to follow lines of transportation, communication, and social interaction.7,8 A recent study by Arroyo et al4 in the Mbeya region of Tanzania reiterated this finding and added another dimension, associating higher prevalence with a higher genetic complexity of strains. Other studies have also found this relation when comparing high- and low-risk cohorts.5,12 This study confirms an association between increased exposure to HIV-1 strains and the development of HIV-1 genetic complexity.
One correlate of these observations could be a decline in the genetic complexity of HIV-1 when prevalence rates fall. Yirrell et al,13 studying another region of Uganda, recently reported a parallel decline in the prevalence and fraction of recombinant strains over 10 years. The low HIV-1 incidence reported in the CHER study was in accordance with previous reports in southwest Uganda.14,15 The challenge of depicting differences in HIV-1 incidence in a low-risk population might have accounted for the lack of statistical significance observed. Despite this fact, we were able to observe subtle differences in HIV-1 incidence in different population strata, and these were associated with differing HIV-1 subtype complexity.
With respect to planning for vaccine trials, the results reported here as well as others reported previously confirm that there is considerable heterogeneity with respect to incidence and subtype complexity within potential source populations for vaccine trials. Finding the optimal balance between high HIV-1 incidence, which simplifies trials, and the complexity of the HIV-1 strains that challenge volunteers, which may complicate the interpretation of trial results, can be facilitated by studies that reveal the complex interactions between gender, social setting, and mobility in potential trial sites.
This study was approved by the following IRBs: the Office of the Deputy Chief of Staff for Regulatory Compliance and Quality, US Army Medical Research and Material Command; the Columbia University IRB; the Johns Hopkins Committee on Human Research; the AIDS Research Subcommittee of the Uganda National Council for Science and Technology; the Uganda Ministry of Health; and the US Army Medical Research Unit in Kampala, Uganda.
1. McCutchan FE. Understanding the genetic diversity of HIV-1. AIDS
. 2000;14(Suppl 3):S31-S44.
2. Peeters M, Sharp PM. Genetic diversity of HIV-1: the moving target. AIDS
. 2000;14(Suppl 3):S129-S140.
3. Suligoi B, Wagner TM, Ciccozzi M, et al. The epidemiological contribution to the preparation of field trials for HIV and STI vaccines: objectives and methods of feasibility studies. Vaccine
4. Arroyo MA, Hoelscher M, Sateren W, et al. HIV-1 diversity and prevalence differ between urban and rural areas in the Mbeya region of Tanzania. AIDS
5. McCutchan FE, Hoelscher M, Tovanabutra S, et al. In-depth analysis of a heterosexually acquired human immunodeficiency virus type 1 superinfection: evolution, temporal fluctuation, and intercompartment dynamics from the seronegative window period through 30 months postinfection. J Virol
6. Hoffmann O, Zaba B, Wolff B, et al. Methodological lessons from a cohort study of high risk women in Tanzania. Sex Transm Infect
. 2004;80(Suppl 2):ii69-ii73.
7. Wawer MJ, Serwadda D, Musgrave SD, et al. Dynamics of spread of HIV-1 infection in a rural district of Uganda. BMJ
8. Serwadda D, Wawer MJ, Musgrave SD, et al. HIV risk factors in three geographic strata of rural Rakai District, Uganda. AIDS
9. Harris ME, Serwadda D, Sewankambo N, et al. Among 46 near full length HIV type 1 genome sequences from Rakai District, Uganda, subtype D and AD recombinants predominate. AIDS Res Hum Retroviruses
10. Arroyo MA, Hoelscher M, Sanders-Buell E, et al. HIV type 1 subtypes among blood donors in the Mbeya region of southwest Tanzania. AIDS Res Hum Retroviruses
11. Hoelscher M, Dowling WE, Sanders-Buell E, et al. Detection of HIV-1 subtypes, recombinants, and dual infections in East Africa by a multi-region hybridization assay. AIDS
12. Herbinger KHGM, Piyasirisilp S, Mloka D, et al. Frequency of HIV multiple infection and HIV diversity: analysis of low and high risk populations in Mbeya region, Tanzania. AIDS Res Hum Retroviruses
13. Yirrell DL, Kaleebu P, Morgan D, et al. HIV-1 subtype dynamics over 10 years in a rural Ugandan cohort. Int J STD AIDS
14. Mbulaiteye SM, Mahe C, Whitworth JA, et al. Declining HIV-1 incidence and associated prevalence over 10 years in a rural population in south-west Uganda: a cohort study. Lancet
15. Kamali A, Carpenter LM, Whitworth JA, et al. Seven-year trends in HIV-1 infection rates, and changes in sexual behaviour, among adults in rural Uganda. AIDS
© 2006 Lippincott Williams & Wilkins, Inc.