Genomic Epidemiology and Lineage Dynamics of SARS-CoV-2 in Bulgaria: Insights from a Three-Year Pandemic Analysis.

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has brought about significant challenges worldwide. In this study, we present a comprehensive analysis of the genomic epidemiology and lineage dynamics of SARS-CoV-2 in Bulgaria over a three-year period. Through extensive genomic sequencing and data analysis, we investigated the evolution of the virus, the emergence of variants of concern (VOCs), and their impact on the country's pandemic trajectory. We also assessed the relationship between viral diversity and COVID-19 morbidity and mortality in Bulgaria. Our findings shed light on the temporal and spatial distribution of SARS-CoV-2 lineages and provide crucial insights into the dynamics of the pandemic in the country. The interplay between international travel and viral transmission plays a significant role in the emergence and dissemination of different SARS-CoV-2 variants. The observed proportions of exportation to various continents provide insights into the potential pathways through which these lineages spread globally. Understanding the genomic epidemiology of SARS-CoV-2 in Bulgaria is essential for formulating targeted public health strategies, enhancing vaccination efforts, and effectively managing future outbreaks.

Early Detection of the Recombinant SARS-CoV-2 XAN Variant in Bulgaria: Initial Genomic Insights into Yet Another Piece of the Growing Puzzle of Recombinant Clades.

The first recombinant SARS-CoV-2 variants were identified in 2022, causing public health concerns. The importance of recombinant variants has increased especially since the WHO designated the recombinant variant XBB and its lineages as subvariants that require monitoring on 20 November 2022. In this study, we provide the first insights into the new SARS-CoV-2 variant named XAN, a recombinant composed of Omicron sub-lineages BA.2 and BA.5. To our knowledge, this is the first report on the recombinant SARS-CoV-2 XAN variant identified in Bulgaria.

Transmitted HIV Drug Resistance in Bulgaria Occurs in Clusters of Individuals from Different Transmission Groups and Various Subtypes (2012-2020).

Transmitted HIV drug resistance in Bulgaria was first reported in 2015 using data from 1988-2011. We determined the prevalence of surveillance drug resistance mutations (SDRMs) and HIV-1 genetic diversity in Bulgaria during 2012-2020 using polymerase sequences from 1053 of 2010 (52.4%) antiretroviral therapy (ART)-naive individuals. Sequences were analyzed for DRM using the WHO HIV SDRM list implemented in the calculated population resistance tool at Stanford University. Genetic diversity was inferred using automated subtyping tools and phylogenetics. Cluster detection and characterization was performed using MicrobeTrace. The overall rate of SDRMs was 5.7% (60/1053), with 2.2% having resistance to nucleoside reverse transcriptase inhibitors (NRTIs), 1.8% to non-nucleoside reverse transcriptase inhibitors (NNRTIs), 2.1% to protease inhibitors (PIs), and 0.4% with dual-class SDRMs. We found high HIV-1 diversity, with the majority being subtype B (60.4%), followed by F1 (6.9%), CRF02_AG (5.2%), A1 (3.7%), CRF12_BF (0.8%), and other subtypes and recombinant forms (23%). Most (34/60, 56.7%) of the SDRMs were present in transmission clusters of different subtypes composed mostly of male-to-male sexual contact (MMSC), including a 14-member cluster of subtype B sequences from 12 MMSC and two males reporting heterosexual contact; 13 had the L90M PI mutation and one had the T215S NRTI SDRM. We found a low SDRM prevalence amid high HIV-1 diversity among ART-naive patients in Bulgaria during 2012-2020. The majority of SDRMs were found in transmission clusters containing MMSC, indicative of onward spread of SDRM in drug-naive individuals. Our study provides valuable information on the transmission dynamics of HIV drug resistance in the context of high genetic diversity in Bulgaria, for the development of enhanced prevention strategies to end the epidemic.

Assessment of hepatitis E seropositivity among HIV-infected patients in Bulgaria.

It is debatable whether HIV-infected patients are at greater risk for hepatitis E virus (HEV) infection compared with healthy subjects. The reported anti-HEV seroprevalence among different groups in Bulgaria varied from 9.04% to 25.9%, but the information regarding the HIV population is still missing. The aim of the present study was to evaluate hepatitis E seroprevalence among HIV-infected patients in Bulgaria and to analyze demographic and immunological factors associated with HEV infection. Serum samples of 312 HIV-infected patients were analyzed retrospectively. Age, sex, residence and laboratory markers for HEV, HBV, HCV and HIV infection, and lymphocytes subpopulations were collected for all patients. None of the tested samples were positive for HEV RNA. HEV seroprevalence among HIV-infected patients was 10.9%. Males were more affected with the highest prevalence of positivity in the age group > 30 to ≤ 40 years. The documented HIV transmission routes in HIV/HEV co-infected group were heterosexual, homosexual, intravenous drug use (IDU), and vertical with predominace of the heterosexual route (z = 0.2; p = 0.804). There was a statistically significant trend of HIV mixed infection with routes of HIV transmission other than homosexual - heterosexual in HIV/HEV group and injection drug use in HIV/HBV/HCV co-infected group. The route of HIV transmission, in contexts of patients' behavior, was associated with HEV prevalence among HIV-infected patients.

Assessment of hepatitis E seropositivity among HIV-infected patients in Bulgaria

Abstract It is debatable whether HIV-infected patients are at greater risk for hepatitis E virus (HEV) infection compared with healthy subjects. The reported anti-HEV seroprevalence among different groups in Bulgaria varied from 9.04% to 25.9%, but the information regarding the HIV population is still missing. The aim of the present study was to evaluate hepatitis E seroprevalence among HIV-infected patients in Bulgaria and to analyze demographic and immunological factors associated with HEV infection. Serum samples of 312 HIV-infected patients were analyzed retrospectively. Age, sex, residence and laboratory markers for HEV, HBV, HCV and HIV infection, and lymphocytes subpopulations were collected for all patients. None of the tested samples were positive for HEV RNA. HEV seroprevalence among HIV-infected patients was 10.9%. Males were more affected with the highest prevalence of positivity in the age group > 30 to ≤ 40 years. The documented HIV transmission routes in HIV/HEV co-infected group were heterosexual, homosexual, intravenous drug use (IDU), and vertical with predominace of the heterosexual route (z = 0.2; p = 0.804). There was a statistically significant trend of HIV mixed infection with routes of HIV transmission other than homosexual - heterosexual in HIV/HEV group and injection drug use in HIV/HBV/HCV co-infected group. The route of HIV transmission, in contexts of patients' behavior, was associated with HEV prevalence among HIV-infected patients.

Molecular Epidemiological Analysis of the Origin and Transmission Dynamics of the HIV-1 CRF01_AE Sub-Epidemic in Bulgaria.

HIV-1 subtype CRF01_AE is the second most predominant strain in Bulgaria, yet little is known about the molecular epidemiology of its origin and transmissibility. We used a phylodynamics approach to better understand this sub-epidemic by analyzing 270 HIV-1 polymerase (pol) sequences collected from persons diagnosed with HIV/AIDS between 1995 and 2019. Using network analyses at a 1.5% genetic distance threshold (d), we found a large 154-member outbreak cluster composed mostly of persons who inject drugs (PWID) that were predominantly men. At d = 0.5%, which was used to identify more recent transmission, the large cluster dissociated into three clusters of 18, 12, and 7 members, respectively, five dyads, and 107 singletons. Phylogenetic analysis of the Bulgarian sequences with publicly available global sequences showed that CRF01_AE likely originated from multiple Asian countries, with Vietnam as the likely source of the outbreak cluster between 1988 and 1990. Our findings indicate that CRF01_AE was introduced into Bulgaria multiple times since 1988, and infections then rapidly spread among PWID locally with bridging to other risk groups and countries. CRF01_AE continues to spread in Bulgaria as evidenced by the more recent large clusters identified at d = 0.5%, highlighting the importance of public health prevention efforts in the PWID communities.

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.