Enhanced detection of viral diversity using partial and near full-length genomes of human immunodeficiency virus Type 1 provirus deep sequencing data from recently infected donors at four blood centers in Brazil

A validated in-house assay for HIV drug resistance mutation surveillance from dried blood spot specimens

Despite increasing scale-up of antiretroviral therapy (ART) coverage, challenges related to adherence and HIV drug resistance (HIVDR) remain. The high cost of HIVDR surveillance is a persistent challenge with implementation in resource-constrained settings. Dried blood spot (DBS) specimens have been demonstrated to be a feasible alternative to plasma or serum for HIVDR genotyping and are more suitable for lower resource settings. There is a need for affordable HIVDR genotyping assays which can amplify HIV-1 sequences from DBS specimens, particularly those with low viral loads, at a low cost. Here, we present an in-house assay capable of reliably amplifying HIV-1 protease and partial reverse transcriptase genes from DBS specimens, which covers the complete World Health Organization 2009 list of drug resistance mutations under surveillance. DBS specimens were prepared using whole blood spiked with HIV-1 at concentrations of 10,000, 5000, 1000, and 500 copies/mL (n=30 for each concentration). Specimens were tested in triplicate. A two-step approach was used consisting of cDNA synthesis followed by nested PCR. The limit of detection of the assay was calculated to be approximately 5000 (95% CI: 3200-10,700) copies/mL for the protease gene and 3600 (95% CI: 2200-10,000) copies/mL for reverse transcriptase. The assay was observed to be most sensitive with higher viral load specimens (97.8% [95% CI: 92.2-99.7]) for both protease and reverse transcriptase at 10,000 copies/mL with performance decreasing with the use of specimens with lower viral loads (46.7% [36.1-57.5] and 60.0% [49.1-70.2] at 500 copies/mL for protease and reverse transcriptase, respectively). Ultimately, this assay presents a promising opportunity for use in resource-constrained settings. Future work should involve validation under field conditions including sub-optimal storage conditions and preparation of DBS with fingerprick blood in order to accurately reflect real-world collection scenarios.

Identification of a HIV-1 circulating BF1 recombinant form (CRF75_BF1) of Brazilian origin that also circulates in Southwestern Europe

Introduction

The high recombinogenic potential of HIV-1 has resulted in the generation of countless unique recombinant forms (URFs) and around 120 reported circulating recombinant forms (CRFs). Here we identify through analyses of near full-length genomes (NFLG) a new HIV-1 CRF derived from subtypes B and F1.

Methods

HIV-1 protease-reverse transcriptase (Pr-RT) sequences were obtained by RT-PCR amplification from plasma RNA. Near full-length genome sequences were obtained after amplification by RT-PCR in 5 overlapping fragments. Phylogenetic sequence analyses were performed via maximum likelihood. Mosaic structures were analyzed by bootscanning and phylogenetic analyses of genome segments. Temporal and geographical estimations of clade emergence were performed with a Bayesian coalescent method.

Results

Through phylogenetic analyses of HIV-1 Pr-RT sequences obtained by us from samples collected in Spain and downloaded from databases, we identified a BF1 recombinant cluster segregating from previously reported CRFs comprising 52 viruses, most from Brazil (n = 26), Spain (n = 11), and Italy (n = 9). The analyses of NFLG genomes of 4 viruses of the cluster, 2 from Spain and 2 from Italy, allowed to identify a new CRF, designated CRF75_BF1, which exhibits a complex mosaic structure with 20 breakpoints. All 4 patients harboring CRF75_BF1 viruses studied by us had CD4+ T-cell lymphocyte counts below 220/mm3 less than one year after diagnosis, a proportion significantly higher (p = 0.0074) than the 29% found in other patients studied in Spain by us during the same period. The origin of the clade comprising CRF75_BF1 and related viruses was estimated around 1984 in Brazil, with subsequent introduction of CRF75_BF1 in Italy around 1992, and migration from Italy to Spain around 1999.

Conclusion

A new HIV-1 CRF, designated CRF75_BF1, has been identified. CRF75_BF1 is the 6th CRF of South American origin initially identified in Western Europe, reflecting the increasing relationship of South American and European HIV-1 epidemics. The finding of low CD4+ T-cell lymphocyte counts early after diagnosis in patients harboring CRF75_BF1 viruses warrants further investigation on the virulence of this variant.

HIV-1 subtype frequency in Northeast Brazil: A systematic review and meta-analysis

Information on human immunodeficiency virus (HIV) molecular epidemiology is required to verify HIV/AIDS (acquired immune deficiency syndrome) epidemic dynamics in different regions, as well as provide support for response to antiretroviral therapy, transmission of resistance mutations, disease progression, and viral spread. The aim of this study was to conduct a systematic review and meta-analysis of the frequency of HIV-1 subtypes in Northeast Brazil. Seventy-six articles that refer to HIV-1 and its subtypes in the Northeast Brazil and published between 1 January 1999 and 31 August 2019 were identified. We included 27 articles for the qualitative synthesis, thus analyzing results from 4466 patients and 4298 genomic sequences. The results showed that subtypes B, F, and C and recombinant BF were responsible for 76% (IC95%: 71-80), 8% (IC95%: 5-11), 2% (IC95%: 2-3), and 7% (IC95%: 4-12) infections, respectively. The highest proportion of subtype B infections (82.2%) was observed in Piauí, while the subtype F had a high frequency in Pernambuco (23.4%). Bahia presented 11.6% of the proportion of recombinant BF. In addition, several recombinants such as AG, BC, BCF, and BD have been identified in the region. This is the first systematic review and meta-analysis on the HIV-1 subtype distribution in Northeast Brazil and has shown a high circulating viral diversity. Although subtype B is predominant in Brazil, a large frequency of non-B subtypes has also been found, which may have consequences for response to antiretroviral therapy, disease progression, and transmission. Thus, HIV molecular epidemiological data are essential for epidemic prevention and control strategies.

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.

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.

HIV-1 Molecular Epidemiology, Transmission Clusters and Transmitted Drug Resistance Mutations in Central Brazil

We aimed to characterize HIV-1 molecular epidemiology and transmission clusters among heterosexual (HET) and men who have sex with men (MSM) individuals, as well as transmitted drug resistance mutations (TDRM) in Central-Western Brazil. This cross-sectional survey was conducted among 190 antiretroviral naïve HIV-1 infected individuals. Proviral DNA was extracted, and nested PCR amplified partial polymerase gene (PR/RT). After sequencing, subtypes were assigned, and the sequences were analyzed for the occurrence of possible transmission networks. Calibrated Population Resistance (CPR) tool from Stanford HIV Database was used to investigate the presence of TDRM. Among 150 individuals whose samples were successfully sequenced, the most prevalent HIV-1 subtype was B, followed by recombinant forms. The occurrence of twenty transmission clusters composed by at least two sequences was verified, suggesting the existence of transmission clusters among individuals from the same or distinct sexual orientations. Intermediate level of TDRM (12%) was found in the study population, and almost half of the subjects with TDRM had more than one resistance mutation. No correlations between sexual orientation and the presence of TDRM, HIV-1 subtypes/recombinants forms were verified. Taken together, the necessity of the continuous monitoring of the TDRM to verify the importance of pre-genotyping and to delineate future strategies in primary antiretroviral therapy. Likewise, the knowledge of the HIV-1 transmission networks in Brazil would allow the implementation of effective HIV-1 prevention strategies in local settings.

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.

Increase in human immunodeficiency virus 1 diversity and detection of various subtypes and recombinants in north-eastern Brazil

PURPOSE

Diverse human immunodeficiency virus 1 (HIV-1) subtypes and circulating recombinant forms are found in Brazil. The majority of HIV-1 molecular epidemiological studies in Brazil have been conducted in the southern and south-eastern regions of the country, although several recent studies in the north-eastern region have addressed this issue. The objective of this study was to molecularly characterize HIV-1 circulating in Pernambuco, north-eastern Brazil.

METHODOLOGY

A total of 64 samples were collected from 2002 to 2003, and another 103 were collected from 2007 to 2009. The protease and partial reverse transcriptase regions of the HIV-1 polymerase-encoding (pol) gene were sequenced, and subtyping, recombination and phylogenetic analyses were performed.Results/Key findings. Subtype B (60.9 %) was found to be predominant, followed by HIV-1 F (31.4 %). Several BF recombinants (4.2 %), and BC and AG recombinants were also identified. The intra-subtype genetic diversity was estimated to be 0.065 (sd±0.004) for HIV-1 B and 0.055 (sd±0.004) for HIV-1 F, reflecting a greater accumulation of mutations in subtype B (P<0.01). More codons were found to be under positive selective pressure in samples collected from 2007 to 2009, from individuals with a T-cell count≥200 cells mm-3 and from women. Coalescence data indicated that the subtype F population has been continuously expanding.

CONCLUSIONS

HIV-1 shows high genetic diversity in the state of Pernambuco. Thus, additional molecular evaluations of circulating strains will provide a better understanding of the epidemic and may lead to more effective preventive strategies.

Genetic Characterization of a Panel of Diverse HIV-1 Isolates at Seven International Sites

HIV-1 subtypes and drug resistance are routinely tested by many international surveillance groups. However, results from different sites often vary. A systematic comparison of results from multiple sites is needed to determine whether a standardized protocol is required for consistent and accurate data analysis. A panel of well-characterized HIV-1 isolates (N = 50) from the External Quality Assurance Program Oversight Laboratory (EQAPOL) was assembled for evaluation at seven international sites. This virus panel included seven subtypes, six circulating recombinant forms (CRFs), nine unique recombinant forms (URFs) and three group O viruses. Seven viruses contained 10 major drug resistance mutations (DRMs). HIV-1 isolates were prepared at a concentration of 107 copies/ml and compiled into blinded panels. Subtypes and DRMs were determined with partial or full pol gene sequences by conventional Sanger sequencing and/or Next Generation Sequencing (NGS). Subtype and DRM results were reported and decoded for comparison with full-length genome sequences generated by EQAPOL. The partial pol gene was amplified by RT-PCR and sequenced for 89.4%-100% of group M viruses at six sites. Subtyping results of majority of the viruses (83%-97.9%) were correctly determined for the partial pol sequences. All 10 major DRMs in seven isolates were detected at these six sites. The complete pol gene sequence was also obtained by NGS at one site. However, this method missed six group M viruses and sequences contained host chromosome fragments. Three group O viruses were only characterized with additional group O-specific RT-PCR primers employed by one site. These results indicate that PCR protocols and subtyping tools should be standardized to efficiently amplify diverse viruses and more consistently assign virus genotypes, which is critical for accurate global subtype and drug resistance surveillance. Targeted NGS analysis of partial pol sequences can serve as an alternative approach, especially for detection of low-abundance DRMs.

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.

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

Background

Here, we aimed to gain a comprehensive picture of the HIV-1 diversity in the northeast and southeast part of Brazil. To this end, a high-throughput sequencing-by-synthesis protocol and instrument were used to characterize the near full length (NFLG) and partial HIV-1 proviral genome in 259 HIV-1 infected blood donors at four major blood centers in Brazil: Pro-Sangue foundation (São Paulo state (SP), n 51), Hemominas foundation (Minas Gerais state (MG), n 41), Hemope foundation (Recife state (PE), n 96) and Hemorio blood bank (Rio de Janeiro (RJ), n 70).

Materials and Methods

A total of 259 blood samples were obtained from 195 donors with long-standing infections and 64 donors with a lack of stage information. DNA was extracted from the peripheral blood mononuclear cells (PBMCs) to amplify the HIV-1 NFLGs from five overlapping fragments. The amplicons were molecularly bar-coded, pooled, and sequenced by Illumina paired-end protocol.

Results

Of the 259 samples studied, 208 (80%) NFLGs and 49 (18.8%) partial fragments were de novo assembled into contiguous sequences and successfully subtyped. Of these 257 samples, 183 (71.2%) were pure subtypes consisting of clade B (n = 167, 65%), C (n = 10, 3.9%), F1 (n = 4, 1.5%), and D (n = 2, 0.7%). Recombinant viruses were detected in 74 (28.8%) samples and consist of unique BF1 (n = 41, 15.9%), BC (n = 7, 2.7%), BCF1 (n = 4, 1.5%), CF1 and CDK (n = 1, 0.4%, each), CRF70_BF1 (n = 4, 1.5%), CRF71_BF1 (n = 12, 4.7%), and CRF72_BF1 (n = 4, 1.5%). Evidence of dual infection was detected in four patients coinfected with the same subtype (n = 3) and distinct subtype (n = 1).

Conclusion

Based on this work, subtype B appears to be the prevalent subtype followed by a high proportion of intersubtype recombinants that appeared to be arising continually in this country. Our study represents the largest analysis of the viral NFLG ever undertaken worldwide and provides insights into the understanding the genesis of the HIV-1 epidemic in this particular area of South America and informs vaccine design and clinical trials.

Frequency of coreceptor tropism in PBMC samples from HIV-1 recently infected blood donors by massively parallel sequencing: the REDS II study

Background

The interaction of HIV-1 and target cells involves sequential binding of the viral gp120 Env protein to the CD4 receptor and a chemokine co-receptor (either CCR5 or CXCR4). CCR5 antagonists have proved to be an effective salvage therapy in patients with CCR5 using variants (R5) but not with variants capable of using CXCR4 (×4) phenotype. Thus, it is critically important to determine cellular tropism of a country’s circulating HIV strains to guide a management decision to improve treatment outcome. In this study, we report the prevalence of R5 and ×4 HIV strains in 45 proviral DNA massively parallel sequencing “MPS” data from recently infected Brazilian blood donors.

Methods

The MPS data encompassing the tropism-related V3 loop region of the HIV‐1 env gene was extracted from our recently published HIV-1 genomes sequenced by a paired-end protocol (Illumina). HIV‐1 tropism was inferred using Geno2pheno_[coreceptor] algorithm (3.5 % false-positive rate). V3 net charge and 11/25 rules were also used for coreceptor prediction.

Results

Among the 45 samples for which tropism were determined, 39 were exclusively R5 variants, 5 ×4 variants, and one dual-tropic or mixed (D/M) populations of R5 and ×4 viruses, corresponding to 86.7, 11.1 and 2.2 %, respectively. Thus, the proportion of all blood donors that harbor CXCR4-using virus was 13.3 % including individuals with D/M-tropic viruses.

Conclusions

The presence of CCR5-tropic variants in more than 85 % of our cohort of antiretroviral-naïve blood donors with recent HIV-1 infection indicates a potential benefit of CCR5 antagonists as a therapeutic option in Brazil. Therefore, determination of viral co-receptor tropism is an important diagnostic prerequisite.