Emerging phylogenetic structure of the SARS-CoV-2 pandemic

Prediction of Prospective Mutational Landscape of SARS-CoV-2 Spike ssRNA and Evolutionary Basis of Its Host Interaction

Booster doses are crucial against severe COVID-19, as rapid virus mutations and variant emergence prolong the pandemic crisis. The virus's quick evolution, short generation-time, and adaptive changes impact virulence and evolvability, helping predictions about variant of concerns' (VOCs') landscapes. Here, in this study, we used a new computational algorithm, to predict the mutational pattern in SARS-CoV-2 ssRNA, proteomics, structural identification, mutation stability, and functional correlation, as well as immune escape mechanisms. Interestingly, the sequence diversity of SARS Coronavirus-2 has demonstrated a predominance of G- > A and C- > U substitutions. The best validation statistics are explored here in seven homologous models of the expected mutant SARS-CoV-2 spike ssRNA and employed for hACE2 and IgG interactions. The interactome profile of SARS-CoV-2 spike with hACE2 and IgG revealed a strong correlation between phylogeny and divergence time. The systematic adaptation of SARS-CoV-2 spike ssRNA influences infectivity and immune escape. Data suggest higher propensity of Adenine rich sequence promotes MHC system avoidance, preferred by A-rich codons. Phylogenetic data revealed the evolution of SARS-CoV-2 lineages' epidemiology. Our findings may unveil processes governing the genesis of immune-resistant variants, prompting a critical reassessment of the coronavirus mutation rate and exploration of hypotheses beyond mechanical aspects.

Identification of differences in the magnitude and specificity of SARS-CoV-2 nucleocapsid antibody responses in naturally infected and vaccinated individuals

As there are limited data on B-cell epitopes for the nucleocapsid protein in SARS-CoV-2, we sought to identify the immunodominant regions within the N protein, recognized by patients with varying severity of natural infection with the Wuhan strain (WT), delta, omicron, and in those who received the Sinopharm vaccines, which is an inactivated, whole virus vaccine. Using overlapping peptides representing the N protein, with an in-house ELISA, we mapped the immunodominant regions within the N protein, in seronegative (n = 30), WT infected (n = 30), delta infected (n = 30), omicron infected + vaccinated (n = 20) and Sinopharm (BBIBP-CorV) vaccinees (n = 30). We then investigated the sensitivity and specificity of these immunodominant regions and analyzed their conservation with other SARS-CoV-2 variants of concern, seasonal human coronaviruses, and bat Sarbecoviruses. We identified four immunodominant regions aa 29-52, aa 155-178, aa 274-297, and aa 365-388, which were highly conserved within SARS-CoV-2 and the bat coronaviruses. The magnitude of responses to these regions varied based on the infecting SARS-CoV-2 variants, >80% of individuals gave responses above the positive cut-off threshold to many of the four regions, with some differences with individuals who were infected with different VoCs. These regions were found to be 100% specific, as none of the seronegative individuals gave any responses. As these regions were highly specific with high sensitivity, they have a potential to be used to develop diagnostic assays and to be used in development of vaccines.

Fitness, growth and transmissibility of SARS-CoV-2 genetic variants

The massive scale of the global SARS-CoV-2 sequencing effort created new opportunities and challenges for understanding SARS-CoV-2 evolution. Rapid detection and assessment of new variants has become one of the principal objectives of genomic surveillance of SARS-CoV-2. Because of the pace and scale of sequencing, new strategies have been developed for characterizing fitness and transmissibility of emerging variants. In this Review, I discuss a wide range of approaches that have been rapidly developed in response to the public health threat posed by emerging variants, ranging from new applications of classic population genetics models to contemporary synthesis of epidemiological models and phylodynamic analysis. Many of these approaches can be adapted to other pathogens and will have increasing relevance as large-scale pathogen sequencing becomes a regular feature of many public health systems.

The tumour is in the detail: Local phylogenetic, population and epidemiological dynamics of a transmissible cancer in Tasmanian devils

Infectious diseases are a major threat for biodiversity conservation and can exert strong influence on wildlife population dynamics. Understanding the mechanisms driving infection rates and epidemic outcomes requires empirical data on the evolutionary trajectory of pathogens and host selective processes. Phylodynamics is a robust framework to understand the interaction of pathogen evolutionary processes with epidemiological dynamics, providing a powerful tool to evaluate disease control strategies. Tasmanian devils have been threatened by a fatal transmissible cancer, devil facial tumour disease (DFTD), for more than two decades. Here we employ a phylodynamic approach using tumour mitochondrial genomes to assess the role of tumour genetic diversity in epidemiological and population dynamics in a devil population subject to 12 years of intensive monitoring, since the beginning of the epidemic outbreak. DFTD molecular clock estimates of disease introduction mirrored observed estimates in the field, and DFTD genetic diversity was positively correlated with estimates of devil population size. However, prevalence and force of infection were the lowest when devil population size and tumour genetic diversity was the highest. This could be due to either differential virulence or transmissibility in tumour lineages or the development of host defence strategies against infection. Our results support the view that evolutionary processes and epidemiological trade-offs can drive host-pathogen coexistence, even when disease-induced mortality is extremely high. We highlight the importance of integrating pathogen and population evolutionary interactions to better understand long-term epidemic dynamics and evaluating disease control strategies.

Identification of differences in the magnitude and specificity of SARS-CoV-2 nucleocapsid antibody responses in naturally infected and vaccinated individuals

Background

As there are limited data on B cell epitopes for the nucleocapsid protein in SARS-CoV-2, we sought to identify the immunodominant regions within the N protein, recognized by patients with varying severity of natural infection with the Wuhan strain (WT), delta, omicron and in those who received the Sinopharm vaccines, which is an inactivated, whole virus vaccine.

Methods

Using overlapping peptides representing the N protein, with an in-house ELISA, we mapped the immunodominant regions within the N protein, in seronegative (n=30), WT infected (n=30), delta infected (n=30), omicron infected+vaccinated (n=20) and Sinopharm (BBIBP-CorV) vaccinees (n=30). We then investigated the sensitivity and specificity of these immunodominant regions and analysed their conservation with other SARS-CoV-2 variants of concern, seasonal human coronaviruses and bat Sarbecoviruses. We then investigated the kinetics of responses to these regions in those with varying severity of acute COVID-19.

Results

We identified four immunodominant regions aa 29-52, aa 155-178, aa 274 to 297 and aa 365 to 388, were highly conserved within SARS-CoV-2 and the bat coronaviruses. The magnitude of responses to these regions varied based on the infecting SARS-CoV-2 variants, with WT infected individuals predominantly recognizing aa155 to 178 regions, delta infected individuals and vaccinated+omicron infected individuals predominantly recognizing regions aa 29 to 52 and aa 274 to 294 regions. Sinopharm vaccinees recognized all four regions, with the magnitude of responses significantly lower than other groups. >80% of individuals gave responses above the positive cut-off threshold to many of the four regions, with some differences with individuals who were infected with different VoCs. These regions were found to be 100% specific, as none of the seronegative individuals gave any responses.

Conclusions

N-protein specific responses appear to be detectable in over 90% of those who were naturally infected or vaccinated with a whole virus inactivated vaccine, with responses mainly directed against four regions of the protein, which were highly conserved. As these regions were highly specific with high sensitivity, they have a potential to be used to develop diagnostic assays and to be used in development of vaccines.

Phylogenetic and phylodynamic approaches to understanding and combating the early SARS-CoV-2 pandemic

Determining the transmissibility, prevalence and patterns of movement of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is central to our understanding of the impact of the pandemic and to the design of effective control strategies. Phylogenies (evolutionary trees) have provided key insights into the international spread of SARS-CoV-2 and enabled investigation of individual outbreaks and transmission chains in specific settings. Phylodynamic approaches combine evolutionary, demographic and epidemiological concepts and have helped track virus genetic changes, identify emerging variants and inform public health strategy. Here, we review and synthesize studies that illustrate how phylogenetic and phylodynamic techniques were applied during the first year of the pandemic, and summarize their contributions to our understanding of SARS-CoV-2 transmission and control.

Comparative Phylodynamics Reveals the Evolutionary History of SARS-CoV-2 Emerging Variants in the Arabian Peninsula

Emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants continue to be responsible for an unprecedented worldwide public health and economic catastrophe. Accurate understanding and comparison of global and regional evolutionary epidemiology of novel SARS-CoV-2 variants are critical to guide current and future interventions. Here, we utilized a Bayesian phylodynamic pipeline to trace and compare the evolutionary dynamics, spatiotemporal origins, and spread of five variants (Alpha, Beta, Delta, Kappa, and Eta) across the Arabian Peninsula. We found variant-specific signatures of evolution and spread that are likely linked to air travel and disease control interventions in the region. Alpha, Beta, and Delta variants went through sequential periods of growth and decline, whereas we inferred inconclusive population growth patterns for the Kappa and Eta variants due to their sporadic introductions in the region. Non-pharmaceutical interventions imposed between mid-2020 and early 2021 likely played a role in reducing the epidemic progression of the Beta and the Alpha variants. In comparison, the combination of the non-pharmaceutical interventions and the rapid rollout of vaccination might have shaped Delta variant dynamics. We found that the Alpha and Beta variants were frequently introduced into the Arab peninsula between mid-2020 and early 2021 from Europe and Africa, respectively, whereas the Delta variant was frequently introduced between early 2021 and mid-2021 from East Asia. For these three variants, we also revealed significant and intense dispersal routes between the Arab region and Africa, Europe, Asia, and Oceania. In contrast, the restricted spread and stable effective population size of the Kappa and the Eta variants suggest that they no longer need to be targeted in genomic surveillance activities in the region. In contrast, the evolutionary characteristics of the Alpha, Beta, and Delta variants confirm the dominance of these variants in the recent outbreaks. Our study highlights the urgent need to establish regional molecular surveillance programs to ensure effective decision making related to the allocation of intervention activities targeted toward the most relevant variants.

On the Origin and Propagation of the COVID-19 Outbreak in the Italian Province of Trento, a Tourist Region of Northern Italy

Background: Trentino is an Italian province with a tourism-based economy, bordering the regions of Lombardy and Veneto, where the two earliest and largest outbreaks of COVID-19 occurred in Italy. The earliest cases in Trentino were reported in the first week of March 2020, with most of the cases occurring in the winter sport areas in the Dolomites mountain range. The number of reported cases decreased over the summer months and was followed by a second wave in the autumn and winter of 2020. Methods: we performed high-coverage Oxford Nanopore sequencing of 253 positive SARS-CoV-2 swabs collected in Trentino between March and December 2020. Results: in this work, we analyzed genome sequences to trace the routes through which the virus entered the area, and assessed whether the autumnal resurgence could be attributed to lineages persisting undetected during summer, or as a consequence of new introductions. Conclusions: Comparing the draft genomes analyzed with a large selection of European sequences retrieved from GISAID we found that multiple introductions of the virus occurred at the early stage of the epidemics; the two epidemic waves were unrelated; the second wave was due to reintroductions of the virus in summer when traveling restrictions were uplifted.

Hunting alters viral transmission and evolution in a large carnivore

Hunting can fundamentally alter wildlife population dynamics but the consequences of hunting on pathogen transmission and evolution remain poorly understood. Here, we present a study that leverages a unique landscape-scale quasi-experiment coupled with pathogen-transmission tracing, network simulation and phylodynamics to provide insights into how hunting shapes feline immunodeficiency virus (FIV) dynamics in puma (Puma concolor). We show that removing hunting pressure enhances the role of males in transmission, increases the viral population growth rate and increases the role of evolutionary forces on the pathogen compared to when hunting was reinstated. Changes in transmission observed with the removal of hunting could be linked to short-term social changes while the male puma population increased. These findings are supported through comparison with a region with stable hunting management over the same time period. This study shows that routine wildlife management can have impacts on pathogen transmission and evolution not previously considered.

Genomic and Epidemiological Analysis of SARS-CoV-2 Viruses in Sri Lanka

Background: In order to understand the molecular epidemiology of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) in Sri Lanka, since March 2020, we carried out genomic sequencing overlaid on available epidemiological data until April 2021.

Methods: Whole genome sequencing was carried out on diagnostic sputum or nasopharyngeal swabs from 373 patients with COVID-19. Molecular clock phylogenetic analysis was undertaken to further explore dominant lineages.

Results: The B.1.411 lineage was most prevalent, which was established in Sri Lanka and caused outbreaks throughout the country until March 2021. The estimated time of the most recent common ancestor (tMRCA) of this lineage was June 1, 2020 (with 95% lower and upper bounds March 30 to July 27) suggesting cryptic transmission may have occurred, prior to a large epidemic starting in October 2020. Returning travellers were identified with infections caused by lineage B.1.258, as well as the more transmissible B.1.1.7 lineage, which has replaced B.1.411 to fuel the ongoing large outbreak in the country.

Conclusions: The large outbreak that started in early October, is due to spread of a single virus lineage, B.1.411 until the end of March 2021, when B.1.1.7 emerged and became the dominant lineage.

Emerging SARS-CoV-2 Diversity Revealed by Rapid Whole-Genome Sequence Typing

Discrete classification of SARS-CoV-2 viral genotypes can identify emerging strains and detect geographic spread, viral diversity, and transmission events. We developed a tool (GNU-based Virus IDentification [GNUVID]) that integrates whole-genome multilocus sequence typing and a supervised machine learning random forest-based classifier. We used GNUVID to assign sequence type (ST) profiles to all high-quality genomes available from GISAID. STs were clustered into clonal complexes (CCs) and then used to train a machine learning classifier. We used this tool to detect potential introduction and exportation events and to estimate effective viral diversity across locations and over time in 16 US states. GNUVID is a highly scalable tool for viral genotype classification (https://github.com/ahmedmagds/GNUVID) that can quickly classify hundreds of thousands of genomes in a way that is consistent with phylogeny. Our genotyping ST/CC analysis uncovered dynamic local changes in ST/CC prevalence and diversity with multiple replacement events in different states, an average of 20.6 putative introductions and 7.5 exportations for each state over the time period analyzed. We introduce the use of effective diversity metrics (Hill numbers) that can be used to estimate the impact of interventions (e.g., travel restrictions, vaccine uptake, mask mandates) on the variation in circulating viruses. Our classification tool uncovered multiple introduction and exportation events, as well as waves of expansion and replacement of SARS-CoV-2 genotypes in different states. GNUVID classification lends itself to measures of ecological diversity, and, with systematic genomic sampling, it could be used to track circulating viral diversity and identify emerging clones and hotspots.

Ongoing global and regional adaptive evolution of SARS-CoV-2

Understanding the trends in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) evolution is paramount to control the COVID-19 pandemic. We analyzed more than 300,000 high-quality genome sequences of SARS-CoV-2 variants available as of January 2021. The results show that the ongoing evolution of SARS-CoV-2 during the pandemic is characterized primarily by purifying selection, but a small set of sites appear to evolve under positive selection. The receptor-binding domain of the spike protein and the region of the nucleocapsid protein associated with nuclear localization signals (NLS) are enriched with positively selected amino acid replacements. These replacements form a strongly connected network of apparent epistatic interactions and are signatures of major partitions in the SARS-CoV-2 phylogeny. Virus diversity within each geographic region has been steadily growing for the entirety of the pandemic, but analysis of the phylogenetic distances between pairs of regions reveals four distinct periods based on global partitioning of the tree and the emergence of key mutations. The initial period of rapid diversification into region-specific phylogenies that ended in February 2020 was followed by a major extinction event and global homogenization concomitant with the spread of D614G in the spike protein, ending in March 2020. The NLS-associated variants across multiple partitions rose to global prominence in March to July, during a period of stasis in terms of interregional diversity. Finally, beginning in July 2020, multiple mutations, some of which have since been demonstrated to enable antibody evasion, began to emerge associated with ongoing regional diversification, which might be indicative of speciation.

The roles of environmental variation and parasite survival in virulence-transmission relationships

Disease outbreaks are a consequence of interactions among the three components of a host-parasite system: the infectious agent, the host and the environment. While virulence and transmission are widely investigated, most studies of parasite life-history trade-offs are conducted with theoretical models or tractable experimental systems where transmission is standardized and the environment controlled. Yet, biotic and abiotic environmental factors can strongly affect disease dynamics, and ultimately, host-parasite coevolution. Here, we review research on how environmental context alters virulence-transmission relationships, focusing on the off-host portion of the parasite life cycle, and how variation in parasite survival affects the evolution of virulence and transmission. We review three inter-related 'approaches' that have dominated the study of the evolution of virulence and transmission for different host-parasite systems: (i) evolutionary trade-off theory, (ii) parasite local adaptation and (iii) parasite phylodynamics. These approaches consider the role of the environment in virulence and transmission evolution from different angles, which entail different advantages and potential biases. We suggest improvements to how to investigate virulence-transmission relationships, through conceptual and methodological developments and taking environmental context into consideration. By combining developments in life-history evolution, phylogenetics, adaptive dynamics and comparative genomics, we can improve our understanding of virulence-transmission relationships across a diversity of host-parasite systems that have eluded experimental study of parasite life history.

Nano-Enabled COVID-19 Vaccines: Meeting the Challenges of Durable Antibody Plus Cellular Immunity and Immune Escape

At the time of preparing this Perspective, large-scale vaccination for COVID-19 is in progress, aiming to bring the pandemic under control through vaccine-induced herd immunity. Not only does this vaccination effort represent an unprecedented scientific and technological breakthrough, moving us from the rapid analysis of viral genomes to design, manufacture, clinical trial testing, and use authorization within the time frame of less than a year, but it also highlights rapid progress in the implementation of nanotechnology to assist vaccine development. These advances enable us to deliver nucleic acid and conformation-stabilized subunit vaccines to regional lymph nodes, with the ability to trigger effective humoral and cellular immunity that prevents viral infection or controls disease severity. In addition to a brief description of the design features of unique cationic lipid and virus-mimicking nanoparticles for accomplishing spike protein delivery and presentation by the cognate immune system, we also discuss the importance of adjuvancy and design features to promote cooperative B- and T-cell interactions in lymph node germinal centers, including the use of epitope-based vaccines. Although current vaccine efforts have demonstrated short-term efficacy and vaccine safety, key issues are now vaccine durability and adaptability against viral variants. We present a forward-looking perspective of how vaccine design can be adapted to improve durability of the immune response and vaccine adaptation to overcome immune escape by viral variants. Finally, we consider the impact of nano-enabled approaches in the development of COVID-19 vaccines for improved vaccine design against other infectious agents, including pathogens that may lead to future pandemics.

Genetic evidence for the association between COVID-19 epidemic severity and timing of non-pharmaceutical interventions

Unprecedented public health interventions including travel restrictions and national lockdowns have been implemented to stem the COVID-19 epidemic, but the effectiveness of non-pharmaceutical interventions is still debated. We carried out a phylogenetic analysis of more than 29,000 publicly available whole genome SARS-CoV-2 sequences from 57 locations to estimate the time that the epidemic originated in different places. These estimates were examined in relation to the dates of the most stringent interventions in each location as well as to the number of cumulative COVID-19 deaths and phylodynamic estimates of epidemic size. Here we report that the time elapsed between epidemic origin and maximum intervention is associated with different measures of epidemic severity and explains 11% of the variance in reported deaths one month after the most stringent intervention. Locations where strong non-pharmaceutical interventions were implemented earlier experienced much less severe COVID-19 morbidity and mortality during the period of study.

Genomic epidemiology of SARS-CoV-2 transmission lineages in Ecuador

Characterisation of SARS-CoV-2 genetic diversity through space and time can reveal trends in virus importation and domestic circulation, and permit the exploration of questions regarding the early transmission dynamics. Here we present a detailed description of SARS-CoV-2 genomic epidemiology in Ecuador, one of the hardest hit countries during the early stages of the COVID-19 pandemic. We generate and analyse 160 whole genome sequences sampled from all provinces of Ecuador in 2020. Molecular clock and phylgeographic analysis of these sequences in the context of global SARS-CoV-2 diversity enable us to identify and characterise individual transmission lineages within Ecuador, explore their spatiotemporal distributions, and consider their introduction and domestic circulation. Our results reveal a pattern of multiple international importations across the country, with apparent differences between key provinces. Transmission lineages were mostly introduced before the implementation of non-pharmaceutical interventions (NPIs), with differential degrees of persistence and national dissemination.

Ongoing Global and Regional Adaptive Evolution of SARS-CoV-2

Understanding the trends in SARS-CoV-2 evolution is paramount to control the COVID-19 pandemic. We analyzed more than 300,000 high quality genome sequences of SARS-CoV-2 variants available as of January 2021. The results show that the ongoing evolution of SARS-CoV-2 during the pandemic is characterized primarily by purifying selection, but a small set of sites appear to evolve under positive selection. The receptor-binding domain of the spike protein and the nuclear localization signal (NLS) associated region of the nucleocapsid protein are enriched with positively selected amino acid replacements. These replacements form a strongly connected network of apparent epistatic interactions and are signatures of major partitions in the SARS-CoV-2 phylogeny. Virus diversity within each geographic region has been steadily growing for the entirety of the pandemic, but analysis of the phylogenetic distances between pairs of regions reveals four distinct periods based on global partitioning of the tree and the emergence of key mutations. The initial period of rapid diversification into region-specific phylogenies that ended in February 2020 was followed by a major extinction event and global homogenization concomitant with the spread of D614G in the spike protein, ending in March 2020. The NLS associated variants across multiple partitions rose to global prominence in March-July, during a period of stasis in terms of inter-regional diversity. Finally, beginning July 2020, multiple mutations, some of which have since been demonstrated to enable antibody evasion, began to emerge associated with ongoing regional diversification, which might be indicative of speciation.

SARS-CoV-2 Infection and Disease Modelling Using Stem Cell Technology and Organoids

In this Review, we briefly describe the basic virology and pathogenesis of SARS-CoV-2, highlighting how stem cell technology and organoids can contribute to the understanding of SARS-CoV-2 cell tropisms and the mechanism of disease in the human host, supporting and clarifying findings from clinical studies in infected individuals. We summarize here the results of studies, which used these technologies to investigate SARS-CoV-2 pathogenesis in different organs. Studies with in vitro models of lung epithelia showed that alveolar epithelial type II cells, but not differentiated lung alveolar epithelial type I cells, are key targets of SARS-CoV-2, which triggers cell apoptosis and inflammation, while impairing surfactant production. Experiments with human small intestinal organoids and colonic organoids showed that the gastrointestinal tract is another relevant target for SARS-CoV-2. The virus can infect and replicate in enterocytes and cholangiocytes, inducing cell damage and inflammation. Direct viral damage was also demonstrated in in vitro models of human cardiomyocytes and choroid plexus epithelial cells. At variance, endothelial cells and neurons are poorly susceptible to viral infection, thus supporting the hypothesis that neurological symptoms and vascular damage result from the indirect effects of systemic inflammatory and immunological hyper-responses to SARS-CoV-2 infection.

Genomic epidemiology of SARS-CoV-2 transmission lineages in Ecuador

Characterisation of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) genetic diversity through space and time can reveal trends in virus importation and domestic circulation and permit the exploration of questions regarding the early transmission dynamics. Here, we present a detailed description of SARS-CoV-2 genomic epidemiology in Ecuador, one of the hardest hit countries during the early stages of the coronavirus-19 pandemic. We generated and analysed 160 whole genome sequences sampled from all provinces of Ecuador in 2020. Molecular clock and phylogeographic analysis of these sequences in the context of global SARS-CoV-2 diversity enable us to identify and characterise individual transmission lineages within Ecuador, explore their spatiotemporal distributions, and consider their introduction and domestic circulation. Our results reveal a pattern of multiple international importations across the country, with apparent differences between key provinces. Transmission lineages were mostly introduced before the implementation of non-pharmaceutical interventions, with differential degrees of persistence and national dissemination.

No evidence for increased transmissibility from recurrent mutations in SARS-CoV-2

COVID-19 is caused by the coronavirus SARS-CoV-2, which jumped into the human population in late 2019 from a currently uncharacterised animal reservoir. Due to this recent association with humans, SARS-CoV-2 may not yet be fully adapted to its human host. This has led to speculations that SARS-CoV-2 may be evolving towards higher transmissibility. The most plausible mutations under putative natural selection are those which have emerged repeatedly and independently (homoplasies). Here, we formally test whether any homoplasies observed in SARS-CoV-2 to date are significantly associated with increased viral transmission. To do so, we develop a phylogenetic index to quantify the relative number of descendants in sister clades with and without a specific allele. We apply this index to a curated set of recurrent mutations identified within a dataset of 46,723 SARS-CoV-2 genomes isolated from patients worldwide. We do not identify a single recurrent mutation in this set convincingly associated with increased viral transmission. Instead, recurrent mutations currently in circulation appear to be evolutionary neutral and primarily induced by the human immune system via RNA editing, rather than being signatures of adaptation. At this stage we find no evidence for significantly more transmissible lineages of SARS-CoV-2 due to recurrent mutations.