Phylodynamics of SARS-CoV-2 in France, Europe and the world in 2020

  1. Romain Coppée  Is a corresponding author
  2. François Blanquart
  3. Aude Jary
  4. Valentin Leducq
  5. Valentine Marie Ferré
  6. Anna Maria Franco Yusti
  7. Léna Daniel
  8. Charlotte Charpentier
  9. Samuel Lebourgeois
  10. Karen Zafilaza
  11. Vincent Calvez
  12. Diane Descamps
  13. Anne-Geneviève Marcelin
  14. Benoit Visseaux
  15. Antoine Bridier-Nahmias  Is a corresponding author
  1. Université Paris Cité and Sorbonne Paris Nord, Inserm, IAME, France
  2. Collège de France, France
  3. Sorbonne Université, Inserm, iPLESP, France
  4. Hôpital Bichat-Claude-Bernard, France

Abstract

Although France was one of the most affected European countries by the COVID-19 pandemic in 2020, the dynamics of 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 1st January to 31th December 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 1st January to 25th July 2020, and from 26th July to 31th December 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, 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.

Data availability

All genome sequences and associated metadata in the dataset are published in GISAID's EpiCoV database. To view the contributors of each individual sequence with details such as accession number, virus name, collection date, originating lab and submitting lab and the list of authors, visit: https://doi.org/10.55876/gis8.230120zd.All the scripts developed for this study were deposited in the following GitHub repository: https://github.com/Rcoppee/PhyloCoV

Article and author information

Author details

  1. Romain Coppée

    Université Paris Cité and Sorbonne Paris Nord, Inserm, IAME, Paris, France
    For correspondence
    romain.coppee@inserm.fr
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3024-5928
  2. François Blanquart

    Centre for Interdisciplinary Research in Biology, Collège de France, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0591-2466
  3. Aude Jary

    Sorbonne Université, Inserm, iPLESP, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
  4. Valentin Leducq

    Sorbonne Université, Inserm, iPLESP, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
  5. Valentine Marie Ferré

    Service de Virologie, Hôpital Bichat-Claude-Bernard, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
  6. Anna Maria Franco Yusti

    Université Paris Cité and Sorbonne Paris Nord, Inserm, IAME, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
  7. Léna Daniel

    Université Paris Cité and Sorbonne Paris Nord, Inserm, IAME, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
  8. Charlotte Charpentier

    Service de Virologie, Hôpital Bichat-Claude-Bernard, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
  9. Samuel Lebourgeois

    Université Paris Cité and Sorbonne Paris Nord, Inserm, IAME, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
  10. Karen Zafilaza

    Sorbonne Université, Inserm, iPLESP, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
  11. Vincent Calvez

    Sorbonne Université, Inserm, iPLESP, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
  12. Diane Descamps

    Service de Virologie, Hôpital Bichat-Claude-Bernard, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
  13. Anne-Geneviève Marcelin

    Sorbonne Université, Inserm, iPLESP, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
  14. Benoit Visseaux

    Service de Virologie, Hôpital Bichat-Claude-Bernard, Paris, France
    Competing interests
    The authors declare that no competing interests exist.
  15. Antoine Bridier-Nahmias

    Université Paris Cité and Sorbonne Paris Nord, Inserm, IAME, Paris, France
    For correspondence
    antoine.bridier-nahmias@inserm.fr
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0376-6840

Funding

No external funding was received for this work.

Copyright

© 2023, Coppée et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 821
    views
  • 117
    downloads
  • 4
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Romain Coppée
  2. François Blanquart
  3. Aude Jary
  4. Valentin Leducq
  5. Valentine Marie Ferré
  6. Anna Maria Franco Yusti
  7. Léna Daniel
  8. Charlotte Charpentier
  9. Samuel Lebourgeois
  10. Karen Zafilaza
  11. Vincent Calvez
  12. Diane Descamps
  13. Anne-Geneviève Marcelin
  14. Benoit Visseaux
  15. Antoine Bridier-Nahmias
(2023)
Phylodynamics of SARS-CoV-2 in France, Europe and the world in 2020
eLife 12:e82538.
https://doi.org/10.7554/eLife.82538

Share this article

https://doi.org/10.7554/eLife.82538

Further reading

    1. Epidemiology and Global Health
    2. Genetics and Genomics
    Tianyu Zhao, Hui Li ... Li Chen
    Research Article

    Alzheimer’s disease (AD) is a complex degenerative disease of the central nervous system, and elucidating its pathogenesis remains challenging. In this study, we used the inverse-variance weighted (IVW) model as the major analysis method to perform hypothesis-free Mendelian randomization (MR) analysis on the data from MRC IEU OpenGWAS (18,097 exposure traits and 16 AD outcome traits), and conducted sensitivity analysis with six models, to assess the robustness of the IVW results, to identify various classes of risk or protective factors for AD, early-onset AD, and late-onset AD. We generated 400,274 data entries in total, among which the major analysis method of the IVW model consists of 73,129 records with 4840 exposure traits, which fall into 10 categories: Disease, Medical laboratory science, Imaging, Anthropometric, Treatment, Molecular trait, Gut microbiota, Past history, Family history, and Lifestyle trait. More importantly, a freely accessed online platform called MRAD (https://gwasmrad.com/mrad/) has been developed using the Shiny package with MR analysis results. Additionally, novel potential AD therapeutic targets (CD33, TBCA, VPS29, GNAI3, PSME1) are identified, among which CD33 was positively associated with the main outcome traits of AD, as well as with both EOAD and LOAD. TBCA and VPS29 were negatively associated with the main outcome traits of AD, as well as with both EOAD and LOAD. GNAI3 and PSME1 were negatively associated with the main outcome traits of AD, as well as with LOAD, but had no significant causal association with EOAD. The findings of our research advance our understanding of the etiology of AD.

    1. Epidemiology and Global Health
    Xiaoning Wang, Jinxiang Zhao ... Dong Liu
    Research Article

    Artificially sweetened beverages containing noncaloric monosaccharides were suggested as healthier alternatives to sugar-sweetened beverages. Nevertheless, the potential detrimental effects of these noncaloric monosaccharides on blood vessel function remain inadequately understood. We have established a zebrafish model that exhibits significant excessive angiogenesis induced by high glucose, resembling the hyperangiogenic characteristics observed in proliferative diabetic retinopathy (PDR). Utilizing this model, we observed that glucose and noncaloric monosaccharides could induce excessive formation of blood vessels, especially intersegmental vessels (ISVs). The excessively branched vessels were observed to be formed by ectopic activation of quiescent endothelial cells (ECs) into tip cells. Single-cell transcriptomic sequencing analysis of the ECs in the embryos exposed to high glucose revealed an augmented ratio of capillary ECs, proliferating ECs, and a series of upregulated proangiogenic genes. Further analysis and experiments validated that reduced foxo1a mediated the excessive angiogenesis induced by monosaccharides via upregulating the expression of marcksl1a. This study has provided new evidence showing the negative effects of noncaloric monosaccharides on the vascular system and the underlying mechanisms.