# Analysis of SARS-CoV-2 Variants From 24,181 Patients Exemplifies the Role of Globalization and Zoonosis in Pandemics

Philippe Colson, Pierre-Edouard Fournier, Hervé Chaudet, Jérémy Delerce, Audrey Giraud-Gatineau, Linda Houhamdi, Claudia Andrieu, Ludivine Brechard, Marielle Bedotto, Elsa Prudent, et al.

► To cite this version:

Philippe Colson, Pierre-Edouard Fournier, Hervé Chaudet, Jérémy Delerce, Audrey Giraud-Gatineau, et al. Analysis of SARS-CoV-2 Variants From 24,181 Patients Exemplifies the Role of Globalization and Zoonosis in Pandemics. Frontiers in Microbiology, 2022, 12, 10.3389/fmicb.2021.786233. hal-03663729

HAL Id: hal-03663729

https://amu.hal.science/hal-03663729

Submitted on 21 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Analysis of SARS-CoV-2 Variants From 24,181 Patients Exemplifies the Role of Globalization and Zoonosis in Pandemics

OPEN ACCESS

**Edited by:**

Ashley C. Banyard,
Animal and Plant Health Agency,
United Kingdom

**Reviewed by:**

Yongming Sang,
Tennessee State University,
United States

M. Nazmul Hoque,
Bangabandhu Sheikh Mujibur
Rahman Agricultural University,
Bangladesh

**\*Correspondence:**

Didier Raoult
didier.raoult@gmail.com

† These authors have contributed equally to this work

**Specialty section:**

This article was submitted to Virology, a section of the journal Frontiers in Microbiology

**Received:** 30 September 2021

**Accepted:** 15 December 2021

**Published:** 07 February 2022

**Citation:**

Colson P, Fournier P-E, Chaudet H, Delerce J, Giraud-Gatineau A, Houhamdi L, Andrieu C, Brechard L, Bedotto M, Prudent E, Gazin C, Beye M, Burel E, Dudouet P, Tissot-Dupont H, Gautret P, Lagier J-C, Million M, Brouqui P, Parola P, Fenollar F, Drancourt M, La Scola B, Levasseur A and Raoult D (2022) Analysis of SARS-CoV-2 Variants From 24,181 Patients Exemplifies the Role of Globalization and Zoonosis in Pandemics. *Front. Microbiol.* 12:786233. doi: 10.3389/fmicb.2021.786233

Philippe Colson¹,²,³†, Pierre-Edouard Fournier¹,²,³†, Hervé Chaudet¹,⁴,⁵, Jérémy Delerce¹, Audrey Giraud-Gatineau¹,³,⁴,⁵, Linda Houhamdi¹, Claudia Andrieu¹, Ludivine Brechard¹, Marielle Bedotto¹, Elsa Prudent¹, Céline Gazin¹, Mamadou Beye¹, Emilie Burel¹, Pierre Dudouet¹,²,³, Hervé Tissot-Dupont¹,²,³, Philippe Gautret¹,³,⁴, Jean-Christophe Lagier¹,²,³, Matthieu Million¹,²,³, Philippe Brouqui¹,²,³, Philippe Parola¹,³,⁴, Florence Fenollar¹,³,⁴, Michel Drancourt¹,²,³, Bernard La Scola¹,²,³, Anthony Levasseur¹,²,³ and Didier Raoult¹,²,³*

¹ IHU Méditerranée Infection, Marseille, France, ² Microbes Evolution Phylogeny and Infections (MEPHI), Institut de Recherche pour le Développement, Aix-Marseille Université, Marseille, France, ³ Assistance Publique-Hôpitaux de Marseille, Marseille, France, ⁴ Vecteurs-Infections Tropicales et Méditerranéennes (VITROME), Institut de Recherche pour le Développement, Aix-Marseille Université, Marseille, France, ⁵ French Armed Forces Center for Epidemiology and Public Health, Marseille, France

After the end of the first epidemic episode of SARS-CoV-2 infections, as cases began to rise again during the summer of 2020, we at IHU Méditerranée Infection in Marseille, France, intensified the genomic surveillance of SARS-CoV-2, and described the first viral variants. In this study, we compared the incidence curves of SARS-CoV-2-associated deaths in different countries and reported the classification of SARS-CoV-2 variants detected in our institute, as well as the kinetics and sources of the infections. We used mortality collected from a COVID-19 data repository for 221 countries. Viral variants were defined based on ≥5 hallmark mutations along the whole genome shared by ≥30 genomes. SARS-CoV-2 genotype was determined for 24,181 patients using next-generation genome and gene sequencing (in 47 and 11% of cases, respectively) or variant-specific qPCR (in 42% of cases). Sixteen variants were identified by analyzing viral genomes from 9,788 SARS-CoV-2-diagnosed patients. Our data show that since the first SARS-CoV-2 epidemic episode in Marseille, importation through travel from abroad was documented for seven of the new variants. In addition, for the B.1.160 variant of Pangolin classification (a.k.a. Marseille-4), we suspect transmission from farm minks. In conclusion, we observed that the successive epidemic peaks of SARS-CoV-2 infections are not linked to rebounds of viral genotypes that are already present but to newly introduced variants. We thus suggest that border control is the best mean of combating this type of introduction, and that intensive control of mink farms is also necessary to prevent the emergence of new variants generated in this animal reservoir.

**Keywords:** SARS-CoV-2, variant, mutant, classification, travel, zoonosis, epidemics, pandemics
# INTRODUCTION

Following its emergence in December 2019 in Wuhan, China, SARS-CoV-2 was declared to be a pandemic in March 2020 and its incidence since then has shown different epidemic patterns according to the country<fnref n="1" />,<fnref n="2" /> (Cucinotta and Vanelli, 2020). In Asia, a bell-shaped curve typical of seasonal viral respiratory infections has been observed. In Western developed countries, successive phases have occurred that have associated an initial bell-shaped curve and additional waves with one or more peaks. In Marseille, France, after a first epidemic period characterized by a bell-shaped curve of incidence followed by a period of almost zero incidence in May and June 2020, an epidemic burst of SARS-CoV-2 cases occurred in July. As this differed from the single bell-shaped curve typically observed for most respiratory viral infections, we suspected this corresponded to the introduction of new SARS-CoV-2 genotypes. Consequently, in the summer of 2020 we expanded our assessment of viral genetic diversity using genome next-generation sequencing. This allowed us to detect and report as early as 7 September 2020 several SARS-CoV-2 lineages that occurred concurrently or successively, and were associated with different genetic, epidemiological and clinical features (Colson et al., 2020b, 2021c; Fournier et al., 2021).

These observations required defining and naming these different lineages to communicate about them, which is the very principle of taxonomy. Classifying and naming emerging infectious agents is essential (Sokal, 1974) but the task is particularly difficult in the case of SARS-CoV-2 for several reasons. First, as for other RNA viruses, SARS-CoV-2 evolves through continuous genetic variation with a high mutation rate (Eigen, 1996; Domingo and Perales, 2019). Second, SARS-CoV-2 has caused a dramatically high number of infections (over 195 million human cases worldwide and 5.9 million cases in France by 31 July 2021 (see footnote 2). It has been estimated that the SARS-CoV-2 genetic diversity in an infected patient may cover every single nucleotide substitutions and every possible pairs of nucleotide substitutions (Sender et al., 2021). In addition, there are animal reservoirs, a major one being mink (Fenollar et al., 2021; Oude Munnink et al., 2021). SARS-CoV-2 is, therefore, currently a rapidly evolving virus with a mutation rate estimated to be 9.8 × 10⁻<fnref n="4" /> substitutions/site/year (van Dorp et al., 2020) and characterized by a high rate of lineage turnover (Rambaut et al., 2020). Several SARS-CoV-2 specific classification and naming systems exist including those of GISAID<fnref n="3" /> (Alm et al., 2020), Nextstrain<fnref n="4" /> (Hadfield et al., 2018; Aksamentov et al., 2021), Pangolin<fnref n="5" /> (Rambaut et al., 2020), or the World Health Organization (WHO)<fnref n="6" />. But, although these are of great interest and usefulness, other classification and naming systems are also valuable, as previously acknowledged (Rambaut et al., 2020), particularly when it comes to better exploring

genomic epidemiology on a local scale through real-time next-generation sequencing. Indeed, some classifications have become complicated with the advent of multiple new variants. In addition, local genetic trends may be underrepresented in the tremendous flow of available genomes. Moreover, correlations between genotypic patterns and epidemiological and clinical characteristics can be accurately performed at institutional or local levels through a large network of people who are responsible for patient diagnosis and clinical management, while remote collection and analysis of data may have poorer performance. Thus, we aimed at using our own definition, classification and naming of SARS-CoV-2 lineages and at correlating these lineages with our local epidemiological data in order to better investigate and monitor their genetic features and their sources and geographical origins. Here, we describe our classification and naming system of SARS-CoV-2 lineages and the history and characteristics of the major lineages detected among patients diagnosed at our institute, which sharply highlights the role of travel and globalization in the current pandemic.

# RESULTS

## Broad Diversity According to Countries of the Patterns of Incidence of SARS-CoV-2-Associated Deaths

Clustering of the incidence curves of SARS-CoV-2-associated deaths for 221 countries delineated ten major groups of countries (Figure 1A). Each group comprised a mean ± standard deviation of 12 ± 4 countries (range, 7–21). A broad range of patterns was revealed through these groups, with countries experiencing one or more epidemic waves of varying intensity and duration combined or not with periods of lower and uneven incidence (Figure 1B). These data indicate that the incidence of SARS-CoV-2-associated deaths did not follow a single bell-shaped curve except in countries from Group 1, which includes China, Australia and a few countries in Central Asia and the Middle East. Instead, the incidence of SARS-CoV-2-associated deaths evolved elsewhere with multiple epidemic patterns according to time and geographic location. The worldwide distribution of these patterns showed that the majority of Western and Northern European countries (including France), the United States and Canada, which are part of Group 10, shared a same pattern (Figures 1A,C). Similarly, several countries in Asia and South America shared the same patterns leading to their classification in Groups 1–4. In contrast, a broad panel of patterns was observed for Africa. Countries that experienced a single peak of SARS-CoV-2-associated deaths (Groups 1, 4, 5, and 7) were mostly in Asia, Oceania, the Middle East, South America, Africa, and Eastern Europe; exceptions were Luxembourg, Austria and Switzerland. Among the countries which experienced several peaks of SARS-CoV-2-associated deaths, those in Group 10 were the developed industrialized countries of Western Europe and North America. These countries are particularly involved in globalization through tourism, business, and immigration (Chinazzi et al., 2020). In Marseille and its geographical area,

<footnote n="1">¹https://ourworldindata.org/</footnote>

<footnote n="2">²https://www.ecdc.europa.eu/en/geographical-distribution-2019-ncov-cases</footnote>

<footnote n="3">³https://www.gisaid.org/</footnote>

<footnote n="4">⁴https://clades.nextstrain.org/</footnote>

<footnote n="5">⁵https://cov-lineages.org/pangolin.html</footnote>

<footnote n="6">⁶https://www.who.int/en/activities/tracking-SARS-CoV-2-variants/</footnote>
the overall pattern of incidence of SARS-CoV-2-associated deaths was similar to that of countries in Group 10, which includes France, with a bell-shaped curve in March–April 2020 then a short period with a very low incidence followed by a longer period with two peaks of incidence and erratic incidence between these peaks (Figure 2A). This pattern of incidence with several peaks suggests multiple introductions of new variants that evolved abroad and caused distinct epidemics.

## Increase of SARS-CoV-2 Genetic Diversity Over Time

A total of 54,703 (10.6%) out of 513,805 patients tested in our institute for SARS-CoV-2 infection between 29 January 2020 and 18 August 2021 were positive by real-time reverse transcription-PCR (qPCR)<fnref n="7" /> (Figure 2A). We observed a bell-shaped curve of incidence between late February and early May 2020. This was followed by a period with almost no diagnoses (between 0 and 8 per day; mean ± standard deviation, 2.5 ± 2.2 per day) which lasted from 9 May to 5 July 2020. We performed whole SARS-CoV-2 genome sequencing using next-generation sequencing procedures since we diagnosed the first infection on 27 February 2020 (Colson et al., 2020a) and primarily obtained viral genomes from 309 patients for the March-April period (Colson et al., 2020b; Levasseur et al., 2020). In July 2020, the incidence of SARS-CoV-2 diagnoses increased again (Figure 2A). As the occurrence of a second epidemic 2 months after the end of the first one was unexpected for a respiratory virus infection, we suspected that it was linked to the introduction of new SARS-CoV-2 genotypes. This led us to quickly expand the next-generation sequencing of SARS-CoV-2 genomes to test this hypothesis. We primarily obtained 382 additional genomes during summer 2020 (Colson et al., 2020b), and ultimately analyzed the viral genomes obtained from a total of 9,788 patients sampled between 27 February 2020 and 18 July 2021. We observed a significant increase in the mean number of mutations in SARS-CoV-2 genomes between April 2020 and July 2020 [mean (±standard deviation) number of mutations per genome: 8.5 ± 2.0 versus 15.0 ± 4.4, respectively; $p < 10^{-3}$ (ANOVA test)] (Figure 2B and Supplementary Table 1). This mean number continued to increase during the summer of 2020, and the increase was significant between July and August 2020 (20.0 ± 5.0). Thereafter, significant increases in the mean number of mutations per genome were observed between September and December 2020, (25.2 ± 4.3), between November 2020 (23.0 ± 4.1) and February 2021 (26.4 ± 5.7), and between February and March 2021 (34.1 ± 4.2). Similar trends were observed for the number of amino acid changes in the spike protein (Supplementary Figure 1). Overall, when analyzing the pairwise genetic distance between genomes obtained over time, two major periods with increased viral genetic diversity compared to previous months were identified, during the summer of 2020 compared to before June 2020, mean ± standard deviation of pairwise genetic distances between genomes being $8.28 \times 10^{-4} \pm 3.50 \times 10^{-4}$ versus $2.34 \times 10^{-4} \pm 1.24 \times 10^{-4}$, respectively ($p < 10^{-3}$), then in January-March 2021 when

the mean ± standard deviation pairwise genetic distances between genomes was $13.93 \times 10^{-4} \pm 4.48 \times 10^{-4}$ compared to $7.14 \times 10^{-4} \pm 4.44 \times 10^{-4}$ in October–December 2020 ($p < 10^{-3}$) (Figure 3).

## The Marseille IHU Méditerranée Infection Classification System of SARS-CoV-2 Variants

Since the summer of 2020, we intensified our activity of SARS-CoV-2 genotyping using next-generation sequencing and also through performing partial spike gene sequencing and implementing variant-specific qPCR (Bedotto et al., 2021a,b). To characterize the increasing SARS-CoV-2 genome diversity observed during the summer of 2020, we promptly implemented a viral genome classification system and a nomenclature of SARS-CoV-2 variants. The term “SARS-CoV-2 variant” had scarcely been used in the scientific literature before September 2020 (Supplementary Figure 2) and no definition had been proposed to differentiate variants from mutants. We first reported the emergence of SARS-CoV-2 variants on 7 September, 2020 (Colson et al., 2020b), which included Marseille-2 (later named B.1.177 according to Pangolin classification) and Marseille-4 (B.1.160) variants that were then reported in October 2020 to have emerged and, for the case of Marseille-2, to have spread across Europe (Hodcroft et al., 2021).

There are no universally recognized or used strategies and criteria to classify viruses below the species level, and only biological classifications and nomenclatures have been proposed (Fauquet et al., 2008; International Committee on Taxonomy of Viruses (ICTV), 2020; Box 1). Classification and naming systems have been specifically proposed for SARS-CoV-2 including those of GISAID (see footnote 3) (Alm et al., 2020), Nextstrain (see footnote 4) (Hadfield et al., 2018; Aksamentov et al., 2021), Pangolin (see footnote 5) (Rambaut et al., 2020), and the WHO (see footnote 6) (Box 2). However, we chose to implement our own classification system and nomenclature and to use names rather than only numbers. This approach fitted the viral genomic epidemiology in our geographical area, and was more understandable and easier to use, which is the very purpose of classification and nomenclature. Our strategy consisted in detecting all mutations in all genomes obtained from patients diagnosed in our institute and in assigning a dynamic nomenclature to track SARS-CoV-2 genetic patterns.

Phylogenetic trees were built that included all genomic sequences of SARS-CoV-2 obtained in our institute (Figure 4 and Supplementary Figure 3). Based on tree topology and to distinguish them from SARS-CoV-2 mutants, we delineated SARS-CoV-2 variants as groups of viral genomes carrying a set of at least five mutations along the whole genome differentiating them from any other viral genomes and obtained from at least 30 (initially 10) different patients. Genomes were, therefore, unambiguously assigned to a given variant based on the presence of particular sets of mutations (Table 1). This differentiated variants from mutants, the genomes of which do not harbor this minimal set of hallmark mutations. We first delineated seven Marseille variants of SARS-CoV-2 as early as 7 September 2020, 2 of which did not eventually reach the threshold number

<footnote n="7">⁷https://www.mediterranee-infection.com/covid-19/</footnote>
of 30 genomes (Colson et al., 2020b). This was the first description of an emergence and expansion of several SARS-CoV-2 variants responsible for distinct epidemics in a same geographical area. Three additional variants were then defined in 2020 (Supplementary Figure 4) and four in 2021, making a total of 12 (Table 1, Figure 4, and Supplementary Figure 4).

In addition, we delineated 9 other viral lineages comprising genomes harboring a specific set of at least five mutations but obtained from < 30 different patients; some of these lineages may grow and become bona fide variants as additional members will be obtained. Besides, we detected four other SARS-CoV-2 variants that were defined as variants of concern (VOC; Box 2)

<figure><img src="image_2.png" /><figcaption>FIGURE 1 | (Continued)</figcaption></figure>
and were named Alpha, Beta, Gamma and Delta by the WHO (see footnote 6). Regarding the nomenclature of the variants we defined, we used our city name, Marseille, followed by a number according to the chronology of their description (e.g., Marseille-1, Marseille-2, ...) or one of their hallmark amino acid substitutions (e.g., Marseille-501). Hereafter, we will use both the Marseille and Pangolin nomenclatures to describe the genotypes of the SARS-CoV-2 genomes obtained in our institute.

# Emergence and Outcome of the SARS-CoV-2 Variants That Circulated in Our Geographical Area

## Overall Prevalence of the Marseille SARS-CoV-2 Variants

SARS-CoV-2 genotyping was performed and successful for 24,181 (44%) of the 54,703 SARS-CoV-2-positive patients

<figure><img src="image_3.png" /></figure>
C

<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | Patterns of incidence curves of SARS-CoV-2-associated deaths according to countries. (A) Hierarchical clustering of incidence curves of SARS-CoV-2-associated deaths per country. (B) Major patterns of incidence curves of SARS-CoV-2-associated deaths defined based on hierarchical clustering, and the example of one country for each pattern. Patterns of incidence curves of SARS-CoV-2-associated deaths are shown in 10 panels numbered from 1 to 10; these numbers correspond to those from the cladogram of panel (A). The example of one country is shown in a panel at the right of each of the ten patterns. X-axis indicates time in numbers of days, and Y-axis indicates normalized daily deaths. (C) Worldwide distribution of the 10 major patterns of incidence curves of SARS-CoV-2-associated deaths defined based on hierarchical clustering. Colors are as used for panel (A).</figcaption></figure>

A

<figure><img src="image_5.png" /></figure>

B

<figure><img src="image_6.png" /><figcaption>**FIGURE 2** | Chronological distribution of SARS-CoV-2 diagnoses by qPCR at IHU Méditerranée Infection institute (A), and of mean (± standard deviation) numbers of mutations in SARS-CoV-2 genomes obtained from patients sampled per month (B). (A) Chronological distribution of SARS-CoV-2 diagnoses is the mean number of diagnoses by qPCR per sliding windows of 7 days and steps of 1 day. (B) Numbers of mutations in SARS-CoV-2 genomes were calculated in reference to the genome of the Wuhan-Hu-1 isolate (GenBank Accession no. NC_045512.2). Whiskers indicate 10–90 percentiles. \*p < 10⁻¹; \*\*p < 10⁻³.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**FIGURE 3** | Chronological distribution of mean (± standard deviation) pairwise genetic distances according to time periods for SARS-CoV-2 genomes obtained from patients SARS-CoV-2-diagnosed at IHU Méditerranée Infection, Marseille. ***p* < 10⁻³.</figcaption></figure>

diagnosed in our institute as of 18 August 2021 (Table 2). SARS-CoV-2 genotyping was performed by genome sequencing in our institute for at least 12% of positive diagnoses of SARS-CoV-2 infection monthly, and for at least one fifth and one third

of positive diagnoses monthly for 14 and 7 of the 18 months from February 2020 to July 2021, respectively (Supplementary Table 2). Genotype was obtained primarily by whole genome next-generation sequencing in 11,387 (47%) cases and by partial genome next-generation sequencing (spike fragment) in 2,621 (11%) cases, and retrospectively or prospectively (since January 2021) by in house variant-specific qPCR assays gradually designed and implemented subsequent to the detection of new variants (Bedotto et al., 2021a,b), in 10,173 (42%) cases. The most frequently detected variants were, in decreasing order, the Alpha (a.k.a. UK) variant, first detected in late December 2020 [Pangolin lineage B.1.1.7; *n* = 9,294 patients (38.4%)]; the Marseille-4 variant, first detected in July 2020 [B.1.160; 6,125 (25.3%)]; the Delta (a.k.a. Indian) variant, first detected in April 2021 [B.1.617.2; 4,113 (17.0%)]; the Marseille-2 variant, first detected in August 2020 [B.1.177; 890 (3.7%)]; and the Beta (a.k.a. South African) variant, first detected in January 2020 [B.1.351; 575 (2.4%)]. Other variants of significant prevalence and/or interest in our geographical area were the Marseille-8 variant [B.1.1.269; 176 (0.7%)]; the Marseille-5 variant [B.1.367; 124 (0.5%)]; the Marseille-1 variant [B.1.416; 123 (0.5%)]; the Marseille-9 variant [B.1.1.241; 117 (0.5%)]; the Marseille-484K.V3 variant [B.1.525; 112 (0.5%)], the Marseille-10 variant [B.1.221; 91 (0.4%)]; the Gamma (a.k.a. Brazilian) variant [P.1; 87 (0.4%)]; and the Marseille-501 variant [A.27; 74 (0.3%)].

## Dynamics of SARS-CoV-2 Variants in the IHU Méditerranée Infection Cohort of SARS-CoV-2-Infected Patients

Based on our definition of variants, we were able to identify several since the summer of 2020, which we would currently call VOC (Box 2). After having defined and named the Marseille variants, we attempted to characterize their dynamic and to analyze their epidemiological features. For some of them, we were able to find their source. During the onset of the SARS-CoV-2 epidemic in late February 2020, we observed two distinct SARS-CoV-2 genotypes that did not make it possible to define variants based on our classification system. They corresponded

**BOX 1** | Definition of virus groups at and below species level.

**Viral species:** Defined according to a polythetic classification as monophyletic groups of viruses whose properties can be distinguished from those of other species by multiple criteria including increasingly genetic criteria (https://talk.ictvonline.org/taxonomy/w/ictv-taxonomy) (Fauquet and Stanley, 2005; Adams et al., 2013).

Below the level of species, there are currently no universally recognized or used strategies and criteria to classify viruses (Fauquet et al., 2008; International Committee on Taxonomy of Viruses (ICTV), 2020).

**Viral clade:** Monophyletic group comprised by a common ancestor and all its descent. These are branches with an independent root that obey a Darwinian evolutionary process.

**Viral isolate:** Viral population obtained by culture, generally linked to a particular virus sample (Fauquet and Stanley, 2005).

**Viral strain:** Virus of a given species with stable and heritable biological, serological, and/or molecular characteristics, and/or links to particular hosts, vectors, or symptoms (Fauquet and Stanley, 2005; Van Regenmortel, 2007).

**Viral quasi-species:** The mutant clouds (or swarms) with closely related genomes generated according to a Darwinian evolutionary process. It includes the accumulation of mutations during the replication of RNA viruses through continuous genetic variation, with possible genetic rearrangement by genetic recombination or genome segment reassortment, followed by competition among the mutants generated and selected in a given environment (Eigen, 1971, 1996; Domingo et al., 1976; Andino and Domingo, 2015).

**Viral mutant:** The result of nucleotide substitutions (non-synonymous or synonymous), deletions or insertions. These nucleotide changes occur randomly, gradually, and spontaneously, according to the Darwinian model (Eigen, 1993; Koonin, 2009; Domingo and Perales, 2019) and possibly at sites critical for viral protein conformation, antigenicity, susceptibility to immune responses, and functionality, including ability to interact with other proteins and enzymatic properties (Rambaut et al., 2004; Harvey et al., 2021).

**Viral variant:** An isolate that differs slightly by its transmission mode or the symptoms it causes from a reference virus (Fauquet and Stanley, 2005; Fauquet et al., 2008) or whose consensus genome sequence differs from that of a reference virus (Fauquet and Stanley, 2005; Van Regenmortel, 2007; Fauquet et al., 2008). Notwithstanding, no definition of the term "variant" is universal or broadly accepted (Fauquet et al., 2008; Kuhn et al., 2013). See Marseille classification of SARS-CoV-2 variants in Box 2.
**BOX 2** | Major SARS-CoV-2 classification systems.

**World Health Organization (WHO) classification (https://www.who.int/en/activities/tracking-SARS-CoV-2-variants/):** Defines SARS-CoV-2 'Variants of Interest' (VOI) and 'Variants of Concern' (VOC). **VOI:** viral genome has mutations with established or suspected phenotypic implications and either has been identified to cause community transmission/multiple COVID-19 cases/clusters or has been detected in multiple countries. **VOC:** a VOI demonstrated through comparative assessment to be associated with (i) increase transmissibility or detrimental change in COVID-19 epidemiology; (ii) increase in virulence or change in clinical disease presentation; or (iii) decrease in effectiveness of public health and social measures or available diagnostics, vaccines, and/or therapeutics. Nomenclature uses letters of the Greek alphabet (currently from Alpha to Theta as of 1 August 2021).

**GISAID classification (https://www.gisaid.org/references/statements-clarifications/clade-and-lineage-nomenclature-aids-in-genomic-epidemiology-of-active-hcov-19-viruses/)** (Alm et al., 2020): Defines high-level phylogenetic clades informed by the statistical distribution of genome distances in phylogenetic clusters followed by merging of smaller lineages into major clades (eight currently) based on shared marker mutations. Nomenclature uses actual letters of marker mutations (S, L, V, G, GH, GR, GV, GRY).

**Nextstrain classification (https://nextstrain.org/blog/2021-01-06-updated-SARS-CoV-2-clade-naming)** (Hadfield et al., 2018; Aksamentov et al., 2021): Defines major clades that reach ≥ 20% global frequency for ≥ two months, ≥ 30% regional frequency for ≥ two months, or correspond to a 'variant of concern' (VOC). Nomenclature uses a year-letter naming system.

**Pangolin (for Phylogenetic Assignment of Named Global Outbreak Lineages) classification (https://cov-lineages.org/pangolin.html)** (Rambaut et al., 2020): Based on a maximum likelihood phylogenetic reconstruction with the following criteria: (i) Phylogenetic evidence of emergence from an ancestral lineage into another geographically distinct population; (ii) ≥ one shared nucleotide difference(s) from the ancestral lineage; (iii) ≥ five genomes; (iv) ≥ one shared nucleotide change(s) among genomes from the lineage; and (v) a bootstrap value > 70% for the lineage node. Regarding the nomenclature: a letter (A for the Wuhan-Hu-1 isolate) is given first, followed by numbers for a maximum of three sublevels after which new descendant lineages are given a new letter.

**Marseille classification of variants:** Defines groups of viral genomes carrying specific sets of ≥ five mutations along the whole genome differentiating them from any other viral genomes and obtained from ≥30 different patients. Nomenclature uses our city name, "Marseille," followed by a number according to the chronology of their description (e.g., Marseille-1, Marseille-2, ...) or to one of their hallmark amino acid substitutions (e.g., Marseille-501).

to Nextstrain clades 19A and 20B that briefly co-circulated in our geographical area during the first period of the SARS-CoV-2 epidemic (Figures 5-7, Supplementary Figure 4, and Supplementary Table 3; IHU Méditerranée Infection, 2021a,b). The first clade was diagnosed from 15 patients in March 2020. The second, which corresponds to viruses harboring the D614G substitution in the spike, was initially retrieved from patients who stayed in northwest Italy and at the French/Italian border, and quickly established in the Marseille area as a majority clade together with clades 20A and 20C until May 2020.

We named the first variant we defined Marseille-1 (corresponds to Pangolin clade B.1.416), which emerged on week 27 of 2020 and predominated in July 2020 before rapidly disappearing at the end of August 2020 (Colson et al., 2021c; Figures 5-7, Supplementary Figure 4, and Supplementary Table 3; IHU Méditerranée Infection, 2021a,b). Its onset

coincided with the resumption on 1 July 2020 of maritime connections between Marseille and the Maghreb. As 16 of the first 20 diagnosed cases (80%) were detected in people who had traveled by boat or on passenger boats circulating between the Maghreb and Marseille, we assumed that this variant was brought to Marseille from the Maghreb. Furthermore, genomic analyses showed that it originated from sub-Saharan Africa, mostly Senegal (Colson et al., 2021c).

The Marseille-4 variant (B.1.160) was the most prevalent SARS-CoV-2 variant in our geographical area in 2020 and was only outdone by the Alpha (UK, B.1.1.7) variant in 2021 (Fournier et al., 2021). It replaced the Marseille-1 variant in August 2020 and faded away in April 2021. It was first detected in south-west France and in Marseille. Its emergence coincided with the re-increase in France of SARS-CoV-2 incidence starting with a department in northwest France (Mayenne) close to the Eure-et-Loir department where a mink farm was found to be massively infected with SARS-CoV-2 (Fournier et al., 2021; Santé Publique France, 2021). The only genome from this farm, obtained from a mink sampled in mid-November 2020 but released in March 2021, is that of a Marseille-4 variant (GISAID Accession no. EPI_ISL_1392906). In addition, the co-occurrence of 13 hallmark mutations in the Marseille-4 variant without any available genome showing a gradual appearance of these mutations suggested that genetic evolution had been overlooked. These data led us to conclude, based on genomic, geographical and temporal coincidences, that the source of the Marseille-4 epidemic was the mink. The Marseille-4 variant was progressively replaced from January 2021 by the epidemic linked to the Alpha (UK, B.1.1.7) variant. Meanwhile, it circulated in 76 countries and accounted for 31,299 genomes in the GISAID database as of 30 June 2021, mostly obtained from European countries [in 29,906 cases (96%)] (Supplementary Figure 5).

During July and August 2020, two additional variants, Marseille-2 and Marseille 5, emerged in our geographical area, which became the second and third most prevalent variants from week 11 of 2020 after the Marseille-4 variant. The Marseille-2 variant (B.1.177) most likely emerged in Spain, where it accounted for two-thirds of released genomes and it spread to our geographical area and through Europe during the summer of 2020 (Hodcroft et al., 2021). The Marseille-5 variant (B.1.367) disappeared in October 2020 and its source and origin are unknown. The Marseille-8 (B.1.1.269) and Marseille-9 (B.1.1.241) variants were first detected in August 2020. They accounted overall for approximately 1.2% of the cases of SARS-CoV-2 infections in our geographical area, but nonetheless circulated until March and January 2021, respectively. The source and origin of these variants are unknown.

The Marseille-501 variant (A.27) emerged in Marseille in January 2021 as the first variant of lineage 19B harboring amino acid substitution N501Y in the spike protein (Colson et al., 2021b). A majority of the genomes primarily available from the GISAID database were linked to Mayotte, which is part of the Comoros archipelago located in the Indian Ocean from where originated some of the patients we diagnosed to be infected with this variant. The Marseille-484K.V3-20A variant (B.1.525) emerged in April 2021 in Marseille and was the majority variant
<figure><img src="image_8.png" /><figcaption>**FIGURE 4** | Phylogenetic tree of 10,773 genomic sequences of SARS-CoV-2 obtained from patients SARS-CoV-2-diagnosed at IHU Méditerranée Infection, Marseille, France. Phylogeny reconstruction was performed using the nextstrain/ncov tool (https://github.com/nextstrain/ncov) then visualized with Auspice (https://docs.nextstrain.org/projects/auspice/en/stable/). The genome of the original Wuhan-Hu-1 coronavirus isolate (GenBank accession no. NC_045512.2) was added as outgroup. Major (most prevalent) variants are labeled. Marseille variants or WHO clades, and Pangolin, or Nextstrain clades are indicated. X-axis shows the number of mutations compared to the genome of the Wuhan-Hu-1 isolate (GenBank accession no. NC_045512.2).</figcaption></figure>

apart from Beta (South African) and Gamma (Brazilian) variants to harbor amino acid substitution E484K in the spike protein. This variant was first described in Nigeria on 11 December 2020, but we have not documented its importation from this country.

The Alpha variant (B.1.1.7) emerged in December 2020 in France and in Marseille and faded away in June 2021 after having been responsible for the majority of SARS-CoV-2 infections from February 2021. The first cases were part of a family cluster during the Christmas holidays after some family members came from England by train. We also had smaller epidemics with the Beta (South African, B.1.351) and the Gamma (Brazilian) (P.1) variants from January 2021. Regarding the Gamma variant, we first detected it among people returning from the Comoros who had traveled through Ethiopia then Tunisia to Marseille. Between February and April 2021, we also detected three cases of a rare spike E484K-harboring SARS-CoV-2, which we named Marseille-484K.V1 (A.23) (Colson et al., 2021a). The three genomes were phylogenetically closest to a genome from the Republic of Guinea (a.k.a. Guinea-Conakry), where one of the patients had traveled.

Most recently, we have mostly detected infections with the Delta (B.617.2) variant, in 4,113 patients. The two first patients diagnosed with this variant in late April 2021 were an Indian sailor re-embarking in Marseille after the resumption of maritime travel (La Scola et al., 2021), and a physiotherapist who had traveled to the Republic of the Maldives, in Southeast Asia, and provided care upon his return to France in an accommodation facility for dependent elderly people, being involved in a cluster among residents and health care workers. Finally, in mid-July 2021, we detected a variant classified in the Pangolin lineage

B.1.621 (a.k.a. Mu), that was first identified in Colombia in January 2021 (Laiton-Donato et al., 2021). Information was available for one of the six patients infected with this variant, a 10-year-old Colombian girl who developed symptoms and was diagnosed in our institute 4 days after she arrived in France.

Overall, we identified several SARS-CoV-2 variants that were responsible as early as the summer of 2020 for distinct epidemics causing several waves of SARS-CoV-2 incidence (Figures 5-7; IHU Méditerranée Infection, 2021a,b). The epidemics caused by these different variants each exhibited a bell-shaped curve, and they occurred successively or concomitantly, therefore contributing to the total case load. We documented that new SARS-CoV-2 variants were introduced in our geographical area by boat in one case [Marseille-1 (B.1.416)], by train in one case [Alpha (B.1.1.7)], and by plane in five cases [Beta (B.1.351), Marseille-501 (A.27), Marseille-484K.V1 (A.23), Delta (B.1.617.2), and B.1.621 variants]. In addition, the importation of Marseille-484K.V3 (B.1.525) from Nigeria through air travel is suspected. Importation from abroad through travel therefore accounted for a substantial source of variants and, subsequently, of cases. The other source identified for the emergence of a new variant is the mink epidemics in Northern France, with the case of Marseille-4 (B.1.160) (Fournier et al., 2021).

DISCUSSION

In this study, we wanted to clarify the classification and naming of SARS-CoV-2 genotypes we implemented in our
<table><caption>TABLE 1 | Nomenclature and description of Marseille (IHU Méditerranée Infection) variants of SARS-CoV-2.</caption><thead><tr><th>Marseille variant (IHU Méditerranée Infection)</th><th>Number of genomes in the IHU Méditerranée Infection dataset</th><th>Nextstrain clade</th><th>Pangolin lineage</th><th>WHO classification</th><th>Mutations compared to the genome of the Wuhan-Hu-1 isolate (GenBank accession no. NC_045512.2)</th></tr></thead><tbody><tr><td colspan="6"><strong>Defined</strong></td></tr><tr><td>Marseille-1</td><td>88</td><td>20A</td><td>B.1.416</td><td>-</td><td>C28833T, G1181T, C1625T, C25886T, G28198T, G28851T</td></tr><tr><td>Marseille-2</td><td>617</td><td>20E (EU1)</td><td>B.1.177</td><td>-</td><td>T445C, C6286T, C22227T, C26801G, C28932T, G29645T, G21255C</td></tr><tr><td>Marseille-3</td><td>43</td><td>20A</td><td>B.1</td><td>-</td><td>C1912T, G5210A, C17470T, C21191T, A23148G, T27125C, C28854T</td></tr><tr><td>Marseille-4</td><td>2,442</td><td>20A.EU2</td><td>B.1.160</td><td>-</td><td>C18877T, C26735T, G5629T, G9526T, G11497T, G13993T, G15766T, G22992A, C25710T, T26876C, G28975C, G29399A, C4543T, A16889G, G17019T</td></tr><tr><td>Marseille-5</td><td>135</td><td>20C</td><td>B.1.367</td><td>-</td><td>C10582T, C28830A, C3099T, G4960T, C4965T, C6070T, C7303T, C7564T, C9246T, C10279T, C10301A, C10525T, G10688T, G11851T, C14230A, G21800T, G27632T, C27804T, G29402T, G29779T</td></tr><tr><td>Marseille-8</td><td>188</td><td>20B</td><td>B.1.1.269</td><td>-</td><td>C5055T, G11851T, A12755G, G24812T, C26895T</td></tr><tr><td>Marseille-9</td><td>114</td><td>20B</td><td>B.1.1.241</td><td>-</td><td>A11782G, T21570G, C21575T, T25473C, C28253T</td></tr><tr><td>Marseille-10</td><td>86</td><td>20A</td><td>B.1.221</td><td>-</td><td>C21855T, A25505G, G25906C, G25996T, C28651T, C28869T, C3602T, C6941T</td></tr><tr><td>Marseille-12</td><td>55</td><td>20A</td><td>-</td><td>-</td><td>G29734C, C17104T, T7767C, C8047T, C22879A</td></tr><tr><td>Marseille-501</td><td>47</td><td>19B</td><td>A.27</td><td>-</td><td>G28878A, G29742A, A361G, C1122T, C2509T, A9204G, A11217G, C16466T, A18366G, A20262G, T22917G, A23063T, C23520T, C23525T, G23948T, G25218T, T25541C, TTAATCCAG26159T, C27247T, TA27385T, AGATTTCC28247A, A28273T</td></tr><tr><td>Marseille-484K.V3</td><td>84</td><td>21D</td><td>B.1.525</td><td>Eta</td><td>C1498T, A1807G, G2659A, C6285T, T8593C, GTCTGGTTTT11287G, C14407T, C18171T, A20724G, A21717G, C21762T, ATACATGT21764AT, TTTATTA21990TTTA, G23012A, G23593C, T24224C, C24748T, C26305T, T26767C, GTTT27204G, TCTG28277T, C28308G, A28699G, C28887T, G29543T, TTTCACCGAGGCCACGCGGAGTACGATCGAGTGTACAGTGAAACAAT GCTAGG29727TGG</td></tr><tr><td>Marseille-452R-19B</td><td>51</td><td>19B</td><td>A.21</td><td>-</td><td>G28878A, G29742A, T1885A, C5907T, C10138T, G11417T, C11824T, C18129T, G18181A, T21982C, G22132T, T22917G, G23402A, G25687T, G26062T, C28253T, C29686T</td></tr><tr><td colspan="6"><strong>Potential</strong></td></tr><tr><td>Marseille-6</td><td>11</td><td>20A</td><td>B.1</td><td>-</td><td>C9430T, A15477T, C18395T, A20622T, G20623T, A20624T, C23730T, A26319G, C28854T, G29044A</td></tr><tr><td>Marseille-7</td><td>17</td><td>20A</td><td>B.1.416.1</td><td>-</td><td>C28833T, C2706T, C25731T, G27463C</td></tr><tr><td>Marseille-11</td><td>19</td><td>20A</td><td>B.1.221</td><td>-</td><td>C21855T, A25505G, G25906C, G25996T, C28651T, C28869T, T445C, G13514T, G25785T, C25810T, G29511T, C29708T</td></tr><tr><td>Marseille-484K.V1-20B</td><td>1</td><td>20B</td><td>R.1</td><td>-</td><td>C14340T, G17551A, C18877T, A19167G, C19274A, T19839C, G22017T, G23012A, G23868T, T26604C, C27213T, TAAAA28270TAAA, C28833T, G29527T</td></tr><tr><td>Marseille-484K.V2-20B</td><td>9</td><td>20B</td><td>B.1.1.318</td><td>-</td><td>C3961T, C9072T, C9891T, C10116T, GTCTGGTTTT11287G, G20578A, C21846T, TTTATTA21990TTTA, G23012A, T23287C, C23604A, G23948C, C24382T, C25276A, T26767C, AAACGAACATGAAATTT27886AT, G28209T, A28271G, GCTA28895G, C29769T</td></tr><tr><td>Marseille-681H-19B</td><td>6</td><td>19B</td><td>A.23</td><td>-</td><td>G28878A, G29742A, G11230T, G22661T, G23401T, G28167A, G569A, C5178T, C7389T, C10029T, A15498T, C21646T, T23030C, C23604A, G25323T, C25603G, G25785T, T26767C</td></tr><tr><td>Marseille-681R-19B</td><td>7</td><td>19B</td><td>A.23.1</td><td>-</td><td>C9430T, A15477T, C18395T, A20622T, G20623T, A20624T, C23730T, A26319G, C28854T, G29044A</td></tr><tr><td>Marseille-4B</td><td>7</td><td>-</td><td>-</td><td>-</td><td>C18877T, C26735T, G5629T, G9526T, G11497T, G13993T, G15766T, G22992A, C25710T, T26876C, G28975C, G29399A, C1059T, C17634T, T19584C, G19684T, T25081C, G28378A, G29553T</td></tr><tr><td>Marseille-20E.2</td><td>14</td><td>-</td><td>-</td><td>-</td><td>C3646T, C6636T, C2862T, G5629T, A10652G, C13275T, C20844T, C21575T, G22104T, C24866T</td></tr></tbody></table>

institute in order to have a simple identification of SARS-CoV-2 variants and to characterize their mutations, their source, and spread. This is essential to get a better understanding of

SARS-CoV-2 epidemiology. Since the SARS-CoV-2 emergence in France, we have performed in our institute extensive surveillance, identification and monitoring of SARS-CoV-2
<table><caption>**TABLE 2** | Numbers of the main SARS-CoV-2 variants and mutants diagnosed among the IHU Méditerranée Infection cohort of SARS-CoV-2-infected patients as of 18 August 2021.</caption><thead><tr><th rowspan="2">Variant/ mutant name</th><th colspan="3">Genotyping approach</th><th rowspan="2">Total</th><th rowspan="2">% of all genotypes</th></tr><tr><th>Next-generation genome sequencing</th><th>Next-generation spike gene fragment sequencing</th><th>qPCR</th></tr></thead><tbody><tr><td>Alpha (B.1.1.7, a.k.a. UK)</td><td>3,506</td><td>2,220</td><td>3,568</td><td>9,294</td><td>38.44</td></tr><tr><td>Marseille-4 (B.1.160)</td><td>2,564</td><td>45</td><td>3,516</td><td>6,125</td><td>25.33</td></tr><tr><td>Delta (B.1.617.2, a.k.a. Indian)</td><td>1,406</td><td>0</td><td>2,707</td><td>4,113</td><td>17.01</td></tr><tr><td>Marseille-2 (B.1.177)</td><td>679</td><td>6</td><td>205</td><td>890</td><td>3.68</td></tr><tr><td>Beta (B.1.351, a.k.a. South African)</td><td>204</td><td>210</td><td>161</td><td>575</td><td>2.38</td></tr><tr><td>20A/8371T (1st period)</td><td>488</td><td>0</td><td>0</td><td>488</td><td>2.02</td></tr><tr><td>20A (1st period)</td><td>323</td><td>44</td><td>0</td><td>367</td><td>1.52</td></tr><tr><td>20C (1st period)</td><td>310</td><td>0</td><td>0</td><td>310</td><td>1.28</td></tr><tr><td>20B (1st period)</td><td>290</td><td>6</td><td>0</td><td>296</td><td>1.22</td></tr><tr><td>20A/15324T (1st period)</td><td>176</td><td>4</td><td>0</td><td>180</td><td>0.74</td></tr><tr><td>Marseille-8 (B.1.1.269)</td><td>174</td><td>2</td><td>0</td><td>176</td><td>0.73</td></tr><tr><td>Marseille-5 (B.1.367)</td><td>124</td><td>0</td><td>0</td><td>124</td><td>0.51</td></tr><tr><td>Marseille-1 (B.1.416)</td><td>123</td><td>0</td><td>0</td><td>123</td><td>0.51</td></tr><tr><td>20A/25563T (1st period)</td><td>123</td><td>0</td><td>0</td><td>123</td><td>0.51</td></tr><tr><td>Marseille-9 (B.1.1.241)</td><td>117</td><td>0</td><td>0</td><td>117</td><td>0.48</td></tr><tr><td>Marseille-484K.V3 (B.1.525)</td><td>112</td><td>0</td><td>0</td><td>112</td><td>0.46</td></tr><tr><td>20A/20268G (1st period)</td><td>95</td><td>0</td><td>0</td><td>95</td><td>0.39</td></tr><tr><td>Marseille-10 (B.1.221)</td><td>89</td><td>2</td><td>0</td><td>91</td><td>0.38</td></tr><tr><td>20A/2416T (1st period)</td><td>89</td><td>0</td><td>0</td><td>89</td><td>0.37</td></tr><tr><td>Gamma (P.1, a.k.a. Brazilian)</td><td>70</td><td>1</td><td>16</td><td>87</td><td>0.36</td></tr><tr><td>Marseille-501 (A.27)</td><td>26</td><td>48</td><td>0</td><td>74</td><td>0.31</td></tr><tr><td>Marseille-452R-19B (A.21)</td><td>43</td><td>2</td><td>0</td><td>45</td><td>0.19</td></tr><tr><td>Marseille-3 (B.1)</td><td>39</td><td>0</td><td>0</td><td>39</td><td>0.16</td></tr><tr><td>19B (1st period)</td><td>13</td><td>23</td><td>0</td><td>36</td><td>0.15</td></tr><tr><td>Marseille-12</td><td>36</td><td>0</td><td>0</td><td>36</td><td>0.15</td></tr></tbody></table>

(Continued)

<table><caption>**TABLE 2** | (Continued)</caption><thead><tr><th rowspan="2">Variant/ mutant name</th><th colspan="3">Genotyping approach</th><th rowspan="2">Total</th><th rowspan="2">% of all genotypes</th></tr><tr><th>Next-generation genome sequencing</th><th>Next-generation spike gene fragment sequencing</th><th>qPCR</th></tr></thead><tbody><tr><td>20D (B.1.1.1)</td><td>31</td><td>2</td><td>0</td><td>33</td><td>0.14</td></tr><tr><td>20E (B.1.177)</td><td>21</td><td>0</td><td>0</td><td>21</td><td>0.09</td></tr><tr><td>Marseille-484K.V2 (B.1.1.318)</td><td>20</td><td>1</td><td>0</td><td>21</td><td>0.09</td></tr><tr><td>Marseille-11 (B.1.221)</td><td>19</td><td>0</td><td>0</td><td>19</td><td>0.08</td></tr><tr><td>19A (1st period)</td><td>15</td><td>0</td><td>0</td><td>15</td><td>0.06</td></tr><tr><td>B.1.214 (a.k.a. Belgian)</td><td>15</td><td>0</td><td>0</td><td>15</td><td>0.06</td></tr><tr><td>Marseille-7 (B.1416.1)</td><td>15</td><td>0</td><td>0</td><td>15</td><td>0.06</td></tr><tr><td>Marseille-6 (B.1)</td><td>9</td><td>0</td><td>0</td><td>9</td><td>0.04</td></tr><tr><td>Marseille-681R-19B (A.23.1)</td><td>7</td><td>0</td><td>0</td><td>7</td><td>0.03</td></tr><tr><td>Marseille-681H-19B (A.23)</td><td>6</td><td>0</td><td>0</td><td>6</td><td>0.02</td></tr><tr><td>Epsilon (B.1.427/B.1.429)</td><td>6</td><td>0</td><td>0</td><td>6</td><td>0.02</td></tr><tr><td>B.1.621</td><td>6</td><td>0</td><td>0</td><td>6</td><td>0.02</td></tr><tr><td>Marseille-484.V1 (A.23)</td><td>3</td><td>0</td><td>0</td><td>3</td><td>0.01</td></tr></tbody><tfoot><tr><td>Total</td><td>11,392</td><td>2,616</td><td>10,173</td><td>24,181</td><td>100.0</td></tr></tfoot></table>

variants. Currently, the viral genotype has been determined for 44% of the patients diagnosed in our institute and the viral genome was obtained in approximatively one fifth of patients diagnosed with SARS-CoV-2. As early as the summer of 2020, we increased our surveillance of SARS-CoV-2 genotypes which allowed us to decipher several features of the current SARS-CoV-2 pandemic on a local scale. Before the summer of 2020, very few studies had pointed out an increase of SARS-CoV-2 genetic and amino acid diversity. Noticeably, Tomaszewski et al. (2020) had reported new pathways of mutational changes in SARS-CoV-2 proteomes. We used several approaches for SARS-CoV-2 classification as it can be determined based on a limited number of genetic markers (in the case of the spike gene fragment sequenced here) (Colson et al., 2021b) and even a single one (in the case of variant-specific qPCR) (Bedotto et al., 2021a,b). Whole genome analysis provides the most precise classification as it considers multiple genetic markers along the genomes and allows characterizing genetic diversity beyond lineage identification. However, variant-specific qPCR allows rapid classification at low cost and high throughput including in case of low viral load.

First, it appears relatively clear that in countries which have been isolated, either geographically (islands that have closed their
access) or politically (China and Korea), there has been only one major epidemic episode with a bell-shaped incidence curve. In Europe and the United States, where such a policy has not been implemented or has been made impossible by political circumstances, several epidemics have occurred successively or concurrently with different viral variants. Marseille, which is a very particular city due to its geographical location, has for centuries been at the forefront of imports of epidemics including of plague and cholera (Barbieri et al., 2021), and operated in an identical way during the SARS-CoV-2 pandemic. Thus, we observed that the first diagnosed SARS-CoV-2 infections were imported cases involving people who had stayed in northwest Italy and at the French/Italian border. This is congruent with findings from a study that inferred the introduction of SARS-CoV-2 from Hubei province into Italy and France before early March 2020 and also estimated that viruses that spread during the epidemic onset in Europe had a common ancestor in Italy (Nadeau et al., 2021). This preceded a first epidemic phase during which international borders were closed. Their reopening in early summer 2020 was followed by a short epidemic. This was linked to a SARS-CoV-2 variant that we named Marseille-1 (B.1.416), which was imported from North Africa by boat and originated from sub-Saharan Africa, and was probably not very transmissible, since this epidemic ended 2 months after its emergence and did not spread beyond the Marseille region. Another variant, which we named Marseille-2 (B.1.177), was introduced to our geographical area during the summer of 2020 from Spain, where it appears to have originated (Hodcroft et al., 2021). Later, in 2021, the Alpha (B.1.1.7, or UK) variant was first imported in Marseille by French people returning from the UK. In addition, in the Marseille population which includes lots of people of Comorian origin, we observed the onset of small epidemics with the Beta (B.1.351, or South African) variant among people traveling by plane from the Comoros, and with the Marseille-501 (A.27) variant that involved people of Comorian origin including from Mayotte (Colson et al., 2021b). Similarly, the reopening of sea cruises in the spring of 2020 was associated with an importation of the Delta (B.1.617.2, or Indian) variant by an Indian sailor who had traveled from India to Marseille (La Scola et al., 2021), and a second case identified of this variant was in a traveler returning from the Republic of the Maldives located in the Indian Ocean. Two other SARS-CoV-2 variants were imported in 2021 from the Republic of Guinea and from Colombia. We therefore observed that the importation from abroad through travel accounted for a substantial source of variants and subsequently of infections, highlighting the role of globalization in the spread of SARS-CoV-2 pandemic. Marseille is a city with many immigrant adults and children, and its geographical location and harbor make it very open to travelers. This explains the importation of infectious agents of various origins, which has been recognized in Marseille since antiquity (Barbieri et al., 2021). The considerable role of travel in the introduction and subsequent spread of new SARS-CoV-2 variants on the national or continental scales has also been reported in a recent phylogeographical study that analyzed genomic, epidemiological, and mobility data (including ours) collected from ten European countries between January and October 2020

(Lemey et al., 2021). This study found that the intensity of international travel predicted the spread of SARS-CoV-2 with the introduction of new viral variants, and it estimated that more than half of the lineages that were spreading in late summer 2020 had been newly introduced to countries since mid-June 2020. Such importations into countries of SARS-CoV-2 variants and their impact on the pattern of incidence of SARS-CoV-2 cases have also been reported in other studies (McLaughlin et al., 2021). This sharply confirms what we described as early as at the beginning of September 2020 (Colson et al., 2020b) with the emergence of seven new SARS-CoV-2 variants including one [Marseille-1 (B.1.416)] for which we had already found a source (Colson et al., 2021c). Thus, we were the first to consider and highlight that the re-increase of SARS-CoV-2 incidence observed during the summer of 2020 in Marseille and in France may have resulted from the importation of new variants rather than from a rebound of previous viral genotypes (Colson et al., 2020b, 2021c). This aligns with the fact that countries for which cross-border travel was the most limited geographically and/or politically by closing the borders experienced a single epidemic peak, whereas countries with high levels of international travel experienced several epidemic peaks and longer periods with a significant incidence of infections. Hence, these findings indicate that while local control measures such as lockdowns and curfews may limit viral spread, they are inefficient to avoid the introduction of new viral genotypes. The role of local lockdowns in controlling the viral spread has been controversial. This measure has been found to be irrelevant by some researchers (Melnick and Ioannidis, 2020). In Japan, lockdown was not associated with a reduction of SARS-CoV-2 incidence, whereas travel restrictions were (Hideki and Skidmore, 2021).

Second, our findings indicate that the role of epizootics in dense animal herds has been overlooked in the generation of variants, their transmission to humans, and their contribution to the SARS-CoV-2 pandemic. Mink was revealed as a major source of SARS-CoV-2 diversity and a key element in the expansion of the human pandemic (Fenollar et al., 2021; Oude Munnink et al., 2021). Data have accumulated indicating that the current pandemic, which is primarily a zoonosis, continued being a zoonosis with the massive involvement of minks, while other animals, including domestic animals such as cats and dogs, have also been affected (Fenollar et al., 2021). It has been demonstrated that minks were a source of new variants that fuelled the pandemic in Denmark (Oude Munnink et al., 2021), and SARS-CoV-2 spike protein gene variants harboring amino acid substitutions N501T and G142D were also found to circulate in both humans and minks in the United States (Cai and Cai, 2021). In Marseille, we accumulated indications that the Marseille-4 variant originated from mink. This hypothesis took us a long time to confirm due to the very late release, in March 2021, of the viral sequence obtained from a mink farm in Eure-et-Loir in mid-November 2020. Interestingly, the detection of two epidemics, in Marseille and retrospectively in Mayenne (Northwest France), appeared to be unrelated as the absence of SARS-CoV-2 genotyping initially carried out in Mayenne and Eure-et-Loir prevented the demonstration that they were due to the same variant
<figure><img src="image_9.png" /></figure>

and that virus transmission could have occurred from farm
minks located in a neighboring department. It is worthy to
note that the Marseille-4 variant was responsible for a large
epidemic in France and also to a lesser extent in other European
countries and in United Kingdom and was probably considerably
overlooked as the global concern about SARS-CoV-2 variants
only emerged several months after we described it. A recent
study of global and regional adaptive evolution of SARS-
CoV-2 delineated a divergent viral clade with a pattern of
signature mutations consistent with the Marseille-4 variant. This
clade was reported to have first appeared in Europe in April
2020 then to have considerably expanded during summer 2020
(Rochman et al., 2021).

These findings show that SARS-CoV-2 genomic surveillance
and the definition of variants are valuable for a better assessment
of SARS-CoV-2 epidemiology. In addition, they allow correlating
viral genotypes with the clinical severity and outcome of
infections. Indeed, several studies have reported different clinical
patterns of infections according to the SARS-CoV-2 variant (Dao
et al., 2021; Domingo et al., 2021). We previously compared
the clinical symptoms and outcomes of patients infected with
different variants we identified. Compared with viruses of the
Nextstrain clade 20A that circulated between February and May
2020 in Marseille, patients infected with the Marseille-1 (B.1.416)
variant were less likely to report dyspnea and rhinitis and to be
hospitalized (Dao et al., 2021). In addition, patients infected with

the Marseille-4 variant (B.1.160) were more likely to exhibit fever and to be hospitalized than those infected with the Marseille-1 variant, and they were also more likely to exhibit fever than those infected with a virus of Nextstrain clade 20A that circulated between February and May 2020 in Marseille. Also, patients infected with N501Y-harboring SARS-CoV-2 variants [namely, Alpha (B.1.1.7), Beta (B.1.351), and Gamma (P.1) variants] were less likely to be hospitalized than those infected with the Marseille-4 variant and with viruses of Nextstrain clade 20A that circulated between February and May 2020 in Marseille (Dao et al., 2021). SARS-CoV-2 genomic surveillance can also provide insights for phenotypic, structural and functional investigations of mutations and variants, as has been already broadly performed (Dong et al., 2021; Fantini et al., 2021; Gupta et al., 2021; Jaafar et al., 2021; Rahman et al., 2021). Particularly, SARS-CoV-2 mutation patterns in the spike protein may affect binding of neutralizing antibodies or T-cell recognition.

Furthermore, a key issue from the present study is highlighting
that the fluctuations in incidence of SARS-CoV-2 infections that
we observed resulted from the succession or combination of
distinct epidemics caused by different viral variants rather than
from the rebound of persisting lineages, which we primarily
described in September 2020 and was confirmed at larger
scale (Lemey et al., 2021). Since September 2020, these genetic
evolutionary and epidemiological patterns continued to be
observed and the surveillance of SARS-CoV-2 variants, the role of
<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>
travel and globalization in their spread, and the turnover of major VOC has come to the frontline, particularly with the Alpha, Beta, Gamma, and Delta variants (Del Rio et al., 2021; Harvey et al., 2021).

Overall, previous findings have revealed that epidemic waves mostly observed in developed Western countries are, in fact, distinct epidemics due to different SARS-CoV-2 variants. Therefore, as the major identified sources for these SARS-CoV-2 variants is cross border travel and animal reservoirs, efficient border controls and surveillance of mustelid farms are critical in order to avoid the introduction of new SARS-CoV-2 variants and the occurrence of new epidemics.

# MATERIALS AND METHODS

## Epidemic Curves for SARS-CoV-2-Associated Deaths in Various Countries and Clustering of the Dynamic of Epidemics at Country Scale

We used the cumulated daily mortality from the COVID-19 data repository operated by Johns Hopkins University (Dong et al., 2020) available at https://github.com/CSSEGISandData/COVID-19 (accessed on 25 June 2021), which covers 221 countries and locations from 22 April 2020 to 24 June 2021. Data mining of time series is a complex scientific field where the analysis process must be tailored to the characteristics of the series and the clustering objectives (Liao, 2005; Fu, 2011). Our objective was to cluster the dynamic of epidemics at country scale as reflected by COVID-19 mortality, taking in account the time of the different series peaks but not country mortality incidences. A threshold of 500 deaths was used for countries to analyze their SARS-CoV-2-associated death curves. We transformed the cumulative death counts into smoothed (7 days centered moving average) daily death rates. The resulting series were then clustered using the Pearson distance after normalization as L2-normed lock-step distance measure (Mori et al., 2016) and the ascendant hierarchical classification algorithm with Ward aggregation. For determining the resulting clusters, we used a dynamic cluster detection method (Langfelder et al., 2008) with a minimum cluster size set at five countries. We then mapped the geographic distribution of cluster members. The data mining process was done using R v4.0.2⁸ with the packages vegan, forecast, dynamictreecut, and rnaturalearth.

## Study Period and Clinical Samples

SARS-CoV-2 genotyping was performed from nasopharyngeal samples tested between 27 February 2020 and 18 August 2021 (18 months) at the IHU Méditerranée Infection institute<fnref n="9" />. Specimens with a cycle threshold value (Ct) of qPCR lower than 20 were selected as a in priority for whole genome sequencing of SARS-CoV-2 genomes, and those with a Ct between 20 and 30 were included secondarily to ensure a more comprehensive coverage of the study period. Partial sequencing of the SARS-CoV-2 spike gene generating a 1,854 nucleotide-long sequence

was performed in 2020 and, mainly, between January and April 2021, as previously described (Colson et al., 2021b). For respiratory specimens with Ct values > 30 or those with Ct values < 30 but from which genome sequences were not obtained, we identified those harboring specific Marseille genotypes using qPCR targeting variant-specific regions, as previously described for the Marseille-1 variant (Colson et al., 2021c) and the Marseille-4 variant (Bedotto et al., 2021b; Fournier et al., 2021). Additional in house and commercial variant-specific qPCR assays were used for other variants. The study was approved by the ethics committee of the University Hospital Institute Méditerranée Infection under No. 2020-016-03.

## SARS-CoV-2 Genome Sequencing

Viral RNA was extracted from 200 µL of nasopharyngeal swab fluid using the EZ1 Virus Mini Kit v2.0 on an EZ1 Advanced XL instrument (Qiagen, Courtaboeuf, France) or the KingFisher Flex system (Thermo Fisher Scientific, Waltham, MA, United States), following the manufacturer's instructions. SARS-CoV-2 genome sequences were obtained by next-generation sequencing with various procedures since February 2020 until August 2021: (i) with the Illumina technology using the Nextera XT paired end strategy on MiSeq instruments (Illumina Inc., San Diego, CA, United States) since February 2020, as previously described (Colson et al., 2021c); (ii) with the Illumina COVIDSeq protocol on a NovaSeq 6000 instrument (Illumina Inc.) since April 2021; or (iii) with Oxford Nanopore technology (ONT) on MinION or GridION instruments (Oxford Nanopore Technologies Ltd., Oxford, United Kingdom), as previously described (Colson et al., 2021c). Next-generation sequencing with ONT was performed without or with (since March 2021) synthesized cDNA amplification using a multiplex PCR protocol with ARTIC nCoV-2019 V3 Panel primers purchased from Integrated DNA technologies (IDT, Coralville, IA, United States) according to the ARTIC procedure<fnref n="10" />. After its extraction, viral RNA was reverse-transcribed using SuperScript IV (ThermoFisher Scientific) prior to cDNA second strand synthesis with Klenow Fragment DNA polymerase (New England Biolabs, Beverly, MA, United States) when performing NGS on the Illumina MiSeq instrument (Illumina Inc.) (Colson et al., 2021c), LunaScript RT SuperMix kit (New England Biolabs) when performing NGS with the ONT, or according to the COVIDSeq protocol (Illumina Inc.) following the manufacturer's recommendations. Generated cDNA was purified using Agencourt AMPure XP beads (Beckman Coulter, Villepinte, France) and quantified using Qubit 2.0 fluorometer (Invitrogen, Carlsbad, CA, United States); fragment sizes were analyzed on an Agilent Fragment analyser 5200 (Agilent Inc., Palo Alto, CA, United States). Next-generation sequencing performed on the NovaSeq 6000 instrument followed the Illumina COVIDSeq protocol (Illumina Inc.), which included first strand cDNA synthesis from extracted viral RNA; cDNA amplification with two COVIDSeq primer pools; fragmentation and tagging of PCR amplicons with adapter sequences; clean up; PCR amplification (seven cycles) of tagmented amplicons; pool, clean up, quantification and normalization of libraries; and

<footnote n="8">⁸https://cran.r-project.org/</footnote>

<footnote n="9">⁹https://www.mediterranee-infection.com/</footnote>

<footnote n="10">¹⁰https://artic.network/</footnote>
library sequencing on a NovaSeq 6000 sequencing system SP flow cell (Illumina Inc.).

## Assembly and Analyses of Genome Sequences

Genome sequences were assembled by mapping on the SARS-CoV-2 genome GenBank accession no. NC_045512.2 (Wuhan-Hu-1 isolate) with the CLC Genomics workbench v.7 (with the following thresholds: 0.8 for coverage and 0.9 for similarity)<fnref n="11" /> as previously described (Colson et al., 2021c), or Minimap2¹² (Li, 2018). Samtools<fnref n="13" /> was used for soft clipping of Artic primers (see footnote 10) and to remove sequence duplicates (Li et al., 2009). Consensus genomes were generated using the CLC Genomics workbench v.7, and Sam2consensus<fnref n="14" /> through a first in-house script written in Python language<fnref n="15" />. Detection of mutations was performed using the Nextclade tool (see footnote 4), and with freebayes<fnref n="16" /> (Garrison and Marth, 2012) using a mapping quality score of 20 and results filtering by the Python script based on majoritary nucleotide frequencies ≥ 70% and nucleotide depths ≥ 10 [when sequence reads were generated on the NovaSeq Illumina instrument (Illumina Inc.)] or ≥5 [when sequences reads were generated on a MiSeq Illumina instrument (Illumina Inc.)]. SARS-CoV-2 genotyping was performed by a second in-house script written in Python language (see footnote 15) by comparing mutation patterns with those of our database of SARS-CoV-2 variants. Nextstrain clades and Pangolin lineages provided in the present study were determined with the Nextclade web application (see footnote 4) (Hadfield et al., 2018; Aksamentov et al., 2021) and the Pangolin web application (see footnote 5) (Rambaut et al., 2020), respectively. Sequences described in the present study have been deposited on the GISAID sequence database (see footnote 3) (Alm et al., 2020) and can be retrieved online using the GISAID online search tool with "Marseille" as keyword, or selecting sequence names containing "IHU" or "MEPHI." In addition, they have been deposited on the IHU Marseille Infection website: https://www.mediterranee-infection.com/tout-sur-le-coronavirus/sequencage-genomique-sars-cov-2/.

## Phylogeny Reconstructions and Definition and Naming of SARS-CoV-2 Variants

Phylogeny reconstructions based on the SARS-CoV-2 genomes obtained in our laboratory were performed for genome sequences larger than 24,000 nucleotides obtained with the Illumina technology using the nextstrain/ncov tool<fnref n="17" /> then visualized with Auspice<fnref n="18" /> or FigTree v1.4.4¹⁹; the limitation of sampling per date

of collection of respiratory specimens was deactivated. Within reconstructed trees, we identified and counted all mutations between nodes and between nodes and leaves. This allowed us to define the set of mutations for the most recent ancestor of the clusters corresponding to each of the SARS-CoV-2 variants.

## Numbers of Mutations in SARS-CoV-2 Genomes and Pairwise Genetic Distances for SARS-CoV-2 Genomes Over Time

To compare the numbers of mutations in SARS-CoV-2 genomes and the pairwise genetic distances for these genomes over time, the numbers of nucleotide changes in the SARS-CoV-2 genomes and of amino acid changes in the spike protein relatively to the genome of the Wuhan-Hu-1 isolate (GenBank accession no. NC_045512.2) were obtained using the Nextclade tool (see footnote 4) (Hadfield et al., 2018; Aksamentov et al., 2021) for SARS-CoV-2 genomes used in the phylogenetic analyses. Pairwise nucleotide distances for SARS-CoV-2 genomes were computed using the MEGA7 software v10.2.5²⁰. The numbers of nucleotide and amino acid changes and the pairwise nucleotide distances were plotted using the GraphPad software v5.01²¹. The mean numbers of mutations in viral genomes and the mean pairwise genetic distances between genomes over time since SARS-CoV-2 emergence in February 2020 were compared using R 4.0.2 (see footnote 8) and GraphPad software v5.01 (see footnote 21). A *p* < 0.05 was considered statistically significant.

## Real-Time Reverse Transcription-PCR Specific of Marseille Variants

For specimens with Ct values > 30 or those with Ct values < 30 but from which genome sequences were not obtained, we identified those harboring specific Marseille genotypes using qPCR targeting variant-specific regions. We previously described qPCR specific of the Marseille-1 variant (Colson et al., 2021c) and the Marseille-4 variant (Fournier et al., 2021). We used other in-house variant-specific qPCR systems for which the sequences of primers and probes are as follows, provided in 5'-3' orientation. For the Beta (B.1.351) variant, qPCR targeted the envelope gene and used forward primer C_SA_3_MBF: TGAATTGCAGACACCTTTTGA, reverse primer C_SA_3_MBR: CAACCCTTGGTTGAATAGTCTTG, and probe C_SA_3_MBP: TGACATCTTCAATGGGGAATGT. For the Gamma (P.1) variant, qPCR targeted the nucleocapsid gene and used forward primer C_Bra_2_MBF: GTCAAGCCTCTTCTCGTTCCT, reverse primer C_Bra_2_MBR: AAAGCAAGAGCAGCATCACC, and probe C_Bra_2_MBP: GCAGCTCTAAACGAACTTCTCCTG. For the Delta (B.1.617.2) variant, two qPCR were used that targeted the envelope gene: the first one used forward primer C_IND_1_MBF:

<footnote n="11">¹¹https://digitalinsights.qiagen.com/</footnote>

<footnote n="12">¹²https://github.com/lh3/minimap2</footnote>

<footnote n="13">¹³https://www.htslib.org/</footnote>

<footnote n="14">¹⁴https://github.com/vbsreenu/Sam2Consensus</footnote>

<footnote n="15">¹⁵https://www.python.org/</footnote>

<footnote n="16">¹⁶https://github.com/freebayes/freebayes</footnote>

<footnote n="17">¹⁷https://github.com/nextstrain/ncov</footnote>

<footnote n="18">¹⁸https://docs.nextstrain.org/projects/auspice/en/stable/</footnote>

<footnote n="19">¹⁹http://tree.bio.ed.ac.uk/software/figtree/</footnote>

<footnote n="20">²⁰https://www.megasoftware.net/</footnote>

<footnote n="21">²¹https://www graphpad com/</footnote>
AATCTTGATTCTAAGGTTGGTGGT, reverse primer C_IND_1_MBR: TGCTACCGGCCTGATAGATT, and probe C_IND_1_MBP: TTACCGGTATAGATTGTTTAGGAAGTCT; the second one used forward primer P681_MBF: TGACATACCCATTGGTCAG, reverse primer P681_MBR: GGCAATGATGGATTGACTAGC, and probe P681_MBP: ATTCTCATCGGCGGCA. Finally, for the Marseille-2 (B.1.177) variant, qPCR targeted the nsp16 gene and used forward primer C2_2_MBF: CATGGTGGACAGCCTTTGTT, reverse primer C2_2_MBR: CATCTATTTGTTCGCGTGGTT, and probe C2_2_MBP: AATGCCTCATCATCTGAAGCAT. PCR conditions were the same than those described previously in Bedotto et al. (2021b) and Fournier et al. (2021). Finally, we used the Applied Biosystems TaqPath COVID-19 kit (Thermo Fisher Scientific, Waltham, United States) then the TaqMan SARS-CoV-2 mutation assay H69-V70 (Thermo Fisher Scientific) for Alpha (B.1.1.7) variant detection, and TaqMan SARS-CoV-2 mutation assays L452R and P681R (Thermo Fisher Scientific) to detect the Delta (B.1.617.2) variant.

## DATA AVAILABILITY STATEMENT

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material.

## ETHICS STATEMENT

This study received approval of the Ethical Committee of the Mediterranee Infection institute under reference N°2020-016-3. Access to the patients' biological and registry data issued from the hospital information system was approved by the data protection committee of Assistance Publique-Hopitaux de Marseille (APHM) and was recorded in the European General Data Protection Regulation registry under number RGPD/APHM 2019-73; the latter authorization allowed us to access the patients' data retrospectively, according to the French law on "Règlement Général sur la Protection des Données (RGPD)".

## REFERENCES

* Adams, M. J., Lefkowitz, E. J., King, A. M., and Carstens, E. B. (2013). Recently agreed changes to the international code of virus classification and nomenclature. *Arch. Virol.* 158, 2633–2639. doi: 10.1007/s00705-013-1749-9

* Aksamentov, I., Roemer, C., Hodcroft, E. B., and Neher, R. A. (2021). Nextclade: clade assignment, mutation calling and quality control for viral genomes. *J. Open Source Softw.* 6:3773. doi: 10.5281/zenodo.5607694

* Alm, E., Broberg, E. K., Connor, T., Hodcroft, E. B., Komissarov, A. B., Maurer-Stroh, S., et al. (2020). Geographical and temporal distribution of SARS-CoV-2 clades in the WHO European Region, January to June 2020. *Euro. Surveill.* 25:2001410. doi: 10.2807/1560-7917.ES.2020.25.32.2001410

* Andino, R., and Domingo, E. (2015). Viral quasispecies. *Virology* 47, 46–51. doi: 10.1016/j.virol.2015.03.022

## AUTHOR CONTRIBUTIONS

DR, PC, and P-EF: conceptualization and original draft preparation. DR, PC, P-EF, and AL: methodology. PC, P-EF, HC, JD, AG-G, LH, CA, LB, MBed, EP, CG, MBey, EB, PD, HT-D, PG, J-CL, MM, PB, PP, FF, MD, BL, and AL: data and investigation. DR: supervision. All authors: review and editing of the manuscript, read and agreed to the published version of the manuscript.

## FUNDING

This work was supported by the French Government under the "Investments for the Future" programme managed by the National Agency for Research (ANR), Méditerranée-Infection 10-IAHU-03 and was also supported by Région Provence Alpes Côte d'Azur and European funding PRIMMI (Plateformes de Recherche et d'Innovation Mutualisées Méditerranée Infection), FEDER PA 0000320 PRIMMI.

## ACKNOWLEDGMENTS

Marseille public hospitals are part of the French Emergen consortium (https://www.santepubliquefrance.fr/dossiers/coronavirus-covid-19/consortium-emergen). We are thankful to the IHU Méditerranée Infection COVID-19 Task Force that notably includes Sophie Amrane, Camille Aubry, Sofiane Bakour, Karim Bendamardji, Cyril Berenger, Nadim Cassir, Anne Darmon, Claire Decoster, Catherine Dhiver, Barbara Doudier, Sophie Edouard, Véronique Filosa, Séverine Guitton, Marie Hocquart, Morgane Mailhe, Florine Pelisson, Coralie Porcheto, Isabelle Ravaux, Piseth Seng, Laurence Thomas, Raphael Tola, Christelle Tomei, and Catherine Triquet. This manuscript has been edited for English language by a native English speaker.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.786233/full#supplementary-material

* Barbieri, R., Colson, P., Raoult, D., and Drancourt, M. (2021). Two-millennia fighting against port-imported epidemics, Marseille. *IHU Preprints* [Preprint] doi: 10.35088/84a4-me41

* Bedotto, M., Fournier, P. E., Houhamdi, L., Colson, P., and Raoult, D. (2021a). Implementation of an in-house real-time reverse transcription-PCR assay to detect the emerging SARS-CoV-2 N501Y variants. *J. Clin. Virol.* 140:104868. doi: 10.1016/j.jcv.2021.104868

* Bedotto, M., Fournier, P. E., Houhamdi, L., Levasseur, A., Delerce, J., Pinault, L., et al. (2021b). Implementation of an in-house real-time reverse transcription-PCR assay for the rapid detection of the SARS-CoV-2 Marseille-4 variant. *J. Clin. Virol.* 139:104814. doi: 10.1016/j.jcv.2021.104814

* Cai, H. Y., and Cai, A. (2021). SARS-CoV2 spike protein gene variants with N501T and G142D mutation-dominated infections in mink in the United States. *J. Vet. Diagn. Invest.* 33, 939–942. doi: 10.1177/10406387211023481
Chinazzi, M., Davis, J. T., Ajelli, M., Gioannini, C., Litvinova, M., Merler, S., et al. (2020). The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak. *Science* 368, 395–400. doi: 10.1126/science. aba9757

Colson, P., Levasseur, A., Delerce, J., Chaudet, H., Bossi, V., Ben Khedher, M., et al. (2020b). Dramatic increase in the SARS-CoV-2 mutation rate and low mortality rate during the second epidemic in summer in Marseille. *IHU Preprints* [Preprint] doi: 10.35088/68c3-ew82

Colson, P., Lagier, J. C., Baudoin, J. P., Bou Khalil, J., La Scola, B., and Raoult, D. (2020a). Ultrarapid diagnosis, microscope imaging, genome sequencing, and culture isolation of SARS-CoV-2. *Eur. J. Clin. Microbiol. Infect. Dis.* 39, 1601–1603. doi: 10.1007/s10096-020-03869-w

Colson, P., Levasseur, A., Gautret, P., Fenollar, F., Thuan, H. V., Delerce, J., et al. (2021c). Introduction into the Marseille geographical area of a mild SARS-CoV-2 variant originating from sub-Saharan Africa: an investigational study. *Travel. Med. Infect. Dis.* 40:101980. doi: 10.1016/j.tmaid.2021.101980

Colson, P., Levasseur, A., Delerce, J., Pinault, L., Dudouet, P., Devaux, C., et al. (2021b). Spreading of a new SARS-CoV-2 N501Y spike variant in a new lineage. *Clin. Microbiol. Infect.* 27, 1352.e1–1352.e5. doi: 10.1016/j.cmi.2021.05.006

Colson, P., Fantini, J., Yahi, N., Delerce, J., Levasseur, A., Fournier, P.-E., et al. (2021a). Limited spread of a rare spike E484K-harboring SARS-CoV-2 in Marseille, France. *IHU Preprints* [Preprints] doi: 10.35088/2ngw-sa32

Cucinotta, D., and Vanelli, M. (2020). WHO Declares COVID-19 a Pandemic. *Acta Biomed.* 91, 157–160.

Dao, T. L., Hoang, V. T., Nguyen, N. N., Delerce, J., Chaudet, H., Levasseur, A., et al. (2021). Clinical outcomes in COVID-19 patients infected with different SARS-CoV-2 variants in Marseille, France. *Clin. Microbiol. Infect.* 10, 1516.e1–1516.e6.

Del Rio, C., Malani, P. N., and Omer, S. B. (2021). Confronting the delta variant of SARS-CoV-2, summer 2021. *JAMA* 326, 1001–1002. doi: 10.1001/jama.2021.14811

Domingo, E., Flavell, R. A., and Weissmann, C. (1976). In vitro site-directed mutagenesis: generation and properties of an infectious extracristonic mutant of bacteriophage Qbeta. *Gene* 1, 3–25. doi: 10.1016/0378-1119(76)90003-2

Domingo, E., and Perales, C. (2019). Viral quasispecies. *PLoS Genet.* 15:e1008271. doi: 10.1371/journal.pgen.1008271

Domingo, P., Pomar, V., Mur, I., Castellvi, I., Corominas, H., and de Benito, N. (2021). Not all COVID-19 pandemic waves are alike. *Clin. Microbiol. Infect.* 27:1040. doi: 10.1016/j.cmi.2021.04.005

Dong, E., Du, H., and Gardner, L. (2020). An interactive web-based dashboard to track COVID-19 in real time. *Lancet Infect. Dis.* 20, 533–534. doi: 10.1016/S1473-3099(20)30120-1

Dong, J., Zost, S. J., Greaney, A. J., Starr, T. N., Dingens, A. S., Chen, E. C., et al. (2021). Genetic and structural basis for SARS-CoV-2 variant neutralization by a two-antibody cocktail. *Nat. Microbiol.* 6, 1233–1244. doi: 10.1038/s41564-021-00972-2

Eigen, M. (1971). Molecular self-organization and the early stages of evolution. *Q. Rev. Biophys.* 4, 149–212. doi: 10.1017/s0033583500000627

Eigen, M. (1993). Viral quasispecies. *Sci. Am.* 269, 42–49. doi: 10.1038/scientificamerican0793-42

Eigen, M. (1996). On the nature of virus quasispecies. *Trends Microbiol.* 4, 216–218. doi: 10.1016/0966-842X(96)20011-3

Fantini, J., Yahi, N., Azzaz, F., and Chahinian, H. (2021). Structural dynamics of SARS-CoV-2 variants: a health monitoring strategy for anticipating Covid-19 outbreaks. *J. Infect.* 83, 197–206. doi: 10.1016/j.jinf.2021.06.001

Fauquet, C. M., Briddon, R. W., Brown, J. K., Moriones, E., Stanley, J., Zerbini, M., et al. (2008). Geminivirus strain demarcation and nomenclature. *Arch. Virol.* 153, 783–821. doi: 10.1007/s00705-008-0037-6

Fauquet, C. M., and Stanley, J. (2005). Revising the way we conceive and name viruses below the species level: a review of geminivirus taxonomy calls for new standardized isolate descriptors. *Arch. Virol.* 150, 2151–2179. doi: 10.1007/s00705-005-0583-0

Fenollar, F., Mediannikov, O., Maurin, M., Devaux, C., Colson, P., Levasseur, A., et al. (2021). Mink, SARS-CoV-2, and the human-animal interface. *Front. Microbiol.* 12:663815. doi: 10.3389/fmicb.2021.663815

Fournier, P. E., Colson, P., Levasseur, A., Devaux, C. A., Gautret, P., Bedotto, M., et al. (2021). Emergence and outcomes of the SARS-CoV-2 'Marseille-4' variant. *Int. J. Infect. Dis.* 106, 228–236. doi: 10.1016/j.ijid.2021.03.068

Fu, T. (2011). A review on time series data mining. *Eng. Appl. Artif. Intell.* 24, 164–181. doi: 10.1016/j.engappai.2010.09.007

Garrison, E., and Marth, G. (2012). *Haplotype-Based Variant Detection From Short-Read Sequencing*. Available online at: https://arxiv.org/abs/1207.3907 (accessed September 9, 2012).

Gupta, D., Sharma, P., Singh, M., Kumar, M., Ethayathulla, A. S., Kaur, P., et al. (2021). Structural and functional insights into the spike protein mutations of emerging SARS-CoV-2 variants. *Cell. Mol. Life Sci.* 78, 7967–7989. doi: 10.1007/s00018-021-04008-0

Hadfield, J., Megill, C., Bell, S. M., Huddleston, J., Potter, B., Callender, C., et al. (2018). Nextstrain: real-time tracking of pathogen evolution. *Bioinformatics* 34, 4121–4123. doi: 10.1093/bioinformatics/bty407

Harvey, W. T., Carabelli, A. M., Jackson, B., Gupta, R. K., Thomson, E. C., Harrison, E. M., et al. (2021). SARS-CoV-2 variants, spike mutations and immune escape. *Nat. Rev. Microbiol.* 19, 409–424. doi: 10.1038/s41579-021-00573-0

Hideki, T., and Skidmore, M. (2021). “A Cross-Country Analysis of the Determinants of Covid-19 Fatalities,” in *Paper Presented at CESifo Working Papers* 9028, (Munich).

Hodcroft, E. B., Zuber, M., Nadeau, S., Vaughan, T. G., Crawford, K. H. D., Althaus, C. L., et al. (2021). Spread of a SARS-CoV-2 variant through Europe in the summer of 2020. *Nature* 595, 707–712. doi: 10.1038/s41586-021-03677-y

IHU Méditerranée Infection (2021a). *Evolution of the Weekly Incidence of SARS-CoV-2 Mutants and Variants Extrapolated to the Total Number of Cases, Based on Their Proportions of Genotyped Cases*. Marseille: HU Méditerranée Infection. Available online at: https://doi.org/10.35081/ffd4-1y77 (accessed September 9, 2021).

IHU Méditerranée Infection (2021b). *Introduction of SARS-CoV-2 Mutants and Variants in the Marseille Geographical Area Through Travel From Abroad*. Marseille: HU Méditerranée Infection. Available online at: https://doi.org/10.35081/nz2g-f980 (accessed September 9, 2021).

International Committee on Taxonomy of Viruses (ICTV) (2020). *Virus Taxonomy: The ICTV Report on Virus Classification and Taxon Nomenclature. The Online 10th Report of the International Committee on Taxonomy of Viruses*. Available online at: https://talkICTVonline.org/ictv1/f/general_ictv_discussions-20/3938/ictv-online-10th-report (accessed September 9, 2021).

Jaafar, R., Boschi, C., Aherfi, S., Bancod, A., Le Bideau, M., Edouard, S., et al. (2021). High individual heterogeneity of neutralizing activities against the original strain and nine different variants of SARS-CoV-2. *Viruses* 13:2177. doi: 10.3390/v13112177

Koonin, E. V. (2009). Darwinian evolution in the light of genomics. *Nucleic Acids Res.* 37, 1011–1034. doi: 10.1093/nar/gkp089

Kuhn, J. H., Bao, Y., Bavari, S., Becker, S., Bradfute, S., Brister, J. R., et al. (2013). Virus nomenclature below the species level: a standardized nomenclature for laboratory animal-adapted strains and variants of viruses assigned to the family Filoviridae. *Arch. Virol.* 158, 1425–1432. doi: 10.1007/s00705-012-1594-2

La Scola, B., Lavrard, P., Fournier, P. E., Colson, P., Lacoste, A., and Raoult, D. (2021). SARS-CoV-2 variant from India to Marseille: the still active role of ports in the introduction of epidemics. *Travel. Med. Infect. Dis.* 42:102085. doi: 10.1016/j.tmaid.2021.102085

Langfelder, P., Zhang, B., and Horvath, S. (2008). Defining clusters from a hierarchical cluster tree: the dynamic tree Cut package for R. *Bioinformatics* 24, 719–720. doi: 10.1093/bioinformatics/btm563

Laiton-Donato, K., Franco-Munoz, C., Alvarez-Diaz, D. A., Ruiz-Moreno, H. A., Usme-Ciro, J. A., Prada, D. A., et al. (2021). Characterization of the emerging B.1.621 variant of interest of SARS-CoV-2. *Infect. Genet. Evol.* 95:105038. doi: 10.1016/j.megid.2021.105038

Lemey, P., Ruktanonchai, N., Hong, S. L., Colizza, V., Poletto, C., Van den Broeck, F., et al. (2021). Untangling introductions and persistence in COVID-19 resurgence in Europe. *Nature* 595, 713–717. doi: 10.1038/s41586-021-03754-2

Levasseur, A., Delerce, J., Caputo, A., Brechard, L., Colson, P., Lagier, J. C., et al. (2020). Genomic diversity and evolution of coronavirus (SARS-CoV-2) in France from 309 COVID-19-infected patients. *bioRxiv* [Preprint] doi: 10.1101/2020.09.04.282616

Li, H. (2018). Minimap2: pairwise alignment for nucleotide sequences. *Bioinformatics* 34, 3094–3100. doi: 10.1093/bioinformatics/bty191

Li, H., Handsaker, B., Wysocker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The sequence alignment/map format and SAMtools. *Bioinformatics* 25, 2078–2079. doi: 10.1093/bioinformatics/btp352
Liao, T. W. (2005). Clustering of time series data—a survey. *Pattern Recogn.* 38, 1857–1874. doi: 10.1016/j.patcog.2005.01.025

McLaughlin, A., Montoya, V., Miller, R. L., Mordecai, G. J., Worobey, M., Poon, A., et al. (2021). Early and ongoing importations of SARS-CoV-2 in Canada. *medRxiv* [Preprint] doi: 10.1101/2021.04.09.21255131

Melnick, E. R., and Ioannidis, J. P. A. (2020). Should governments continue lockdown to slow the spread of covid-19? *BMJ* 369:m1924. doi: 10.1136/bmj.m1924

Mori, U., Mendiburu, A., and Lozano, J. A. (2016). Similarity measure selection for clustering time series databases. *IEEE. Trans. Knowl. Data.* 28, 181–194. doi: 10.1109/TKDE.2015.2462369

Nadeau, S. A., Vaughan, T. G., Scire, J., Huisman, J. S., and Stadler, T. (2021). The origin and early spread of SARS-CoV-2 in Europe. *Proc. Natl. Acad. Sci. U.S.A.* 118:e2012008118. doi: 10.1073/pnas.201208118

Oude Munnink, B. B., Sikkema, R. S., Nieuwenhuijse, D. F., Molenaar, R. J., Munger, E., Molenkamp, R., et al. (2021). Transmission of SARS-CoV-2 on mink farms between humans and mink and back to humans. *Science* 371, 172–177. doi: 10.1126/science.abe5901

Rahman, M. S., Islam, M. R., Alam, A. S. M. R. U., Islam, I., Hoque, M. N., Akter, S., et al. (2021). Evolutionary dynamics of SARS-CoV-2 nucleocapsid protein and its consequences. *J. Med. Virol.* 93, 2177–2195. doi: 10.1002/jmv.26626

Rambaut, A., Holmes, E. C., O'Toole, Å, Hill, V., McCrone, J. T., Ruis, C., et al. (2020). A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology. *Nat. Microbiol.* 5, 1403–1407. doi: 10.1038/s41564-020-0770-5

Rambaut, A., Posada, D., Crandall, K. A., and Holmes, E. C. (2004). The causes and consequences of HIV evolution. *Nat. Rev. Genet.* 5, 52–61. doi: 10.1038/nrg1246

Rochman, N. D., Wolf, Y. I., Faure, G., Mutz, P., Zhang, F., and Koonin, E. V. (2021). Ongoing global and regional adaptive evolution of SARS-CoV-2. *Proc. Natl. Acad. Sci. U.S.A.* 118:e2104241118. doi: 10.1073/pnas.2104241118

Santé Publique France (2021). *Incidence Rate of Sars-CoV-2 Infections by French Department Per Sliding Week*. Saint-Maurice: Santé Publique France. Available online at: https://geodes.santepubliquefrance.fr/ (accessed September 9, 2021).

Sender, R., Bar-On, Y. M., Gleizer, S., Bernshtein, B., Flamholz, A., Phillips, R., et al. (2021). The total number and mass of SARS-CoV-2 virions. *Proc. Natl. Acad. Sci. U.S.A.* 118:e2024815118. doi: 10.1073/pnas.2024815118

Sokal, R. R. (1974). Classification: purposes, principles, progress, prospects. *Science* 185, 1115–1123. doi: 10.1126/science.185.4157.1115

Tomaszewski, T., DeVries, R. S., Dong, M., Bhatia, G., Norsworthy, M. D., Zheng, X., et al. (2020). New pathways of mutational change in SARS-CoV-2 proteomes involve regions of intrinsic disorder important for virus replication and release. *Evol. Bioinform. Online* 16:1176934320965149. doi: 10.1177/1176934320965149

van Dorp, D. L., Richard, D., Tan, C. C. S., Shaw, L. P., Acman, M., and Balloux, F. (2020). No evidence for increased transmissibility from recurrent mutations in SARS-CoV-2. *Nat. Commun.* 11:5986. doi: 10.1038/s41467-020-19818-2

Van Regenmortel, M. H. (2007). Virus species and virus identification: past and current controversies. *Infect. Genet. Evol.* 7, 133–144. doi: 10.1016/j.meegid.2006.04.002

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2022 Colson, Fournier, Chaudet, Delerce, Giraud-Gatineau, Houhamdi, Andrieu, Brechard, Bedotto, Prudent, Gazin, Beye, Burel, Dudouet, Tissot-Dupont, Gautret, Lagier, Million, Brouqui, Parola, Fenollar, Drancourt, La Scola, Levasseur and Raoult. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.