The majority of SARS-CoV-2 genomes obtained during the pandemic were derived by amplifying overlapping windows of the genome (‘tiled amplicons’), reconstructing their sequences and fitting them together. This leads to systematic errors in genomes unless the software is both aware of the amplicon scheme and of the error modes of amplicon sequencing. Additionally, over time, amplicon schemes need to be updated as new mutations in the virus interfere with the primer binding sites at the end of amplicons. Thus, waves of variants swept the world during the pandemic and were followed by waves of systematic errors in the genomes, which had significant impacts on the inferred phylogenetic tree. Here we reconstruct the genomes from all public data as of June 2024 using an assembly tool called Viridian ( https://github.com/iqbal-lab-org/viridian ), developed to rigorously process amplicon sequence data. With these high-quality consensus sequences we provide a global phylogenetic tree of 4,471,579 samples, viewable at https://viridian.taxonium.org . We provide simulation and empirical validation of the methodology, and quantify the improvement in the phylogeny.

Objective: To identify predictors of HIV acquisition in Botswana. Design: We applied machine learning approaches to identify HIV risk predictors using existing data from a large, well characterized HIV incidence cohort. Methods: We applied machine learning (randomForestSRC) to analyze data from a large population-based HIV incidence cohort enrolled in a cluster-randomized HIV prevention trial in 30 communities across Botswana. We sought to identify the most important risk factors for HIV acquisition, starting with 110 potential predictors. Results: During a median 29-month follow-up of 8551 HIV-negative adults, 147 (1.7%) acquired HIV. Our machine learning analysis found that for females, the most important variables for predicting HIV acquisition were the use of injectable hormonal contraception, frequency of sex in the prior 3 months with the most recent partner and residing in a community with HIV prevalence of 29% or higher. For the small proportion (0.3%) of females who had all three risk factors, their estimated probability of acquiring HIV during 29 months of follow-up was 34% (approximate annual incidence of 14%). For males, nonlong-term relationships with the most recent partner and community HIV prevalence of 34% or higher were the most important HIV risk predictors. The 6% of males who had both risk factors had a 5.1% probability of acquiring HIV during the follow-up period (approximate annual incidence of 2.1%). Conclusion: Machine learning approaches allowed us to analyze a large number of variables to efficiently identify key factors strongly predictive of HIV risk. These factors could help target HIV prevention interventions in Botswana. Clinical Trials Registration: NCT01965470

West Nile virus (WNV) is a priority pathogen that poses a high risk for public health emergencies of global concern. Although WNV is endemic to Africa, only few (n=63) whole genomic sequences are available from the continent. In this Review, we examined the status of the molecular testing and genomic sequencing of WNV across Africa and mapped its global spatiotemporal spread. WNV has been detected in 39 African countries, the Canary Islands, and Réunion Island. Although publications, including those with molecular data, originated from 24 of these countries, genomic sequences were available from only 16 countries. Our analysis identified regions with detected viral circulation but without molecular surveillance. The current literature has substantial knowledge gaps in terms of the disease burden, molecular epidemiology, and distribution of WNV in Africa. Addressing these gaps requires an integrated One Health surveillance approach, which is challenging to establish. We propose three key surveillance needs that could improve the current understanding of the WNV disease burden in Africa, to strengthen the global public health response to this vector-borne disease.

Background West Nile Virus (WNV) is a zoonotic flavivirus of significant One Health relevance and is classified as a priority pathogen with a high-risk of causing public health emergencies of global concern. WNV is endemic to Africa; however, the availability of genomic sequences from the continent remains limited. Methods We review the extent of polymerase chain reaction testing and genomic sequencing of WNV conducted across Africa. Using phylogeographic methods, we map the spatiotemporal spread of the virus across the continent and globally. Findings Our study shows that WNV has been detected in 39 African countries (including Comoros, Seychelles, and Mauritius), the Canary Islands, and Reunion Island. Publications including molecular data originate from 24 countries; however, genomic sequences are publicly available for only 16 countries. We identify regions with detected viral circulation but lacking molecular surveillance. Further, we list such regions that overlap with Key Biodiversity Areas (sites harbouring significant bird diversity) as they may host high viral circulation, and high human population density that may be susceptible to spillover. Interpretation We recognise significant knowledge gaps on the true disease burden, molecular epidemiology, and distribution of WNV in Africa. Addressing these gaps requires an integrated One Health surveillance approach which is challenging to establish. We propose three key surveillance needs as potential starting points to improve our understanding of the virus in Africa to strengthen the global public health response to this disease. ### Competing Interest Statement The authors have declared no competing interest.

Botswana, like the rest of the world, has been significantly impacted by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In December 2022, we detected a monophyletic cluster of genomes comprising a sublineage of the Omicron variant of concern (VOC) designated as B.1.1.529.5.3.1.1.1.1.1.1.74.1 (alias FN.1, clade 22E). These genomes were sourced from both epidemiologically linked and unlinked samples collected in three close locations within the district of Greater Gaborone. In this study, we assessed the worldwide prevalence of the FN.1 lineage, evaluated its mutational profile, and conducted a phylogeographic analysis to reveal its global dispersal dynamics. Among approximately 16 million publicly available SARS-CoV-2 sequences generated by 30 September 2023, only 87 were of the FN.1 lineage, including 22 from Botswana, 6 from South Africa, and 59 from the UK. The estimated time to the most recent common ancestor of the 87 FN.1 sequences was 22 October 2022 [95% highest posterior density: 2 September 2022-24 November 2022], with the earliest of the 22 Botswana sequences having been sampled on 7 December 2022. Discrete trait reconstruction of FN.1 identified Botswana as the most probable place of origin. The FN.1 lineage is derived from the BQ.1.1 lineage and carries two missense variants in the spike protein, S:K182E in NTD and S:T478R in RDB. Among the over 90 SARS-CoV-2 lineages circulating in Botswana between September 2020 and July 2023, FN.1 was most closely related to BQ.1.1.74 based on maximum likelihood phylogenetic inference, differing only by the S:K182E mutation found in FN.1. Given the early detection of numerous novel variants from Botswana and its neighbouring countries, our study underscores the necessity of continuous surveillance to monitor the emergence of potential VOCs, integrating molecular and spatial data to identify dissemination patterns enhancing preparedness efforts.

Objectives: We characterized the molecular and mutational landscape of SARS-CoV-2 recombinant strains in Botswana. Methods: We performed genomic, phylogenetic, and immunoinformatic analyses of 5254 near-complete genomes from 2020 to 2023. We assessed the presence of mutations of interested (MutOI) that may be associated with immune escape in silico . Results: We observed a few recombinant strains in Botswana, with the majority being descendants of Omicron (XBB* ), except for XV and XM. Most recombinant sequences corresponded to transmission clusters. Most recombination events occurred within the receptor-binding domain (RDB) of the spike (S) protein. We identified 16 MutOI among different proteins, with the majority occurring at a very low global prevalence ( < 4.8 x 10(- 5 )). We also observed S:Q474K, a MutOI in the RBD, that was predicted to escape human leukocyte antigen class I-mediated immune responses. Molecular surveillance is vital to inform early detection and response to potential variants with heightened immune and vaccine breakthrough properties. Conclusions: These results underscore the need for continued molecular surveillance to map the evolutionary landscape of SARS-CoV-2.

Introduction: The ‘sterile womb paradigm’ is currently under debate and the advent of next generation 16S rRNA gene sequencing is driving the characterization of microbes associated with the amniotic cavity during pregnancy. Objective: To characterize the bacterial diversity in placentas from preterm and term births using next generation 16S rRNA gene sequencing in association with adverse pregnancy outcomes and histopathology studies. Methods: In this prospective study, placentas were collected consecutively from patients attending a public tertiary referral hospital in South Africa, delivering preterm (n=42; 28-34 weeks gestational age) and term (n=20; >37 weeks gestational age). Placentas underwent histopathology tests and next generation 16S rRNA gene sequencing. Results: Analysis of microbial diversity in the placenta showed a significantly higher alpha diversity in term placentas compared to preterm placentas (P=0.000075), in placentas with acute chorioamnionitis compared to placentas without acute chorioamnionitis (P=0.0061), between HIV negative term births and both HIV negative and HIV positive preterm births (P=0.0118 and P=0.0008), respectively. Beta diversity was significantly different between preterm and term births (unweighted UniFrac distance, P=0.003996; Jaccard distance, P=0.03696) and Escherichia/Shigella, Shuttleworthia, Anaeroglobus and Megasphaera were differentially expressed. Conclusion: This is the first South African study to characterize the bacterial diversity in placentas from complicated pregnancies using next generation 16S rRNA gene sequencing in conjunction with full placental histology. Microbial diversity differs between preterm and term placentas where HIV may act as a cofactor associated with decreased bacterial alpha diversity in placentas from preterm birth.

Background Virome studies on birds, including chickens are relatively scarce, particularly from the African continent. Despite the continuous evolution of RNA viruses and severe losses recorded in poultry from seasonal viral outbreaks, the information on RNA virome composition is even scantier as a result of their highly unstable nature, genetic diversity, and difficulties associated with characterization. Also, information on factors that may modulate the occurrence of some viruses in birds is limited, particularly for domesticated birds. Viral metagenomics through advancements in sequencing technologies, has enabled the characterization of the entire virome of diverse host species using various samples. Methods The complex RNA viral constituents present in 27 faecal samples of asymptomatic chickens from a South African farm collected at 3-time points from two independent seasons were determined, and the impact of the chicken’s age and collection season on viral abundance and diversity was further investigated. The study utilized the non-invasive faecal sampling method, mRNA viral targeted enrichment steps, a whole transcriptome amplification strategy, Illumina sequencing, and bioinformatics tools. Results The results obtained revealed a total of 48 viral species spanning across 11 orders, 15 families and 21 genera. Viral RNA families such as Coronaviridae, Picornaviridae, Reoviridae, Astroviridae, Caliciviridae, Picorbirnaviridae and Retroviridae were abundant, among which picornaviruses, demonstrated a 100% prevalence across the three age groups (2, 4 and 7 weeks) and two seasons (summer and winter) of the 27 faecal samples investigated. A further probe into the extent of variation between the different chicken groups investigated indicated that viral diversity and abundance were significantly influenced by age ( P = 0.01099) and season ( P = 0.00099) between chicken groups, while there was no effect on viral shedding within samples in a group (alpha diversity) for age ( P = 0.146) and season ( P = 0.242). Conclusion The presence of an exceedingly varied chicken RNA virome, encompassing avian, mammalian, fungal, and dietary-associated viruses, underscores the complexities inherent in comprehending the causation, dynamics, and interspecies transmission of RNA viruses within the investigated chicken population. Hence, chickens, even in the absence of discernible symptoms, can harbour viruses that may exhibit opportunistic, commensal, or pathogenic characteristics.

From the perspective of developing relevant interventions for treating HIV and controlling its spread, it is particularly important to comprehensively understand the underlying diversity of the virus, especially in countries where the virus has been present and evolving since the cross-species transmission event that triggered the global pandemic. Here, we generate and phylogenetically analyse sequences derived from the gag-protease (2010 bp; n = 115), partial integrase (345 bp; n = 36), and nef (719 bp; n = 321) genes of HIV-1 group M (HIV-1M) isolates sampled between 2000 and 2022 from two cosmopolitan cities and 40 remote villages of Cameroon. While 52.4% of all sequenced viruses belonged to circulating recombinant form (CRF) 02_AG (CRF02_AG), the remainder were highly diverse, collectively representing seven subtypes and sub-subtypes, eight CRFs, and 36 highly divergent lineages that fall outside the established HIV-1M classification. Additionally, in 77 samples for which at least two genes were typed, 31% of the studied viruses apparently had fragments from viruses belonging to different clades. Furthermore, we found that the distribution of HIV-1M populations is similar between different regions of Cameroon. In contrast, HIV-1M demographics in Cameroon differ significantly from those in its neighbouring countries in the Congo Basin (CB). In phylogenetic trees, viral sequences cluster according to the countries where they were sampled, suggesting that while there are minimal geographical or social barriers to viral dissemination throughout Cameroon, there is strongly impeded dispersal of HIV-1M lineages between Cameroon and other locations of the CB. This suggests that the apparent stability of highly diverse Cameroonian HIV-1M populations may be attributable to the extensive mixing of human populations within the country and the concomitant trans-national movements of major lineages with very similar degrees of fitness; coupled with the relatively infrequent inter-national transmission of these lineages from neighbouring countries in the CB.

Health consequences arising from climate change are threatening to offset advances made to reduce the damage of infectious diseases, which vary by region and the resilience of the local health system. Here we discuss how climate change-related migrations and infectious disease burden are linked through various processes, such as the expansion of pathogens into non-endemic areas, overcrowding in new informal settlements, and the increased proximity of disease vectors and susceptible human populations. Countries that are predicted to have the highest burden are those that have made the least contribution to climate change. Further studies are needed to generate robust evidence on the potential consequences of climate change-related human movements and migration, as well as identify effective and bespoke short- and long-term interventions. Both extreme weather events and long-term gradual changes drive human migration, which could aggravate the burden of infectious diseases. This Perspective examines the complex interplay between climate change, migration and infectious diseases then advocates for context-specific adaptations.

The Chikungunya virus (CHIKV) poses a significant global public health concern, especially in Africa. Since its first isolation in Tanzania in 1953, CHIKV has caused recurrent outbreaks, challenging healthcare systems in low-resource settings. Recent outbreaks in Africa highlight the dynamic nature of CHIKV transmission and the challenges of underreporting and underdiagnosis. Here, we review the literature and analyse publicly available cases, outbreaks, and genomic data, providing insights into the epidemiology, genetic diversity, and transmission dynamics of CHIKV in Africa. Our analyses reveal the circulation of geographically distinct CHIKV genotypes, with certain regions experiencing a disproportionate burden of disease. Phylogenetic analysis of sporadic outbreaks in West Africa suggests repeated emergence of the virus through enzootic spillover, which is markedly different from inferred transmission dynamics in East Africa, where the virus is often introduced from Asian outbreaks, including the recent reintroduction of the Indian Ocean lineage from the Indian subcontinent to East Africa. Furthermore, there is limited evidence of viral movement between these two regions. Understanding the history and transmission dynamics of outbreaks is crucial for effective public health planning. Despite advances in surveillance and research, diagnostic and surveillance challenges persist. This review and secondary analysis highlight the importance of ongoing surveillance, research, and collaboration to mitigate the burden of CHIKV in Africa and improve public health outcomes.

Brazil, Russia, India, China, and South Africa (BRICS) are a group of developing countries with shared economic, healthcare, and scientific interests. These countries navigate multiple syndemics, and the COVID-19 pandemic placed severe strain on already burdened BRICS’ healthcare systems, hampering effective pandemic interventions. Genomic surveillance and molecular epidemiology remain indispensable tools for facilitating informed pandemic intervention. To evaluate the combined manner in which the pandemic unfolded in BRICS countries, we reviewed the BRICS pandemic epidemiological and genomic milestones, which included the first reported cases and deaths, and pharmaceutical and non-pharmaceutical interventions implemented in these countries. To assess the development of genomic surveillance capacity and efficiency over the pandemic, we analyzed the turnaround time from sample collection to data availability and the technologies used for genomic analysis. This data provided information on the laboratory capacities that enable the detection of emerging SARS-CoV-2 variants and highlight their potential for monitoring other pathogens in ongoing public health efforts. Our analyses indicated that BRICS suffered >105.6M COVID-19 infections, resulting in >1.7M deaths. BRICS countries detected intricate genetic combinations of SARS-CoV-2 variants that fueled country-specific pandemic waves. BRICS’ genomic surveillance programs enabled the identification and characterization of the majority of globally circulating Variants of Concern (VOCs) and their descending lineages. Pandemic intervention strategies first implemented by BRICS countries included non-pharmaceutical interventions during the onset of the pandemic, such as nationwide lockdowns, quarantine procedures, the establishment of fever clinics, and mask mandates- which were emulated internationally. Vaccination rollout strategies complemented this, some representing the first of their kind. Improvements in BRICS sequencing and data generation turnaround time facilitated quicker detection of circulating and emerging variants, supported by investments in sequencing and bioinformatic infrastructure. Intra-BRICS cooperation contributed to the ongoing intervention in COVID-19 and other pandemics, enhancing collective capabilities in addressing these health challenges. The data generated continues to inform BRICS-centric pandemic intervention strategies and influences global health matters. The increased laboratory and bioinformatic capacity post-COVID-19 will support the detection of emerging pathogens.

Chronic Hepatitis B Virus (HBV) infection remains a significant public health concern, particularly in Africa, where the burden is substantial. HBV is an enveloped virus, classified into ten phylogenetically distinct genotypes (A-J). Tests to determine HBV genotypes are based on full-genome sequencing or reverse hybridization. In practice, both approaches have limitations. Whereas diagnostic sequencing, generally using the Sanger approach, tends to focus only on the S-gene and yields little or no information on intra-patient HBV genetic diversity, reverse hybridization detects only known genotype-specific mutations. To resolve these limitations, we developed an Oxford Nanopore Technology (ONT)-based HBV diagnostic sequencing protocol suitable for clinical virology that yields both complete genome sequences and extensive intra-patient HBV diversity data. Specifically, the protocol involves tiling-based PCR amplification of HBV sequences, library preparation using the ONT Rapid Barcoding Kit (Oxford nanopore Technologies, Oxford, OX4 4DQ, UK), ONT GridION sequencing, genotyping using genome detective software v1.132/1.133, a recombination analysis using jpHMM (26 October 2011 version) and RDP5.61 software, and drug resistance profiling using Geno2pheno v2.0 software. We prove the utility of our protocol by efficiently generating and characterizing high-quality near full-length HBV genomes from 148 residual diagnostic samples from HBV-infected patients in the Western Cape province of South Africa, providing valuable insights into the genetic diversity and epidemiology of HBV in this region of the world.

Improved RNA virus understanding is critical to studying animal and plant health, and environmental processes. However, the continuous and rapid RNA virus evolution makes their identification and characterization challenging. While recent sequence-based advances have led to extensive RNA virus discovery, there is growing variation in how RNA viruses are identified, analyzed, characterized, and reported. To this end, an RdRp Summit was organized and a hybrid meeting took place in Valencia, Spain in May 2023 to convene leading experts with emphasis on early career researchers (ECRs) across diverse scientific communities. Here we synthesize key insights and recommendations and offer these as a first effort to establish a consensus framework for advancing RNA virus discovery. First, we need interoperability through standardized methodologies, data-sharing protocols, metadata provision and interdisciplinary collaborations and offer specific examples as starting points. Second, as an emergent field, we recognize the need to incorporate cutting-edge technologies and knowledge early and often to improve omic-based viral detection and annotation as novel capabilities reveal new biology. Third, we underscore the significance of ECRs in fostering international partnerships to promote inclusivity and equity in virus discovery efforts. The proposed consensus framework serves as a roadmap for the scientific community to collectively contribute to the tremendous challenge of unveiling the RNA virosphere.

We investigated intra-host genetic evolution using two SARS-CoV-2 isolates from a fully vaccinated (primary schedule x2 doses of AstraZeneca plus a booster of Pfizer), >70-year-old woman with a history of lymphoma and hypertension who presented a SARS-CoV-2 infection for 3 weeks prior to death due to COVID-19. Two full genome sequences were determined from samples taken 13 days apart with both belonging to Pango lineage FL.2: the first detection of this Omicron sub-variant in Botswana. FL.2 is a sub-lineage of XBB.1.9.1. The repertoire of mutations and minority variants in the Spike protein differed between the two time points. Notably, we also observed deletions within the ORF1a and Membrane proteins; both regions are associated with high T-cell epitope density. The internal milieu of immune-suppressed individuals may accelerate SARS-CoV-2 evolution; hence, close monitoring is warranted.

The SARS-CoV-2 genome occupies a unique place in infection biology - it is the most highly sequenced genome on earth (making up over 20% of public sequencing datasets) with fine scale information on sampling date and geography, and has been subject to unprecedented intense analysis. As a result, these phylogenetic data are an incredibly valuable resource for science and public health. However, the vast majority of the data was sequenced by tiling amplicons across the full genome, with amplicon schemes that changed over the pandemic as mutations in the viral genome interacted with primer binding sites. In combination with the disparate set of genome assembly workflows and lack of consistent quality control (QC) processes, the current genomes have many systematic errors that have evolved with the virus and amplicon schemes. These errors have significant impacts on the phylogeny, and therefore over the last few years, many thousands of hours of researchers time has been spent in "eyeballing" trees, looking for artefacts, and then patching the tree. Given the huge value of this dataset, we therefore set out to reprocess the complete set of public raw sequence data in a rigorous amplicon-aware manner, and build a cleaner phylogeny. Here we provide a global tree of 4,471,579 samples, built from a consistently assembled set of high quality consensus sequences from all available public data as of June 2024, viewable at https://viridian.taxonium.org. Each genome was constructed using a novel assembly tool called Viridian (https://github.com/iqbal-lab-org/viridian), developed specifically to process amplicon sequence data, eliminating artefactual errors and mask the genome at low quality positions. We provide simulation and empirical validation of the methodology, and quantify the improvement in the phylogeny. We hope the tree, consensus sequences and Viridian will be a valuable resource for researchers.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants with increased transmissibility, virulence and immune escape abilities have heavily altered the COVID-19 pandemic's course. Deciphering local and global transmission patterns of those variants is thus key in building a profound understanding of the virus' spread around the globe. In the present study, we investigate SARS-CoV-2 variant epidemiology in Côte d'Ivoire, Western sub-Saharan Africa. We therefore generated 234 full SARS-CoV-2 genomes stemming from Central and Northern Côte d'Ivoire. Covering the first and second pandemic wave the country had been facing, we identified 20 viral lineages and showed that in Côte d'Ivoire the second pandemic wave in 2021 was driven by the spread of the Alpha (B.1.1.7) and Eta (B.1.525) variant. Our analyses are consistent with a limited number of international introductions of Alpha and Eta into Côte d'Ivoire, and those introduction events mostly stemmed from within the West African subregion. This suggests that subregional travel to Côte d'Ivoire had more impact on local pandemic waves than direct intercontinental travel.

Ethiopia is the second most populous country in Africa and the sixth most affected by COVID-19 on the continent. Despite having experienced five infection waves, >499,000 cases, and ~7500 COVID-19-related deaths as of January 2023, there is still no detailed genomic epidemiological report on the introduction and spread of SARS-CoV-2 in Ethiopia. In this study, we reconstructed and elucidated the COVID-19 epidemic dynamics. Specifically, we investigated the introduction, local transmission, ongoing evolution, and spread of SARS-CoV-2 during the first four infection waves using 353 high-quality near-whole genomes sampled in Ethiopia. Our results show that whereas viral introductions seeded the first wave, subsequent waves were seeded by local transmission. The B.1.480 lineage emerged in the first wave and notably remained in circulation even after the emergence of the Alpha variant. The B.1.480 was outcompeted by the Delta variant. Notably, Ethiopia's lack of local sequencing capacity was further limited by sporadic, uneven, and insufficient sampling that limited the incorporation of genomic epidemiology in the epidemic public health response in Ethiopia. These results highlight Ethiopia's role in SARS-CoV-2 dissemination and the urgent need for balanced, near-real-time genomic sequencing.

The global pandemic caused by SARS-CoV-2 has increased the demand for scalable sequencing and diagnostic methods, especially for genomic surveillance. Although next-generation sequencing has enabled large-scale genomic surveillance, the ability to sequence SARS-CoV-2 in some settings has been limited by the cost of sequencing kits and the time-consuming preparations of sequencing libraries. We compared the sequencing outcomes, cost and turn-around times obtained using the standard Illumina DNA Prep kit protocol to three modified protocols with fewer clean-up steps and different reagent volumes (full volume, half volume, one-tenth volume). We processed a single run of 47 samples under each protocol and compared the yield and mean sequence coverage. The sequencing success rate and quality for the four different reactions were as follows: the full reaction was 98.2%, the one-tenth reaction was 98.0%, the full rapid reaction was 97.5% and the half-reaction, was 97.1%. As a result, uniformity of sequence quality indicated that libraries were not affected by the change in protocol. The cost of sequencing was reduced approximately seven-fold and the time taken to prepare the library was reduced from 6.5 hours to 3 hours. The sequencing results obtained using the miniaturised volumes showed comparability to the results obtained using full volumes as described by the manufacturer. The adaptation of the protocol represents a lower-cost, streamlined approach for SARS-CoV-2 sequencing, which can be used to produce genomic data quickly and more affordably, especially in resource-constrained settings.

The Alpha, Beta, and Gamma SARS-CoV-2 variants of concern (VOCs) co-circulated globally during 2020 and 2021, fueling waves of infections. They were displaced by Delta during a third wave worldwide in 2021, which, in turn, was displaced by Omicron in late 2021. In this study, we use phylogenetic and phylogeographic methods to reconstruct the dispersal patterns of VOCs worldwide. We find that source-sink dynamics varied substantially by VOC and identify countries that acted as global and regional hubs of dissemination. We demonstrate the declining role of presumed origin countries of VOCs in their global dispersal, estimating that India contributed <15% of Delta exports and South Africa <1%-2% of Omicron dispersal. We estimate that >80 countries had received introductions of Omicron within 100 days of its emergence, associated with accelerated passenger air travel and higher transmissibility. Our study highlights the rapid dispersal of highly transmissible variants, with implications for genomic surveillance along the hierarchical airline network.