Next Article in Journal
Implications of COVID-19 Pandemic on the Emergence of Antimicrobial Resistance: Adjusting the Response to Future Outbreaks
Next Article in Special Issue
Self-Reported Real-World Safety and Reactogenicity of COVID-19 Vaccines: A Vaccine Recipient Survey
Previous Article in Journal
Endocannabinoid System as Therapeutic Target of PTSD: A Systematic Review
Previous Article in Special Issue
Animal Coronaviruses Induced Apoptosis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Life
Volume 11
Issue 3
10.3390/life11030219
life-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

SARS-CoV-2 Molecular Transmission Clusters and Containment Measures in Ten European Regions during the First Pandemic Wave

by
Maria Bousali
1,
Aristea Dimadi
1,
Evangelia-Georgia Kostaki
2,
Sotirios Tsiodras
3,
Georgios K. Nikolopoulos
4,
Dionyssios N. Sgouras
5,
Gkikas Magiorkinis
2,
George Papatheodoridis
7,
Vasiliki Pogka
1,5,
Giota Lourida
8,
Aikaterini Argyraki
8,
Emmanouil Angelakis
5,6,
George Sourvinos
9,
Apostolos Beloukas
10,11,*,
Dimitrios Paraskevis
2,* and
Timokratis Karamitros
1,5,*
1
Bioinformatics and Applied Genomics Unit, Department of Microbiology, Hellenic Pasteur Institute, 11521 Athens, Greece
2
Department of Hygiene Epidemiology and Medical Statistics, School of Medicine, National and Kapodistrian University of Athens, 15772 Athens, Greece
3
4th Department of Internal Medicine & Infectious Diseases, School of Medicine, National and Kapodistrian University of Athens, 15772 Athens, Greece
4
Medical School, University of Cyprus, 2029 Nicosia, Cyprus
5
Laboratory of Medical Microbiology, Department of Microbiology, Hellenic Pasteur Institute, 11521 Athens, Greece
6
IRD, APHM, VITROME, IHU-Mediterranean Infections, Aix Marseille University, 13005 Marseille, France
7
Department of Gastroenterology, Medical School of National and Kapodistrian University of Athens, “Laiko” General Hospital of Athens, 11527 Athens, Greece
8
Infectious Diseases Clinic A, Sotiria Chest Diseases Hospital, 11527 Athens, Greece
9
Laboratory of Clinical Virology, School of Medicine, University of Crete, 71500 Heraklion, Greece
10
Department of Biomedical Sciences, University of West Attica, 12243 Athens, Greece
11
Institute of Infection and Global Health, University of Liverpool, Liverpool L69 7BE, UK
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Authors to whom correspondence should be addressed.
Life 2021, 11(3), 219; https://doi.org/10.3390/life11030219
Submission received: 18 January 2021 / Revised: 12 February 2021 / Accepted: 3 March 2021 / Published: 9 March 2021
(This article belongs to the Special Issue Ecology, Evolution and Epidemiology of Coronaviruses)
Browse Figures
Versions Notes

Abstract

:
Background: The spatiotemporal profiling of molecular transmission clusters (MTCs) using viral genomic data can effectively identify transmission networks in order to inform public health actions targeting SARS-CoV-2 spread. Methods: We used whole genome SARS-CoV-2 sequences derived from ten European regions belonging to eight countries to perform phylogenetic and phylodynamic analysis. We developed dedicated bioinformatics pipelines to identify regional MTCs and to assess demographic factors potentially associated with their formation. Results: The total number and the scale of MTCs varied from small household clusters identified in all regions, to a super-spreading event found in Uusimaa-FI. Specific age groups were more likely to belong to MTCs in different regions. The clustered sequences referring to the age groups 50–100 years old (y.o.) were increased in all regions two weeks after the establishment of the lockdown, while those referring to the age group 0–19 y.o. decreased only in those regions where schools’ closure was combined with a lockdown. Conclusions: The spatiotemporal profiling of the SARS-CoV-2 MTCs can be a useful tool to monitor the effectiveness of the interventions and to reveal cryptic transmissions that have not been identified through contact tracing.

1. Introduction

Coronaviruses (CoVs) belong to the Nidovirales order, Coronaviridae family and Coronavirinae subfamily and are the largest known group of viruses. Within the past two decades, two newly emerged coronaviruses, Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) and Middle East Respiratory Syndrome Coronavirus (MERS-CoV), have caused serious respiratory and intestinal infections in humans [1]. In December 2019, a novel Coronavirus was firstly reported in the city of Wuhan, Hubei province in China. This novel coronavirus was further named as “2019-nCoV” by WHO [2] and later as Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) by the International Committee on Taxonomy of Viruses. On 24 January 2020, the first European case was reported in France while four days later, Germany confirmed its first case. After six days, on 30 January 2020, WHO declared the 2019-nCoV outbreak to be a Public Health Emergency of International Concern (PHEIC) under International Health Regulations. As of 10 January 2021, almost 88.4 million infections were recorded worldwide, including 1.9 million deaths (https://www.who.int/publications/m/item/weekly-epidemiological-update---12-January-2021, accessed on 16 October 2020). As soon as, the 11 January 2020, the first whole-genome sequence of SARS-CoV-2 was available and became the baseline for researchers to track SARS-CoV-2, as it spread through the world [3]. Based on phylogenetic studies, SARS-CoV-2 belongs to the same lineage as SARS-CoV and MERS-CoV [4,5,6]. An unprecedented number of full genome sequences have become available thanks to the worldwide effort of scientists and to the GISAID consortium [7].
SARS-CoV-2 exhibits a high potential to undergo human-to-human transmission, while the three major factors that are involved in spreading, are the source of infection, the route of transmission and population susceptibility coupled with viral latency [8,9]. The primary transmission mode is person-to-person contact through respiratory droplets, as well as direct contact with an infected subject or indirect contact, through hand-mediated transfer of the virus from contaminated fomites to the mouth, nose, or eyes. The novel coronavirus SARS-CoV-2 is less deadly, but far more transmissible than MERS-CoV or SARS [10], while it has been found that transmission of SARS-CoV-2 is mostly driven by clusters in close contacts, particularly family clusters, and less so by community transmission (https://www.who.int/docs/default-source/coronaviruse/who-china-joint-mission-on-covid-19-final-report.pdf, accessed on 16 October 2020) although super-spreading events continue to occur in the pandemic.
Based on up to date data, there are seven clades of SARS-CoV-2 in total, denoted as G, GH, GR, L, O, S, and V) [11], while the evolutionary rate has been estimated at approximately 10−3–10−4 substitutions per site per year [12,13,14,15,16], which is broadly in line with those estimated for SARS-CoV [17] and Middle East respiratory syndrome [18], about a third of that estimated for influenza B and in general lower than other RNA viruses [19,20]. Phylogenetic analysis of SARS-CoV-2 data is challenging, not only because of the magnitude of available data, but also because of the way the virus has spread in the population in such a short period [21]. As a result, dedicated bioinformatics pipelines and big data analysis tools are of great value in order to rapidly evaluate the factors that are associated with the spread and the transmissibility of the virus. Phylogenomics and phylodynamics analysis in almost real time after the collection and sequencing of the samples may assist in focusing the prevention efforts after the identification of transmission clusters in the communities.
Transmission clusters are groups of infected individuals who are connected with SARS-COV-2 transmission and potentially represent a subset of a risk network. The identification of molecular transmission clusters (MTCs) provides a tool to identify transmission clusters and risk networks and has been previously used for various pathogens, including HIV, Influenza A, and Mycobacterium tuberculosis [22,23,24,25,26]. A time-space cluster occurs by focusing on the identification of MTCs in a particular geographic area while using time-reversible phylodynamics approaches. MTCs have been used to examine the likely impact of genetic (mutations and subtypes of viruses) [27,28,29], demographic and clinical [22,25,26] factors on regional phylogenetic clustering. Identification of MTCs for HIV has resulted in the characterization of the risk factors that are associated with the spread of the virus in different countries [22,30,31]. As for SARS-CoV-2, an analysis of 3184 sequences from Japan with complete metadata, revealed that the primary source of clusters are healthcare facilities such as hospitals and care facilities such as nursing homes [32].
Outbreak containment through testing, case isolation, contact tracing, and quarantine is often the first line of defense against a novel emerging infectious disease [33,34]. During the first pandemic wave of SARS-CoV-2, efforts to contain the outbreak of the virus failed and in order to flatten the epidemic curve supplemental containment measures were taken in almost all countries worldwide [35]. These measures rely on the capacity to control viral transmission from person-to-person and their prioritization can be determined by understanding the virus’ transmission patterns [36]. In Europe, the containment measures taken included large-scale physical distancing measures and movement restrictions—stay at home orders often referred as lockdowns, non-essential shop closure, national events stop, school, nurseries, kindergartens and educational facilities closure, and national and international movement and flights restrictions [37,38,39]. Containment measures may have different efficacy in diverse age groups [40], and the groups for which the measures are more effective may vary across populations [41]. Moreover, demographic, socioeconomic, climatic and cultural factors have been linked to the distribution of COVID-19 cases across the world and they have a potential impact on the effectiveness of the containment measures, too [42].
In this study, we developed dedicated bioinformatics pipelines to conduct a massive identification and spatiotemporal characterization of SARS-CoV-2 MTCs in 10 European geographical regions. We also examine the likely impact of demographic and epidemiological factors, as well as the effect of the containment measures taken, on the regional phylogenetic clustering of the pandemic.

2. Materials and Methods

2.1. Raw Data Pre-Processing and Filtering

All of the available SARS-CoV-2 sequencing data and metadata derived from European samples were downloaded from GISAID [7] on 16 October 2020. From the 82,572 raw sequences originated from the European continent, 5352 sequences were removed from the dataset as they were either incomplete (<29,000 bp) or of low genomic coverage (sequences containing >5% Ns), resulting in 77,220 sequences originated from 1376 different European geographical regions, in total. In order to ensure a maximum genetic diversity, a threshold was set in order to keep only those geographical regions with more than 80 sequences available. Geographical regions were further filtered based on the existence of more than three different SARS-CoV-2 clades spread in the population as well as on the continuity of samples for a minimum period of four weeks after the establishment or after the end of a containment measure. The geographical regions that met all of the above-mentioned criteria were Munich (DE), La Rioja (ES), Navarra (ES), Madrid (ES), Liege (BE), Vienna (AT), Lombardy (IT), Reykjavik (IS), Uusimaa (FI), and Saint Petersburg (RU).

2.2. Reconstruction of Phylogenetic Trees and Phylodynamics Analysis

For each geographical region a separate phylogenetic and phylodynamics analysis was undertaken in order to perform spatiotemporal characterization of MTCs. Firstly, each geographical region’s data and metadata were joined with world’s formatted data and metadata and a FASTA formatted file was generated (sequences.fasta) as well as a TSV file with the corresponding metadata (metadata.tsv). Those files were used as input in Nextstrain’s “augur” pipeline [43] after excluding duplicates. The pipeline involves sequence alignment with MAFFT [44] using NC_045512.2 (also referred to as SARS-CoV-2 isolate Wuhan-Hu-1 GenBank ID: MN908947.3) as the reference sequence and a subsequent subsampling step of sequences that are not referring to the geographical region of interest, based on genetic (alignment), time (collection date) and space (region of exposure) criteria. The NC_045512.2 sequence, was also used during the phylogeny reconstruction along with “Wuhan/WH01/2019” sequence (MT291826.1). The sub-sampled sequences were filtered/processed in parallel with the study sequences. The phylogenetic tree reconstruction was performed using IQ-TREE [45] and the general time-reversible substitution model with gamma rate heterogeneity among sites (GTR+G substitution model), as suggested by “ModelFinder” (implemented in IQ-TREE) [46], and as described in previous studies on SARS-CoV-2 phylogeny [47,48,49]. Re-rooting, resolving of polytomies, inferring of internal node dates, labeling of internal nodes, pruning of sequences and maximum-likelihood phylodynamics analysis was performed with TreeTime [50] using the default parameters that are pre-optimized for SARS-CoV-2 by Nextstrain.

2.3. Identification of Molecular Transmission Clusters (MTCs)

In order to identify MTCs, a two-step approach was followed as previously described by Paraskevis et al. [22]. Given the rapid spread in the population and the relatively slow mutation rate of SARS-CoV-2, which results in low overall genetic variability [51,52], the identification of MTCs from the generated phylogenetic trees was performed using two different criteria: clusters with maximum genetic distance ≤0.005 [16,53] as well as Shimodaira–Hasegawa support (SH-Test) ≥0.75 [54] (phylogenetic confidence criterion) and clusters consisting of geographical region’s specific sequences at a proportion of >65% compared to the total number of sequences within the cluster (geographic criterion). The filtering of MTCs in the trees was performed using in-house developed scripts written in R programming language v4.0.3 (https://www.R-project.org/, accessed on 16 October 2020), utilizing the “ape” [55] and the “phangorn” [56] libraries for phylogenetics and evolutionary analysis, the tidyverse R package for data handling and the “ggtree” library for the visualization of the phylogenies [57] as well as Bash sh-compatible command language.

2.4. Statistical Analysis

All variables were examined as categorical. The continuous variables of age and collection date that were retrieved from GISAID, were transformed into categorical. More specifically, we created five age groups based on social activity (0–19, 20–34, 35–49, 50–64, 65+ years old) as described by Pollan et al. [58]. As for the collection date, we created 32 classes that are referring to weeks, starting from 2020-02-20 when the first sample in the dataset was collected. For simple comparisons of distributions at different levels of categorical variables we performed pairwise chi-square tests, as well as Fisher’s exact test, while the statistical significance threshold was set at 5% (p-value < 0.05). All calculations were performed using R programming language v4.0.3 (https://www.R-project.org/, accessed on 16 October 2020).

3. Results

3.1. Phylogenetic Analysis and Characterization of Molecular Transmission Clusters

For each geographical region a phylogenetic tree was reconstructed using the filtered sequence data of the geographical region and the sub-sampled sequences through the world (Figure 1). Reykjavik, Madrid and Liege (Figure 1g,i,j) account for more than 500 geographical region’s specific sequences while Munich, Vienna, and Navarra (Figure 1a–c) account for less than 200 sequences. After the subsampling, however, all trees were generated with more than 2950 sequences (Table 1). The collection dates in our dataset ranged as the first sequences were reported in Lombardy on 2020-02-20 while the last ones were from Liege on 2020-09-25. Madrid and Lombardy were found to have the most MTCs (20 and 14 respectively), while the other geographical regions had less than 10 MTCs. At the same time, the percentage of clustered sequences was higher in Uusimaa and La Rioja (88.9% and 67.5%, respectively), while Saint Petersburg and Liege were found to have high percentage of sequences not belonging to MTCs (unclustered). Large community clusters and were observed in all geographical regions with the exception of Saint Petersburg. Especially in Uusimaa, a super-spreading event was found (>20 geographical region’s specific sequences [59]), possibly occurred as an MTC consisting of 156 sequences. Smaller clusters were found in all geographical regions. Moreover, patterns based on the clade of the clustered sequences were observed. All geographical regions had great proportion of clustered sequences belonging to the GR clade, except for the Spanish geographical regions of Madrid, La Rioja as well as Uusimaa for which the GR clade was present only in unclustered sequences (25 in Madrid, 3 in La Rioja, and 6 in Uusimaa). In general, geographical regions with fewer clustered sequences showed greater diversity in the distribution of clades, while the combination of distributions of O and V clades were rare and were observed only in Reykjavik.

3.2. Demographic Characteristics and Molecular Transmission Clusters

We investigated the hypothesis that cases belonging to particular age groups might be more important contributors in the formation of MTCs, as their social activity profiles differ substantially. Interestingly, as presented in Table 2, specific age groups were more likely to belong to MTCs in different regions. In detail, in Reykjavik, cases referring to the age group 0–19 were more likely to belong to MTCs (p ≤ 0.05). In Liege, cases referring to the age group 20–34 were more likely to belong to MTCs (p ≤ 0.05) while in Saint Petersburg cases referring to the age group 35–49 were found more likely to belong to MTCs (p < 0.01). Although no specific age groups were directly associated with molecular clustering across all regions analyzed, it is of great interest that these three particular regions presented the lower levels of total clustering, with only 4.8%, 6.7%, and 20.2% of total sequences in MTCs for Saint Petersburg, Liege and Reykjavik, respectively (Figure 2, outer five-color circles). Simultaneously, the age groups of the sampled sequences were evenly distributed in these regions (Figure 2. inner five-color circles), indicating no sampling bias. These findings suggest that in settings with lower levels of clustered dispersal, younger (<50 y.o.) and consequently socially more active individuals drive the formation of MTCs. As expected, gender was not found to be associated with the formation of MTCs in none of the geographical regions studied.

3.3. Containment Measures Taken and Molecular Transmission Clusters

In order to investigate the effect of the containment measures taken on the individual population groups for each geographical region we analyzed the dynamics of both the clustered (in MTCs) and the unclustered sequences on the time scale, with regards to each of the five age groups (0–19, 20–34, 35–49, 50–64, 65+) (Figure 3). In order to examine the effects of each measure in the clustering of each age group, we compared the number of clustered sequences belonging to particular age groups (where available) before the establishment of a measure and one/two weeks after the establishment—or the end—of the measure, as SARS-CoV-2 has an average incubation period of 6–7 days in average.
We studied the effect of the containment measures in each of the ten European regions on the total number of clusters (MTCs, N), as well as the total number of sequences that were found to be clustered (Number of sequences in each MTC). In all regions that chose lockdown as a containment measure, we observed an increase of the clustered sequences referring to the age groups 50–64 and 65+, two weeks after the establishment of the measure. Especially in Munich, La Rioja, Vienna, Navarra, and Lombardy, MTCs referring to the age group 65+ were de novo formed two or more weeks after the establishment of the measure. At the same time, in Madrid, an increase of MTCs is observed in the age group 50–64. Moreover, in Munich (DE), clustered sequences referring to cases at the age group 20–34 were eliminated two weeks after implementation of movement restrictions, while those referring to age groups 20–34 and 35–49 were increased after the end of this measure (Figure 3).
As for the closure of the schools, in Vienna (AT) and La Rioja (ES), we observed a decrease (75% and 100% decrease, respectively) in the number of pre-existing clustered sequences referring to the age group 0–19 two weeks after the establishment of the measure, indicating a successful containment of the transmissions within this group. MTCs referring to age group 20–34 were also shrunk. The decrease of the MTCs belonging to age group 20–34 in Munich (DE) coincided and may also be associated with schools’ closure, as MTCs were formed one week after school closure, thereby suggesting that the corresponding infections were traced at one week before the establishment of the measure. This trend was not obvious in Uusimaa (FI) and Reykjavik (IS), where the pre-existing clusters of younger infected individuals persisted two weeks after the establishment of the measure (Figure 3). Of note, in these geographical regions, school closure was not combined with a generalized lockdown.
The national borders closure was generally applied for longer periods and data availability after the end of the measure was limited for most of the geographical regions evaluated in this study. Where available, data from Liege (BE) and Madrid (ES) suggested that after the reopening of the borders the expected increase of the unclustered sequences was followed by new MTCs referring to all age groups and particularly to the age group of 20–34 (100% and 50% increase, respectively), yet excluding young individuals 0–19 y.o. (Figure 3).

4. Discussion

In the present study, we focused on the spatio-temporal characterization of SARS-CoV-2 molecular transmission clusters in ten geographical regions of the European continent and we evaluated the effect of individual containment measures taken in each geographical region on the clustering patterns of the sequenced cases and the virus dispersal profiles. One major factor we focused on was the role of the different age groups of the infected individuals as the contribution of this characteristic in the spread of COVID-19 remains unclear [60]. The containment measures we focused on were the lockdowns, the school closures, and border closures, as these measures have different impacts on different parts of the population, but also directly affect all aspects of the society. To this end, each country should undertake all containment measures needed towards protecting health, both physical and mental [61,62,63], while preventing economic and social disruption and respecting human rights (https://www.who.int/director-general/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19---11-march-2020, accessed on 16 October 2020). At certain points, some countries have had no other option but to issue stay-at-home orders in order to buy time and prepare their health-care systems for the pandemic, test the population in a wide-scale and trace and quarantine contacts (https://www.who.int/news-room/q-a-detail/herd-immunity-lockdowns-and-covid-19, accessed on 16 October 2020).
The social activity profiles of the different age groups are expected to differ. Younger individuals are usually more socially active, thus are expected to be more important contributors in the formation of MTCs [64]. Our findings suggest that this pattern did not apply uniformly across the different settings. Analyzing the distribution of clustered sequences according to the age group they belonged to, we found that younger individuals were more likely to belong to MTCs (Table 2) but only in settings with lower overall levels of total clustering (Saint Petersburg, Liege and Reykjavik) (Figure 2, outer five-color circles). At the same time, the age distribution of the sampled cases in these regions was balanced, as in all regions included in the study (Figure 2, inner five-color circles). Multiple factors could have been associated with the formation of MTCs but by analyzing the demographic and socio-economic characteristics of these as well as of all regions studied, we did not identify other factors that could have driven this trend (Table S1). The fact that especially during the initial phase of the pandemic there was a lot of fear for the elderly—thus they were more protected—could explain the observation that MTCs were primarily formed amongst younger individuals, in these particular settings [65].
Across Europe, the kind and the timing of the containment measures taken as a response to the coronavirus pandemic has differed from country to country as presented in detail in Supplementary Figure S1. Among the first measures taken were travel bans from specific locations. This was followed by local or regional quarantines, calls for self-isolation and social distancing. In some countries, extensive closures and quarantines followed. Based on the results of the present study, in most of the geographical regions that established lockdowns during the first pandemic wave, an increase of the MTCs that refer to the older population of the geographical regions evaluated was observed two weeks after the establishment of the measure (age classes 50–64 and 65+ yrs) (Figure 3). However, it is of higher importance to note that sampling bias may apply on this type of analysis and have an impact on the conclusions drawn, especially since different countries may have drastically divergent sampling and sequencing approaches, which may also change over time.
Taking these limitations into account, our results suggest that the generalized lockdown might be effective in containing the epidemic within the most socially active groups of the population but the impact of this measure on the transmissions associated with older adults warrants further investigation. Similarly, Salazar et al. analyzed the daily number of COVID-19 cases for different age groups before and after the lockdown period, detecting a relative increase in the incidence of detected SARS-CoV-2 infection in age groups 50–54 and 55–59 years, after the initial implementation of the measure [41]. A possible explanation could be the household transmissions that occurred because of the lockdown, especially in geographical regions or areas where most of the households are multi-generational (all Spanish geographical regions and Lombardy in Italy). Another explanation could be the higher employment rates in middle-aged adults compared with the younger or even that younger workforce is more related to employment that could be performed digitally/from home compared to middle-aged population. Adam et al. observed that transmission within family households were more frequent than within social and work settings. They highlighted that social settings were associated with both younger cases and more secondary cases compared to households, although this was not the case for households versus work setting [65]. These trends can only be confirmed by carefully designed controlled studies that specifically measure the changes in the MTCs’ distribution over time.
School closures have been widespread in some countries during influenza pandemics, and many studies report important effects on reducing transmission [66]. Expecting that the benefits observed in influenza outbreaks would also apply to SARS-CoV-2, many countries instituted large-scale or national closure of schools by March, 2020. Previous studies regarding the transmission dynamics and clinical characteristics of SARS-CoV-2 infection in children have suggested that children may become a significant spreader at the explosion stage of the outbreak [66], while at the same time the risk of diffusion of infection from child to child is low and even the transmission from a child to an adult is uncommon [67,68,69]. Moreover, as Piovani et al. suggest [70], combination of mass gathering bans and school closures is associated with the reduction in COVID-19 mortality. In this study, we sought to investigate the effect of school closure in MTCs patterns. Based on our analyses, and for those geographical regions with adequate supporting data, school closure was possibly related to the decrease of clustered sequences referring to the age group 0–19 and/or 20–34 in Vienna (AT), La Rioja (ES), and Munich (DE) (Figure 3). This leads to the hypothesis that clustered spread of the virus was decreased due to the limited interactions of these age groups, but also probably due to limited movements and interactions of young parents taking their children to schools. At the same time, in Uusimaa (FI) and Reykjavik (IS), the pre-existing clusters of younger infected individuals remained the same or increased after the establishment of school closure, in the absence of a generalized lockdown. This finding suggests that the combination of the two measures is more effective than school closures alone, towards the containment of the clustered spread of the virus between younger individuals.
Regarding the closure of borders, the impact towards the containment of the epidemic was not reflected on the clustering levels of the sequenced cases, as this measure is mainly effective in controlling the incoming new cases. Observations could be performed for Madrid and Liege where sequence data were available before the establishment and after the end of the measure. The observation that, after the reopening of the borders, there was an increase of the unclustered sequences, which gradually transformed into MTCs specifically in the age group of 20–34 and generally in the age groups referring to adults, suggests that such measures may have a positive impact in controlling the epidemic, while the reopening of the borders can lead to new introductions (unclustered sequences), which are gradually converted into local spreads, reflected in increased MTCs (Figure 3). This finding is in concordance with the studies performed by Wells et al. [71], John et al. [72] and Wilder-Smith et al. [73], suggesting that this measure has a great impact if enacted during the early stages of the epidemic. The estimation of the average time it takes for an incoming case to become domestic and create new MTCs, could be very useful in the context of designing public health interventions as a respond to these introductions. It is very important to rely on molecular tests for the detection of SARS-CoV-2 in the borders after their reopening. Additionally, it is important to highlight the value of continuing molecular testing via real time PCR methods at the border level for Public Health purposes, which should not be ultimately replaced by the recently introduced rapid Ag testing; the isolation of the viral RNA is a prerequisite for PCR testing, providing the necessary material for monitoring the introduction of newly emerged virus variants, like the B.1.1.7, which was recently emerged in the UK. Emerging variants play a pivotal role in the dispersal patterns of cross-borders transmissions, since their virulence and transmissibility might differ [74].
In this study, we relied on data mainly generated during the first pandemic wave, thus some of our inferences may change as more SARS-CoV-2 complete sequences become available. Moreover, because the number of tests for the detection of SARS-CoV-2 that are being held and the proportion of sequenced cases are dissimilar between the geographic locations included in the study, the total number and the profile of the MTCs isolated is not necessarily comparable between these regions. Sequences available on GISAID are unrelated to clinical data, thus the sampling dates do not necessarily reflect the actual infection dates. Other socioeconomic and environmental factors could have also affected the spread of the virus and the formation of MTCs (Figure 2, Table S1). At the same time, other events, such as concerts, parties, political gatherings and sports events could possibly be linked to MTCs, extensive testing, contact tracing, and complete metadata are needed in order to conclude on such associations. However, despite essential differences between countries’ conditions and policies, common trends emerge that prove MTCs are a valuable tool for virus spread surveillance.
Spatiotemporal analysis of MTCs may reflect virus transmission that has not yet been identified through contact tracing, thus cryptic transmissions can be revealed through this process. The investigation of SARS-CoV-2 MTCs via robust bioinformatics pipelines can be a useful tool in order to focus prevention efforts. Routine use of this systematic method in near real-time can automate the detection of SARS-CoV-2 transmission and merits further investigation regarding guidance of Public Health efforts to contain the spread of the pandemic virus.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-1729/11/3/219/s1, Table S1: Demographic and socio-economic factors in 10 European regions, Figure S1: Containment measures in 10 European regions.

Author Contributions

Conceptualization, M.B., A.D., E.-G.K., D.P. and T.K.; Data curation, M.B., A.D., G.L. and A.A.; Formal analysis, M.B., A.D., E.-G.K., E.A. and T.K.; Methodology, M.B., D.P. and T.K.; Software, M.B., G.M. and T.K.; Supervision, A.B., D.P. and T.K.; Validation, M.B., V.P.; Visualization, M.B., A.D. and T.K.; Writing—original draft, M.B., A.D., E.-G.K., S.T., G.K.N., D.N.S., A.B., D.P. and T.K.; Writing—review and editing, M.B., A.D., E.-G.K., S.T., G.K.N., D.N.S., G.M., V.P., G.P., G.L., A.A., E.A., G.S., A.B., D.P. and T.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Publicly available sequencing data and metadata were downloaded from www.gisaid.org (accessed on 16 October 2020).

Acknowledgments

We gratefully acknowledge the authors, and the originating and submitting laboratories of the sequences from GISAID’s Database on which this research is based.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cui, J.; Li, F.; Shi, Z.-L. Origin and evolution of pathogenic coronaviruses. Nat. Rev. Microbiol. 2018, 17, 181–192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Zhou, H.; Chen, X.; Hu, T.; Li, J.; Song, H.; Liu, Y.; Wang, P.; Liu, D.; Yang, J.; Holmes, E.C.; et al. A novel bat coronavirus closely related to SARS-CoV-2 contains natural insertions at the S1/S2 cleavage site of the spike protein. Curr. Biol. 2020, 30, 2196–2203. [Google Scholar] [CrossRef]
  3. Wu, F.; Zhao, S.; Yu, B.; Chen, Y.M.; Wang, W.; Song, Z.G.; Hu, Y.; Tao, Z.W.; Tian, J.H.; Pei, Y.Y.; et al. A new coronavirus associated with human respiratory disease in China. Nature 2020, 579, 265–269. [Google Scholar] [CrossRef] [Green Version]
  4. Drosten, C.; Preiser, W.; Günther, S.; Schmitz, H.; Doerr, H.W. Severe acute respiratory syndrome: Identification of the etiological agent. Trends Mol. Med. 2003, 9, 325–327. [Google Scholar] [CrossRef] [Green Version]
  5. Gandhi, R.T.; Lynch, J.B.; del Rio, C. Mild or moderate Covid-19. N. Engl. J. Med. 2020, 383, 1757–1766. [Google Scholar] [CrossRef]
  6. de Groot, R.J.; Baker, S.C.; Baric, R.S.; Brown, C.S.; Drosten, C.; Enjuanes, L.; Fouchier, R.A.; Galiano, M.; Gorbalenya, A.E.; Memish, Z.A.; et al. Middle East respiratory syndrome coronavirus (MERS-CoV): Announcement of the coronavirus study group. J. Virol. 2013, 87, 7790–7792. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Shu, Y.; McCauley, J. GISAID: Global initiative on sharing all influenza data from vision to reality. Eurosurveillance 2017, 22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Barreto, M.L. Infectious diseases epidemiology. J. Epidemiol. Community Health 2006, 60, 192–195. [Google Scholar] [CrossRef] [Green Version]
  9. Wang, L.; Wang, Y.; Ye, D.; Liu, Q. Review of the 2019 novel coronavirus (SARS-CoV-2) based on current evidence. Int. J. Antimicrob. Agents 2020, 55, 105948. [Google Scholar] [CrossRef]
  10. Petersen, E.; Koopmans, M.; Go, U.; Hamer, D.H.; Petrosillo, N.; Castelli, F.; Storgaard, M.; Al Khalili, S.; Simonsen, L. Comparing SARS-CoV-2 with SARS-CoV and influenza pandemics. Lancet Infect. Dis. 2020, 20, e238–e244. [Google Scholar] [CrossRef]
  11. Alm, E.; Broberg, E.K.; Connor, T.; Hodcroft, E.B.; Komissarov, A.B.; Maurer-Stroh, S.; Melidou, A.; Neher, R.A.; O’Toole, Á.; Pereyaslov, D.; et al. Geographical and temporal distribution of SARS-CoV-2 clades in the WHO European region, January to June 2020. Eurosurveillance 2020, 25, 1–8. [Google Scholar] [CrossRef]
  12. van Dorp, L.; Acman, M.; Richard, D.; Shaw, L.P.; Ford, E.F.; Ormond, L.; Owen, C.J.; Pang, J.; Tan, C.C.S.; Boshier, F.A.T.; et al. Emergence of genomic diversity and recurrent mutations in SARS-CoV-2. Infect. Genet. E 2020, 83. [Google Scholar] [CrossRef]
  13. Bai, Y.; Jiang, D.; Lon, J.R.; Chen, X.; Hu, M.; Lin, S.; Chen, Z.; Wang, X.; Meng, Y.; Du, H.; et al. Comprehensive evolution and molecular characteristics of a large number of SARS-CoV-2 genomes reveal its epidemic trends. Int. J. Infect. Dis. 2020, 100, 164–173. [Google Scholar] [CrossRef] [PubMed]
  14. Su, Y.C.F.; Anderson, D.E.; Young, B.E.; Barnaby, E.; Linster, M.; Zhu, F.; Jayakumar, J.; Zhuang, Y.; Kalimuddin, S.L.; Jenny, G.H.; et al. Discovery and genomic characterization of a 382-nucleotide deletion in ORF7B and orf8 during the early evolution of SARS-CoV-2. mBio 2020, 11, 1–9. [Google Scholar] [CrossRef]
  15. Alinda, A.K.; Ambutsi, M.J.; Kakai, S.K.; Okoth, P.K. SARS-CoV-2 molecular clock and zoonosis. J. Gen. Mol. Virol. 2020, 10, 1–8. [Google Scholar] [CrossRef]
  16. Li, X.; Wang, W.; Zhao, X.; Zai, J.; Zhao, Q.; Li, Y.; Chaillon, A. Transmission dynamics and evolutionary history of 2019-nCoV. J. Med. Virol. 2020, 92, 501–511. [Google Scholar] [CrossRef] [PubMed]
  17. Zhao, Z.; Li, H.; Wu, X.; Zhong, Y.; Zhang, K.; Zhang, Y.P.; Boerwinkle, E.; Fu, Y.X. Moderate mutation rate in the SARS coronavirus genome and its implications. BMC Evol. Biol. 2004, 4, 1–9. [Google Scholar] [CrossRef] [Green Version]
  18. Dudas, G.; Carvalho, L.M.; Rambaut, A.; Bedford, T. MERS-CoV spillover at the camel-human interface. eLife 2018, 7, e31257. [Google Scholar] [CrossRef] [PubMed]
  19. Virk, R.K.; Jayakumar, J.; Mendenhall, I.H.; Moorthy, M.; Lam, P.; Linster, M.; Lim, J.; Lin, C.; Oon, L.L.E.; Lee, H.K.; et al. Divergent evolutionary trajectories of influenza B viruses underlie their contemporaneous epidemic activity. Proc. Natl. Acad. Sci. USA 2020, 117, 619–628. [Google Scholar] [CrossRef] [Green Version]
  20. Williams, T.C.; Burgers, W.A. SARS-CoV-2 evolution and vaccines: Cause for concern? Lancet Respir. Med. 2021, 2600, 19–21. [Google Scholar] [CrossRef]
  21. Morel, B.; Barbera, P.; Czech, L.; Bettisworth, B.; Hubner, L.; Lutteropp, S.; Serdari, D.; Kostaki, E.G.; Mamais, I.; Kozlov, A.M.; et al. Phylogenetic analysis of SARS-CoV-2 data is difficult. Mol. Biol. Evol. 2020, 1–16. [Google Scholar] [CrossRef]
  22. Paraskevis, D.; Beloukas, A.; Stasinos, K.; Pantazis, N.; de Mendoza, C.; Bannert, N.; Meyer, L.; Zangerle, R.; Gill, J.; Prins, M.; et al. HIV-1 molecular transmission clusters in nine European countries and Canada: Association with demographic and clinical factors. BMC Med. 2019, 17. [Google Scholar] [CrossRef] [Green Version]
  23. Rich, S.N.; Richards, V.L.; Mavian, C.N.; Switzer, W.M.; Magalis, B.R.; Poschman, K.; Geary, S.; Broadway, S.E.; Bennett, S.B.; Blanton, J.; et al. Employing molecular phylodynamic methods to identify and forecast HIV transmission clusters in public health settings: A qualitative study. Viruses 2020, 12, 921. [Google Scholar] [CrossRef] [PubMed]
  24. Wertheim, J.O.; Panneer, N.; France, A.M.; Saduvala, N.; Oster, A.M. Incident infection in high-priority HIV molecular transmission clusters in the United States. AIDS 2020, 34, 1187–1193. [Google Scholar] [CrossRef]
  25. Wertheim, J.O.; Murrell, B.; Mehta, S.R.; Forgione, L.A.; Pond, S.L.K.; Smith, D.M.; Torian, L.V. Growth of HIV-1 Molecular Transmission Clusters in New York City. J. Infect. Dis. 2018, 218, 1943–1953. [Google Scholar] [CrossRef] [Green Version]
  26. Huang, L.; Abe, E.M.; Li, X.X.; Bergquist, R.; Xu, L.; Xue, J.B.; Ruan, Y.; Cao, C.L.; Li, S.Z. Space-time clustering and associated risk factors of pulmonary tuberculosis in southwest China. Infect. Dis. Poverty 2018, 7, 91. [Google Scholar] [CrossRef] [PubMed]
  27. Wertheim, J.O.; Oster, A.M.; Switzer, W.M.; Zhang, C.; Panneer, N.; Campbell, E.; Saduvala, N.; Johnson, J.A.; Heneine, W. Natural selection favoring more transmissible HIV detected in United States molecular transmission network. Nat. Commun. 2019, 10, 1–10. [Google Scholar] [CrossRef] [PubMed]
  28. di Giallonardo, F.; Pinto, A.N.; Keen, P.; Shaik, A.; Carrera, A.; Salem, H.; Selvey, C.; Nigro, S.J.; Frases, N.; Price, K.; et al. Subtype-specific differences in transmission cluster dynamics of HIV-1 B and CRF01_AE in New South Wales, Australia. J. Int. AIDS Soc. 2019, 24, e25655. [Google Scholar] [CrossRef]
  29. Kostaki, E.G.; Gova, M.; Adamis, G.; Xylomenos, G.; Chini, M.; Mangafas, M.; Metallidis, S.; Tsachouridou, O.; Papastamopoulos, V.; Chatzidimitriou, D.; et al. A nationwide study about the dispersal patterns of the predominant HIV-1 subtypes A1 and B in Greece: Inference of the molecular transmission clusters. Viruses 2020, 12, 1183. [Google Scholar] [CrossRef]
  30. Dennis, A.M.; Volz, E.; Frost, S.; Hossain, M.; Poon, A.F.Y.; Rebeiro, P.F.; Vermund, S.H.; Sterling, T.R.; Kalish, M.L. HIV-1 Transmission clustering and phylodynamics highlight the important role of young men who have sex with men. AIDS Res. Hum. Retrovir. 2018, 34, 879–888. [Google Scholar] [CrossRef] [Green Version]
  31. Petersen, A.; Cowan, S.A.; Nielsen, J.; Fischer, T.K.; Fonager, J. Characterisation of HIV-1 transmission clusters and drug-resistant mutations in Denmark, 2004 to 2016. Eurosurveillance 2018, 23, 44. [Google Scholar] [CrossRef]
  32. Furuse, Y.; Sando, E.; Tsuchiya, N.; Miyahara, R.; Yasuda, I.; Ko, Y.K.; Saito, M.; Morimoto, K.; Imamura, T.; Shobugawa, Y.; et al. Clusters of coronavirus disease in communities, Japan, January–April 2020. Emerg. Infect. Dis. 2020, 26, 2176–2179. [Google Scholar] [CrossRef] [PubMed]
  33. Longini, I.M.; Nizam, A.; Xu, S.; Ungchusak, K.; Hanshaoworakul, W.; Cummings, D.A.T.; Halloran, E.M. Containing pandemic influenza at the source. Science 2005, 309, 1083–1087. [Google Scholar] [CrossRef] [Green Version]
  34. Bauch, C.T.; Lloyd-Smith, J.O.; Coffee, M.P.; Galvani, A.P. Dynamically modeling SARS and other newly emerging respiratory illnesses: Past, present, and future. Epidemiology 2005, 16, 791–801. [Google Scholar] [CrossRef] [Green Version]
  35. Karatayev, V.A.; Anand, M.; Bauch, C.T. Local lockdowns outperform global lockdown on the far side of the COVID-19 epidemic curve. Proc. Natl. Acad. Sci. USA 2020, 117, 24575–24580. [Google Scholar] [CrossRef]
  36. Koh, W.C.; Naing, L.; Chaw, L.; Rosledzana, M.A.; Alikhan, M.F.; Jamaludin, S.A.; Amin, F.; Omar, A.; Shazli, A.; Griffith, M.; et al. What do we know about SARS-CoV-2 transmission? A systematic review and meta-analysis of the secondary attack rate and associated risk factors. PLoS ONE 2020, 15, e0240205. [Google Scholar] [CrossRef]
  37. Li, C.; Romagnani, P.; von Brunn, A.; Anders, H.-J. SARS-CoV-2 and Europe: Timing of containment measures for outbreak control. Infection 2020, 48, 483–486. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Primc, K.; Slabe-Erker, R. The success of public health measures in Europe during the COVID-19 pandemic. Sustainability 2020, 12, 4321. [Google Scholar] [CrossRef]
  39. Meier, K.; Glatz, T.; Guijt, M.C.; Piccininni, M.; van der Meulen, M.; Atmar, K.; Jolink, A.-T.C.; Kurth, T.; Rohmann, J.L.; Najafabadi, A.H.Z. Public perspectives on protective measures during the COVID-19 pandemic in the Netherlands, Germany and Italy: A survey study. PLoS ONE 2020, 15, e0236917. [Google Scholar] [CrossRef] [PubMed]
  40. Zhang, J.; Litvinova, M.; Liang, Y.; Wang, Y.; Wang, W.; Zhao, S.; Wu, Q.; Merler, S.; Viboud, C.; Vespignani, A.; et al. Changes in contact patterns shape the dynamics of the COVID-19 outbreak in China. Science 2020, 368, 1481–1486. [Google Scholar] [CrossRef] [PubMed]
  41. de Salazar, P.M.; Gomez-Barroso, D.; Pampaka, D.; Gill, J.M.; Penalver, B.; Fernandez-Escobar, C.; Lipsitch, M.; Larrauri, A.; Goldstein, E.; Hernan, M.A. Lockdown measures and relative changes in the age-specific incidence of SARS-CoV-2 in Spain. Epidemiol. Infect. 2020, 148, e268. [Google Scholar] [CrossRef]
  42. Sannigrahi, S.; Pilla, F.; Basu, B.; Basu, A.S.; Molter, A. Examining the association between socio-demographic composition and COVID-19 fatalities in the European region using spatial regression approach. Sustain. Cities Soc. 2020, 62, 102418. [Google Scholar] [CrossRef]
  43. Hadfield, J.; Megill, C.; Bell, S.M.; Huddleston, J.; Potter, B.; Callender, C.; Sagulenko, P.; Bedford, T.; Neher, R.A. Nextstrain: Real-time tracking of pathogen evolution. Bioinformatics 2018, 34, 4121–4123. [Google Scholar] [CrossRef]
  44. Katoh, K.; Asimenos, G.; Toh, H. Multiple alignment of DNA sequences with MAFFT. Methods Mol Biol. 2009, 537, 39–64. [Google Scholar] [CrossRef]
  45. Nguyen, L.-T.; Schmidt, H.A.; von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2014, 32, 268–274. [Google Scholar] [CrossRef]
  46. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.F.; von Haeseler, A.; Jermiin, L.S. Model finder: Fast model selection for accurate phylogenetic estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Turakhia, Y.; De Maio, N.; Thornlow, B.; Gozashti, L.; Lanfear, R.; Walker, C.R.; Hinrichs, A.S.; Fernandes, J.D.; Borges, R.; Slodkowicz, G.; et al. Stability of SARS-CoV-2 phylogenies. PLoS Genet. 2020. [Google Scholar] [CrossRef] [PubMed]
  48. Pipes, L.; Wang, H.; Huelsenbeck, J.P.; Nielsen, R. Assessing uncertainty in the rooting of the SARS-CoV-2 phylogeny. Mol. Biol. Evol. 2020. [Google Scholar] [CrossRef]
  49. Boni, M.F.; Lemey, P.; Jiang, X.; Lam, T.T.-Y.; Perry, B.W.; Castoe, T.A.; Rambaut, A.; Robertson, D.L. Evolutionary origins of the SARS-CoV-2 sarbecovirus lineage responsible for the COVID-19 pandemic. Nat. Microbiol. 2020, 5, 1408–1417. [Google Scholar] [CrossRef] [PubMed]
  50. Sagulenko, P.; Puller, V.; Neher, R.A. TreeTime: Maximum-likelihood phylodynamic analysis. Virus Evol. 2018, 4. [Google Scholar] [CrossRef]
  51. Jones, L.R.; Manrique, J.M. Assessing SARS-CoV-2 spatial phylogenetic structure: Evidence from RNA and protein sequences. bioRxiv 2020. [Google Scholar] [CrossRef]
  52. Li, T.; Liu, D.; Yang, Y.; Guo, J.; Feng, Y.; Zhang, X.; Cheng, S.; Feng, J. Phylogenetic supertree reveals detailed evolution of SARS-CoV-2. Sci. Rep. 2020, 10, 1–9. [Google Scholar] [CrossRef] [PubMed]
  53. Chalmet, K.; Staelens, D.; Blot, S.; Dinakis, S.; Pelgrom, J.; Plum, J.; Vogalaers, D.; Vandekerckhove, L.; Verhofstede, C. Epidemiological study of phylogenetic transmission clusters in a local HIV-1 epidemic reveals distinct differences between subtype B and non-B infections. BMC Infect. Dis. 2010, 10, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Paraskevis, D.; Nikolopoulos, G.; Sypsa, V.; Psichigiou, M.; Pantavou, K.; Kostaki, E.; Karamitros, T.; Paraskeva, D.; Schneider, J.; Malliori, M.; et al. Molecular investigation of HIV-1 cross-group transmissions during an outbreak among people who inject drugs (2011–2014) in Athens, Greece. Infect. Genet. Evol. 2018, 62, 11–16. [Google Scholar] [CrossRef]
  55. Paradis, E.; Schliep, K. Ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 2019, 35, 526–528. [Google Scholar] [CrossRef]
  56. Schliep, K.P. Phangorn: Phylogenetic analysis in R. Bioinformatics 2011, 27, 592–593. [Google Scholar] [CrossRef] [Green Version]
  57. Yu, G. Using ggtree to Visualize Data on Tree-Like Structures. Curr. Protoc. Bioinform. 2020, 69. [Google Scholar] [CrossRef] [PubMed]
  58. Pollán, M.; Pérez-Gómez, B.; Pastor-Barriuso, R.; Oteo, J.; Hernán, M.A.; Pérez-Olmeda, M.; Sanmartín, J.L.; Fernández-García, A.; Cruz, I.; de Larrea, N.F.; et al. Prevalence of SARS-CoV-2 in Spain (ENE-COVID): A nationwide, population-based seroepidemiological study. Lancet 2020, 396, 535–544. [Google Scholar] [CrossRef]
  59. Lemieux, J.E.; Siddle, K.J.; Shaw, B.M.; Loreth, C.; Schaffner, S.F.; Gladden-Young, A.; Adams, G.; Fink, T.; Tomkins-Tinch, C.H.; Krasilnikova, L.A.; et al. Phylogenetic analysis of SARS-CoV-2 in Boston highlights the impact of superspreading events. Am. Assoc. Adv. Sci. 2021, 371. [Google Scholar] [CrossRef]
  60. Goldstein, E.; Lipsitch, M. Temporal rise in the proportion of younger adults and older adolescents among coronavirus disease (COVID-19) cases following the introduction of physical distancing measures, Germany, March to April 2020. Eurosurveillance 2020, 25. [Google Scholar] [CrossRef]
  61. Singh, S.; Roy, D.; Sinha, K.; Parveen, S.; Sharma, G.; Joshi, G. Impact of COVID-19 and lockdown on mental health of children and adolescents: A narrative review with recommendations. Psychiatry Res. 2020, 293, 113429. [Google Scholar] [CrossRef]
  62. Pierce, M.; Hope, H.; Ford, T.; Hatch, S.; Hotopf, M.; John, A.; Kontopantelis, E.; Webb, R.; Wessely, S.; McManus, S.; et al. Mental health before and during the COVID-19 pandemic: A longitudinal probability sample survey of the UK population. Lancet Psychiatry 2020, 7, 883–892. [Google Scholar] [CrossRef]
  63. Townsend, E. COVID-19 policies in the UK and consequences for mental health. Lancet Psychiatry 2020, 7, 1014–1015. [Google Scholar] [CrossRef]
  64. Leclerc, Q.J.; Fuller, N.M.; Knight, L.E.; Funk, S.; Knight, G.M. What settings have been linked to SARS-CoV-2 transmission clusters? Wellcome Open Res. 2020, 5. [Google Scholar] [CrossRef]
  65. Adam, D.C.; Wu, P.; Wong, J.Y.; Lau, E.H.; Tsang, T.K.; Cauchemez, S.; Leung, G.M.; Cowling, B.J. Clustering and superspreading potential of SARS-CoV-2 infections in Hong Kong. Nat. Med. 2020, 26, 1714–1719. [Google Scholar] [CrossRef] [PubMed]
  66. Viner, R.M.; Russell, S.J.; Croker, H.; Packer, J.; Ward, J.; Stansfield, C.; Mytton, O.; Bonell, C.; Booy, R. School closure and management practices during coronavirus outbreaks including COVID-19: A rapid systematic review. Lancet Child Adolesc. Health 2020, 4, 397–404. [Google Scholar] [CrossRef]
  67. Götzinger, F.; Santiago-García, B.; Noguera-Julián, A.; Lanaspa, M.; Lancella, L.; Carducci, F.I.; Gabrovska, N.; Velizarova, S.; Prunk, P.; Osterman, V.; et al. COVID-19 in children and adolescents in Europe: A multinational, multicentre cohort study. Lancet Child Adolesc. Health 2020, 4, 653–661. [Google Scholar] [CrossRef]
  68. Esposito, S.; Principi, N. Debates Around the Role of School Closures in the Coronavirus 2019 Pandemic-Reply. J. Assoc. Med. Am. Pediatr. 2020. [Google Scholar] [CrossRef]
  69. Cao, Q.; Chen, Y.-C.; Chen, C.-L.; Chiu, C.-H. SARS-CoV-2 infection in children: Transmission dynamics and clinical characteristics. J. Formos. Med. Assoc. 2020, 119, 670–673. [Google Scholar] [CrossRef]
  70. Piovani, D.; Christodoulou, M.N.; Hadjidemetriou, A.; Pantavou, K.; Zaza, P.; Bagos, P.G.; Bonovas, S.; Nikolopoulos, G.K. Effect of early application of social distancing interventions on COVID-19 mortality over the first pandemic wave: An analysis of longitudinal data from 37 countries. J. Infect. 2020. [Google Scholar] [CrossRef]
  71. Wells, C.R.; Sah, P.; Moghadas, S.M.; Pandey, A.; Shoukat, A.; Wang, Y.; Wang, Z.; Meyers, L.A.; Singer, B.H.; Galvani, A.P. Impact of international travel and border control measures on the global spread of the novel 2019 coronavirus outbreak. Proc. Natl. Acad. Sci. USA 2020, 117, 7504–7509. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. St. John, R.K.; King, A.; de Jong, D.; Bodie-Collins, M.; Squires, S.G.; Tam, T.W.S. Border screening for SARS. Emerg. Infect. Dis. 2005, 11, 6–10. [Google Scholar] [CrossRef]
  73. Wilder-Smith, A.; Kee, T.G.; Paton, N.I. Experience of severe acute respiratory syndrome in Singapore: Importation of cases, and defense strategies at the airport. J. Travel Med. 2006, 10, 259–262. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Korber, B.; Fischer, W.M.; Gnanakaran, S.; Yoon, H.; Theiler, J.; Abfalterer, W.; Hengartner, N.; Giorgi, E.E.; Bhattacharya, T.; Foley, B.; et al. Tracking changes in SARS-CoV-2 spike: Evidence that D614G increases infectivity of the COVID-19 virus. Cell 2020, 182, 812–827. [Google Scholar] [CrossRef] [PubMed]
Life 11 00219 g001 550
Figure 1. Rooted phylogenetic inference of full genome SARS-CoV-2 sequences from 10 European geographical regions (a–j) and a globally sampled dataset (in grey). Nodes significantly enriched in geographical region’s specific sequences representing molecular transmission clusters (MTCs) are colored in orange, while nodes with geographical region’s specific unclustered sequences are in blue.
Figure 1. Rooted phylogenetic inference of full genome SARS-CoV-2 sequences from 10 European geographical regions (a–j) and a globally sampled dataset (in grey). Nodes significantly enriched in geographical region’s specific sequences representing molecular transmission clusters (MTCs) are colored in orange, while nodes with geographical region’s specific unclustered sequences are in blue.
Life 11 00219 g001
Life 11 00219 g002 550
Figure 2. Distribution of total sequenced cases (five-color inner circle) and clustered sequences (belonging to MTCs—outer circle) according to age group. Density of population, containment measures, total levels of clustering and averages temperature range for each geographical region are also presented.
Figure 2. Distribution of total sequenced cases (five-color inner circle) and clustered sequences (belonging to MTCs—outer circle) according to age group. Density of population, containment measures, total levels of clustering and averages temperature range for each geographical region are also presented.
Life 11 00219 g002
Life 11 00219 g003 550
Figure 3. Distribution of clustered and unclustered sequences for each age group during the containment measures of schools’ closure (blue), the lockdown (red) and closure of national borders (orange). The virus symbol refers to the first case reported in each of the geographical regions.
Figure 3. Distribution of clustered and unclustered sequences for each age group during the containment measures of schools’ closure (blue), the lockdown (red) and closure of national borders (orange). The virus symbol refers to the first case reported in each of the geographical regions.
Life 11 00219 g003
Table
Table 1. Sequences analyzed and MTCs isolated in 10 European geographical regions.
Table 1. Sequences analyzed and MTCs isolated in 10 European geographical regions.
EU
Geographical Region		Region’s
Specific Sequences, N		Total Sequences in Tree, NSamples Collection PeriodMTCs, NSequences in MTCs, N (%)Number of Sequences in Each MTC, NRegion’s Specific Sequences Enrichment, N (%)MTC Clades
Munich (DE)19580322020-03-02 to 2020-05-26551 (26.0)28
17
12
6
5		19 (67.9)
16 (94.1)
8 (66.7)
4 (66.7)
4 (80.0)		G 32/51
GR 19/51
Vienna (AT)14929902020-02-26 to 2020-04-14882 (55.0)44
10
9
9
8
6
6
5		38 (86.4)
9 (90.0)
8 (88.9)
8 (88.9)
6 (75.0)
5 (83.3)
4 (66.7)
4 (80.0)		G 9/82
GR 59/82
S 14/82
Navarra (ES)10936462020-03-07 to 2020-03-29333 (30.2)24
11
7		17 (70.8)
11 (100.0)
5 (71.4)		GR 17/33
S 16/33
La Rioja (ES)25637792020-02-29 to 2020-04-049173 (67.5)68
51
22
21
16
9
8
6
6		62 (91.2)
35 (68.6)
20 (90.9)
19 (90.5)
15 (93.7)
8 (88.9)
6 (75.0)
4 (66.7)
4 (66.7)		GH 74/173
S 99/173
Lombardy (IT)41233332020-02-20 to 2020-05-1014101 (24.5)16
12
12
12
11
9
8
8
7
6
6
5
5
5		12 (75.0)
12 (100.0)
11 (91.7)
8 (66.7)
8 (72.7)
8 (88.9)
7 (87.5)
6 (75.0)
6 (85.7)
5 (83.3)
5 (83.3)
5 (100.0)
4 (80.0)
4 (80.0)		G 81/101
GR 20/101
Uusimaa (FI)22729792020-03-13 to 2020-05-167202 (88.9)156
66
40
6
6
6
6
5		109 (69.9)
44 (66.7)
31 (77.5)
5 (83.3)
5 (83.3)
4 (66.7)
4 (66.7)
4 (80.0)		G 17/202
GH 185/202
Madrid (ES)58240252020-02-25 to 2020-08-3020163 (28.0)26
20
19
18
16
14
13
13
9
8
7
6
6
6
6
5
5
5
5
5		17 (65.4)
14 (70.0)
13 (68.4)
13 (72.2)
13 (81.2)
13 (92.9)
11 (84.6)
10 (76.9)
8 (88.9)
7 (87.5)
6 (85.7)
5 (83.3)
5 (83.3)
4 (66.7)
4 (66.7)
4 (80.0)
4 (80.0)
4 (80.0)
4 (80.0)
4 (80.0)		G 108/163
S 50/163
V 5/163
Saint Petersburg (RU)26729672020-03-13 to 2020-06-16213 (4.8)11
5		9 (81.8)
4 (80.0)		G 4/13
GR 9/13
Liege (BE)53533492020-03-05 to 2020-09-25336 (6.7)27
10
6		21 (77.8)
10 (100.0)
5 (83.3)		G 10/36
GR 26/36
Reykjavik (IS)60130442020-02-27 to 2020-03-296122 (20.2)62
47
14
12
9
7		55 (88.7)
31 (66.0)
12 (85.7)
12 (100.0)
6 (66.7)
6 (85.7)		G 62/122
GR 12/122
L 6/122
O 2/122
S 6/122
V 34/122
Table
Table 2. Distribution of clustered (belonging to MTCs) and unclustered sequences, according to age and gender of sequenced cases in ten European regions.
Table 2. Distribution of clustered (belonging to MTCs) and unclustered sequences, according to age and gender of sequenced cases in ten European regions.
Munich (DE)Vienna (AT)Navarra (ES)La Rioja (ES)Lombardy (IT)Uusimaa (FI)Madrid (ES)Saint Petersburg (RU)Liege (BE)Reykjavik (IS)Total CL/UN
Clustered/Unclustered sequences
AGE0–191/45/40/111/31/317/31/80/231/1518/3368/133
20–349/288/78/618/75/1060/427/600/4418/9529/113194/405
35–499/2617/2415/2136/2017/3938/225/708/495/8741/157233/594
50–6413/3729/77/2640/2119/5730/436/1103/746/7025/131243/615
65–10019/4833/245/2027/3354/18441/570/1700/475/798/21295/693
Overall51/14392/6635/7413/7496/293186/18159/41811/23735/346121/4551033/2440
Total1941581092063892045772483815763473
Chi’s p-valueNSNSNSNSNSNSNS****NS
GENDERFemale25/5442/2621/4165/2940/13290/984/2066/14311/15161/214503/1144
Male26/9040/4012/3461/3556/16295/979/2137/9819/16360/241512/1262
Overall51/14482/6633/75126/6496/294185/18163/41913/24130/314121/4551015/2406
Total1951481081903902035822543445763421
Chi’s p-valueNSNSNSNSNSNSNSNSNSNSNS
Note: NS: not significant (p > 0.05), *: p ≤ 0.05, **: p < 0.01. Statistically significant associations in bold.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bousali, M.; Dimadi, A.; Kostaki, E.-G.; Tsiodras, S.; Nikolopoulos, G.K.; Sgouras, D.N.; Magiorkinis, G.; Papatheodoridis, G.; Pogka, V.; Lourida, G.; et al. SARS-CoV-2 Molecular Transmission Clusters and Containment Measures in Ten European Regions during the First Pandemic Wave. Life 2021, 11, 219. https://doi.org/10.3390/life11030219

AMA Style

Bousali M, Dimadi A, Kostaki E-G, Tsiodras S, Nikolopoulos GK, Sgouras DN, Magiorkinis G, Papatheodoridis G, Pogka V, Lourida G, et al. SARS-CoV-2 Molecular Transmission Clusters and Containment Measures in Ten European Regions during the First Pandemic Wave. Life. 2021; 11(3):219. https://doi.org/10.3390/life11030219

Chicago/Turabian Style

Bousali, Maria, Aristea Dimadi, Evangelia-Georgia Kostaki, Sotirios Tsiodras, Georgios K. Nikolopoulos, Dionyssios N. Sgouras, Gkikas Magiorkinis, George Papatheodoridis, Vasiliki Pogka, Giota Lourida, and et al. 2021. "SARS-CoV-2 Molecular Transmission Clusters and Containment Measures in Ten European Regions during the First Pandemic Wave" Life 11, no. 3: 219. https://doi.org/10.3390/life11030219

APA Style

Bousali, M., Dimadi, A., Kostaki, E. -G., Tsiodras, S., Nikolopoulos, G. K., Sgouras, D. N., Magiorkinis, G., Papatheodoridis, G., Pogka, V., Lourida, G., Argyraki, A., Angelakis, E., Sourvinos, G., Beloukas, A., Paraskevis, D., & Karamitros, T. (2021). SARS-CoV-2 Molecular Transmission Clusters and Containment Measures in Ten European Regions during the First Pandemic Wave. Life, 11(3), 219. https://doi.org/10.3390/life11030219

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop