Mutational signature dynamics indicate SARS-CoV-2's evolutionary capacity is driven by host antiviral molecules.

The COVID-19 pandemic has been characterised by sequential variant-specific waves shaped by viral, individual human and population factors. SARS-CoV-2 variants are defined by their unique combinations of mutations and there has been a clear adaptation to more efficient human infection since the emergence of this new human coronavirus in late 2019. Here, we use machine learning models to identify shared signatures, i.e., common underlying mutational processes and link these to the subset of mutations that define the variants of concern (VOCs). First, we examined the global SARS-CoV-2 genomes and associated metadata to determine how viral properties and public health measures have influenced the magnitude of waves, as measured by the number of infection cases, in different geographic locations using regression models. This analysis showed that, as expected, both public health measures and virus properties were associated with the waves of regional SARS-CoV-2 reported infection numbers and this impact varies geographically. We attribute this to intrinsic differences such as vaccine coverage, testing and sequencing capacity and the effectiveness of government stringency. To assess underlying evolutionary change, we used non-negative matrix factorisation and observed three distinct mutational signatures, unique in their substitution patterns and exposures from the SARS-CoV-2 genomes. Signatures 1, 2 and 3 were biased to C→T, T→C/A→G and G→T point mutations. We hypothesise assignments of these mutational signatures to the host antiviral molecules APOBEC, ADAR and ROS respectively. We observe a shift amidst the pandemic in relative mutational signature activity from predominantly Signature 1 changes to an increasingly high proportion of changes consistent with Signature 2. This could represent changes in how the virus and the host immune response interact and indicates how SARS-CoV-2 may continue to generate variation in the future. Linkage of the detected mutational signatures to the VOC-defining amino acids substitutions indicates the majority of SARS-CoV-2's evolutionary capacity is likely to be associated with the action of host antiviral molecules rather than virus replication errors.

Genomic epidemiology of the rotavirus G2P[4] strains in coastal Kenya pre- and post-rotavirus vaccine introduction, 2012-8.

The introduction of rotavirus vaccines into the national immunization programme in many countries has led to a decline in childhood diarrhoea disease burden. Coincidentally, the incidence of some rotavirus group A (RVA) genotypes has increased, which may result from non-vaccine-type replacement. Here, we investigate the evolutionary genomics of rotavirus G2P[4] which has shown an increase in countries that introduced the monovalent Rotarix® vaccine. We examined sixty-three RVA G2P[4] strains sampled from children (aged below 13 years) admitted to Kilifi County Hospital, coastal Kenya, pre- (2012 to June 2014) and post-(July 2014 to 2018) rotavirus vaccine introduction. All the sixty-three genome sequences showed a typical DS-1-like genome constellation (G2-P[4]-I2-R2-C2-M2-A2-N2-T2-E2-H2). Pre-vaccine G2 sequences predominantly classified as sub-lineage IVa-3 and co-circulated with low numbers of sub-lineage IVa-1 strains, whereas post-vaccine G2 sequences mainly classified into sub-lineage IVa-3. In addition, in the pre-vaccine period, P[4] sub-lineage IVa strains co-circulated with low numbers of P[4] lineage II strains, but P[4] sub-lineage IVa strains predominated in the post-vaccine period. On the global phylogeny, the Kenyan pre- and post-vaccine G2P[4] strains clustered separately, suggesting that different virus populations circulated in the two periods. However, the strains from both periods exhibited conserved amino acid changes in the known antigenic epitopes, suggesting that replacement of the predominant G2P[4] cluster was unlikely a result of immune escape. Our findings demonstrate that the pre- and post-vaccine G2P[4] strains circulating in Kilifi, coastal Kenya, differed genetically but likely were antigenically similar. This information informs the discussion on the consequences of rotavirus vaccination on rotavirus diversity.

Rapid, early, and potent Spike-directed IgG, IgM, and IgA distinguish asymptomatic from mildly symptomatic COVID-19 in Uganda, with IgG persisting for 28 months.

Introduction: Understanding how spike (S)-, nucleoprotein (N)-, and RBD-directed antibody responses evolved in mild and asymptomatic COVID-19 in Africa and their interactions with SARS-CoV-2 might inform development of targeted treatments and vaccines. Methods: Here, we used a validated indirect in-house ELISA to characterise development and persistence of S- and N-directed IgG, IgM, and IgA antibody responses for 2430 SARS-CoV-2 rt-PCR-diagnosed Ugandan specimens from 320 mild and asymptomatic COVID-19 cases, 50 uninfected contacts, and 54 uninfected non-contacts collected weekly for one month, then monthly for 28 months. Results: During acute infection, asymptomatic patients mounted a faster and more robust spike-directed IgG, IgM, and IgA response than those with mild symptoms (Wilcoxon rank test, p-values 0.046, 0.053, and 0.057); this was more pronounced in males than females. Spike IgG antibodies peaked between 25 and 37 days (86.46; IQR 29.47-242.56 BAU/ml), were significantly higher and more durable than N- and RBD IgG antibodies and lasted for 28 months. Anti-spike seroconversion rates consistently exceeded RBD and nucleoprotein rates. Spike- and RBD-directed IgG antibodies were positively correlated until 14 months (Spearman's rank correlation test, p-values 0.0001 to 0.05), although RBD diminished faster. Significant anti-spike immunity persisted without RBD. 64% and 59% of PCR-negative, non-infected non-contacts and suspects, exhibited baseline SARS-CoV-2 N-IgM serological cross-reactivity, suggesting undetected exposure or abortive infection. N-IgG levels waned after 787 days, while N-IgM levels remained undetectable throughout. Discussion: Lower N-IgG seroconversion rates and the absence of N-IgM indicate that these markers substantially underestimate the prior exposure rates. Our findings provide insights into the development of S-directed antibody responses in mild and asymptomatic infections, with varying degrees of symptoms eliciting distinct immune responses, suggesting distinct pathogenic pathways. These longer-lasting data inform vaccine design, boosting strategies, and surveillance efforts in this and comparable settings.

Viral epidemic preparedness: a perspective from five clinical microbiology laboratories in Europe.

BACKGROUND: Pandemic preparedness is critical to respond effectively to existing and emerging/new viral pathogens. Important lessons have been learned during the last pandemic at various levels. This revision discusses some of the major challenges and potential ways to address them in the likely event of future pandemics. OBJECTIVES: To identify critical points of readiness that may help us accelerate the response to future pandemics from a clinical microbiology laboratory perspective with a focus on viral diagnostics and genomic sequencing. The potential areas of improvement identified are discussed from the sample collection to information reporting. SOURCES: Microbiologists and researchers from five countries reflect on challenges encountered during the COVID-19 pandemic, review published literature on prior and current pandemics, and suggest potential solutions in preparation for future outbreaks. CONTENT: Major challenges identified in the pre-analytic and post-analytic phases from sample collection to result reporting are discussed. From the perspective of clinical microbiology laboratories, the preparedness for a new pandemic should focus on zoonotic viruses. Laboratory readiness for scalability is critical and should include elements related to material procurement, training personnel, specific funding programmes, and regulatory issues to rapidly implement "in-house" tests. Laboratories across various countries should establish (or re-use) operational networks to communicate to respond effectively, ensuring the presence of agile circuits with full traceability of samples. IMPLICATIONS: Laboratory preparedness is paramount to respond effectively to emerging and re-emerging viral infections and to limit the clinical and societal impact of new potential pandemics. Agile and fully traceable methods for sample collection to report are the cornerstone of a successful response. Expert group communication and early involvement of information technology personnel are critical for preparedness. A specific budget for pandemic preparedness should be ring-fenced and added to the national health budgets.

Optimisation and Validation of a conventional ELISA and cut-offs for detecting and quantifying anti-SARS-CoV-2 Spike, RBD, and Nucleoprotein IgG, IgM, and IgA antibodies in Uganda.

There is an urgent need for better immunoassays to measure antibody responses as part of immune-surveillance activities and to profile immunological responses to emerging SARS-CoV-2 variants. We optimised and validated an in-house conventional ELISA to identify and quantify SARS-CoV-2 spike- (S-), receptor binding domain- (RBD-), and nucleoprotein- (N-) directed IgG, IgM, and IgA binding antibodies in the Ugandan population and similar settings. Pre- and post-pandemic specimens were used to compare the utility of mean ± 2SD, mean ± 3SD, 4-fold above blanks, bootstrapping, and receiver operating characteristic (ROC) analyses in determining optimal cut-off optical densities at 450 nm (OD) for discriminating between antibody positives and negatives. "Limits of detection" (LOD) and "limits of quantitation" (LOQ) were validated alongside the assay's uniformity, accuracy, inter-assay and inter-operator precision, and parallelism. With spike-directed sensitivity and specificity of 95.33 and 94.15%, respectively, and nucleoprotein sensitivity and specificity of 82.69 and 79.71%, ROC was chosen as the best method for determining cutoffs. Accuracy measurements were within the expected CV range of 25%. Serum and plasma OD values were highly correlated (r = 0.93, p=0.0001). ROC-derived cut-offs for S-, RBD-, and N-directed IgG, IgM, and IgA were 0.432, 0.356, 0.201 (S), 0.214, 0.350, 0.303 (RBD), and 0.395, 0.229, 0.188 (N). The sensitivity and specificity of the S-IgG cut-off were equivalent to the WHO 20/B770-02 S-IgG reference standard at 100% level. Spike negative IgG, IgM, and IgA ODs corresponded to median antibody concentrations of 1.49, 3.16, and 0 BAU/mL, respectively, consistent with WHO low titre estimates. Anti-spike IgG, IgM, and IgA cut-offs were equivalent to 18.94, 20.06, and 55.08 BAU/mL. For the first time, we provide validated parameters and cut-off criteria for the in-house detection of subclinical SARS-CoV-2 infection and vaccine-elicited binding antibodies in the context of Sub-Saharan Africa and populations with comparable risk factors.

Evolution of increased positive charge on the SARS-CoV-2 spike protein may be adaptation to human transmission.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to evolve and infect individuals. The exterior surface of the SARS-CoV-2 virion is dominated by the spike protein, and the current work examined spike protein biochemical features that have changed during the 3 years in which SARS-CoV-2 has infected humans. Our analysis identified a striking change in spike protein charge, from -8.3 in the original Lineage A and B viruses to -1.26 in most of the current Omicron viruses. We conclude that in addition to immune selection pressure, the evolution of SARS-CoV-2 has also altered viral spike protein biochemical properties, which may influence virion survival and promote transmission. Future vaccine and therapeutic development should also exploit and target these biochemical properties.

SARS-CoV-2 Omicron BA.5 Infections in Vaccinated Persons, Rural Uganda.

We describe a cluster of COVID-19 breakthrough infections after vaccination in Kyamulibwa, Kalungu District, Uganda. All but 1 infection were from SARS-CoV-2 Omicron strain BA.5.2.1. We identified 6 distinct genotypes by genome sequencing. Infections were mild, suggesting vaccination is not protective against infection but may limit disease severity.

The evolving SARS-CoV-2 epidemic in Africa: Insights from rapidly expanding genomic surveillance.

Investment in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequencing in Africa over the past year has led to a major increase in the number of sequences that have been generated and used to track the pandemic on the continent, a number that now exceeds 100,000 genomes. Our results show an increase in the number of African countries that are able to sequence domestically and highlight that local sequencing enables faster turnaround times and more-regular routine surveillance. Despite limitations of low testing proportions, findings from this genomic surveillance study underscore the heterogeneous nature of the pandemic and illuminate the distinct dispersal dynamics of variants of concern-particularly Alpha, Beta, Delta, and Omicron-on the continent. Sustained investment for diagnostics and genomic surveillance in Africa is needed as the virus continues to evolve while the continent faces many emerging and reemerging infectious disease threats. These investments are crucial for pandemic preparedness and response and will serve the health of the continent well into the 21st century.

Establishing farm dust as a useful viral metagenomic surveillance matrix.

Farm animals may harbor viral pathogens, some with zoonotic potential which can possibly cause severe clinical outcomes in animals and humans. Documenting the viral content of dust may provide information on the potential sources and movement of viruses. Here, we describe a dust sequencing strategy that provides detailed viral sequence characterization from farm dust samples and use this method to document the virus communities from chicken farm dust samples and paired feces collected from the same broiler farms in the Netherlands. From the sequencing data, Parvoviridae and Picornaviridae were the most frequently found virus families, detected in 85-100% of all fecal and dust samples with a large genomic diversity identified from the Picornaviridae. Sequences from the Caliciviridae and Astroviridae familes were also obtained. This study provides a unique characterization of virus communities in farmed chickens and paired farm dust samples and our sequencing methodology enabled the recovery of viral genome sequences from farm dust, providing important tracking details for virus movement between livestock animals and their farm environment. This study serves as a proof of concept supporting dust sampling to be used in viral metagenomic surveillance.

Near-Complete Genome Sequences of Measles Virus Strains from 10 Years of Uganda Country-wide Surveillance.

Measles remains a global health challenge despite the availability of a safe and effective vaccine. Sporadic outbreaks of measles virus infections continue in Uganda. We report eight near-complete genome sequences of measles virus strains from Uganda cases from 2011 to 2020, providing useful data for assessing vaccine escape and local/regional transmission.

Transmission networks of SARS-CoV-2 in Coastal Kenya during the first two waves: A retrospective genomic study.

Background: Detailed understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) regional transmission networks within sub-Saharan Africa is key for guiding local public health interventions against the pandemic. Methods: Here, we analysed 1139 SARS-CoV-2 genomes from positive samples collected between March 2020 and February 2021 across six counties of Coastal Kenya (Mombasa, Kilifi, Taita Taveta, Kwale, Tana River, and Lamu) to infer virus introductions and local transmission patterns during the first two waves of infections. Virus importations were inferred using ancestral state reconstruction, and virus dispersal between counties was estimated using discrete phylogeographic analysis. Results: During Wave 1, 23 distinct Pango lineages were detected across the six counties, while during Wave 2, 29 lineages were detected; 9 of which occurred in both waves and 4 seemed to be Kenya specific (B.1.530, B.1.549, B.1.596.1, and N.8). Most of the sequenced infections belonged to lineage B.1 (n = 723, 63%), which predominated in both Wave 1 (73%, followed by lineages N.8 [6%] and B.1.1 [6%]) and Wave 2 (56%, followed by lineages B.1.549 [21%] and B.1.530 [5%]). Over the study period, we estimated 280 SARS-CoV-2 virus importations into Coastal Kenya. Mombasa City, a vital tourist and commercial centre for the region, was a major route for virus imports, most of which occurred during Wave 1, when many Coronavirus Disease 2019 (COVID-19) government restrictions were still in force. In Wave 2, inter-county transmission predominated, resulting in the emergence of local transmission chains and diversity. Conclusions: Our analysis supports moving COVID-19 control strategies in the region from a focus on international travel to strategies that will reduce local transmission. Funding: This work was funded by The Wellcome (grant numbers: 220985, 203077/Z/16/Z, 220977/Z/20/Z, and 222574/Z/21/Z) and the National Institute for Health and Care Research (NIHR), project references: 17/63/and 16/136/33 using UK Aid from the UK government to support global health research, The UK Foreign, Commonwealth and Development Office. The views expressed in this publication are those of the author(s) and not necessarily those of the funding agencies.

Spike Protein Cleavage-Activation in the Context of the SARS-CoV-2 P681R Mutation: an Analysis from Its First Appearance in Lineage A.23.1 Identified in Uganda.

Based on its predicted ability to affect transmissibility and pathogenesis, surveillance studies have highlighted the role of a specific mutation (P681R) in the S1/S2 furin cleavage site of the SARS-CoV-2 spike protein. Here we analyzed A.23.1, first identified in Uganda, as a P681R-containing virus several months prior to the emergence of B.1.617.2 (Delta variant). We performed assays using peptides mimicking the S1/S2 from A.23.1 and B.1.617 and observed significantly increased cleavability with furin compared to both an original B lineage (Wuhan-Hu1) and B.1.1.7 (Alpha variant). We also performed cell-cell fusion and functional infectivity assays using pseudotyped particles and observed an increase in activity for A.23.1 compared to an original B lineage spike. However, these changes in activity were not reproduced in the B lineage spike bearing only the P681R substitution. Our findings suggest that while A.23.1 has increased furin-mediated cleavage linked to the P681R substitution, this substitution needs to occur on the background of other spike protein changes to enable its functional consequences. IMPORTANCE During the course of the SARS-CoV-2 pandemic, viral variants have emerged that often contain notable mutations in the spike gene. Mutations that encode changes in the spike S1/S2 (furin) activation site have been considered especially impactful. The S1/S2 change from proline to arginine at position 681 (P681R) first emerged in the A.23.1 variant in Uganda, and subsequently occurred in the more widely transmitted Delta variant. We show that the A.23.1 spike is more readily activated by the host cell protease furin, but that this is not reproduced in an original SARS-CoV-2 spike containing the P681R mutation. Changes to the S1/S2 (furin) activation site play a role in SARS-CoV-2 infection and spread, but successful viruses combine these mutations with other less well identified changes, occurring as part of natural selection.

Field evaluation of the performance of seven Antigen Rapid diagnostic tests for the diagnosis of SARs-CoV-2 virus infection in Uganda.

OBJECTIVE: The objective of this study was to evaluate the performance of seven antigen rapid diagnostic tests (Ag RDTs) in a clinical setting to identify those that could be recommended for use in the diagnosis of SARS-CoV-2 infection in Uganda. METHODS: This was a cross-sectional prospective study. Nasopharyngeal swabs were collected consecutively from COVID-19 PCR positive and COVID-19 PCR negative participants at isolation centers and points of entry, and tested with the SARS-CoV-2 Ag RDTs. Test sensitivity and specificity were generated by comparing results against qRT-PCR results (Berlin Protocol) at a cycle threshold (Ct) cut-off of ≤39. Sensitivity was also calculated at Ct cut-offs ≤29 and ≤33. RESULTS: None of the Ag RDTs had a sensitivity of ≥80% at Ct cut-off values ≤33 and ≤39. Two kits, Panbio™ COVID-19 Ag and VivaDiag™ SARS-CoV-2 Ag had a sensitivity of ≥80% at a Ct cut-off value of ≤29. Four kits: BIOCREDIT COVID -19 Ag, COVID-19 Ag Respi-Strip, MEDsan® SARS-CoV-2 Antigen Rapid Test and Panbio™ COVID-19 Ag Rapid Test had a specificity of ≥97%. CONCLUSIONS: This evaluation identified one Ag RDT, Panbio™ COVID-19 Ag with a performance at high viral load (Ct value ≤29) reaching that recommended by WHO. This kit was recommended for screening of patients with COVID -19-like symptoms presenting at health facilities.

Identification of missed viruses by metagenomic sequencing of clinical respiratory samples from Kenya.

Pneumonia remains a major cause of mortality and morbidity. Most molecular diagnoses of viruses rely on polymerase chain reaction (PCR) assays that however can fail due to primer mismatch. We investigated the performance of routine virus diagnostics in Kilifi, Kenya, using random-primed viral next generation sequencing (viral NGS) on respiratory samples which tested negative for the common viral respiratory pathogens by a local standard diagnostic panel. Among 95 hospitalised pneumonia patients and 95 household-cohort individuals, analysis of viral NGS identified at least one respiratory-associated virus in 35 (37%) and 23 (24%) samples, respectively. The majority (66%; 42/64) belonged to the Picornaviridae family. The NGS data analysis identified a number of viruses that were missed by the diagnostic panel (rhinovirus, human metapneumovirus, respiratory syncytial virus and parainfluenza virus), and these failures could be attributed to PCR primer/probe binding site mismatches. Unexpected viruses identified included parvovirus B19, enterovirus D68, coxsackievirus A16 and A24 and rubella virus. The regular application of such viral NGS could help evaluate assay performance, identify molecular causes of missed diagnoses and reveal gaps in the respiratory virus set used for local screening assays. The results can provide actionable information to improve the local pneumonia diagnostics and reveal locally important viral pathogens.

Spike protein cleavage-activation mediated by the SARS-CoV-2 P681R mutation: a case-study from its first appearance in variant of interest (VOI) A.23.1 identified in Uganda.

The African continent like all other parts of the world with high infection/low vaccination rates can, and will, be a source of novel SARS-CoV-2 variants. The A.23 viral lineage, characterized by three spike mutations F157L, V367F and Q613H, was first identified in COVID-19 cases from a Ugandan prison in July 2020, and then was identified in the general population with additional spike mutations (R102I, L141F, E484K and P681R) to comprise lineage A.23.1 by September 2020, with this virus being designated a variant of interest (VOI) in Africa and with subsequent spread to 26 other countries. The P681R spike substitution of the A.23.1 VOI is of note as it increases the number of basic residues in the sub-optimal SARS-CoV-2 spike protein furin cleavage site; as such, this substitution may affect viral replication, transmissibility or pathogenic properties. The same P681R substitution has also appeared in B.1.617 variants, including B.1.617.2 (Delta). Here, we performed assays using fluorogenic peptides mimicking the S1/S2 sequence from A.23.1 and B.1.617.2 and observed significantly increased cleavability with furin, compared to sequences derived from the original Wuhan-Hu1 S1/S2. We performed functional infectivity assays using pseudotyped MLV particles harboring SARS-CoV-2 spike proteins and observed an increase in transduction for A.23.1-pseudotyped particles compared to Wuhan-Hu-1 in Vero-TMPRSS2 and Calu-3 cells (with a presumed early entry pathway), although lowered infection in Vero E6 cells (with a presumed late entry pathway). However, these changes in infectivity were not reproduced in the original Wuhan-Hu-1 spike bearing only the P681R substitution. Our findings suggest that while A.23.1 has increased furin-mediated cleavage linked to the P681R substitution, which may affect viral infection and transmissibility, this substitution alone is not sufficient and needs to occur on the background of other spike protein changes to enable its full functional consequences.

SARS-CoV-2 Variants, South Sudan, January-March 2021.

As the coronavirus pandemic continues, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequence data are required to inform vaccine efforts. We provide SARS-CoV-2 sequence data from South Sudan and document the dominance of SARS-CoV-2 lineage B.1.525 (Eta variant) during the country's second wave of infection.

Evaluation of the performance of 25 SARS-CoV-2 serological rapid diagnostic tests using a reference panel of plasma specimens at the Uganda Virus Research Institute.

INTRODUCTION: Serological testing is needed to better understand the epidemiology of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus. Rapid diagnostic tests (RDTs) have been developed to detect specific antibodies, IgM and IgG, to the virus. The performance of 25 of these RDTs was evaluated. METHODS: A serological reference panel of 50 positive and 100 negative plasma specimens was developed from SARS-CoV-2 PCR and antibody positive patients and pre-pandemic SARS-CoV-2-negative specimens collected in 2016. Test performance of the 25 RDTs was evaluated against this panel. RESULTS: A total of 10 RDTs had a sensitivity ≥98%, while 13 RDTs had a specificity ≥98% to anti-SARS-CoV-2 IgG antibodies. Four RDTs (Boson, MultiG, Standard Q, and VivaDiag) had both sensitivity and specificity ≥98% to anti-SARS-CoV-2 IgG antibodies. Only three RDTs had a sensitivity ≥98%, while 10 RDTs had a specificity ≥98% to anti-SARS-CoV-2 IgM antibodies. Three RDTs (Autobio, MultiG, and Standard Q) had sensitivity and specificity ≥98% to combined IgG/IgM. The RDTs that performed well also had perfect or almost perfect inter-reader agreement. CONCLUSIONS: This evaluation identified three RDTs with a sensitivity and specificity to IgM/IgG antibodies of ≥98% with the potential for widespread antibody testing in Uganda.

Unique protein features of SARS-CoV-2 relative to other Sarbecoviruses.

Defining the unique properties of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) protein sequences has potential to explain the range of Coronavirus Disease 2019 severity. To achieve this we compared proteins encoded by all Sarbecoviruses using profile Hidden Markov Model similarities to identify protein features unique to SARS-CoV-2. Consistent with previous reports, a small set of bat- and pangolin-derived Sarbecoviruses show the greatest similarity to SARS-CoV-2 but are unlikely to be the direct source of SARS-CoV-2. Three proteins (nsp3, spike, and orf9) showed regions differing between the bat Sarbecoviruses and SARS-CoV-2 and indicate virus protein features that might have evolved to support human infection and/or transmission. Spike analysis identified all regions of the protein that have tolerated change and revealed that the current SARS-CoV-2 variants of concern have sampled only a fraction (∼31 per cent) of the possible spike domain changes which have occurred historically in Sarbecovirus evolution. This result emphasises the evolvability of these coronaviruses and the potential for further change in virus replication and transmission properties over the coming years.