Origin and evolutionary history of HIV-1 subtype C in Brazil

Soluble CD14 is subtype-dependent in serum but not in cerebrospinal fluid in people with HIV

Monocytes and macrophages activation are crucial in human immunodeficiency virus (HIV) central nervous system (CNS) infection and HIV associated neurocognitive disorders (HAND) pathogenesis. The soluble form of CD14 (sCD14) is a marker of monocyte activation. We hypothesized that sCD14 levels would be lower in people with HIV-1 subtype C (HIV-1C) than in HIV-1B owing to a variant Tat cysteine dimotif (C30S31) with reduced chemotactic activity. A total of 68 paired cerebrospinal fluid (CSF) and blood samples from people with HIV (PWH); 27 samples of the HIV-1B subtype and 40 of the non-B HIV-1 subtypes (including 26,HIV-1C), and 18 HIV-negative controls were included. sCD14 levels were quantified using a high-sensitivity enzyme-linked immunosorbent assay. sCD14 increase in serum, but not in CSF, was higher in samples from HIV-1B than HIV-1C (p = 0.002; Cohen's d, 0.7). CSF or serum sCD14 values were not correlated with global deficit score or specific cognitive domains. The impact of HIV-1 on monocyte stimulation biomarkers evaluated by sCD14 in serum was subtype-dependent, higher in HIV-1B than HIV-1C, consistent with reduced chemotactic activity as hypothesized.

Dynamic Dispersion of HIV-1 Subtype C Toward Brazilian Northeastern Region

The subtype C accounts for >50% of HIV type 1 (HIV-1) infections worldwide and it is currently the predominant viral form in South Brazil. Subtype C has been reported in all Brazilian regions; however, the phylogenetic relationship among strains circulating in those regions still remains unclear. This study aimed to investigate the origin and dynamic dispersion of HIV-1 subtype C toward Northeast Brazil. Our phylogenetic analysis suggests that most subtype C strains circulating in Brazil (99%) are descendant from the main lineage whose entrance in the country was previously described in the 1970s. According to the literature, additional introductions of subtype C were reported in the country through the Southeast region and in this study we identified another entry event that occurred most likely through the North region. Furthermore, our analysis suggests that the spread of subtype C to Brazilian Northeastern states occurred through multiple independent introductions of the main lineage that originated in South Brazil between mid-1980s and late 1990s. Despite the observation of eventual new HIV-1 subtype C introductions, our results highlight the predominance of a single lineage of this subtype in Brazil and the importance of South region in its dissemination throughout the country.

Evolutionary dynamics of HIV-1 subtype C in Brazil

The extensive genetic diversity of HIV-1 is a major challenge for the prevention and treatment of HIV-1 infections. Subtype C accounts for most of the HIV-1 infections in the world but has been mainly localized in Southern Africa, Ethiopia and India. For elusive reasons, South Brazil harbors the largest HIV-1 subtype C epidemic in the American continent that is elsewhere dominated by subtype B. To investigate this topic, we collected clinical data and viral sequences from 2611 treatment-naïve patients diagnosed with HIV-1 in Brazil. Molecular epidemiology analysis supported 35 well-delimited transmission clusters of subtype C highlighting transmission within South Brazil but also from the South to all other Brazilian regions and internationally. Individuals infected with subtype C had lower probability to be deficient in CD4⁺ T cells when compared to subtype B. The HIV-1 epidemics in the South was characterized by high female-to-male infection ratios and women-to-child transmission. Our results suggest that HIV-1 subtype C probably takes advantage of longer asymptomatic periods to maximize transmission and is unlikely to outcompete subtype B in settings where the infection of women is relatively less relevant. This study contributes to elucidate factors possibly underlying the geographical distribution and expansion patterns of the most spread HIV-1 subtypes.

HIV-1C and HIV-1B Tat protein polymorphism in Southern Brazil

The transactivator of transcription (Tat) is a key HIV regulatory protein. We aimed to identify the frequency of key polymorphisms in HIV-1C compared with HIV-1B Tat protein, chiefly in the cysteine-, arginine-, and glutamine-rich domains and identify novel point mutations in HIV-1B and C sequences from Southern Brazil. This study was the first to investigate the genetic diversity and point mutations within HIV-1 Tat C in a Brazilian cohort. This was an observational, cross-sectional study, which included sequences of HIV-1B (n = 20) and HIV-1C (n = 21) from Southern Brazil. Additionally, 344 HIV-1C sequences were obtained from the Los Alamos database: 29 from Brazil and 315 from Africa, Asia, and Europe. The frequency of C31S substitution on HIV-1 Tat C in Brazil was 82% vs. 10% in the HIV-1B group (p < 0.0001). The frequency of the R57S substitution among the HIV-1C sequences from Brazil was 74% vs. 20% in HIV-1B (p = 0.004), and that of substitution Q63E in HIV-1C was 80% and 20% in HIV-1B (p < 0.0001). The mutation P60Q was more frequent in HIV-1B than in HIV-1C (55% and 6.12%, respectively, p < 0.0001)). Novel point mutations in the HIV-1C and B Tat functional domains were described. The frequency of C31S and other key point mutations in HIV-1 Tat C in Brazil were similar to those described in Africa, although lower than those in India. The Tat-B and C sequences found in Southern Brazil are consistent with biological differences and have potential implications for HIV-1 subtype pathogenesis.

PROFILE OF HCV GENOTYPES AND HIV-SUBTYPES AMONG HIV-COINFECTED PATIENTS IN SOUTHERN BRAZIL

BACKGROUND

Hepatitis B and C virus (HBV and HCV) are the two most common infections among human immunodeficiency virus (HIV)-infected patients.

OBJECTIVE

To identify the frequency of HIV subtypes and HCV genotypes in HIV-coinfected patients.

METHODS

A cross-sectional and retrospective study was carried out into two reference centers in Southern Brazil between January 1, 2002 and June 30, 2016. The Abbott Real Time HCV Genotype II system was used for routine diagnostics to determine the HCV genotype based on dual-target real-time PCR. Proviral HIV-1 RNA was extracted from serum samples and fragments of the pol gene were generated by PCR. The HIV-1 PT and RT gene sequences were submitted to Maximum Likelihood Phylogenetic analysis by collecting reference sequences from the HIV-1 group M subtype of the Los Alamos database.

RESULTS

During the study period, 3340 patients with HIV were diagnosed at both referral centers, of which 4.97% (166/3340) had HBV and/or HCV coinfection. Seroprevalence of HIV-HBV, HIV-HCV and HIV-HBV-HCV was 37.4%, 58.4%, and 4.2%, respectively. HIV-HCV-coinfected patients had a lower median nadir CD4+ T-cell count when compared to HIV-HBV-coinfected patients (P=0.01). Among those coinfected with HCV, HCV-1 (HCV-1) and HCV-3 (HCV-3) genotypes were the most prevalent, being detected in 73.8% and 21.4%, respectively. Among the HCV-1 coinfected patients, 79.3% and 20.1% had subtypes 1a and 1b, respectively. HIV subtype B was the most prevalent in HIV-coinfected patients. There was no significant difference regarding nadir CD4+ T-cell count and HIV viral load when compared to coinfected with HCV-1 with HCV-3, as well as those co-infected with HCV-1a with HCV-1b.

CONCLUSION

In the present study, a higher frequency of subtype B of HIV and HCV-1 were found in HIV-coinfected patients. Further larger-scale and long-term studies are needed to better understand the effect of HCV genotypes in HIV-infected patients.

Inferring population dynamics of HIV-1 subtype C epidemics in Eastern Africa and Southern Brazil applying different Bayesian phylodynamics approaches

The subtype C Eastern Africa clade (C_EA), a particularly successful HIV-1 subtype C lineage, has seeded several sub-epidemics in Eastern African countries and Southern Brazil during the 1960s and 1970s. Here, we characterized the past population dynamics of the major C_EA sub-epidemics in Eastern Africa and Brazil by using Bayesian phylodynamic approaches based on coalescent and birth-death models. All phylodynamic models support similar epidemic dynamics and exponential growth rates until roughly the mid-1980s for all the C_EA sub-epidemics. Divergent growth patterns, however, were supported afterwards. The Bayesian skygrid coalescent model (BSKG) and the birth-death skyline model (BDSKY) supported longer exponential growth phases than the Bayesian skyline coalescent model (BSKL). The BDSKY model uncovers patterns of a recent decline for the C_EA sub-epidemics in Burundi/Rwanda and Tanzania (R_e < 1) and a recent growth for Southern Brazil (R_e > 1); whereas coalescent models infer an epidemic stabilization. To the contrary, the BSKG model captured a decline of Ethiopian C_EA sub-epidemic between the mid-1990s and mid-2000s that was not uncovered by the BDSKY model. These results underscore that the joint use of different phylodynamic approaches may yield complementary insights into the past HIV population dynamics.

C. Pardini MI, Feiteiro MVC, Tolentino FM, Diaz RS, Lopes GISL, Francisco RBL, Véras NMC, Pires AF, Franchini M, Mesquita F, Tanuri A. Brazilian network for HIV Drug Resistance Surveillance (HIV-BresNet): a survey of treatment-naive individuals

INTRODUCTION

In Brazil, more than 487,450 individuals are currently undergoing antiretroviral treatment. In order to monitor the transmission of drug-resistant strains and HIV subtype distribution in the country, this work aimed to estimate its prevalence and to characterize the nationwide pretreatment drug resistance in individuals recently diagnosed with HIV between 2013 and 2015.

METHODS

The HIV threshold survey methodology (HIV-THS, WHO) targeting antiretroviral-naive individuals with recent HIV diagnosis was utilized, and subjects were selected from 51 highly populated cities in all five Brazilian macroregions. The HIV pol genotypic test was performed by genomic sequencing.

RESULTS

We analysed samples from 1568 antiretroviral-naive individuals recently diagnosed with HIV, and the overall transmitted drug resistance (TDR) prevalence was 9.5% (150 sequences). The regional prevalence of resistance according to Brazilian geographical regions was 9.4% in the northeast, 11.2% in the southeast, 6.8% in the central region, 10.2% in the north and 8.8% in the south. The inhibitor-specific TDR prevalence was 3.6% for nucleoside reverse transcriptase inhibitors (NRTIs), 5.8% for non-nucleoside reverse transcriptase inhibitors (NNRTIs) and 1.6% for protease inhibitors (PIs); 1.0% of individuals presented resistance to more than one class of inhibitors. Overall, subtype B was more prevalent in every region except for the southern, where subtype C prevails.

CONCLUSIONS

To the best of our knowledge, this is the first TDR study conducted in Brazil with nationwide representative sampling. The TDR prevalence revealed a moderate rate in the five Brazilian geographical regions, although some cities presented higher TDR prevalence rates, reaching 14% in São Paulo, for example. These results further illustrate the importance of surveillance studies for designing future strategies in primary antiretroviral therapy, aiming to mitigate TDR, as well as for predicting future trends in other regions of the globe where mass antiretroviral (ARV) treatment was implemented.

Divergent HIV-1 strains (CRF92_C2U and CRF93_cpx) co-circulating in the Democratic Republic of the Congo: Phylogenetic insights on the early evolutionary history of subtype C

Molecular epidemiological studies revealed that the epicenter of the HIV pandemic was Kinshasa, the capital city of the Democratic Republic of the Congo (DRC) in Central Africa. All known subtypes and numerous complex recombinant strains co-circulate in the DRC. Moreover, high intra-subtype diversity has been also documented. During two previous surveys on HIV-1 antiretroviral drug resistance in the DRC, we identified two divergent subtype C lineages in the protease and partial reverse transcriptase gene regions. We sequenced eight near full-length genomes and classified them using bootscanning and likelihood-based phylogenetic analyses. Four strains are more closely related to subtype C although within the range of inter sub-subtype distances. However, these strains also have small unclassified fragments and thus were named CRF92_C2U. Another strain is a unique recombinant of CRF92_C2U with an additional small unclassified fragment and a small divergent subtype A fragment. The three remaining strains represent a complex mosaic named CRF93_cpx. CRF93_cpx have two fragments of divergent subtype C sequences, which are not conventional subtype C nor the above described C2, and multiple divergent subtype A-like fragments. We then inferred the time-scaled evolutionary history of subtype C following a Bayesian approach and a partitioned analysis using major genomic regions. CRF92_C2U and CRF93_cpx had the most recent common ancestor with conventional subtype C around 1932 and 1928, respectively. A Bayesian demographic reconstruction corroborated that the subtype C transition to a faster phase of exponential growth occurred during the 1950s. Our analysis showed considerable differences between the newly discovered early-divergent strains and the conventional subtype C and therefore suggested that this virus has been diverging in humans for several decades before the HIV/M diversity boom in the 1950s.

Profile of HIV subtypes in HIV/HBV- and HIV/HCV-coinfected patients in Southern Brazil

INTRODUCTION:

HIV and viral hepatitis infections are major causes of chronic disease worldwide and have some similarities with regard to routes of transmission, epidemiology, front barriers faced during access of treatment, and strategies for a global public health response. The objective was to describe the HIV-1 subtypes, viral tropism and single-nucleotide polymorphisms (SNPs) of interleukin 28B (IL28B) from a case series of HIV/viral hepatitis coinfected patients from southern Brazil.

METHODS:

Clinical and epidemiological data were evaluated by a review of medical records. Periodic blood draws were taken to determine the viral and host characteristics.

RESULTS:

This study included 38 patients with HIV/HBV or HIV/HCV coinfection; the median age was 49 years. Thirty-seven (97.4%) were on antiretroviral therapy, 32 (84.2%) had an undetectable viral load, a median CD4+ T-cell count of 452 cells/mm3. HIV-1 subtyping showed 47.4 and 31.6% of patients with subtypes C and B, respectively. Analysis of viral co-receptor usage showed a predominance of the R5 variant (64.7%), with no significant difference between the subtypes. Twenty patients with HIV/HCV coinfection were eligible to receive HCV therapy with pegylated-interferon-alpha plus ribavirin, and 10/20 (50%) of them achieved sustained virological response. SNPs of IL28B were evaluated in 93.3% of patients with HIV/HCV coinfection, and 17 (60.7%) presented the CC genotype.

CONCLUSIONS:

In the present case series, a higher frequency of HIV subtype C was found in coinfected patients. However such findings need to be prospectively evaluated with the inclusion of data from regional multicenter analyses.

HIV-1 genetic diversity and antiretroviral drug resistance among individuals from Roraima state, northern Brazil

The HIV-1 epidemic in Brazil has spread towards the Northern country region, but little is known about HIV-1 subtypes and prevalence of HIV strains with resistance mutations to antiretrovirals in some of the Northern states. HIV-1 protease (PR) and reverse transcriptase (RT) sequences were obtained from 73 treatment-naive and -experienced subjects followed between 2013 and 2014 at a public health reference unit from Roraima, the northernmost Brazilian state. The most prevalent HIV-1 clade observed in the study population was the subtype B (91%), followed by subtype C (9%). Among 12 HIV-1 strains from treatment-naïve patients, only one had a transmitted drug resistance mutation for NNRTI. Among 59 treatment-experienced patients, 12 (20%) harbored HIV-1 strains with acquired drug resistance mutations (ADRM) that reduce the susceptibility to two classes of antiretroviral drugs (NRTI and NNRTI or NRTI and PI), and five (8%) harbored HIV-1 strains with ADRM that reduced susceptibility to only one class of antiretroviral drugs (NNRTI or PI). No patients harboring HIV strains with reduced susceptibility to all three classes of antiretroviral drugs were detected. A substantial fraction of treatment-experienced patients with (63%) and without (70%) ADRM had undetectable plasma viral loads (<40 copies/ml) at the time of sampling. Among treatment-experienced with plasma viral loads above 2,000 copies/ml, 44% displayed no ADRM. This data showed that the HIV-1 epidemic in Roraima displayed a much lower level of genetic diversity and a lower prevalence of ADRM than that described in other Brazilian states.

Comprehensive Characterization of HIV-1 Molecular Epidemiology and Demographic History in the Brazilian Region Most Heavily Affected by AIDS

The high incidence of AIDS cases and the dominance of HIV-1 subtype C infections are two features that distinguish the HIV-1 epidemic in the two southernmost Brazilian states (Rio Grande do Sul [RS] and Santa Catarina [SC]) from the epidemic in other parts of the country. Nevertheless, previous studies on HIV molecular epidemiology were conducted mainly in capital cities, and a more comprehensive understanding of factors driving this unique epidemic in Brazil is necessary. Blood samples were collected from individuals in 13 municipalities in the Brazilian southern region. HIV-1 env and pol genes were submitted to phylogenetic analyses for assignment of subtype, and viral population phylodynamics were reconstructed by applying Skygrid and logistic coalescent models in a Bayesian analysis. A high prevalence of subtype C was observed in all sampled locations; however, an increased frequency of recombinant strains was found in RS, with evidence for new circulating forms (CRFs). In the SC state, subtype B and C epidemics were associated with distinct exposure groups. Although logistic models estimated similar growth rates for HIV-1 subtype C (HIV-1C) and HIV-1B, a Skygrid plot reveals that the former epidemic has been expanding for a longer time. Our results highlight a consistent expansion of HIV-1C in south Brazil, and we also discuss how heterosexual and men who have sex with men (MSM) transmission chains might have impacted the current prevalence of HIV-1 subtypes in this region.

IMPORTANCE

The AIDS epidemic in south Brazil is expanding rapidly, but the circumstances driving this condition are not well known. A high prevalence of HIV-1 subtype C was reported in the capital cities of this region, in contrast to the subtype B dominance in the rest of the country. This study sought to comparatively investigate the HIV-1 subtype B and C epidemics by sampling individuals from several cities in the two states with the highest AIDS incidences in Brazil. Our analyses showed distinct epidemic growth curves for the two epidemics, and we also found evidence suggesting that separate transmission chains may be impacting the viral phylodynamics and the emergence of new recombinant forms.

The application of HIV molecular epidemiology to public health

HIV is responsible for one of the largest viral pandemics in human history. Despite a concerted global response for prevention and treatment, the virus persists. Thus, urgent public health action, utilizing novel interventions, is needed to prevent future transmission events, critical to eliminating HIV. For public health planning to prove effective and successful, we need to understand the dynamics of regional epidemics and to intervene appropriately. HIV molecular epidemiology tools as implemented in phylogenetic, phylodynamic and phylogeographic analyses have proven to be powerful tools in public health planning across many studies. Numerous applications with HIV suggest that molecular methods alone or in combination with mathematical modelling can provide inferences about the transmission dynamics, critical epidemiological parameters (prevalence, incidence, effective number of infections, Re, generation times, time between infection and diagnosis), or the spatiotemporal characteristics of epidemics. Molecular tools have been used to assess the impact of an intervention and outbreak investigation which are of great public health relevance. In some settings, molecular sequence data may be more readily available than HIV surveillance data, and can therefore allow for molecular analyses to be conducted more easily. Nonetheless, classic methods have an integral role in monitoring and evaluation of public health programmes, and should supplement emerging techniques from the field of molecular epidemiology. Importantly, molecular epidemiology remains a promising approach in responding to viral diseases.

Spatio-Temporal History of HIV-1 CRF35_AD in Afghanistan and Iran

HIV-1 Circulating Recombinant Form 35_AD (CRF35_AD) has an important position in the epidemiological profile of Afghanistan and Iran. Despite the presence of this clade in Afghanistan and Iran for over a decade, our understanding of its origin and dissemination patterns is limited. In this study, we performed a Bayesian phylogeographic analysis to reconstruct the spatio-temporal dispersion pattern of this clade using eligible CRF35_AD gag and pol sequences available in the Los Alamos HIV database (432 sequences available from Iran, 16 sequences available from Afghanistan, and a single CRF35_AD-like pol sequence available from USA). Bayesian Markov Chain Monte Carlo algorithm was implemented in BEAST v1.8.1. Between-country dispersion rates were tested with Bayesian stochastic search variable selection method and were considered significant where Bayes factor values were greater than three. The findings suggested that CRF35_AD sequences were genetically similar to parental sequences from Kenya and Uganda, and to a set of subtype A1 sequences available from Afghan refugees living in Pakistan. Our results also showed that across all phylogenies, Afghan and Iranian CRF35_AD sequences formed a monophyletic cluster (posterior clade credibility> 0.7). The divergence date of this cluster was estimated to be between 1990 and 1992. Within this cluster, a bidirectional dispersion of the virus was observed across Afghanistan and Iran. We could not clearly identify if Afghanistan or Iran first established or received this epidemic, as the root location of this cluster could not be robustly estimated. Three CRF35_AD sequences from Afghan refugees living in Pakistan nested among Afghan and Iranian CRF35_AD branches. However, the CRF35_AD-like sequence available from USA diverged independently from Kenyan subtype A1 sequences, suggesting it not to be a true CRF35_AD lineage. Potential factors contributing to viral exchange between Afghanistan and Iran could be injection drug networks and mass migration of Afghan refugees and labours to Iran, which calls for extensive preventive efforts.

Tracing the origin of a singular HIV-1 CRF45_cpx clade identified in Brazil

The HIV-1 epidemiology has changed over the past decade toward a marked increase in the circulation of strains previously restricted to local epidemics. Recent molecular epidemiological surveys identified some HIV-1 strains of probable African origin circulating in Brazil, including the Circulating Recombinant Form (CRF) 45_cpx, a complex A1/K/U recombinant that circulates in Central Africa. Here, we characterize partial genomic sequences and reconstruct the evolutionary history of HIV-1 CRF45_cpx-related recombinant samples identified in independent studies carried out with HIV+ individuals in Brazil. The sequences were obtained by overlapping PCR amplifications followed by direct sequencing. Recombination profiles were determined by phylogenetic and bootscaning analyses. The evolutionary history was estimated by a Bayesian coalescent-based method using datasets representing the gag, pol and env gene fragments. Six of the 10 samples isolated in Rio de Janeiro showed a CRF45_cpx-like pattern throughout the sequenced genome. The remaining were classified as second-generation recombinants, showing the mosaic patterns: CRF45_cpx/B/D/F1/U, CRF45_cpx/B/F1/U, CRF45_cpx/B/U and CRF45_cpx/F1. All Brazilian CRF45_cpx sequences, except one, formed a monophyletic clade (CRF45-BR), which seems to be the result of a single introduction event that has spread to the Rio de Janeiro, São Paulo and Minas Gerais states and is related to sequences from Argentina, Italy and Belgium. The Bayesian analyses pointed out quite consistent onset dates for CRF45-BR clade (~1984: 1976-1996) in the three gene datasets. These results indicate that the CRF45-BR clade has been circulating in the Southeastern Brazilian region for about 30years, although its presence was not detected until recently due to its very low prevalence. This reinforces the relevance of large-scale molecular surveillance data to identify the emergence of new HIV variants and their impact on local epidemics.

HIV-1 Genetic Diversity and Transmitted Drug Resistance in Antiretroviral Treatment-Naive Individuals from Amapá State, Northern Brazil

The pattern of HIV-1 subtype distribution and prevalence of transmitted drug resistance mutations (TDRM) is heterogeneous across different Brazilian regions. Little information is available about the molecular epidemiologic profile in Northern Brazil. HIV-1 protease (PR) and reverse transcriptase (RT) sequences were obtained from 97 drug-naive HIV-1-infected individuals from Amapá, one of the most isolated Northern Brazilian states, for subtype determination and analysis of drug resistance mutations. The most prevalent HIV-1 clade observed in Amapá was subtype B (74%), followed by subtype F1 (14%), BF1 recombinants (8%), subtype C (1%), CRF31_BC (1%), and CRF02_AG (1%). Only one TDRM (K103N) was detected in a single patient from our study population. This study reveals that the HIV-1 epidemic in Amapá is characterized by a high level of genetic diversity comparable to that observed in major Brazilian cities, but a much lower rate of TDRM (1%).

HIV-1 epidemiology and circulating subtypes in the countryside of South Brazil

INTRODUCTION

Human immunodeficiency virus type 1 (HIV-1) has spread worldwide, with several subtypes and circulating recombinant forms. Brazil has an incidence of 20.5 HIV-1/acquired immunodeficiency syndrome (AIDS) patients per 100,000 inhabitants; however, the Southernmost State of Rio Grande do Sul (RS) has more than twice the number of HIV-1-infected people (41.3/100,000 inhabitants) and a different pattern of subtype frequencies, as previously reported in studies conducted in the capital (Porto Alegre) and its metropolitan region. This study examined HIV-1/AIDS epidemiological and molecular aspects in the countryside of Rio Grande do Sul.

METHODS

Socio-demographic, clinical and risk behavioral characteristics were obtained from HIV-1-positive adult patients using a structured questionnaire. HIV-1 subtypes were determined by nested-polymerase chain reaction (PCR) and sequencing of the pol and env genes.

RESULTS

The study sample included 149 (55% women) patients with a mean age of 41.8 ± 11.9 years. Most (73.8%) patients had a low education level and reported heterosexual practices as the most (91.9%) probable transmission route. HIV-1 subtypes were detected in 26 patients: 18 (69.2%) infected with subtype C, six (23.1%) infected with subtype B and two (7.7%) infected with BC recombinant forms.

CONCLUSIONS

These data highlight the increasing number of HIV-1 subtype C infections in the countryside of South Brazil.

Short Communication: Reassessing the Origin of the HIV-1 CRF02_AG Lineages Circulating in Brazil

HIV-1 CRF02_AG is responsible for at least 8% of the HIV-1 infections worldwide and is distributed mainly in West Africa. CRF02_AG has recently been reported in countries where it is not native, including Brazil. In a previous study including 10 CRF02_AG Brazilian samples, we found at least four independent introductions and two autochthonous transmission networks of this clade in Brazil. As more CRF02_AG samples have been identified in Brazil, we performed a new phylogeographic analysis using a larger dataset than before. A total of 20 Brazilian (18 from Rio de Janeiro and two from São Paulo) and 1,485 African HIV-1 CRF02_AG pol sequences were analyzed using maximum likelihood (ML). The ML tree showed that the Brazilian sequences were distributed in five different lineages. The Bayesian phylogeographic analysis of the Brazilian and their most closely related African sequences (n = 212) placed the origin of all Brazilian lineages in West Africa, probably Ghana, Senegal, and Nigeria. Two monophyletic clades were identified, comprising only sequences from Rio de Janeiro, and their date of origin was estimated at around 1985 (95% highest posterior density: 1979-1992). These results support the existence of at least five independent introductions of the CRF02_AG lineage from West Africa into Brazil and further indicate that at least two of these lineages have been locally disseminated in the Rio de Janeiro state over the past 30 years.

History and origin of the HIV-1 subtype C epidemic in South Africa and the greater southern African region

HIV has spread at an alarming rate in South Africa, making it the country with the highest number of HIV infections. Several studies have investigated the histories of HIV-1 subtype C epidemics but none have done so in the context of social and political transformation in southern Africa. There is a need to understand how these processes affects epidemics, as socio-political transformation is a common and on-going process in Africa. Here, we genotyped strains from the start of the epidemic and applied phylodynamic techniques to determine the history of the southern Africa and South African epidemic from longitudinal sampled data. The southern African epidemic's estimated dates of origin was placed around 1960 (95% HPD 1956-64), while dynamic reconstruction revealed strong growth during the 1970s and 80s. The South African epidemic has a similar origin, caused by multiple introductions from neighbouring countries, and grew exponentially during the 1980s and 90s, coinciding with socio-political changes in South Africa. These findings provide an indication as to when the epidemic started and how it has grown, while the inclusion of sequence data from the start of the epidemic provided better estimates. The epidemic have stabilized in recent years with the expansion of antiretroviral therapy.

Assessing the HIV-1 Epidemic in Brazilian Drug Users: A Molecular Epidemiology Approach

Person who inject illicit substances have an important role in HIV-1 blood and sexual transmission and together with person who uses heavy non-injecting drugs may have less than optimal adherence to anti-retroviral treatment and eventually could transmit resistant HIV variants. Unfortunately, molecular biology data on such key population remain fragmentary in most low and middle-income countries. The aim of the present study was to assess HIV infection rates, evaluate HIV-1 genetic diversity, drug resistance, and to identify HIV transmission clusters in heavy drug users (DUs). For this purpose, DUs were recruited in the context of a Respondent-Driven Sampling (RDS) study in different Brazilian cities during 2009. Overall, 2,812 individuals were tested for HIV, and 168 (6%) of them were positive, of which 19 (11.3%) were classified as recent seroconverters, corresponding to an estimated incidence rate of 1.58%/year (95% CI 0.92–2.43%). Neighbor joining phylogenetic trees from env and pol regions and bootscan analyses were employed to subtype the virus from132 HIV-1-infected individuals. HIV-1 subtype B was prevalent in most of the cities under analysis, followed by BF recombinants (9%-35%). HIV-1 subtype C was the most prevalent in Curitiba (46%) and Itajaí (86%) and was also detected in Brasília (9%) and Campo Grande (20%). Pure HIV-1F infections were detected in Rio de Janeiro (9%), Recife (6%), Salvador (6%) and Brasília (9%). Clusters of HIV transmission were assessed by Maximum likelihood analyses and were cross-compared with the RDS network structure. Drug resistance mutations were verified in 12.2% of DUs. Our findings reinforce the importance of the permanent HIV-1 surveillance in distinct Brazilian cities due to viral resistance and increasing subtype heterogeneity all over Brazil, with relevant implications in terms of treatment monitoring, prophylaxis and vaccine development.

Fast Dating Using Least-Squares Criteria and Algorithms

Phylogenies provide a useful way to understand the evolutionary history of genetic samples, and data sets with more than a thousand taxa are becoming increasingly common, notably with viruses (e.g., human immunodeficiency virus (HIV)). Dating ancestral events is one of the first, essential goals with such data. However, current sophisticated probabilistic approaches struggle to handle data sets of this size. Here, we present very fast dating algorithms, based on a Gaussian model closely related to the Langley–Fitch molecular-clock model. We show that this model is robust to uncorrelated violations of the molecular clock. Our algorithms apply to serial data, where the tips of the tree have been sampled through times. They estimate the substitution rate and the dates of all ancestral nodes. When the input tree is unrooted, they can provide an estimate for the root position, thus representing a new, practical alternative to the standard rooting methods (e.g., midpoint). Our algorithms exploit the tree (recursive) structure of the problem at hand, and the close relationships between least-squares and linear algebra. We distinguish between an unconstrained setting and the case where the temporal precedence constraint (i.e., an ancestral node must be older that its daughter nodes) is accounted for. With rooted trees, the former is solved using linear algebra in linear computing time (i.e., proportional to the number of taxa), while the resolution of the latter, constrained setting, is based on an active-set method that runs in nearly linear time. With unrooted trees the computing time becomes (nearly) quadratic (i.e., proportional to the square of the number of taxa). In all cases, very large input trees (>10,000 taxa) can easily be processed and transformed into time-scaled trees. We compare these algorithms to standard methods (root-to-tip, r8s version of Langley–Fitch method, and BEAST). Using simulated data, we show that their estimation accuracy is similar to that of the most sophisticated methods, while their computing time is much faster. We apply these algorithms on a large data set comprising 1194 strains of Influenza virus from the pdm09 H1N1 Human pandemic. Again the results show that these algorithms provide a very fast alternative with results similar to those of other computer programs. These algorithms are implemented in the LSD software (least-squares dating), which can be downloaded from http://www.atgc-montpellier.fr/LSD/, along with all our data sets and detailed results. An Online Appendix, providing additional algorithm descriptions, tables, and figures can be found in the Supplementary Material available on Dryad at http://dx.doi.org/10.5061/dryad.968t3.

Spatiotemporal dynamics of dissemination of non-pandemic HIV-1 subtype B clades in the Caribbean region

The Human immunodeficiency virus type-1 (HIV-1) epidemic in the Caribbean region is mostly driven by subtype B; but information about the pattern of viral spread in this geographic region is scarce and different studies point to quite divergent models of viral dissemination. In this study, we reconstructed the spatiotemporal and population dynamics of the HIV-1 subtype B epidemic in the Caribbean. A total of 1,806 HIV-1 subtype B pol sequences collected from 17 different Caribbean islands between 1996 and 2011 were analyzed together with sequences from the United States (n = 525) and France (n = 340) included as control. Maximum Likelihood phylogenetic analyses revealed that HIV-1 subtype B infections in the Caribbean are driven by dissemination of the pandemic clade (BPANDEMIC) responsible for most subtype B infections across the world, and older non-pandemic lineages (BCAR) characteristics of the Caribbean region. The non-pandemic BCAR strains account for >40% of HIV-1 infections in most Caribbean islands; with exception of Cuba and Puerto Rico. Bayesian phylogeographic analyses indicate that BCAR strains probably arose in the island of Hispaniola (Haiti/Dominican Republic) around the middle 1960s and were later disseminated to Trinidad and Tobago and to Jamaica between the late 1960s and the early 1970s. In the following years, the BCAR strains were also disseminated from Hispaniola and Trinidad and Tobago to other Lesser Antilles islands at multiple times. The BCAR clades circulating in Hispaniola, Jamaica and Trinidad and Tobago appear to have experienced an initial phase of exponential growth, with mean estimated growth rates of 0.35-0.45 year(-1), followed by a more recent stabilization since the middle 1990s. These results demonstrate that non-pandemic subtype B lineages have been widely disseminated through the Caribbean since the late 1960s and account for an important fraction of current HIV-1 infections in the region.

Tracing the origin and northward dissemination dynamics of HIV-1 subtype C in Brazil

Previous studies indicate that the HIV-1 subtype C epidemic in southern Brazil was initiated by the introduction of a single founder strain probably originating from east Africa. However, the exact country of origin of such a founder strain as well as the origin of the subtype C viruses detected outside the Brazilian southern region remains unknown. HIV-1 subtype C pol sequences isolated in the southern, southeastern and central-western Brazilian regions (n = 209) were compared with a large number (n ~ 2,000) of subtype C pol sequences of African origin. Maximum-likelihood analyses revealed that most HIV-1 subtype C Brazilian sequences branched in a single monophyletic clade (CBR-I), nested within a larger monophyletic lineage characteristic of east Africa. Bayesian analyses indicate that the CBR-I clade most probably originated in Burundi and was introduced into the Paraná state (southern region) around the middle 1970s, after which it rapidly disseminated to neighboring regions. The states of Paraná and Santa Catarina have been the most important hubs of subtype C dissemination, and routine travel and spatial accessibility seems to have been the major driving forces of this process. Five additional introductions of HIV-1 subtype C strains probably originated in eastern (n = 2), southern (n = 2) and central (n = 1) African countries were detected in the Rio de Janeiro state (southeastern region). These results indicate a continuous influx of HIV-1 subtype C strains of African origin into Brazil and also unveil the existence of unrecognized transmission networks linking this country to east Africa.

HIV behind bars: human immunodeficiency virus cluster analysis and drug resistance in a reference correctional unit from southern Brazil

People deprived of liberty in prisons are at higher risk of infection by the human immunodeficiency virus (HIV) due to their increased exposure through intravenous drug use, unprotected sexual activity, tattooing in prison and blood exposure in fights and rebellions. Yet, the contribution of intramural HIV transmission to the epidemic is scarcely known, especially in low- and middle-income settings. In this study, we surveyed 1,667 inmates incarcerated at Presídio Central de Porto Alegre, located in southern Brazil, for HIV infection and molecular characterization. The HIV seroprevalence was 6.6% (110/1,667). Further analyses were carried out on 40 HIV-seropositive inmates to assess HIV transmission clusters and drug resistance within the facility with the use of molecular and phylogenetic techniques. The molecular epidemiology of HIV-1 subtypes observed was similar to the one reported for the general population in southern Brazil, with the predominance of HIV-1 subtypes C, B, CRF31_BC and unique BC recombinants. In particular, the high rate (24%) of URF_BC found here may reflect multiple exposures of the population investigated to HIV infection. We failed to find HIV-infected inmates sharing transmission clusters with each other. Importantly, the analysis of HIV-1 pol genomic fragments evidenced high rates of HIV primary and secondary (acquired) drug resistance and an alarming proportion of virologic failure among patients under treatment, unveiling suboptimal access to antiretroviral therapy (ARV), low ARV adherence and dissemination of drug resistant HIV strains in primary infections. Our results call for immediate actions of public authority to implement preventive measures, serological screening and, for HIV-seropositive subjects, clinical and treatment follow-up in order to control HIV infection and limit the spread of drug resistance strains in Brazilian prisons.

Sexually-transmitted/founder HIV-1 cannot be directly predicted from plasma or PBMC-derived viral quasispecies in the transmitting partner

Objective

Characterization of HIV-1 sequences in newly infected individuals is important for elucidating the mechanisms of viral sexual transmission. We report the identification of transmitted/founder viruses in eight pairs of HIV-1 sexually-infected patients enrolled at the time of primary infection (“recipients”) and their transmitting partners (“donors”).

Methods

Using a single genome-amplification approach, we compared quasispecies in donors and recipients on the basis of 316 and 376 C2V5 env sequences amplified from plasma viral RNA and PBMC-associated DNA, respectively.

Results

Both DNA and RNA sequences indicated very homogeneous viral populations in all recipients, suggesting transmission of a single variant, even in cases of recent sexually transmitted infections (STIs) in donors (n = 2) or recipients (n = 3). In all pairs, the transmitted/founder virus was derived from an infrequent variant population within the blood of the donor. The donor variant sequences most closely related to the recipient sequences were found in plasma samples in 3/8 cases and/or in PBMC samples in 6/8 cases. Although donors were exclusively (n = 4) or predominantly (n = 4) infected by CCR5-tropic (R5) strains, two recipients were infected with highly homogeneous CXCR4/dual-mixed-tropic (X4/DM) viral populations, identified in both DNA and RNA. The proportion of X4/DM quasispecies in donors was higher in cases of X4/DM than R5 HIV transmission (16.7–22.0% versus 0–2.6%), suggesting that X4/DM transmission may be associated with a threshold population of X4/DM circulating quasispecies in donors.

Conclusions

These suggest that a severe genetic bottleneck occurs during subtype B HIV-1 heterosexual and homosexual transmission. Sexually-transmitted/founder virus cannot be directly predicted by analysis of the donor’s quasispecies in plasma and/or PBMC. Additional studies are required to fully understand the traits that confer the capacity to transmit and establish infection, and determine the role of concomitant STIs in mitigating the genetic bottleneck in mucosal HIV transmission.

Peak HIV prevalence: a useful outcome variable for ecological studies

A key question for ecological studies with HIV as the outcome variable is what measure of HIV prevalence to use. In this study we compared the strengths and weaknesses of a variety of measures of HIV prevalence, focusing on peak HIV prevalence and HIV prevalence measured at the same time as the exposure variable. We explored the theoretical problems with each of the two measures of HIV prevalence. We then investigated the difference that substituting one variable for the other made to two published ecological studies. One published study found a strong relationship between migration intensity and HIV prevalence measured at the time the migration was measured. When we repeated the analysis using peak HIV prevalence as the outcome variable, there was no evidence of an association. The second study found evidence of a strong relationship between concurrency and peak HIV prevalence. On repetition of the analysis (but utilizing HIV prevalence at the time the concurrency was measured as the outcome variable) there was no longer a significant association. The choice of HIV measure as outcome variable in ecological studies makes a large difference to the study results. The choice of peak HIV prevalence as outcome variable offers the advantage of avoiding the HIV introduction time bias.

Temporal dynamics of HIV-1 circulating subtypes in distinct exposure categories in southern Brazil

Background

The HIV-1 epidemic in Brazil is predominantly driven by subtype B. However, in Brazilian Southern region subtype C prevails and a relatively high AIDS incidence rate is observed. The aim of the present study was to assess the temporal dynamics of HIV-1 subtypes circulating in patients from distinct exposure categories in Southern Brazil. For this purpose 166 HIV-1 samples collected at the years of 1998 (group I) and 2005–2008 (group II) were analyzed.

Results

Analysis of group I revealed statistically significant (p < 0.05) associations between MSM and subtype B as well as between IDU and subtype C; while no statistical significant association between HIV-1 subtypes and exposure category was verified for group II. An overall temporal increase in the prevalence of subtype C and BC recombinants was observed in both HET and MSM populations, accompanied by a proportional decrease in the prevalence of the pure subtype B.

Conclusions

The present study shows an association between HIV subtypes and exposure categories at the middle 1990s in Southern Brazil. Our findings suggest that MSM and IDU populations might have played a major role in the introduction and initial dissemination of subtypes B and C, respectively, in Southern Brazil. This study also suggests a trend towards homogenization of HIV-1 strains across distinct exposure categories as a consequence of an overall increase in the prevalence of subtype C and BC recombinants in both HET and MSM populations.

Genotypical diversity of HIV clades and central nervous system impairment

The central nervous system (CNS) and the immune system are considered major target organs for HIV infection. The neurological manifestations directly related to HIV are acute viral meningitis, chronic meningitis, HIV associated dementia, vacuolar myelopathy and involvement of the peripheral nervous system. Changes in diagnosis and clinical management have changed the aspect of HIV infection so that it is no longer a fatal disease, and has become a chronic disease requiring sustained medical management. After HAART the incidence of most opportunistic infections, including those affecting the CNS, has dropped markedly. Some studies suggest that neurological involvement of infected patient occur with different frequency, depending on HIV subtype involved in the infection. Subtype C may have reduced neuroinvasive capacity, possibly due to its different primary conformation of HIV transactivating regulatory protein (Tat), involved in monocyte chemotaxis. This review focus on physiopathologic aspects of HIV infection in CNS and its correlation with HIV clades.

The increasing prevalence of HIV-1 subtype C in Southern Brazil and its dispersion through the continent

The HIV-1 has evolved swiftly and the scenario of HIV-1 genetic diversity is constantly changing. In South America, recombinant forms of subtypes B, F1, and BF1 have historically driven the HIV-1 epidemic. In recent years, however, infection with subtype C has gained prominence as its prevalence increased in Southern Brazil as well as neighboring countries. Current studies point to a single introduction of closely related strains as the beginning of the Brazilian subtype C epidemic. However, the place of origin of these strains, date, and route of introduction are under continuous debate as well as the clinical outcomes of the emergence of subtype C. Therefore, this paper reviews the history of the HIV-1 subtype C in Brazil, particularly in the Southern region, covering its demographic and evolutionary history and the possible implications to the Brazilian AIDS epidemic as well as to neighboring countries.

HIV-1 genetic diversity and transmitted drug resistance mutations among patients from the North, Central and South regions of Angola

Background

Angola presents a very complex HIV-1 epidemic characterized by the co-circulation of several HIV-1 group M subtypes, intersubtype recombinants and unclassified (U) variants. The viral diversity outside the major metropolitan regions (Luanda and Cabinda) and the prevalence of transmitted drug resistance mutations (DRM) since the introduction of HAART in 2004, however, has been barely studied.

Methods

One hundred and one individuals from the Central (n = 44), North (n = 35), and South (n = 22) regions of Angola were diagnosed as HIV-1 positive and had their blood collected between 2008 and 2010, at one of the National Referral Centers for HIV diagnosis, the Kifangondo Medical Center, located in the border between the Luanda and Bengo provinces. Angolan samples were genotyped based on phylogenetic and bootscanning analyses of the pol (PR/RT) gene and their drug resistance profile was analyzed.

Results

Among the 101 samples analyzed, 51% clustered within a pure group M subtype, 42% were classified as intersubtype recombinants, and 7% were denoted as U. We observed an important variation in the prevalence of different HIV-1 genetic variants among country regions, with high frequency of subtype F1 in the North (20%), intersubtype recombinants in the Central (42%), and subtype C in the South (45%). Statistically significant difference in HIV-1 clade distribution was only observed in subtype C prevalence between North vs South (p = 0.0005) and Central vs South (p = 0.0012) regions. DRM to NRTI and/or NNRTI were detected in 16.3% of patients analyzed.

Conclusions

These results demonstrate a heterogeneous distribution of HIV-1 genetic variants across different regions in Angola and also revealed an unexpected high frequency of DRM to RT inhibitors in patients that have reported no antiretroviral usage, which may decrease the efficiency of the standard first-line antiretroviral regimens currently used in the country.

Phylodynamics of HIV-1 subtype C epidemic in east Africa

The HIV-1 subtype C accounts for an important fraction of HIV infections in east Africa, but little is known about the genetic characteristics and evolutionary history of this epidemic. Here we reconstruct the origin and spatiotemporal dynamics of the major HIV-1 subtype C clades circulating in east Africa. A large number (n = 1,981) of subtype C pol sequences were retrieved from public databases to explore relationships between strains from the east, southern and central African regions. Maximum-likelihood phylogenetic analysis of those sequences revealed that most (>70%) strains from east Africa segregated in a single regional-specific monophyletic group, here called C_EA. A second major Ethiopian subtype C lineage and a large collection of minor Kenyan and Tanzanian subtype C clades of southern African origin were also detected. A Bayesian coalescent-based method was then used to reconstruct evolutionary parameters and migration pathways of the C_EA African lineage. This analysis indicates that the C_EA clade most probably originated in Burundi around the early 1960s, and later spread to Ethiopia, Kenya, Tanzania and Uganda, giving rise to major country-specific monophyletic sub-clusters between the early 1970s and early 1980s. The results presented here demonstrate that a substantial proportion of subtype C infections in east Africa resulted from dissemination of a single HIV local variant, probably originated in Burundi during the 1960s. Burundi was the most important hub of dissemination of that subtype C clade in east Africa, fueling the origin of new local epidemics in Ethiopia, Kenya, Tanzania and Uganda. Subtype C lineages of southern African origin have also been introduced in east Africa, but seem to have had a much more restricted spread.

Molecular epidemiology of HIV-1 subtypes in India: origin and evolutionary history of the predominant subtype C

Background

India has the third largest HIV-1 epidemic with 2.4 million infected individuals. Molecular epidemiological analysis has identified the predominance of HIV-1 subtype C (HIV-1C). However, the previous reports have been limited by sample size, and uneven geographical distribution. The introduction of HIV-1C in India remains uncertain due to this lack of structured studies. To fill the gap, we characterised the distribution pattern of HIV-1 subtypes in India based on data collection from nationwide clinical cohorts between 2007 and 2011. We also reconstructed the time to the most recent common ancestor (tMRCA) of the predominant HIV-1C strains.

Methodology/Principal Findings

Blood samples were collected from 168 HIV-1 seropositive subjects from 7 different states. HIV-1 subtypes were determined using two or three genes, gag, pol, and env using several methods. Bayesian coalescent-based approach was used to reconstruct the time of introduction and population growth patterns of the Indian HIV-1C. For the first time, a high prevalence (10%) of unique recombinant forms (BC and A1C) was observed when two or three genes were used instead of one gene (p<0.01; p = 0.02, respectively). The tMRCA of Indian HIV-1C was estimated using the three viral genes, ranged from 1967 (gag) to 1974 (env). Pol-gene analysis was considered to provide the most reliable estimate [1971, (95% CI: 1965–1976)]. The population growth pattern revealed an initial slow growth phase in the mid-1970s, an exponential phase through the 1980s, and a stationary phase since the early 1990s.

Conclusions/Significance

The Indian HIV-1C epidemic originated around 40 years ago from a single or few genetically related African lineages, and since then largely evolved independently. The effective population size in the country has been broadly stable since the 1990s. The evolving viral epidemic, as indicated by the increase of recombinant strains, warrants a need for continued molecular surveillance to guide efficient disease intervention strategies.

Evidence of multiple introductions of HIV-1 subtype C in Angola

HIV-1 subtype C is the most prevalent group M clade in southern Africa and some eastern African countries. Subtype C is also the most frequent subtype in Angola (southwestern Africa), with an estimated prevalence of 10-20%. In order to better understand the origin of the HIV-1 subtype C strains circulating in Angola, 31 subtype C pol sequences of Angolan origin were compared with 1950 subtype C pol sequences sampled in other African countries. Phylogenetic analyses reveal that the Angolan subtype C sequences were distributed in 16 different lineages that were widely dispersed among other African strains. Ten subtype C Angolan lineages were composed by only one sequence, while the remaining six clades contain between two and seven sequences. Bayesian phylogeographic analysis indicates that most Angolan clades probably originated in different southern African countries with the exception of one lineage that most likely originated in Burundi. Evolutionary analysis suggests that those Angolan subtype C clades composed by ≥ 2 sequences were introduced into the country between the late 1970s and the mid 2000s. The median estimated time frame for the origin of those Angolan lineages coincides with periods of positive migration influx in Angola that were preceded by phases of negative migratory outflow. These results demonstrate that the Angolan subtype C epidemic resulted from multiple introductions of subtype C viruses mainly imported from southern African countries over the last 30years, some of which have been locally disseminated establishing several autochthonous transmission networks. This study also suggests that population mobility between Angola and southern African countries during civil war (1974-2002) may have played a key role in the emergence of the Angolan subtype C epidemic.

Global analysis of sequence diversity within HIV-1 subtypes across geographic regions

AIMS: HIV-1 sequence diversity can affect host immune responses and phenotypic characteristics such as antiretroviral drug resistance. Current HIV-1 sequence diversity classification uses phylogeny-based methods to identify subtypes and recombinants, which may overlook distinct subpopulations within subtypes. While local epidemic studies have characterized sequence-level clustering within subtypes using phylogeny, identification of new genotype - phenotype associations are based on mutational correlations at individual sequence positions. We perform a systematic, global analysis of position-specific pol gene sequence variation across geographic regions within HIV-1 subtypes to characterize subpopulation differences that may be missed by standard subtyping methods and sequence-level phylogenetic clustering analyses. MATERIALS #ENTITYSTARTX00026; METHODS: Analysis was performed on a large, globally diverse, cross-sectional pol sequence dataset. Sequences were partitioned into subtypes and geographic subpopulations within subtypes. For each subtype, we identified positions that varied according to geography using VESPA (viral epidemiology signature pattern analysis) to identify sequence signature differences and a likelihood ratio test adjusted for multiple comparisons to characterize differences in amino acid (AA) frequencies, including minority mutations. Synonymous nonsynonymous analysis program (SNAP) was used to explore the role of evolutionary selection witihin subtype C. RESULTS: In 7693 protease (PR) and reverse transcriptase (RT) sequences from untreated patients in multiple geographic regions, 11 PR and 11 RT positions exhibited sequence signature differences within subtypes. Thirty six PR and 80 RT positions exhibited within-subtype geography-dependent differences in AA distributions, including minority mutations, at both conserved and variable loci. Among subtype C samples from India and South Africa, nine PR and nine RT positions had significantly different AA distributions, including one PR and five RT positions that differed in consensus AA between regions. A selection analysis of subtype C using SNAP demonstrated that estimated rates of nonsynonymous and synonymous mutations are consistent with the possibility of positive selection across geographic subpopulations within subtypes. CONCLUSION: We characterized systematic genotypic pol differences across geographic regions within subtypes that are not captured by the subtyping nomenclature. Awareness of such differences may improve the interpretation of future studies determining the phenotypic consequences of genetic backgrounds.

Phylogeographic analysis of HIV-1 subtype C dissemination in Southern Brazil

The HIV-1 subtype C has spread efficiently in the southern states of Brazil (Rio Grande do Sul, Santa Catarina and Paraná). Phylogeographic studies indicate that the subtype C epidemic in southern Brazil was initiated by the introduction of a single founder virus population at some time point between 1960 and 1980, but little is known about the spatial dynamics of viral spread. A total of 135 Brazilian HIV-1 subtype C pol sequences collected from 1992 to 2009 at the three southern state capitals (Porto Alegre, Florianópolis and Curitiba) were analyzed. Maximum-likelihood and Bayesian methods were used to explore the degree of phylogenetic mixing of subtype C sequences from different cities and to reconstruct the geographical pattern of viral spread in this country region. Phylogeographic analyses supported the monophyletic origin of the HIV-1 subtype C clade circulating in southern Brazil and placed the root of that clade in Curitiba (Paraná state). This analysis further suggested that Florianópolis (Santa Catarina state) is an important staging post in the subtype C dissemination displaying high viral migration rates from and to the other cities, while viral flux between Curitiba and Porto Alegre (Rio Grande do Sul state) is very low. We found a positive correlation (r² = 0.64) between routine travel and viral migration rates among localities. Despite the intense viral movement, phylogenetic intermixing of subtype C sequences from different Brazilian cities is lower than expected by chance. Notably, a high proportion (67%) of subtype C sequences from Porto Alegre branched within a single local monophyletic sub-cluster. These results suggest that the HIV-1 subtype C epidemic in southern Brazil has been shaped by both frequent viral migration among states and in situ dissemination of local clades.

Phylodynamics of HIV-1 subtype F1 in Angola, Brazil and Romania

The HIV-1 subtype F1 is exceptionally prevalent in Angola, Brazil and Romania. The epidemiological context in which the spread of HIV occurred was highly variable from one country to another, mainly due to the existence of a long-term civil war in Angola and the contamination of a large number of children in Romania. Here we apply phylogenetic and Bayesian coalescent-based methods to reconstruct the phylodynamic patterns of HIV-1 subtype F1 in such different epidemiological settings. The phylogenetic analyses of HIV-1 subtype F1 pol sequences sampled worldwide confirmed that most sequences from Angola, Brazil and Romania segregated in country-specific monophyletic groups, while most subtype F1 sequences from Romanian children branched as a monophyletic sub-cluster (Romania-CH) nested within sequences from adults. The inferred time of the most recent common ancestor of the different subtype F1 clades were as follow: Angola=1983 (1978-1989), Brazil=1977 (1972-1981), Romania adults=1980 (1973-1987), and Romania-CH=1985 (1978-1989). All subtype F1 clades showed a demographic history best explained by a model of logistic population growth. Although the expansion phase of subtype F1 epidemic in Angola (mid 1980s to early 2000s) overlaps with the civil war period (1975-2002), the mean estimated growth rate of the Angolan F1 clade (0.49 year(-1)) was not exceptionally high, but quite similar to that estimated for the Brazilian (0.69 year(-1)) and Romanian adult (0.36 year(-1)) subtype F1 clades. The Romania-CH subtype F1 lineage, by contrast, displayed a short and explosive dissemination phase, with a median growth rate (2.47 year(-1)) much higher than that estimated for adult populations. This result supports the idea that the AIDS epidemic that affected the Romanian children was mainly caused by the spread of the HIV through highly efficient parenteral transmission networks, unlike adult populations where HIV is predominantly transmitted through sexual route.

The origin and evolutionary history of HIV-1 subtype C in Senegal

Background

The classification of HIV-1 strains in subtypes and Circulating Recombinant Forms (CRFs) has helped in tracking the course of the HIV pandemic. In Senegal, which is located at the tip of West Africa, CRF02_AG predominates in the general population and Female Sex Workers (FSWs). In contrast, 40% of Men having Sex with Men (MSM) in Senegal are infected with subtype C. In this study we analyzed the geographical origins and introduction dates of HIV-1 C in Senegal in order to better understand the evolutionary history of this subtype, which predominates today in the MSM population

Methodology/Principal Findings

We used a combination of phylogenetic analyses and a Bayesian coalescent-based approach, to study the phylogenetic relationships in pol of 56 subtype C isolates from Senegal with 3,025 subtype C strains that were sampled worldwide. Our analysis shows a significantly well supported cluster which contains all subtype C strains that circulate among MSM in Senegal. The MSM cluster and other strains from Senegal are widely dispersed among the different subclusters of African HIV-1 C strains, suggesting multiple introductions of subtype C in Senegal from many different southern and east African countries. More detailed analyses show that HIV-1 C strains from MSM are more closely related to those from southern Africa. The estimated date of the MRCA of subtype C in the MSM population in Senegal is estimated to be in the early 80's.

Conclusions/Significance

Our evolutionary reconstructions suggest that multiple subtype C viruses with a common ancestor originating in the early 1970s entered Senegal. There was only one efficient spread in the MSM population, which most likely resulted from a single introduction, underlining the importance of high-risk behavior in spread of viruses.

Genetic diversity and molecular epidemiology of HIV transmission

The high genetic diversity of HIV is one of its most significant features, as it has consequences in global distribution, vaccine design, therapy success, disease progression, transmissibility and viral load testing. Studying HIV diversity helps to understand its origins, migration patterns, current distribution and transmission events. New advances in sequencing technologies based on the parallel acquisition of data are now used to characterize within-host and population processes in depth. Additionally, we have seen similar advances in statistical methods designed to model the past history of lineages (the phylodynamic framework) to ultimately gain better insights into the evolutionary history of HIV. We can, for example, estimate population size changes, lineage dispersion over geographic areas and epidemiological parameters solely from sequence data. In this article, we review some of the evolutionary approaches used to study transmission patterns and processes in HIV and the insights gained from such studies.

Naturally occurring resistance mutations to HIV-1 entry inhibitors in subtypes B, C, and CRF31_BC

BACKGROUND

Entry inhibitors are a class of antiretroviral (ARV) drugs that prevent HIV replication by blocking viral entry into the host cell. The investigation of naturally occurring mutations associated with entry inhibitors across subtypes is required because genetic differences between HIV-1 variants may influence the emergence of drug resistance. Despite the importance of subtype C, which predominates globally, the majority of studies include only subtype B strains.

OBJECTIVES

To investigate the presence of natural resistance mutations to entry inhibitors in HIV-1 subtypes B, C, and CRF31_BC strains.

STUDY DESIGN

Eighty samples were collected from antiretroviral-naïve patients. The gp41 gene from 67 patients and the gp120 gene from 65 patients were partially sequenced. Resistance mutations to entry inhibitors Enfuvirtide, Maraviroc, and Vicriviroc were screened.

RESULTS

ENF resistance-associated mutations of HR1 and HR2 on gp41 were not associated with any subtype. However, the major polymorphisms detected in HR1: N42S, L54M, and A67T were most prevalent in subtype C (p<0.001). Mutations A316T and R315Q in gp120, which are related to MVC and VCV reduced susceptibility respectively, were predominant in subtype C (p<0.05).

CONCLUSIONS

This study shows that many more resistance-associated mutations to entry inhibitors in ARV-naïve patients occur in subtype C compared with subtype B strains. However, further studies will be necessary to elucidate if the differential genetic background of HIV subtypes can affect the efficacy of treatment with entry inhibitors.

Heterosexual transmission of human immunodeficiency virus type 1 subtype C in southern Brazil

BACKGROUND

Human immunodeficiency virus type 1 (HIV-1) subtype B predominates in Brazil, but in the southern region subtype C is the most frequent, followed by subtypes B, F1 and recombinant forms. In southern Brazil, these subtypes co-circulate in subjects with homogeneous demographic and clinical features, enabling a better understanding of the role of HIV-1 subtypes on the characteristics of infection.

OBJECTIVES

To evaluate the prevalence of different HIV-1 subtypes in subjects with recent diagnosis for HIV infection in the extreme south of Brazil, and to study their association with demographic, behavioral, clinical and laboratorial characteristics.

STUDY DESIGN

We have determined the genetic sequence of viral protease and reverse transcriptase (polymerase, connection and RNase H domains) isolated from studied subjects. Viral subtype was inferred by comparison with reference HIV sequences, and recombination was determined with Simplot analysis. The association of HIV-1 subtypes with studied characteristics was evaluated by chi-square, Fisher's exact, Student's t and Kruskal-Wallis tests.

RESULTS

Two hundred and forty-five HIV isolates were molecularly characterized, and the association with variables was studied for 233 (95.1%) patients. Of those, 46.8% followed AIDS defining criteria. HIV-1C was responsible for 56.3% of infections, and was associated with heterosexual transmission (p=0.001) and with higher CD4(+) T-cell counts (p=0.02).

CONCLUSIONS

The molecular epidemiology of HIV-1 in the southernmost Brazil is currently steady with predominance of HIV-1C. This is the first study showing a robust association of the infection by this subtype and heterosexual transmission in the state of Rio Grande do Sul, Brazil.

Co-circulation HIV-1 subtypes B, C, and CRF31_BC in a drug-naïve population from Southernmost Brazil: analysis of primary resistance mutations

In Southernmost Brazil HIV-1 subtypes B, C, and CRF31_BC co-circulates and, since 1996 with the implementation of free access to highly active antiretroviral treatment (HAART), this epidemic is under a quite characteristic selective pressure. The profile of mutations and polymorphisms in the protease (PR) and reverse transcriptase (RT) genes of HIV-1 from untreated patients living in Porto Alegre, Southernmost Brazil were evaluated in order to identify the subtypes and circulating drug resistant genotypes. Blood samples from 99 HIV-1 positive drugs-naïve patients were collected from 2006 to 2007 in Porto Alegre, Brazil. HIV PR and RT genes were amplified, sequenced, and subtyped. The HIV-1 genotyping was performed by partial sequence analysis of the pol in the HIV Drug Resistance Database of Stanford University. Phylogenetic analyses allowed to classify the HIV samples according to their subtypes: B (26.2%), C (39.4%), F (1.1%), CRF31_CB (19.2%), and URF (14.1%). Eight (8.1%) samples showed primary resistance mutations according to the Calibrated Population Resistance tool based in the 2009 Surveillance Drug Resistance Mutation list. Two samples presented resistance mutations to PI, three NRTI and three NNRTI. There was no significant association between presence of resistant genotypes and subtypes, but resistance mutations seem to be less frequent in the subtype C. In addition, this study describes for the first time the mutational profile of CRF31_BC to PI, NRTI, and NNRTI. Genetic analyses of HIV-1 from naïve patients are a promising and important method for surveillance of HIV infection.

The HIV-1 epidemic in Bolivia is dominated by subtype B and CRF12_BF "family" strains

BACKGROUND

Molecular epidemiological studies of HIV-1 in South America have revealed the occurrence of subtypes B, F1 and BF1 recombinants. Even so, little information concerning the HIV-1 molecular epidemiology in Bolivia is available. In this study we performed phylogenetic analyses from samples collected in Bolivia at two different points in time over a 10 year span. We analyzed these samples to estimate the trends in the HIV subtype and recombinant forms over time.

MATERIALS AND METHODS

Fifty one HIV-1 positive samples were collected in Bolivia over two distinct periods (1996 and 2005). These samples were genetically characterized based on partial pol protease/reverse transcriptase (pr/rt) and env regions. Alignment and neighbor-joining (NJ) phylogenetic analyses were established from partial env (n = 37) and all pol sequences using Mega 4. The remaining 14 env sequences from 1996 were previously characterized based on HMA-env (Heteroduplex mobility assay). The Simplot v.3.5.1 program was used to verify intragenic recombination, and SplitsTree 4.0 was employed to confirm the phylogenetic relationship of the BF1 recombinant samples.

RESULTS

Phylogenetic analysis of both env and pol regions confirmed the predominance of "pure" subtype B (72.5%) samples circulating in Bolivia and revealed a high prevalence of BF1 genotypes (27.5%). Eleven out of 14 BF1 recombinants displayed a mosaic structure identical or similar to that described for the CRF12_BF variant, one sample was classified as CRF17_BF, and two others were F1pol/Benv. No "pure" HIV-1 subtype F1 or B" variant of subtype B was detected in the present study. Of note, samples characterized as CRF12_BF-related were depicted only in 2005.

CONCLUSION

HIV-1 genetic diversity in Bolivia is mostly driven by subtype B followed by BF1 recombinant strains from the CRF12_BF "family". No significant temporal changes were detected between the mid-1990s and the mid-2000s for subtype B (76.2% vs 70.0%) or BF1 recombinant (23.8% vs 30.0%) samples from Bolivia.

The origin and diversity of the HIV-1 pandemic

This review examines the enormous progress that has been made in the past decade in understanding the origin of HIV, HIV genetic variability, and the impact of global HIV diversity on the pandemic. Multiple zoonotic transmissions of simian immunodeficiency virus (SIV) have resulted in different HIV lineages in humans. In addition, the high mutation and recombination rates during viral replication result in a great genetic variability of HIV within individuals, as well as within populations, upon which evolutionary selection pressures act. The global HIV pandemic is examined in the context of HIV evolution, and the global diversity of HIV subtypes and recombinants is discussed in detail. Finally, the impact of HIV diversity on pathogenesis, transmission, diagnosis, treatment, the immune response, and vaccine development is reviewed.

Case report of a haemovigilance investigation using phylogenetic analysis of HIV-1 in Brazil

OBJECTIVES

The aim of this work is to provide the first report of a transfusion-acquired HIV-1 infection and to verify transmission from the donor to the recipients using phylogenetic analysis of HIV-1 DNA sequences in a Brazilian blood bank.

BACKGROUND

Although haemovigilance procedures based on phylogenetic analysis of HIV have been reported in several countries, this type of study has yet to be conducted in Latin America.

MATERIALS AND METHODS

Upon identifying a HIV-1-positive repeat blood donor by enzyme immunoassay (EIA) blood screening, all recipients of the donor's previous donation were identified and tested for HIV-1 by EIA, nucleic acid amplification test and HIV-1 DNA sequencing and phylogenetic analysis.

RESULTS

One of the recipients tested positive for HIV-1. The phylogenetic analysis showed a high genetic similarity among the viruses, thus supporting the hypothesis of transmission from the donor to the recipient.

CONCLUSIONS

Phylogenetic analysis of HIV-1 DNA sequences has been a decisive tool in verifying suspected transmission of the virus from blood donor to recipient in Brazil.

Reviewing the history of HIV-1: spread of subtype B in the Americas

The dispersal of HIV-1 subtype B (HIV-1B) is a reflection of the movement of human populations in response to social, political, and geographical issues. The initial dissemination of HIV-1B outside Africa seems to have included the passive involvement of human populations from the Caribbean in spreading the virus to the United States. However, the exact pathways taken during the establishment of the pandemic in the Americas remain unclear. Here, we propose a geographical scenario for the dissemination of HIV-1B in the Americas, based on phylogenetic and genetic statistical analyses of 313 available sequences of the pol gene from 27 countries. Maximum likelihood and Bayesian inference methods were used to explore the phylogenetic relationships between HIV-1B sequences, and molecular variance estimates were analyzed to infer the genetic structure of the viral population. We found that the initial dissemination and subsequent spread of subtype B in the Americas occurred via a single introduction event in the Caribbean around 1964 (1950–1967). Phylogenetic trees present evidence of several primary outbreaks in countries in South America, directly seeded by the Caribbean epidemic. Cuba is an exception insofar as its epidemic seems to have been introduced from South America. One clade comprising isolates from different countries emerged in the most-derived branches, reflecting the intense circulation of the virus throughout the American continents. Statistical analysis supports the genetic compartmentalization of the virus among the Americas, with a close relationship between the South American and Caribbean epidemics. These findings reflect the complex establishment of the HIV-1B pandemic and contribute to our understanding between the migration process of human populations and virus diffusion.

The Use of Bioinformatics for Studying HIV Evolutionary and Epidemiological History in South America

The South American human immunodeficiency virus type 1 (HIV-1) epidemic is driven by several subtypes (B, C, and F1) and circulating and unique recombinant forms derived from those subtypes. Those variants are heterogeneously distributed around the continent in a country-specific manner. Despite some inconsistencies mainly derived from sampling biases and analytical constrains, most of studies carried out in the area agreed in pointing out specificities in the evolutionary dynamics of the circulating HIV-1 lineages. In this paper, we covered the theoretical basis, and the application of bioinformatics methods to reconstruct the HIV spatial-temporal dynamics, unveiling relevant information to understand the origin, geographical dissemination and the current molecular scenario of the HIV epidemic in the continent, particularly in the countries of Southern Cone.

Origin and dynamics of HIV-1 subtype C infection in India

Objective

To investigate the geographical origin and evolution dynamics of HIV-1 subtype C infection in India.

Design

Ninety HIV-1 subtype C env gp120 subtype C sequences from India were compared with 312 env gp120 reference subtype C sequences from 27 different countries obtained from Los Alamos HIV database. All the HIV-1 subtype C env gp120 sequences from India were used for the geographical origin analysis and 61 subtype C env gp120 sequences with known sampling year (from 1991 to 2008) were employed to determine the origin of HIV infection in India.

Methods

Phylogenetic analysis of HIV-1 env sequences was used to investigate the geographical origin and tMRCA of Indian HIV-1 subtype C. Evolutionary parameters including origin date and demographic growth patterns of Indian subtype C were estimated using a Bayesian coalescent-based approach under relaxed molecular clock models.

Findings

The majority of the analyzed Indian and South African HIV-1 subtype C sequences formed a single monophyletic cluster. The most recent common ancestor date was calculated to be 1975.56 (95% HPD, 1968.78–1981.52). Reconstruction of the effective population size revealed three phases of epidemic growth: an initial slow growth, followed by exponential growth, and then a plateau phase approaching present time. Stabilization of the epidemic growth phase correlated with the foundation of National AIDS Control Organization in India.

Interpretation

Indian subtype C originated from a single South African lineage in the middle of 1970s. The current study emphasizes not only the utility of HIV-1 sequence data for epidemiological studies but more notably highlights the effectiveness of community or government intervention strategies in controlling the trend of the epidemic.

Identification of a current hot spot of HIV type 1 transmission in Mongolia by molecular epidemiological analysis

We investigated the current molecular epidemiological status of HIV-1 in Mongolia, a country with very low incidence of HIV-1 though with rapid expansion in recent years. HIV-1 pol (1065 nt) and env (447 nt) genes were sequenced to construct phylogenetic trees. The evolutionary rates, molecular clock phylogenies, and other evolutionary parameters were estimated from heterochronous genomic sequences of HIV-1 subtype B by the Bayesian Markov chain Monte Carlo method. We obtained 41 sera from 56 reported HIV-1-positive cases as of May 2009. The main route of infection was men who have sex with men (MSM). Dominant subtypes were subtype B in 32 cases (78%) followed by subtype CRF02_AG (9.8%). The phylogenetic analysis of the pol gene identified two clusters in subtype B sequences. Cluster 1 consisted of 21 cases including MSM and other routes of infection, and cluster 2 consisted of eight MSM cases. The tree analyses demonstrated very short branch lengths in cluster 1, suggesting a surprisingly active expansion of HIV-1 transmission during a short period with the same ancestor virus. Evolutionary analysis indicated that the outbreak started around the early 2000s. This study identified a current hot spot of HIV-1 transmission and potential seed of the epidemic in Mongolia. Comprehensive preventive measures targeting this group are urgently needed.

Molecular epidemiology of HIV-1 clades in Southern Brazil

Human immunodeficiency virus (HIV) clades B and C account for more than 60% of the HIV-1 infections worldwide. In this paper, we describe the profiles of patients infected with subtypes of HIV-1 from the state of Paraná, Southern Brazil, and correlate them with demographic and epidemiological findings. A retrospective analysis of HIV cases reported from 1999-2007 was also performed. Data from 293 patients were reviewed and 245 were older than 13 (58% female). The distribution of clades was as follows: B 140 (57%), C 67 (23%), F 24 (10%) and mosaic or unique recombinant forms (URFs) 24 (10%). Of the 48 patients younger than 13 years of age (62.5% male), vertical transmission occurred in 46 and the distribution of clades was as follows: B 14 (29%), C 24 (50%), F 7 (15%) and URFs 6 (13%). There was no significant difference in mortality between HIV-1 subtypes. In both groups, patients infected with clade C tended to have higher rates of injection drug use exposure risk.