Identification of an HIV-1 BG Intersubtype Recombinant Form (CRF73_BG), Partially Related to CRF14_BG, Which Is Circulating in Portugal and Spain

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.

Identification of CRF89_BF, a new member of an HIV-1 circulating BF intersubtype recombinant form family widely spread in South America

Circulating recombinant forms (CRFs) contribute substantially to the HIV-1 pandemic. Among 105 CRFs described in the literature, 16 are BF intersubtype recombinants, most of South American origin, of which CRF12_BF is the most widely spread. A BF recombinant cluster identified in Bolivia was suggested to represent a new CRF_BF. Here we find that it belongs to a larger cluster incorporating 39 viruses collected in 7 countries from 3 continents, 22 of them in Spain, most from Bolivian or Peruvian individuals, and 12 in South America (Bolivia, Argentina, and Peru). This BF cluster comprises three major subclusters, two associated with Bolivian and one with Peruvian individuals. Near full-length genome sequence analyses of nine viruses, collected in Spain, Bolivia, and Peru, revealed coincident BF mosaic structures, with 13 breakpoints, 6 and 7 of which coincided with CRF12_BF and CRF17_BF, respectively. In a phylogenetic tree, they grouped in a clade closely related to these CRFs, and more distantly to CRF38_BF and CRF44_BF, all circulating in South America. These results allowed to identify a new HIV-1 CRF, designated CRF89_BF. Through phylodynamic analyses, CRF89_BF emergence was estimated in Bolivia around 1986. CRF89_BF is the fifth CRF member of the HIV-1 recombinant family related to CRF12_BF.

Identification of a New HIV-1 BC Intersubtype Circulating Recombinant Form (CRF108_BC) in Spain

The extraordinary genetic variability of human immunodeficiency virus type 1 (HIV-1) group M has led to the identification of 10 subtypes, 102 circulating recombinant forms (CRFs) and numerous unique recombinant forms. Among CRFs, 11 derived from subtypes B and C have been identified in China, Brazil, and Italy. Here we identify a new HIV-1 CRF_BC in Northern Spain. Originally, a phylogenetic cluster of 15 viruses of subtype C in protease-reverse transcriptase was identified in an HIV-1 molecular surveillance study in Spain, most of them from individuals from the Basque Country and heterosexually transmitted. Analyses of near full-length genome sequences from six viruses from three cities revealed that they were BC recombinant with coincident mosaic structures different from known CRFs. This allowed the definition of a new HIV-1 CRF designated CRF108_BC, whose genome is predominantly of subtype C, with four short subtype B fragments. Phylogenetic analyses with database sequences supported a Brazilian ancestry of the parental subtype C strain. Coalescent Bayesian analyses estimated the most recent common ancestor of CRF108_BC in the city of Vitoria, Basque Country, around 2000. CRF108_BC is the first CRF_BC identified in Spain and the second in Europe, after CRF60_BC, both phylogenetically related to Brazilian subtype C strains.

Prevalence of Different Genotypes of HIV-1 in Injection Drug Users in China: A Systematic Review and Meta-Analysis

BACKGROUND

Since 1981, an increasing trend in HIV has been observed for transmission via injection drug users (IDUs), sexual transmission and mother-to-child transmission. The IDUs are blamed for early increases in HIV-positive cases in China.

OBJECTIVE

HIV genotypes of IDUs were comprehensively analysed to trace the source and relationships of the AIDS epidemic in China.

METHODS

Relevant databases written in English and Chinese were searched. Overall, 7,149 publications were identified in six databases. After screening 7,104 articles according to the inclusion and exclusion criteria, 45 studies consisting of 2,765 cases were finally identified. A meta-analysis was conducted using R MATLAB software, RevMan and SPSS. Subgroup analyses focused on time frame, region, and location of different genotypes of IDUs in China.

RESULTS

There were five dominant HIV-1 genotypes among the 2,765 IDU cases. The proportions of CRF07_BC, CRF01_AE, CRF08_BC, subtype B/B', and subtype C were 45.18% (95% CI: 33.55-57.08%), 16.00% (95% CI: 9.39-23.82%), 13.43% (95% CI: 7.32-20.84%), 3.58% (95% CI: 1.52-6.24%), and 0.90% (95% CI: 0.04-2.43%), respectively. HIV genotypes transmitted among IDUs in China are primarily CRF07-BC, followed by CRF01-AE and CRF08-BC. Across the different time frames and regions, CRF07_BC was the most prevalent HIV-1 genotype among IDUs, while CRF08_BC was the most prevalent genotype in the southwest region.

CONCLUSION

Our study reveals that CRF07-BC was the dominant prevalent strain among IDUs from 1991 to 2015 in China, while CRF08-BC was the dominant prevalent strain among IDUs in southwestern China. This systematic review and meta-analysis shows evidence of the comprehensive prevalence of different genotypes, data and characteristics of HIV among IDUs in China.

HIV-1 Unique Recombinant Forms Identified in Slovenia and Their Characterization by Near Full-Length Genome Sequencing

Surveillance of HIV circulating recombinant forms (CRFs) is important because HIV diversity can affect various aspects of HIV infection from prevention to diagnosis and patient management. A comprehensive collection of pol sequences obtained from individuals diagnosed with HIV-1 from 2000 to 2016 in Slovenia was subtyped to identify possible unique recombinant forms (URFs). Selected samples were subjected to near full-length genome (NFLG) sequencing and detailed recombination analyses. Discordant subtyping results were observed for 68/387 (17.6%) sequences and 20 sequences were identified as the most probable URFs and selected for NFLG characterization. Further, 11 NFLGs and two sequences of >7000 base pairs were obtained. Seven sequences were identified as “pure” subtypes or already characterized CRFs: subtype B (n = 5), sub-subtype A6 (n = 1), and CRF01_AE (n = 1). The remaining six sequences were determined to be URFs; four displayed a single recombination event and two exhibited a complex recombination pattern involving several subtypes or CRFs. Finally, three HIV strains were recognized as having epidemic potential and could be further characterized as new CRFs. Our study shows that the identification of new CRFs is possible, even in countries where HIV diversity is considered limited, emphasizing the importance of the surveillance of HIV recombinant forms.

Characterization of New Recombinant Forms of HIV-1 From the Comunitat Valenciana (Spain) by Phylogenetic Incongruence

Recombination is one of the main processes shaping the evolution of HIV-1, with relevant consequences for its epidemiology. In fact, Circulating and Unique Recombinant Forms (CRFs and URFs) cause 23% of current infections. The routine analyses of antiretroviral resistance yield partial pol gene sequences that can be exploited for molecular epidemiology surveillance but also to study viral diversity and to detect potential recombinant samples. Among the pol sequences derived from a large sample dataset from the Comunitat Valenciana (Spain), we identified nine putative recombinant samples. We aimed at fully characterizing these samples and performing a detailed analysis of the corresponding recombination events. We obtained nearly full-genome sequences and used jpHMM and RDP4 to detect and characterize recombinant fragments. We assessed the confidence of these inferences by likelihood mapping and phylogenetic placement with topology congruence tests. Next, we performed a phylogenetic analysis of each putative recombinant fragment to determine its relationships to previously described recombinant forms. We found that two samples related to CRF44_BF whereas the rest corresponded to new URFs (two URF_AD, one URF_BG that can constitute a new CRF resulting from subtype B and CRF24_BG, and two URF_cpx composed of A, G, K, H, and J subtypes). These URFs have a complex recombination pattern that cannot be determined accurately. They seem to have arisen by successive recombination events among lineages, including other CRFs. Our results highlight the usefulness of routine surveillance analysis for the detection of new HIV-1 recombination forms and, at the same time, the need for full-genome sequencing and recombination detection guidelines to properly characterize this complex process.

Molecular Epidemiology of HIV-1 Subtype G in the Russian Federation

Although HIV-1 subtype A has predominated in Russia since the end of the 20th century, other viral variants also circulate in this country. The dramatic outbreak of HIV-1 subtype G in 1988-1990 represents the origin of this variant spreading in Russia. However, full genome sequencing of the nosocomial viral variant and an analysis of the current circulating variants have not been conducted. We performed near full-length genome sequencing and phylogenetic and recombination analyses of 11 samples; the samples were determined to be subtype G based on an analysis of the pol region. Three samples were reliably obtained from patients infected during the nosocomial outbreak. The other 8 samples were obtained from patients who were diagnosed in 2010-2015. Phylogenetic analysis confirmed that a man from the Democratic Republic of the Congo was the origin of the outbreak. We also found that currently circulating viral variants that were genotyped as subtype G according to their pol region are in fact unique recombinant forms. These recombinant forms are similar to the BG-recombinants from Western Europe, particularly Spain and Portugal. The limitations of subtyping based on the pol region suggest that these viral variants are more widespread in Europe than is currently supposed.

Local Epidemics Gone Viral: Evolution and Diffusion of the Italian HIV-1 Recombinant Form CRF60_BC

The molecular epidemiology of HIV-1 in Italy is becoming increasingly complex, mainly due to the spread of non-B subtypes and the emergence of new recombinant forms. We previously characterized the outbreak of the first Italian circulating recombinant form (CRF60_BC), occurring among young MSM living in Apulia between the years 2009 and 2011. Here we show a 5-year follow-up surveillance to trace the evolution of CRF60_BC and to investigate its further spread in Italy. We collected additional sequences and clinical data from patients harboring CRF60_BC, enrolled at the Infectious Diseases Clinic of the University of Bari. In addition to the 24 previously identified sequences, we retrieved 27 CRF60_BC sequences from patients residing in Apulia, whose epidemiological and clinical features did not differ from those of the initial outbreak, i.e., the Italian origin, young age at HIV diagnosis (median: 24 years; range: 18–37), MSM risk factor (23/25, 92%) and recent infection (from 2008 to 2017). Sequence analysis revealed a growing overall nucleotide diversity, with few nucleotide changes that were fixed over time. Twenty-seven additional sequences were detected across Italy, spanning multiple distant regions. Using a BLAST search, we also identified a CRF60_BC sequence isolated in United Kingdom in 2013. Three patients harbored a unique second generation recombinant form in which CRF60_BC was one of the parental strains. Our data show that CRF60_BC gained epidemic importance, spreading among young MSM in multiple Italian regions and increasing its population size in few years, as the number of sequences identified so far has triplicated since our first report. The observed further divergence of CRF60_BC is likely due to evolutionary bottlenecks and host adaptation during transmission chains. Of note, we detected three second-generation recombinants, further supporting a widespread circulation of CRF60_BC and the increasing complexity of the HIV-1 epidemic in Italy.

A New HIV-1 Circulating Recombinant Form (CRF98_cpx) Between CRF06_cpx and Subtype B Identified in Southwestern France

During a recent study on the sequencing data of our database between 2012 and 2016 in Southwestern France, we observed that eight patients harbored what seemed to be the same virus. Indeed, routine genotyping at the time of HIV diagnosis showed that protease and reverse transcriptase were related to CRF06_cpx and subtype B, respectively. The integrase sequences (available for three patients) were clustering with CRF06_cpx and envelope (Env) gp120 sequences (available for two patients) with subtype B. Since such a recombinant has not been recorded in the Los Alamos database, we decided to characterize the full-length genome of this virus. The data suggest the identification of a new circulating recombinant form (CRF) between CRF06_cpx and subtype B, the structure of which is very complex with multiple breakpoints. We will refer this CRF as CRF98_cpx.

Trends in Drug Resistance Prevalence, HIV-1 Variants and Clinical Status in HIV-1-infected Pediatric Population in Madrid: 1993 to 2015 Analysis

BACKGROUND

The expanded use of long-term antiretroviral treatments in infected children may exacerbate the problem of drug resistance mutations selection, which can compromise treatment efficiency.

OBJECTIVE

We describe the temporal trends of HIV drug resistance mutations and the HIV-1 variants during 23 years (1993 to March 2016) in the Madrid cohort of HIV-infected children and adolescents.

METHODS

We selected patients with at least one available HIV-1 pol sequence/genotypic resistance profile, establishing different groups according to the sampling year of first resistance data. We determined the prevalence of transmitted drug resistance mutations or acquired drug resistance mutations (DRM), the drug susceptibility among resistant viruses and HIV-1 variants characterized by phylogeny across time.

RESULTS

A total of 245 pediatric patients were selected, being mainly female, Spanish native, perinatally infected and carrying HIV-1 subtype B. At first sampling, most pediatric patients were on antiretroviral therapy and heavily pretreated. During 1993 to 2016, transmitted drug resistance mutations was found in 13 (26%) of 50 naive children [non-nucleoside reverse transcriptase inhibitors (NNRTI), 14.6%; nucleoside reverse transcriptase inhibitors (NRTI), 10.4%; protease inhibitors, 8.7%]. DRM appeared in 139 (73.2%) of 190 pretreated patients (NRTI, 64.5%; NNRTI, 36%; protease inhibitors, 35.1%). DRM to NNRTI was higher in last 5 years. Non-B variants infected 14.5% of children and adolescents of the Madrid Cohort, being mainly intersubtype recombinants (76.5%), including complex unique recombinant strains. They caused 3.4% infections before 2000, rising to 85.7% during 2011 to 2016.

CONCLUSIONS

Periodic surveillance resistance and molecular epidemiology studies in long-term pretreated HIV-infected pediatric populations are required to optimize treatment regimens. Results will permit a better understanding of long-time dynamics of viral resistance and HIV-1 variants in Spain.

A high HIV-1 strain variability in London, UK, revealed by full-genome analysis: Results from the ICONIC project

BACKGROUND & METHODS

The ICONIC project has developed an automated high-throughput pipeline to generate HIV nearly full-length genomes (NFLG, i.e. from gag to nef) from next-generation sequencing (NGS) data. The pipeline was applied to 420 HIV samples collected at University College London Hospitals NHS Trust and Barts Health NHS Trust (London) and sequenced using an Illumina MiSeq at the Wellcome Trust Sanger Institute (Cambridge). Consensus genomes were generated and subtyped using COMET, and unique recombinants were studied with jpHMM and SimPlot. Maximum-likelihood phylogenetic trees were constructed using RAxML to identify transmission networks using the Cluster Picker.

RESULTS

The pipeline generated sequences of at least 1Kb of length (median = 7.46Kb, IQR = 4.01Kb) for 375 out of the 420 samples (89%), with 174 (46.4%) being NFLG. A total of 365 sequences (169 of them NFLG) corresponded to unique subjects and were included in the down-stream analyses. The most frequent HIV subtypes were B (n = 149, 40.8%) and C (n = 77, 21.1%) and the circulating recombinant form CRF02_AG (n = 32, 8.8%). We found 14 different CRFs (n = 66, 18.1%) and multiple URFs (n = 32, 8.8%) that involved recombination between 12 different subtypes/CRFs. The most frequent URFs were B/CRF01_AE (4 cases) and A1/D, B/C, and B/CRF02_AG (3 cases each). Most URFs (19/26, 73%) lacked breakpoints in the PR+RT pol region, rendering them undetectable if only that was sequenced. Twelve (37.5%) of the URFs could have emerged within the UK, whereas the rest were probably imported from sub-Saharan Africa, South East Asia and South America. For 2 URFs we found highly similar pol sequences circulating in the UK. We detected 31 phylogenetic clusters using the full dataset: 25 pairs (mostly subtypes B and C), 4 triplets and 2 quadruplets. Some of these were not consistent across different genes due to inter- and intra-subtype recombination. Clusters involved 70 sequences, 19.2% of the dataset.

CONCLUSIONS

The initial analysis of genome sequences detected substantial hidden variability in the London HIV epidemic. Analysing full genome sequences, as opposed to only PR+RT, identified previously undetected recombinants. It provided a more reliable description of CRFs (that would be otherwise misclassified) and transmission clusters.