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.

Circulation of Human Immunodeficiency Virus 1 A6 Variant in the Eastern Border of the European Union-Dynamics of the Virus Transmissions Between Poland and Ukraine.

BACKGROUND: The human immunodeficiency virus (HIV) type 1 A6 variant is dominating in high-prevalence Eastern European countries, with increasing prevalence over the remaining regions of Europe. The recent war in Ukraine may contribute to further introductions of this A6 lineage. Our aim was to model the transmission dynamics of the HIV-1 A6 variant between Poland and Ukraine. METHODS: HIV-1 A6 partial pol sequences originating from Poland (n = 1185) and Ukraine (n = 653) were combined with publicly available sequences (n = 7675) from 37 other countries. We used maximum likelihood-based tree estimation followed by a bayesian inference strategy to characterize the putative transmission clades. Asymmetric discrete phylogeographic analysis was used to identify the best-supported virus migration events across administrative regions of Poland and Ukraine. RESULTS: We identified 206 clades (n = 1362 sequences) circulating in Poland or Ukraine (63 binational clades, 79 exclusively Polish, and 64 exclusively Ukrainian). Cross-border migrations were almost exclusively unidirectional (from Ukraine to Poland, 99.4%), mainly from Eastern and Southern Ukraine (Donetsk, 49.7%; Odesa, 17.6% regions) to the Central (Masovian, 67.3%; Lodz, 18.2%) and West Pomeranian (10.1%) districts of Poland. The primary sources of viral dispersal were the Eastern regions of Ukraine, long affected by armed conflict, and large population centers in Poland. CONCLUSIONS: The Polish outbreak of the A6 epidemic was fueled by complex viral migration patterns across the country, together with cross-border transmissions from Ukraine. There is an urgent need to include war-displaced people in the national HIV prevention and treatment programs to reduce the further spread of transmission networks.

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.

Remdesivir Reduces Mortality in Hemato-Oncology Patients with COVID-19.

Introduction: Remdesivir is the first agent with proven clinical efficacy against coronavirus disease 2019 (COVID-19); however, its benefit is associated with early use, and its efficacy has been poorly studied in patients with hemato-oncological diseases, who have an increased risk of a severe course of infection. This study aimed to assess the effects of remdesivir on mortality, mechanical ventilation, and the duration of hospitalization in both the general population and in patients with hemato-oncological diseases. Materials and Methods: Longitudinal data for 4287 patients with confirmed COVID-19 were analyzed, including a subset of 200 individuals with hemato-oncological diseases. In total, 1285 (30.0%) patients received remdesivir, while the remaining patients were treated with other methods. Survival statistics for the 14- and 30-day observation time points were calculated using non-parametric and multivariate Cox models. Results: Mortality for the 14- and 30-day observation time points was notably lower among patients receiving remdesivir (7.2% vs 11.6%, p < 0.001 and 12.7% vs 16.0, p = 0.005, respectively); however, in multivariate models adjusted for age, sex, lung involvement, and lactate dehydrogenase and interleukin-6 levels, the administration of remdesivir did not reduce patient mortality at either the 14-day or 30-day time points. Among patients with haemato-oncological disease, significant survival benefit was observed at 14 and 30 days for patients treated with remdesivir (11.3% vs.16.7% and 24.2% vs 26.1%, respectively; p < 0.001). A favorable effect of remdesivir was also noted for the 14-day time point in multivariate survival analysis (HR:4.03 [95% confidence interval:1.37-11.88]; p = 0.01). Conclusion: Remdesivir significantly reduced the early mortality rate in COVID-19 patients with comorbid hemato-oncological disease, which emphasizes the need to administer this agent to immunosuppressed patients.

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.

Factors Influencing Immune Restoration in People Living with HIV/AIDS.

Introduction: Immune restoration is a key clinical aspect that is pursued in the care of human immunodeficiency virus (HIV)-infected patients. Despite effective antiretroviral treatment and undetectable viremia, immune recovery is often incomplete. Materials and methods: Data from 311 Caucasian patients were collected. SNP in CCR2(rs1799864), CX3CR1(rs3732378), HLAC-35(rs9264942), and CCR5(promoter, rs1799988); a 32bp deletion(Δ32) in CCR5; and HLA-B*5701 genotypes were correlated with clinical data and selected endpoints. Kaplan−Meier and Cox proportional hazards models were used to analyze the effects of genetic factors over time. Results: For HLA-B*5701, the effect on the CD4+/CD8+ >0.8 cell ratio was lost within 48 months (HR = 2.04, 95% CI: 1.04−4.03), and the effect on the CD4+ cell count >500 cells/µL was lost within 12 months (HR = 2.12, CI: 1.11−4.04). The effect of CCR2 GG on the CD4+/CD8+ >0.8 cell ratio was lost within 36 months (HR = 1.7, CI: 1.05−2.75). For CCR5 wt/Δ32, the effect on the CD4+/CD8+ >1.0 cell ratio was lost within 24 months (HR = 2.0, CI: 1.08−3.69), and the effect on the CD4+ >800 cells/µL cell count was lost within 18 months (HR = 1.98, CI: 1.14−4.73). Conclusions: Selected genetic polymorphisms, namely CCR2 GG and CCR5 Δ32, and the presence of the HLA-B*5701 allele positively influenced immune restoration in cART-treated patients with HIV/AIDS.

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

The occurrence of HIV-1 subtypes differs worldwide and within Europe, with non-B variants mainly found across different exposure groups. In this study, we investigated the distribution and temporal trends in HIV-1 subtype variability across Poland between 2015 and 2019. Sequences of the pol gene fragment from 2518 individuals were used for the analysis of subtype prevalence. Subtype B was dominant (n = 2163, 85.90%). The proportion of subtype B-infected individuals decreased significantly, from 89.3% in 2015 to 80.3% in 2019. This was related to the increasing number of subtype A infections. In 355 (14.10%) sequences, non-B variants were identified. In 65 (2.58%) samples, recombinant forms (RFs) were noted. Unique recombinant forms (URFs) were found in 30 (1.19%) sequences. Three A/B recombinant clusters were identified of which two were A6/B mosaic viruses not previously described. Non-B clades were significantly more common among females (n = 81, 22.8%, p = 0.001) and heterosexually infected individuals (n = 45, 32.4%, p = 0.0031). The predominance of subtype B is evident, but the variability of HIV-1 in Poland is notable. Almost half of RFs (n = 65, 2.58%) was comprised of URFs (n = 30, 1.19%); thus those forms were common in the analyzed population. Hence, molecular surveillance of identified variants ensures recognition of HIV-1 evolution in Poland.

HCV resistance-associated substitutions following direct-acting antiviral therapy failure - Real-life data from Poland.

BACKGROUND: This study analysed the NS3 and NS5A mutation frequencies, persistence and drug susceptibility in a cohort of real-life patients, with failed hepatitis C virus (HCV) therapy following directly acting antiviral (DAA) treatment. METHODS: NS3/NS5A Sanger sequences from 105 patients infected with HCV genotype (G) 1a (6,5.7%), G1b (94,89.5%), G3a (4,3.8%), and G4 (1,1.0%) post DAA treatment failure were analysed. NS3 and NS5A resistance-associated substitutions (RASs) were identified using the geno2pheno algorithm and associated with clinical variables. Time trends were examined using logistic regression. RESULTS: NS5A RAS were found in 87.9% of sequences derived from patients exposed to this class of agents, whereas NS3 RAS was found in 59.1% of HCV protease-exposed subjects. The frequency of the NS3 RAS increased with fibrosis stage, from 40.0% among F0/F1 individuals to 81.8% among patients with liver cirrhosis (F4, p = 0.094). NS5A mutation frequencies were 7.6% for 28A/V/M, 10.6% for 30 K/Q/R, 42.4% for 31I/F/M/V, and 75.8% for 93H. For NS3, the most common RASs were 56F-23.7%, 168A/E/I/Y/T/V-14.0%, and 117H-5.4%. Susceptibility to glecaprevir/pibrentasvir, velpatasvir/voxlaprevir, and elbasvir/grazoprevir was retained in 92.9%, 43.4%, and, 25.3% of patients, respectively. The frequency of NS3 RAS decreased with time elapsed from failure to sampling (p = 0.034 for trend). NS5A RAS frequency remained stable over the 24-months. CONCLUSIONS: Following DAA treatment failure, NS5A and NS3 RASs were common with increasing frequency among patients with advanced liver disease. In most cases, despite the presence of RASs, susceptibility to DAA combinations with higher genetic barrier was retained.

Transmitted HIV drug resistance and subtype patterns among blood donors in Poland.

Surveillance on the HIV molecular variability, risk of drug resistance transmission and evolution of novel viral variants among blood donors remains an understudied aspect of hemovigilance. This nationwide study analyses patterns of HIV diversity and transmitted resistance mutations. Study included 185 samples from the first time and repeat blood donors with HIV infection identified by molecular assay. HIV protease, reverse transcriptase and integrase were sequenced using population methods. Drug resistance mutation (DRM) patterns were analyzed based on the Stanford Interpretation Algorithm and standardized lists of transmitted mutations. Phylogeny was used to investigate subtyping, clustering and recombination patterns. HIV-1 subtype B (89.2%) followed by subtype A6 (7.6%) were predominant, while in three (1.6%) cases, novel recombinant B/A6 variants were identified. Non-B variants were more common among repeat donors (14.5%) compared to the first time ones (1.8%), p = 0.011, with higher frequency (9.9%) of A6 variant in the repeat donor group, p = 0.04. Major NRTI DRMs were observed in 3.8%, NNRTI and PI in 0.6% and INSTI 1.1% of cases. Additionally, E157Q polymorphism was observed in 9.8% and L74I in 11.5% of integrase sequences. Transmission of drug resistance among blood donors remains infrequent. Subtype patters increase in complexity with emergence of novel intersubtype A6B recombinants.

Infection with HIV-1 subtype D adversely affects the live expectancy independently of antiretroviral drug use.

INTRODUCTION: HIV-1 subtypes have been associated with less favourable clinical profiles, differences in disease progression and higher risk of neurocognitive deficit. In this study we aimed to analyse the long term survival disparities between patients infected with the most common HIV-1 variants observed in Poland. METHODS: For the study data from 518 Caucasian non-immigrant patients of Polish origin infected with divergent HIV subtypes and variants [subtype A (n = 35, 6.8%), subtype B (n = 386, 74.5%), subtype C (n = 13, 2.5%), subtype D (n = 58, 11.19%) or other non-A,B,C,D (n = 26, 5.01%)variants] were analysed. Subtyping was performed using the partial pol (reverse transcriptase and protease) sequencing. HIV variant was coupled with clinical, virologic and survival data censored at 20 years of observation. Overall survival and on antiretroviral treatment survival was analysed using Kaplan-Meyer as well as unadjusted and multivariate Cox proportional hazards models. RESULTS: Significantly higher mortality was observed among subtype D (28.8%) infected subjects compared to subtype B (11.7%, p = 0.0004). Increased risk of death among subtype D cases remained significant when cART treated individuals were analysed, with on-treatment mortality of 26.9% for subtype D (p = 0.006) compared to 10.73% in subtype B infected cases. Kaplan-Meyer survival estimates differed significantly across all investigated HIV-1 variant groups when overall 20 year mortality was analysed (log rank p = 0.029), being non-significant for the cART treated group. In multivariate model of overall 20 year survival, adjusted for age at diagnosis, gender, HCV and AIDS status, lymphocyte CD4 count, transmission route and HIV viral load, only age and subtype D were independently associated with higher likelihood of death [HR: 1.08 (95%CI: 1.03-1.14, p = 0.002) and HR: 7.91 (95%CI:2.33-26.86), p < 0.001, respectively]. In the on-treatment (cART) multivariate model of 20 year survival adjusted for the same parameters only subtype D remained as the independent factor associated with higher mortality risk [HR: 4.24 (95%CI:1.31-13.7), p = 0.02]. CONCLUSIONS: Subtype D has an independent deleterious effect of survival, even in the setting of antiretroviral treatment. Observed effect indicated higher clinical vigilance for patients infected with this subtype even after long time of stable antiretroviral treatment.

Low prevalence of doravirine-associated resistance mutations among polish human immunodeficiency-1 (HIV-1)-infected patients.

INTRODUCTION: Doravirine (DOR) is a novel non-nucleoside reverse transcriptase inhibitor (NNRTI) that retains activity against common NNRTI resistance mutations. In this study, we aimed to investigate the prevalence of DOR resistance mutations compared with that of resistance mutations for other NNRTIs among HIV-1-infected treatment-experienced and -naïve patients from Poland. METHODS: Resistance to DOR and other NNRTIs was assessed in two datasets: 1760 antiretroviral treatment-naïve HIV-1 patients and 200 treatment-experienced patients. All 1960 sequences were derived from the patients using bulk sequencing. For resistance analyses, Stanford HIV drug resistance database scores were used. RESULTS: Overall, DOR resistance was present in 32 patients (1.62%), of whom 13 (0.74%) were naïve and 19 (9.50%) were treatment-experienced. The most common DOR resistance mutations observed among the naïve patients were A98G and K101E (0.2% each), and those among cART-experienced patients were L100I (2.0%), K101E, V108I, H221Y, and P225H (1.5% each). Furthermore, among the naïve patients, less common resistance to DOR (0.7%) compared with that to nevirapine (NVP) (2.1%; p = 0.0013) and rilpivirine (5.40%; p < 0.0001) was observed. For sequences obtained from treatment-experienced patients, the frequency of resistance to DOR (9.5%) was lower than that for efavirenz (25.5%; p < 0.0001) and NVP (26.0%; p < 0.0001). CONCLUSIONS: The frequency of transmitted drug resistance to DOR is low, allowing for effective treatment of antiretroviral treatment-naïve patients and rapid treatment initiation. In cART-experienced patients, this agent remains an attractive NNRTI option with a higher genetic barrier to resistance.