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

HIV-1 genetic diversity and transmitted drug resistance to integrase strand transfer inhibitors among recently diagnosed adults in Porto Alegre, South Brazil

OBJECTIVES

The HIV-1 genetic diversity and the presence of transmitted drug resistance mutations (TDRMs) against integrase strand transfer inhibitors (INSTIs) were assessed sequencing samples of antiretroviral (ARV)-naive HIV-1-infected individuals from South Brazil.

METHODS

Viral RNA from 42 ART-naive individuals was submitted to complete HIV-1 integrase gene amplification by RT-PCR and sequencing.

RESULTS

Viral strains carrying TDRMs against INSTIs were not detected in the present study. However, the polymorphisms L74M and L74I were each observed in 4.8% of the individuals. These accessory mutations have been reported as putative causes of TDRMs in ART with raltegravir, but only when associated with additional major mutations. When submitted to HIV-1 subtyping, 50% were classified as subtype C, 21% as recombinant BC, 19% as subtype B, 4.8% as subtype F1 and 4.8% as recombinant CF1.

CONCLUSIONS

All 42 ARV-naive individuals were apparently susceptible to INSTIs, included in the Brazilian therapeutic guideline since 2009. To the best of our knowledge, this is the first study to evaluate TDRMs against INSTIs in Brazil. The most prevalent HIV-1 subtypes were subtype C, followed by the recombinant BC and subtype B, which is in agreement with previous studies. However, the presence of subtype F1 and recombinant CF1 reported herein was not observed in previous studies.

Multi-arm RNA junctions encoding molecular logic unconstrained by input sequence for versatile cell-free diagnostics

Applications of RNA-based molecular logic have been hampered by sequence constraints imposed on the input and output of the circuits. Here we show that the sequence constraints can be substantially reduced by appropriately encoded multi-arm junctions of single-stranded RNA structures. To conditionally activate RNA translation, we integrated multi-arm junctions, self-assembled upstream of a regulated gene and designed to unfold sequentially in response to different RNA inputs, with motifs of loop-initiated RNA activators that function independently of the sequence of the input RNAs and that reduce interference with the output gene. We used the integrated RNA system and sequence-independent input RNAs to execute two-input and three-input OR and AND logic in Escherichia coli, and designed paper-based cell-free colourimetric assays that accurately identified two human immunodeficiency virus (HIV) subtypes (by executing OR logic) in amplified synthetic HIV RNA as well as severe acute respiratory syndrome coronavirus-2 (via two-input AND logic) in amplified RNA from saliva samples. The sequence-independent molecular logic enabled by the integration of multi-arm junction RNAs with motifs for loop-initiated RNA activators may be broadly applicable in biotechnology.

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.

HIV-1 molecular diversity in Brazil unveiled by 10 years of sampling by the national genotyping network

HIV-1 has diversified into several subtypes and recombinant forms that are heterogeneously spread around the world. Understanding the distribution of viral variants and their temporal dynamics can help to design vaccines and monitor changes in viral transmission patterns. Brazil has one of the largest HIV-1 epidemics in the western-world and the molecular features of the virus circulating in the country are still not completely known. Over 50,000 partial HIV-1 genomes sampled between 2008 and 2017 by the Brazilian genotyping network (RENAGENO) were analyzed. Sequences were filtered by quality, duplicate sequences per patient were removed and subtyping was performed with online tools and molecular phylogeny. Association between patients’ demographic data and subtypes were performed by calculating the relative risk in a multinomial analysis and trends in subtype prevalence were tested by Pearson correlation. HIV-1B was found to be the most prevalent subtype throughout the country except in the south, where HIV-1C prevails. An increasing trend in the proportion of HIV-1C and F1 was observed in several regions of the country, while HIV-1B tended to decrease. Men and highly educated individuals were more frequently infected by HIV-1B and non-B variants were more prevalent among women with lower education. Our results suggest that socio-demographic factors partially segregate HIV-1 diversity in Brazil while shaping viral transmission networks. Historical events could explain a preferential circulation of HIV-1B among men who have sex with men (MSM) and non-B variants among heterosexual individuals. In view of an increasing male/female ratio of AIDS cases in Brazil in the last 10–15 years, the decrease of HIV-1B prevalence is surprising and suggests a greater penetrance of non-B subtypes in MSM transmission chains.

Transmission dynamics and molecular characterization of HIV-1 epidemic among therapeutic failure patients in Santa Catarina state, southern Brazil

The HIV-1 epidemic in southern Brazil is mostly caused by subtype C, which contrasts the dominance of subtype B in the other regions of the country. Santa Catarina (SC), although the smallest state in the southern region, presents one of the highest incidences and mortality rates in Brazil due to AIDS. This work investigated the HIV-1 molecular diversity and phylogenetic transmission networks in SC state by analyzing a database of 3070 sequences of the national genotyping service. HIV-1C proved to be the most frequent subtype, with a significant increase in prevalence over time. HIV-1B was observed to be associated with highly educated men, suggesting a compartmentalization from other subtypes. Such observation was confirmed by the high frequency of HIV-1B circulating in MSM transmission networks. Identified transmission clusters were majority composed by individuals living up to 25 km away and interstate linkages were mainly between southern neighbor states. In general, individuals between 25 and 40 years old and sequences sampled after 2014 were more likely to be in transmission chains, in agreement with the universal treatment protocol launched in 2014. The present study brings new insights about HIV-1 transmission dynamics in southern Brazil.

Gut Microbiome Profiles and Associated Metabolic Pathways in HIV-Infected Treatment-Naïve Patients

The normal composition of the intestinal microbiota is a key factor for maintaining healthy homeostasis, and accordingly, dysbiosis is well known to be present in HIV-1 patients. This article investigates the gut microbiota profile of antiretroviral therapy-naive HIV-1 patients and healthy donors living in Latin America in a cohort of 13 HIV positive patients (six elite controllers, EC, and seven non-controllers, NC) and nine healthy donors (HD). Microbiota compositions in stool samples were determined by sequencing the V3-V4 region of the bacterial 16S rRNA, and functional prediction was inferred using PICRUSt. Several taxa were enriched in EC compared to NC or HD groups, including Acidaminococcus, Clostridium methylpentosum, Barnesiella, Eubacterium coprostanoligenes, and Lachnospiraceae UCG-004. In addition, our data indicate that the route of infection is an important factor associated with changes in gut microbiome composition, and we extend these results by identifying several metabolic pathways associated with each route of infection. Importantly, we observed several bacterial taxa that might be associated with different viral subtypes, such as Succinivibrio, which were more abundant in patients infected by HIV subtype B, and Streptococcus enrichment in patients infected by subtype C. In conclusion, our data brings a significant contribution to the understanding of dysbiosis-associated changes in HIV infection and describes, for the first time, differences in microbiota composition according to HIV subtypes. These results warrant further confirmation in a larger cohort of patients.

Molecular Surveillance of HIV-1 Infection in Krasnoyarsk Region, Russia: Epidemiology, Phylodynamics and Phylogeography

BACKGROUND

The information about the dynamics of the viral population and migration events that affect the epidemic in different parts of the Russia is insufficient. Possibly, the huge size of the country and limited transport accessibility to certain territories may determine unique traits of the HIV-1 evolutionary history in different regions.

OBJECTIVE

The aim of this study was to explore the genetic diversity of HIV-1 in the Krasnoyarsk region and reconstruct spatial-temporal dynamics of the infection in the region.

METHODS

The demographic and virologic data from 281 HIV-infected individuals in Krasnoyarsk region collected during 2011-2016 were analyzed. The time to the most recent common ancestor, evolutionary rates, population growth, and ancestral geographic movements was estimated using Bayesian coalescent-based methods.

RESULTS

The study revealed moderate diversity of the HIV-1 subtypes found in the region, which included A6 (92.3%), CRF063_02A (4.3%), B (1.1%), and unique recombinants (2.5%). Phylogenetic reconstruction revealed that the A6 subtype was introduced into Krasnoyarsk region by one viral lineage, which arose around 1996.9 (1994.5-1999.5). The phylogeography analysis pointed to Krasnoyarsk city as the geographical center of the epidemic, which further spread to central neighboring districts of the region. At least two epidemic growth phases of subtype A6 were identified which included exponential growth in early-2000s followed by the decline in the mid/late 2010s.

CONCLUSION

This study demonstrates a change in the genetic diversity of HIV-1 in the Krasnoyarsk region. At the beginning of the epidemic, subtype A6 prevailed, subtypes B and CRF063_02A appeared in the region later.

Estimating HIV-1 Genetic Diversity in Brazil Through Next-Generation Sequencing

Approximately 36.7 million people were living with the human immunodeficiency virus (HIV) at the end of 2016 according to UNAIDS, representing a global prevalence rate of 0.8%. In Brazil, an HIV prevalence of 0.24% has been estimated, which represents approximately 830,000 individuals living with the virus. As a touristic and commercial hub in Latin America, Brazil harbors an elevated HIV genetic variability, further contributed by the selective pressure exerted by the host immune system and by antiretroviral treatment. Through the progress of the next-generation sequencing (NGS) techniques, it has been possible to expand the study of HIV genetic diversity, evolutionary, and epidemic processes, allowing the generation of HIV complete or near full-length genomes (NFLG) and improving the characterization of intra- and interhost diversity of viral populations. Greater sensitivity in the detection of viral recombinant forms represents one of the major improvements associated with this development. It is possible to identify unique or circulating recombinant forms using the near full-length viral genomes with increasing accuracy. It also permits the characterization of multiple viral infections within individual hosts. Previous Brazilian studies using NGS to analyze HIV diversity were able to identify several distinct unique and circulating recombinant forms and evidenced dual infections. These data unveiled unprecedented high rates of viral recombination and highlighted that novel recombinants are continually arising in the Brazilian epidemic. In the pooled analysis depicted in this report, HIV subtypes have been determined from HIV-positive patients in five states of Brazil with some of the highest HIV prevalence, three in the Southeast (Rio de Janeiro, São Paulo, and Minas Gerais), one in the Northeast (Pernambuco) and one in the South (Rio Grande do Sul). Combined data analysis showed a significant prevalence of recombinant forms (29%; 101/350), and a similar 26% when only NFLGs were considered. Moreover, the analysis was able to evidence the occurrence of multiple infections in some individuals. Our data highlight the great HIV genetic diversity found in Brazil and unveils a more accurate scenario of the HIV evolutionary dynamics in the region.

Identification of New HIV-1 Circulating Recombinant Forms CRF81_cpx and CRF99_BF1 in Central Western Brazil and of Unique BF1 Recombinant Forms

Intersubtype recombinants classified as circulating recombinant forms (CRFs) or unique recombinant forms (URFs) have been shown to play an important role in the complex and dynamic Brazilian HIV/AIDS epidemic. Previous pol region studies (2003-2013) in 828 patients from six states from Central Western, Northern and Northeastern Brazil reported variable rates of BF1, F1CB, BC, and CF1 mosaics. In this study HIV-1 subtype diversity BF1, F1CB, BC, and CF1 recombinants in pol were analyzed. Full/near-full/partial genome sequences were generated from F1CB and BF1 recombinants. Genomic DNA extracted from whole blood was used in nested-PCR to amplify four overlapping fragments encompassing the full HIV-1 genome. Phylogenetic trees were generated using the neighbor-joining/NJ method (MEGA 6.0). The time of the most recent common ancestor (TMRCA) of F1CB and BF1 clades was estimated using a Bayesian Markov Chain Monte Carlo approach (BEAST v1.8; BEAGLE). Bootscanning was used for recombination analyses (Simplot v3.5.1); separate NJ phylogenetic analysis of fragments confirmed subtypes. The phylogenetic analyses of protease/reverse-transcriptase sequences in 828 patients revealed 76% subtype B (n = 629), 6.4% subtype C (n = 53), 4.2% subtype F1 (n = 35), 13.4% intersubtype recombinants: 10.5% BF1 (n = 87), 2.3% BC (n = 19), 0.4% F1CB (n = 3), and 0.2% CF1 (n = 2). Two full and one partial BF1C genomes allowed the characterization of the CRF81_cpx that has 9 breakpoints dividing the genome into 10 subregions. Basic Local Alignment Search Tool searches (Los Alamos HIV Sequence Database) identified six other sequences with the same recombination profile in pol, five from Brazil, and one from Italy. The estimated median TMRCA of CRF81_cpx was 1999 (1992-2003). CRF60_BC-like sequences, originally described in Italy, were also found. Two full and one near full-length BF1 genomes led to the characterization of the new CRF99_BF1 that has six recombination breakpoints dividing the genome into seven subregions. Two new URFs BF1, with six recombination breakpoints and seven subregions were also characterized. The description of the first Brazilian BF1C CRF81_cpx and of the new CRF99_BF1 corroborate the important role of CRFs in the HIV/AIDS epidemic throughout Brazil. Our data also highlight the value of HIV-1 full-genome sequence studies in order to fully reveal the complexity of the epidemic in a huge country as Brazil.

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.

Molecular Epidemiology of the HIV Epidemic in Three German Metropolitan Regions - Cologne/Bonn, Munich and Hannover, 1999-2016

Using HIV sequence data to characterize clusters of HIV transmission may provide insight into the epidemic. Phylogenetic and network analyses were performed to infer putative relationships between HIV-1 partial pol sequences from 2,774 individuals receiving care in three German regions between 1999-2016. The regions have in common that they host some of the largest annual festivals in Europe (Carnival and Oktoberfest). Putative links with sequences (n = 150,396) from the Los Alamos HIV Sequence database were evaluated. A total of 595/2,774 (21.4%) sequences linked with at least one other sequence, forming 184 transmission clusters. Clustering individuals were significantly more likely to be younger, male, and report sex with men as their main risk factor (p < 0.001 each). Most clusters (77.2%) consisted exclusively of men; 41 (28.9%) of these included men reporting sex with women. Thirty-two clusters (17.4%) contained sequences from more than one region; clustering men were significantly more likely to be in a position bridging regional HIV epidemics than clustering women (p = 0.027). We found 236 clusters linking 547 sequences from our sample with sequences from the Los Alamos database (n = 1407; 31% from other German centres). These results highlight the pitfalls of focusing HIV prevention efforts on specific risk groups or specific locales.

Spatiotemporal and demographic history of the HIV-1 circulating recombinant form CRF31_BC in Brazil

CRF31_BC is an HIVs-1 recombinant form very prevalent in the southernmost capital city of Brazil, Porto Alegre. Recent studies have been describing a growing number of cases of infection by CRF31_BC in other Brazilian cities and countries, suggesting a process of expansion of this strain. Aiming to describe the city of origin, dispersion routes and demographic history of CRF31_BC, this study analyzed all HIV-1 CRF31_BC and Brazilian BC mosaic publicly available sequences. CRF31_BC classification was performed by bootscanning and tree reconstruction methods. Bayesian phylogeographic and phylodynamic model approaches were used to reconstruct the spatiotemporal and demographic history of 95 sequences identified as CRF31_BC-like. Porto Alegre was estimated to be the origin and center of the dispersion of the CRF31_BC for most of the analyzed locations. However, some viral transitions independent from Porto Alegre were observed in other cities from the Rio Grande do Sul state and also in other Brazilian states. The estimated CRF31_BC epidemic growth rate was similar to subtype C and B in Brazil. Our findings suggest that CRF31_BC, although mostly prevalent in south region, is circulating nation-wide with some localities presenting autochthonous transmissions.

Time-scale of minor HIV-1 complex circulating recombinant forms from Central and West Africa

Background

Several HIV-1 circulating recombinant forms with a complex mosaic structure (CRFs_cpx) circulate in central and western African regions. Here we reconstruct the evolutionary history of some of these complex CRFs (09_cpx, 11_cpx, 13_cpx and 45_cpx) and further investigate the dissemination dynamic of the CRF11_cpx clade by using a Bayesian coalescent-based method.

Results

The analysis of two HIV-1 datasets comprising 181 pol (36 CRF09_cpx, 116 CRF11_cpx, 20 CRF13_cpx and 9 CRF45_cpx) and 125 env (12 CRF09_cpx, 67 CRF11_cpx, 17 CRF13_cpx and 29 CRF45_cpx) sequences pointed to quite consistent onset dates for CRF09_cpx (~1966: 1958–1979), CRF11_cpx (~1957: 1950–1966) and CRF13_cpx (~1965: 1958–1973) clades; while some divergence was found for the estimated date of origin of CRF45_cpx clade [pol = 1970 (1964–1976); env = 1960 (1952–1969)]. Phylogeographic reconstructions indicate that the HIV-1 CRF11_cpx clade most probably emerged in Cameroon and from there it was first disseminated to the Central Africa Republic and Chad in the early 1970s and to other central and western African countries from the early 1980s onwards. Demographic reconstructions suggest that the CRF11_cpx epidemic grew between 1960 and 1990 with a median exponential growth rate of 0.27 year⁻¹, and stabilized after.

Conclusions

These results reveal that HIV-1 CRFs_cpx clades have been circulating in Central Africa for a period comparable to other much more prevalent HIV-1 group M lineages. Cameroon was probably the epicenter of dissemination of the CRF11_cpx clade that seems to have experienced a long epidemic growth phase before stabilization. The epidemic growth of the CRF11_cpx clade was roughly comparable to other HIV-1 group M lineages circulating in Central Africa.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0824-8) contains supplementary material, which is available to authorized users.