Variability of HIV-1 Genomes among Children and Adolescents from São Paulo, Brazil

Abstract

Background

Genetic variability is a major feature of the human immunodeficiency virus type 1 (HIV-1) and considered the key factor to frustrating efforts to halt the virus epidemic. In this study, we aimed to investigate the genetic variability of HIV-1 strains among children and adolescents born from 1992 to 2009 in the state of Sao Paulo, Brazil.

Methodology

Plasma and peripheral blood mononuclear cells (PBMC) were collected from 51 HIV-1-positive children and adolescents on ART followed between September 1992 and July 2009. After extraction, the genetic materials were used in a polymerase chain reaction (PCR) to amplify the viral near full length genomes (NFLGs) from 5 overlapped fragments. NFLGs and partial amplicons were directly sequenced and data were phylogenetically inferred.

Results

Of the 51 samples studied, the NFLGs and partial fragments of HIV-1 from 42 PBMCs and 25 plasma were successfully subtyped. Results based on proviral DNA revealed that 22 (52.4%) patients were infected with subtype B, 16 (38.1%) were infected with BF1 mosaic variants and 4 (9.5%) were infected with sub-subtype F1. All the BF1 recombinants were unique and distinct from any previously identified unique or circulating recombinant forms in South America. Evidence of dual infections was detected in 3 patients coinfected with the same or distinct HIV-1 subtypes. Ten of the 31 (32.2%) and 12 of the 21 (57.1%) subjects with recovered proviral and plasma, respectively, protease sequences were infected with major mutants resistant to protease inhibitors. The V3 sequences of 14 patients with available sequences from PBMC/or plasma were predicted to be R5-tropic virus except for two patients who harbored an X4 strain.

Conclusions

The high proportion of HIV-1 BF1 recombinant, coinfection rate and vertical transmission in Brazil merits urgent attention and effective measures to reduce the transmission of HIV among spouses and sex partners.

About the Authors:

Sabri Saeed Sanabani

Contributed equally to this work with: Sabri Saeed Sanabani, Rodrigo Pessôa

* E-mail: [email protected]

Affiliations Clinical and Research Laboratory (LIM 03), School of Medicine, University of São Paulo, São Paulo, Brazil, Virology Laboratory (LIM 52-HCFMUSP), Institute of Tropical Medicine, University of São Paulo, São Paulo, Brazil

Rodrigo Pessôa

Contributed equally to this work with: Sabri Saeed Sanabani, Rodrigo Pessôa

Affiliation: Virology Laboratory (LIM 52-HCFMUSP), Institute of Tropical Medicine, University of São Paulo, São Paulo, Brazil

Ana Carolina Soares de Oliveira

Affiliation: Virology Laboratory (LIM 52-HCFMUSP), Institute of Tropical Medicine, University of São Paulo, São Paulo, Brazil

Vanessa Pouza Martinez

Affiliation: Virology Laboratory (LIM 52-HCFMUSP), Institute of Tropical Medicine, University of São Paulo, São Paulo, Brazil

Maria Teresa Maidana Giret

Affiliation: Two Story Lab, Miller School of Medicine, University of Florida, Miami, Florida, United States of America

Regina Célia de Menezes Succi

Affiliation: Department of Pediatrics, Paulista School of Medicine, Federal University of São Paulo, São Paulo, Brazil

Karina Carvalho

Affiliation: Division of Clinical Immunology and Allergy, School of Medicine, University of São Paulo, São Paulo, Brazil

Claudia Satiko Tomiyama

Affiliation: Division of Clinical Immunology and Allergy, School of Medicine, University of São Paulo, São Paulo, Brazil

Douglas F. Nixon

Affiliation: Division of Experimental Medicine, Department of Medicine, University of California San Francisco, San Francisco, California, United States of America

Ester Cerdeira Sabino

Affiliation: Virology Laboratory (LIM 52-HCFMUSP), Institute of Tropical Medicine, University of São Paulo, São Paulo, Brazil

Esper Georges Kallas

Affiliation: Division of Clinical Immunology and Allergy, School of Medicine, University of São Paulo, São Paulo, Brazil

Introduction

Since the beginning of the HIV/AIDS pandemic, until the end of 2010, over 3 million children under 15 years of age have been infected with HIV-1 and 390,000 new viral infections each year (most recent data from UNAIDS/WHO; http://www.who.int/hiv/pub/progress_report2011/hiv_full_report_2011) have been reported. Up until 2009, the Centers for Disease Control and Prevention (CDC) estimates that perinatal transmission of the infection by the mother accounts for 91% of all AIDS cases among children under the age of 13 (CDC- Basic Statistics. Available at: http://www.cdc.gov/hiv/topics/surveillance/basic.htm. Accessed November 21, 2011). While mother to child transmission (MTCT) has been drastically reduced (1–2%) in rich countries, pregnant women living with HIV in poorer countries still have limited access to the same quality of counseling and antiretroviral therapy (ART) [1]. Therefore, these women are at higher risk of transmitting the virus to their offspring during pregnancy, labor or after childbirth via breastfeeding. Risk factors associated with MTCT include lack of receipt of prenatal ART, advanced maternal clinical status, detectable maternal viral load at delivery, low maternal CD4 T cell counts, immunogenetic host factors, and a high viral heterogeneity in the mother [2], [3], [4], [5], [6]. Without treatment, the chance of transmitting HIV from a mother to a baby is somewhere between 12% and 25% in resource rich settings, and between 20% and 45% in resource poor settings [7]. Although MTCT is being addressed by interventions of highly active ART (HAART), which usually comprises three drugs, the ultimate solution to HIV/AIDS will be a globally effective vaccine to curb HIV from spreading further. However, the development of such vaccines requires an in-depth knowledge of the virus strains being transmitted in the target population.

One of the most prominent features of HIV-1 is the remarkable accumulation of genetic diversity in its population during the course of infection. This diversity reflects the high mutation rate of reverse transcriptase (3×10−5 substitutions per site per generation) [8], rapid viral turnover (10−8 to 10−9 virions per day) [9], large number of infected cells (107 to 108 cells) [10], and recombination [11]. Consequently, the HIV-1 population is composed of a swarm of highly genetically related variants, i.e. a quasispecies, capable of rapidly adapting to various selective pressures. This diversity has been shown to have an impact not only on viral phenotypes at the level of transmission patterns, pathogenicity and immunology but also in responses to ART and vaccines [12], [13], [14]. Nine distinct genetic subtypes, (A–D, F–H, J and K) are joined in the pandemic today by more than 50 major circulating recombinant forms (CRFs) [http://www.hiv.lanl.gov/content/sequence/HIV/CRFs/CRFs.html] and numerous unique recombinant forms (URFs) have been isolated from individual patients [15], [16].

Brazil, the most populous country in the Latin America, is home to about one third of the people living with HIV (608,230) in Central and South America (UNAIDS. 2010–2011 Report on the Brazilian response to HIV/AIDS; http://www.unaids.org/en/dataanalysis/knowyourresponse/countryprogressreports/2012countries/UNGASS_2012_ingles_rev_08jun.pdf). According to 2008 estimates from UNAIDS (most recent data), 13,728 cases of HIV-infected Brazilian children were notified and 84.5% of them were infected by vertical transmission. Despite the availability of free HIV diagnostic tests and ART, perinatal transmission of HIV-1 infection by the mother remains a national public health challenge in some areas, where difficulties in providing the requisite antenatal HIV screening exist. Like in other European countries and in North America, HIV-1 subtype B is a major genetic clade circulating in the country but the overall prevalence of non-B strains, particularly URF BF1, C and URF BC, is increasing [17], [18], [19], [20]. Data from recent studies of the viral near full length genomes (NFLG) have provided evidence of Brazilian CRF strains designated as CRF28_BF, CRF29_BF, CRF39_BF, CRF40_BF, CRF46_BF and CRF31_BC [18], [19], [21], [22] (http://www.hiv.lanl.gov/content/sequence/HIV/CRFs/CRFs.html.). However, most HIV molecular variability studies in Brazil so far have been limited to adult patients and/or are based on sequence information derived from small genetic fragments that may not predict the subtype classification of other regions of the viral genomes.

Unlike HIV-1 infection in adults, multiple HIV-1 strains are not expected to co-circulate in children and the precise timing of the transmission event can be traced to the time of birth, providing a unique opportunity to explore the viral evolution [23].

The present study was undertaken to describe the genetic variability and identify resistance-associated mutations in HIV-1 isolates recovered from children, adolescents, and young adults with majority perinatally acquired HIV infection.

Materials and Methods

Study Population

A cross sectional study was conducted among 51 HIV-1 infected children (ages 0–14 years) and adolescents (aged 15–20 years) on ART followed between September 1992 and July 2009 at the Division of Pediatric Infectious Disease Clinic (CEADIPe), at the Federal University of São Paulo (UNIFESP), Brazil. All, but three of the participants had been born to HIV-1 seropositive mothers. From their medical records, almost all of the seropositive mothers had multiple sexual partners or had an exclusive sexual relationship with a partner whom they knew either to have had other sexual partners or, less commonly, to have a history of injection drug use. Three participants were born to women seronegative for HIV-1 and thus had an unknown transmission mode of their infection. All study participants gave written informed consent. Parents or legal guardians provided written informed consent on behalf of the children. The study plan and consent procedures were approved by the ethics committee of the federal University of São Paulo.

Amplification & Sequencing

Proviral DNA and RNA were extracted from peripheral blood mononuclear cells (PBMC) and plasma with commercial kits (QIAamp DNA Blood mini Kit and QIAamp Viral RNA Kit, QIAGEN, Germany) according to the manufacturer’s instructions. To make complementary DNAs, the extracted RNA samples were subjected to reverse transcription PCR using SuperScript III (Invitrogen, Carlsbad, CA). Both cDNAs and proviral DNAs were used as the PCR template, as this allowed amplification of the NFLGs from five overlapping fragments as previously described [24], [25], [26]. Amplification reactions were done in duplicate to eliminate PCR artifacts, ensuring that sequenced full-length genomes were not assembled from heterogeneous DNA targets. The expected sizes of the amplified products were verified using ethidium-bromide staining after agarose gel electrophoresis. Each PCR included a known HIV-1 subtype B positive control and an interspersed no DNA template negative controls. Strict laboratory precautions were taken to avoid cross contamination.

Both DNA complementary strands were sequenced directly from purified PCR products in an overlapping fragment of 400 nucleotides by using a variety of sequence specific primers, fluorescent-dye terminators, and Taq polymerase on an automated sequencer (ABI 3100, Applied Biosystems Inc., Foster City, CA). The data from the sequenced fragments were edited, assembled into contiguous sequences on a minimum overlap of 20 bp with a 85%–90% minimal mismatch and a consensus of both strands was formed by the Sequencher program (Gene Code Corp., Ann Arbor, MI). Such assembly criteria would prevent any fragment from overlapping if it is not derived from the same variant. All the sequences were checked for contamination by BLAST search against HIV-1 sequence database and among themselves [27].

Phylogenetic Analysis

Full genome sequences were aligned with reference sequences representing subtypes A–D, F–H, J and K obtained from the Los Alamos database (http://hiv-web.lanl.gov) using the CLUSTAL X program [28] with the “slow-accurate” default alignment parameters and IUB DNA weight matrix. Aligned sequences were manually edited and trimmed to the minimal shared length in the BioEdit Sequence Alignment Editor Program. The gap-stripped aligned sequences were screened for the presence of recombination patterns by the bootscan methods implemented in the SIMPLOT program v3.2 beta 29,30 and the jumping profile Hidden Markov Model [31]. For the bootscan method, nucleotide distances were calculated in a sliding window of 300 bp moving in steps of 30 bp by the F84 model of evolution, transition\transversion ratio of 2.0. Recombinant regions of the alignment as determined by the crossover points from bootscanning were analyzed separately by phylogenetic analysis. Maximum likelihood phylogenies were constructed using the GTR+I+G substitution model and a BIONJ starting tree. Heuristic tree searches under the ML optimality criterion were performed using the NNI branch-swapping algorithm. The approximate likelihood ratio test (aLRT) based on a Shimodaira-Hasegawa-like procedure was used as a statistical test to calculate branch support. The maximum composite likelihood in MEGA version 4 [32] was used to calculate the genetic distances between and within isolates.

Genotyping Analysis

All amino acid positions associated with ART resistance, in the protease (Pro) and reverse transcription (RT) regions, according to IAS-USA 2011 and Stanford HIV drug resistance database were evaluated on both plasma and blood samples.

Measurement of HIV RNA and Cell Count

The viral load was measured using the Roche Amplicor HIV-1 Monitor test (Roche, Branchburg, NJ; lower limit of detection 50 copies per ml). CD4+ and CD8+ T cell counts were performed using a lymphocyte staining panel containing CD3, CD4, CD8, and CD45 conjugated monoclonal antibodies (BD Biosciences, San Diego, California, USA).

Genotypic Tropism Analysis

For the predictions of HIV tropism, the env region identified in the NFLGs and partial-length env sequences that would encompass the V3 region were analyzed using a tropism prediction algorithm implemented as the web-based service geno2pheno [coreceptor] http://www.geno2pheno.org. To minimize the number of false predictions of CXC chemokine receptor 4 (CXCR4 or X4) tropic sequences as C–C chemokine receptor 5 (CCR5 or R5) tropic, a conservative false-positive rate (FPR) of 20% was used as a cutoff. Therefore, X4 or X4 dual/mixed-tropic viruses (X4/DM) were reported positive if their sequences had a prediction result FPR of ≤20% or the 11/25 rule predicted a X4 virus, otherwise, they were considered R5-tropic viruses.

All nucleotide sequences obtained during our study were reported to GenBank (Accession numbers pending).

Results

Samples

In total, 49 paired samples of whole blood and plasma and 2 unpaired samples of each type from an additional 4 different patients were subjected to NFLG amplification and sequencing. Of these 51 subjects, 21 (41.2%) were males and 30 (58.8%) were females. The participants were predominantly white (69.6%), were followed from birth and reached a median age of 11.5 years (range between 4 and 20 years). The median HIV-1 viral load and CD4 cell count, as judged by levels at the time of inclusion, were 6.34×102 copies/ml (range,<49−7.5×105) and 640 cells/mm3 (range, 18–1821 cells/mm3), respectively. Four (8%) patients were naïve and an additional 2 (4%) patients were not taking therapy at the time of enrollment, although they were drug-experienced. The therapeutic status was not known for one subject. Among the 43 subjects, 69.7% had received ≥3 different ART regimens during their follow-up. The median duration of ART at the time of genotyping was 27.5 months. The main characteristics and ART regimens of the study population are given in Table 1.

[Figure omitted. See PDF.]

Table 1. Patient characteristics for the study population.

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

NFLG and Partial Amplification of HIV-1 from both PBMC and Plasma Specimens

Sequences were obtained for all five overlapped fragments that cover the NFLGs of 4 PBMC DNA and one plasma RNA virus. Partial sequences were obtained from at least one fragment derived from 38 blood and 24 plasma samples as shown in Table 2. Of the 26 plasma samples for which partial and NFLGs failed, 21 (80.7%) had a viral load under 500 copies/ml and the remaining 5 (19.3%) RNA viruses had multiple peaks present in the sequencing chromatogram probably indicating different quasispecies in the same sample or HIV-1 dual infections. On the other hand, partial amplification of 1126 bp of fragment B1(Nucleotide position from start of HXB2 genome 2196–3322) and 494 bp stretch of fragment D (Nucleotide position from start of HXB2 genome; 8997–9491) were subtype B positive for isolates 010BR_IMT_010 and 010BR_IMT_051, respectively, and both patients had viral loads below 50 copies/ml (Table 2). These results may suggest an underestimation of the measured viremia or high efficacy of our nested PCR approach in some patients. On the other hand, our results among paired samples demonstrated that 20 patients had detectable HIV proviral DNA and undetectable viral RNA, 21 were dually positive for viral RNA and DNA, and 4 patients were dually negative. It is unclear why we were unable to amplify more plasma RNA viruses, particularly for patients 010_BR_IMT_05, 010_BR_IMT_12, 010_BR_IMT_54, and 010_BR_IMT_58 (median viral load 1.5×104, range 1.6×103–2.9×104) using our fragment-based amplification strategy, although RNA degradation may account for this finding.

[Figure omitted. See PDF.]

Table 2. The near full-length genomic (NFLG) and partial fragments subtyping of HIV-from plasma and blood samples.

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

HIV Variants and Sequences

Based on phylogenetic analysis, the NFLGs and partial proviral nucleotide sequences (n = 42) of the clinical HIV-1 isolates indicated that 22 (52.4%) patients were infected with HIV-1 subtype B, 16 (38.1%) were infected with a mosaic consisting of subtype BF1 and 4 (9.5%) were infected with sub-subtype F1 (Table 2). Of the total 25 plasma samples for which viral subtype was determined, 17 (68%) were classified as subtype B, 3 (12%) were sub-subtype F1, and 5 (20%) were BF1 recombinant viruses (Table 2). All chimeric viruses were unique according to their recombination profile, i.e., not assigned to any subtype or CRF (Figure 1). The relationships of the viral sequences from patients’ PBMCs to the sequences obtained from the corresponding RNA virus within the same regions were examined for each patient to assess the viral diversity in both compartments. The results revealed that all but one patient, 010BR_IMT_020, had plasma RNA and proviral DNA variation only ranging between 0–2.7% (Figure 2). These relations were further confirmed by phylogenetic analysis, which showed close branching as demonstrated in Figure 3. These findings may indicate that the primary infected PBMCs of these patients were likely the source of plasma circulating viral sequences however; more sophisticated genetic tests able to detect viral population structure are needed to confirm this conclusion. The observed differences in the percent nucleotide variations between proviruses and plasma free viruses in this group may reflect evolution that occurs during the initial phase of acute infection, before the therapeutic control of HIV-1 replication is established. Surprisingly, the intra-individual plasma and proviral sequence variation for patient 010BR_IMT_020 in the overlapped regions depicted in Figure 4 were 9.8% and 6.5%, respectively, indicating that the plasma viruses were derived from a population significantly distinct from those of the cellular sources in this 13 years old asymptomatic patient. This result is consistent with dual distinct variants of the same subtype being involved in establishing infection. Dual infection with subclade F1 and BF1 recombinant was observed in patient 010BR_IMT_041 plasma sample (Figure 5). This patient was a nine year old child who diagnosed in February 2005 and until the sampling period had been asymptomatic. The patient had been receiving ART since September 2005. This observation of dual infection occurred accidentally during assembling of the generated data, in which some sequences failed to assemble to other overlapping stretches of fragment B1. As a result, we sought to compare this stretch to HIV sequences available from public databases. Upon analysis with the basic local alignment search tool (BLAST) available from GenBank, the stretch (010BR_IMT_041_PL- REC; 548 bp) from plasma revealed high percentages of nucleotide sequence identity to the BF1 isolate 99JY-TRA0133 (Genbank accession: JN235964), whereas the other larger fragment (010BR_IMT_041_pl; 5720 bp) revealed high homology to subclade F1 isolate 02BR082 (Genbank accession: FJ771006) at the nucleotide levels. To ensure that the generation of the two consensus sequences from patient 010BR_IMT_041 plasma sample was not the result of sample contamination, repeat sequence analysis using the purified B1 amplicon was performed and revealed identical findings. These results possibly indicate that some internal sequencing primers of fragment B1 preferentially annealed to the BF1 string during sequencing reaction. Regions that were the same F1 subclade in the two pols were then compared to determine whether the 010BR_IMT_041_PL viruses were the actual parents of the recombinant fragment or if an infection in this patient was acquired with two genetically distinct viruses (Figure 5B). While both partial pol genes were sub-subtype F1 fragments, these were from different subclade F1 isolates, since the sequences from the two plasma demonstrated high nucleotide divergence (up to 6.8%). Moreover, as shown in Figure 5B, both F1 non-recombinant sequences recovered from plasma and PBMC clustered separately (aLRT 100%) and the branch lengths separating them from the F1 fragment involved in the recombination event were typical for other sequences of unrelated F1 variants. The analysis was then extended to include isolates with non-overlapping fragments, namely 010BR_IMT_013 and 010BR_IMT_027, to determine whether the PBMC viruses were truly parental strains to those recovered from the plasma. For this purpose, the phylogenetic clustering profile of the non-overlapped fragments from both compartments were compared to a number of additional Brazilin subtype B and other HIV-1 reference sequences to increase our confidence in the analyses and provide a broader perspective. These results revealed the magnitude of aLRT value supporting the identical clustering of the plasma isolate 010BR_IMT_013_pl (Figure 6A) and the proviral 010BR_IMT_013_pr strains (Figure 6B) with the Brazilian subtype B BREPM 1040 and 05BR 1092 subtype B sequences (branch marked with green color). Based on these results, it is possible to assume that the primary infected PBMCs in this patient were likely the source of the plasma circulating viral sequences. However, this interpretation does not hold true when the analysis was applied to the plasma and proviral non-overlapping fragments of patient 010BR_IMT_027. The clustering profile to subtype B references and genetic distances as shown in Figure S1 were significantly different between both fragments in this patient indicating dual infection with distinct subtype B variants.

[Figure omitted. See PDF.]

Figure 1. Schematic representation of the NFLG, partial structure and breakpoint profiles of the BF1 sequences identified in this study both from HIV RNA and proviral DNA.

Samples that were identified in this study to host distinct viruses are indicated with the star symbol. The region of subclade F1 and subtypes B are indicated at the bottom.

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

[Figure omitted. See PDF.]

Figure 2. Genetic distances of overlapping regions between isolates recovered from patients with paired plasma and PBMC samples.

Concatenated sequences are indicated with the star symbol. The region of subclade F1 and subtypes B are indicated at the bottom.

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

[Figure omitted. See PDF.]

Figure 3. Maximum-likelihood phylogenetic trees form each non-recombinant fragment were constructed using all available sequences from proviral DNA (indicated by black circles) and plasma isolate (indicated by an empty circles) along with HIV-1 reference sequences from the Los Alamos HIV-1 database representing 11 genetic subtypes.

The numbering for the HIV-1 fragment A, B1, B2 and C sequences corresponds to the HXB2 reference sequence. For purposes of clarity, the tree was midpoint rooted. The approximate likelihood ratio test (aLRT) values of ≥70% are indicated at nodes. The scale bar represents 0.05 nucleotide substitutions per site.

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

[Figure omitted. See PDF.]

Figure 4. Genetic distances and phylogenetic tree constructed using a maximum-likelihood method from concatenated regions of HIV RNA (indicated by black circles) and proviral DNA (indicated by black squares) marked with Arabic numbers 1 and 2 from patient 010BR_IMT_020 along with HIV-1 reference sequences from the Los Alamos HIV-1 database representing 11 genetic subtypes.

For purposes of clarity, the tree was midpoint rooted. The approximate likelihood ratio test (aLRT) values of ≥70% are indicated at nodes. The scale bar represents 0.05 nucleotide substitutions per site.

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

[Figure omitted. See PDF.]

Figure 5. Genetic distances and phylogenetic tree constructed using a maximum-likelihood method of HIV isolates recovered from patient 01BR_IMT_041.

A: Phylogenetic tree from concatenated regions assigned as subtype B from the BF1 recombinant isolate. B: Phylogenetic tree showing the clustering pattern of F1 sequences (marked by dotted box). F1 region from genuine F1 sequence recovered from plasma and PBMCs are marked by a black circle and square, respectively while the F1 region from the BF1 recombinant sequence is marked by an empty square. The approximate likelihood ratio test (aLRT) values of ≥70% are indicated at nodes. The scale bar represents 0.05 nucleotide substitutions per site.

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

[Figure omitted. See PDF.]

Figure 6. Phylogenetic clustering profile of the non-overlapped fragments assigned as subtype B from both plasma and provirus isolate 01BR_IMT_013 were compared to a number of additional Brazilian subtype B sequences and other HIV-1 reference sequences from the Los Alamos HIV-1 database representing 11 genetic subtypes.

For purposes of clarity, the tree was midpoint rooted. The approximate likelihood ratio test (aLRT) values of ≥70% are indicated at nodes. The scale bar represents 0.05 nucleotide substitutions per site.

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

Genotypic Drug Resistance Test Results

The results of the genotyping analysis from both plasma RNA and whole blood DNA are presented in Table 3 and 4. Resistance analyses were performed from proviral DNA in 32 patients (19 PR/RT and 13 PR) and from plasma-associated RNA in 21 patients (8 PR/RT and 13 PR). Regardless of paired or unpaired samples, 10 of the 31 (32.2%) and 12 of the 21 (57.1%) subjects with recovered proviral and plasma PR sequences, respectively, were on ART at the time of specimen collection and were infected with major mutants resistant to protease inhibitors (PIs). Regarding the primary resistant mutations for PIs among the naive patients with available PR sequences from PBMCs (n = 4) and plasma (n = 3), various mutations were detected in only one patient (010BR_IMT_034; Table 3 and 4). The RT region of the provirus of the same patient displayed some major transmitted mutations both for the NRTIs and NNRTIs. Detailed frequency of single Pro and RT mutations detected in patients on ART both in plasma RNA and whole blood DNA is also illustrated in Table 3 and 4.

[Figure omitted. See PDF.]

Table 3. Drug-resistance mutations detected in plasma.

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

[Figure omitted. See PDF.]

Table 4. Drug-resistance mutations detected in PBMC.

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

V3 Sequence Analysis and Viral Tropism

An evaluation of the V3 loop amino acids and predictions of viral tropism were performed for patients with available sequences from PBMCs (n = 10) and plasma (n = 7) of the derived fragment C intact frame sequences. The inferred HIV tropism in paired samples of plasma and PBMCs was successfully determined in 3 samples and all concordant with the R5 virus. The inferred HIV tropism study in the other 3 plasma demonstrated that 2 patients harbored the R5 virus. The V3 sequences of the 7 patients with available sequences from only PBMC were predicted to be R5-tropic virus except for patient 01BR_IMT_035 who harbored an X4 strain.

Discussion

This study describes the genetic variability and the prevalence of drug resistance mutations and co-receptor usage of HIV-1 variants in a small, well sampled group of children and adolescents. The majority of these patients acquired their infection through vertical transmission during the period 1992–2007. The results presented confirmed that subtype B is still the main HIV-1 variant and concordant with data from other studies on adult and children populations from Brazil [24], [33], [34], [35]. The most remarkable observations in this study are that at least 38.1% of the 42 patients with proviral DNA sequences are infected with HIV-1 BF1 recombinant variants, which is relatively much higher if compared to earlier studies on children and adolescent patients in Brazil [33], [34], [35], [36], [37], [38]. This difference is not surprising, because small fragments from different regions of HIV genomes were characterized in the previous studies while we used larger overlapped fragments to sequence the full-length genome, which undoubtedly provides efficient discrimination of HIV subtypes and the recombinant forms. Thus, the earlier study is likely to have missed some recombinants. Despite the high rate of recombination in our study, it is probable that our results have also underestimated the true rate of infection with BF1 recombinant viruses, particularly among patients with partially sequenced viral fragments. Thus, it is possible that the BF1 infection in this group may be higher than what was observed if we had sequenced the virus NFLG in all samples. Our attempts to amplify the NFLG or additional larger fragments for some samples to determine if recombination had occurred were unsuccessful. Other likely explanations for underestimation of BF1 recombination rate is that some isolates could have been undetected by our PCR strategy because of a mismatch at the primer binding sites, low proviral load, employment of consensus sequences or that the BF1 isolates were maintained in another reservoir other than the CD4-positive compartment that was sampled in the peripheral blood. The results that indicate none of the BF1 recombinant structures identified in this study showed any similarity to the known CRFs or other recombinants strongly suggests that new recombinants are arising continually in São Paulo, Brazil.

Additional observations of this study are the description of the high level of intra-host diversity with evidence of mixed infections with the same or distinct HIV-1 subtypes. The observation that patients may be simultaneously infected with different HIV-1 subtypes has been reported in numerous cases and considered of significant interest. For instance, the first documented dual infection of two distinct HIV-1 subtypes B and E (later designated as CRF01_AE) was reported in Thailand [39]. Janini et al. [40] reported the first case of both horizontal and subsequent vertical transmission of 2 distinct HIV-1 subtypes from 1 dually infected person to another. In the present study, dual infections were evident in three patients, lending further support to previous studies [41], [42], [43], as this event is far more common in Brazil where both subtypes co-circulate. The fact that existence of dual infection in some patients contrasts with the hypothesis that an initial viral infection produces some degree of protection against a second infecting HIV subtype. If we assume that super-infection occurs in these cases, then it is conceivable that antiviral immunity evoked by one subtype had insufficient broad protection at the time of primary infection against a second infecting virus. Indeed this assumption has been challenged previously by convincing findings revealing that a second super-infection with a different HIV-1 strain can occur long after an initial infection is established and can hasten the disease process [44], [45]. Alternatively, the subjects 010BR_IMT_041 and 010BR_IMT_027 (confirmed MTCT) may have been vertically and concomitantly infected with different HIV strains at the same time. On the basis of this assumption, our results may suggest that despite the genetic bottleneck occurring upon vertical transmission of HIV-1, the replication capacity of transmitted variants is not necessarily reduced. This interpretation is in line with previous studies that provided evidences of multiple-variant transmission in MTCT, and also agreed with the conclusion that in a majority of cases the infant is infected with a single isolate [46], [47], [48]. By the lack of mode of transmission in patient 010BR_IMT_020, it was therefore not possible to interpret the simultaneous detection of both viruses.

Our results on genotype resistance mutations are consistent with previous studies with similar subjects indicating that prevalence of major mutations conferring ART resistance in viral DNA/RNA of such chronically infected groups is common (20/38, 52.6%) [34]. In patients with available sequences from the Pro and/or RT, the mutations found in PBMCs were generally also found in the plasma, although some of the patients showed few differences between the two compartments, while in one patient (010BR_IMT_011) the 69 insertion in the protease region was found in PBMC, but not in plasma.

Regardless of the sample compartment, the analysis of HIV tropism revealed two patients with X4 viruses and both with CDC class “C3” reflecting advanced disease. The assessment of HIV tropism in our study was limited to sequence- based algorithms rather than using phenotypic methods. Although phenotypic assays still have an edge over genotypic methods, genotypic predictors prove to be highly concordant with phenotype data and can reliably be used to determine viral tropism with better results in PBMC than in plasma samples [49]. In this study, we used geno2pheno, which has shown a similar performance to the Trofile phenotypic assay, the most often used tropism method [50]. Moreover, the method has been shown to achieve higher sensitivity while retaining high level of specificity when compared with the performance of different algorithms [51], [52].

We are aware that the demonstration of the high recombination rate and evidence of double infections and their association with virological response and viral tropism should be based on a larger dataset to establish statistical influence of these factors in determining the outcome. Another limitation of this study is that direct bulk sequencing and genotyping of HIV-1 in plasma and whole blood might underestimate low-level minority species present as quasispecies which could be evidenced by more complex methods, such as massive parallel pyrosequencing [53], [54], [55]. Despite these limitations, the results of this analysis indicate that HIV-1 recombination and dual infections are much more frequent than thought previously among children and adolescents in this region. Evidently, more extensive studies with large sample sizes are required to unravel the mechanisms underlying the emergence of these recombinants and their implications for HIV control.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Comparison of phylogenetic clustering profile of the fragments assigned as subtype B from both plasma and provirus isolate 01BR_IMT_027 were compared to a number of additional Brazilian subtype B sequences and other HIV-1 reference sequences from the Los Alamos HIV-1 database representing 11 genetic subtypes. For purposes of clarity, the tree was midpoint rooted. The approximate likelihood ratio test (aLRT) values of ≥70% are indicated at nodes. The scale bar represents 0.05 nucleotide substitutions per site.

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

(TIFF)

Author Contributions

Attending physician for sample collection: RCMS. Conceived and designed the experiments: SSS RP. Performed the experiments: SSS RP ACSO VPM. Analyzed the data: SSS RP ACSO. Contributed reagents/materials/analysis tools: MTMG KC CST DFN ECS EGK. Wrote the paper: SSS.

Citation: Sanabani SS, Pessôa R, Soares de Oliveira AC, Martinez VP, Giret MTM, de Menezes Succi RC, et al. (2013) Variability of HIV-1 Genomes among Children and Adolescents from São Paulo, Brazil. PLoS ONE 8(5): e62552. https://doi.org/10.1371/journal.pone.0062552

References

1. Templeton AR, Kramer MG, Jarvis J, Kowalski J, Gange S, et al. (2009) Multiple-infection and recombination in HIV-1 within a longitudinal cohort of women. Retrovirology 6: 54.

2. Bitnun A, King SM, Arneson C, Read SE (2002) Failure to prevent perinatal HIV infection. CMAJ 166: 904–905.

3. Peters V, Liu KL, Dominguez K, Frederick T, Melville S, et al. (2003) Missed opportunities for perinatal HIV prevention among HIV-exposed infants born 1996–2000, pediatric spectrum of HIV disease cohort. Pediatrics 111: 1186–1191.

4. Sperling RS, Shapiro DE, Coombs RW, Todd JA, Herman SA, et al. (1996) Maternal viral load, zidovudine treatment, and the risk of transmission of human immunodeficiency virus type 1 from mother to infant. Pediatric AIDS Clinical Trials Group Protocol 076 Study Group. N Engl J Med 335: 1621–1629.

5. Pitt J, Brambilla D, Reichelderfer P, Landay A, McIntosh K, et al. (1997) Maternal immunologic and virologic risk factors for infant human immunodeficiency virus type 1 infection: findings from the Women and Infants Transmission Study. J Infect Dis 175: 567–575.

6. Wolinsky SM, Wike CM, Korber BT, Hutto C, Parks WP, et al. (1992) Selective transmission of human immunodeficiency virus type-1 variants from mothers to infants. Science 255: 1134–1137.

7. Wiktor SZ, Ekpini E, Nduati RW (1997) Prevention of mother-to-child transmission of HIV-1 in Africa. AIDS 11 Suppl B: S79–87.

8. Mansky LM, Temin HM (1995) Lower in vivo mutation rate of human immunodeficiency virus type 1 than that predicted from the fidelity of purified reverse transcriptase. J Virol 69: 5087–5094.

9. Wei X, Ghosh SK, Taylor ME, Johnson VA, Emini EA, et al. (1995) Viral dynamics in human immunodeficiency virus type 1 infection. Nature 373: 117–122.

10. Chun TW, Carruth L, Finzi D, Shen X, DiGiuseppe JA, et al. (1997) Quantification of latent tissue reservoirs and total body viral load in HIV-1 infection. Nature 387: 183–188.

11. Sabino EC, Shpaer EG, Morgado MG, Korber BT, Diaz RS, et al. (1994) Identification of human immunodeficiency virus type 1 envelope genes recombinant between subtypes B and F in two epidemiologically linked individuals from Brazil. J Virol 68: 6340–6346.

12. Spira S, Wainberg MA, Loemba H, Turner D, Brenner BG (2003) Impact of clade diversity on HIV-1 virulence, antiretroviral drug sensitivity and drug resistance. J Antimicrob Chemother 51: 229–240.

13. McBurney SP, Ross TM (2008) Viral sequence diversity: challenges for AIDS vaccine designs. Expert Rev Vaccines 7: 1405–1417.

14. Rambaut A, Posada D, Crandall KA, Holmes EC (2004) The causes and consequences of HIV evolution. Nat Rev Genet 5: 52–61.

15. Peeters M, Toure-Kane C, Nkengasong JN (2003) Genetic diversity of HIV in Africa: impact on diagnosis, treatment, vaccine development and trials. AIDS 17: 2547–2560.

16. Thomson MM, Najera R (2005) Molecular epidemiology of HIV-1 variants in the global AIDS pandemic: an update. AIDS Rev 7: 210–224.

17. Barreto CC, Nishyia A, Araujo LV, Ferreira JE, Busch MP, et al. (2006) Trends in antiretroviral drug resistance and clade distributions among HIV-1–infected blood donors in Sao Paulo, Brazil. J Acquir Immune Defic Syndr 41: 338–341.

18. Sanabani S, Neto WK, de Sa Filho DJ, Diaz RS, Munerato P, et al. (2006) Full-length genome analysis of human immunodeficiency virus type 1 subtype C in Brazil. AIDS Res Hum Retroviruses 22: 171–176.

19. Sanabani S, Kleine Neto W, Kalmar EM, Diaz RS, Janini LM, et al. (2006) Analysis of the near full length genomes of HIV-1 subtypes B, F and BF recombinant from a cohort of 14 patients in Sao Paulo, Brazil. Infect Genet Evol 6: 368–377.

20. Passaes CP, Guimaraes ML, Bello G, Morgado MG (2009) Near full-length genome characterization of HIV type 1 unique BC recombinant forms from Southern Brazil. AIDS Res Hum Retroviruses 25: 1339–1344.

21. Guimaraes ML, Eyer-Silva WA, Couto-Fernandez JC, Morgado MG (2008) Identification of two new CRF_BF in Rio de Janeiro State, Brazil. AIDS 22: 433–435.

22. Santos AF, Sousa TM, Soares EA, Sanabani S, Martinez AM, et al. (2006) Characterization of a new circulating recombinant form comprising HIV-1 subtypes C and B in southern Brazil. AIDS 20: 2011–2019.

23. Aulicino PC, Bello G, Guimaraes ML, Ruchansky D, Rocco C, et al. (2011) Longitudinal analysis of HIV-1 BF1 recombinant strains in vertically infected children from Argentina reveals a decrease in CRF12_BF pol gene mosaic patterns and high diversity of BF unique recombinant forms. Infect Genet Evol 11: 349–357.

24. Sanabani SS, Pastena ER, da Costa AC, Martinez VP, Kleine-Neto W, et al. (2011) Characterization of partial and near full-length genomes of HIV-1 strains sampled from recently infected individuals in Sao Paulo, Brazil. PLoS One 6: e25869.

25. Sanabani SS, Pastena ER, Neto WK, Martinez VP, Sabino EC (2010) Characterization and frequency of a newly identified HIV-1 BF1 intersubtype circulating recombinant form in Sao Paulo, Brazil. Virol J 7: 74.

26. Sanabani SS, Pastena ER, Kleine Neto W, Barreto CC, Ferrari KT, et al. (2009) Near full-length genome analysis of low prevalent human immunodeficiency virus type 1 subclade F1 in Sao Paulo, Brazil. Virol J 6: 78.

27. Korber B, Gaschen B, Yusim K, Thakallapally R, Kesmir C, et al. (2001) Evolutionary and immunological implications of contemporary HIV-1 variation. Br Med Bull 58: 19–42.

28. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25: 4876–4882.

29. Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, et al. (1999) Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol 73: 152–160.

30. Salminen MO, Carr JK, Burke DS, McCutchan FE (1995) Identification of breakpoints in intergenotypic recombinants of HIV type 1 by bootscanning. AIDS Res Hum Retroviruses 11: 1423–1425.

31. Schultz AK, Zhang M, Leitner T, Kuiken C, Korber B, et al. (2006) A jumping profile Hidden Markov Model and applications to recombination sites in HIV and HCV genomes. BMC Bioinformatics 7: 265.

32. Kumar S, Tamura K, Nei M (2004) MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform 5: 150–163.

33. de Oliveira CM, Almeida FJ, Rodrigues R, Crozatti M, Vazquez CM, et al. (2008) High frequency of BF mosaic genomes among HIV-1-infected children from Sao Paulo, Brazil. Arch Virol 153: 1799–1806.

34. Almeida FJ, Berezin EN, Rodrigues R, Safadi MA, Arnoni MV, et al. (2009) Diversity and prevalence of antiretroviral genotypic resistance mutations among HIV-1-infected children. J Pediatr (Rio J) 85: 104–109.

35. da Silveira AA, Cardoso LP, Francisco RB, de Araujo Stefani MM (2012) HIV type 1 molecular epidemiology in pol and gp41 genes among naive patients from Mato Grosso do Sul State, central western Brazil. AIDS Res Hum Retroviruses 28: 304–307.

36. Alcantara KC, Lins JB, Albuquerque M, Aires LM, Cardoso LP, et al. (2012) HIV-1 mother-to-child transmission and drug resistance among Brazilian pregnant women with high access to diagnosis and prophylactic measures. J Clin Virol 54: 15–20.

37. Brindeiro PA, Brindeiro RM, Mortensen C, Hertogs K, De Vroey V, et al. (2002) Testing genotypic and phenotypic resistance in human immunodeficiency virus type 1 isolates of clade B and other clades from children failing antiretroviral therapy. J Clin Microbiol 40: 4512–4519.

38. Machado ES, Lambert JS, Watson DC, Afonso AO, da Cunha SM, et al. (2004) Genotypic resistance and HIV-1 subtype in Brazilian children on dual and triple combination therapy. J Clin Virol 30: 24–31.

39. Artenstein AW, VanCott TC, Mascola JR, Carr JK, Hegerich PA, et al. (1995) Dual infection with human immunodeficiency virus type 1 of distinct envelope subtypes in humans. J Infect Dis 171: 805–810.

40. Janini LM, Tanuri A, Schechter M, Peralta JM, Vicente AC, et al. (1998) Horizontal and vertical transmission of human immunodeficiency virus type 1 dual infections caused by viruses of subtypes B and C. J Infect Dis. 177: 227–231.

41. Soares de Oliveira AC, Pessoa de Farias R, da Costa AC, Sauer MM, Bassichetto KC, et al. (2012) Frequency of subtype B and F1 dual infection in HIV-1 positive, Brazilian men who have sex with men. Virol J 9: 223.

42. Ramos A, Tanuri A, Schechter M, Rayfield MA, Hu DJ, et al. (1999) Dual and recombinant infections: an integral part of the HIV-1 epidemic in Brazil. Emerg Infect Dis 5: 65–74.

43. Diaz RS, Sabino EC, Mayer A, Mosley JW, Busch MP (1995) Dual human immunodeficiency virus type 1 infection and recombination in a dually exposed transfusion recipient. The Transfusion Safety Study Group. J Virol 69: 3273–3281.

44. Altfeld M, Allen TM, Yu XG, Johnston MN, Agrawal D, et al. (2002) HIV-1 superinfection despite broad CD8+ T-cell responses containing replication of the primary virus. Nature 420: 434–439.

45. Koelsch KK, Smith DM, Little SJ, Ignacio CC, Macaranas TR, et al. (2003) Clade B HIV-1 superinfection with wild-type virus after primary infection with drug-resistant clade B virus. AIDS 17: F11–16.

46. Samleerat T, Braibant M, Jourdain G, Moreau A, Ngo-Giang-Huong N, et al. (2008) Characteristics of HIV type 1 (HIV-1) glycoprotein 120 env sequences in mother-infant pairs infected with HIV-1 subtype CRF01_AE. J Infect Dis 198: 868–876.

47. Russell ES, Kwiek JJ, Keys J, Barton K, Mwapasa V, et al. (2011) The genetic bottleneck in vertical transmission of subtype C HIV-1 is not driven by selection of especially neutralization-resistant virus from the maternal viral population. J Virol 85: 8253–8262.

48. Verhofstede C, Demecheleer E, De Cabooter N, Gaillard P, Mwanyumba F, et al. (2003) Diversity of the human immunodeficiency virus type 1 (HIV-1) env sequence after vertical transmission in mother-child pairs infected with HIV-1 subtype A. J Virol. 77: 3050–3057.

49. Skrabal K, Low AJ, Dong W, Sing T, Cheung PK, et al. (2007) Determining human immunodeficiency virus coreceptor use in a clinical setting: degree of correlation between two phenotypic assays and a bioinformatic model. J Clin Microbiol 45: 279–284.

50. McGovern RA, Thielen A, Mo T, Dong W, Woods CK, et al. (2010) Population-based V3 genotypic tropism assay: a retrospective analysis using screening samples from the A4001029 and MOTIVATE studies. AIDS 24: 2517–2525.

51. Chueca N, Garrido C, Alvarez M, Poveda E, de Dios Luna J, et al. (2009) Improvement in the determination of HIV-1 tropism using the V3 gene sequence and a combination of bioinformatic tools. J Med Virol 81: 763–767.

52. Raymond S, Delobel P, Mavigner M, Cazabat M, Souyris C, et al. (2008) Correlation between genotypic predictions based on V3 sequences and phenotypic determination of HIV-1 tropism. AIDS 22: F11–16.

53. Bushman FD, Hoffmann C, Ronen K, Malani N, Minkah N, et al. (2008) Massively parallel pyrosequencing in HIV research. AIDS 22: 1411–1415.

54. Dudley DM, Chin EN, Bimber BN, Sanabani SS, Tarosso LF, et al. (2012) Low-cost ultra-wide genotyping using Roche/454 pyrosequencing for surveillance of HIV drug resistance. PLoS One 7: e36494.

55. Bimber BN, Dudley DM, Lauck M, Becker EA, Chin EN, et al. (2010) Whole-genome characterization of human and simian immunodeficiency virus intrahost diversity by ultradeep pyrosequencing. J Virol 84: 12087–12092.