Genomic epidemiology of SARS-CoV-2 in Mauritius reveals a new wave of infections dominated by the B.1.1.318, a variant under investigation ========================================================================================================================================== * Houriiyah Tegally * Magalutcheemee Ramuth * Daniel Amoaka * Cathrine Scheepers * Eduan Wilkinson * Marta Giovanetti * Richard J Lessells * Jennifer Giandhari * Arshad Ismail * Darren Martin * Emmanuel James San * Margaret Crawford * Rodney S Daniels * Ruth Harvey * Somduthsingh Bahadoor * Janaki Sonoo * Myriam Timol * Lovena Veerapa-Mangroo * Anne von Gottberg * Jinal N. Bhiman * Tulio de Oliveira * Shyam Manraj ## Abstract Mauritius, a small island in the Indian Ocean, has had a unique experience of the SARS-CoV-2 pandemic. In March 2020, Mauritius endured a small first wave and quickly implemented control measures which allowed elimination of local transmission of SARS-CoV-2. When borders to the island reopened, it was accompanied by mandatory quarantine and testing of incoming passengers to avoid reintroduction of the virus into the community. As variants of concern (VOCs) emerged elsewhere in the world, Mauritius began using genomic surveillance to keep track of quarantined cases of these variants. In March 2021, another local outbreak occurred, and sequencing was used to investigate this new wave of local infections. Here, we analyze 154 SARS-CoV-2 viral genomes from Mauritius, which represent 12% of all the infections seem in Mauritius, these were both from specimens of incoming passengers before March 2021 and those of cases during the second wave. Our findings indicate that despite the presence of known VOCs Beta (B.1.351) and Alpha (B.1.1.7) among quarantined passengers, the second wave of local SARS-CoV-2 infections in Mauritius was caused by a single introduction and dominant circulation of the B.1.1.318 virus. The B.1.1.318 variant is characterized by fourteen non-synonymous mutations in the S-gene, with five encoded amino acid substitutions (T95I, E484K, D614G, P681H, D796H) and one deletion (Y144del) in the Spike glycoprotein. This variant seems to be increasing in prevalence and it is now present in 34 countries. This study highlights that despite having stopped the introduction of more transmissible VOCs by travel quarantines, a single undetected introduction of a B.1.1.318 lineage virus was enough to initiate a large local outbreak in Mauritius and demonstrated the need for continuous genomic surveillance to fully inform public health decisions. ## Introduction Genomic surveillance has been an important tool for characterizing the COVID-19 pandemic, through mapping dissemination dynamics of the virus across the world or tracking the local movement and population dynamics of the virus within countries, cities and communities (1– 6). Consequently, by the end of May 2021, the GISAID database contained more than 1.6 million SARS-CoV-2 near-full genome sequences (7). In November and December 2020, after months of ongoing transmission in many areas of the world, genomic surveillance in the UK, South Africa, and Brazil detected three SARS-CoV-2 variants of concern (VOCs): Alpha (B.1.1.7) (8), Beta (501Y.V2/B.1.351) (9) and Gamma (P.1) (10). These VOCs each have multiple mutations, predominantly in the S-gene, that have subsequently been found to be associated with greater resistance to immunity induced by infection/vaccination, enhanced transmissibility, and/or increased pathogenicity (11–14). As these VOCs began spreading throughout the world via air travel (15), the impact of travel restrictions and traveler quarantines were again debated. A few countries which had been successful at near-elimination of COVID-19, turned to the sequencing of variants within SARS-CoV-2 positive travelers to inform public health measures aimed at keeping VOCs at bay. Sydney, Australia, for example, first detected an outbreak of the Beta variant in a travel quarantine facility and quickly followed through with a city-wide lockdown for a few days (16). Mauritius, a small island in the Indian Ocean with 1.2 million inhabitants, has had a unique experience of, and response to, the COVID-19 pandemic. Since early 2020, before the identification of the first cases in the country, Mauritius began requiring 14-day quarantine for incoming passengers arriving from the most affected countries of Asia and Europe (17). On 18 March 2020, the country recorded its first three positive cases in travelers who had not been required to quarantine (not arriving from countries defined to be most at risk), and the following day the country went into a strict national lockdown, consisting of the closure of all non-essential economic activities, educational institutions and a ban of entry into the country. The first wave of infections lasted for approximately one month with a total of around 340 local cases across the island. With zero ongoing community cases following this, the country lifted all local restrictions on 15 June 2020 but borders remained closed. Total travel bans were eased on 1 October 2020 but a strict government managed 14-day hotel quarantine and regular RT-PCR tests for any incoming passengers were imposed by the government to prevent new infections from entering the community. As VOC emerged, travel bans were again imposed on passengers coming directly from the UK and South Africa. Between 30 April 2020 to February 2021, 23,471 passengers arrived in Mauritius through this quarantine system. Of these, 281 tested positive for COVID-19 during the said quarantine period. From June 2020 to February 2021, no local COVID-19 cases were recorded in the country. This unfortunately changed in March 2021 when a cluster of local cases was detected followed by a second wave of more than 650 local cases. So far, this second wave has been controlled with rigorous testing of symptomatic individuals, their contacts a significant number of individuals within targeted areas as well as by the isolation of cases, contact tracing and testing the contacts of those in quarantine. Here we describe the use of SARS-CoV-2 sequencing for the genomic surveillance of both community and travel-associated COVID-19 cases in Mauritius to help investigate this second wave and to inform the public health response to avert future outbreaks. We demonstrate how coupling quarantining of incoming travelers with SARS-CoV-2 testing has helped to prevent multiple introductions of two known VOCs into the country. Genomic surveillance also revealed that the second wave was dominated by a single introduction to the country of a potential variant of interest (VOI), B.1.1.318. This VOI has fourteen amino acid mutations in the Spike protein, including the E484K and P681H, which are also found in other VOCs. ## Results ### Travel quarantine control of VOCs From October 2020 to February 2021, Mauritius re-opened its borders to commercial flights and allowed the entry of both citizens and tourists under a strict testing and (managed) quarantine protocol. All visitors to Mauritius were required to provide confirmation of a negative PCR test 48 hours prior to flight-boarding and had to quarantine for 14 days in a hotel, undergoing three more RT-PCR tests at days 0, 7 and 14 of the quarantine period, all at their own cost. During this time, the country recorded a low and stable number of positive cases in quarantine facilities (Fig 1A), from which infected individuals were immediately transferred to treatment centers and isolated until testing negative again. Sequencing of residual diagnostic samples from these individuals revealed multiple cases of infection with VOCs. There were six independent introductions (including two paired travelers) of the Beta (501Y.V2/B.1.351) variant from Africa and Europe, and two independent introductions of the Alpha (B.1.1.7) from Europe (Fig 1B, Supp Fig S1). Sequencing of specimens from these travelers also detected one case of the A.23.1 variant, three cases of the European lineage B.1.177.4 and an infection incident of one health worker in a quarantine facility by two travel-related cases carrying the B.1.36 lineage (Supp Figure S1). The B.1.36 lineage was common in India, where these travelers had come from at the beginning of 2021. Subsequent sequencing of community cases during the second wave of local infections in March 2021 (124 of 432 positive cases sequenced in March 2021 suggested that the VOCs detected in quarantine did not spill over into the community (Fig 1B). ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/16/2021.06.16.21259017/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2021/06/16/2021.06.16.21259017/F1) Figure 1. B.1.1.318 variants dominating second wave of infections in Mauritius. A) Epidemiological curve showing daily number of COVID-19 cases (blue) and associated deaths (orange) in Mauritius from October 2020 to May 2021, alongside a timeline of the status of lockdown and border control status. The dates of genome sampling are also shown, distinguishing genomes from incoming passengers (n=34) (pink) or local cases (n=120) (green). B) ML tree of 4484 genomes, including 154 genomes from Mauritius, colored to show VOCs Alpha (501Y.V2/B.1.351) and Beta (B.1.1.7) and the VOI B.1.1.318. C) Protein map of amino acid mutations in the B.1.1.318 VOI, with a zoomed focus on the Spike glycoprotein. ### Second wave dominated by B.1.1.318, potential VOI On 5 March 2021, Mauritius detected a first community case of COVID-19 in what would become the country’s second wave of infections. Epidemiologically, the second wave was characterized by multiple clusters in family, workplace and health facility settings and the country went into its second lockdown on 9 March 2021. At the peak of the second wave, Mauritius recorded 68 cases in one day on 11 March 2021 with 2555 tests carried out that day (positivity rate of 2.6%) following an average of approximately 1500 tests per day during this second wave (Source: Government Information Services Mauritius). Unfortunately, the country also recorded seven new COVID-related deaths during this time. To understand the origins and dynamics of this new wave of infections and to assess the potential impact of the lineages on the planned vaccine roll-out, virus-positive samples from 120 cases from the March 2021 outbreak were sequenced. The geographical sampling of these sequenced cases closely matched the concentration of cases during this outbreak within the districts of Plaines Wilhems and Black River (Supp Fig S2 – genomes sampling map and cases map). These were in addition to the successful sequencing of 34 samples from quarantined incoming passengers, placing Mauritius in fifth position worldwide for the proportion of cases (n=1287) sequenced (12%) at this timepoint (Supp Table S1). A maximum likelihood (ML) tree representing the evolutionary relationships of the 154 Mauritian genomes with a representative subset of 4330 genomes sampled globally indicated that all 120 of the genomes sampled from local community-derived infections fall within a cluster of sequences belonging to variant lineage B.1.1.318 (Fig 1B). This lineage is characterized by fourteen non-synonymous mutations in the S-gene, with five encoded amino acid substitutions (T95I, E484K, D614G, P681H, D796H) and one deletion (Y144del) in the Spike glycoprotein (Fig 1C). While this lineage has not been defined as a VOI by the World Health Organization (WHO), several of the Spike glycoprotein mutations are common to other VOCs and VOIs (Supp Fig S3), including mutations associated with immune escape (E484K, Y144del and likely T95I) and increased transmissibility (P681H). Moreover, all five of the spike codon sites where lineage-defining mutations have occurred are detectably evolving under positive selection (p <0.05 with both the MEME and IFEL positive selection detection methods) (Supp Table S2). Overall, therefore, phylogenetic analysis of this outbreak cluster in Mauritius suggests that it was initiated by a single undetected introduction of a B.1.1.318 variant virus. ### International epidemiology of the B.1.1.318 variant lineage To contextualize the B.1.1.318 variant outbreak in Mauritius, we downloaded all sequences belonging to the B.1.1.318 lineage on GISAID as of 17 May 2021. We identified 34 countries in total where cases of infection with B.1.1.318 variants had been sequenced (Fig 2A). Although Mauritius is the only country where B.1.1.318 has accounted for the majority of infections (76% of all sequenced cases), this variant has been classified as a variant under investigation (VUI-202102/04) by Public Health England since February 2021 due to it carrying the Spike glycoprotein E484K substitution (18). While this variant has been consistently sampled since the beginning of 2021 in several countries, such as the UK, the US, Germany and Ireland, with first detection being in Nigeria (West Africa), it has not been shown to circulate elsewhere at such high proportions as was seen in Mauritius (Fig 2A). Other countries with the highest prevalence of B.1.1.318 variants are Gabon (accounting for 19% of all sequenced cases), Togo (7%), and Greece (3%) (Fig 2A). ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/16/2021.06.16.21259017/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2021/06/16/2021.06.16.21259017/F2) Figure 2. Global distribution and spread of the B.1.1.318 lineage. A) Prevalence of the B.1.1.318 lineage worldwide shown over time and as a proportion of all sequenced cases in each country. B) MCC tree of the B.1.1.318 lineage worldwide discretized by continental regions. C) International phylogeography of the B.1.1.318 lineage (lines showing spread in anti-clockwise direction). A time-resolved Maximum Clade Credibility (MCC) analysis was performed with all B.1.1.318 variant sequences in the EpiCov™ database of GISAID. The resulting phylogeny indicates a potential origin of this lineage in Nigeria in late December 2020 (HPD: late December 2020 – early January 2021) (Fig 2B). The results show evidence of multiple introductions to Europe in January 2021 and onward transmission within Europe (Fig 2B). The lineage was introduced to North America both from Africa in late January 2021 and from Europe multiple times between mid-January 2021 and March 2021, with evidence of onward transmission in North America. Introduction of the lineage to Asia happened in late January 2021 with no evidence of onward transmission (Fig 2B, 2C). The introduction of the B.1.1.318 to Mauritius seems to have occurred via a single source around mid-February 2021 HPD: mid-February 2021 – late February 2021) (Fig 2B) and spread as small onward clusters within the country. The virus seems to have been imported to Mauritius from Europe with the closest relatives in the phylogeny being sequences of viruses from Ireland, Finland, and Norway (Fig 2C, Supp Fig S4). ## Discussion Being a small island state, Mauritius was able to respond rapidly and aggressively to the first wave of SARS-CoV-2 infections in the country in March 2020. By closing borders and instituting an immediate strict lockdown, Mauritius was able to achieve zero transmission of the virus in the community and with ongoing public health and social measures, life resumed normally a couple of months later. In an attempt to revive tourism, borders reopened to both Mauritian citizens and foreigners in October 2020, but a strict testing and managed quarantine protocol was established to prevent any re-introduction of SARS-CoV-2 into the local community. As in Australia and New Zealand, this system proved to be effective in preventing new virus introductions and local outbreaks. In fact, in this study we show that Mauritius stopped the introduction of more transmissible VOCs at least nine times due to this approach. For nine months, Mauritius was considered to be COVID-19 safe with minimal restrictions and COVID-19 protocols imposed on local residents. Unfortunately, the evidence presented here suggests that despite these effective border controls, a single undetected introduction of a B.1.1.318 lineage virus was enough to initiate a large local outbreak. It is as yet unknown how this variant was introduced to Mauritius but it is likely to have occurred through a porous layer of the quarantine system for travelers; sequencing of specimens from travelers from January to March 2021 is ongoing in an effort to trace a probable source. Given that there were no direct flights between Finland, Ireland and Norway to support the introduction of B.1.1.318 directly from these countries, it is plausible that the introduction may have occurred through intermediate links (such as travel hubs) and not be reflected in this genomic surveillance data. The current data demonstrates the dominance of B.1.1.318 infections in the second wave. However, ongoing sequencing of cases from later in the outbreak will rule out the possibility of the co-circulation of other variants. While B.1.1.318 has been classified as a VUI by Public Health England (PHE), largely because of the presence of E484K, not much is known about this lineage other than it having three lineage-defining changes in the spike gene (Y144 Del, E484K and P681H) and one in nsp6 (106-108Del) that have convergently arisen in one or more of the better-known VOCs. All four of the defining amino acid substitutions in Spike (T95I, E484K, P681H and D796H) are at genome sites that have been detectably evolving under positive selection within the global SARS-CoV2 dataset since December 2020 when this lineage likely emerged. Further, the D796H mutation has been detected in a long-term SARS-CoV-2 infection, suggesting that it too might be significantly adaptive (19). Deletions within the NTD region of the Spike glycoprotein, especially in position 141-146, like the Y144Del seen in the B.1.1.318 virus, have been associated with resistance to several neutralizing antibodies but do not seem to confer significant reduction in vaccine efficacy (20). The P681H mutation, commonly present in the B.1.1.318 and the Alpha (B.1.1.7) variants, is proximal to the S1/S2 furin cleavage site of the Spike glycoprotein, a region known to impact transmissibility of the virus (12,21). Lastly, the E484K, found in the Receptor Binding Domain (RBD) of the Spike glycoprotein and also present in the B.1.1.318, the B.1.351, the P.1 and other variants, is known to cause escape to polyclonal antibodies and show decreased susceptibility both to convalescent and vaccine plasmas (14,22–24). Although the B.1.1.318 has been detected in more than 30 countries on four continents and is currently still being detected, this lineage had not been able to outcompete VOCs in other regions of the world until recently when it became dominant in Greece. Although it seeded a large outbreak in Mauritius, it was introduced at a time when no other variants were locally transmitting, so it is unclear whether specific biological properties of B.1.1.318 viruses (such as increased transmissibility) may have contributed to initiation and propagation of the outbreak. Immune evasion properties introduced by Spike glycoprotein substitutions (e.g. E484K and Y144del) are unlikely to be relevant in explaining the outbreak as there was presumably a low level of population immunity following a small first wave. This may be relevant when weighing the importance of controlling the outbreak, to ensure the best chance of success of the local vaccination program. The vaccination campaign in Mauritius started in January 2021 for front-line workers and vulnerable individuals and in March 2021 was extended to the general adult population. Vaccines being used are Astrazeneca/Covishield, Covaxin and Sinopharm. At the end of May 2021, close to 20% of the population of Mauritius had received the first dose and half of those had also received the second dose of one of the vaccines, but vaccination is not yet thought to have played a major role in slowing down the second wave in March and April 2021. In conclusion, SARS-CoV-2 genomic surveillance in Mauritius demonstrated a second wave that was dominated by a VUI, despite PCR testing and quarantine effectively preventing local transmission of other variants. This report also highlights the need for continuous genomic surveillance to fully to understand the transmission dynamics of the VUI in Mauritius. ## Supporting information Supplementary Table [[supplements/259017_file04.pdf]](pending:yes) ## Data Availability All of the genomes are available in GISAID. ## Methods ### Epidemiological data We analyzed COVID-19 cases and death counts in Mauritius from publicly released data up to 17 May 2021 from the Our World in Data (OWID) data repository ([https://github.com/owid/covid-19-data](https://github.com/owid/covid-19-data)). For correlation with government epidemic control measures, information was extracted from government press releases and speech transcripts. ### SARS-CoV-2 samples and metadata We obtained deidentified remnant nasopharyngeal and oropharyngeal swab samples from patients testing positive for SARS-CoV-2 by RT-qPCR from the Ministry of Health and Wellness of Mauritius primary testing laboratory, the Central Health Laboratory at Victoria Hospital. These were used for SARS-CoV-2 Whole Genome Sequencing (WGS). RNA was extracted using the MagMAX™ Viral/Pathogen II Nucleic Acid Isolation Kit (Cat. No. A48383) using the automated nucleic acid extractor KingFisher™ Apex Magnetic Particle Processor with 96 Deep-Well Head using sample volume input of 200 µL. Associated metadata for the samples included date and location (district) of sampling, and sex and age of the patients. ### Real Time RT-PCR In order to detect the SARS-CoV-2 virus by RT-PCR, the TaqPath COVID-19 CE-IVD RT-PCR Kit (Life Technologies, Carlsbad, CA) was used according to the manufacturer’s instructions. The assays target genomic regions (ORF1ab, S protein and N protein) of the SARS-CoV-2 genome. RT-PCR was performed on a QuantStudio 7 Flex Real-Time PCR instrument (Life Technologies, Carlsbad, CA). Cycle threshold (Ct) values were analyzed using auto-analysis settings with the threshold lines falling within the exponential phase of the fluorescence curves and above any background signal. ### WGS and genome assembly RNA was extracted manually using the Qiagen Viral RNA Mini Kit (QIAGEN, California, USA). cDNA synthesis was performed on the RNA using random primers followed by gene specific multiplex PCR using the ARTIC protocol (25). Briefly, extracted RNA was converted to cDNA using the Superscript IV First Strand synthesis system (Life Technologies, Carlsbad, CA) and random hexamer primers. SARS-CoV-2 whole genome amplification by multiplex PCR was carried out using primers designed on Primal Scheme ([http://primal.zibraproject.org/](http://primal.zibraproject.org/)) to generate 400bp amplicons with an overlap of 70bp that covers the 30Kb SARS-CoV-2 genome. PCR products were cleaned up using AmpureXP purification beads (Beckman Coulter, High Wycombe, UK) and quantified using the Qubit dsDNA High Sensitivity assay on the Qubit 4.0 instrument (Life Technologies Carlsbad, CA). The Illumina® Nextera Flex DNA Library Prep kit was used according to the manufacturer’s protocol to prepare uniquely indexed paired end libraries of genomic DNA. Sequencing libraries were normalized to 2nM, pooled and denatured with 0.2N sodium hydroxide. 1.3pM sample library was spiked with 1% PhiX (PhiX Control v3 adapter-ligated library used as a control). Libraries were loaded onto a 300-cycle NextSeq 500/550 Mid Output Kit v2. and run on the Illumina NextSeq550 instrument (Illumina, San Diego, CA, USA). Raw reads coming from Illumina sequencing were assembled using Exatype ([https://sars-cov-2.exatype.com/](https://sars-cov-2.exatype.com/)), Galaxy ARTIC pipeline (2. [https://covid19.galaxyproject.org/artic/](https://covid19.galaxyproject.org/artic/)), Genome Detective 1.126 ([https://www.genomedetective.com/](https://www.genomedetective.com/)) and the Coronavirus Typing Tool (26,27). The initial assembly obtained from Genome Detective was polished by aligning mapped reads to the references and filtering out low-quality mutations using bcftools 1.7-2 mpileup method. Seventeen of the samples in this study were sequenced at the Francis Crick Institute on ONT GridION following the ARTIC LoCost V3 protocol (28). Data was demultiplexed and processed using v1.1.1 of the ncov2019-artic-nf pipeline [[https://github.com/connor-lab/ncov2019-artic-nf](https://github.com/connor-lab/ncov2019-artic-nf)] with the additional parameters “--schemeVersion V3 -- minReadsPerBarcode 1 --minReadsArticGuppyPlex 1”. All of the sequences were deposited in GISAID ([https://www.gisaid.org/](https://www.gisaid.org/)) (29) and the GISAID accession included as part of the Supplementary Data (Table S3). ### Global reference dataset Mauritian sequences were analyzed against a backdrop of globally representative SARS-CoV-2 genotypes. At the time of sequence analysis, more than 1.5 million SARS-CoV-2 genomes had been shared publicly. Due to the sheer size of this dataset and over sampling in specific countries, a systematic down-sampling was carried out. Important lineage defining genomes along with randomly sampled genotypes per location were included in the phylogenetic reconstruction. The dataset was also enriched with sequences from the B.1.1.318 lineage. The final reference set of 4484 genomes contained 635 from Africa, 582 from Asia, 1690 from Europe, 735 and 465 from North and South America respectively and 377 genotypes from Oceania. ### Phylogenetic analysis of SARS-CoV-2 in Mauritius Mauritian sequences were analyzed against the global reference dataset using a custom build of the SARS-CoV-2 NextStrain pipeline ([https://github.com/nextstrain/ncov](https://github.com/nextstrain/ncov)) (30). The pipeline contains several python scripts that manage the analysis workflow. In short it allows for the filtering of genotypes, the alignment of genotypes in MAFFT (31), phylogenetic tree inference in IQ-Tree(32), tree dating and ancestral state construction and annotation. The resulting time scaled phylogeny can be viewed interactively and has been shared publicly on the NGS-SA NextStrain page ([https://nextstrain.org/groups/ngs-sa/COVID19-Mauritius-2021.05.10](https://nextstrain.org/groups/ngs-sa/COVID19-Mauritius-2021.05.10)). It was also visualized using R package ggtree (33). ### Lineage & Clade classification We used the dynamic lineage classification method proposed in this study via the Phylogenetic Assignment of named Global Outbreak LINeages (PANGOLin) software suite ([https://github.com/hCoV-2019/pangolin](https://github.com/hCoV-2019/pangolin)) (34). This was aimed at identifying the most epidemiologically important lineages of SARS-CoV-2 at the time of analysis, allowing researchers to monitor the epidemic in a specific geographical region. Accordingly, with this recently proposed dynamic lineage classification, many factors might suggest a new lineage including: i) monophyletic clusters on a global tree; ii) the presence of a statistically significant support (bootstrap/ultrafast bootstrap) on the node of the new lineages; iii) introduction into a novel geographic region; iv) epidemiological support (location; travel history) and v) characteristic Single Nucleotide Polymorphisms. Accordingly, three main SARS-CoV-2 lineages with those characteristics are currently recognized; lineage A, defined by Wuhan/WH04/2020, lineage B, defined by Wuhan-Hu-1 strain, and lineage C, a sub-classification from the B lineage. The BLAST function in GISAID (7) was also used to find sequences most similar to the Mauritius dataset and infer belonging to the B.1.1.318 lineage. We also classified the SARS-CoV-2 genomes in our dataset using the clade classification proposed by Nextstrain, divided into 19A, 19B, 20A, 20B, and 20C clades (35). ### Selection analysis To identify which, if any, of the observed mutations in the Spike protein was most likely to increase viral fitness, we used the natural selection analysis of SARS-CoV-2 pipeline ([https://observablehq.com/@spond/revised-sars-cov-2-analytics-page](https://observablehq.com/@spond/revised-sars-cov-2-analytics-page)) (36,37). This pipeline examines the entire global SARS-CoV-2 nucleotide sequence dataset for evidence of:(*) (i) polymorphisms having arisen in multiple epidemiologically unlinked lineages, (ii) the likelihood of these polymorphisms being associated with a greater than expected ratio of non-synonymous:synonymous nucleotide substitution rates, and (iii) whether these polymorphisms have increased in frequency in the regions of the world where they have occurred. ### Dated phylogenetics To estimate time-calibrated phylogenies dated from time-stamped genome data, we conducted phylogenetic analysis using the Bayesian software package BEASTv.1.10.4 (38), on a subset of data belonging to the B.1.1.318 lineage (n = 280) identified in the ML phylogeny to contain isolates from Mauritius(n=134). ML trees from this data subset was inspected in TempEst v1.5.3 (39) for the presence of a temporal (i.e. molecular clock) signal. Linear regression of root-to-tip genetic distances against sampling dates indicated that the SARS-CoV-2 sequences in this cluster evolve in a relatively-strong clock-like manner (r = 0.51; r2 = 0·3) (Supplementary Fig S4). For this analysis we employed the strict molecular clock model, the HKY+I, nucleotide substitution model and the exponential growth coalescent model. We computed MCMC (Markov chain Monte Carlo) triplicate runs of 100 million states each, sampling every 10.000 steps for each data set. Convergence of MCMC chains was checked using Tracer v.1.7.141 (40). Maximum clade credibility trees were summarised from the MCMC samples using TreeAnnotator after discarding 10% as burn-in. We discretized sequence sampling locations as North America, Europe, Asia, West Africa and Mauritius. ### Phylogeographic analysis To model phylogenetic diffusion of the B.1.1.318 cluster in Mauritius, we used a flexible relaxed random walk (RRW) diffusion model that accommodates branch-specific variation in rates of dispersal with a Cauchy distribution. For each sequence, latitude and longitude were attributed to a point randomly sampled within the patient’s district of residence. MCMC chains were run for >100 million generations and sampled every 10000th step, with convergence assessed using Tracer v1.7 (40). Maximum clade credibility trees were summarized using TreeAnnotator after discarding 10% as burn-in. We used the R package “seraphim” to extract and map spatiotemporal information embedded in posterior trees (41). ## Supplementary Figures ![Supplementary Fig S1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/16/2021.06.16.21259017/F3.medium.gif) [Supplementary Fig S1.](http://medrxiv.org/content/early/2021/06/16/2021.06.16.21259017/F3) Supplementary Fig S1. ML Timetree for three clusters of interest for introductions of SARS-CoV-2 to Mauritius (into quarantine centres): B.1.1.7 (A), 501Y.V2 (B), B.1.36 (C). ![Supplementary Fig S2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/16/2021.06.16.21259017/F4.medium.gif) [Supplementary Fig S2.](http://medrxiv.org/content/early/2021/06/16/2021.06.16.21259017/F4) Supplementary Fig S2. A) Geographical distribution of cases during the second wave in March 2021 (Picture property of Ministry of Health and Wellness of Mauritius and published by the Government Information Services on 2 April 2021) B) Location of sampling of genomes from the local (2nd) second wave in March 2021 (n=120). ![Supplementary Fig S3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/16/2021.06.16.21259017/F5.medium.gif) [Supplementary Fig S3.](http://medrxiv.org/content/early/2021/06/16/2021.06.16.21259017/F5) Supplementary Fig S3. Genome maps showing a comparison of mutations between the B.1.1.318 variant and two known VOCs, the B.1.1.7 and the B.1.351 ![Supplementary Fig S4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/16/2021.06.16.21259017/F6.medium.gif) [Supplementary Fig S4.](http://medrxiv.org/content/early/2021/06/16/2021.06.16.21259017/F6) Supplementary Fig S4. MCC tree of global B.1.1.318 sequences coloured by country ![Supplementary Fig S4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/16/2021.06.16.21259017/F7.medium.gif) [Supplementary Fig S4.](http://medrxiv.org/content/early/2021/06/16/2021.06.16.21259017/F7) Supplementary Fig S4. Root-to-tip regression plot of global B.1.1.318 lineage obtained through TempEst analysis showing the evolution of this lineage over time through accumulation of mutation (Distance) respective to the SARS-CoV-2 reference. ## Supplementary Tables View this table: [Supplementary Table S1.](http://medrxiv.org/content/early/2021/06/16/2021.06.16.21259017/T1) Supplementary Table S1. Ranking of the top 20 countries having sequenced the greatest proportions of SARS-CoV-2 cases relative to the total number of positive cases in the country. View this table: [Supplementary Table S2.](http://medrxiv.org/content/early/2021/06/16/2021.06.16.21259017/T2) Supplementary Table S2. Sites of B.1.1.318 evolving under positive selection View this table: [Supplementary Table S3.](http://medrxiv.org/content/early/2021/06/16/2021.06.16.21259017/T3) Supplementary Table S3. Genomes from Mauritius on GISAID. ## Acknowledgments We are grateful to the Ministry of Health and Wellness of Mauritius for providing samples for this study. Equally, we are grateful to Mungra Priscilla, Mathur Hari, Sujeewon Chitra and all the staff of the Virology and Molecular Biology departments of the Victoria Hospital Mauritius for their dedication and hard work. We would also like to thank our colleagues Professor Aris Katzourakis from the University of Oxford and Dr. Gkikas Magiorkinis from National and Kapodistrian University of Athens for insightful discussions on the epidemic in Greece. We are immensely grateful to the laboratories and institutions around the world which have sampled and sequenced B.1.1.318 variants in their respective countries and made the data publicly available on GISAID (Supp Table S4), and to Public Health England for initially flagging this variant. These were all crucial in contextualizing our study. We also acknowledge Harshil Patel, Robert Goldstone, John Mccauley and Jerome Nicod at the Francis Crick Institute for assistance in sample management and data processing. ## Footnotes * * joint first authors * ¶ jointly supervised this work * Received June 16, 2021. * Revision received June 16, 2021. * Accepted June 16, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Alteri C, Cento V, Piralla A, Costabile V, Tallarita M, Colagrossi L, et al. Genomic epidemiology of SARS-CoV-2 reveals multiple lineages and early spread of SARS- CoV-2 infections in Lombardy, Italy. Nat Commun [Internet]. 2021;12(1):434. Available from: [https://doi.org/10.1038/s41467-020-20688-x](https://doi.org/10.1038/s41467-020-20688-x) 2. 2.Gonzalez-Reiche AS, Hernandez MM, Sullivan M, Ciferri B, Alshammary H, Obla A, et al. Introductions and early spread of SARS-CoV-2 in the New York City area. medRxiv. 2020 Apr 16;2020.04.08.20056929. 3. 3.Deng X, Gu W, Federman S, du Plessis L, Pybus OG, Faria N, et al. Genomic surveillance reveals multiple introductions of SARS-CoV-2 into Northern California. Science (80-) [Internet]. 2020 Jun 8;eabb9263. Available from: [http://science.sciencemag.org/content/early/2020/06/05/science.abb9263.abstract](http://science.sciencemag.org/content/early/2020/06/05/science.abb9263.abstract) 4. 4.du Plessis L, McCrone JT, Zarebski AE, Hill V, Ruis C, Gutierrez B, et al. Establishment and lineage dynamics of the SARS-CoV-2 epidemic in the UK. Science (80-) [Internet]. 2021 Jan 8;eabf2946. Available from: [http://science.sciencemag.org/content/early/2021/01/07/science.abf2946.abstract](http://science.sciencemag.org/content/early/2021/01/07/science.abf2946.abstract) 5. 5.Tegally H, Wilkinson E, Lessells RJ, Giandhari J, Pillay S, Msomi N, et al. Sixteen novel lineages of SARS-CoV-2 in South Africa. Nat Med [Internet]. 2021;27(3):440–6. Available from: [https://doi.org/10.1038/s41591-021-01255-3](https://doi.org/10.1038/s41591-021-01255-3) 6. 6.Lu J, du Plessis L, Liu Z, Hill V, Kang M, Lin H, et al. Genomic Epidemiology of SARS-CoV-2 in Guangdong Province, China. Cell [Internet]. 2020 Apr 30 [cited 2020 May 8]; Available from: [https://linkinghub.elsevier.com/retrieve/pii/S0092867420304864](https://linkinghub.elsevier.com/retrieve/pii/S0092867420304864) 7. 7.Jasper J. Koehorst, Jesse C. J. van Dam, Edoardo Saccenti, Vitor A. P. Martins dos Santos MS-D and PJS. GISAID Global Initiative on Sharing All Influenza Data. Phylogeny of SARS-like betacoronaviruses including novel coronavirus (nCoV). Oxford [Internet]. 2017 [cited 2020 May 8];34(8):1401–3. Available from: [https://www.researchgate.net/publication/321428915\_SAPP\_functional\_genome\_annotation\_and\_analysis\_through\_a\_semantic\_framework\_using\_FAIR\_principles](https://www.researchgate.net/publication/321428915\_SAPP\_functional\_genome\_annotation\_and\_analysis\_through\_a\_semantic\_framework_using_FAIR_principles) 8. 8.Preliminary genomic characterisation of an emergent SARS-CoV-2 lineage in the UK defined by a novel set of spike mutations - SARS-CoV-2 coronavirus / nCoV-2019 Genomic Epidemiology - Virological [Internet]. [cited 2021 May 25]. Available from: [https://virological.org/t/preliminary-genomic-characterisation-of-an-emergent-sars-cov-2-lineage-in-the-uk-defined-by-a-novel-set-of-spike-mutations/563](https://virological.org/t/preliminary-genomic-characterisation-of-an-emergent-sars-cov-2-lineage-in-the-uk-defined-by-a-novel-set-of-spike-mutations/563) 9. 9.Tegally H, Wilkinson E, Giovanetti M, Iranzadeh A, Fonseca V, Giandhari J, et al. Detection of a SARS-CoV-2 variant of concern in South Africa. Nature [Internet]. 2021;592(7854):438–43. Available from: [https://doi.org/10.1038/s41586-021-03402-9](https://doi.org/10.1038/s41586-021-03402-9) 10. 10.Faria NR, Mellan TA, Whittaker C, Claro IM, Candido D da S, Mishra S, et al. Genomics and epidemiology of the P.1 SARS-CoV-2 lineage in Manaus, Brazil. Science (80-) [Internet]. 2021 May 21;372(6544):815 LP–821. Available from: [http://science.sciencemag.org/content/372/6544/815.abstract](http://science.sciencemag.org/content/372/6544/815.abstract) 11. 11.Abdool Karim SS, de Oliveira T. New SARS-CoV-2 Variants — Clinical, Public Health, and Vaccine Implications. N Engl J Med [Internet]. 2021 Mar 24;384(19):1866–8. Available from: [https://doi.org/10.1056/NEJMc2100362](https://doi.org/10.1056/NEJMc2100362) 12. 12.Davies NG, Abbott S, Barnard RC, Jarvis CI, Kucharski AJ, Munday JD, et al. Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England. Science (80-) [Internet]. 2021 Apr 9;372(6538):eabg3055. Available from: [http://science.sciencemag.org/content/372/6538/eabg3055.abstract](http://science.sciencemag.org/content/372/6538/eabg3055.abstract) 13. 13.Challen R, Brooks-Pollock E, Read JM, Dyson L, Tsaneva-Atanasova K, Danon L. Risk of mortality in patients infected with SARS-CoV-2 variant of concern 202012/1: matched cohort study. BMJ [Internet]. 2021 Mar 10;372:579. Available from: [http://www.bmj.com/content/372/bmj.n579.abstract](http://www.bmj.com/content/372/bmj.n579.abstract) 14. 14.Cele S, Gazy I, Jackson L, Hwa S-H, Tegally H, Lustig G, et al. Escape of SARS- CoV-2 501Y.V2 from neutralization by convalescent plasma. medRxiv [Internet]. 2021 Jan 1;2021.01.26.21250224. Available from: [http://medrxiv.org/content/early/2021/02/27/2021.01.26.21250224.abstract](http://medrxiv.org/content/early/2021/02/27/2021.01.26.21250224.abstract) 15. 15.O’Toole LJ, Hill V, Pybus OG, Watts A, Bogoch II, Khan K, et al. Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2 [version 1; peer review: awaiting peer review]. Wellcome Open Res [Internet]. 2021;6(121). Available from: [https://wellcomeopenresearch.org/articles/6-121/v1](https://wellcomeopenresearch.org/articles/6-121/v1) 16. 16.Grout L, Katar A, Ouakrim DA, Summers JA, Kvalsvig A, Baker MG, et al. Estimating the Failure Risk of Quarantine Systems for Preventing COVID-19 Outbreaks in Australia and New Zealand. medRxiv [Internet]. 2021 Jan 1;2021.02.17.21251946. Available from: [http://medrxiv.org/content/early/2021/04/30/2021.02.17.21251946.abstract](http://medrxiv.org/content/early/2021/04/30/2021.02.17.21251946.abstract) 17. 17.Kowlessur S, Ori B, Zimmet P, Tuomilehto J, Chitson P, Ramphul Y. Tackling the COVID-19 pandemic in paradise: the Mauritian experience. Lancet Diabetes Endocrinol [Internet]. 2020 Nov 1;8(11):878–9. Available from: [https://doi.org/10.1016/S2213-8587(20)30336-3](https://doi.org/10.1016/S2213-8587(20)30336-3) 18. 18.SARS-CoV-2 variants of concern and variants under investigation in England. 19. 19.Kemp SA, Collier DA, Datir RP, Ferreira IATM, Gayed S, Jahun A, et al. SARS-CoV- 2 evolution during treatment of chronic infection. Nature [Internet]. 2021;592(7853):277–82. Available from: [https://doi.org/10.1038/s41586-021-03291-y](https://doi.org/10.1038/s41586-021-03291-y) 20. 20.McCallum M, Marco A De, Lempp F, Tortorici MA, Pinto D, Walls AC, et al. N- terminal domain antigenic mapping reveals a site of vulnerability for SARS-CoV-2. bioRxiv [Internet]. 2021 Jan 1;2021.01.14.426475. Available from: [http://biorxiv.org/content/early/2021/01/14/2021.01.14.426475.abstract](http://biorxiv.org/content/early/2021/01/14/2021.01.14.426475.abstract) 21. 21.Peacock TP, Goldhill DH, Zhou J, Baillon L, Frise R, Swann OC, et al. The furin cleavage site of SARS-CoV-2 spike protein is a key determinant for transmission due to enhanced replication in airway cells. bioRxiv [Internet]. 2020 Jan 1;2020.09.30.318311. Available from: [http://biorxiv.org/content/early/2020/09/30/2020.09.30.318311.abstract](http://biorxiv.org/content/early/2020/09/30/2020.09.30.318311.abstract) 22. 22.Jangra S, Ye C, Rathnasinghe R, Stadlbauer D, Krammer F, Simon V, et al. The E484K mutation in the SARS-CoV-2 spike protein reduces but does not abolish neutralizing activity of human convalescent and post-vaccination sera. medRxivLJ: the preprint server for health sciences. 2021. 23. 23.Garcia-Beltran WF, Lam EC, St. Denis K, Nitido AD, Garcia ZH, Hauser BM, et al. Multiple SARS-CoV-2 variants escape neutralization by vaccine-induced humoral immunity. Cell [Internet]. 2021; Available from: [https://www.sciencedirect.com/science/article/pii/S0092867421002981](https://www.sciencedirect.com/science/article/pii/S0092867421002981) 24. 24.Greaney AJ, Starr TN, Barnes CO, Weisblum Y, Schmidt F, Caskey M, et al. Mutational escape from the polyclonal antibody response to SARS-CoV-2 infection is largely shaped by a single class of antibodies. bioRxiv [Internet]. 2021 Jan 1;2021.03.17.435863. Available from: [http://biorxiv.org/content/early/2021/03/18/2021.03.17.435863.abstract](http://biorxiv.org/content/early/2021/03/18/2021.03.17.435863.abstract) 25. 25.Quick J. ARTIC Coronavirus Method Development Community 1 more workspace [Internet]. 2020 Aug [cited 2020 Nov 19]. Available from: [https://protocols.io/view/ncov-2019-sequencing-protocol-v3-locost-bh42j8ye](https://protocols.io/view/ncov-2019-sequencing-protocol-v3-locost-bh42j8ye) 26. 26.Cleemput S, Dumon W, Fonseca V, Abdool Karim W, Giovanetti M, Alcantara LC, et al. Genome Detective Coronavirus Typing Tool for rapid identification and characterization of novel coronavirus genomes. Bioinformatics. 2020;36(11):3552–5. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btaa145&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32108862&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F16%2F2021.06.16.21259017.atom) 27. 27.Vilsker M, Moosa Y, Nooij S, Fonseca V, Ghysens Y, Dumon K, et al. Genome Detective: an automated system for virus identification from high-throughput sequencing data. Bioinformatics. 2019;35(5):871–3. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bty695&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30124794&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F16%2F2021.06.16.21259017.atom) 28. 28.Pliner HA, Packer JS, McFaline-Figueroa JL, Cusanovich DA, Daza RM, Aghamirzaie D, et al. Cicero Predicts cis-Regulatory DNA Interactions from Single- Cell Chromatin Accessibility Data. Mol Cell [Internet]. 2018 Sep 6;71(5):858-871.e8. Available from: [https://doi.org/10.1016/j.molcel.2018.06.044](https://doi.org/10.1016/j.molcel.2018.06.044) 29. 29.Shu Y, McCauley J. GISAID: Global initiative on sharing all influenza data–from vision to reality. Eurosurveillance. 2017;22(13):30494. 30. 30. James Hadfield, Colin Megill, Sidney M Bell, John Huddleston, Barney Potter, Charlton Callender, Pavel Sagulenko, Trevor Bedford RAN. Nextstrain: real-time tracking of pathogen evolution. Bioinformatics [Internet]. 2018 [cited 2020 May 8];Volume 34(Issue 23):4121–4123. Available from: [https://www.researchgate.net/publication/325344695\_Nextstrain\_real-time\_tracking\_of\_pathogen\_evolution](https://www.researchgate.net/publication/325344695\_Nextstrain\_real-time\_tracking_of_pathogen_evolution) 31. 31.Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol [Internet]. 2013/01/16. 2013 Apr;30(4):772–80. Available from: [https://pubmed.ncbi.nlm.nih.gov/23329690](https://pubmed.ncbi.nlm.nih.gov/23329690) 32. 32.Nguyen L-T, Schmidt HA, Von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol. 2015;32(1):268–74. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msu300&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25371430&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F16%2F2021.06.16.21259017.atom) 33. 33.Yu G. Using ggtree to Visualize Data on Tree-Like Structures. Curr Protoc Bioinforma. 2020 Mar;69(1):e96. 34. 34.Rambaut A, Holmes EC, O’Toole Á, Hill V, McCrone JT, Ruis C, et al. A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology. Nat Microbiol [Internet]. 2020 Jul 15 [cited 2020 Oct 1];1–5. Available from: [https://doi.org/10.1038/s41564-020-0770-5](https://doi.org/10.1038/s41564-020-0770-5) 35. 35.Year-letter Genetic Clade Naming for SARS-CoV-2 on Nextstain.org [Internet]. [cited 2020 Oct 4]. Available from: [https://nextstrain.org/blog/2020-06-02-SARSCoV2-clade-naming](https://nextstrain.org/blog/2020-06-02-SARSCoV2-clade-naming) 36. 36.Murrell B, Wertheim JO, Moola S, Weighill T, Scheffler K, Kosakovsky Pond SL. Detecting Individual Sites Subject to Episodic Diversifying Selection. PLOS Genet [Internet]. 2012 Jul 12;8(7):e1002764. Available from: [https://doi.org/10.1371/journal.pgen.1002764](https://doi.org/10.1371/journal.pgen.1002764) 37. 37.Kosakovsky Pond SL, Frost SDW. Not So Different After All: A Comparison of Methods for Detecting Amino Acid Sites Under Selection. Mol Biol Evol [Internet]. 2005 May 1;22(5):1208–22. Available from: [https://doi.org/10.1093/molbev/msi105](https://doi.org/10.1093/molbev/msi105) 38. 38.Suchard MA, Lemey P, Baele G, Ayres DL, Drummond AJ, Rambaut A. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol. 2018 Jan;4(1):vey016. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ve/vey016&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29942656&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F16%2F2021.06.16.21259017.atom) 39. 39.Rambaut A, Lam TT, Max Carvalho L, Pybus OG. Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen). Virus Evol. 2016 Jan;2(1):vew007. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ve/vew007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27774300&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F16%2F2021.06.16.21259017.atom) 40. 40.Rambaut A, Drummond AJ, Xie D, Baele G, Suchard MA. Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7. Syst Biol. 2018 Sep;67(5):901–4. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/sysbio/syy032&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29718447&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F16%2F2021.06.16.21259017.atom) 41. 41.Dellicour S, Rose R, Faria NR, Lemey P, Pybus OG. SERAPHIM: studying environmental rasters and phylogenetically informed movements. Bioinformatics. 2016 Oct;32(20):3204–6. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btw384&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27334476&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F16%2F2021.06.16.21259017.atom)