Early phylogenetic estimate of the effective reproduction number of SARS-CoV-2

Fulminant myocarditis induced by SARS-CoV-2 infection without severe lung involvement: insights into COVID-19 pathogenesis

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection often leads to pulmonary complications. Cardiovascular sequelae, including myocarditis and heart failure, have also been reported. Here, the study presents two fulminant myocarditis cases infected by SARS-CoV-2 exhibiting remarkable elevation of cardiac biomarkers without significant pulmonary injury, as determined by imaging examinations. Immunohistochemical staining reveals the viral antigen within cardiomyocytes, indicating that SARS-CoV-2 could directly infect the myocardium. The full viral genomes from respiratory, anal, and myocardial specimens are obtained via next-generation sequencing. Phylogenetic analyses of the whole genome and spike gene indicate that viruses in the myocardium/pericardial effusion and anal swabs are closely related and cluster together yet diverge from those in the respiratory samples. In addition, unique mutations are found in the anal/myocardial strains compared to the respiratory strains, suggesting tissue-specific virus mutation and adaptation. These findings indicate genetically distinct SARS-CoV-2 variants have infiltrated and disseminated within myocardial tissues, independent of pulmonary injury, and point to different infection routes between the myocardium and respiratory tract, with myocardial infections potentially arising from intestinal infection. These findings highlight the potential for systemic SARS-CoV-2 infection and the importance of a thorough multi-organ assessment in patients for a comprehensive understanding of the pathogenesis of COVID-19.

Host genetics and COVID-19 severity: increasing the accuracy of latest severity scores by Boolean quantum features

The impact of common and rare variants in COVID-19 host genetics has been widely studied. In particular, in Fallerini et al. (Human genetics, 2022, 141, 147-173), common and rare variants were used to define an interpretable machine learning model for predicting COVID-19 severity. First, variants were converted into sets of Boolean features, depending on the absence or the presence of variants in each gene. An ensemble of LASSO logistic regression models was used to identify the most informative Boolean features with respect to the genetic bases of severity. After that, the Boolean features, selected by these logistic models, were combined into an Integrated PolyGenic Score (IPGS), which offers a very simple description of the contribution of host genetics in COVID-19 severity.. IPGS leads to an accuracy of 55%-60% on different cohorts, and, after a logistic regression with both IPGS and age as inputs, it leads to an accuracy of 75%. The goal of this paper is to improve the previous results, using not only the most informative Boolean features with respect to the genetic bases of severity but also the information on host organs involved in the disease. In this study, we generalize the IPGS adding a statistical weight for each organ, through the transformation of Boolean features into "Boolean quantum features," inspired by quantum mechanics. The organ coefficients were set via the application of the genetic algorithm PyGAD, and, after that, we defined two new integrated polygenic scores (IPGS p h 1 and IPGS p h 2 ). By applying a logistic regression with both IPGS, (IPGS p h 2 (or indifferently IPGS p h 1 ) and age as inputs, we reached an accuracy of 84%-86%, thus improving the results previously shown in Fallerini et al. (Human genetics, 2022, 141, 147-173) by a factor of 10%.

Estimates of early outbreak-specific SARS-CoV-2 epidemiological parameters from genomic data

We estimate the basic reproductive number and case counts for 15 distinct Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) outbreaks, distributed across 11 populations (10 countries and one cruise ship), based solely on phylodynamic analyses of genomic data. Our results indicate that, prior to significant public health interventions, the reproductive numbers for 10 (out of 15) of these outbreaks are similar, with median posterior estimates ranging between 1.4 and 2.8. These estimates provide a view which is complementary to that provided by those based on traditional line listing data. The genomic-based view is arguably less susceptible to biases resulting from differences in testing protocols, testing intensity, and import of cases into the community of interest. In the analyses reported here, the genomic data primarily provide information regarding which samples belong to a particular outbreak. We observe that once these outbreaks are identified, the sampling dates carry the majority of the information regarding the reproductive number. Finally, we provide genome-based estimates of the cumulative number of infections for each outbreak. For 7 out of 11 of the populations studied, the number of confirmed cases is much bigger than the cumulative number of infections estimated from the sequence data, a possible explanation being the presence of unsequenced outbreaks in these populations.

Development and efficacy assessment of hand sanitizers and polylactic acid films incorporating caffeic acid and vanillin for enhanced antiviral properties against HCoV-229E

BACKGROUND

Although three years after the outbreak of SARS-CoV-2, the virus is still having a significant impact on human health and the global economy. Infection through respiratory droplets is the main transmission route, but the transmission of the virus by surface contact cannot be ignored. Hand sanitizers and antiviral films can be applied to control SARS-CoV-2, but sanitizers and films show drawbacks such as resistance of the virus against ethanol and environmental problems including the overuse of plastics. Therefore, this study suggested applying natural substrates to hand sanitizers and antiviral films made of biodegradable plastic (PLA). This approach is expected to provide advantages for the easy control of SARS-CoV-2 through the application of natural substances.

METHODS

Antiviral disinfectants and films were manufactured by adding caffeic acid and vanillin to ethanol, isopropyl alcohol, benzalkonium chloride, and PLA. Antiviral efficacies were evaluated with slightly modified international standard testing methods EN 14,476 and ISO 21,702.

RESULTS

In suspension, all the hand sanitizers evaluated in this study showed a reduction of more than 4 log within 2 min against HCoV-229E. After natural substances were added to the hand sanitizers, the time needed to reach the detection limit of the viral titer was shortened both in suspension and porcine skin. However, no difference in the time needed to reach the detection limit of the viral titer was observed in benzalkonium chloride. In the case of antiviral films, those made using both PLA and natural substances showed a 1 log reduction of HCoV-229E compared to the neat PLA film for all treatment groups. Furthermore, the influence of the organic load was evaluated according to the number of contacts of the antiviral products with porcine skin. Ten rubs on the skin resulted in slightly higher antiviral activity than 50 rubs.

CONCLUSION

This study revealed that caffeic acid and vanillin can be effectively used to control HCoV-229E for hand sanitizers and antiviral films. In addition, it is recommended to remove organic matter from the skin for maintaining the antiviral activity of hand sanitizer and antiviral film as the antiviral activity decreased as the organic load increased in this study.

SARS-CoV-2 Recombinants: Genomic Comparison between XBF and Its Parental Lineages

Recombination events are very common and represent one of the primary drivers of RNA virus evolution. The XBF SARS-CoV-2 lineage is one of the most recently generated recombinants during the COVID-19 pandemic. It is a recombinant of BA.5.2.3 and BA.2.75.3, both descendants of lineages that caused many concerns (BA.5 and BA.2.75, respectively). Here, we performed a genomic survey focused on comparing the recombinant XBF with its parental lineages to provide a comprehensive assessment of the evolutionary potential, epidemiological trajectory, and potential risks. Genetic analyses indicated that although XBF initially showed the typical expansion depicted by a steep curve, causing several concerns, currently there is no indication of significant expansion potential or a contagion rate surpassing that of other currently active or previously prevalent lineages. BSP indicated that the peak has been reached around 19 October 2022 and then the genetic variability suffered slight oscillations until early 5 March 2023 when the population size reduced for the last time starting its last plateau that is still lasting. Structural analyses confirmed its reduced potential, also indicating that properties of NTDs and RBDs of XBF and its parental lineages present no significant difference. Of course, cautionary measures must still be taken and genome-based monitoring remains the best tool for detecting any important changes in viral genome composition.

Molecular In-Depth on the Epidemiological Expansion of SARS-CoV-2 XBB.1.5

Since the beginning of the pandemic, the generation of new variants periodically recurs. The XBB.1.5 SARS-CoV-2 variant is one of the most recent. This research was aimed at verifying the potential hazard of this new subvariant. To achieve this objective, we performed a genome-based integrative approach, integrating results from genetic variability/phylodynamics with structural and immunoinformatic analyses to obtain as comprehensive a viewpoint as possible. The Bayesian Skyline Plot (BSP) shows that the viral population size reached the plateau phase on 24 November 2022, and the number of lineages peaked at the same time. The evolutionary rate is relatively low, amounting to 6.9 × 10−4 subs/sites/years. The NTD domain is identical for XBB.1 and XBB.1.5 whereas their RBDs only differ for the mutations at position 486, where the Phe (in the original Wuhan) is replaced by a Ser in XBB and XBB.1, and by a Pro in XBB.1.5. The variant XBB.1.5 seems to spread more slowly than sub-variants that have caused concerns in 2022. The multidisciplinary molecular in-depth analyses on XBB.1.5 performed here does not provide evidence for a particularly high risk of viral expansion. Results indicate that XBB.1.5 does not possess features to become a new, global, public health threat. As of now, in its current molecular make-up, XBB.1.5 does not represent the most dangerous variant.

Molecular Evolution of SARS-CoV-2 during the COVID-19 Pandemic

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) produced diverse molecular variants during its recent expansion in humans that caused different transmissibility and severity of the associated disease as well as resistance to monoclonal antibodies and polyclonal sera, among other treatments. In order to understand the causes and consequences of the observed SARS-CoV-2 molecular diversity, a variety of recent studies investigated the molecular evolution of this virus during its expansion in humans. In general, this virus evolves with a moderate rate of evolution, in the order of 10-3-10-4 substitutions per site and per year, which presents continuous fluctuations over time. Despite its origin being frequently associated with recombination events between related coronaviruses, little evidence of recombination was detected, and it was mostly located in the spike coding region. Molecular adaptation is heterogeneous among SARS-CoV-2 genes. Although most of the genes evolved under purifying selection, several genes showed genetic signatures of diversifying selection, including a number of positively selected sites that affect proteins relevant for the virus replication. Here, we review current knowledge about the molecular evolution of SARS-CoV-2 in humans, including the emergence and establishment of variants of concern. We also clarify relationships between the nomenclatures of SARS-CoV-2 lineages. We conclude that the molecular evolution of this virus should be monitored over time for predicting relevant phenotypic consequences and designing future efficient treatments.

Long-term SARS-CoV-2 surveillance in the wastewater of Stockholm: What lessons can be learned from the Swedish perspective?

Wastewater-based epidemiology (WBE) can be used to track the spread of SARS-CoV-2 in a population. This study presents the learning outcomes from over two-year long monitoring of SARS-CoV-2 in Stockholm, Sweden. The three main wastewater treatment plants in Stockholm, with a total of six inlets, were monitored from April 2020 until June 2022 (in total 600 samples). This spans five major SARS-CoV-2 waves, where WBE data provided early warning signals for each wave. Further, the measured SARS-CoV-2 content in the wastewater correlated significantly with the level of positive COVID-19 tests (r = 0.86; p << 0.0001) measured by widespread testing of the population. Moreover, as a proof-of-concept, six SARS-CoV-2 variants of concern were monitored using hpPCR assay, demonstrating that variants can be traced through wastewater monitoring. During this long-term surveillance, two sampling protocols, two RNA concentration/extraction methods, two calculation approaches, and normalization to the RNA virus Pepper mild mottle virus (PMMoV) were evaluated. In addition, a study of storage conditions was performed, demonstrating that the decay of viral RNA was significantly reduced upon the addition of glycerol to the wastewater before storage at -80 °C. Our results provide valuable information that can facilitate the incorporation of WBE as a prediction tool for possible future outbreaks of SARS-CoV-2 and preparations for future pandemics.

Tracing the First Days of the SARS-CoV-2 Pandemic in Greece and the Role of the First Imported Group of Travelers

The first SARS-CoV-2 case in Greece was confirmed on February 26, 2020, and since then, multiple strains have circulated the country, leading to regional and country-wide outbreaks. Our aim is to enlighten the events that took place during the first days of the SARS-CoV-2 pandemic in Greece, focusing on the role of the first imported group of travelers. We used whole-genome SARS-CoV-2 sequences obtained from the infected travelers of the group as well as Greece-derived and globally subsampled sequences and applied dedicated phylogenetics and phylodynamics tools as well as in-house-developed bioinformatics pipelines. Our analyses reveal the genetic variants circulating in Greece during the first days of the pandemic and the role of the group's imported strains in the course of the first pandemic wave in Greece. The strain that dominated in Greece throughout the first wave, bearing the D614G mutation, was primarily imported from a certain group of travelers, while molecular and clinical data suggest that the infection of the travelers occurred in Egypt. Founder effects early in the pandemic are important for the success of certain strains, as those arriving early, several times, and to diverse locations lead to the formation of large transmission clusters that can be estimated using molecular epidemiology approaches and can be a useful surveillance tool for the prioritization of nonpharmaceutical interventions and combating present and future outbreaks. IMPORTANCE The strain that dominated in Greece during the first pandemic wave was primarily imported from a group of returning travelers in February 2020, while molecular and clinical data suggest that the origin of the transmission was Egypt. The observed molecular transmission clusters reflect the transmission dynamics of this particular strain bearing the D614G mutation while highlighting the necessity of their use as a surveillance tool for the prioritization of nonpharmaceutical interventions and combating present and future outbreaks.

Genetic and Structural Data on the SARS-CoV-2 Omicron BQ.1 Variant Reveal Its Low Potential for Epidemiological Expansion

The BQ.1 SARS-CoV-2 variant, also known as Cerberus, is one of the most recent Omicron descendant lineages. Compared to its direct progenitor BA.5, BQ.1 has some additional spike mutations in some key antigenic sites, which confer further immune escape ability over other circulating lineages. In such a context, here, we perform a genome-based survey aimed at obtaining a complete-as-possible nuance of this rapidly evolving Omicron subvariant. Genetic data suggest that BQ.1 represents an evolutionary blind background, lacking the rapid diversification that is typical of a dangerous lineage. Indeed, the evolutionary rate of BQ.1 is very similar to that of BA.5 (7.6 × 10-4 and 7 × 10-4 subs/site/year, respectively), which has been circulating for several months. The Bayesian Skyline Plot reconstruction indicates a low level of genetic variability, suggesting that the peak was reached around 3 September 2022. Concerning the affinity for ACE2, structure analyses (also performed by comparing the properties of BQ.1 and BA.5 RBD) indicate that the impact of the BQ.1 mutations may be modest. Likewise, immunoinformatic analyses showed moderate differences between the BQ.1 and BA5 potential B-cell epitopes. In conclusion, genetic and structural analyses on SARS-CoV-2 BQ.1 suggest no evidence of a particularly dangerous or high expansion capability. Genome-based monitoring must continue uninterrupted for a better understanding of its descendants and all other lineages.

Structural analysis of a simplified model reproducing SARS-CoV-2 S RBD/ACE2 binding site

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an RNA virus identified as the cause of the coronavirus outbreak in December 2019 (COVID-19). Like all the RNA viruses, SARS-CoV-2 constantly evolves through mutations in its genome, accumulating 1–2 nucleotide changes every month, giving the virus a selective advantage through enhanced transmissibility, greater pathogenicity, and the possibility of circumventing immunity previously acquired by an individual either by natural infection or by vaccination. Several SARS-CoV-2 variants of concern (VoC) have been identified, among which we find Alpha (Lineage B.1.1.7), Beta (Lineage B.1.351), and Gamma (Lineage P.1) variants. Most of the mutations occur in the spike (S) protein, a surface glycoprotein that plays a crucial role in viral infection; the S protein binds the host cell receptor, the angiotensin-converting enzyme of type 2 (ACE2) via the receptor binding domain (RBD) and catalyzes the fusion of the viral membrane with the host cell. In this work, we present the development of a simplified system that would afford to study the change in the SARS-CoV-2 S RBD/ACE2 binding related to the frequent mutations. In particular, we synthesized and studied the structure of short amino acid sequences, mimicking the two proteins’ critical portions. Variations in the residues were easily managed through the one-point alteration of the sequences. Nuclear magnetic resonance (NMR) and circular dichroism (CD) spectroscopies provide insights into ACE2 and SARS-CoV-2 S RBD structure with its related three variants (Alpha, Beta, and Gamma). Spectroscopy data supported by molecular dynamics lead to the description of an ACE2/RBD binding model in which the effect of a single amino acid mutation in changing the binding of S protein to the ACE2 receptor is predictable.

Mathematical modeling to study the impact of immigration on the dynamics of the COVID-19 pandemic: A case study for Venezuela

We propose two different mathematical models to study the effect of immigration on the COVID-19 pandemic. The first model does not consider immigration, whereas the second one does. Both mathematical models consider five different subpopulations: susceptible, exposed, infected, asymptomatic carriers, and recovered. We find the basic reproduction number R₀ using the next-generation matrix method for the mathematical model without immigration. This threshold parameter is paramount because it allows us to characterize the evolution of the disease and identify what parameters substantially affect the COVID-19 pandemic outcome. We focus on the Venezuelan scenario, where immigration and emigration have been important over recent years, particularly during the pandemic. We show that the estimation of the transmission rates of the SARS-CoV-2 are affected when the immigration of infected people is considered. This has an important consequence from a public health perspective because if the basic reproduction number is less than unity, we can expect that the SARS-CoV-2 would disappear. Thus, if the basic reproduction number is slightly above one, we can predict that some mild non-pharmaceutical interventions would be enough to decrease the number of infected people. The results show that the dynamics of the spread of SARS-CoV-2 through the population must consider immigration to obtain better insight into the outcomes and create awareness in the population regarding the population flow.

In silico evaluation of Philippine Natural Products against SARS-CoV-2 Main Protease

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, is a novel strain of coronavirus first reported in December 2019 which rapidly spread throughout the world and was subsequently declared a pandemic by the World Health Organization (WHO) in March 2020. Although vaccines, as well as treatments, have been rapidly developed and deployed, these are still spread thin, especially in the developing world. There is also a continuing threat of the emergence of mutated variants which may not be as responsive to available vaccines and drugs. Accessible and affordable sources of antiviral drugs against SARS-CoV-2 offer wider options for the clinical treatment of populations at risk for severe COVID-19. Using in silico methods, this study identified potential inhibitors against the SARS-CoV-2 main protease (Mpro), the protease directly responsible for the activation of the viral replication enzyme, from a consolidated database of 1516 Philippine natural products. Molecular docking experiments, along with in silico ADME predictions, determined top ligands from this database with the highest potential inhibitory effects against Mpro. Molecular dynamic trajectories of the apo and diosmetin-7-O-b-D-glucopyranoside (DG) in complex with the protein predicted potential mechanisms of action for the ligand-by separating the Cys145-His41 catalytic dyad and by influencing the protein network through key intra-signaling residues within the Mpro binding site. These findings show the inhibitory potential of DG against the SARS-CoV-2 Mpro, and further validation is recommended through in vitro or in vivo experimentation.

Estimation of the Serial Interval and the Effective Reproductive Number of COVID-19 Outbreak Using Contact Data in Burkina Faso, a Sub-Saharan African Country

The COVID-19 outbreak has spread all around the world in less than four months. However, the pattern of the epidemic was different according to the countries. We propose this paper to describe the transmission network and to estimate the serial interval and the reproductive number of the novel coronavirus disease (COVID-19) in Burkina Faso, a Sub-Saharan African country. Data from the COVID-19 response team was analyzed. Information on the 804 first detected cases were pulled together. From contact tracing information, 126 infector-infectee pairs were built. The principal infection clusters with their index cases were observed, principally the two major identified indexes in Burkina. However, the generations of infections were usually short (less than four). The serial interval was estimated to follow a gamma distribution with a shape parameter 1.04 (95% credibility interval: 0.69–1.57) and a scale parameter of 5.69 (95% credibility interval: 3.76–9.11). The basic reproductive number was estimated at 2.36 (95% confidence interval: 1.46–3.26). However, the effective reproductive number decreases very quickly, reaching a minimum value of 0.20 (95% confidence interval: 0.06–0.34). Estimated parameters are made available to monitor the outbreak in Sub-Saharan African countries. These show serial intervals like in the other continents but less infectiousness.

Commentary on the use of the reproduction number R during the COVID-19 pandemic

Since the beginning of the COVID-19 pandemic, the reproduction number [Formula: see text] has become a popular epidemiological metric used to communicate the state of the epidemic. At its most basic, [Formula: see text] is defined as the average number of secondary infections caused by one primary infected individual. [Formula: see text] seems convenient, because the epidemic is expanding if [Formula: see text] and contracting if [Formula: see text]. The magnitude of [Formula: see text] indicates by how much transmission needs to be reduced to control the epidemic. Using [Formula: see text] in a naïve way can cause new problems. The reasons for this are threefold: (1) There is not just one definition of [Formula: see text] but many, and the precise definition of [Formula: see text] affects both its estimated value and how it should be interpreted. (2) Even with a particular clearly defined [Formula: see text], there may be different statistical methods used to estimate its value, and the choice of method will affect the estimate. (3) The availability and type of data used to estimate [Formula: see text] vary, and it is not always clear what data should be included in the estimation. In this review, we discuss when [Formula: see text] is useful, when it may be of use but needs to be interpreted with care, and when it may be an inappropriate indicator of the progress of the epidemic. We also argue that careful definition of [Formula: see text], and the data and methods used to estimate it, can make [Formula: see text] a more useful metric for future management of the epidemic.

matOptimize: a parallel tree optimization method enables online phylogenetics for SARS-CoV-2

MOTIVATION

Phylogenetic tree optimization is necessary for precise analysis of evolutionary and transmission dynamics, but existing tools are inadequate for handling the scale and pace of data produced during the coronavirus disease 2019 (COVID-19) pandemic. One transformative approach, online phylogenetics, aims to incrementally add samples to an ever-growing phylogeny, but there are no previously existing approaches that can efficiently optimize this vast phylogeny under the time constraints of the pandemic.

RESULTS

Here, we present matOptimize, a fast and memory-efficient phylogenetic tree optimization tool based on parsimony that can be parallelized across multiple CPU threads and nodes, and provides orders of magnitude improvement in runtime and peak memory usage compared to existing state-of-the-art methods. We have developed this method particularly to address the pressing need during the COVID-19 pandemic for daily maintenance and optimization of a comprehensive SARS-CoV-2 phylogeny. matOptimize is currently helping refine on a daily basis possibly the largest-ever phylogenetic tree, containing millions of SARS-CoV-2 sequences.

AVAILABILITY AND IMPLEMENTATION

The matOptimize code is freely available as part of the UShER package (https://github.com/yatisht/usher) and can also be installed via bioconda (https://bioconda.github.io/recipes/usher/README.html). All scripts we used to perform the experiments in this manuscript are available at https://github.com/yceh/matOptimize-experiments.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Variational phylodynamic inference using pandemic-scale data

The ongoing global pandemic has sharply increased the amount of data available to researchers in epidemiology and public health. Unfortunately, few existing analysis tools are capable of exploiting all of the information contained in a pandemic-scale data set, resulting in missed opportunities for improved surveillance and contact tracing. In this paper, we develop the variational Bayesian skyline (VBSKY), a method for fitting Bayesian phylodynamic models to very large pathogen genetic data sets. By combining recent advances in phylodynamic modeling, scalable Bayesian inference and differentiable programming, along with a few tailored heuristics, VBSKY is capable of analyzing thousands of genomes in a few minutes, providing accurate estimates of epidemiologically relevant quantities such as the effective reproduction number and overall sampling effort through time. We illustrate the utility of our method by performing a rapid analysis of a large number of SARS-CoV-2 genomes, and demonstrate that the resulting estimates closely track those derived from alternative sources of public health data.

Early Emergence and Dispersal of Delta SARS-CoV-2 Lineage AY.99.2 in Brazil

The year of 2021 was marked by the emergence and dispersal of a number of SARS-CoV-2 lineages, resulting in the “third wave” of COVID-19 in several countries despite the level of vaccine coverage. Soon after the first confirmed cases of COVID-19 by the Delta variant in Brazil, at least seven Delta sub-lineages emerged, including the globally spread AY.101 and AY.99.2. In this study we performed a detailed analysis of the COVID-19 scenario in Brazil from April to December 2021 by using data collected by the largest private medical diagnostic company in Latin America (Dasa), and SARS-CoV-2 genomic sequences generated by its SARS-CoV-2 genomic surveillance project (GENOV). For phylogenetic and Bayesian analysis, SARS-CoV-2 genomes available at GISAID public database were also retrieved. We confirmed that the Brazilian AY.99.2 and AY.101 were the most prevalent lineages during this period, overpassing the Gamma variant in July/August. We also estimated that AY.99.2 likely emerged a few weeks after the entry of the B.1.617.2 in the country, at some point between late April and May and rapidly spread to other countries. Despite no increased fitness described for the AY.99.2 lineage, a rapid shift in the composition of Delta SARS-CoV-2 lineages prevalence in Brazil took place. Understanding the reasons leading the AY.99.2 to become the dominant lineage in the country is important to understand the process of lineage competitions that may inform future control measures.

SARS-CoV-2 Mutant Spectra at Different Depth Levels Reveal an Overwhelming Abundance of Low Frequency Mutations

Populations of RNA viruses are composed of complex and dynamic mixtures of variant genomes that are termed mutant spectra or mutant clouds. This applies also to SARS-CoV-2, and mutations that are detected at low frequency in an infected individual can be dominant (represented in the consensus sequence) in subsequent variants of interest or variants of concern. Here we briefly review the main conclusions of our work on mutant spectrum characterization of hepatitis C virus (HCV) and SARS-CoV-2 at the nucleotide and amino acid levels and address the following two new questions derived from previous results: (i) how is the SARS-CoV-2 mutant and deletion spectrum composition in diagnostic samples, when examined at progressively lower cut-off mutant frequency values in ultra-deep sequencing; (ii) how the frequency distribution of minority amino acid substitutions in SARS-CoV-2 compares with that of HCV sampled also from infected patients. The main conclusions are the following: (i) the number of different mutations found at low frequency in SARS-CoV-2 mutant spectra increases dramatically (50- to 100-fold) as the cut-off frequency for mutation detection is lowered from 0.5% to 0.1%, and (ii) that, contrary to HCV, SARS-CoV-2 mutant spectra exhibit a deficit of intermediate frequency amino acid substitutions. The possible origin and implications of mutant spectrum differences among RNA viruses are discussed.

Le virus SARS-CoV-2 n’a pas « d’origine »

Depuis le début de la pandémie de COVID-19, la question de l’origine de ce virus fait l’objet d’une vive polémique. C’est la question de « l’origine » qui est biaisée. Darwin a montré qu’il n’y a pas d’origine déterminée à aucune espèce animale ou végétale, simplement un processus évolutif et sélectif. Il en est de même pour les virus, il n’y a pas d’origine, mais un processus évolutif. Les virus circulent d’hôte à hôte, animaux ou humains. Les virus pandémiques circulent déjà chez l’homme et évoluent avant l’apparition d’une maladie. Ce processus évolutif se poursuit et donne naissance à des variants successifs. La solution n’est pas de cibler la maladie ou le possible agent causal, mais plutôt de cibler le processus d’émergence de la maladie lui-même.

Estimating the effect of non-pharmaceutical interventions on US SARS-CoV-2 infections in the first year of the pandemic

SARS-CoV-2 emerged in late 2019 as a zoonotic infection of humans, and proceeded to cause a worldwide pandemic of historic magnitude. Here, we use a simple epidemiological model and consider the full range of initial estimates from published studies for infection and recovery rates, seasonality, changes in mobility, the effectiveness of masks and the fraction of people wearing them. Monte Carlo simulations are used to simulate the progression of possible pandemics and we show a match for the real progression of the pandemic during 2020 with an R 2 of 0.91. The results show that the combination of masks and changes in mobility avoided approximately 248.3 million (σ = 31.2 million) infections in the US before vaccinations became available.

The Occurrence Birth-Death Process for combined-evidence analysis in macroevolution and epidemiology

Phylodynamic models generally aim at jointly inferring phylogenetic relationships, model parameters, and more recently, the number of lineages through time, based on molecular sequence data. In the fields of epidemiology and macroevolution, these models can be used to estimate, respectively, the past number of infected individuals (prevalence) or the past number of species (paleodiversity) through time. Recent years have seen the development of “total-evidence” analyses, which combine molecular and morphological data from extant and past sampled individuals in a unified Bayesian inference framework. Even sampled individuals characterized only by their sampling time, that is, lacking morphological and molecular data, which we call occurrences, provide invaluable information to estimate the past number of lineages. Here, we present new methodological developments around the fossilized birth–death process enabling us to (i) incorporate occurrence data in the likelihood function; (ii) consider piecewise-constant birth, death, and sampling rates; and (iii) estimate the past number of lineages, with or without knowledge of the underlying tree. We implement our method in the RevBayes software environment, enabling its use along with a large set of models of molecular and morphological evolution, and validate the inference workflow using simulations under a wide range of conditions. We finally illustrate our new implementation using two empirical data sets stemming from the fields of epidemiology and macroevolution. In epidemiology, we infer the prevalence of the coronavirus disease 2019 outbreak on the Diamond Princess ship, by taking into account jointly the case count record (occurrences) along with viral sequences for a fraction of infected individuals. In macroevolution, we infer the diversity trajectory of cetaceans using molecular and morphological data from extant taxa, morphological data from fossils, as well as numerous fossil occurrences. The joint modeling of occurrences and trees holds the promise to further bridge the gap between traditional epidemiology and pathogen genomics, as well as paleontology and molecular phylogenetics. [Birth–death model; epidemiology; fossils; macroevolution; occurrences; phylogenetics; skyline.]

Identification of vaccine targets & design of vaccine against SARS-CoV-2 coronavirus using computational and deep learning-based approaches

An unusual pneumonia infection, named COVID-19, was reported on December 2019 in China. It was reported to be caused by a novel coronavirus which has infected approximately 220 million people worldwide with a death toll of 4.5 million as of September 2021. This study is focused on finding potential vaccine candidates and designing an in-silico subunit multi-epitope vaccine candidates using a unique computational pipeline, integrating reverse vaccinology, molecular docking and simulation methods. A protein named spike protein of SARS-CoV-2 with the GenBank ID QHD43416.1 was shortlisted as a potential vaccine candidate and was examined for presence of B-cell and T-cell epitopes. We also investigated antigenicity and interaction with distinct polymorphic alleles of the epitopes. High ranking epitopes such as DLCFTNVY (B cell epitope), KIADYNKL (MHC Class-I) and VKNKCVNFN (MHC class-II) were shortlisted for subsequent analysis. Digestion analysis verified the safety and stability of the shortlisted peptides. Docking study reported a strong binding of proposed peptides with HLA-A*02 and HLA-B7 alleles. We used standard methods to construct vaccine model and this construct was evaluated further for its antigenicity, physicochemical properties, 2D and 3D structure prediction and validation. Further, molecular docking followed by molecular dynamics simulation was performed to evaluate the binding affinity and stability of TLR-4 and vaccine complex. Finally, the vaccine construct was reverse transcribed and adapted for E. coli strain K 12 prior to the insertion within the pET-28-a (+) vector for determining translational and microbial expression followed by conservancy analysis. Also, six multi-epitope subunit vaccines were constructed using different strategies containing immunogenic epitopes, appropriate adjuvants and linker sequences. We propose that our vaccine constructs can be used for downstream investigations using in-vitro and in-vivo studies to design effective and safe vaccine against different strains of COVID-19.

There is no "origin" to SARS-CoV-2

Since the beginning of the COVID-19 pandemic in 2020 caused by SARS-CoV-2, the question of the origin of this virus has been a highly debated issue. Debates have been, and are still, very disputed and often violent between the two main hypotheses: a natural origin through the "spillover" model or a laboratory-leak origin. Tenants of these two options are building arguments often based on the discrepancies of the other theory. The main problem is that it is the initial question of the origin itself which is biased. Charles Darwin demonstrated in 1859 that all species are appearing through a process of evolution, adaptation and selection. There is no determined origin to any animal or plant species, simply an evolutionary and selective process in which chance and environment play a key role. The very same is true for viruses. There is no determined origin to viruses, simply also an evolutionary and selective process in which chance and environment play a key role. However, in the case of viruses the process is slightly more complex because the "environment" is another living organism. Pandemic viruses already circulate in humans prior to the emergence of a disease. They are simply not capable of triggering an epidemic yet. They must evolve in-host, i.e. in-humans, for that. The evolutionary process which gave rise to SARS-CoV-2 is still ongoing with regular emergence of novel variants more adapted than the previous ones. The real relevant question is how these viruses can emerge as pandemic viruses and what the society can do to prevent the future emergence of pandemic viruses.

SARS-CoV-2 Point Mutation and Deletion Spectra and Their Association with Different Disease Outcomes

Mutant spectra of RNA viruses are important to understand viral pathogenesis and response to selective pressures. There is a need to characterize the complexity of mutant spectra in coronaviruses sampled from infected patients. In particular, the possible relationship between SARS-CoV-2 mutant spectrum complexity and disease associations has not been established. In the present study, we report an ultradeep sequencing (UDS) analysis of the mutant spectrum of amplicons from the nsp12 (polymerase)- and spike (S)-coding regions of 30 nasopharyngeal isolates (diagnostic samples) of SARS-CoV-2 of the first COVID-19 pandemic wave (Madrid, Spain, April 2020) classified according to the severity of ensuing COVID-19. Low-frequency mutations and deletions, counted relative to the consensus sequence of the corresponding isolate, were overwhelmingly abundant. We show that the average number of different point mutations, mutations per haplotype, and several diversity indices was significantly higher in SARS-CoV-2 isolated from patients who developed mild disease than in those associated with moderate or severe disease (exitus). No such bias was observed with RNA deletions. Location of amino acid substitutions in the three-dimensional structures of nsp12 (polymerase) and S suggest significant structural or functional effects. Thus, patients who develop mild symptoms may be a richer source of genetic variants of SARS-CoV-2 than patients with moderate or severe COVID-19. IMPORTANCE The study shows that mutant spectra of SARS-CoV-2 from diagnostic samples differ in point mutation abundance and complexity and that significantly larger values were observed in virus from patients who developed mild COVID-19 symptoms. Mutant spectrum complexity is not a uniform trait among isolates. The nature and location of low-frequency amino acid substitutions present in mutant spectra anticipate great potential for phenotypic diversification of SARS-CoV-2.

Phylogeography and genomic epidemiology of SARS-CoV-2 in Italy and Europe with newly characterized Italian genomes between February-June 2020

The aims of this study were to characterize new SARS-CoV-2 genomes sampled all over Italy and to reconstruct the origin and the evolutionary dynamics in Italy and Europe between February and June 2020. The cluster analysis showed only small clusters including < 80 Italian isolates, while most of the Italian strains were intermixed in the whole tree. Pure Italian clusters were observed mainly after the lockdown and distancing measures were adopted. Lineage B and B.1 spread between late January and early February 2020, from China to Veneto and Lombardy, respectively. Lineage B.1.1 (20B) most probably evolved within Italy and spread from central to south Italian regions, and to European countries. The lineage B.1.1.1 (20D) developed most probably in other European countries entering Italy only in the second half of March and remained localized in Piedmont until June 2020. In conclusion, within the limitations of phylogeographical reconstruction, the estimated ancestral scenario suggests an important role of China and Italy in the widespread diffusion of the D614G variant in Europe in the early phase of the pandemic and more dispersed exchanges involving several European countries from the second half of March 2020.

Analysis of Genomic Characteristics of SARS-CoV-2 in Italy, 29 January to 27 March 2020

We performed next-generation sequencing (NGS), phylogenetic analysis, gene flows, and N- and O-glycosylation prediction on SARS-CoV-2 genomes collected from lab-confirmed cases from different Italian regions. To this end, a total of 111 SARS-CoV-2 genomes collected in Italy between 29 January and 27 March 2020 were investigated. The majority of the genomes belonged to lineage B.1, with some descendant lineages. The gene flow analysis showed that the spread occurred mainly from the north to the center and to the south of Italy, as confirmed by epidemiological data. The mean evolutionary rate estimated here was 8.731 × 10⁻⁴ (95% highest posterior density, HPD intervals 5.809 × 10⁻⁴ to 1.19 × 10⁻³), in line with values reported by other authors. The dated phylogeny suggested that SARS-CoV-2 lineage B.1 probably entered Italy between the end of January and early February 2020. Continuous molecular surveillance is needed to trace virus circulation and evolution.

Evolutionary dynamics of the severe acute respiratory syndrome coronavirus 2 genomes

The coronavirus disease 2019 (COVID-19) pandemic has caused immense losses in human lives and the global economy and posed significant challenges for global public health. As severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, has evolved, thousands of single nucleotide variants (SNVs) have been identified across the viral genome. The roles of individual SNVs in the zoonotic origin, evolution, and transmission of SARS-CoV-2 have become the focus of many studies. This review summarizes recent comparative genomic analyses of SARS-CoV-2 and related coronaviruses (SC2r-CoVs) found in non-human animals, including delineation of SARS-CoV-2 lineages based on characteristic SNVs. We also discuss the current understanding of receptor-binding domain (RBD) evolution and characteristic mutations in variants of concern (VOCs) of SARS-CoV-2, as well as possible co-evolution between RBD and its receptor, angiotensin-converting enzyme 2 (ACE2). We propose that the interplay between SARS-CoV-2 and host RNA editing mechanisms might have partially resulted in the bias in nucleotide changes during SARS-CoV-2 evolution. Finally, we outline some current challenges, including difficulty in deciphering the complicated relationship between viral pathogenicity and infectivity of different variants, and monitoring transmission of SARS-CoV-2 between humans and animals as the pandemic progresses.

Epidemiology and Genetic Analysis of SARS-CoV-2 in Myanmar during the Community Outbreaks in 2020

We aimed to analyze the situation of the first two epidemic waves in Myanmar using the publicly available daily situation of COVID-19 and whole-genome sequencing data of SARS-CoV-2. From March 23 to December 31, 2020, there were 33,917 confirmed cases and 741 deaths in Myanmar (case fatality rate of 2.18%). The first wave in Myanmar from March to July was linked to overseas travel, and then a second wave started from Rakhine State, a western border state, leading to the second wave spreading countrywide in Myanmar from August to December 2020. The estimated effective reproductive number (R_t) nationwide reached 6–8 at the beginning of each wave and gradually decreased as the epidemic spread to the community. The whole-genome analysis of 10 Myanmar SARS-CoV-2 strains together with 31 previously registered strains showed that the first wave was caused by GISAID clade O or PANGOLIN lineage B.6 and the second wave was changed to clade GH or lineage B.1.36.16 with a close genetic relationship with other South Asian strains. Constant monitoring of epidemiological situations combined with SARS-CoV-2 genome analysis is important for adjusting public health measures to mitigate the community transmissions of COVID-19.

OUP accepted manuscript

Transmission chains within small urban areas (accommodating ∼30 per cent of the European population) greatly contribute to case burden and economic impact during the ongoing coronavirus pandemic and should be a focus for preventive measures to achieve containment. Here, at very high spatio-temporal resolution, we analysed determinants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission in a European urban area, Basel-City (Switzerland). We combined detailed epidemiological, intra-city mobility and socio-economic data sets with whole-genome sequencing during the first SARS-CoV-2 wave. For this, we succeeded in sequencing 44 per cent of all reported cases from Basel-City and performed phylogenetic clustering and compartmental modelling based on the dominating viral variant (B.1-C15324T; 60 per cent of cases) to identify drivers and patterns of transmission. Based on these results we simulated vaccination scenarios and corresponding healthcare system burden (intensive care unit (ICU) occupancy). Transmissions were driven by socio-economically weaker and highly mobile population groups with mostly cryptic transmissions which lacked genetic and identifiable epidemiological links. Amongst more senior population transmission was clustered. Simulated vaccination scenarios assuming 60–90 per cent transmission reduction and 70–90 per cent reduction of severe cases showed that prioritising mobile, socio-economically weaker populations for vaccination would effectively reduce case numbers. However, long-term ICU occupation would also be effectively reduced if senior population groups were prioritised, provided there were no changes in testing and prevention strategies. Reducing SARS-CoV-2 transmission through vaccination strongly depends on the efficacy of the deployed vaccine. A combined strategy of protecting risk groups by extensive testing coupled with vaccination of the drivers of transmission (i.e. highly mobile groups) would be most effective at reducing the spread of SARS-CoV-2 within an urban area.

Timeline of SARS-CoV-2 Spread in Italy: Results from an Independent Serological Retesting

The massive emergence of COVID-19 cases in the first phase of pandemic within an extremely short period of time suggest that an undetected earlier circulation of SARS-CoV-2 might have occurred. Given the importance of this evidence, an independent evaluation was recommended by the World Health Organization (WHO) to test a subset of samples selected on the level of positivity in ELISA assays (positive, low positive, negative) detected in our previous study of prepandemic samples collected in Italy. SARS-CoV-2 antibodies were blindly retested by two independent centers in 29 blood samples collected in the prepandemic period in Italy, 29 samples collected one year before and 11 COVID-19 control samples. The methodologies used included IgG-RBD/IgM-RBD ELISA assays, a qualitative micro-neutralization CPE-based assay, a multiplex IgG protein array, an ELISA IgM kit (Wantai), and a plaque-reduction neutralization test. The results suggest the presence of SARS-CoV-2 antibodies in some samples collected in the prepandemic period, with the oldest samples found to be positive for IgM by both laboratories collected on 10 October 2019 (Lombardy), 11 November 2019 (Lombardy) and 5 February 2020 (Lazio), the latter with neutralizing antibodies. The detection of IgM and/or IgG binding and neutralizing antibodies was strongly dependent on the different serological assays and thresholds employed, and they were not detected in control samples collected one year before. These findings, although gathered in a small and selected set of samples, highlight the importance of harmonizing serological assays for testing the spread of the SARS-CoV-2 virus and may contribute to a better understanding of future virus dynamics.

Pandemic-scale phylogenetics

Phylogenetics has been central to the genomic surveillance, epidemiology and contact tracing efforts during the COVD-19 pandemic. But the massive scale of genomic sequencing has rendered the pre-pandemic tools inadequate for comprehensive phylogenetic analyses. Here, we discuss the phylogenetic package that we developed to address the needs imposed by this pandemic. The package incorporates several pandemic-specific optimization and parallelization techniques and comprises four programs: UShER, matOptimize, RIPPLES and matUtils. Using high-performance computing, UShER and matOptimize maintain and refine daily a massive mutation-annotated phylogenetic tree consisting of all SARS-CoV-2 sequences available in online repositories. With UShER and RIPPLES, individual labs - even with modest compute resources - incorporate newly-sequenced SARS-CoV-2 genomes on this phylogeny and discover evidence for recombination in real-time. With matUtils, they rapidly query and visualize massive SARS-CoV-2 phylogenies. These tools have empowered scientists worldwide to study the SARS-CoV-2 evolution and transmission at an unprecedented scale, resolution and speed.

Temporal spread and evolution of SARS-CoV-2 in the second pandemic wave in Brazil

Severe acute respiratory syndrome coronavirus‐2 (SARS‐CoV‐2) pandemic spread rapidly and this scenario is concerning in South America, mainly in Brazil that presented more than 21 million coronavirus disease 2019 cases and 590 000 deaths. The recent emergence of novel lineages carrying several mutations in the spike protein has raised additional public health concerns worldwide. The present study describes the temporal spreading and evolution of SARS‐CoV2 in the beginning of the second pandemic wave in Brazil, highlighting the fast dissemination of the two major concerning variants (P.1 and P.2). A total of 2507 SARS‐CoV‐2 whole‐genome sequences (WGSs) with available information from the country (Brazil) and sampling date (July 2020–February 2021), were obtained and the frequencies of the lineages were evaluated in the period of the growing second pandemic wave. The results demonstrated the increasing prevalence of P.1 and P.2 lineages in the period evaluated. P.2 lineage was first detected in the middle of 2020, but a high increase occurred only in the last trimester of this same year and the spreading to all Brazilian regions. P.1 lineage emerged even later, first in the North region in December 2020 and really fast dissemination to all other Brazilian regions in January and February 2021. All SARS‐CoV‐2 WGSs of P.1 and P.2 were further separately evaluated with a Bayesian approach. The rates of nucleotide and amino acid substitutions were statistically higher in P.1 than P.2 (p < 0.01). The phylodynamic analysis demonstrated that P.2 gradually spread in all the country from September 2020 to January 2021, while P.1 disseminated even faster from December 2020 to February 2021. Skyline plots of both lineages demonstrated a slight rise in the spreading for P.2 and exponential growth for P.1. In conclusion, these data demonstrated that the P.1 (recently renamed as Gamma) and P.2 lineages have predominated in the second pandemic wave due to the very high spreading across all geographic regions in Brazil at the end of 2020 and beginning of 2021.

In Brazil, P.1 (Gamma) and P.2 lineages have predominated in the second pandemic wave.

The Bayesian approach showed very high spreading for both lineages across all geographic regions at the end of 2020 and the beginning of 2021.

P.2 increased only in the last trimester of 2020 and the spreading to all Brazilian regions.

P.1 (Gamma) emerged even later with fast dissemination to all Brazilian regions in January and February 2021.

Early reports of epidemiological parameters of the COVID-19 pandemic

Background

The emergence of a new pathogen requires a rapid assessment of its transmissibility, to inform appropriate public health interventions.

Methods

The peer-reviewed literature published between 1 January and 30 April 2020 on COVID-19 in PubMed was searched. Estimates of the incubation period, serial interval and reproduction number for COVID-19 were obtained and compared.

Results

A total of 86 studies met the inclusion criteria. Of these, 33 estimated the mean incubation period (4–7 days) and 15 included estimates of the serial interval (mean 4–8 days; median length 4–5 days). Fifty-two studies estimated the reproduction number. Although reproduction number estimates ranged from 0.3 to 14.8, in 33 studies (63%), they fell between 2 and 3.

Discussion

Studies calculating the incubation period and effective reproduction number were published from the beginning of the pandemic until the end of the study period (30 April 2020); however, most of the studies calculating the serial interval were published in April 2020. The calculated incubation period was similar over the study period and in different settings, whereas estimates of the serial interval and effective reproduction number were setting-specific. Estimates of the serial interval were shorter at the end of the study period as increasing evidence of pre-symptomatic transmission was documented and as jurisdictions enacted outbreak control measures. Estimates of the effective reproduction number varied with the setting and the underlying model assumptions. Early analysis of epidemic parameters provides vital information to inform the outbreak response.

Early Spatiotemporal Patterns and Population Characteristics of the COVID-19 Pandemic in Southeast Asia

This observational study aims to investigate the early disease patterns of coronavirus disease 2019 (COVID-19) in Southeast Asia, consequently providing historical experience for further interventions. Data were extracted from official websites of the WHO and health authorities of relevant countries. A total of 1346 confirmed cases of COVID-19, with 217 recoveries and 18 deaths, were reported in Southeast Asia as of 16 March 2020. The basic reproductive number (R0) of COVID-19 in the region was estimated as 2.51 (95% CI:2.31 to 2.73), and there were significant geographical variations at the subregional level. Early transmission dynamics were examined with an exponential regression model: y = 0.30e0.13x (p < 0.01, R2 = 0.96), which could help predict short-term incidence. Country-level disease burden was positively correlated with Human Development Index (r = 0.86, p < 0.01). A potential early shift in spatial diffusion patterns and a spatiotemporal cluster occurring in Malaysia and Singapore were detected. Demographic analyses of 925 confirmed cases indicated a median age of 44 years and a sex ratio (male/female) of 1.25. Age may play a significant role in both susceptibilities and outcomes. The COVID-19 situation in Southeast Asia is challenging and unevenly geographically distributed. Hence, enhanced real-time surveillance and more efficient resource allocation are urgently needed.

SARS-CoV-2 transmission dynamics in Belarus in 2020 revealed by genomic and incidence data analysis

Background

Non-pharmaceutical interventions (NPIs) have been implemented worldwide to curb COVID-19 spread. Belarus is a rare case of a country with a relatively modern healthcare system, where highly limited NPIs have been enacted. Thus, investigation of Belarusian COVID-19 dynamics is essential for the local and global assessment of the impact of NPI strategies.

Methods

We integrate genomic epidemiology and surveillance methods to investigate the spread of SARS-CoV-2 in Belarus in 2020. We utilize phylodynamics, phylogeography, and probabilistic bias inference to study the virus import and export routes, the dynamics of the effective reproduction number, and the incidence of SARS-CoV-2 infection.

Results

Here we show that the estimated cumulative number of infections by June 2020 exceeds the confirmed case number by a factor of ~4 (95% confidence interval (2; 9)). Intra-country SARS-CoV-2 genomic diversity originates from at least 18 introductions from different regions, with a high proportion of regional transmissions. Phylodynamic analysis indicates a moderate reduction of the effective reproductive number after the introduction of limited NPIs, but its magnitude is lower than for developed countries with large-scale NPIs. On the other hand, the effective reproduction number estimate is comparable with that for the neighboring Ukraine, where NPIs were broader.

Conclusions

The example of Belarus demonstrates how countries with relatively low outward population mobility continue to be integral parts of the global epidemiological environment. Comparison of the effective reproduction number dynamics for Belarus and other countries reveals the effect of different NPI strategies but also emphasizes the role of regional Eastern European sociodemographic factors in the virus spread.

Belarus is one of few European countries that has enacted limited measures to contain SARS-CoV-2, the virus that causes COVID-19. We study the genetic sequences of the SARS-CoV-2 virus circulating in Belarus and other countries in 2020 to investigate how it might have been imported into the country and spread there. We show that the virus was repeatedly imported from and exported to different regions, including a large portion of regional transmissions that occurred despite stricter measures implemented by Belarus’ neighbors. There was a moderate reduction of the virus reproductive number—a measure of virus transmission speed—after April 2020, but its magnitude was lower than for developed countries with more stringent epidemiological interventions. These findings shed light on the COVID-19 spread in Eastern Europe and highlight the impact of public health policies and of regional factors on this spread.

Nemira et al. study the genomic epidemiology and phylodynamics of SARS-CoV-2 in Belarus. They identify potential introduction routes of the virus from other countries, determine that during the first wave of the pandemic the number of infections was likely several times higher than reported case numbers, and estimate the impact of early non-pharmaceutical interventions on SARS-CoV-2 transmission.

COVID-19 pandemic spread against countries' non-pharmaceutical interventions responses: a data-mining driven comparative study

Background

The first half of 2020 has been marked as the era of COVID-19 pandemic which affected the world globally in almost every aspect of the daily life from societal to economical. To prevent the spread of COVID-19, countries have implemented diverse policies regarding Non-Pharmaceutical Intervention (NPI) measures. This is because in the first stage countries had limited knowledge about the virus and its contagiousness. Also, there was no effective medication or vaccines. This paper studies the effectiveness of the implemented policies and measures against the deaths attributed to the virus between January and May 2020.

Methods

Data from the European Centre for Disease Prevention and Control regarding the identified cases and deaths of COVID-19 from 48 countries have been used. Additionally, data concerning the NPI measures related policies implemented by the 48 countries and the capacity of their health care systems was collected manually from their national gazettes and official institutes. Data mining, time series analysis, pattern detection, machine learning, clustering methods and visual analytics techniques have been applied to analyze the collected data and discover possible relationships between the implemented NPIs and COVID-19 spread and mortality. Further, we recorded and analyzed the responses of the countries against COVID-19 pandemic, mainly in urban areas which are over-populated and accordingly COVID-19 has the potential to spread easier among humans.

Results

The data mining and clustering analysis of the collected data showed that the implementation of the NPI measures before the first death case seems to be very effective in controlling the spread of the disease. In other words, delaying the implementation of the NPI measures to after the first death case has practically little effect on limiting the spread of the disease. The success of implementing the NPI measures further depends on the way each government monitored their application. Countries with stricter policing of the measures seems to be more effective in controlling the transmission of the disease.

Conclusions

The conducted comparative data mining study provides insights regarding the correlation between the early implementation of the NPI measures and controlling COVID-19 contagiousness and mortality. We reported a number of useful observations that could be very helpful to the decision makers or epidemiologists regarding the rapid implementation and monitoring of the NPI measures in case of a future wave of COVID-19 or to deal with other unknown infectious pandemics. Regardless, after the first wave of COVID-19, most countries have decided to lift the restrictions and return to normal. This has resulted in a severe second wave in some countries, a situation which requires re-evaluating the whole process and inspiring lessons for the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12889-021-11251-4.

Fundamental Identifiability Limits in Molecular Epidemiology

Viral phylogenies provide crucial information on the spread of infectious diseases, and many studies fit mathematical models to phylogenetic data to estimate epidemiological parameters such as the effective reproduction ratio (Re) over time. Such phylodynamic inferences often complement or even substitute for conventional surveillance data, particularly when sampling is poor or delayed. It remains generally unknown, however, how robust phylodynamic epidemiological inferences are, especially when there is uncertainty regarding pathogen prevalence and sampling intensity. Here, we use recently developed mathematical techniques to fully characterize the information that can possibly be extracted from serially collected viral phylogenetic data, in the context of the commonly used birth-death-sampling model. We show that for any candidate epidemiological scenario, there exists a myriad of alternative, markedly different, and yet plausible "congruent" scenarios that cannot be distinguished using phylogenetic data alone, no matter how large the data set. In the absence of strong constraints or rate priors across the entire study period, neither maximum-likelihood fitting nor Bayesian inference can reliably reconstruct the true epidemiological dynamics from phylogenetic data alone; rather, estimators can only converge to the "congruence class" of the true dynamics. We propose concrete and feasible strategies for making more robust epidemiological inferences from viral phylogenetic data.

An Evolutionary Portrait of the Progenitor SARS-CoV-2 and Its Dominant Offshoots in COVID-19 Pandemic

Global sequencing of genomes of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continued to reveal new genetic variants that are the key to unraveling its early evolutionary history and tracking its global spread over time. Here we present the heretofore cryptic mutational history and spatiotemporal dynamics of SARS-CoV-2 from an analysis of thousands of high-quality genomes. We report the likely most recent common ancestor of SARS-CoV-2, reconstructed through a novel application and advancement of computational methods initially developed to infer the mutational history of tumor cells in a patient. This progenitor genome differs from genomes of the first coronaviruses sampled in China by three variants, implying that none of the earliest patients represent the index case or gave rise to all the human infections. However, multiple coronavirus infections in China and the United States harbored the progenitor genetic fingerprint in January 2020 and later, suggesting that the progenitor was spreading worldwide months before and after the first reported cases of COVID-19 in China. Mutations of the progenitor and its offshoots have produced many dominant coronavirus strains that have spread episodically over time. Fingerprinting based on common mutations reveals that the same coronavirus lineage has dominated North America for most of the pandemic in 2020. There have been multiple replacements of predominant coronavirus strains in Europe and Asia as well as continued presence of multiple high-frequency strains in Asia and North America. We have developed a continually updating dashboard of global evolution and spatiotemporal trends of SARS-CoV-2 spread (http://sars2evo.datamonkey.org/).

Analysis of Delayed Vaccination Regimens: A Mathematical Modeling Approach

The first round of vaccination against coronavirus disease 2019 (COVID-19) began in early December of 2020 in a few countries. There are several vaccines, and each has a different efficacy and mechanism of action. Several countries, for example, the United Kingdom and the USA, have been able to develop consistent vaccination programs where a great percentage of the population has been vaccinated (May 2021). However, in other countries, a low percentage of the population has been vaccinated due to constraints related to vaccine supply and distribution capacity. Countries such as the USA and the UK have implemented different vaccination strategies, and some scholars have been debating the optimal strategy for vaccine campaigns. This problem is complex due to the great number of variables that affect the relevant outcomes. In this article, we study the impact of different vaccination regimens on main health outcomes such as deaths, hospitalizations, and the number of infected. We develop a mathematical model of COVID-19 transmission to focus on this important health policy issue. Thus, we are able to identify the optimal strategy regarding vaccination campaigns. We find that for vaccines with high efficacy (>70%) after the first dose, the optimal strategy is to delay inoculation with the second dose. On the other hand, for a low first dose vaccine efficacy, it is better to use the standard vaccination regimen of 4 weeks between doses. Thus, under the delayed second dose option, a campaign focus on generating a certain immunity in as great a number of people as fast as possible is preferable to having an almost perfect immunity in fewer people first. Therefore, based on these results, we suggest that the UK implemented a better vaccination campaign than that in the USA with regard to time between doses. The results presented here provide scientific guidelines for other countries where vaccination campaigns are just starting, or the percentage of vaccinated people is small.

Dynamics of SARS-CoV-2 with waning immunity in the UK population

The dynamics of immunity are crucial to understanding the long-term patterns of the SARS-CoV-2 pandemic. Several cases of reinfection with SARS-CoV-2 have been documented 48-142 days after the initial infection and immunity to seasonal circulating coronaviruses is estimated to be shorter than 1 year. Using an age-structured, deterministic model, we explore potential immunity dynamics using contact data from the UK population. In the scenario where immunity to SARS-CoV-2 lasts an average of three months for non-hospitalized individuals, a year for hospitalized individuals, and the effective reproduction number after lockdown ends is 1.2 (our worst-case scenario), we find that the secondary peak occurs in winter 2020 with a daily maximum of 387 000 infectious individuals and 125 000 daily new cases; threefold greater than in a scenario with permanent immunity. Our models suggest that longitudinal serological surveys to determine if immunity in the population is waning will be most informative when sampling takes place from the end of the lockdown in June until autumn 2020. After this period, the proportion of the population with antibodies to SARS-CoV-2 is expected to increase due to the secondary wave. Overall, our analysis presents considerations for policy makers on the longer-term dynamics of SARS-CoV-2 in the UK and suggests that strategies designed to achieve herd immunity may lead to repeated waves of infection as immunity to reinfection is not permanent. This article is part of the theme issue 'Modelling that shaped the early COVID-19 pandemic response in the UK'.

Molecular Evolution and Epidemiological Characteristics of SARS COV-2 in (Northwestern) Poland

The emergence of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) evolved into a worldwide outbreak, with the first Polish cases in February/March 2020. This study aimed to investigate the molecular epidemiology of the circulating virus lineages between March 2020 and February 2021. We performed variant identification, spike mutation pattern analysis, and phylogenetic and evolutionary analyses for 1106 high-coverage whole-genome sequences, implementing maximum likelihood, multiple continuous-time Markov chain, and Bayesian birth-death skyline models. For time trends, logistic regression was used. In the dataset, virus B.1.221 lineage was predominant (15.37%), followed by B.1.258 (15.01%) and B.1.1.29 (11.48%) strains. Three clades were identified, being responsible for 74.41% of infections over the analyzed period. Expansion in variant diversity was observed since September 2020 with increasing frequency of the number in spike substitutions, mainly H69V70 deletion, P681H, N439K, and S98F. In population dynamics inferences, three periods with exponential increase in infection were observed, beginning in March, July, and September 2020, respectively, and were driven by different virus clades. Additionally, a notable increase in infections caused by the B.1.1.7 lineage since February 2021 was noted. Over time, the virus accumulated mutations related to optimized transmissibility; therefore, faster dissemination is reflected by the second wave of epidemics in Poland.

COVID-19 Community Temporal Visualizer: a new methodology for the network-based analysis and visualization of COVID-19 data

Understanding the evolution of the spread of the COVID-19 pandemic requires the analysis of several data at the spatial and temporal levels. Here, we present a new network-based methodology to analyze COVID-19 data measures containing spatial and temporal features and its application on a real dataset. The goal of the methodology is to analyze sets of homogeneous datasets (i.e. COVID-19 data taken in different periods and in several regions) using a statistical test to find similar/dissimilar datasets, mapping such similarity information on a graph and then using a community detection algorithm to visualize and analyze the spatio-temporal evolution of data. We evaluated diverse Italian COVID-19 data made publicly available by the Italian Protezione Civile Department at https://github.com/pcm-dpc/COVID-19/. Furthermore, we considered the climate data related to two periods and we integrated them with COVID-19 data measures to detect new communities related to climate changes. In conclusion, the application of the proposed methodology provides a network-based representation of the COVID-19 measures by highlighting the different behaviour of regions with respect to pandemics data released by Protezione Civile and climate data. The methodology and its implementation as R function are publicly available at https://github.com/mmilano87/analyzeC19D.

The role of disinfectants and sanitizers during COVID-19 pandemic: advantages and deleterious effects on humans and the environment

Disinfectants and sanitizers are essential preventive agents against the coronavirus disease 2019 (COVID-19) pandemic; however, the pandemic crisis was marred by undue hype, which led to the indiscriminate use of disinfectants and sanitizers. Despite demonstrating a beneficial role in the control and prevention of COVID-19, there are crucial concerns regarding the large-scale use of disinfectants and sanitizers, including the side effects on human and animal health along with harmful impacts exerted on the environment and ecological balance. This article discusses the roles of disinfectants and sanitizers in the control and prevention of the current pandemic and highlights updated disinfection techniques against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This article provides evidence of the deleterious effects of disinfectants and sanitizers exerted on humans, animals, and the environment as well as suggests mitigation strategies to reduce these effects. Additionally, potential technologies and approaches for the reduction of these effects and the development of safe, affordable, and effective disinfectants are discussed, particularly, eco-friendly technologies using nanotechnology and nanomedicine.

Phylodynamics reveals the role of human travel and contact tracing in controlling the first wave of COVID-19 in four island nations

New Zealand, Australia, Iceland, and Taiwan all saw success in controlling their first waves of Coronavirus Disease 2019 (COVID-19). As islands, they make excellent case studies for exploring the effects of international travel and human movement on the spread of COVID-19. We employed a range of robust phylodynamic methods and genome subsampling strategies to infer the epidemiological history of Severe acute respiratory syndrome coronavirus 2 in these four countries. We compared these results to transmission clusters identified by the New Zealand Ministry of Health by contact tracing strategies. We estimated the effective reproduction number of COVID-19 as 1-1.4 during early stages of the pandemic and show that it declined below 1 as human movement was restricted. We also showed that this disease was introduced many times into each country and that introductions slowed down markedly following the reduction of international travel in mid-March 2020. Finally, we confirmed that New Zealand transmission clusters identified via standard health surveillance strategies largely agree with those defined by genomic data. We have demonstrated how the use of genomic data and computational biology methods can assist health officials in characterising the epidemiology of viral epidemics and for contact tracing.

Analysis of Key Factors of a SARS-CoV-2 Vaccination Program: A Mathematical Modeling Approach

The administration of vaccines against the coronavirus disease 2019 (COVID-19) started in early December of 2020. Currently, there are only a few approved vaccines, each with different efficacies and mechanisms of action. Moreover, vaccination programs in different regions may vary due to differences in implementation, for instance, simply the availability of the vaccine. In this article, we study the impact of the pace of vaccination and the intrinsic efficacy of the vaccine on prevalence, hospitalizations, and deaths related to the SARS-CoV-2 virus. Then we study different potential scenarios regarding the burden of the COVID-19 pandemic in the near future. We construct a compartmental mathematical model and use computational methodologies to study these different scenarios. Thus, we are able to identify some key factors to reach the aims of the vaccination programs. We use some metrics related to the outcomes of the COVID-19 pandemic in order to assess the impact of the efficacy of the vaccine and the pace of the vaccine inoculation. We found that both factors have a high impact on the outcomes. However, the rate of vaccine administration has a higher impact in reducing the burden of the COVID-19 pandemic. This result shows that health institutions need to focus on increasing the vaccine inoculation pace and create awareness in the population about the importance of COVID-19 vaccines.

Clinical features of COVID-19 and SARS epidemics. A literature review

SARS-CoV-2, responsible for the current pandemic, is a novel strain of the Coronaviridae family, which has infected humans as a result of the leap to a new species. It causes an atypical pneumonia similar to that caused by SARS-CoV in 2003. SARS-CoV-2 has currently infected more than 9,200,000 people and caused almost 480,000 deaths worldwide. Although SARS-CoV-2 and SARS-CoV have similar phylogenetic and pathogenetic characteristics, they show important differences in clinical manifestations.

We have reviewed the recent literature comparing the characteristics of the two epidemics and highlight their peculiar aspects. An analysis of all signs and symptoms of 3,365 SARS patients and 23,280 COVID-19 patients as well as of the comorbidities has been carried out. A total of 17 and 75 studies regarding patients with SARS and COVID-19, respectively, were included in the analysis. The analysis revealed an overlap of some symptoms between the two infections. Unlike SARS patients, COVID-19 patients have developed respiratory, neurological and gastrointestinal symptoms, and, in a limited number of subjects, symptoms involving organs such as skin and subcutaneous tissue, kidneys, cardiovascular system, liver and eyes.

This analysis was conducted in order to direct towards an early identification of the infection, a suitable diagnostic procedure and the adoption of appropriate containment measures.

SARS-CoV-2 transmission dynamics in Belarus revealed by genomic and incidence data analysis

Since the emergence of COVID-19, a series of non-pharmaceutical interventions (NPIs) has been implemented by governments and public health authorities world-wide to control and curb the ongoing pandemic spread. From that perspective, Belarus is one of a few countries with a relatively modern healthcare system, where much narrower NPIs have been put in place. Given the uniqueness of this Belarusian experience, the understanding its COVID-19 epidemiological dynamics is essential not only for the local assessment, but also for a better insight into the impact of different NPI strategies globally. In this work, we integrate genomic epidemiology and surveillance methods to investigate the emergence and spread of SARS-CoV-2 in the country. The observed Belarusian SARS-CoV-2 genetic diversity originated from at least eighteen separate introductions, at least five of which resulted in on-going domestic transmissions. The introduction sources represent a wide variety of regions, although the proportion of regional virus introductions and exports from/to geographical neighbors appears to be higher than for other European countries. Phylodynamic analysis indicates a moderate reduction in the effective reproductive number ℛ e after the introduction of limited NPIs, with the reduction magnitude generally being lower than for countries with large-scale NPIs. On the other hand, the estimate of the Belarusian ℛ e at the early epidemic stage is comparable with this number for the neighboring ex-USSR country of Ukraine, where much broader NPIs have been implemented. The actual number of cases by the end of May, 2020 was predicted to be 2-9 times higher than the detected number of cases.