Molecular epidemiology and HIV-1 variant evolution in Poland between 2015 and 2019

Identifying the unknown: Application of molecular epidemiology tools to identify clustering and HIV transmission routes in Poland

BACKGROUND

Understanding the dynamics of HIV-1 transmission is essential for developing effective screening and intervention strategies. Viral genetic sequences provide valuable information that can be used to infer the history and patterns of viral transmission.

PURPOSE

Our study explores the structure and dynamics of HIV transmissions in Poland from 1999 to 2022 to elucidate key patterns related with national epidemics.

METHODS

To understand the temporal dynamics of transmission routes we examined HIV pol sequence data from 5705 Polish PWH. The HIV-TRAnsmission Cluster Engine (HIV-TRACE) was utilized to identify potential links between different risk groups and putative links to individuals with unreported transmission risk.

RESULTS

Our analyses generated 503 clusters, containing 3942 individuals, and identified 13,917 putative links. Approximately 69.1 % of the sequences formed clusters. In the dataset 32.2 % of individuals were reported MSM transmission route, 7.9 % by heterosexual, and 5.6 % by PWID transmissions. The transmission route was unknown for 54.2 % of patients. Putative transmissions from MSM to all other groups revealed that 45.1 % of links lead to people with unregistered transmission mode. For heterosexual patients, 40.2 % of connections were directed to patients lacking information on infection routes and 30.5 % to MSM individuals. Our analysis unveiled that 45.1 % of cases with unreported transmission routes may be identified as MSM, while 3.5 % might be potential non-disclosed MSM.

CONCLUSIONS

Genetic linkages can provide valuable insights into the transmission dynamics among individuals, even in cases where transmission risk information is missing or unreported. The observed association between MSM and unreported cases highlights the potential of molecular epidemiology to complete missing patient data.

No Association Between HIV-1 Subtype and Primary Resistance Mutations with CD4 Reconstitution During Effective Antiretroviral Treatment: An Observational, Cohort Study

Some people with Human Immunodeficiency Virus (HIV) on effective antiretroviral therapy have persistent low lymphocyte CD4 counts and remain at an increased risk of Acquired Immunodeficiency Syndrome (AIDS). We investigated whether primary drug resistance mutations (DRMs) and HIV-1 subtype could be related to immunologic reconstitution in these people. In a multicenter, observational cohort study among treatment-naïve patients, we analyzed HIV-1 subtype, primary drug resistance mutations, CD4 counts, and CD4:CD8 ratios during effective antiretroviral therapy. We compared these variables between patients with different HIV subtypes and between those with or without drug-resistance mutations up to 48 weeks post-baseline. In 156 patients, CD4 count normalization (≥500 cells/µL) was observed in 39% of patients, while CD4:CD8 ratio ≥ 1 in 27% after treatment implementation. HIV-1 subtype B was present in 75% of the patients and subtype A in 22%. Primary resistance mutations were found in 57% of the individuals. The percentage of immunological nonrespondents did not differ significantly between those with different HIV subtypes or between those with or without primary resistance mutations (p > 0.05). In conclusion, there was no significant coincidence between the HIV subtype and primary drug resistance mutations with immunological reconstitution in patients receiving effective antiretroviral therapy.

Drug Resistance Mutations (DRMs) for Long-Acting Injectable Cabotegravir and Rilpivirine (CAB/RPV LAI) in the HIV-1 Subtype A6 Epidemic in Poland

HIV subtype A6 with the L74I polymorphism, which increases the risk of cabotegravir/rilpivirine treatment failure, causes more and more infections in Poland. In this multicenter, observational, cross-sectional study (2023-2024), we analyzed viral subtypes and drug-resistance mutations to drugs used for long-acting injectable antiretroviral treatment and pre-exposure prophylaxis. Among 357 people with HIV, 247 (69%) were Polish nationals, and 102 (29%) were from former Soviet Union countries. Of the 357 people included, 159 (45%) had subtype B, and 177 (50%) had subtype A6 infections, with 165 (87%) of the latter characterized by the L74I polymorphism. Subtype A6 was more frequent in women (66% vs. 46% in men, p < 0.05) and among people from former Soviet countries (77% vs. 39% in Polish nationals, p < 0.05). About 40% of people had either drug-resistance mutations for cabotegravir/rilpivirine or HIV A6 subtype with the L74I polymorphism; 4.5% had both of these conditions. Compared to subtype B infections, subtype A6 infections were characterized by more frequent major transmitted drug-resistance mutations for non-nucleoside reverse transcriptase inhibitors (8.5% vs. 1.9%, p = 0.007) and rilpivirine (5.1% vs. 0.6%, p = 0.016). Due to the frequent occurrence of the L74I polymorphism and drug-resistance mutations in HIV A6 subtype infection, about 40% of people with HIV in Poland may be at risk of long-acting injectable treatment failure.

Frequency of Major Transmitted Integrase Resistance in Poland Remains Low Despite Change in Subtype Variability

With the widespread use of integrase inhibitors and the expanding use of long-acting cabotegravir in both pre-exposure prophylaxis and antiretroviral treatment, molecular surveillance on the transmission of integrase resistance has regained clinical significance. This study aimed to determine the frequency of INSTI-transmitted drug resistance mutations (DRMs) among treatment-naïve individuals in Poland from 2016 to 2023. INSTI resistance was analyzed in 882 antiretroviral treatment-naïve individuals using Sanger sequencing. Integrase DRMs were defined based on the Stanford HIV drug resistance database scores. Phylogeny was used to investigate subtyping and clustering. For the analysis of time-trends, logistic regression was used. Major (E138K and R263K) integrase mutations were detected in 0.45% of cases with minor resistance observed in 14.85%, most commonly (13.95%) E157Q. Overall, no major clusters of transmitted drug resistance were identified, and the transmission of E157Q showed a decreasing trend (p < 0.001). While the frequency of sub-subtype A6 increased, it was predominantly found among migrants and associated with L74 mutations. The frequency of major integrase-transmitted DRMs remains low, despite the changes in subtype variability. Surveillance of changing HIV molecular variation patterns is vital from the perspective of the optimal use of integrase inhibitors, especially due to expanding long-acting cabotegravir implementation.

Investigating alignment-free machine learning methods for HIV-1 subtype classification

Motivation

Many viruses are organized into taxonomies of subtypes based on their genetic similarities. For human immunodeficiency virus 1 (HIV-1), subtype classification plays a crucial role in infection management. Sequence alignment-based methods for subtype classification are impractical for large datasets because they are costly and time-consuming. Alignment-free methods involve creating numerical representations for genetic sequences and applying statistical or machine learning methods. Despite their high overall accuracy, existing models perform poorly on less common subtypes. Furthermore, there is limited work investigating the impact of sequence vectorization methods, in particular natural language-inspired embedding methods, on HIV-1 subtype classification.

Results

We present a comprehensive analysis of sequence vectorization methods across machine learning methods. We report a k-mer-based XGBoost model with a balanced accuracy of 0.84, indicating that it has good overall performance for both common and uncommon HIV-1 subtypes. We also report a Word2Vec-based support vector machine that achieves promising results on precision and balanced accuracy. Our study sheds light on the effect of sequence vectorization methods on HIV-1 subtype classification and suggests that natural language-inspired encoding methods show promise. Our results could help to develop improved HIV-1 subtype classification methods, leading to improved individual patient outcomes, and the development of subtype-specific treatments.

Availability and implementation

Source code is available at https://www.github.com/kwade4/HIV_Subtypes.

Phylodynamic evolution of HIV-1 A6 sub-subtype epidemics in Poland

The human immunodeficiency virus type 1 (HIV-1) A6 sub-subtype is highly prevalent in Eastern Europe. Over the past decade, the dissemination of the A6 lineage has been expanding in Poland. The recent Russian invasion of Ukraine may further escalate the spread of this sub-subtype. While evolutionary studies using viral sequences have been instrumental in identifying the HIV epidemic patterns, the origins, and dynamics of the A6 sub-subtype in Poland remain to be explored. We analyzed 1185 HIV-1 A6 pol sequences from Poland, along with 8318 publicly available sequences from other countries. For analyses, phylogenetic tree construction, population dynamics inference, Bayesian analysis, and discrete phylogeographic modeling were employed. Of the introduction events to Poland, 69.94% originated from Ukraine, followed by 29.17% from Russia. Most A6 sequences in Poland (53.16%) formed four large clades, with their introductions spanning 1993-2008. Central and Southern Polish regions significantly influenced migration events. Transmissions among men who have sex with men (MSM) emerged as the dominant risk group for virus circulation, representing 72.92% of migration events. Sequences from migrants were found primarily outside the large clades. Past migration from Ukraine has fueled the spread of the A6 sub-subtype and the current influx of war-displaced people maintains the growing national epidemic.

Female Sex Hormones Upregulate the Replication Activity of HIV-1 Sub-Subtype A6 and CRF02_AG but Not HIV-1 Subtype B

More than 50% of all people living with HIV worldwide are women. Globally, HIV/AIDS is the leading cause of death among women aged 15 to 44. The safe and effective methods of hormonal contraception are an essential component of preventive medical care in order to reduce maternal and infant mortality. However, there is limited knowledge regarding the effect of hormones on the rate of viral replication in HIV infection, especially non-B subtypes. The goal of the present work was to study in vitro how the female hormones β-estradiol and progesterone affect the replication of the HIV-1 subtypes A6, CRF02_AG, and B. The findings show that high doses of hormones enhanced the replication of HIV-1 sub-subtype A6 by an average of 1.75 times and the recombinant variant CRF02_AG by 1.4 times but did not affect the replication of HIV-1 subtype B. No difference was detected in the expression of CCR5 and CXCR4 co-receptors on the cell surface, either in the presence or absence of hormones. However, one of the reasons for the increased viral replication could be the modulated TLRs secretion, as it was found that high doses of estradiol and progesterone upregulated, to varying degrees, the expression of TLR2 and TLR9 genes in the PBMCs of female donors infected with HIV-1 sub-subtype A6.

Clinical Perspective on Human Immunodeficiency Virus Care of Ukrainian War Refugees in Poland

BACKGROUND

The Russian invasion of Ukraine forced migration for safety, protection, and assistance. Poland is the primary sheltering country for Ukrainian refugees, providing support including medical care, which resulted in the rapid ∼15% increase in the number of followed-up people with human immunodeficiency virus (HIV) (PWH) in the country. Here, we present the national experience on HIV care provided for refugees from Ukraine.

METHODS

Clinical, antiretroviral, immunological, and virologic data from 955 Ukrainian PWH entering care in Poland since February 2022 were analyzed. The dataset included both antiretroviral-treated (n = 851) and newly diagnosed (n = 104) patients. In 76 cases, protease/reverse transcriptase/integrase sequencing was performed to identify drug resistance and subtype.

RESULTS

Most (70.05%) of the patients were female, with a predominance of heterosexual (70.3%) transmissions. Anti-hepatitis C antibody and hepatitis B antigen were present in 28.7% and 2.9% of the patients, respectively. A history of tuberculosis was reported in 10.1% of cases. Among previously treated patients, the viral suppression rate was 89.6%; 77.3% of newly HIV diagnosed cases were diagnosed late (with lymphocyte CD4 count <350 cells/μL or AIDS). The A6 variant was observed in 89.0% of sequences. Transmitted mutations in the reverse transcriptase were found in 15.4% treatment-naive cases. Two patients with treatment failure exhibited multiclass drug resistance.

CONCLUSIONS

Migration from Ukraine influences the characteristics of HIV epidemics in Europe, with an increase in the proportion of women and hepatitis C coinfected patients. Antiretroviral treatment efficacy among previously treated refugees was high, with new HIV cases frequently diagnosed late. The A6 subtype was the most common variant.

Detection of Polish cases of highly virulent subtype B of HIV-1 originating in the Netherlands

Infection with the human immunodeficiency virus type 1 (HIV-1) subtype B is most commonly acquired in Poland through men who have sex with men (MSM) comparable to the HIV epidemic in the Netherlands. Following a paper by Chris Wymant et al. on February 4, 2022 in Science on a highly virulent variant of HIV-1 subtype-B (VB-variant) in the Netherlands raised concerns about the possibility of the variant dissemination to other European countries. We aim to report the spread of HIV-1 VB-variant, recently identified in the Netherlands, into other European regions. Subtype B pol gene fragments of protease (P), reverse transcriptase (RT), and integrase (IN) from our laboratory supplemented with publicly available sequences were inferred with VB samples from the Netherlands. For positively clustering samples, clinical observations were compiled. Between May 2009 and August 2014, three cases of VB sequences of Polish origin and one additional from Belgium were identified. Patients presented with elevated viral loads and fast CD4 decline as original characteristics. The mean number of base substitutions per site within the clade versus interclade variability showed a high intragroup sequence similarity, reflecting an ongoing MSM transmission cluster for Polish sequences. The sampling period coincides with the ongoing Dutch VB-variant spread reported between 2003 and 2014. This study informs on phylogenetic descriptions, and clinical symptoms from the rare and emerging VBs placed in Poland. VB is not expanding since 2014 and the Inviduals infected with the VB virus can be treated successfully. Studies on the propagation of novel and potentially virulent virus variants in the undersampled regions add to the understanding of the pan-European HIV-1 transmission dynamics.

Frequency of genotypic factors possibly associated with cabotegravir/rilpivirine failure in antiretroviral treatment-naïve and -experienced HIV-1- infected population

OBJECTIVES

The long-acting injectable (LAI) cabotegravir (CAB) and rilpivirine (RPV) treatment offers important advantages over oral ART (antiretroviral therapy), however baseline factors possibly contributing to the CAB/RPV treatment failure were identified. The purpose of this study was to describe the frequency of virologic factors previously influencing efficacy of this treatment, namely RPV and CAB resistance-associated mutations (RAMs) and A1/A6 subtype among naïve and treatment-experienced HIV-1-infected patients from Poland.

METHODS

The following datasets of HIV-1 sequences were analysed: 4809 protease and reverse transcriptase (PR/RT) sequences obtained from 4649 Polish Caucasian patients (4122 naive and 687 non-naïve) supplemented with integrase (PR/RT/INT) sequences in 1217 cases (942 naïve and 275 non-naïve). Sub-subtypes A were assigned by phylogenetic methods. Major and minor CAB and RPV RAMs were determined according to the IAS-USA 2019 list, while minor RAMs were additionally defined based on the Stanford database algorithm.

RESULTS

Subtype A1/A6 frequency ranged from 6.11% in ART failing cases with PR/RT sequences only, to 15.92% for the PR/RT/INT treatment-naïve dataset, while RPV RAMs were found in up to 5.89% of treatment-naïve and 14.56% of ART failing cases. Regardless treatment history, only <1% sequences had combination of two factors (RPV RAMs and A1/A6 subtype). Furthermore, CAB RAMs were found in 1.27% of treatment-naïve and 14.54% of experienced patients.

CONCLUSIONS

Despite notable frequency of subtype A1/A6 or CAB/RPV RAMs analysed separately, combination of at least two factors previously associated with failure or this treatment is rare. As subtype A1/A6 becomes more common across real-life cohorts continued subtyping and RAM screening will remain of key importance for LAI treatment implementation. Sequence data from this article have been deposited in GenBank under accession numbers: GU906860, GU906864, GU906871-GU906874, JQ305750-JQ305791, KC409134-KC409222, KM057341-KM057362, KM283892-KM284490, KT340108-KT340205, MZ468643-MZ468894, MZ671788-MZ671823, OP298017-OP302727.