Université Paris Cité and Sorbonne Paris Nord, Inserm, IAME, Paris, France; Centre for Interdisciplinary Research in Biology (CIRB), Collège de France, CNRS Inserm, PSL Research University, Paris, France
Aude Jary
Sorbonne Université, Inserm, iPLESP, Paris, France
Valentin Leducq
Sorbonne Université, Inserm, iPLESP, Paris, France
Valentine Marie Ferré
Université Paris Cité and Sorbonne Paris Nord, Inserm, IAME, Paris, France; Service de Virologie, Hôpital Bichat Claude Bernard, AP-HP, Paris, France
Anna Maria Franco Yusti
Université Paris Cité and Sorbonne Paris Nord, Inserm, IAME, Paris, France; Service de Virologie, Hôpital Bichat Claude Bernard, AP-HP, Paris, France
Léna Daniel
Université Paris Cité and Sorbonne Paris Nord, Inserm, IAME, Paris, France; Service de Virologie, Hôpital Bichat Claude Bernard, AP-HP, Paris, France
Charlotte Charpentier
Université Paris Cité and Sorbonne Paris Nord, Inserm, IAME, Paris, France; Service de Virologie, Hôpital Bichat Claude Bernard, AP-HP, Paris, France
Samuel Lebourgeois
Université Paris Cité and Sorbonne Paris Nord, Inserm, IAME, Paris, France; Service de Virologie, Hôpital Bichat Claude Bernard, AP-HP, Paris, France
Karen Zafilaza
Sorbonne Université, Inserm, iPLESP, Paris, France; Service de Virologie, Hôpital de la Pitié-Salpêtrière, AP-HP, Paris, France
Vincent Calvez
Sorbonne Université, Inserm, iPLESP, Paris, France; Service de Virologie, Hôpital de la Pitié-Salpêtrière, AP-HP, Paris, France
Diane Descamps
Université Paris Cité and Sorbonne Paris Nord, Inserm, IAME, Paris, France; Service de Virologie, Hôpital Bichat Claude Bernard, AP-HP, Paris, France
Anne-Geneviève Marcelin
Sorbonne Université, Inserm, iPLESP, Paris, France; Service de Virologie, Hôpital de la Pitié-Salpêtrière, AP-HP, Paris, France
Benoit Visseaux
Université Paris Cité and Sorbonne Paris Nord, Inserm, IAME, Paris, France; Service de Virologie, Hôpital Bichat Claude Bernard, AP-HP, Paris, France
Although France was one of the most affected European countries by the COVID-19 pandemic in 2020, the dynamics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) movement within France, but also involving France in Europe and in the world, remain only partially characterized in this timeframe. Here, we analyzed GISAID deposited sequences from January 1 to December 31, 2020 (n = 638,706 sequences at the time of writing). To tackle the challenging number of sequences without the bias of analyzing a single subsample of sequences, we produced 100 subsamples of sequences and related phylogenetic trees from the whole dataset for different geographic scales (worldwide, European countries, and French administrative regions) and time periods (from January 1 to July 25, 2020, and from July 26 to December 31, 2020). We applied a maximum likelihood discrete trait phylogeographic method to date exchange events (i.e., a transition from one location to another one), to estimate the geographic spread of SARS-CoV-2 transmissions and lineages into, from and within France, Europe, and the world. The results unraveled two different patterns of exchange events between the first and second half of 2020. Throughout the year, Europe was systematically associated with most of the intercontinental exchanges. SARS-CoV-2 was mainly introduced into France from North America and Europe (mostly by Italy, Spain, the United Kingdom, Belgium, and Germany) during the first European epidemic wave. During the second wave, exchange events were limited to neighboring countries without strong intercontinental movement, but Russia widely exported the virus into Europe during the summer of 2020. France mostly exported B.1 and B.1.160 lineages, respectively, during the first and second European epidemic waves. At the level of French administrative regions, the Paris area was the main exporter during the first wave. But, for the second epidemic wave, it equally contributed to virus spread with Lyon area, the second most populated urban area after Paris in France. The main circulating lineages were similarly distributed among the French regions. To conclude, by enabling the inclusion of tens of thousands of viral sequences, this original phylodynamic method enabled us to robustly describe SARS-CoV-2 geographic spread through France, Europe, and worldwide in 2020.