Next Article in Journal
A Chimeric Sudan Virus-Like Particle Vaccine Candidate Produced by a Recombinant Baculovirus System Induces Specific Immune Responses in Mice and Horses
Next Article in Special Issue
The Impact of Human Mobility on Regional and Global Efforts to Control HIV Transmission
Previous Article in Journal
Prevalence and Genetic Diversity of Rotaviruses among under-Five Children in Ethiopia: A Systematic Review and Meta-Analysis
Previous Article in Special Issue
Exploring HIV-1 Transmission Dynamics by Combining Phylogenetic Analysis and Infection Timing
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 1
10.3390/v12010063
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Maja M. Lunar
,
Jana Mlakar
,
Tomaž Mark Zorec
and
Mario Poljak
*
Institute of Microbiology and Immunology, Faculty of Medicine, University of Ljubljana, 1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Viruses 2020, 12(1), 63; https://0-doi-org.brum.beds.ac.uk/10.3390/v12010063
Submission received: 3 December 2019 / Revised: 28 December 2019 / Accepted: 31 December 2019 / Published: 3 January 2020
(This article belongs to the Special Issue HIV Molecular Epidemiology for Prevention)
Browse Figures
Versions Notes

Abstract

:
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.

1. Introduction

Even though a successful antiretroviral treatment (ART) has been available for over two decades, human immunodeficiency virus (HIV) remains one of the major global disease burdens, causing 770,000 deaths in 2018 [1]. New HIV infections have decreased considerably since the peak in 1997, from 2.9 million to 1.7 million in 2018; however, we are still far from ending the epidemic and all measures need to be addressed to do so [1]. Two major hurdles are the viruses ability to integrate into the human genome and its high genetic diversity, hampering the extensive efforts to design a drug to eradicate HIV from its reservoirs and to design an effective prophylactic vaccine. Even if a vaccine becomes available and adopted, continuous surveillance of HIV subtypes and recombinant forms will remain crucial [2].
The major pandemic HIV-1 group M is represented by 10 subtypes (A–D, F–H, J, K, and the recently established subtype L) and 98 characterized circulating recombinant forms (CRFs), found in at least three epidemiologically unlinked individuals [3,4,5]. Some of the CRFs are so-called second-generation recombinants (SGRs), composed of two or more previously characterized CRFs [6]. A recent systematic review of all published and unpublished HIV-1 subtyping data collected between 1990 and 2015 showed that subtype C is currently responsible for almost half of HIV-1 infections globally, followed by subtype B and subtype A. The global prevalence of all recombinant forms was estimated at 22.8% between 2010 and 2015, with CRFs accounting for 16.7% and unique recombinant forms (URFs) for 6.1% of all global HIV-1 infections. A steady increase in the prevalence of recombinants over time has been noted [7]. HIV diversity can affect various aspects of HIV infection, from transmission and pathogenesis to diagnostic tests and viral load measurements, ART effectiveness and drug resistance, immune response, and vaccine development [6]. Thus, surveillance of circulating recombinant forms is considered essential.
The Slovenian HIV-1 epidemic is mainly driven by subtype B; however, 22.2%, 10.7%, and 17.3% of non-B subtypes were identified in Slovenia in the 2000–2004, 2005–2010, and 2011–2016 study periods, respectively [8,9,10]. In addition, 5.2% of CRFs (02_AG) and 7.8% of URFs in the 2000–2004 study period, 0.7% of CRFs and 2.0% of URFs in 2005–2010 period, and 2.4% of CRFs and 3.6% of URFs in the 2011–2016 study period were observed among the subtyped sequences [8,9,10].
In this study, we have characterized potential HIV-1 URFs identified in Slovenia by obtaining near full-length genome (NFLG) sequences and performing detailed recombination analyses.

2. Materials and Methods

A total of 634 persons were diagnosed with HIV-1 in Slovenia from 2000 to 2016. Among them, 387 (61.0%) partial pol sequences were obtained previously to study the prevalence of pre-treatment HIV-1 drug resistance in Slovenia as part of the European Commission—sponsored SPREAD program [8,9,10], approved by the national Medical Ethics Committee at the Slovenian Ministry of Health (approval reference number 126/12/03). This study was conducted in accordance with the Code of Ethics of the World Medical Association (Declaration of Helsinki) and all the patients gave consent to participate in the study.
Sequences were analyzed for evidence of recombination with the following subtyping tools: Rega HIV-1 & 2 automated subtyping tool, version 2.0 (Rega 2) [11], Rega HIV-1 subtyping tool, version 3.0 (Rega 3) [12], COMET HIV-1, version 1.0 (Comet) [13], jpHMM [14], SCUEAL [15], Geno2pheno [16], and Stanford HIVdb Program [17].
Next, phylogenetic analysis was performed for further inspection of Slovenian sequences that yielded at least one discordant subtyping result with the seven subtyping tools. These sequences were aligned together with the Compendium alignment which was obtained from the Los Alamos database (all subtypes and CRFs of group M available at that time and three sequences of group O to be used as the outgroup) and a set of the most similar sequences to serve as controls [18]. Control sequences were obtained from a Los Alamos database blast search, in which 10 sequences per each Slovenian sequence were downloaded [18,19]. Duplicated sequences were removed from the dataset and aligned by MUSCLE available in Mega version X [20]. Maximum likelihood phylogeny was inferred by employing PhyML with automatic model selection according to the Akaike information criterion [21]. Nearest neighbor interchange was employed as a tree improvement method and aLRT (approximate likelihood ratio test) was used to assess the significance of transmission clusters. The phylogenetic tree was visualized in Figtree, version 1.4.3 [22].
Samples were selected for further NFLG analyses according to the results of the subtyping tools and the phylogenetic analysis; specifically, if the sequences were positioned: (i) at the outskirts of subtype clusters or (ii) alongside previously determined URFs on the phylogenetic tree.
NFLG sequences were obtained according to the previously published protocol [23]. Briefly, four overlapping PCR amplicons were obtained and joined in equimolar amounts by using a Qubit dsDNA HS Assay kit (Thermo Fisher Scientific, Waltham, MA, USA). A library was prepared using Nextera DNA Flex and sequenced on an Illumina MiSeq platform (Illumina, San Diego, CA, USA). The FastQ files that were obtained were assembled with a SHIVER pipeline [24]. Briefly, contigs were assembled de novo from paired-end reads using IVA [25] and compared to the reference set of HIV sequences, Compendium HIV-1 genome alignment (2017), which were retrieved from the Los Alamos database using BLASTN [18,19]. Contigs representing HIV sequences were corrected and aligned to the references for visual inspection. Then, a mapping reference was created from the contigs and the gaps between the contigs were filled with the closest reference sequence. Reads were trimmed of primer and adapter sequences and of low-quality reads using Trimmomatic [26] and Fastaq [27] and then mapped to the compiled reference using Smalt [28]. Finally, a consensus was called using SAMtools [29].
The NFLG sequences that were obtained were first analyzed with the following subtyping tools: Rega 3, Comet, and jpHMM. SCUEAL can only be used to determine the subtype for pol sequences and was therefore not employed. Similarly, the primary aim of Geno2pheno and the Stanford HIVdb Program is genotypic HIV drug resistance interpretation and they are therefore designed to analyze genome regions targeted by antiretroviral drugs and were also not employed.
Further recombination analyses were carried out using SimPlot 3.5.1 [30]. The NFLG sequences that were obtained were first examined by performing exploratory similarity plot analyses available in SimPlot, with all HIV subtype reference sequences available from the SimPlot 3.5.1 Los Alamos HIV database [18]. Next, the bootscan method was employed with query sequences corresponding to subtypes A–D, F–H, J, K, and CRF01_AE or with a corresponding subset of subtype references. The bootscan procedure was run with 500 bootstrap replicates and a 400 base-pair (bp) sliding window with a 10 bp step size, without gap stripping. If the sequence was suspected to be a secondary recombinant (with recombination among CRFs), CRFs that were identified in a similarity plot were also included in the bootscanning. Finally, recombination breakpoints were identified and slices of the alignments were saved for further phylogenetic analyses, performed as described above. Final unique recombinant forms were visualized with the Recombinant HIV-1 Drawing Tool [31].
An additional BLAST search of NFLG sequences determined to be URFs was employed to search for sequences with the highest identity based on the pol region. The BLAST hits that were obtained were further analyzed with the selected subtyping tools, which displayed a recombination signal when query URF sequences were first analyzed [19]. If two or more sequences were identified with a similar recombination pattern, the URF was identified as a candidate for a new CRF.

3. Results

3.1. Subtyping of Pol Region

The subtyping tools yielded at least one discordant result in 68/387 (17.6%) partial pol sequences obtained from individuals with an HIV-1 diagnosis from 2000 to 2016 in Slovenia (Table 1). The most discordant results were obtained with Rega 2 (n = 32), followed by Rega 3 (n = 27), Geno2pheno (n = 23), SCUEAL (n = 21), Stanford (n = 21), jpHMM (n = 12), and Comet (n = 11). On the other hand, the subtyping tool that gave an indication of a possible URF in the most cases was, again, Rega 2 (n = 28), followed by Rega 3 (n = 22), SCUEAL (n = 21), Stanford (n = 20), jpHMM (n = 13), and Comet (n = 11). Geno2pheno did not report any cases of possible recombination (Table 1).
All 68 potential recombinant sequences were included in the phylogenetic analysis together with 365 control sequences, all available reference sequences of HIV-1 group M, and three sequences of HIV-1 group O as the outgroup. The alignment was carried out on 614 sequences altogether (positioning according to the HXB2 reference sequence: 2268–3269 bp) with a length of 1002 bp, and phylogeny was inferred (Figure 1). According to the results of the subtyping tools and the pol phylogenetic analysis, 20 sequences were identified as the most probable unique recombinant forms and were selected for further NFLG characterization (Figure 1, Table 1).

3.2. Near Full-Length HIV Genome Sequencing and Characterization

We successfully obtained all four overlapping amplicons for 13/20 (65%) of the samples, which were afterwards sequenced using the Illumina MiSeq platform. Due to low read coverage of the fourth PCR amplicon, the expected length of the NFLG sequence was not obtained in two cases and ended at nucleotide positions HXB2: 7842 bp and HXB2: 7843 bp, respectively.
The NFLG sequences obtained were first analyzed with Rega 3, Comet, and jpHMM (Table 2). Second, an exploratory similarity plot analysis and bootscanning, available in SimPlot, were performed (Table 2, Figure 2 and Figure S1). Recombination breakpoints were identified and slices of the alignments were saved for further phylogenetic analyses. Subtypes/CRFs were assigned and possible unique recombinants were characterized in seven and six cases, respectively (Table 2, Figure 3).
To observe the sensitivity and specificity of subtyping tools based on the pol region, we inspected the subtyping results among the six sequences determined to be URFs. SCUEAL performed best, with five results indicative of recombination, followed by Rega 2 (n = 4) and jpHMM (n = 4). Rega 3, Comet, and Stanford all detected recombination in half of the cases, whereas Geno2pheno did not detect them at all. On the other hand, a false-positive recombination signal in the sequences, characterized as non-recombinant by NFLG analysis, was detected with Rega 3 (5/7), Rega 2 (4/7), and SCUEAL (1/7).
No clinically relevant drug resistance mutations were identified in the obtained NFLG sequences by the Stanford HIVdb Program.

3.2.1. Sample 11

The sequence of 7320 bp obtained was determined to be subtype A1/K recombinant with jpHMM, whereas Comet and Rega 3 assigned subtype A1 to it. SimPlot analysis characterized it as A1/K recombinant or CRF45_cpx (Table 2, Figure 2). Phylogenetic analysis was performed on nine separate slices of alignment and displayed clustering with subtype A in 7/9 occurrences with aLRT values of >0.85 (Figure S2). Phylogenetically, the closest sequence was CRF45_cpx in 3/9 segments, followed by subtype A1, CRF09_cpx, and CRF63_02A in one segment each. The remaining three alignments of the query sequence did not cluster within transmission pairs. Alignment slice HXB2: 5178–5436 bp clustered within the subtype A cluster, but could not be distinguished from CRF01_AE, and thus, this region was designated as an area of less certainty (A/unclassified (U)). Slice HXB2: 6031–6367 bp clustered with subtype K sequences (aLRT = 0.93). Both partitions clustered in a transmission pair with CRF45_cpx (Figure S2). Finally, Sample 11 was characterized as a URF of subtypes A and K with an area of uncertainty at position HXB2: 5307 ± 129 bp (Figure 3).
The sample was obtained in 2008 from an individual who reported that he most likely acquired HIV from a male of African origin residing in Slovenia. A BLAST search of the pol region revealed another similar Slovenian sequence without a known epidemiological link and several sequences from Nigeria and Cameroon. It is apparent that this variant is composed of subtype K and of a divergent subtype A1, as shown in Figure 1, which is possibly related to the ancestor of CRF01_AE. A highly significant phylogenetic cluster (aLRT = 0.99) encompassing the Sample 11 sequence and several sequences from Africa was observed on the pol phylogenetic tree. Thus, this variant could be a candidate CRF when two additional NFLGs obtained from unlinked individuals are sequenced.

3.2.2. Sample 12

We obtained a sequence of 9014 bp and the three subtyping tools determined it to be of subtype B (Table 2). This was confirmed by the detailed SimPlot analysis; however, two regions exhibiting less significant bootstrap support with subtype B were identified (HXB2: 3797–4218 bp and 8882–9183 bp; Figure S1). Further phylogenetic analyses of the two segments confirmed clustering within subtype B sequences with aLRT of 0.73 and 0.96, respectively.

3.2.3. Sample 25

We successfully sequenced 9053 bp of the HIV genome. The three subtyping tools determined this sequence to be of sub-subtype A1; however, SimPlot analyses determined it to be of sub-subtype A6 (Table 2, Figure S1). The sample was obtained from a person originating from Ukraine. Three additional individuals were identified with this subtype in Slovenia, according to the pol phylogenetic tree.

3.2.4. Sample 29

We sequenced 8964 bp of the genome and assigned it subtype B by the subtyping tools as well as SimPlot analyses (Table 2, Figure S1).

3.2.5. Sample 31

Only 7308 bp were obtained due to poor read coverage of the fourth PCR amplicon. jpHMM and Comet determined the sequence to be a complex recombinant, whereas Rega 3 assigned subtype C to it. SimPlot analysis suggested it to be CRF60_BC with an area of poor bootstrap support of the segment HXB2: 4811–5293 bp (Figure 2). This region was identified as CRF01_AE by phylogenetic analysis with aLRT of 0.76 (Figure S2). The surrounding two regions were confirmed as CRF60_BC with aLRT of 1.00 and therefore, this sequence was proposed to be a URF consisting of CRF60_BC and CRF01_AE (Table 2, Figure 3).
A BLAST search identified sequences of the recently characterized Italian CRF60_BC recombinant as sequences with the highest identity scores [32]. The phylogenetic tree of the pol region displayed a Sample 11 sequence on the root of the cluster of these CRF60_BC sequences from Italy, with aLRT = 1.00. Thus far, we do not have data for additional samples with this CRF60_BC and CRF01_AE recombinant and therefore, we do not have any indications that this variant is still circulating and we propose that it is a unique recombinant.

3.2.6. Sample 34

We sequenced 9030 bp of the genome and all three subtyping tools identified it as a complex recombinant (Table 2). This was also evidenced in the further SimPlot analysis, which indicated a recombinant of subtypes A, G, J, C, and CRF01_AE (Figure 2). Namely, 17 separate slices of alignment were saved according to bootscanning and corresponding phylogenetic analyses were performed (Figure S2). A subtype could not be assigned to seven segments; however, the remaining segments were assigned subtypes A1, G, J, and CRF01_AE. CRF13_cpx was found to be most closely related to this sequence on 5/17 phylogenetic trees, followed by CRF11_cpx (3/17; Figure S2). Finally, a complex URF of subtypes A, G, J, CRF01_AE, and U was described, as shown in Table 2 and with HXB2 breakpoints in Figure 3.
The individual originates from Gambia and a BLAST search revealed the sequence to be most similar to another Slovenian sequence obtained from a person with an epidemiological link, as we know from the patient’s records. Other similar sequences were from Liberia and Nigeria; however, those sequences did not show the same recombination signal when the pol region was analyzed and therefore, we do not have data that more sequences are circulating and we determined that this variant is a URF and not a CRF candidate.

3.2.7. Sample 37

The NFLG sequence of 8979 bp displayed a complex recombination pattern with all three subtyping tools (Table 2). SimPlot proposed a recombinant of CRF56_cpx and subtype B. Five different segments of the alignment were generated and phylogenetic analysis confirmed CRF56_cpx as the most closely related sequence in 4/5 segments, whereas subtype B was determined to be the origin of the remaining segment of the genome (HXB2: 7853–8362 bp; Figure S2). CRF56_cpx clustered with 02_AG sequences on the phylogenetic tree of this segment (Figure S2).
A BLAST search identified several sequences from Turkey, displaying high identity. This was also observed on the pol phylogenetic tree, where the query sequence was in a transmission cluster with Turkish sequences (aLRT = 0.99). Comet subtyped the most similar sequence as an unassigned sequence composed of CRF02_AG, CRF56_cpx, and subtype B, and therefore it is possible that this recombinant is in circulation and could thus be a CRF candidate.
In addition, we obtained this sample in 2012 from a person who reported that he most probably acquired HIV in Turkey.

3.2.8. Sample 39

Comet detected possible recombination between CRF01_AE, subtype B, and subtype D in the sequence of 9048 bp; however, Rega 3 and jpHMM identified it as CRF01_AE. SimPlot analysis showed high bootstrap support for CRF01_AE throughout the sequence and thus, we concluded that this sequence is CRF01_AE (Figure S1).

3.2.9. Sample 43

Only Comet identified this sequence of 8910 bp as a potential recombinant of subtypes B, C, and G, whereas jpHMM and Rega 3 determined it to be subtype B (Table 2). SimPlot, on the other hand, displayed three regions of low bootstrap support with subtype B and indicated a recombinant of subtype D (Figure S1). Two of the segments were clearly within the subtype B cluster on the phylogenetic tree. For the last segment (HXB2: 8929–9412 bp), the subtype could not be determined, possibly due to the high similarity between subtype B and subtype D because subtype D sequences were nested within the subtype B cluster. We could not confirm any recombination event in this sequence and determined it to be of subtype B (Table 2).

3.2.10. Sample 54

Rega 3 and Comet subtyped the sequence of 8970 bp as CRF02_AG and jpHMM as subtype A1/G recombinant (Table 2). Bootscan identified it as CRF02_AG with two sections of lower bootstrap support, which were further analyzed phylogenetically (Figure 2 and Figure S2). Segment HXB2: 6580–6971 bp clustered with sub-subtype A1 sequences (aLRT = 0.87), distinctly separated from the CRF02_AG cluster (Figure S2). This sequence was therefore characterized as CRF02_AG and subtype A recombinant (Figure 3).
The sequence with the highest identity retrieved from a BLAST search was a sequence from Cameroon (MH705136) and was subtyped as CRF02_AG. Therefore, we do not have any evidence to suggest that this variant is still in circulation and so we propose that it is a URF.

3.2.11. Sample 56

A sequence of 9010 bp was obtained and identified as subtype B by the three subtyping tools as well as SimPlot analysis (Table 2, Figure S1).

3.2.12. Sample 65

The sequence of 9039 bp obtained was determined to be subtype B by the subtyping tools and SimPlot (Table 2, Figure S1).

3.2.13. Sample 68

The NFLG sequence of 9048 bp was recognized as a complex recombinant by the three subtyping tools. Bootscan and SimPlot displayed the highest similarity with CRF19_cpx and subtype B/G recombinants (the group of CRF20_BG, CRF23_BG, and CRF24_BG). Separate phylogenetic analyses were performed on 15 sections of the alignment (Figure S2). Query sequence was found in a transmission pair with CRF19_cpx in 5/15 segments, in a cluster with CRF19_cpx and CRF37_cpx in two segments, and as a parental strain to CRF19_cpx in another segment. In another 5/15 segments, the sequence was observed in a phylogenetic cluster with subtype B/G recombinants, mostly in a close relationship with CRF20_BG (3/5). In one segment, phylogenetic analysis determined a close relationship with sub-subtype A1, and in two other segments, a decision could not be made and it was therefore determined to be U (Figure S2). Because we could not assign a clear BG recombinant for the non-CRF19_cpx regions of the genome, we propose that the recombinant is probably between CRF19_cpx and a still-undescribed subtype B/G/A recombinant and we have characterized the sequence as a recombinant of CRF19_cpx, subtype G, subtype B, subtype A, and U (Table 2, Figure 3).
The sequence was obtained in 2016 from a person of Slovenian origin who reported that he probably acquired HIV in Slovenia. A BLAST search showed the sequence to be most similar to two Slovenian sequences. The first was collected in 2008 from an individual who reported that he acquired the HIV infection in Italy. The second was collected in 2016. For now, we do not have data on whether these three sequences were obtained from epidemiologically linked individuals, apart from them being men who have sex with men (MSM). Another sequence exhibiting high identity was a sequence collected in 2008 in Spain (GU830956) and submitted under an unpublished study called “Identification of a new HIV-1 circulating recombinant form in Spain.” These four sequences formed a phylogenetic cluster (aLRT = 0.996) nested within the CRF19_cpx cluster on the phylogenetic tree of the pol region (Figure 1). Thus, this variant seems to have spread between countries, possibly continents, and it is therefore a good candidate for a new CRF.

4. Discussion

A comprehensive collection of sequences obtained from individuals diagnosed with HIV-1 from 2000 to 2016 in Slovenia was subtyped to study possible URFs. Discordant subtyping results were observed for 18% (68/387) of sequences and 20 samples were selected for NFLG sequencing. Eleven complete NFLGs and an additional two sequences of over 7000 bp of HIV-1 possible unique recombinants were obtained and characterized in this study. 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 unique recombinants. Among these, four displayed a single recombination event between either “pure” subtypes or CRFs and three were so-called SGRs (CRF60_BC|CRF01_AE|CRF60_BC, CRF56_cpx|B|CRF56_cpx, and CRF02_AG|A|CRF02_AG). On the other hand, two sequences exhibited a complex recombination pattern involving several subtypes or CRFs as well as regions of unknown origin. A search for additional sequences with a similar recombination profile revealed that additional samples of unlinked individuals could be sequenced for three of the URFs identified and could therefore lead to the characterization of new CRFs.
The HIV epidemic in Slovenia is mainly driven by subtype B; therefore, it is not surprising that 25 of the 68 divergently subtyped pol sequences were subtyped as potential recombinants of subtype B by the subtyping tools. Five of the samples were subsequently sequenced and no apparent recombination was determined in the NFLG sequences—all sequences were determined to be subtype B.
An additional two samples were also subsequently not characterized as unique recombinants (CRF01_AE and sub-subtype A6). Interestingly, Sample 25 was chosen because it clustered in a divergent A1 clade on the pol phylogenetic tree, together with the sub-subtype A6 reference sequence. It had a discordant subtyping score because Stanford and Geno2pheno both subtyped it as CRF01_AE. However, all of the subtyping tools designed specifically for subtyping characterized it as sub-subtype A1 because sub-subtype A6 is not included in their reference alignments. Thus, none of the subtyping tools were able to correctly assign sub-subtype A6. This sub-subtype is still the predominant sub-subtype in Russia, hence its former name subtype AFSU [33]. Sub-subtype A6 was shown to have the most recent common ancestor (tMRCA) around 1984 in the city of Odessa, Ukraine, from where it most likely spread to other cities of Ukraine and other FSU countries [34]. We obtained Sample 25 from a person who most probably acquired HIV in Ukraine. Forward transmission in Slovenia was observed, with three additional sequences identified.
We employed a subsequent BLAST search of those sequences recognized as potential URFs to look for other related sequences. It was carried out on the pol region only because of the most comprehensive database and it is therefore possible that some similar sequences were overlooked. We identified Sample 68 as having the most potential for a new CRF. This variant was detected in three individuals from Slovenia and it could have originated in Italy, based on an individual report of one patient. Another highly similar sequence was obtained in Spain. We determined it to be a recombinant of CRF19_cpx and most probably of a still-undescribed B/G/A recombinant, similar to that of CRF20_BG. We can speculate that the recombination occurred in Cuba and spread from there, based on the origin of CRF19_cpx, CRF20_BG, CRF23_BG, and CRF24_BG. A high prevalence of CRF19_cpx and CRF20, 23, 24_BG was observed among HIV-positive Cuban adults in a recent study: 24% and 23%, respectively [35]. This makes recombination events more plausible. On the other hand, CRF19_cpx spread was observed in southern Spain, where several transmission networks of young native Spanish MSM were determined, with tMRCAs in the period from 2007 to 2010 [36,37,38]. This, coupled with the recent characterization of new recombinant forms of subtypes B and G described in Spain and Portugal, could suggest that the URF identified in this study might also have originated in Europe [39,40].
Worryingly, CRF19_cpx has been associated with rapid disease progression among Cuban patients. Patients had a higher HIV viral load at diagnosis, they presented with AIDS within three years, and the variant was observed to be more frequently CXCR4-tropic and had better fitness of the protease [41]. We do not have data regarding disease progression for all three Slovenian individuals infected with this strain because two were diagnosed in 2016 and started treatment shortly after diagnosis, following the European HIV treatment guidelines [42]. On the other hand, the person diagnosed in 2008 was certain to have acquired HIV infection in Italy in 2007. He had a CD4 cell count of 445 cells/mm3 at diagnosis, three months later his CD4 cell count dropped to 351 cells/mm3, and afterward, antiretroviral treatment was initiated. This could indicate rapid disease progression, as mentioned above. Indeed, the protease region of Sample 68 partially consisted of CRF19_cpx.
Another potential candidate for CRF characterization is Sample 11. Five of the nine genome segments were observed in a close phylogenetic relationship with CRF09_cpx and/or CRF45_cpx. A study examining the evolutionary history of these complex recombinants found that CRF09_cpx and CRF45_cpx were positioned at the root of subtype A/A1 on the phylogenetic tree of the env region. This indicates that these two CRFs originate from variants that diverged from the ancestral lineage at a similar time as the ancestor of subtype A. Thus, the tMRCA of the env region was estimated to be around 1965 and 1960 for CRF09_cpx and CRF45_cpx, respectively. The tMRCA of the pol region was estimated to be a decade later [43]. The sequence of Sample 11 shows a less complex structure and yet the subtype A part of the sequence, which corresponds to the env part of that study, shows a high similarity to that of CRF09_cpx and CRF45_cpx. Similarly, we observed that these sequences clustered basally at the root of CRF01_AE on the pol phylogenetic tree. Interestingly, a NFLG sequence collected in 2010 in Nigeria (KX389608) was found among the clustered sequences. It was subtyped as subtype A by the authors; however, our jpHMM analysis displayed a nearly identical recombination profile to that of the Sample 11 sequence, with additional complexity at the 3′-end, which we failed to sequence. We obtained the Sample 11 sequence from a person who probably acquired the infection from an HIV-positive person of African origin and we identified another similar sequence from Slovenia and several African sequences. We can speculate that the variant is an old lineage that is still in circulation and it could be one of the progenitors of other CRFs.
Sample 37 was characterized as a CRF56_cpx|subtype B|CRF56_cpx recombinant. We identified several sequences with high similarity from Turkey, where the person whose sample was obtained also reported having been infected with this particular strain. Interestingly, CRF56_cpx was first described among young MSM from France. The tMRCA of around 2007 was estimated for this CRF, suggesting a recent diffusion. The authors described related strains found in Cyprus and Greece, and they proposed that they probably share a common ancestor [44]. This falls in line with our findings, implying that the origin of these strains could have been in the eastern Mediterranean region. In addition, CRF56_cpx has been observed in Pakistan and the Philippines, thus reaching global diffusion [45].
Sample 31 was recognized as a recombinant of CRF60_BC and CRF01_AE. The sequence was obtained from a person who reported having acquired HIV in Slovenia. We do not have further evidence that this variant has epidemic potential. However, a few SGRs composed of CRF60_BC have been described in Italy, where the CRF60_BC epidemic grew exponentially between 2006 and 2011 [32]. Because the CRF01_AE section of the Sample 31 sequence is outside the routinely genotyped pol region for resistance detection, further detection of this URF will be more challenging.
A complex recombination pattern was observed in the Sample 34 sequence, probably originating from Gambia. We could not assign a subtype to several regions of the sequence that was obtained. This indicates that these parts of the sequence are too divergent from the known subtypes and therefore this could represent a recombinant of an ancestral lineage or a lineage that subsequently died out. In fact, the sequence was found to be related in parts to CRF13_cpx and CRF11_cpx. These CRFs were shown to probably originate from Central Africa, where they emerged in the late 1950s to late 1960s [43]. Tongo et al. (2015) proposed that parts of genomes of the unique recombinant sequences from Cameroon that clustered outside of all identified subtypes and CRF clades are in fact historical lineages that diverged early in the HIV epidemic but did not reach the magnitude of an epidemic like other subtypes did [46].
The main limitation of this study is the selection of samples for URF characterization, which was based on subtyping and phylogenetic analysis of the pol region only, therefore possibly missing URFs that do not display recombination breakpoints within this region. However, the pol region was employed because the sequences were available from previous studies analyzing transmitted drug resistance in Slovenia [8,9,10]. Most Western European countries perform routine baseline drug resistance testing based on the pol region before treatment initiation and therefore, these generated sequences could additionally be employed for surveillance of HIV subtypes and CRFs and detection of emerging new CRFs with possible implications for HIV prevention, diagnosis, and management of HIV-positive individuals [7].
Seven different subtyping tools were employed for the initial subtyping of the pol region. The principal role of five of the aforementioned tools is the determination of the HIV-1 subtype; however, Stanford and Geno2pheno are used foremost for genotypic drug resistance interpretation, yet do allow for additional information of HIV-1 subtype. The two were included to compare the results with other subtyping tools and indeed Geno2pheno did not report any evidence of recombination, probably because this tool was not designed to do so in the first place. The Stanford algorithm showed recombination in half of the URFs subsequently characterized by NFLG sequencing and was thus more in line with other more complex subtyping algorithms. Because all the subtyping tools available seem to have some limitations, a combination of tools should preferably be used when looking for recombination.
Unfortunately, we were not able to obtain all four amplicons for seven of the samples selected for NFLG sequencing using the protocol by Gall et al. (2012) [23]. This protocol faces some problems, especially when older samples with a lower viral load were processed. For example, Sample 19 was detected in a transmission pair with CRF90_BF and Sample 22 was observed basal to the CRF02_AG cluster on the phylogenetic tree of the pol region and therefore, both were indicated for NFLG sequencing. Sadly, we could not obtain all the PCR products and further analyze these two samples of interest. In addition, read coverage of the fourth amplicon was too low for two samples and we were therefore not able to obtain the complete sequences.
The highest HIV-1 genetic diversity was observed in the Democratic Republic of Congo, where the HIV-1 epidemic was shown to originate [47,48]. Due to error prone reverse transcriptase enzyme and high replication rates, HIV has the ability to generate genetically diverse viral variants within a single individual that can differ by up to 10% [49,50,51]. Variation within a subtype and between subtypes is even greater and can be as high as 30% and 42% at the amino acid level, respectively [6,51]. That is why the HIV pandemic is characterized by the remarkable global genetic diversity, reflected in the abundance of HIV subtypes and recombinant forms circulating in the world. Apart from HIV subtypes, several CRFs have obtained epidemic proportions and have become drivers of the HIV epidemic in various regions, such as CRF01_AE in Southeast Asia and CRF02_AG in West Africa. Furthermore, CRFs and URFs are continuously increasing in prevalence [7]. However, what specific characteristics of CRFs contribute to their dispersal is still to be elucidated. Propagation of a newly generated URF into a CRF with epidemic significance is probably a result of a mosaic of factors. Genetic characteristics of the virus itself and immunological traits of the population at risk of infection both contribute to transmissibility and disease progression [7,52,53]. These biological properties coupled with socioeconomic factors, such as transportation links, migration, population growth, access to treatment, and to some extent, stochastic effect, all add to the effective spread of a particular viral variant [7,48,54,55].
Since HIV diversity plays a central role in the HIV pandemic, today it is imperative that global molecular epidemiology surveillance is continued and improved using rigorous sampling strategies in different parts of the world. With the implementation of next-generation sequencing in smaller-scale laboratories, this technique allows identification of new CRFs, even in countries where HIV diversity is substantially limited, like Slovenia. For this reason, Rodgers et al. (2017) propose making changes to HIV-1 nomenclature and suggest that it should only include CRFs of epidemiological significance [47]. Indeed, as mentioned above, it is difficult to predict which of the CRFs identified will reach epidemic potential, thus their surveillance is of utmost importance. For example, CRFs such as the previously mentioned CRF19_cpx and CRF14_BG were associated with rapid disease progression to AIDS and death, suggesting high pathogenicity and virulence of these strains [41,56]. In addition, HIV variants can differ in their rate of transmission and can influence the sensitivity of diagnostic tests and viral load measurements and the effectiveness of ART and vaccines, highlighting the importance of continuous surveillance of circulating recombinant forms [6].

5. Conclusions

Six unique recombinant forms were identified among individuals diagnosed with HIV-1 between 2000 and 2016 in Slovenia, and 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 limited, emphasizing the importance of the surveillance of HIV-1 recombinant forms.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/12/1/63/s1, Figure S1: SimPlot bootscan analysis of the Slovenian near full-length genome sequences determined to be non-recombinant, Figure S2: Phylogenetic analyses of the segments of the alignment of the near full-length genome sequences determined to be recombinant.

Author Contributions

Conceptualization, M.M.L. and M.P.; methodology, M.M.L.; formal analysis, M.M.L. and T.M.Z.; investigation, M.M.L. and J.M.; resources, M.P.; data curation, M.M.L. and T.M.Z.; writing—original draft preparation, M.M.L.; writing—review and editing, T.M.Z. and M.P.; visualization, M.M.L.; supervision, M.P.; project administration, M.M.L.; funding acquisition, M.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank Ana B. Abecasis for her kind guidance on the conceptualization of the research.

Conflicts of Interest

The authors declare no potential conflict of interest. M.M.L. received speaker’s honoraria and travel support from Gilead Sciences.

References

  1. UNAIDS. Global HIV & AIDS Statistics—2019 Fact Sheet. Available online: https://www.unaids.org/en/resources/fact-sheet (accessed on 16 October 2019).
  2. Ghosn, J.; Taiwo, B.; Seedat, S.; Autran, B.; Katlama, C. HIV. Lancet 2018, 392, 685–697. [Google Scholar] [CrossRef]
  3. HIV Circulating Recombinant Forms (CRFs). Available online: https://www.hiv.lanl.gov/content/sequence/HIV/CRFs/CRFs.html (accessed on 18 October 2019).
  4. Robertson, D.L.; Anderson, J.P.; Bradac, J.A.; Carr, J.K.; Foley, B.; Funkhouser, R.K.; Gao, F.; Hahn, B.H.; Kalish, M.L.; Kuiken, C.; et al. HIV-1 nomenclature proposal. Science 2000, 288, 55–56. [Google Scholar] [CrossRef] [PubMed]
  5. Yamaguchi, J.; McArthur, C.; Vallari, A.; Sthreshley, L.; Cloherty, G.A.; Berg, M.G.; Rodgers, M.A. Complete genome sequence of CG-0018a-01 establishes HIV-1 subtype L. J. Acquir. Immune Defic. Syndr. 2019. [Google Scholar] [CrossRef] [PubMed]
  6. Hemelaar, J. The origin and diversity of the HIV-1 pandemic. Trends Mol. Med. 2012, 18, 182–192. [Google Scholar] [CrossRef] [PubMed]
  7. Hemelaar, J.; Elangovan, R.; Yun, J.; Dickson-Tetteh, L.; Fleminger, I.; Kirtley, S.; Williams, B.; Gouws-Williams, E.; Ghys, P.D.; WHO–UNAIDS Network for HIV Isolation Characterisation. Global and regional molecular epidemiology of HIV-1, 1990–2015: A systematic review, global survey, and trend analysis. Lancet Infect. Dis. 2019, 19, 143–155. [Google Scholar] [CrossRef]
  8. Babič, D.Z.; Zelnikar, M.; Seme, K.; Vandamme, A.M.; Snoeck, J.; Tomažič, J.; Vidmar, L.; Karner, P.; Poljak, M. Prevalence of antiretroviral drug resistance mutations and HIV-1 non-B subtypes in newly diagnosed drug-naïve patients in Slovenia, 2000–2004. Virus Res. 2006, 118, 156–163. [Google Scholar] [CrossRef]
  9. Lunar, M.M.; Židovec Lepej, S.; Abecasis, A.B.; Tomažič, J.; Vidmar, L.; Karner, P.; Vovko, T.D.; Pečavar, B.; Maver, P.J.; Seme, K.; et al. Short communication: Prevalence of HIV type 1 transmitted drug resistance in Slovenia: 2005–2010. AIDS Res. Hum. Retrovir. 2013, 29, 343–349. [Google Scholar] [CrossRef] [Green Version]
  10. Lunar, M.M.; Židovec Lepej, S.; Tomažič, J.; Vovko, T.D.; Pečavar, B.; Turel, G.; Maver, M.; Poljak, M. HIV-1 transmitted drug resistance in Slovenia and its impact on predicted treatment effectiveness: 2011–2016 update. PLoS ONE 2018, 13, e0196670. [Google Scholar] [CrossRef] [Green Version]
  11. de Oliveira, T.; Deforche, K.; Cassol, S.; Salminen, M.; Paraskevis, D.; Seebregts, C.; Snoeck, J.; van Rensburg, E.J.; Wensing, A.M.; van de Vijver, D.A.; et al. An automated genotyping system for analysis of HIV-1 and other microbial sequences. Bioinformatics 2005, 21, 3797–3800. [Google Scholar] [CrossRef] [Green Version]
  12. Pineda-Peña, A.C.; Faria, N.R.; Imbrechts, S.; Libin, P.; Abecasis, A.B.; Deforche, K.; Gómez-López, A.; Camacho, R.J.; de Oliveira, T.; Vandamme, A.M. Automated subtyping of HIV-1 genetic sequences for clinical and surveillance purposes: Performance evaluation of the new REGA version 3 and seven other tools. Infect. Genet. Evol. 2013, 19, 337–348. [Google Scholar] [CrossRef] [Green Version]
  13. Struck, D.; Lawyer, G.; Ternes, A.M.; Schmit, J.C.; Bercoff, D.P. COMET: Adaptive context-based modeling for ultrafast HIV-1 subtype identification. Nucleic Acids Res. 2014, 42, e144. [Google Scholar] [CrossRef] [PubMed]
  14. Schultz, A.K.; Zhang, M.; Bulla, I.; Leitner, T.; Korber, B.; Morgenstern, B.; Stanke, M. jpHMM: Improving the reliability of recombination prediction in HIV-1. Nucleic Acids Res. 2009, 37, W647–W651. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Kosakovsky Pond, S.L.; Posada, D.; Stawiski, E.; Chappey, C.; Poon, A.F.; Hughes, G.; Fearnhill, E.; Gravenor, M.B.; Leigh Brown, A.J.; Frost, S.D. An evolutionary model-based algorithm for accurate phylogenetic breakpoint mapping and subtype prediction in HIV-1. PLoS Comput. Biol. 2009, 5, e1000581. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Beerenwinkel, N.; Däumer, M.; Oette, M.; Korn, K.; Hoffmann, D.; Kaiser, R.; Lengauer, T.; Selbig, J.; Walter, H. Geno2pheno: Estimating phenotypic drug resistance from HIV-1 genotypes. Nucleic Acids Res. 2003, 31, 3850–3855. [Google Scholar] [CrossRef] [Green Version]
  17. Rhee, S.Y.; Gonzales, M.J.; Kantor, R.; Betts, B.J.; Ravela, J.; Shafer, R.W. Human immunodeficiency virus reverse transcriptase and protease sequence database. Nucleic Acids Res. 2003, 31, 298–303. [Google Scholar] [CrossRef] [Green Version]
  18. HIV Sequence Database. Available online: https://www.hiv.lanl.gov/content/sequence/HIV/mainpage.html (accessed on 23 February 2018).
  19. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  20. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  21. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef] [Green Version]
  22. FigTree. Available online: http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 11 December 2018).
  23. Gall, A.; Ferns, B.; Morris, C.; Watson, S.; Cotton, M.; Robinson, M.; Berry, N.; Pillay, D.; Kellam, P. Universal amplification, next-generation sequencing, and assembly of HIV-1 genomes. J. Clin. Microbiol. 2012, 50, 3838–3844. [Google Scholar] [CrossRef] [Green Version]
  24. Wymant, C.; Blanquart, F.; Golubchik, T.; Gall, A.; Bakker, M.; Bezemer, D.; Croucher, N.J.; Hall, M.; Hillebregt, M.; Ong, S.H.; et al. Easy and accurate reconstruction of whole HIV genomes from short-read sequence data with shiver. Virus Evol. 2018, 4, vey007. [Google Scholar] [CrossRef] [Green Version]
  25. Hunt, M.; Gall, A.; Ong, S.H.; Brener, J.; Ferns, B.; Goulder, P.; Nastouli, E.; Keane, J.A.; Kellam, P.; Otto, T.D. IVA: Accurate de novo assembly of RNA virus genomes. Bioinformatics 2015, 31, 2374–2376. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  27. Fastaq. Available online: https://github.com/sanger-pathogens/Fastaq (accessed on 8 May 2019).
  28. SMALT. Available online: https://www.sanger.ac.uk/science/tools/smalt-0 (accessed on 8 May 2019).
  29. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; 1000 Genome Project Data Processing Subgroup. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Lole, K.S.; Bollinger, R.C.; Paranjape, R.S.; Gadkari, D.; Kulkarni, S.S.; Novak, N.G.; Ingersoll, R.; Sheppard, H.W.; Ray, S.C. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J. Virol. 1999, 73, 152–160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Recombinant HIV-1 Drawing Tool. Available online: https://www.hiv.lanl.gov/content/sequence/DRAW_CRF/recom_mapper.html (accessed on 3 September 2019).
  32. Lai, A.; Simonetti, F.R.; Brindicci, G.; Bergna, A.; Di Giambenedetto, S.; Sterrantino, G.; Mussini, C.; Menzo, S.; Bagnarelli, P.; Zazzi, M.; et al. Local epidemics gone viral: Evolution and diffusion of the Italian HIV-1 recombinant form CRF60_BC. Front. Microbiol. 2019, 10, 769. [Google Scholar] [CrossRef] [PubMed]
  33. Lebedev, A.; Lebedeva, N.; Moskaleychik, F.; Pronin, A.; Kazennova, E.; Bobkova, M. Human immunodeficiency virus-1 diversity in the Moscow region, Russia: Phylodynamics of the most common subtypes. Front. Microbiol. 2019, 10, 320. [Google Scholar] [CrossRef]
  34. Díez-Fuertes, F.; Cabello, M.; Thomson, M.M. Bayesian phylogeographic analyses clarify the origin of the HIV-1 subtype A variant circulating in former Soviet Union’s countries. Infect. Genet. Evol. 2015, 33, 197–205. [Google Scholar] [CrossRef]
  35. Machado, L.Y.; Blanco, M.; López, L.S.; Díaz, H.M.; Dubed, M.; Valdés, N.; Noa, E.; Martínez, L.; Pérez, M.T.; Romay, D.M.; et al. National survey of pre-treatment HIV drug resistance in Cuban patients. PLoS ONE 2019, 14, e0221879. [Google Scholar] [CrossRef]
  36. Galindo, J.A.P.; Torres-Puente, M.; Gimeno, C.; Ortega, E.; Navarro, D.; Galindo, M.J.; Navarro, L.; Navarro, V.; Juan, A.; Belda, J.; et al. Expansion of the CRF19_cpx variant in Spain. J. Clin. Virol. 2015, 69, 146–149. [Google Scholar] [CrossRef]
  37. González-Domenech, C.M.; Viciana, I.; Delaye, L.; Mayorga, M.L.; Palacios, R.; de la Torre, J.; Jarilla, F.; Castaño, M.; Del Arco, A.; Clavijo, E.; et al. Emergence as an outbreak of the HIV-1 CRF19_cpx variant in treatment-naïve patients in southern Spain. PLoS ONE 2018, 13, e0190544. [Google Scholar] [CrossRef] [Green Version]
  38. Pérez-Parra, S.; Álvarez, M.; Fernandez-Caballero, J.A.; Pérez, A.B.; Santos, J.; Bisbal, O.; Aguilera, A.; Rivero, M.; García-Fraile, L.; García, F. Continued propagation of the CRF19_cpx variant among HIV-positive MSM patients in Spain. J. Antimicrob. Chemother. 2018, 73, 1031–1038. [Google Scholar] [CrossRef] [PubMed]
  39. Fernández-García, A.; Delgado, E.; Cuevas, M.T.; Vega, Y.; Montero, V.; Sánchez, M.; Carrera, C.; López-Álvarez, M.J.; Miralles, C.; Pérez-Castro, S.; et al. Identification of an HIV-1 BG intersubtype recombinant form (CRF73_BG), partially related to CRF14_BG, which is circulating in Portugal and Spain. PLoS ONE 2016, 11, e0148549. [Google Scholar] [CrossRef] [PubMed]
  40. Beamud, B.; Bracho, M.A.; González-Candelas, F. Characterization of new recombinant forms of HIV-1 from the Comunitat Valenciana (Spain) by phylogenetic incongruence. Front. Microbiol. 2019, 10, 1006. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Kouri, V.; Khouri, R.; Alemán, Y.; Abrahantes, Y.; Vercauteren, J.; Pineda-Peña, A.C.; Theys, K.; Megens, S.; Moutschen, M.; Pfeifer, N.; et al. CRF19_cpx is an evolutionary fit HIV-1 variant strongly associated with rapid progression to AIDS in Cuba. EBioMedicine 2015, 2, 244–254. [Google Scholar] [CrossRef] [Green Version]
  42. EACS Guidelines Version 8.1, October 2016. Available online: https://www.eacsociety.org/files/guidelines_8.1-english.pdf (accessed on 5 November 2019).
  43. Delatorre, E.; Bello, G. Time-scale of minor HIV-1 complex circulating recombinant forms from Central and West Africa. BMC Evol. Biol. 2016, 16, 249. [Google Scholar] [CrossRef] [Green Version]
  44. Leoz, M.; Feyertag, F.; Charpentier, C.; Delaugerre, C.; Wirden, M.; Lemee, V.; Plantier, J.C. Characterization of CRF56_cpx, a new circulating B/CRF02/G recombinant form identified in MSM in France. AIDS 2013, 27, 2309–2312. [Google Scholar] [CrossRef]
  45. Tariq, U.; Parveen, A.; Akhtar, F.; Mahmood, F.; Ali, S.; Abidi, S.H. Emergence of circulating recombinant form 56_cpx in Pakistan. AIDS Res. Hum. Retrovir. 2018, 34, 1002–1004. [Google Scholar] [CrossRef]
  46. Tongo, M.; Dorfman, J.R.; Abrahams, M.R.; Mpoudi-Ngole, E.; Burgers, W.A.; Martin, D.P. Near full-length HIV type 1M genomic sequences from Cameroon: Evidence of early diverging under-sampled lineages in the country. Evol. Med. Public Health 2015, 2015, 254–265. [Google Scholar] [CrossRef] [Green Version]
  47. Rodgers, M.A.; Wilkinson, E.; Vallari, A.; McArthur, C.; Sthreshley, L.; Brennan, C.A.; Cloherty, G.; de Oliveira, T. Sensitive next-generation sequencing method reveals deep genetic diversity of HIV-1 in the Democratic Republic of the Congo. J. Virol. 2017, 91, e01841-16. [Google Scholar] [CrossRef] [Green Version]
  48. Faria, N.R.; Rambaut, A.; Suchard, M.A.; Baele, G.; Bedford, T.; Ward, M.J.; Tatem, A.J.; Sousa, J.D.; Arinaminpathy, N.; Pépin, J.; et al. HIV epidemiology. The early spread and epidemic ignition of HIV-1 in human populations. Science 2014, 346, 56–61. [Google Scholar] [CrossRef] [Green Version]
  49. Roberts, J.D.; Bebenek, K.; Kunkel, T.A. The accuracy of reverse transcriptase from HIV-1. Science 1988, 242, 1171–1173. [Google Scholar] [CrossRef] [PubMed]
  50. Ho, D.D.; Neumann, A.U.; Perelson, A.S.; Chen, W.; Leonard, J.M.; Markowitz, M. Rapid turnover of plasma virions and CD4 lymphocytes in HIV-1 infection. Nature 1995, 373, 123–126. [Google Scholar] [CrossRef] [PubMed]
  51. Korber, B.; Gaschen, B.; Yusim, K.; Thakallapally, R.; Kesmir, C.; Detours, V. Evolutionary and immunological implications of contemporary HIV-1 variation. Br. Med. Bull. 2001, 58, 19–42. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Kiwanuka, N.; Laeyendecker, O.; Quinn, T.C.; Wawer, M.J.; Shepherd, J.; Robb, M.; Kigozi, G.; Kagaayi, J.; Serwadda, D.; Makumbi, F.E.; et al. HIV-1 subtypes and differences in heterosexual HIV transmission among HIV-discordant couples in Rakai, Uganda. AIDS 2009, 23, 2479–2484. [Google Scholar] [CrossRef] [PubMed]
  53. Venner, C.M.; Nankya, I.; Kyeyune, F.; Demers, K.; Kwok, C.; Chen, P.L.; Rwambuya, S.; Munjoma, M.; Chipato, T.; Byamugisha, J.; et al. Infecting HIV-1 subtype predicts disease progression in women of Sub-Saharan Africa. EBioMedicine 2016, 13, 305–314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Tebit, D.M.; Arts, E.J. Tracking a century of global expansion and evolution of HIV to drive understanding and to combat disease. Lancet Infect. Dis. 2011, 11, 45–56. [Google Scholar] [CrossRef]
  55. Oberle, C.S.; Joos, B.; Rusert, P.; Campbell, N.K.; Beauparlant, D.; Kuster, H.; Weber, J.; Schenkel, C.D.; Scherrer, A.U.; Magnus, C.; et al. Tracing HIV-1 transmission: Envelope traits of HIV-1 transmitter and recipient pairs. Retrovirology 2016, 13, 62. [Google Scholar] [CrossRef] [Green Version]
  56. Bártolo, I.; Camacho, R.; Barroso, H.; Bezerra, V.; Taveira, N. Rapid clinical progression to AIDS and death in a persistently seronegative HIV-1 infected heterosexual young man. AIDS 2009, 23, 2359–2362. [Google Scholar] [CrossRef]
Viruses 12 00063 g001 550
Figure 1. Maximum likelihood phylogenetic tree of the pol region of the divergent Slovenian sequences. Sample numbers selected for further near full-length genome characterization are depicted at the tips of the branches, which are colored red if the near full-length genome sequence was obtained, yellow if the sequencing was not successful, and blue if it was not initially selected for sequencing. Subtype reference sequences of subtypes A–D, F–H, J, K, CRF01_AE, and CRF02_AG are depicted in green, other CRF references are shown in purple, and other control sequences are in black.
Figure 1. Maximum likelihood phylogenetic tree of the pol region of the divergent Slovenian sequences. Sample numbers selected for further near full-length genome characterization are depicted at the tips of the branches, which are colored red if the near full-length genome sequence was obtained, yellow if the sequencing was not successful, and blue if it was not initially selected for sequencing. Subtype reference sequences of subtypes A–D, F–H, J, K, CRF01_AE, and CRF02_AG are depicted in green, other CRF references are shown in purple, and other control sequences are in black.
Viruses 12 00063 g001
Viruses 12 00063 g002 550
Figure 2. SimPlot bootscan analysis of the Slovenian near full-length HIV-1 genome sequences characterized as recombinants, where the x-axis represents the nucleotide position in the alignment.
Figure 2. SimPlot bootscan analysis of the Slovenian near full-length HIV-1 genome sequences characterized as recombinants, where the x-axis represents the nucleotide position in the alignment.
Viruses 12 00063 g002
Viruses 12 00063 g003 550
Figure 3. Final unique recombinant forms visualized with the Recombinant HIV-1 Drawing Tool, where recombination breakpoints are depicted according to nucleotide position in the reference HXB2.
Figure 3. Final unique recombinant forms visualized with the Recombinant HIV-1 Drawing Tool, where recombination breakpoints are depicted according to nucleotide position in the reference HXB2.
Viruses 12 00063 g003
Table
Table 1. Subtyping results of sequences that had at least one divergent result with the seven subtyping tools. Samples subsequently selected for further near full-length genome sequencing according to phylogenetic analysis are shown in bold.
Table 1. Subtyping results of sequences that had at least one divergent result with the seven subtyping tools. Samples subsequently selected for further near full-length genome sequencing according to phylogenetic analysis are shown in bold.
No.Accession NumberRega 2.0Rega 3.0Comet 1.0jpHMMSCUEALStanfordGeno2pheno 3.3
1AJ971094NA02_AG02_AGA1, G02_AG-like02_AG/02_AG02_AG
2AJ971095NA02_AG02_AGA1, G02_AG-like/complex02_AG/02_AG02_AG
3AJ971100F (F1)F (F1)F1F1F1D/FF1
4AJ971096NA02_AG02_AGA1, G02_AG-like02_AG/02_AG02_AG
5AJ971133NAB-likeBBB, D recombinant/BB/BB
6GQ399318BRecombinant of B, DBBBB/BB
7AM113750NA02_AG02_AGA1, Gcomplex02_AG/02_AG02_AG
8JX046417A (A1)A (A1)A1A1A-ancestral, A1 rec/A, A1 recombinant/A1A/CRF01_AE01_AE
9JX046416A (A1)A (01_AE)01_AE01_AEAE01_AE/01_AE01_AE
10JX046415NA02_AG02_AGA1, Gcomplex02_AG/02_AG02_AG
11JX046413A (A1)A (A1)A1A1complex/A, A-ancestral recombinantA/01_AE01_AE
12JX028352BNA (recombinant of B, D)BBBB/BB
13JX04641419_cpx19_cpxunassigned; D, 11_cpx, G, 20_BGDcomplexD/DD
14JX028338NABBBBB/BB
15JX046406DDDB, DDD/DD
16JX046405DDDB, DDD/DD
17JX028323NABBBBB/BB
18JX046407NAB, potential recombinantBBBB/BB
19JX046408NARecombinant of F1, 05_DFunassigned; B, F1B, F1B, F1 recombinantB/FF1
20JX046409A (A1)A (A1)A1A1A1A/01_AE01_AE
21JX046410A (A1)A (A1)A1A1A1A/01_AE01_AE
22KP013668GGunassigned; 18_cpx, 02_AG, GGGG/G02_AG
23JX046412A (A1)A (A1)A1A1A1A/01_AE01_AE
24JX046404A (A1)A (A1)A1A1A1A/01_AE01_AE
25JX046411A (A1)A (A1)A1A1A101_AE/01_AE01_AE
26KP013669A (A1)A (A1)A1A1A1A/01_AE01_AE
27JX028308BNA, Recombinant of B, DBBBB/BB
28JX046403A (A1)A (A1)A1A1A1A/01_AE01_AE
29KF753740NARecombinant of B, DBBBB/BB
30KF753737NARecombinant of B, DBBBB/BB
31KF753747NAC, potential recombinantunassigned; F1, A1, 01_AE, B, C, 60_BCA1, B, CcomplexC/DC
32KF753751A (A1)A (A1)A1A1A101_AE/01_AE01_AE
33KF753741A (A1)A (A1)A1A1A1A/01_AE01_AE
34KF753731GGunassigned; D, J, GG, JAE, G recombinantA/GG
35KF753711BBunassigned; B-DBBB/BB
36KP013665BBBBB, D recombinant/BB/BB
37KF753704NANA (G, B)G (check for 02_AG)A1, B, Gcomplex02_AG/B02_AG
38KF753709NABBBBB/BB
39KF753699A (A1)A (01_AE)01_AE01_AE15_01B/A1, AE recombinant/AE, F-ancestral recombinant01_AE/01_AE01_AE
40KP013656BBB (check for 29_BF)BBB/BB
41KP013648A (A1)A (A1)A1A1A101_AE/AA1
42KP013643NAB-likeBBBB/BB
43KP013657NARecombinant of B, DBBBB/BB
44KP013641A (A1)A (A1)A1A1A1A/01_AE01_AE
45KP013660A (A1)A (A1)unassigned; D, 10_CD, 01_AE, B01_AEA1, AE recombinant01_AE/01_AE01_AE
46KP013644NABBBBB/BB
47KY656620BB-likeBBBBB
48KY656621NAB-likeBBBBB
49MN736711BBBBB, C recombinantBB
50KY656626A (A1)A (A1)A1A1A1A (A_FSU)01_AE
51KY656630NA02_AG02_AGA1, GA4, G recombinant02_AG02_AG
52KY656636A (A1)A (A1)A1A1A1A (A_FSU)01_AE
53KY656616NAB-likeBBBBB
54KY656639NAG (02_AG)02_AGF1, GG02_AG02_AG
55KY656643A (A1)A (A1)A1A1A1A01_AE
56KY656648NARecombinant of B, DBBBBB
57KY656649BBBBB, F1 recombinantBB
58KY656650BRecombinant of B, DBBBBB
59KY656640A (01_AE)A (01_AE)01_AE (check for 15_01B)01_AEAE01_AE01_AE
60KY656656NAB-likeBBBBB
61KY656657A (A1)A (A1)A1A1A1A, 01_AE01_AE
62KY656658A (A1)A (A1)A1A1A1A, 01_AE01_AE
63KY656659A (A1)A (A1)A1A1A1A, 01_AE01_AE
64MF987697NARecombinant of B, DBBBBB
65KY656670NARecombinant of B, DBBBBB
66KY656671A (A1)A (A1)A1A1A1A/01_AE01_AE
67MF987719NANAunassigned, B–GDcomplexDD
68MF987720NANAunassigned, B–GDcomplexDD
NA = not assigned; FSU = former Soviet Union.
Table
Table 2. Subtyping results of the near full-length genome sequences.
Table 2. Subtyping results of the near full-length genome sequences.
No.Accession NumberGenome Sequence Range 1Rega 3.0CometjpHMMSimPlotFinal Result
11MN736709497–7842A (A1)A1A1, K, A1A1/K or 45_cpxA|A/U|A|K|A
12MN736698497–9496BBBBB
25MN736699497–9487A (A1)A1A1A6A6
29MN736700497–9496BBBBB
31MN736710497–7843Cunassigned; 60_BC, B, 82_cpx, 01_AE, 58_01B, 01_AE, 58_01B, 01_AE, 58_01B, 59_01B, 58_01B, 59_01B, 58_01B, 59_01B, 01_AE, 59_01B, 58_01B, 02_AG, 58_01B, 02_AG, A1, 69_01B, 59_01B, 69_01B, 53_01B, 69_01B, 53_01B, 69_01B, 55_01[..]C, B, C, 01_AE, C60_BC60_BC|01_AE|60_BC
34MN736701497–9496recombinant of 13_cpx, G, A1, J, Hunassigned; G, A1, 01_AE, D, J, C, A2, C, B, D, B, D, B, D, J, C, A1, 01_AEA1, G, D, J, B, J, G, A1, 01_AEcomplexU|A|G|U|G|J|
G|U|J|U|G|
A|01_AE
37MN736702497–9433recombinant of 02_AG, B, G, A1unassigned; 02_AG, 20_BG, 56_cpx, 20_BG, 56_cpx, 90_BF1, B, 56_cpx, B, 56_cpx, B, 71_BF1, 51_01B, 71_BF1, 39_BF, B, 69_01B, B, 69_01B, B, 69_01B, B, 56_cpx, 19_cpx, 24_BG, 19_cpx, 56_cpx, 24_BG, 56_cpx, G, 14_BG, 43_02G, G, 43_02G[..]A1, B, A1, G, B, A1, G, B, A1, B, G56_cpx, B56_cpx|B|56_cpx
39MN736703497–949601_AEunassigned; 01_AE, B, D01_AE01_AE01_AE
43MN736704497–9412Bunassigned; B, C, GBB, B/D, B, B/DB
54MN736705497–944502_AG02_AGA1, G, A1, G, A1, G, A1, G, A102_AG02_AG|A|02_AG
56MN736706497–9496BBBBB
65MN736707497–9496BBBBB
68MN736708497–9496recombinant of G, A1, D, B, F1unassigned; G, D, B, D, B, D, A1A1, G, D, B, A1, G, A2, G, A1, G, A1, GcomplexU|19_cpx|G|19_cpx|
B|19_cpx|A|19_cpx|U|19_cpx|G|19_cpx|G|19_cpx|G
1 nucleotide position according to HXB2, U = unknown.

Share and Cite

MDPI and ACS Style

Lunar, M.M.; Mlakar, J.; Zorec, T.M.; Poljak, M. HIV-1 Unique Recombinant Forms Identified in Slovenia and Their Characterization by Near Full-Length Genome Sequencing. Viruses 2020, 12, 63. https://0-doi-org.brum.beds.ac.uk/10.3390/v12010063

AMA Style

Lunar MM, Mlakar J, Zorec TM, Poljak M. HIV-1 Unique Recombinant Forms Identified in Slovenia and Their Characterization by Near Full-Length Genome Sequencing. Viruses. 2020; 12(1):63. https://0-doi-org.brum.beds.ac.uk/10.3390/v12010063

Chicago/Turabian Style

Lunar, Maja M., Jana Mlakar, Tomaž Mark Zorec, and Mario Poljak. 2020. "HIV-1 Unique Recombinant Forms Identified in Slovenia and Their Characterization by Near Full-Length Genome Sequencing" Viruses 12, no. 1: 63. https://0-doi-org.brum.beds.ac.uk/10.3390/v12010063

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop