Molecular characterization of HIV type 1 strains from newly diagnosed patients in Cyprus (2007-2009) recovers multiple clades including unique recombinant strains and lack of transmitted drug resistance

Molecular Epidemiology of HIV-1 Subtypes and Primary Antiretroviral Resistance Profiles in Northern Cyprus: First Data Series

The distribution of human immunodeficiency virus-1 (HIV-1) subtypes indicates difference from region to region and in risk groups acquiring the disease worldwide. Although subtype C is more in terms of total cases, subtype B is dominant in certain regions, especially in western and central Europe. Molecular epidemiological studies are essential for the control, effective treatment, and understanding in transmission routes of HIV-1 infection. This study aims to determine the molecular epidemiology and antiretroviral drug resistance profiles of HIV-1 in northern Cyprus. The study involved 71 naive HIV-positive patients diagnosed in northern Cyprus between 2016 and 2022. HIV-1 subtypes and circulating recombinant forms (CRFs) were identified by phylogenetic analysis (neighbor-joining method) of pol gene sequences. Drug resistance mutations were analyzed using the World Health Organization (WHO) lists of mutations for surveillance. The Stanford University HIVdb program was used to interpret drug resistance mutations. In our study, 40 of 71 samples were successfully sequenced. Subtype B of HIV-1 was dominant with a rate of 52.5%, followed by CRF02_AG (20%) and G (7.5%) subtypes. The rate of subtype B (71.4%) in northern Cyprus was significantly higher than in the other country of origin (p = .028). Antiretroviral drug resistance was found in 15% of the sequenced serum samples. Nucleoside/nucleotide reverse transcriptase inhibitor (NRTI), non-nucleoside nucleotide reverse transcriptase inhibitor (NNRTI), and protease inhibitor (PI) resistance rates were 10% (4/40), 7.5% (3/40), and 2.5% (1/40), respectively. According to the results, it is noteworthy that the dominant subtype circulating in northern Cyprus is the B subtype, and CRFs were detected at a higher rate than expected.

Comprehensive Genetic Characterization of Four Novel HIV-1 Circulating Recombinant Forms (CRF129_56G, CRF130_A1B, CRF131_A1B, and CRF138_cpx): Insights from Molecular Epidemiology in Cyprus

Molecular investigations of the HIV-1 pol region (2253-5250 in the HXB2 genome) were conducted on sequences obtained from 331 individuals infected with HIV-1 in Cyprus between 2017 and 2021. This study unveiled four distinct HIV-1 putative transmission clusters, encompassing 19 previously unidentified HIV-1 recombinants. These recombinants, each comprising eight, three, four, and four sequences, respectively, did not align with previously established Circulating Recombinant Forms (CRFs). To characterize these novel HIV-1 recombinants, near-full-length genome sequences were successfully obtained for 16 of the 19 recombinants (790-8795 in the HXB2 genome) using an in-house-developed RT-PCR assay. Phylogenetic analyses, employing MEGAX and Cluster-Picker, along with confirmatory neighbor-joining tree analyses of subregions, were conducted to identify distinct clusters and determine subtypes. The uniqueness of the HIV-1 recombinants was evident in their exclusive clustering within generated maximum likelihood trees. Recombination analyses highlighted the distinct chimeric nature of these recombinants, with consistent mosaic patterns observed across all sequences within each of the four putative transmission clusters. Conclusive genetic characterization identified four novel HIV-1 CRFs: CRF129_56G, CRF130_A1B, CRF131_A1B, and CRF138_cpx. CRF129_56G exhibited two recombination breakpoints and three fragments of subtypes CRF56_cpx and G. Both CRF130_A1B and CRF131_A1B featured seven recombination breakpoints and eight fragments of subtypes A1 and B. CRF138_cpx displayed five recombination breakpoints and six fragments of subtypes CRF22_01A1 and F2, along with an unclassified fragment. Additional BLAST analyses identified a Unique Recombinant Form (URF) of CRF138_cpx with three additional recombination sites, involving subtype F2, a fragment of unknown subtype origin, and CRF138_cpx. Post-identification, all putative transmission clusters remained active, with CRF130_A1B, CRF131_A1B, and CRF138_cpx clusters exhibiting further growth. Furthermore, international connections were identified through BLAST analyses, linking one sequence from the USA to the CRF130_A1B strain, and three sequences from Belgium and Cameroon to the CRF138_cpx strain. This study contributes valuable insights into the dynamic landscape of HIV-1 diversity and transmission patterns, emphasizing the need for ongoing molecular surveillance and global collaboration in tracking emerging viral variants.

Characterization of a novel HIV-1 circulating recombinant form, CRF91_cpx, comprising CRF02_AG, G, J, and U, mostly among men who have sex with men

Prospective molecular studies of HIV-1 pol region (2253-5250 in HXB2 genome) sequences from sequenced samples of 269 HIV-1-infected patients in Cyprus (2017-2021) revealed a transmission cluster of 14 unknown HIV-1 recombinants that were not classified as previously established CRFs. The earliest recombinant was collected in September 2017, and the transmission cluster continued to grow until November 2020. Near full-length HIV-1 genome sequences of the 11 of the 14 recombinants were successfully obtained (790-8795 in HXB2 genome) and aligned against a reference dataset of HIV-1 subtypes and CRFs. We employed MEGAX for maximum-likelihood tree construction (GTR model, 1000 bootstrap replicates), Cluster-Picker for phylogenetic clustering analysis (genetic distance ≤0.045, bootstrap support value ≥70%), and REGA-3.0 for subtype determination. Bootscan and similarity plot analyses (sliding window of 400 nucleotides overlapped by 40 nucleotides) were conducted using SimPlot-v3.5.1, and subregion confirmatory neighbour-joining tree analyses were conducted using MEGAX (Kimura two-parameter model, 1000 bootstrap replicates, ≥70% bootstrap-support value). Exclusive clustering of the HIV-1 recombinants revealed their uniqueness. The recombination analyses illustrated the same unique mosaic pattern with six putative intersubtype recombination breakpoints, seven fragments of subtypes CRF02_AG, G, J and an unclassified fragment. We conclusively characterized the mosaic structure of the novel HIV-1 CRF, named CRF91_cpx, by the Los Alamos HIV Sequence Database. Additionally, we identified a URF of CRF91_cpx with two additional recombination sites, generated by a recombination event between subtype B and CRF91_cpx. Since the identification of CRF91_cpx, two additional patient samples have been entered into the CRF91_cpx transmission cluster, demonstrating active growth.

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.

Molecular Epidemiology of the HIV-1 Subtype B Sub-Epidemic in Bulgaria

HIV-1 subtype B is the predominant strain in Bulgaria, yet little is known about the molecular epidemiology of these infections, including its origin and transmissibility. We used a phylodynamics approach by combining and analyzing 663 HIV-1 polymerase (pol) sequences collected from persons diagnosed with HIV/AIDS between 1988-2018 and associated epidemiologic data to better understand this sub-epidemic in Bulgaria. Using network analyses at a 1.5% genetic distance threshold (d) we found several large phylogenetic clusters composed mostly of men who have sex with men (MSM) and male heterosexuals (HET). However, at d = 0.5%, used to identify more recent transmission, the largest clusters dissociated to become smaller in size. The majority of female HET and persons with other transmission risks were singletons or pairs in the network. Phylogenetic analysis of the Bulgarian pol sequences with publicly available global sequences showed that subtype B was likely introduced into Bulgaria from multiple countries, including Israel and several European countries. Our findings indicate that subtype B was introduced into Bulgaria multiple times since 1988 and then infections rapidly spread among MSM and non-disclosed MSM. These high-risk behaviors continue to spread subtype B infection in Bulgaria as evidenced by the large clusters at d = 0.5%. Relatively low levels of antiretroviral drug resistance were observed in our study. Prevention strategies should continue to include increased testing and linkage to care and treatment, as well as expanded outreach to the MSM communities.

HIV-1 transmission networks across Cyprus (2010-2012)

A molecular epidemiology study of HIV-1 infection was conducted in one hundred diagnosed and untreated HIV-1-infected patients in Cyprus between 2010 and 2012, representing 65.4% of all the reported HIV-1 infections in Cyprus in this three-year period, using a previously defined enrolment strategy. Eighty-two patients were newly diagnosed (genotypic drug resistance testing within six months from diagnosis), and eighteen patients were HIV-1 diagnosed for a longer period or the diagnosis date was unknown. Phylogenetic trees of the pol sequences obtained in this study with reference sequences indicated that subtypes B and A1 were the most common subtypes present and accounted for 41.0 and 19.0% respectively, followed by subtype C (7.0%), F1 (8.0%), CRF02_AG (4.0%), A2 (2.0%), other circulating recombinant forms (CRFs) (7.0%) and unknown recombinant forms (URFs) (12%). Most of the newly-diagnosed study subjects were Cypriots (63%), males (78%) with median age 39 (Interquartile Range, IQR 33–48) reporting having sex with other men (MSM) (51%). A high rate of clustered transmission of subtype B drug-sensitive strains to reverse transcriptase and protease inhibitors was observed among MSM, twenty-eight out of forty-one MSM study subjects (68.0%) infected were implicated in five transmission clusters, two of which are sub-subtype A1 and three of which are subtype B strains. The two largest MSM subtype B clusters included nine and eight Cypriot men, respectively, living in all major cities in Cyprus. There were only three newly diagnosed patients with transmitted drug resistant HIV-1 strains, one study subject from the United Kingdom infected with subtype B strain and one from Romania with sub-subtype A2 strain, both with PI drug resistance mutation M46L and one from Greece with sub-subtype A1 with non-nucleoside reverse transcriptase inhibitors (NNRTI) drug resistance mutation K103N.

HIV-1 Infection in Cyprus, the Eastern Mediterranean European Frontier: A Densely Sampled Transmission Dynamics Analysis from 1986 to 2012

Since HIV-1 treatment is increasingly considered an effective preventionstrategy, it is important to study local HIV-1 epidemics to formulate tailored preventionpolicies. The prevalence of HIV-1 in Cyprus was historically low until 2005. To investigatethe shift in epidemiological trends, we studied the transmission dynamics of HIV-1 in Cyprususing a densely sampled Cypriot HIV-1 transmission cohort that included 85 percent ofHIV-1-infected individuals linked to clinical care between 1986 and 2012 based on detailedclinical, epidemiological, behavioral and HIV-1 genetic information. Subtyping andtransmission cluster reconstruction were performed using maximum likelihood and Bayesianmethods, and the transmission chain network was linked to the clinical, epidemiological andbehavioral data. The results reveal that for the main HIV-1 subtype A1 and B sub-epidemics,young and drug-naïve HIV-1-infected individuals in Cyprus are driving the dynamics of thelocal HIV-1 epidemic. The results of this study provide a better understanding of thedynamics of the HIV-1 infection in Cyprus, which may impact the development of preventionstrategies. Furthermore, this methodology for analyzing densely sampled transmissiondynamics is applicable to other geographic regions to implement effective HIV-1 preventionstrategies in local settings.

Origin and spread of HIV-1 in persons who inject drugs in Bulgaria

Increased HIV transmission in persons who inject drugs (PWIDs) has led to subepidemics and outbreaks in several countries in Europe, including Bulgaria. In this study in Bulgaria, we investigate the origin and spatiotemporal evolutionary history of HIV-1 infections in PWIDs and the distribution of antiretroviral resistance mutations and hepatitis co-infections in these populations. We analyzed HIV-1 polymerase sequences available from 117 of 359 PWIDs diagnosed with HIV/AIDS from 1999 to 2011. Of these, 50 (42.7%) were classified as CRF02_AG, 41 (35.0%) CRF01_AE, 12 (10.3%) URFs, ten (8.5%) subtype B, two (1.7%) subtype F1 and two (1.7%) CRF14_BG. Most recent common ancestor dating suggests that CRF01_AE was likely first introduced from Southeast Asia into persons reporting heterosexual infection in Bulgaria in 1992 and spread subsequently to PWIDs in the capital city of Sofia around 2003. Conversely, CRF02_AG in Bulgaria was likely first introduced into PWID from Germany in 2000 and later entered heterosexual populations around 2009. The overall prevalence of resistance mutations was 6.8% (8/117), of which 5.1% (5/117) was observed in patients on antiretroviral therapy and 1.7% (2/117) was from transmitted drug resistance mutations in drug-naïve individuals. 189/204 (92.6%) PWIDs were also co-infected with hepatitis C (HCV) and 31/183 (16.9%) were co-infected with hepatitis B (HBV). Our study provides valuable molecular epidemiological information on the introduction and distribution of the main HIV-1 subtypes, resistance mutations and hepatitis co-infections among PWIDs with HIV-1 in Bulgaria which can be used to target prevention efforts.

Phylogeographic analysis on the travel-related introduction of HIV-1 non-B subtypes to Northern Poland

Phylodynamic, sequence data based reconstructions for the surveillance of the geographic spatial spread are a powerful tool in molecular epidemiology. In this study region of origin for the set of 57 partial pol sequences derived from the patients the history of travel-related HIV transmission was analyzed using phylogeographic approach. Maximum likelihood trees based on the sets of country-annotated reference sequences were inferred for identified non-B variants. Region of sequence import was assigned using on the highest approximate likelihood ratios. Import of the A1 clades was traced to the Eastern Europe and associated with immigration from this region. Subtype C infections clustered most frequently with sequences of the South African origin while majority of subtype Ds were similar to the European clades. Subtype G sequences clustered with Portuguese lineage, CRF01_AE with Eastern or South-Eastern Asian. Eastern European, Middle African or Western African lineage was assigned for the CFR02_AG. Rare circulating recombinants originated either from Central Africa (CRF11_cpx - Democratic Republic of Congo, CRF13_cpx - Central African Republic, CRF37_cpx - Cameroon) or South America (CRF28_BF and CRF46_BF - Brazil). Import of the HIV-1 non-B variants, including recombinant forms previously rarely found in Poland and Europe is frequent among travelers. Observed founder events result in the heterosexually-driven introduction of the novel HIV-1 variants into the population.

Ten years survey of primary HIV-1 resistance in Serbia: the occurrence of multiclass resistance

In Serbia, the first cases of HIV infection were reported in 1985, whereas antiretroviral (ARV) therapy has been in use since 1987. With this study we aimed to assess the occurrence and pattern of HIV resistance mutations among newly diagnosed patients in the period 2002-2011. The study prospectively included 181 adult patients. Genotypic HIV-1 drug resistance testing was performed and drug resistance was scored according to the 2009 WHO list for surveillance of drug resistance mutations (SDRMs). A bioinformatic approach was used to estimate the duration of infection by calculating the percentage of ambiguous basecalls per sequence, with a cutoff of 0.47% as the delimiter for recent infection. The overall prevalence of transmitted drug resistance (TDR) found in the study was 8.8% (16/181, 95% CI=5.5-13.8). Thirty-one percent of resistant samples contained multiple SDRMs. In particular, 5/16 patients with resistance carried viral strains with SDRMs to multiple ARV classes, hence one-third of resistant strains were multiclass resistant, including non-B strains. A total of 51.9% of samples (94/181) were classified as recent infection, with a significant increase in the second part of the study period. However, the prevalence of TDR in recent infection was 6.4% (6/94, 95% CI=2.9-13.2), not statistically different from that found in nonrecent infection. We showed a changing pattern of TDR mutations over the study period, with a substantial occurrence of multiclass resistance, across different HIV subtypes. Our results highlight the need for continued surveillance of primary resistance.

Cellular HIV type 1 DNA levels are equivalent among drug-sensitive and drug-resistant strains in newly diagnosed and antiretroviral naive patients

The emergence of resistance against current antiretroviral drugs to human immunodeficiency virus type 1 (HIV-1) is an increasingly important concern to the continuous success of antiretroviral therapy to HIV-1-infected patients. In the past decade, a number of studies reported that the prevalence of transmitted drug resistance among newly diagnosed patients has reached an overall 9% prevalence worldwide. Also, a number of studies using longitudinal HIV-1 patient study cohorts demonstrated that the cellular HIV-1 DNA level in peripheral blood mononuclear cells (PBMCs) has a prognostic value for the progression of HIV-1 disease independently of plasma HIV-1 RNA load and CD4 count. Using a previously established molecular-beacon-based real-time PCR methodology, cellular HIV-1 DNA levels were quantified in newly diagnosed and antiretroviral-naive patients in Northern Greece recruited between 2009 and 2010 using a predefined enrolling strategy, in an effort to investigate whether there is any relationship between cellular HIV-1 DNA levels and HIV-1 transmitted drug resistance. As part of the same study, DNA sequences encoding the env (C2-C5 region of gp120) were also amplified from PBMC-extracted DNA in order to determine the genotypic coreceptor tropism and genetic subtype. Cellular HIV-1 DNA levels had a median of 3.309 log10 HIV-1 copies per 10(6) PBMCs and demonstrated no correlation between cellular HIV-1 DNA levels and HIV-1 transmitted drug resistance. An absence of association between cellular HIV-1 DNA levels with plasma viral HIV-1 RNA load and CD4 levels was also found reconfirming the previously published study. Genotypic analysis of coreceptor tropism indicated that 96% of samples, independently of the presence or not of genotypic drug resistance, were CCR5-tropic. Overall, the findings reconfirmed the previously proposed proposition that transmitted drug resistance does not have an impact on disease progression in HIV-1-infected individuals. Also, CCR5 coreceptor tropism dominance suggests that both drug-resistant and drug-sensitive strains behave similarly in early infection in newly diagnosed patients.

Short communication: molecular epidemiology of HIV type 1 infection in northern Greece (2009-2010): evidence of a transmission cluster of HIV type 1 subtype A1 drug-resistant strains among men who have sex with men

A prospective molecular epidemiology study of HIV-1 infection was conducted in newly diagnosed and antiretroviral-naive patients in Northern Greece between 2009 and 2010 using a predefined enrolling strategy. Phylogenetic trees of the pol sequences obtained in this study with reference sequences indicated that subtypes B and A1 were the most common subtypes present and accounted for 44.9% and 42.9%, respectively, followed by subtype C (3.1%), CRF02_AG (4.1%), CRF04_cpx (2.0%), and subtypes CRF01_01, F1, and G (1.0%). A high rate of clustered transmission of subtype A1-resistant strains to reverse transcriptase (RT) inhibitors was observed among men having sex with men. Indeed, 15 out of 17 study subjects (88.2%) infected with transmitted drug resistance (TDR) strains were implicated in transmission clusters, 10 of whom (66.7%) were men who have sex with men (MSM), and were also infected with subsubtype A1 strains. The main cluster within subtype A1 (I) included eight men reporting having sex with men from Thessaloniki infected with dual-class RT-resistant strains carrying both T215C and Y181C mutations.

Characterization of CRF56_cpx, a new circulating B/CRF02/G recombinant form identified in MSM in France

Several B/CRF02_AG Unique Recombinant Forms (URFs) have previously been identified in France. Here we show that one of them (URF5_B/02/G) is emerging in MSM, a high-risk population where HIV incidence and number of superinfections are increasing. We describe this new Circulating Recombinant Form, CRF56_cpx, estimate the time to its most recent common ancestor, investigate its origins and show that it probably shares common ancestors with strains from the East Mediterranean.

The emerging and re-emerging human immunodeficiency virus epidemics in Europe

We provide a summary of the current status of the human immunodeficiency virus (HIV) epidemic, including data on the transmission of drug-resistant virus in the European Region of the WHO. The review was conducted by searching the reports of the European Centre for Disease Prevention and Control and the United Nations General Assembly Special Session country reports to identify the number of HIV cases reported in 2002-2011, the number of HIV tests performed, and the results of the most recent HIV surveys in at-risk groups. In the West, a 5-year linear trend analysis suggests an increase in the number of reported HIV cases in Belgium, Greece, and Iceland, and a decline in The Netherlands, Ireland, Portugal, Switzerland, and the UK. In the Centre, the number of reported cases increased in Bulgaria and Turkey, and in the East, the number of reported cases increased in Armenia, Georgia, and Ukraine. Estonia and Latvia reported a significant downward trend. HIV transmission in men who have sex with men (MSM) and heterosexual transmission accounted for, respectively, 40% and 38% of newly reported HIV cases in the West in 2011, whereas the epidemic in the Centre is predominantly concentrated in MSM. Although injection drug use remains the key driving force of the epidemic in the East, there is increasing evidence of sexual transmission from the core groups of injection drug users, and increasing spread within MSM. The pattern of transmitted drug resistance (TDR) is heterogeneous accross Europe; a significant number of clusters of TDR have been recently observed in several European countries.

The analysis of near full-length genome sequences of HIV type 1 subtype A viruses from Russia supports the monophyly of major intrasubtype clusters

The HIV-1 epidemic in Russia has been insufficiently studied, with only 11 complete genome sequences from this country currently available, only three of which are of the locally predominant genetic form, the former Soviet Union (FSU) subtype A variant (A(FSU)). Here we analyze 10 newly derived A(FSU) near full-length genome sequences from Russia. Samples were selected based on phylogenetic clustering in protease-reverse transcriptase in two of the major A(FSU) clusters, V77I(PR) (n=6), widely circulating in Russia and other FSU countries, and A(SP1) (n=4), predominant in St. Petersburg. The phylogenetic analysis shows that the V77I(PR) genomes group in a monophyletic cluster together with 10 previously obtained A(FSU) genome sequences from Uzbekistan, Kazakhstan, Russia, and Cyprus, all bearing the V77I substitution in protease. Similarly, the four A(SP1) genomes group in a monophyletic cluster. These results therefore show that the monophyly of V77I(PR) and A(SP1) A(FSU) clusters is supported in near complete genomes.