Association of vaccination, international travel, public health and social measures with lineage dynamics of SARS-CoV-2.

Continually emerging SARS-CoV-2 variants of concern that can evade immune defenses are driving recurrent epidemic waves of COVID-19 globally. However, the impact of measures to contain the virus and their effect on lineage diversity dynamics are poorly understood. Here, we jointly analyzed international travel, public health and social measures (PHSM), COVID-19 vaccine rollout, SARS-CoV-2 lineage diversity, and the case growth rate (GR) from March 2020 to September 2022 across 63 countries. We showed that despite worldwide vaccine rollout, PHSM are effective in mitigating epidemic waves and lineage diversity. An increase of 10,000 monthly travelers in a single country-to-country route between endemic countries corresponds to a 5.5% (95% CI: 2.9 to 8.2%) rise in local lineage diversity. After accounting for PHSM, natural immunity from previous infections, and waning immunity, we discovered a negative association between the GR of cases and adjusted vaccine coverage (AVC). We also observed a complex relationship between lineage diversity and vaccine rollout. Specifically, we found a significant negative association between lineage diversity and AVC at both low and high levels but not significant at the medium level. Our study deepens the understanding of population immunity and lineage dynamics for future pandemic preparedness and responsiveness.

Bidirectional Movement of Emerging H5N8 Avian Influenza Viruses Between Europe and Asia via Migratory Birds Since Early 2020.

Migratory birds play a critical role in the rapid spread of highly pathogenic avian influenza (HPAI) H5N8 virus clade 2.3.4.4 across Eurasia. Elucidating the timing and pattern of virus transmission is essential therefore for understanding the spatial dissemination of these viruses. In this study, we surveyed >27,000 wild birds in China, tracked the year-round migration patterns of 20 bird species across China since 2006, and generated new HPAI H5N8 virus genomic data. Using this new data set, we investigated the seasonal transmission dynamics of HPAI H5N8 viruses across Eurasia. We found that introductions of HPAI H5N8 viruses to different Eurasian regions were associated with the seasonal migration of wild birds. Moreover, we report a backflow of HPAI H5N8 virus lineages from Europe to Asia, suggesting that Europe acts as both a source and a sink in the global HPAI virus transmission network.

Dispersal dynamics of SARS-CoV-2 lineages during the first epidemic wave in New York City.

During the first phase of the COVID-19 epidemic, New York City rapidly became the epicenter of the pandemic in the United States. While molecular phylogenetic analyses have previously highlighted multiple introductions and a period of cryptic community transmission within New York City, little is known about the circulation of SARS-CoV-2 within and among its boroughs. We here perform phylogeographic investigations to gain insights into the circulation of viral lineages during the first months of the New York City outbreak. Our analyses describe the dispersal dynamics of viral lineages at the state and city levels, illustrating that peripheral samples likely correspond to distinct dispersal events originating from the main metropolitan city areas. In line with the high prevalence recorded in this area, our results highlight the relatively important role of the borough of Queens as a transmission hub associated with higher local circulation and dispersal of viral lineages toward the surrounding boroughs.

A systematic review on global RSV genetic data: Identification of knowledge gaps.

Respiratory syncytial virus (RSV) is a major health problem. A better understanding of the geographical and temporal dynamics of RSV circulation will assist in tracking resistance against therapeutics currently under development. Since 2015, the field of RSV molecular epidemiology has evolved rapidly with around 20-30 published articles per year. The objective of this systematic review is to identify knowledge gaps in recent RSV genetic literature to guide global molecular epidemiology research. We included 78 studies published between 2015 and 2020 describing 12,998 RSV sequences of which 8,233 (63%) have been uploaded to GenBank. Seventeen (22%) studies were performed in low- and middle-income countries (LMICs), and seven (9%) studies sequenced whole-genomes. Although most reported polymorphisms for monoclonal antibodies in clinical development (nirsevimab, MK-1654) have not been tested for resistance in neutralisation essays, known resistance was detected at low levels for the nirsevimab and palivizumab binding site. High resistance was found for the suptavumab binding site. We present the first literature review of an enormous amount of RSV genetic data. The need for global monitoring of RSV molecular epidemiology becomes increasingly important in evaluating the effectiveness of monoclonal antibody candidates as they reach their final stages of clinical development. We have identified the following three knowledge gaps: whole-genome data to study global RSV evolution, data from LMICs and data from global surveillance programs.

Dynamics and Dispersal of Local Human Immunodeficiency Virus Epidemics Within San Diego and Across the San Diego-Tijuana Border.

BACKGROUND: Evolutionary analyses of well-annotated human immunodeficiency virus (HIV) sequence data can provide insights into viral transmission patterns and associated factors. Here, we explored the transmission dynamics of the HIV-1 subtype B epidemic across the San Diego (US) and Tijuana (Mexico) border region to identify factors that could help guide public health policy. METHODS: HIV pol sequences were collected from people with HIV in San Diego County and Tijuana between 1996-2018. A multistep phylogenetic approach was used to characterize the dynamics of spread. The contributions of geospatial factors and HIV risk group to the local dynamics were evaluated. RESULTS: Phylogeographic analyses of the 2034 sequences revealed an important contribution of local transmission in sustaining the epidemic, as well as a complex viral migration network across the region. Geospatial viral dispersal between San Diego communities occurred predominantly among men who have sex with men, with central San Diego being the main source (34.9%) and recipient (39.5%) of migration events. HIV migration was more frequent from San Diego county towards Tijuana than vice versa. Migrations were best explained by the driving time between locations. CONCLUSIONS: The US-Mexico border may not be a major barrier to the spread of HIV, which may stimulate coordinated transnational intervention approaches. Whereas a focus on central San Diego has the potential to avert most spread, the substantial viral migration independent of central San Diego shows that county-wide efforts will be more effective. Combined, this work shows that epidemiological information gleaned from pathogen genomes can uncover mechanisms that underlie sustained spread and, in turn, can be a building block of public health decision-making.

Comparative Circulation Dynamics of the Five Main HIV Types in China.

The HIV epidemic in China accounts for 3% of the global HIV incidence. We compared the patterns and determinants of interprovincial spread of the five most prevalent circulating types. HIV pol sequences sampled across China were used to identify relevant transmission networks of the five most relevant HIV-1 types (B and circulating recombinant forms [CRFs] CRF01_AE, CRF07_BC, CRF08_BC, and CRF55_01B) in China. From these, the dispersal history across provinces was inferred. A generalized linear model (GLM) was used to test the association between migration rates among provinces and several measures of human mobility. A total of 10,707 sequences were collected between 2004 and 2017 across 26 provinces, among which 1,962 are newly reported here. A mean of 18 (minimum and maximum, 1 and 54) independent transmission networks involving up to 17 provinces were identified. Discrete phylogeographic analysis largely recapitulates the documented spread of the HIV types, which in turn, mirrors within-China population migration flows to a large extent. In line with the different spatiotemporal spread dynamics, the identified drivers thereof were also heterogeneous but are consistent with a central role of human mobility. The comparative analysis of the dispersal dynamics of the five main HIV types circulating in China suggests a key role of large population centers and developed transportation infrastructures as hubs of HIV dispersal. This advocates for coordinated public health efforts in addition to local targeted interventions.IMPORTANCE While traditional epidemiological studies are of great interest in describing the dynamics of epidemics, they struggle to fully capture the geospatial dynamics and factors driving the dispersal of pathogens like HIV as they have difficulties capturing linkages between infections. To overcome this, we used a discrete phylogeographic approach coupled to a generalized linear model extension to characterize the dynamics and drivers of the across-province spread of the five main HIV types circulating in China. Our results indicate that large urbanized areas with dense populations and developed transportation infrastructures are facilitators of HIV dispersal throughout China and highlight the need to consider harmonized country-wide public policies to control local HIV epidemics.

Genomic Epidemiology of 2015-2016 Zika Virus Outbreak in Cape Verde.

During 2015-2016, Cape Verde, an island nation off the coast of West Africa, experienced a Zika virus (ZIKV) outbreak involving 7,580 suspected Zika cases and 18 microcephaly cases. Analysis of the complete genomes of 3 ZIKV isolates from the outbreak indicated the strain was of the Asian (not African) lineage. The Cape Verde ZIKV sequences formed a distinct monophylogenetic group and possessed 1-2 (T659A, I756V) unique amino acid changes in the envelope protein. Phylogeographic and serologic evidence support earlier introduction of this lineage into Cape Verde, possibly from northeast Brazil, between June 2014 and August 2015, suggesting cryptic circulation of the virus before the initial wave of cases were detected in October 2015. These findings underscore the utility of genomic-scale epidemiology for outbreak investigations.

Evaluating predictive markers for viral rebound and safety assessment in blood and lumbar fluid during HIV-1 treatment interruption.

BACKGROUND: Validated biomarkers to evaluate HIV-1 cure strategies are currently lacking, therefore requiring analytical treatment interruption (ATI) in study participants. Little is known about the safety of ATI and its long-term impact on patient health. OBJECTIVES: ATI safety was assessed and potential biomarkers predicting viral rebound were evaluated. METHODS: PBMCs, plasma and CSF were collected from 11 HIV-1-positive individuals at four different timepoints during ATI (NCT02641756). Total and integrated HIV-1 DNA, cell-associated (CA) HIV-1 RNA transcripts and restriction factor (RF) expression were measured by PCR-based assays. Markers of neuroinflammation and neuronal injury [neurofilament light chain (NFL) and YKL-40 protein] were measured in CSF. Additionally, neopterin, tryptophan and kynurenine were measured, both in plasma and CSF, as markers of immune activation. RESULTS: Total HIV-1 DNA, integrated HIV-1 DNA and CA viral RNA transcripts did not differ pre- and post-ATI. Similarly, no significant NFL or YKL-40 increases in CSF were observed between baseline and viral rebound. Furthermore, markers of immune activation did not increase during ATI. Interestingly, the RFs SLFN11 and APOBEC3G increased after ATI before viral rebound. Similarly, Tat-Rev transcripts were increased preceding viral rebound after interruption. CONCLUSIONS: ATI did not increase viral reservoir size and it did not reveal signs of increased neuronal injury or inflammation, suggesting that these well-monitored ATIs are safe. Elevation of Tat-Rev transcription and induced expression of the RFs SLFN11 and APOBEC3G after ATI, prior to viral rebound, indicates that these factors could be used as potential biomarkers predicting viral rebound.

Transmission dynamics of re-emerging rabies in domestic dogs of rural China.

Despite ongoing efforts to control transmission, rabies prevention remains a challenge in many developing countries, especially in rural areas of China where re-emerging rabies is under-reported due to a lack of sustained animal surveillance. By taking advantage of detailed genomic and epidemiological data for the re-emerging rabies outbreak in Yunnan Province, China, collected between 1999 and 2015, we reconstruct the demographic and dispersal history of domestic dog rabies virus (RABV) as well as the dynamics of dog-to-dog and dog-to-human transmission. Phylogeographic analyses reveal a lower diffusion coefficient than previously estimated for dog RABV dissemination in northern Africa. Furthermore, epidemiological analyses reveal transmission rates between dogs, as well as between dogs and humans, lower than estimates for Africa. Finally, we show that reconstructed epidemic history of RABV among dogs and the dynamics of rabid dogs are consistent with the recorded human rabies cases. This work illustrates the benefits of combining phylogeographic and epidemic modelling approaches for uncovering the spatiotemporal dynamics of zoonotic diseases, with both approaches providing estimates of key epidemiological parameters.

Increasing importance of European lineages in seeding the hepatitis C virus subtype 1a epidemic in Spain.

BackgroundReducing the burden of the hepatitis C virus (HCV) requires large-scale deployment of intervention programmes, which can be informed by the dynamic pattern of HCV spread. In Spain, ongoing transmission of HCV is mostly fuelled by people who inject drugs (PWID) infected with subtype 1a (HCV1a).AimOur aim was to map how infections spread within and between populations, which could help formulate more effective intervention programmes to halt the HCV1a epidemic in Spain.MethodsEpidemiological links between HCV1a viruses from a convenience sample of 283 patients in Spain, mostly PWID, collected between 2014 and 2016, and 1,317, 1,291 and 1,009 samples collected abroad between 1989 and 2016 were reconstructed using sequences covering the NS3, NS5A and NS5B genes. To efficiently do so, fast maximum likelihood-based tree estimation was coupled to a flexible Bayesian discrete phylogeographic inference method.ResultsThe transmission network structure of the Spanish HCV1a epidemic was shaped by continuous seeding of HCV1a into Spain, almost exclusively from North America and European countries. The latter became increasingly relevant and have dominated in recent times. Export from Spain to other countries in Europe was also strongly supported, although Spain was a net sink for European HCV1a lineages. Spatial reconstructions showed that the epidemic in Spain is diffuse, without large, dominant within-country networks.ConclusionTo boost the effectiveness of local intervention efforts, concerted supra-national strategies to control HCV1a transmission are needed, with a strong focus on the most important drivers of ongoing transmission, i.e. PWID and other high-risk populations.

Discovery and genome characterization of three new Jeilongviruses, a lineage of paramyxoviruses characterized by their unique membrane proteins.

BACKGROUND: In the past decade, many new paramyxoviruses that do not belong to any of the seven established genera in the family Paramyxoviridae have been discovered. Amongst them are J-virus (JPV), Beilong virus (BeiPV) and Tailam virus (TlmPV), three paramyxovirus species found in rodents. Based on their similarities, it has been suggested that these viruses should compose a new genus, tentatively called 'Jeilongvirus'. RESULTS: Here we present the complete genomes of three newly discovered paramyxoviruses, one found in a bank vole (Myodes glareolus) from Slovenia and two in a single, co-infected Rungwe brush-furred rat (Lophuromys machangui) from Mozambique, that represent three new, separate species within the putative genus 'Jeilongvirus'. The genome organization of these viruses is similar to other paramyxoviruses, but like JPV, BeiPV and TlmPV, they possess an additional open reading frame, encoding a transmembrane protein, that is located between the F and G genes. As is the case for all Jeilongviruses, the G genes of the viruses described here are unusually large, and their encoded proteins are characterized by a remarkable amino acid composition pattern that is not seen in other paramyxoviruses, but resembles certain motifs found in Orthopneumovirus G proteins. CONCLUSIONS: The phylogenetic clustering of JPV, BeiPV and TlmPV with the viruses described here, as well as their shared features that set them apart from other paramyxoviruses, provide additional support for the recognition of the genus 'Jeilongvirus'.

Implications of hepatitis C virus subtype 1a migration patterns for virus genetic sequencing policies in Italy.

BACKGROUND: In-depth phylogeographic analysis can reveal migration patterns relevant for public health planning. Here, as a model, we focused on the provenance, in the current Italian HCV subtype 1a epidemic, of the NS3 resistance-associated variant (RAV) Q80K, known to interfere with the action of NS3/4A protease inhibitor simeprevir. HCV1a migration patterns were analysed using Bayesian phylodynamic tools, capitalising on newly generated and publicly available time and geo-referenced NS3 encoding virus genetic sequence data. RESULTS: Our results showed that both immigration and local circulation fuel the current Italian HCV1a epidemic. The United States and European continental lineages dominate import into Italy, with the latter taking the lead from the 1970s onwards. Since similar migration patterns were found for Q80K and other lineages, no clear differentiation of the risk for failing simeprevir can be made between patients based on their migration and travel history. Importantly, since HCV only occasionally recombines, these results are readily transferable to the genetic sequencing policy concerning NS5A RAVs. CONCLUSIONS: The patient migration and travel history cannot be used to target only part of the HCV1a infected population for drug resistance testing before start of antiviral therapy. Consequently, it may be cost-effective to expand genotyping efforts to all HCV1a infected patients eligible for simeprevir-based therapies.

SpreaD3: Interactive Visualization of Spatiotemporal History and Trait Evolutionary Processes.

Model-based phylogenetic reconstructions increasingly consider spatial or phenotypic traits in conjunction with sequence data to study evolutionary processes. Alongside parameter estimation, visualization of ancestral reconstructions represents an integral part of these analyses. Here, we present a complete overhaul of the spatial phylogenetic reconstruction of evolutionary dynamics software, now called SpreaD3 to emphasize the use of data-driven documents, as an analysis and visualization package that primarily complements Bayesian inference in BEAST (http://beast.bio.ed.ac.uk, last accessed 9 May 2016). The integration of JavaScript D3 libraries (www.d3.org, last accessed 9 May 2016) offers novel interactive web-based visualization capacities that are not restricted to spatial traits and extend to any discrete or continuously valued trait for any organism of interest.

Cross-border spread, lineage displacement and evolutionary rate estimation of rabies virus in Yunnan Province, China.

BACKGROUND: Rabies is an important but underestimated threat to public health, with most cases reported in Asia. Since 2000, a new epidemic wave of rabies has emerged in Yunnan Province, southwestern China, which borders three countries in Southeast Asia. METHOD: We estimated gene-specific evolutionary rates for rabies virus using available data in GenBank, then used this information to calibrate the timescale of rabies virus (RABV) spread in Asia. We used 452 publicly available geo-referenced complete nucleoprotein (N) gene sequences, including 52 RABV sequences that were recently generated from samples collected in Yunnan between 2008 and 2012. RESULTS: The RABV N gene evolutionary rate was estimated to be 1.88 × 10-4 (1.37-2.41 × 10-4, 95% Bayesian credible interval, BCI) substitutions per site per year. Phylogenetic reconstructions show that the currently circulating RABV lineages in Yunnan result from at least seven independent introductions (95% BCI: 6-9 introductions) and represent each of the three main Asian RABV lineages, SEA-1, -2 and -3. We find that Yunnan is a sink location for the domestic spread of RABV and connects RABV epidemics in North China, South China, and Southeast Asia. Cross-border spread from southeast Asia (SEA) into South China, and intermixing of the North and South China epidemics is also well supported. The influx of RABV into Yunnan from SEA was not well-supported, likely due to the poor sampling of SEA RABV diversity. We found evidence for a lineage displacement of the Yunnan SEA-2 and -3 lineages by Yunnan SEA-1 strains, and considered whether this could be attributed to fitness differences. CONCLUSION: Overall, our study contributes to a better understanding of the spread of RABV that could facilitate future rabies virus control and prevention efforts.

Distinct Effects of T-705 (Favipiravir) and Ribavirin on Influenza Virus Replication and Viral RNA Synthesis.

T-705 (favipiravir) is a new antiviral agent in advanced clinical development for influenza therapy. It is supposed to act as an alternative substrate for the viral polymerase, causing inhibition of viral RNA synthesis or virus mutagenesis. These mechanisms were also proposed for ribavirin, an established and broad antiviral drug that shares structural similarity with T-705. We here performed a comparative analysis of the effects of T-705 and ribavirin on influenza virus and host cell functions. Influenza virus-infected cell cultures were exposed to T-705 or ribavirin during single or serial virus passaging. The effects on viral RNA synthesis and infectious virus yield were determined and mutations appearing in the viral genome were detected by whole-genome virus sequencing. In addition, the cellular nucleotide pools as well as direct inhibition of the viral polymerase enzyme were quantified. We demonstrate that the anti-influenza virus effect of ribavirin is based on IMP dehydrogenase inhibition, which results in fast and profound GTP depletion and an imbalance in the nucleotide pools. In contrast, T-705 acts as a potent and GTP-competitive inhibitor of the viral polymerase. In infected cells, viral RNA synthesis is completely inhibited by T-705 or ribavirin at ≥50 µM, whereas exposure to lower drug concentrations induces formation of noninfectious particles and accumulation of random point mutations in the viral genome. This mutagenic effect is 2-fold higher for T-705 than for ribavirin. Hence, T-705 and ribavirin both act as purine pseudobases but profoundly differ with regard to the mechanism behind their antiviral and mutagenic effects on influenza virus.

Contribution of Epidemiological Predictors in Unraveling the Phylogeographic History of HIV-1 Subtype C in Brazil.

UNLABELLED: The phylogeographic history of the Brazilian HIV-1 subtype C (HIV-1C) epidemic is still unclear. Previous studies have mainly focused on the capital cities of Brazilian federal states, and the fact that HIV-1C infections increase at a higher rate than subtype B infections in Brazil calls for a better understanding of the process of spatial spread. A comprehensive sequence data set sampled across 22 Brazilian locations was assembled and analyzed. A Bayesian phylogeographic generalized linear model approach was used to reconstruct the spatiotemporal history of HIV-1C in Brazil, considering several potential explanatory predictors of the viral diffusion process. Analyses were performed on several subsampled data sets in order to mitigate potential sample biases. We reveal a central role for the city of Porto Alegre, the capital of the southernmost state, in the Brazilian HIV-1C epidemic (HIV-1C_BR), and the northward expansion of HIV-1C_BR could be linked to source populations with higher HIV-1 burdens and larger proportions of HIV-1C infections. The results presented here bring new insights to the continuing discussion about the HIV-1C epidemic in Brazil and raise an alternative hypothesis for its spatiotemporal history. The current work also highlights how sampling bias can confound phylogeographic analyses and demonstrates the importance of incorporating external information to protect against this. IMPORTANCE: Subtype C is responsible for the largest HIV infection burden worldwide, but our understanding of its transmission dynamics remains incomplete. Brazil witnessed a relatively recent introduction of HIV-1C compared to HIV-1B, but it swiftly spread throughout the south, where it now circulates as the dominant variant. The northward spread has been comparatively slow, and HIV-1B still prevails in that region. While epidemiological data and viral genetic analyses have both independently shed light on the dynamics of spread in isolation, their combination has not yet been explored. Here, we complement publically available sequences and new genetic data from 13 cities with epidemiological data to reconstruct the history of HIV-1C spread in Brazil. The combined approach results in more robust reconstructions and can protect against sampling bias. We found evidence for an alternative view of the HIV-1C spatiotemporal history in Brazil that, contrary to previous explanations, integrates seamlessly with other observational data.

Detection of Natural Resistance-Associated Substitutions by Ion Semiconductor Technology in HCV1b Positive, Direct-Acting Antiviral Agents-Naïve Patients.

Naturally occurring resistance-associated substitutions (RASs) can negatively impact the response to direct-acting antivirals (DAAs) agents-based therapies for hepatitis C virus (HCV) infection. Herein, we set out to characterize the RASs in the HCV1b genome from serum samples of DAA-naïve patients in the context of the SINERGIE (South Italian Network for Rational Guidelines and International Epidemiology, 2014) project. We deep-sequenced the NS3/4A protease region of the viral population using the Ion Torrent Personal Genome Machine, and patient-specific majority rule consensus sequence summaries were constructed with a combination of freely available next generation sequencing data analysis software. We detected NS3/4A protease major and minor variants associated with resistance to boceprevir (V36L), telaprevir (V36L, I132V), simeprevir (V36L), and grazoprevir (V36L, V170I). Furthermore, we sequenced part of HCV NS5B polymerase using Sanger-sequencing and detected a natural RAS for dasabuvir (C316N). This mutation could be important for treatment strategies in cases of previous therapy failure.

The genealogical population dynamics of HIV-1 in a large transmission chain: bridging within and among host evolutionary rates.

Transmission lies at the interface of human immunodeficiency virus type 1 (HIV-1) evolution within and among hosts and separates distinct selective pressures that impose differences in both the mode of diversification and the tempo of evolution. In the absence of comprehensive direct comparative analyses of the evolutionary processes at different biological scales, our understanding of how fast within-host HIV-1 evolutionary rates translate to lower rates at the between host level remains incomplete. Here, we address this by analyzing pol and env data from a large HIV-1 subtype C transmission chain for which both the timing and the direction is known for most transmission events. To this purpose, we develop a new transmission model in a Bayesian genealogical inference framework and demonstrate how to constrain the viral evolutionary history to be compatible with the transmission history while simultaneously inferring the within-host evolutionary and population dynamics. We show that accommodating a transmission bottleneck affords the best fit our data, but the sparse within-host HIV-1 sampling prevents accurate quantification of the concomitant loss in genetic diversity. We draw inference under the transmission model to estimate HIV-1 evolutionary rates among epidemiologically-related patients and demonstrate that they lie in between fast intra-host rates and lower rates among epidemiologically unrelated individuals infected with HIV subtype C. Using a new molecular clock approach, we quantify and find support for a lower evolutionary rate along branches that accommodate a transmission event or branches that represent the entire backbone of transmitted lineages in our transmission history. Finally, we recover the rate differences at the different biological scales for both synonymous and non-synonymous substitution rates, which is only compatible with the 'store and retrieve' hypothesis positing that viruses stored early in latently infected cells preferentially transmit or establish new infections upon reactivation.

Mapping Transmission Dynamics and Drug Resistance Surveillance in the Cyprus HIV-1 Epidemic (2017-2021).

The human immunodeficiency virus type 1 (HIV-1) epidemic has been a major public health threat on a global scale since the early 1980s. Despite the introduction of combination antiretroviral therapy (cART), the incidence of new HIV-1 infections continues to rise in some regions around the world. Thus, with the continuous transmission of HIV-1 and the lack of a cure, it is imperative for molecular epidemiological studies to be performed, to monitor the infection and ultimately be able to control the spread of this virus. This work provides a comprehensive molecular epidemiological analysis of the HIV-1 infection in Cyprus, through examining 305 HIV-1 sequences collected between 9 March 2017 and 14 October 2021. Employing advanced statistical and bioinformatic techniques, the research delved deeply into understanding the transmission dynamics of the HIV-1 epidemic in Cyprus, as well as the monitoring of HIV-1's genetic diversity and the surveillance of transmitted drug resistance. The characterization of Cyprus's HIV-1 epidemic revealed a diverse landscape, comprising 21 HIV-1 group M pure subtypes and circulating recombinant forms (CRFs), alongside numerous uncharacterized recombinant strains. Subtypes A1 and B emerged as the most prevalent strains, followed by CRF02_AG. The findings of this study also revealed high levels of transmitted drug resistance (TDR) patterns, raising concerns for the efficacy of cART. The demographic profiles of individuals involved in HIV-1 transmission underscored the disproportionate burden borne by young to middle-aged Cypriot males, particularly those in the MSM community, who reported contracting the virus in Cyprus. An assessment of the spatiotemporal evolutionary dynamics illustrated the global interconnectedness of HIV-1 transmission networks, implicating five continents in the dissemination of strains within Cyprus: Europe, Africa, Asia, North America, and Oceania. Overall, this study advances the comprehension of the HIV-1 epidemic in Cyprus and highlights the importance of understanding HIV-1's transmission dynamics through continuous surveillance efforts. Furthermore, this work emphasizes the critical role of state-of-the-art bioinformatics analyses in addressing the challenges posed by HIV-1 transmission globally, laying the groundwork for public health interventions aimed at curbing its spread and improving patient outcomes.