Ultra-Deep Sequencing of HIV-1 near Full-Length and Partial Proviral Genomes Reveals High Genetic Diversity among Brazilian Blood Donors

Description of the new HIV-1 intersubtype B/C circulating recombinant form (CRF146_BC) detected in Brazil

BACKGROUND

The human immunodeficiency virus 1 (HIV-1) infections in Brazil are predominantly caused by two subtypes, B and C.

OBJECTIVES

Here we present the characterisation of a novel HIV-1 recombinant form, indicating a new Brazilian CRF_BC, named CRF146_BC.

METHODS

RDP, JphMM and Simplot recombination tools were used to evaluate the mosaic pattern.

FINDINGS

In this work, we identified three HIV-1 nucleotide sequences previously classified as unique recombinant forms (URFs), plus one new partial genome sharing the same BC recombination pattern. The mosaic genome is almost entirely represented by the subtype C sequence, with a small subtype B recombination region in the pol gene, at the Integrase level. The phylogenetic analyses strongly indicate a common origin between the strains, which were isolated in Rio Grande do Sul, Rio de Janeiro and Bahia states.

MAIN CONCLUSIONS

Thus, the new HIV-1 CRF146_BC is circulating in three different Brazilian regions: South, Southeast and Northeast.

Genetic diversity of hepatitis B virus quasispecies in different biological compartments reveals distinct genotypes

The selection pressure imposed by the host immune system impacts hepatitis B virus (HBV) quasispecies variability. This study evaluates HBV genetic diversity in different biological fluids. Twenty paired serum, oral fluid, and DBS samples from chronic HBV carriers were analyzed using both Sanger and next generation sequencing (NGS). The mean HBV viral load in serum was 5.19 ± 4.3 log IU/mL (median 5.29, IQR 3.01-7.93). Genotype distribution was: HBV/A1 55% (11/20), A2 15% (3/20), D3 10% (2/20), F2 15% (3/20), and F4 5% (1/20). Genotype agreement between serum and oral fluid was 100% (genetic distances 0.0-0.006), while that between serum and DBS was 80% (genetic distances 0.0-0.115). Two individuals presented discordant genotypes in serum and DBS. Minor population analysis revealed a mixed population. All samples displayed mutations in polymerase and/or surface genes. Major population analysis of the polymerase pointed to positions H122 and M129 as the most polymorphic (≥ 75% variability), followed by V163 (55%) and I253 (50%). Neither Sanger nor NGS detected any antiviral primary resistance mutations in the major populations. Minor population analysis, however, demonstrated the rtM204I resistance mutation in all individuals, ranging from 2.8 to 7.5% in serum, 2.5 to 6.3% in oral fluid, and 3.6 to 7.2% in DBS. This study demonstrated that different fluids can be used to assess HBV diversity, nonetheless, genotypic differences according to biological compartments can be observed.

Identification of a New HIV-1 Circulating Recombinant Form CRF112_01B Strain in Baoding City, Hebei Province, China

BACKGROUND

A large number of HIV-1 recombinants that originated from CRF01_AE and B strains are constantly emerging in men who have sex with men populations in China and deserve more attention and further monitoring.

OBJECTIVE

To analyze the near-full-length genome structure and recombination characteristics of a new HIV-1 strain (BD226AJ) detected in Baoding City and determine its subtype.

CASE REPRESENTATION

Viral RNA was extracted from a blood sample collected from an infected individual and reverse transcribed to cDNA. Two overlapping segments of the HIV-1 genome were amplified using a near-endpoint dilution method and sequenced. Recombinant breakpoints were determined using RIP, jpHMM, and SimPlot 3.5.1 software. MEGA v6.0 was used to construct a neighbor-joining phylogenetic tree to determine the homology relationships of this strain.

RESULTS AND DISCUSSION

We obtained 8830 nucleotides (nt) of the HIV-1 genome sequence by amplification and sequencing, and four recombinant fragments were identified by recombination analysis, namely CRF01_AE (HXB2, 823-4224 nt), subtype B (HXB2, 4225-5991 nt), CRF01_AE (HXB2, 5992-9295 nt), and subtype B (HXB2, 9296-9406 nt). The BLAST results showed that 96% of the sequence was similar to CRF112_01B. The jpHMM results confirmed that BD226AJ was the CRF112_01B strain.

CONCLUSION

Our results confirm the first epidemic of CRF112_01B in Hebei Province. This finding suggests that HIV-1 CRF112_01B may have been introduced into Hebei by men who have sex with men and indicates that the epidemic trend of this strain should be closely monitored.

Viruses Previously Identified in Brazil as Belonging to HIV-1 CRF72_BF1 Represent Two Closely Related Circulating Recombinant Forms, One of Which, Designated CRF122_BF1, Is Also Circulating in Spain

Circulating recombinant forms (CRFs) are important components of the HIV-1 pandemic. Those derived from recombination between subtype B and subsubtype F1, with 18 reported, most of them of South American origin, are among the most diverse. In this study, we identified a HIV-1 BF1 recombinant cluster that is expanding in Spain, transmitted mainly via heterosexual contact, which, analyzed in near full-length genomes in four viruses, exhibited a coincident BF1 mosaic structure, with 12 breakpoints, that fully coincided with that of two viruses (10BR_MG003 and 10BR_MG005) from Brazil, previously classified as CRF72_BF1. The three remaining Brazilian viruses (10BR_MG002, 10BR_MG004, and 10BR_MG008) previously identified as CRF72_BF1 exhibited mosaic structures highly similar, but not identical, to that of the Spanish viruses and to 10BR_MG003 and 10BR_MG005, with discrepant subtypes in two short genome segments, located in pol and gp120^env. Based on these results, we propose that the five viruses from Brazil previously identified as CRF72_BF1 actually belong to two closely related CRFs, one comprising 10BR_MG002, 10BR_MG004, and 10BR_MG008, which keep their CRF72_BF1 designation, and the other, designated CRF122_BF1, comprising 10BR_MG003, 10BR_MG005, and the viruses of the identified Spanish cluster. Three other BF1 recombinant genomes, two from Brazil and one from Italy, previously identified as unique recombinant forms, were classified as CRF72_BF1. CRF122_BF1, but not CRF72_BF1, was associated with protease L89M substitution, which was reported to contribute to antiretroviral drug resistance. Phylodynamic analyses estimate the emergence of CRF122_BF1 in Brazil around 1987. Given their close phylogenetic relationship and similar structures, the grouping of CRF72_BF1 and CRF122_BF1 in a CRF family is proposed.

HIV-1 genetic diversity and drug resistance mutations in the northern Brazilian region

Brazil is a huge continental country with striking geographic differences which are well illustrated in the HIV/AIDS epidemic. Contrasting with the significant decline in the national AIDS detection rate in the last decade, a linear growth has been reported in the Northern region. Despite its public health and epidemiologic importance, there is scarce HIV-1 molecular data from Northern Brazil. This scoping review summarizes recent epidemiologic data with special emphasis on HIV-1 genetic diversity and antiretroviral drug resistance mutations in patients from the seven Northern states of Brazil. Studies from the Northern Brazil on different HIV-1 genomic regions, mostly pol (protease/reverse transcriptase) sequences of naïve/antiretroviral treated adults/children were retrieved from PubMed/MEDLINE electronic database. These studies indicate a consistent molecular profile largely dominated by HIV-1 subtype B with minor contribution of subtypes F1 and C and infrequent detection of other subtypes (A1, D, K), recombinants (BF1, BC), circulating recombinant forms (CRF) as the new CRF90_BF1 and CRF02_AG-like, CRF28–29_BF-like, CRF31_BC-like, and a potential new CRF_BF1. This pattern indicates a founder effect of subtype B and the introduction of non-B-subtypes and recombinants probably generated in the Southern/Southeastern regions. In naïve populations transmitted drug resistance (TDR) can impact the outcome of first-line antiretroviral treatment and prophylactic/preventive regimens. In the Northern region TDR rates are moderate while patients failing highly active antiretroviral therapy (HAART) showed high prevalence of acquired drug resistance mutations. The limited HIV-1 molecular data from Northern Brazil reflects the great challenges to generate comprehensive scientific data in isolated, underprivileged areas. It also highlights the need to invest in local capacity building which supported by adequate infrastructure and funding can promote robust research activities to help reduce the scientific asymmetries in the Northern region. Currently the impacts of the overwhelming COVID-19 pandemic on the expanding HIV/AIDS epidemic in Northern Brazil deserves to be closely monitored.

Validation of Variant Assembly Using HAPHPIPE with Next-Generation Sequence Data from Viruses

Next-generation sequencing (NGS) offers a powerful opportunity to identify low-abundance, intra-host viral sequence variants, yet the focus of many bioinformatic tools on consensus sequence construction has precluded a thorough analysis of intra-host diversity. To take full advantage of the resolution of NGS data, we developed HAplotype PHylodynamics PIPEline (HAPHPIPE), an open-source tool for the de novo and reference-based assembly of viral NGS data, with both consensus sequence assembly and a focus on the quantification of intra-host variation through haplotype reconstruction. We validate and compare the consensus sequence assembly methods of HAPHPIPE to those of two alternative software packages, HyDRA and Geneious, using simulated HIV and empirical HIV, HCV, and SARS-CoV-2 datasets. Our validation methods included read mapping, genetic distance, and genetic diversity metrics. In simulated NGS data, HAPHPIPE generated pol consensus sequences significantly closer to the true consensus sequence than those produced by HyDRA and Geneious and performed comparably to Geneious for HIV gp120 sequences. Furthermore, using empirical data from multiple viruses, we demonstrate that HAPHPIPE can analyze larger sequence datasets due to its greater computational speed. Therefore, we contend that HAPHPIPE provides a more user-friendly platform for users with and without bioinformatics experience to implement current best practices for viral NGS assembly than other currently available options.

Measurement error and variant-calling in deep Illumina sequencing of HIV

MOTIVATION

Next-generation deep sequencing of viral genomes, particularly on the Illumina platform, is increasingly applied in HIV research. Yet, there is no standard protocol or method used by the research community to account for measurement errors that arise during sample preparation and sequencing. Correctly calling high and low-frequency variants while controlling for erroneous variants is an important precursor to downstream interpretation, such as studying the emergence of HIV drug-resistance mutations, which in turn has clinical applications and can improve patient care.

RESULTS

We developed a new variant-calling pipeline, hivmmer, for Illumina sequences from HIV viral genomes. First, we validated hivmmer by comparing it to other variant-calling pipelines on real HIV plasmid datasets. We found that hivmmer achieves a lower rate of erroneous variants, and that all methods agree on the frequency of correctly called variants. Next, we compared the methods on an HIV plasmid dataset that was sequenced using Primer ID, an amplicon-tagging protocol, which is designed to reduce errors and amplification bias during library preparation. We show that the Primer ID consensus exhibits fewer erroneous variants compared to the variant-calling pipelines, and that hivmmer more closely approaches this low error rate compared to the other pipelines. The frequency estimates from the Primer ID consensus do not differ significantly from those of the variant-calling pipelines.

AVAILABILITY AND IMPLEMENTATION

hivmmer is freely available for non-commercial use from https://github.com/kantorlab/hivmmer.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Characterization of HIV-1 genetic diversity and antiretroviral resistance in the state of Maranhão, Northeast Brazil

The HIV-1 epidemic in Brazil has been growing in northeast and north regions, particularly an increase in AIDS cases among the younger male population has been observed. This study aims to characterize the HIV-1 genetic diversity and to evaluate its antiretroviral resistance profile among individuals presenting virological failure in the state of Maranhão—Brazil. HIV-1 pol gene sequences from 633 patients on antiretroviral therapy were obtained from the Department of Surveillance, Prevention and Control of Sexually Transmitted Infections, HIV/AIDS and Viral Hepatitis of the Brazilian Ministry of Health. Phylogenetic and recombination analyses were performed to characterize viral genetic diversity. The presence of antiretroviral resistance mutations was assessed using the HIV Drug Resistance Database online platform of Stanford University. A predominance of subtype B (84.5%) was observed, followed by recombinant BF (9.5%), where more than half of the sequences were dispersed in 3 clusters. Antiretroviral resistance was detected in 74.1% of the sequences, and it was significantly higher for nucleoside analogue reverse-transcriptase inhibitors (NRTIs) than for non-nucleoside analogue reverse-transcriptase inhibitors (NNRTIs) and protease inhibitors (PIs). Inference of putative transmissions clusters identified 11 clusters with 22 query sequences (22/633, 3.5%). Thus, we conclude that continuous monitoring of the molecular epidemiology of HIV-1 is essential for prevention strategies, epidemic control, and treatment adequacy.

Homogenous HIV-1 subtype B from the Brazilian Amazon with infrequent diverse BF1 recombinants, subtypes F1 and C among blood donors

In the last decade a growing HIV/AIDS epidemic with increased incidence and AIDS-related mortality has been reported in Northern Brazil from which molecular data are scarce. Also, apparently healthy, adult blood donors, recently diagnosed with HIV-1 represent important sentinel populations for molecular studies. This cross-sectional study describes HIV-1 subtypes in blood donors from three reference public blood centers located in three States in Northern Brazil. HIV-1 pol sequencing (protease/PR, reverse transcriptase/RT) was performed on plasma samples of HIV-1 positive donors from HEMOAM, Manaus, Amazonas (n = 198), HEMERON, Porto Velho, Rondônia (n = 20) and HEMORAIMA, Boa Vista, Roraima (n = 9) collected from 2011-2017. HIV-1 subtypes were identified by REGA, phylogenetic inference; recombinant viruses were characterized by SIMPLOT. Young, single, males predominated, around half was first-time donors. Syphilis co-infection was detected in 17% (39 out of 227), 8% (18 out of 227) was anti-HBc positive. Subtype B represented ≥ 90% in Amazonas, Rondônia and Roraima, subtype C (3.1%) was found in Amazonas and Rondônia; subtype F1 (0.9%) and BF1 recombinants (5.3%) were only detected in Amazonas. Subtype B sequences from Amazonas (n = 179), Rondônia (n = 18) and Roraima (n = 9) were combined with viral strains representative of the BPANDEMIC (n = 300) and BCARIBBEAN/BCAR (n = 200) lineages. The BPANDEMIC lineage predominated (78%) although BCAR lineages were frequent in Roraima (56%) and Amazonas (22%). Subtype C and subtype F1 sequences identified here clustered within Brazilian CBR and F1BR lineages, respectively. Twelve BF1 mosaics showed 11 different recombination profiles: six were singleton unique-recombinant-forms/URFs, one displays a CRF28/29_BF-like recombinant pattern and the remaining four BF1 isolates branched with other Brazilian BF1 viruses previously described and may represent putative new CRF_BF1 from Northern Brazil. Our study shows a highly homogeneous molecular pattern with prevalent subtype B, followed by BF1, and sporadic subtype C and F1 in blood donors from the Northern region. Surveillance studies are important to monitor HIV-1 diversity which can reveal patterns of viral dissemination, especially in a highly endemic, remote and geographically isolated region as Northern Brazil.

Molecular and geographic characterization of hiv-1 bf recombinant viruses

The Human Immunodeficiency Virus Type 1 (HIV-1) presents a wide genetic variability, which is represented by four groups, nine subtypes of group M and several recombinant forms. Among these, the BF recombinants have been distinguished by a high global dispersion and an increase in number and diversity. To date, 15 BF Circulating Recombinant Forms (CRFs) and diverse BF Unique Recombinant Forms (URFs) have been described. In Brazil, nine CRF_BF have been identified. The aim of this work was to perform molecular and geographic characterization of HIV-1 BF recombinant strains. Near full-length genomes of 265 BF recombinant viruses were collected from public databases and molecular analyses were performed. These sequences were originally retrieved between 1993-2006 and isolated from 16 countries (51.3% from Brazil). Diagnostic's year analysis showed that BF recombinants circulate in Brazil since at least 1985. Most sequences displayed recombination in the pol (84.9%), gag (69.3%) and env (51.4%) regions. The subtype B predominated in all accessory and regulatory genes, except in vif, in which the F subtype was predominant (40.4%). Twelve regions with a recombination rate higher than 10% were identified, especially one region inside p24 gene (1359-1397) whose recombination was present in more than 30% of the sequences. Coreceptor usage prediction during viral entry showed that BF recombinants preferentially use CCR5 (67.2%) and the most frequent tetrapeptides found in the V3 loop were GPGR (47.9%) and GPGQ (21.1%). The frequency of X4/dual viruses was lower amongst F subtype (25.8%) V3 sequences, compared with B subtype (43%). In addition, mutations associated with intermediate or high resistance levels to PI (10.6%), NRTI (15.0%), NNRTI (14.0%) and INSTI (2.6%) were identified. The great diversity of the recombination patterns evidences that the recombination between the subtypes B and F is frequent, reflecting a probable high rate of dual infection and the acquisition of advantageous characteristics for viral fitness.

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.

Bioinformatic data processing pipelines in support of next-generation sequencing-based HIV drug resistance testing: the Winnipeg Consensus

INTRODUCTION

Next-generation sequencing (NGS) has several advantages over conventional Sanger sequencing for HIV drug resistance (HIVDR) genotyping, including detection and quantitation of low-abundance variants bearing drug resistance mutations (DRMs). However, the high HIV genomic diversity, unprecedented large volume of data, complexity of analysis and potential for error pose significant challenges for data processing. Several NGS analysis pipelines have been developed and used in HIVDR research; however, the absence of uniformity in data processing strategies results in lack of consistency and comparability of outputs from different pipelines. To fill this gap, an international symposium on bioinformatic strategies for NGS-based HIVDR testing was held in February 2018 in Winnipeg, Canada, convening laboratory scientists, bioinformaticians and clinicians involved in four recently developed, publicly available NGS HIVDR pipelines. The goal of this symposium was to establish a consensus on effective bioinformatic strategies for NGS data management and its use for HIVDR reporting.

DISCUSSION

Essential functionalities of an NGS HIVDR pipeline were divided into five analytic blocks: (1) NGS read quality control (QC)/quality assurance (QA); (2) NGS read alignment and reference mapping; (3) HIV variant calling and variant QC; (4) NGS HIVDR reporting; and (5) extended data applications and additional considerations for data management. The consensuses reached among the participants on all major aspects of these blocks are summarized here. They encompass not only recommended data management and analysis strategies, but also detailed bioinformatic approaches that help ensure accuracy of the derived HIVDR analysis outputs for both research and potential clinical use.

CONCLUSIONS

While NGS is being adopted more broadly in HIVDR testing laboratories, data processing is often a bottleneck hindering its generalized application. The proposed standardization of NGS read QC/QA, read alignment and reference mapping, variant calling and QC, HIVDR reporting and relevant data management strategies in this "Winnipeg Consensus" may serve as a starting guideline for NGS HIVDR data processing that informs the refinement of existing pipelines and those yet to be developed. Moreover, the bioinformatic strategies presented here may apply more broadly to NGS data analysis of microbes harbouring significant genomic diversity.

HIV-1 infection among crack cocaine users in a region far from the epicenter of the HIV epidemic in Brazil: Prevalence and molecular characteristics

Brazil has the largest cocaine market in South America, and crack cocaine use is closely associated with HIV-1 infection. This study investigated the prevalence, risk factors, and HIV-1 subtypes, including recombinant forms and mutations associated with drug resistance, among crack cocaine users in Central-West Brazil. We recruited 600 crack cocaine users admitted to a referral hospital in Goiânia for psychiatric disorders. The participants were interviewed; blood samples were collected for anti-HIV-1/2 serological screening. HIV-1 pol gene sequences (entire protease [PR] and partial reverse transcriptase [RT]) were obtained from plasma RNA. HIV-1 subtypes, recombinant viruses, transmitted drug resistance (TDR), and secondary drug resistance mutations were investigated. The median participant age was 30 years (range, 18-68 years); most were male, single, unemployed, and of mixed races. Among them, 2.8% (17/600) were HIV-1 positive: 2.2% of men (11/507) and 6.5% of women (6/93). The main predictors of HIV-1 seropositivity were a sexual partner with HIV infection, irregular condom use, and previous homelessness. HIV-1 pol sequences (12/17) indicated the predominance of subtype B (n = 7), followed by recombinant forms FPR/BRT (n = 1) and BPR/FRT (n = 2) and subtypes F1 (n = 1) and C (n = 1). TDR prevalence was 58.3% (7/12). Isolates from two participants showed mutations associated with resistance to nucleoside reverse transcriptase inhibitors (NRTI) only (M41L, T125C, T125F, M184V), while an isolate from one patient who had received antiretroviral therapy (ART) since 2008 had a mutation associated with resistance to non-NRTI (G190S). Five isolates had secondary mutations to protease inhibitors (K20M, L10V, L33I, A71T, A71V). In conclusion, the findings of HIV-1 circulation, TDR to NRTI, and secondary mutations to protease inhibitors in ART-naïve crack cocaine users support the importance of monitoring this population in regions far from the epicenter of the HIV epidemic.

Viral diversity is an obligate consideration in CRISPR/Cas9 designs for targeting the HIV reservoir

BACKGROUND

RNA-guided CRISPR/Cas9 systems can be designed to mutate or excise the integrated HIV genome from latently infected cells and have therefore been proposed as a curative approach for HIV. However, most studies to date have focused on molecular clones with ideal target site recognition and do not account for target site variability observed within and between patients. For clinical success and broad applicability, guide RNA (gRNA) selection must account for circulating strain diversity and incorporate the within-host diversity of HIV.

RESULTS

We identified a set of gRNAs targeting HIV LTR, gag, and pol using publicly available sequences for these genes and ranked gRNAs according to global conservation across HIV-1 group M and within subtypes A-C. By considering paired and triplet combinations of gRNAs, we found triplet sets of target sites such that at least one of the gRNAs in the set was present in over 98% of all globally available sequences. We then selected 59 gRNAs from our list of highly conserved LTR target sites and evaluated in vitro activity using a loss-of-function LTR-GFP fusion reporter. We achieved efficient GFP knockdown with multiple gRNAs and found clustering of highly active gRNA target sites near the middle of the LTR. Using published deep-sequence data from HIV-infected patients, we found that globally conserved sites also had greater within-host target conservation. Lastly, we developed a mathematical model based on varying distributions of within-host HIV sequence diversity and enzyme efficacy. We used the model to estimate the number of doses required to deplete the latent reservoir and achieve functional cure thresholds. Our modeling results highlight the importance of within-host target site conservation. While increased doses may overcome low target cleavage efficiency, inadequate targeting of rare strains is predicted to lead to rebound upon cART cessation even with many doses.

CONCLUSIONS

Target site selection must account for global and within host viral genetic diversity. Globally conserved target sites are good starting points for design, but multiplexing is essential for depleting quasispecies and preventing viral load rebound upon therapy cessation.

Characterization of HIV-1 Near Full-Length Proviral Genome Quasispecies from Patients with Undetectable Viral Load Undergoing First-Line HAART Therapy

Increased access to highly active antiretroviral therapy (HAART) by human immunodeficiency virus postive (HIV⁺) individuals has become a reality worldwide. In Brazil, HAART currently reaches over half of HIV-infected subjects. In the context of a remarkable HIV-1 genetic variability, highly related variants, called quasispecies, are generated. HIV quasispecies generated during infection can influence virus persistence and pathogenicity, representing a challenge to treatment. However, the clinical relevance of minority quasispecies is still uncertain. In this study, we have determined the archived proviral sequences, viral subtype and drug resistance mutations from a cohort of HIV⁺ patients with undetectable viral load undergoing HAART as first-line therapy using next-generation sequencing for near full-length virus genome (NFLG) assembly. HIV-1 consensus sequences representing NFLG were obtained for eleven patients, while for another twelve varying genome coverage rates were obtained. Phylogenetic analysis showed the predominance of subtype B (83%; 19/23). Considering the minority variants, 18 patients carried archived virus harboring at least one mutation conferring antiretroviral resistance; for six patients, the mutations correlated with the current ARVs used. These data highlight the importance of monitoring HIV minority drug resistant variants and their clinical impact, to guide future regimen switches and improve HIV treatment success.

High prevalence of HIV-1 transmitted drug-resistance mutations from proviral DNA massively parallel sequencing data of therapy-naïve chronically infected Brazilian blood donors

BACKGROUND

An improved understanding of the prevalence of low-abundance transmitted drug-resistance mutations (TDRM) in therapy-naïve HIV-1-infected patients may help determine which patients are the best candidates for therapy. In this study, we aimed to obtain a comprehensive picture of the evolving HIV-1 TDRM across the massive parallel sequences (MPS) of the viral entire proviral genome in a well-characterized Brazilian blood donor naïve to antiretroviral drugs.

MATERIALS AND METHODS

The MPS data from 128 samples used in the analysis were sourced from Brazilian blood donors and were previously classified by less-sensitive (LS) or "detuned" enzyme immunoassay as non-recent or longstanding HIV-1 infections. The Stanford HIV Resistance Database (HIVDBv 6.2) and IAS-USA mutation lists were used to interpret the pattern of drug resistance. The minority variants with TDRM were identified using a threshold of ≥ 1.0% and ≤ 20% of the reads sequenced. The rate of TDRM in the MPS data of the proviral genome were compared with the corresponding published consensus sequences of their plasma viruses.

RESULTS

No TDRM were detected in the integrase or envelope regions. The overall prevalence of TDRM in the protease (PR) and reverse transcriptase (RT) regions of the HIV-1 pol gene was 44.5% (57/128), including any mutations to the nucleoside analogue reverse transcriptase inhibitors (NRTI) and non-nucleoside analogue reverse transcriptase inhibitors (NNRTI). Of the 57 subjects, 43 (75.4%) harbored a minority variant containing at least one clinically relevant TDRM. Among the 43 subjects, 33 (76.7%) had detectable minority resistant variants to NRTIs, 6 (13.9%) to NNRTIs, and 16 (37.2%) to PR inhibitors. The comparison of viral sequences in both sources, plasma and cells, would have detected 48 DNA provirus disclosed TDRM by MPS previously missed by plasma bulk analysis.

CONCLUSION

Our findings revealed a high prevalence of TDRM found in this group, as the use of MPS drastically increased the detection of these mutations. Sequencing proviral DNA provided additional information about TDRM, which may impact treatment decisions. The overall results emphasize the importance of continuous monitoring.

HIV-1 Genetic Diversity in Northeastern Brazil: High Prevalence of Non-B Subtypes

The Northeastern Brazilian region has experienced a constant increase in the number of newly reported AIDS cases over the last decade, but the genetic diversity of HIV-1 strains currently disseminated in this region remains poorly explored. HIV-1 pol sequences were obtained from 140 patients followed at outpatient clinics from four Northeastern Brazilian states (Alagoas, Bahia, Ceará, and Piauí) between 2014 and 2015. Subtype B was the most prevalent HIV-1 clade (72%) detected in the Northeastern region, followed by subtypes F1 (6%), C (5%), and D (1%). The remaining strains (16%) displayed a recombinant structure and were classified as follows: BF1 (11%), BC (4%), BCF1 (1%), and CRF02_AG like (1%). The 20 HIV-1 BF1 and BC recombinant sequences detected were distributed among 11 lineages classified as follows: CRF28/29_BF like (n = 5), CRF39_BF like (n = 1), URF_BF (n = 9), and URF_BC (n = 5). Non-B subtypes were detected in all Northeastern Brazilian states, but with variable prevalence, ranging from 16% in Ceará to 55% in Alagoas. Phylogenetic analyses support that subtype D and CRF02_AG strains detected in the Northeastern region resulted from the expansion of autochthonous transmission networks, rather than from exogenous introductions from other countries. These results reveal that HIV-1 epidemic spreading in the Northeastern Brazilian region comprised by multiple subtypes and recombinant strains and the molecular epidemiologic pattern in this Brazilian region is much more complex than originally estimated.

Characterization of HIV-1 CRF90_BF1 and putative novel CRFs_BF1 in Central West, North and Northeast Brazilian regions

The Brazilian AIDS epidemic has been characterized by an increasing rate of BF1 recombinants and so far eight circulating recombinant forms/CRFs_BF1 have been described countrywide. In this study, pol sequences (protease/PR, reverse transcriptase/RT) of 87 BF1 mosaic isolates identified among 828 patients living in six Brazilian States from three geographic regions (Central West, North, Northeast) were analyzed. Phylogenetic and bootscan analyses were performed to investigate the evolutionary relationship and mosaic structure of BF1 isolates. Those analyses showed that 20.7% of mosaics (18 out of 87) were CRFs-like isolates, mostly represented by CRF28/CRF29_BF-like viruses (14 out of 18). We also identified five highly supported clusters that together comprise 42 out of 87 (48.3%) BF1 sequences, each cluster containing at least five sequences sharing a similar mosaic structure, suggesting possible new unidentified CRFs_BF1. The divergence time of these five potential new CRFs_BF1 clusters was estimated using a Bayesian approach and indicate that they probably originated between the middle 1980s and the middle 1990s. DNA was extracted from whole blood and four overlapping fragments were amplified by PCR providing full/near full length genomes (FLG/NFLG) and partial genomes. Eleven HIV-1 isolates from Cluster # 5 identified in epidemiologically unlinked individuals living in Central West and North regions provided FLG/NFLG/partial genome sequences with identical mosaic structure. These viruses differ from any known CRF_BF1 reported to date and were named CRF90_BF1 by the Los Alamos National Laboratory. This is the 9^th CRF_BF1 described in Brazil and the first one identified in Central West and North regions. Our results highlight the importance of continued molecular screening and surveillance studies, especially of full genome sequences to understand the evolutionary dynamics of the HIV-1 epidemic in a country of continental dimensions as Brazil.

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.

Understanding the complex evolution of rapidly mutating viruses with deep sequencing: Beyond the analysis of viral diversity

With the advent of affordable deep sequencing technologies, detection of low frequency variants within genetically diverse viral populations can now be achieved with unprecedented depth and efficiency. The high-resolution data provided by next generation sequencing technologies is currently recognised as the gold standard in estimation of viral diversity. In the analysis of rapidly mutating viruses, longitudinal deep sequencing datasets from viral genomes during individual infection episodes, as well as at the epidemiological level during outbreaks, now allow for more sophisticated analyses such as statistical estimates of the impact of complex mutation patterns on the evolution of the viral populations both within and between hosts. These analyses are revealing more accurate descriptions of the evolutionary dynamics that underpin the rapid adaptation of these viruses to the host response, and to drug therapies. This review assesses recent developments in methods and provide informative research examples using deep sequencing data generated from rapidly mutating viruses infecting humans, particularly hepatitis C virus (HCV), human immunodeficiency virus (HIV), Ebola virus and influenza virus, to understand the evolution of viral genomes and to explore the relationship between viral mutations and the host adaptive immune response. Finally, we discuss limitations in current technologies, and future directions that take advantage of publically available large deep sequencing datasets.

High Prevalence and Onward Transmission of Non-Pandemic HIV-1 Subtype B Clades in Northern and Northeastern Brazilian Regions

The Human immunodeficiency virus type-1 (HIV-1) epidemic in Brazil is mainly driven by the subtype B pandemic lineage (BPANDEMIC), while Caribbean non-pandemic subtype B clades (BCAR) seem to account for a very low fraction of HIV-infections in this country. The molecular characteristics of the HIV-1 subtype B strains disseminated in the Northern and Northeastern Brazilian regions, however, have not been explored so far. In this study, we estimate the prevalence of the HIV-1 BPANDEMIC and BCAR clades across different Brazilian regions and we reconstruct the spatiotemporal dynamics of dissemination of the major Brazilian BCAR clades. A total of 2,682 HIV-1 subtype B pol sequences collected from 21 different Brazilian states from the five country regions between 1998 and 2013 were analyzed. Maximum Likelihood phylogenetic analyses revealed that the BCAR strains reached 16 out 21 Brazilian states here analyzed. The BCAR clades comprise a low fraction (<10%) of subtype B infections in most Brazilian states analyzed, with exception of Roraima (41%), Amazonas (14%) and Maranhão (14%). Bayesian phylogeographic analyses indicate that BCAR strains originally from the Hispaniola and Trinidad and Tobago were introduced at multiple times into different states from all Brazilian regions and a few of those strains, probably introduced into Roraima, Maranhão and São Paulo between the late 1970s and the early 1980s, established secondary outbreaks in the Brazilian population. These results support that the HIV-1 subtype B epidemics in some Brazilian states from the Northern and Northeastern regions display a unique molecular pattern characterized by the high prevalence of BCAR lineages, which probably reflects a strong epidemiological link with the HIV-1 epidemics in the Caribbean region.