SARS-CoV-2 exhibits intra-host genomic plasticity and low-frequency polymorphic quasispecies

Point mutations at specific sites of the nsp12-nsp8 interface dramatically affect the RNA polymerization activity of SARS-CoV-2

In a recent characterization of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variability present in 30 diagnostic samples from patients of the first COVID-19 pandemic wave, 41 amino acid substitutions were documented in the RNA-dependent RNA polymerase (RdRp) nsp12. Eight substitutions were selected in this work to determine whether they had an impact on the RdRp activity of the SARS-CoV-2 nsp12-nsp8-nsp7 replication complex. Three of these substitutions were found around the polymerase central cavity, in the template entry channel (D499G and M668V), and within the motif B (V560A), and they showed polymerization rates similar to the wild type RdRp. The remaining five mutations (P323L, L372F, L372P, V373A, and L527H) were placed near the nsp12-nsp8F contact surface; residues L372, V373, and L527 participated in a large hydrophobic cluster involving contacts between two helices in the nsp12 fingers and the long α-helix of nsp8F. The presence of any of these five amino acid substitutions resulted in important alterations in the RNA polymerization activity. Comparative primer elongation assays showed different behavior depending on the hydrophobicity of their side chains. The substitution of L by the bulkier F side chain at position 372 slightly promoted RdRp activity. However, this activity was dramatically reduced with the L372P, and L527H mutations, and to a lesser extent with V373A, all of which weaken the hydrophobic interactions within the cluster. Additional mutations, specifically designed to disrupt the nsp12-nsp8F interactions (nsp12-V330S, nsp12-V341S, and nsp8-R111A/D112A), also resulted in an impaired RdRp activity, further illustrating the importance of this contact interface in the regulation of RNA synthesis.

The host-targeted antiviral drug Zapnometinib exhibits a high barrier to the development of SARS-CoV-2 resistance

Host targeting antiviral drugs (HTA) are directed against cellular mechanisms which can be exploited by viruses. These mechanisms are essential for viral replication, because missing functions cannot be compensated by the virus. However, this assumption needs experimental proof. Here we compared the HTA Zapnometinib (ZMN), with direct acting antivirals (DAA) (Remdesivir (RDV), Molnupiravir (MPV), Nirmatrelvir (NTV), Ritonavir (RTV), Paxlovid PAX)), in terms of their potency to induce reduced drug susceptibilities in SARS-CoV-2. During serial passage of δ-B1.617.2 adaptation to all DAAs occurred, while the inhibitory capacity of ZMN was not altered. Known single nucleotide polymorphisms (SNPs) responsible for partial resistances were found for RDV, NTV and PAX. Additionally, the high mutagenic potential of MPV was confirmed and decreased drug efficacies were found for the first time. Reduced DAA efficacy did not alter the inhibitory potential of ZMN. These results show that ZMN confers a high barrier towards the development of viral resistance and has the potential to act against partially DAA-insensitive viruses.

Incipient functional SARS-CoV-2 diversification identified through neural network haplotype maps

Since its introduction in the human population, SARS-CoV-2 has evolved into multiple clades, but the events in its intrahost diversification are not well understood. Here, we compare three-dimensional (3D) self-organized neural haplotype maps (SOMs) of SARS-CoV-2 from thirty individual nasopharyngeal diagnostic samples obtained within a 19-day interval in Madrid (Spain), at the time of transition between clades 19 and 20. SOMs have been trained with the haplotype repertoire present in the mutant spectra of the nsp12- and spike (S)-coding regions. Each SOM consisted of a dominant neuron (displaying the maximum frequency), surrounded by a low-frequency neuron cloud. The sequence of the master (dominant) neuron was either identical to that of the reference Wuhan-Hu-1 genome or differed from it at one nucleotide position. Six different deviant haplotype sequences were identified among the master neurons. Some of the substitutions in the neural clouds affected critical sites of the nsp12-nsp8-nsp7 polymerase complex and resulted in altered kinetics of RNA synthesis in an in vitro primer extension assay. Thus, the analysis has identified mutations that are relevant to modification of viral RNA synthesis, present in the mutant clouds of SARS-CoV-2 quasispecies. These mutations most likely occurred during intrahost diversification in several COVID-19 patients, during an initial stage of the pandemic, and within a brief time period.

SARS-CoV-2 rapidly evolves lineage-specific phenotypic differences when passaged repeatedly in immune-naïve mice

The persistence of SARS-CoV-2 despite the development of vaccines and a degree of herd immunity is partly due to viral evolution reducing vaccine and treatment efficacy. Serial infections of wild-type (WT) SARS-CoV-2 in Balb/c mice yield mouse-adapted strains with greater infectivity and mortality. We investigate if passaging unmodified B.1.351 (Beta) and B.1.617.2 (Delta) 20 times in K18-ACE2 mice, expressing the human ACE2 receptor, in a BSL-3 laboratory without selective pressures, drives human health-relevant evolution and if evolution is lineage-dependent. Late-passage virus causes more severe disease, at organism and lung tissue scales, with late-passage Delta demonstrating antibody resistance and interferon suppression. This resistance co-occurs with a de novo spike S371F mutation, linked with both traits. S371F, an Omicron-characteristic mutation, is co-inherited at times with spike E1182G per Nanopore sequencing, existing in different within-sample viral variants at others. Both S371F and E1182G are linked to mammalian GOLGA7 and ZDHHC5 interactions, which mediate viral-cell entry and antiviral response. This study demonstrates SARS-CoV-2's tendency to evolve with phenotypic consequences, its evolution varying by lineage, and suggests non-dominant quasi-species contribution.

Tracking the evolution of the SARS-CoV-2 Delta variant of concern: analysis of genetic diversity and selection across the whole viral genome

Background

Since 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has diversified extensively, producing five highly virulent lineages designated as variants of concern (VOCs). The Delta VOC emerged in India with increased transmission, immune evasion, and mortality, causing a massive global case surge in 2021. This study aims to understand how the Delta VOC evolved by characterizing mutation patterns in the viral population before and after its emergence. Furthermore, we aim to identify the influence of positive and negative selection on VOC evolution and understand the prevalence of different mutation types in the viral genome.

Methods

Three groups of whole viral genomes were retrieved from GISAID, sourced from India, with collection periods as follows: Group A-during the initial appearance of SARS-CoV-2; Group B-just before the emergence of the Delta variant; Group C-after the establishment of the Delta variant in India. Mutations in >1% of each group were identified with BioEdit to reveal differences in mutation quantity and type. Sites under positive or negative selection were identified with FUBAR. The results were compared to determine how mutations correspond with selective pressures and how viral mutation profiles changed to reflect genetic diversity before and after VOC emergence.

Results

The number of mutations increased progressively in Groups A-C, with Group C reporting a 2.2- and 1.9-fold increase from Groups A and B, respectively. Among all the observed mutations, Group C had the highest percentage of deletions (22.7%; vs. 4.2% and 2.6% in Groups A and B, respectively), and most mutations altered the final amino acid code, such as non-synonymous substitutions and deletions. Conversely, Group B had the most synonymous substitutions that are effectively silent. The number of sites experiencing positive selection increased in Groups A-C, but Group B had 2.4- and 2.6 times more sites under negative selection compared to Groups A and C, respectively.

Conclusion

Our findings demonstrated that viral genetic diversity continuously increased during and after the emergence of the Delta VOC. Despite this, Group B reports heightened negative selection, which potentially preserves important gene regions during evolution. Group C contains an unprecedented quantity of mutations and positively selected sites, providing strong evidence of active viral adaptation in the population.

Assessing the evolution of SARS-CoV-2 lineages and the dynamic associations between nucleotide variations

Despite seminal advances towards understanding the infection mechanism of SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2), it continues to cause significant morbidity and mortality worldwide. Though mass immunization programmes have been implemented in several countries, the viral transmission cycle has shown a continuous progression in the form of multiple waves. A constant change in the frequencies of dominant viral lineages, arising from the accumulation of nucleotide variations (NVs) through favourable selection, is understandably expected to be a major determinant of disease severity and possible vaccine escape. Indeed, worldwide efforts have been initiated to identify specific virus lineage(s) and/or NVs that may cause a severe clinical presentation or facilitate vaccination breakthrough. Since host genetics is expected to play a major role in shaping virus evolution, it is imperative to study the role of genome-wide SARS-CoV-2 NVs across various populations. In the current study, we analysed the whole genome sequence of 3543 SARS-CoV-2-infected samples obtained from the state of Telangana, India (including 210 from our previous study), collected over an extended period from April 2020 to October 2021. We present a unique perspective on the evolution of prevalent virus lineages and NVs during this period. We also highlight the presence of specific NVs likely to be associated favourably with samples classified as vaccination breakthroughs. Finally, we report genome-wide intra-host variations at novel genomic positions. The results presented here provide critical insights into virus evolution over an extended period and pave the way to rigorously investigate the role of specific NVs in vaccination breakthroughs.

ACE2 receptor polymorphism in humans and animals increases the risk of the emergence of SARS-CoV-2 variants during repeated intra- and inter-species host-switching of the virus

Like other coronaviruses, SARS-CoV-2 has ability to spread through human-to-human transmission and to circulate from humans to animals and from animals to humans. A high frequency of SARS-CoV-2 mutations has been observed in the viruses isolated from both humans and animals, suggesting a genetic fitness under positive selection in both ecological niches. The most documented positive selection force driving SARS-CoV-2 mutations is the host-specific immune response. However, after electrostatic interactions with lipid rafts, the first contact between the virus and host proteins is the viral spike-cellular receptor binding. Therefore, it is likely that the first level of selection pressure impacting viral fitness relates to the virus's affinity for its receptor, the angiotensin I converting enzyme 2 (ACE2). Although sufficiently conserved in a huge number of species to support binding of the viral spike with enough affinity to initiate fusion, ACE2 is highly polymorphic both among species and within a species. Here, we provide evidence suggesting that when the viral spike-ACE2 receptor interaction is not optimal, due to host-switching, mutations can be selected to improve the affinity of the spike for the ACE2 expressed by the new host. Notably, SARS-CoV-2 is mutation-prone in the spike receptor binding domain (RBD), allowing a better fit for ACE2 orthologs in animals. It is possibly that this may also be true for rare human alleles of ACE2 when the virus is spreading to billions of people. In this study, we present evidence that human subjects expressing the rare E₃₂₉G allele of ACE2 with higher allele frequencies in European populations exhibit a improved affinity for the SARS-CoV-2 spike N₅₀₁Y variant of the virus. This may suggest that this viral N₅₀₁Y variant emerged in the human population after SARS-CoV-2 had infected a human carrying the rare E₃₂₉G allele of ACE2. In addition, this viral evolution could impact viral replication as well as the ability of the adaptive humoral response to control infection with RBD-specific neutralizing antibodies. In a shifting landscape, this ACE2-driven genetic drift of SARS-CoV-2 which we have named the 'boomerang effect', could complicate the challenge of preventing COVID with a SARS-CoV-2 spike-derived vaccine.

CoVigator-A Knowledge Base for Navigating SARS-CoV-2 Genomic Variants

BACKGROUND

The outbreak of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) resulted in the global COVID-19 pandemic. The urgency for an effective SARS-CoV-2 vaccine has led to the development of the first series of vaccines at unprecedented speed. The discovery of SARS-CoV-2 spike-glycoprotein mutants, however, and consequentially the potential to escape vaccine-induced protection and increased infectivity, demonstrates the persisting importance of monitoring SARS-CoV-2 mutations to enable early detection and tracking of genomic variants of concern.

RESULTS

We developed the CoVigator tool with three components: (1) a knowledge base that collects new SARS-CoV-2 genomic data, processes it and stores its results; (2) a comprehensive variant calling pipeline; (3) an interactive dashboard highlighting the most relevant findings. The knowledge base routinely downloads and processes virus genome assemblies or raw sequencing data from the COVID-19 Data Portal (C19DP) and the European Nucleotide Archive (ENA), respectively. The results of variant calling are visualized through the dashboard in the form of tables and customizable graphs, making it a versatile tool for tracking SARS-CoV-2 variants. We put a special emphasis on the identification of intrahost mutations and make available to the community what is, to the best of our knowledge, the largest dataset on SARS-CoV-2 intrahost mutations. In the spirit of open data, all CoVigator results are available for download. The CoVigator dashboard is accessible via covigator.tron-mainz.de.

CONCLUSIONS

With increasing demand worldwide in genome surveillance for tracking the spread of SARS-CoV-2, CoVigator will be a valuable resource of an up-to-date list of mutations, which can be incorporated into global efforts.

Sequential intrahost evolution and onward transmission of SARS-CoV-2 variants

Persistent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections have been reported in immune-compromised individuals and people undergoing immune-modulatory treatments. Although intrahost evolution has been documented, direct evidence of subsequent transmission and continued stepwise adaptation is lacking. Here we describe sequential persistent SARS-CoV-2 infections in three individuals that led to the emergence, forward transmission, and continued evolution of a new Omicron sublineage, BA.1.23, over an eight-month period. The initially transmitted BA.1.23 variant encoded seven additional amino acid substitutions within the spike protein (E96D, R346T, L455W, K458M, A484V, H681R, A688V), and displayed substantial resistance to neutralization by sera from boosted and/or Omicron BA.1-infected study participants. Subsequent continued BA.1.23 replication resulted in additional substitutions in the spike protein (S254F, N448S, F456L, M458K, F981L, S982L) as well as in five other virus proteins. Our findings demonstrate not only that the Omicron BA.1 lineage can diverge further from its already exceptionally mutated genome but also that patients with persistent infections can transmit these viral variants. Thus, there is, an urgent need to implement strategies to prevent prolonged SARS-CoV-2 replication and to limit the spread of newly emerging, neutralization-resistant variants in vulnerable patients.

A comprehensive insight on the challenges for COVID-19 vaccine: A lesson learnt from other viral vaccines

The aim of this study is to comprehensively analyze previous viral vaccine programs and identify potential challenges and effective measures for the COVID-19 vaccine program. Previous viral vaccine programs, such as those for HIV, Zika, Influenza, Ebola, Dengue, SARS, and MERS, were evaluated. Paramount challenges were identified, including quasi-species, cross-reactivity, duration of immunity, revaccination, mutation, immunosenescence, and adverse events related to viral vaccines. Although a large population has been vaccinated, mutations in SARS-CoV-2 and adverse events related to vaccines pose significant challenges. Previous vaccine programs have taught us that predicting the final outcome of the current vaccine program for COVID-19 cannot be determined at a given state. Long-term follow-up studies are essential. Validated preclinical studies, long-term follow-up studies, alternative therapeutic approaches, and alternative vaccines are necessary.

Unravelling Antigenic Cross-Reactions toward the World of Coronaviruses: Extent of the Stability of Shared Epitopes and SARS-CoV-2 Anti-Spike Cross-Neutralizing Antibodies

The human immune repertoire retains the molecular memory of a very great diversity of target antigens (epitopes) and can recall this upon a second encounter with epitopes against which it has previously been primed. Although genetically diverse, proteins of coronaviruses exhibit sufficient conservation to lead to antigenic cross-reactions. In this review, our goal is to question whether pre-existing immunity against seasonal human coronaviruses (HCoVs) or exposure to animal CoVs has influenced the susceptibility of human populations to SARS-CoV-2 and/or had an impact upon the physiopathological outcome of COVID-19. With the hindsight that we now have regarding COVID-19, we conclude that although antigenic cross-reactions between different coronaviruses exist, cross-reactive antibody levels (titers) do not necessarily reflect on memory B cell frequencies and are not always directed against epitopes which confer cross-protection against SARS-CoV-2. Moreover, the immunological memory of these infections is short-term and occurs in only a small percentage of the population. Thus, in contrast to what might be observed in terms of cross-protection at the level of a single individual recently exposed to circulating coronaviruses, a pre-existing immunity against HCoVs or other CoVs can only have a very minor impact on SARS-CoV-2 circulation at the level of human populations.

Assessing the hidden diversity underlying consensus sequences of SARS-CoV-2 using VICOS, a novel bioinformatic pipeline for identification of mixed viral populations

INTRODUCTION

Coinfection with two SARS-CoV-2 viruses is still a very understudied phenomenon. Although next generation sequencing methods are very sensitive to detect heterogeneous viral populations in a sample, there is no standardized method for their characterization, so their clinical and epidemiological importance is unknown.

MATERIAL AND METHODS

We developed VICOS (Viral COinfection Surveillance), a new bioinformatic algorithm for variant calling, filtering and statistical analysis to identify samples suspected of being mixed SARS-CoV-2 populations from a large dataset in the framework of a community genomic surveillance. VICOS was used to detect SARS-CoV-2 coinfections in a dataset of 1,097 complete genomes collected between March 2020 and August 2021 in Argentina.

RESULTS

We detected 23 cases (2%) of SARS-CoV-2 coinfections. Detailed study of VICOS's results together with additional phylogenetic analysis revealed 3 cases of coinfections by two viruses of the same lineage, 2 cases by viruses of different genetic lineages, 13 were compatible with both coinfection and intra-host evolution, and 5 cases were likely a product of laboratory contamination.

DISCUSSION

Intra-sample viral diversity provides important information to understand the transmission dynamics of SARS-CoV-2. Advanced bioinformatics tools, such as VICOS, are a necessary resource to help unveil the hidden diversity of SARS-CoV-2.

Quasispecies Nature of RNA Viruses: Lessons from the Past

Viral quasispecies are distinct but closely related mutants formed by the disparity in viral genomes due to recombination, mutations, competition, and selection pressure. Theoretical derivation for the origin of a quasispecies is owed to the error-prone replication by polymerase and mutants of RNA replicators. Here, we briefly addressed the theoretical and mathematical origin of quasispecies and their dynamics. The impact of quasispecies for major salient human pathogens is reviewed. In the current global scenario, rapid changes in geographical landscapes favor the origin and selection of mutants. It comes as no surprise that a cauldron of mutants poses a significant risk to public health, capable of causing pandemics. Mutation rates in RNA viruses are magnitudes higher than in DNA organisms, explaining their enhanced virulence and evolvability. RNA viruses cause the most devastating pandemics; for example, members of the Orthomyxoviridae family caused the great influenza pandemic (1918 flu or Spanish flu), the SARS (severe acute respiratory syndrome) and MERS (Middle East respiratory syndrome) outbreak, and the human immunodeficiency viruses (HIV), lentiviruses of the Retroviridae family, caused worldwide devastation. Rapidly evolving RNA virus populations are a daunting challenge for the designing of effective control measures like vaccines. Developing awareness of the evolutionary dispositions of RNA viral mutant spectra and what influences their adaptation and virulence will help curtail outbreaks of past and future pathogens.

Viral Fitness, Population Complexity, Host Interactions, and Resistance to Antiviral Agents

Fitness of viruses has become a standard parameter to quantify their adaptation to a biological environment. Fitness determinations for RNA viruses (and some highly variable DNA viruses) meet with several uncertainties. Of particular interest are those that arise from mutant spectrum complexity, absence of population equilibrium, and internal interactions among components of a mutant spectrum. Here, concepts, fitness measurements, limitations, and current views on experimental viral fitness landscapes are discussed. The effect of viral fitness on resistance to antiviral agents is covered in some detail since it constitutes a widespread problem in antiviral pharmacology, and a challenge for the design of effective antiviral treatments. Recent evidence with hepatitis C virus suggests the operation of mechanisms of antiviral resistance additional to the standard selection of drug-escape mutants. The possibility that high replicative fitness may be the driver of such alternative mechanisms is considered. New broad-spectrum antiviral designs that target viral fitness may curtail the impact of drug-escape mutants in treatment failures. We consider to what extent fitness-related concepts apply to coronaviruses and how they may affect strategies for COVID-19 prevention and treatment.

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.

Role of spike compensatory mutations in the interspecies transmission of SARS-CoV-2

SARS-CoV-2, the virus responsible for COVID-19 in humans, can efficiently infect a large number of animal species. Like any virus, and particularly RNA viruses, SARS-CoV-2 undergoes mutations during its life cycle some of which bring a selective advantage, leading to the selection of a given lineage. Minks are very susceptible to SARS-CoV-2 and owing to their presence in mass rearing, they make a good model for studying the relative importance of mutations in viral adaptation to host species. Variants, such as the mink-selected SARS-CoV-2 Y453F and D614G or H69del/V70del, Y453F, I692V and M1229I were identified in humans after spreading through densely caged minks. However, not all mink-specific mutations are conserved when the virus infects human populations back. Many questions remain regarding the interspecies evolution of SARS-CoV-2 and the dynamics of transmission leading to the emergence of new variant strains. We compared the human and mink ACE2 receptor structures and their interactions with SARS-CVoV-2 variants. In minks, ACE2 presents a Y34 amino acid instead of the H34 amino acid found in the human ACE2. H34 is essential for the interaction with the Y453 residue of the SARS-CoV-2 Spike protein. The Y453F mink mutation abolishes this conflict. A series of 18 mutations not involved in the direct ACE2 interaction was observed in addition to the Y453F and D614G in 16 different SARS-CoV-2 strains following bidirectional infections between humans and minks. These mutations were not random and were distributed into five different functional groups having an effect on the kinetics of ACE2-RD interaction. The interspecies transmission of SARS-CoV-2 from humans to minks and back to humans, generated specific mutations in each species which improved the affinity for the ACE2 receptor either by direct mutation of the core 453 residue or by associated compensatory mutations.

As the SARS-CoV-2 virus evolves, should Omicron subvariant BA.2 be subjected to quarantine, or should we learn to live with it?

It has been nearly 35 months since the COVID-19 outbreak. The pathogen SARS-CoV-2 has evolved into several variants. Among them, Omicron is the fifth variant of concern which have rapidly spread globally during the past 8 months. Omicron variant shows different characteristics from previous variants, which is highly infectious, highly transmissible, minimally pathogenic, vaccine and antibody tolerant; however, it is less likely to cause severe illness, resulting in fewer deaths. Omicron has evolved into five main lineages, including BA.1, BA.2, BA.3, BA.4, and BA.5. Before BA.5, Omicron BA.2 sublineage was the dominant strain all over the world for several months. The experience of prevention and treatment against BA.2 is worth studying and learning for overcoming other Omicron subvariants. Although the Omicron subvariant BA.2 is significantly less severe than that caused by ancestral strains, it is still far more dangerous than influenza, and its long-term sequelae are unknown. Effective treatments are currently limited; therefore, effective defense may be the key to controlling the epidemic today, rather than just "living with" the virus.

SARS-CoV-2 quasi-species analysis from patients with persistent nasopharyngeal shedding

At the time of a new and unprecedented viral pandemic, many questions are being asked about the genomic evolution of SARS-CoV-2 and the emergence of different variants, leading to therapeutic and immune evasion and survival of this genetically highly labile RNA virus. The nasopharyngeal persistence of infectious virus beyond 17 days proves its constant interaction with the human immune system and increases the intra-individual mutational possibilities. We performed a prospective high-throughput sequencing study (ARTIC Nanopore) of SARS-CoV-2 from so-called "persistent" patients, comparing them with a non-persistent population, and analyzing the quasi-species present in a single sample at time t. Global intra-individual variability in persistent patients was found to be higher than in controls (mean 5.3%, Standard deviation 0.9 versus 4.6% SD 0.3, respectively, p < 0.001). In the detailed analysis, we found a greater difference between persistent and non-persistent patients with non-severe COVID 19, and between the two groups infected with clade 20A. Furthermore, we found minority N501Y and P681H mutation clouds in all patients, with no significant differences found both groups. The question of the SARS-CoV-2 viral variants' genesis remains to be further investigated, with the need to prevent new viral propagations and their consequences, and quasi-species analysis could be an important key to watch out.

SARS-CoV-2 Variants Identification: Overview of Molecular Existing Methods

Since the beginning of COVID-19 pandemic the Real Time sharing of genome sequences of circulating virus supported the diagnostics and surveillance of SARS-CoV-2 and its transmission dynamics. SARS-CoV-2 straightaway showed its tendency to mutate and adapt to the host, culminating in the emergence of variants; so it immediately became of crucial importance to be able to detect them quickly but also to be able to monitor in depth the changes on the whole genome to early identify the new possibly emerging variants. In this scenario, this manuscript aims to provide an overview of the existing methods for the identification of SARS-CoV-2 variants (from rapid method based on identification of one or more specific mutations to Whole Genome sequencing approach-WGS), taking into account limitations, advantages and applications of them in the field of diagnosis and surveillance of SARS-CoV-2.

Comparative analysis of SARS-CoV-2 quasispecies in the upper and lower respiratory tract shows an ongoing evolution in the spike cleavage site

Studies are needed to better understand the genomic evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This study aimed to describe viral quasispecies population of upper and lower respiratory tract by next-generation sequencing in patients admitted to intensive care unit. A deep sequencing of the S gene of SARS-CoV-2 from 109 clinical specimens, sampled from the upper respiratory tract (URT) and lower respiratory tract (LRT) of 77 patients was performed. A higher incidence of non-synonymous mutations and indels was observed in the LRT among minority variants. This might be explained by the ability of the virus to invade cells without interacting with ACE2 (e.g. exploiting macrophage phagocytosis). Minority variants are highly concentrated around the gene portion encoding for the Spike cleavage site, with a higher incidence in the URT; four mutations are highly recurring among samples and were found associated with the URT. Interestingly, 55.8% of minority variants detected in this locus were T>G and G>T transversions. Results from this study evidenced the presence of selective pressure and suggest that an evolutionary process is still ongoing in one of the crucial sites of spike protein associated with the spillover to humans.

Targeted Virome Sequencing Enhances Unbiased Detection and Genome Assembly of Known and Emerging Viruses-The Example of SARS-CoV-2

Targeted virome enrichment and sequencing (VirCapSeq-VERT) utilizes a pool of oligos (baits) to enrich all known—up to 2015—vertebrate-infecting viruses, increasing their detection sensitivity. The hybridisation of the baits to the target sequences can be partial, thus enabling the detection and genomic reconstruction of novel pathogens with <40% genetic diversity compared to the strains used for the baits’ design. In this study, we deploy this method in multiplexed mixes of viral extracts, and we assess its performance in the unbiased detection of DNA and RNA viruses after cDNA synthesis. We further assess its efficiency in depleting various background genomic material. Finally, as a proof-of-concept, we explore the potential usage of the method for the characterization of unknown, emerging human viruses, such as SARS-CoV-2, which may not be included in the baits’ panel. We mixed positive samples of equimolar DNA/RNA viral extracts from SARS-CoV-2, coronavirus OC43, cytomegalovirus, influenza A virus H3N2, parvovirus B19, respiratory syncytial virus, adenovirus C and coxsackievirus A16. Targeted virome enrichment was performed on a dsDNA mix, followed by sequencing on the NextSeq500 (Illumina) and the portable MinION sequencer, to evaluate its usability as a point-of-care (PoC) application. Genome mapping assembly was performed using viral reference sequences. The untargeted libraries contained less than 1% of total reads mapped on most viral genomes, while RNA viruses remained undetected. In the targeted libraries, the percentage of viral-mapped reads were substantially increased, allowing full genome assembly in most cases. Targeted virome sequencing can enrich a broad range of viruses, potentially enabling the discovery of emerging viruses.

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.

SARS-CoV-2 Within-Host and in vitro Genomic Variability and Sub-Genomic RNA Levels Indicate Differences in Viral Expression Between Clinical Cohorts and in vitro Culture

Background

Low frequency intrahost single nucleotide variants (iSNVs) of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) have been increasingly recognised as predictive indicators of positive selection. Particularly as growing numbers of SARS-CoV-2 variants of interest (VOI) and concern (VOC) emerge. However, the dynamics of subgenomic RNA (sgRNA) expression and its impact on genomic diversity and infection outcome remain poorly understood. This study aims to investigate and quantify iSNVs and sgRNA expression in single and longitudinally sampled cohorts over the course of mild and severe SARS-CoV-2 infection, benchmarked against an in vitro infection model.

Methods

Two clinical cohorts of SARS-CoV-2 positive cases in New South Wales, Australia collected between March 2020 and August 2021 were sequenced. Longitudinal samples from cases hospitalised due to SARS-CoV-2 infection (severe) (n = 16) were analysed and compared with cases that presented with SARS-CoV-2 symptoms but were not hospitalised (mild) (n = 23). SARS-CoV-2 genomic diversity profiles were also examined from daily sampling of culture experiments for three SARS-CoV-2 variants (Lineage A, B.1.351, and B.1.617.2) cultured in VeroE6 C1008 cells (n = 33).

Results

Intrahost single nucleotide variants were detected in 83% (19/23) of the mild cohort cases and 100% (16/16) of the severe cohort cases. SNP profiles remained relatively fixed over time, with an average of 1.66 SNPs gained or lost, and an average of 4.2 and 5.9 low frequency variants per patient were detected in severe and mild infection, respectively. sgRNA was detected in 100% (25/25) of the mild genomes and 92% (24/26) of the severe genomes. Total sgRNA expressed across all genes in the mild cohort was significantly higher than that of the severe cohort. Significantly higher expression levels were detected in the spike and the nucleocapsid genes. There was significantly less sgRNA detected in the culture dilutions than the clinical cohorts.

Discussion and Conclusion

The positions and frequencies of iSNVs in the severe and mild infection cohorts were dynamic overtime, highlighting the importance of continual monitoring, particularly during community outbreaks where multiple SARS-CoV-2 variants may co-circulate. sgRNA levels can vary across patients and the overall level of sgRNA reads compared to genomic RNA can be less than 1%. The relative contribution of sgRNA to the severity of illness warrants further investigation given the level of variation between genomes. Further monitoring of sgRNAs will improve the understanding of SARS-CoV-2 evolution and the effectiveness of therapeutic and public health containment measures during the pandemic.

The third international hackathon for applying insights into large-scale genomic composition to use cases in a wide range of organisms

In October 2021, 59 scientists from 14 countries and 13 U.S. states collaborated virtually in the Third Annual Baylor College of Medicine & DNANexus Structural Variation hackathon. The goal of the hackathon was to advance research on structural variants (SVs) by prototyping and iterating on open-source software. This led to nine hackathon projects focused on diverse genomics research interests, including various SV discovery and genotyping methods, SV sequence reconstruction, and clinically relevant structural variation, including SARS-CoV-2 variants. Repositories for the projects that participated in the hackathon are available at https://github.com/collaborativebioinformatics.

Vaccine breakthrough infections with SARS-CoV-2 Alpha mirror mutations in Delta Plus, Iota, and Omicron

Replication of SARS-CoV-2 in the human population is defined by distributions of mutants that are present at different frequencies within the infected host and can be detected by ultra-deep sequencing techniques. In this study, we examined the SARS-CoV-2 mutant spectra of amplicons from the spike-coding (S-coding) region of 5 nasopharyngeal isolates derived from patients with vaccine breakthrough. Interestingly, all patients became infected with the Alpha variant, but amino acid substitutions that correspond to the Delta Plus, Iota, and Omicron variants were present in the mutant spectra of the resident virus. Deep sequencing analysis of SARS-CoV-2 from patients with vaccine breakthrough revealed a rich reservoir of mutant types and may also identify tolerated substitutions that can be represented in epidemiologically dominant variants.

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.

Molecular and Clinical Prognostic Biomarkers of COVID-19 Severity and Persistence

The coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), poses several challenges to clinicians, due to its unpredictable clinical course. The identification of laboratory biomarkers, specific cellular, and molecular mediators of immune response could contribute to the prognosis and management of COVID-19 patients. Of utmost importance is also the detection of differentially expressed genes, which can serve as transcriptomic signatures, providing information valuable to stratify patients into groups, based on the severity of the disease. The role of biomarkers such as IL-6, procalcitonin, neutrophil–lymphocyte ratio, white blood cell counts, etc. has already been highlighted in recently published studies; however, there is a notable amount of new evidence that has not been summarized yet, especially regarding transcriptomic signatures. Hence, in this review, we assess the latest cellular and molecular data and determine the significance of abnormalities in potential biomarkers for COVID-19 severity and persistence. Furthermore, we applied Gene Ontology (GO) enrichment analysis using the genes reported as differentially expressed in the literature in order to investigate which biological pathways are significantly enriched. The analysis revealed a number of processes, such as inflammatory response, and monocyte and neutrophil chemotaxis, which occur as part of the complex immune response to SARS-CoV-2.

SARS-CoV-2 Mutations and Their Impact on Diagnostics, Therapeutics and Vaccines

With the high rate of COVID-19 infections worldwide, the emergence of SARS-CoV-2 variants was inevitable. Several mutations have been identified in the SARS-CoV-2 genome, with the spike protein as one of the mutational hot spots. Specific amino acid substitutions such as D614G and N501Y were found to alter the transmissibility and virulence of the virus. The WHO has classified the variants identified with fitness-enhancing mutations as variants of concern (VOC), variants of interest (VOI) or variants under monitoring (VUM). The VOCs pose an imminent threat as they exhibit higher transmissibility, disease severity and ability to evade vaccine-induced and natural immunity. Here we review the mutational landscape on the SARS-CoV-2 structural and non-structural proteins and their impact on diagnostics, therapeutics and vaccines. We also look at the effectiveness of approved vaccines, antibody therapy and convalescent plasma on the currently prevalent VOCs, which are B.1.17, B.1.351, P.1, B.1.617.2 and B.1.1.529. We further discuss the possible factors influencing mutation rates and future directions.

Efficacy of mRNA, adenoviral vector, and perfusion protein COVID-19 vaccines

Coronavirus disease 2019 (COVID-19) has a devastating impact on global populations triggered by a highly infectious viral sickness, produced by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The third major cause of mortality in the United States, following heart disease and cancer in 2020, was undoubtedly COVID-19. The centers for disease control and prevention (CDC) and the world health organization (WHO) separately developed a categorization system for differentiating new strains of SARS-CoV-2 into variants of concern (VoCs) and variants of interest (VoIs) with the continuing development of various strains SARS-CoV-2. By December 2021, five of the SARS-CoV-2 VoCs were discovered from the onset of the pandemic depending on the latest epidemiologic report by the WHO: Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2), and Omicron (B.1.1.529). Mutations in the receptor-binding domain (RBD) and n-terminal domain (NTD) have been found throughout all five identified VoCs. All strains other than the delta mutant are often found with the N501Y mutation situated on the RBD, resulting in higher binding between the spike protein and angiotensin-converting enzyme 2 (ACE2) receptors, enhanced viral adhesion, and following the entrance to host cells. The introduction of these new strains of SRAS-CoV-2 is likely to overcome the remarkable achievements gained in restricting this viral disease to the point where it is presented with remarkable vaccine developments against COVID-19 and strong worldwide mass immunization initiatives. Throughout this literature review, the effectiveness of current COVID-19 vaccines for managing and prohibiting SARS-CoV-2 strains is thoroughly described.

Space and Genotype-Dependent Virus Distribution during Infection Progression

The paper is devoted to a nonlocal reaction-diffusion equation describing the development of viral infection in tissue, taking into account virus distribution in the space of genotypes, the antiviral immune response, and natural genotype-dependent virus death. It is shown that infection propagates as a reaction-diffusion wave. In some particular cases, the 2D problem can be reduced to a 1D problem by separation of variables, allowing for proof of wave existence and stability. In general, this reduction provides an approximation of the 2D problem by a 1D problem. The analysis of the reduced problem allows us to determine how viral load and virulence depend on genotype distribution, the strength of the immune response, and the level of immunity.

Jumping a Moving Train: SARS-CoV-2 Evolution in Real Time

The field of molecular epidemiology responded to the SARS-CoV-2 pandemic with an unrivaled amount of whole viral genome sequencing. By the time this sentence is published we will have well surpassed 1.5 million whole genomes, more than 4 times the number of all microbial whole genomes deposited in GenBank and 35 times the total number of viral genomes. This extraordinary dataset that accrued in near real time has also given us an opportunity to chart the global and local evolution of a virus as it moves through the world population. The data itself presents challenges that have never been dealt with in molecular epidemiology, and tracking a virus that is changing so rapidly means that we are often running to catch up. Here we review what is known about the evolution of the virus, and the critical impact that whole genomes have had on our ability to trace back and track forward the spread of lineages of SARS-CoV-2. We then review what whole genomes have told us about basic biological properties of the virus such as transmissibility, virulence, and immune escape with a special emphasis on pediatric disease. We couch this discussion within the framework of systematic biology and phylogenetics, disciplines that have proven their worth again and again for identifying and deciphering the spread of epidemics, though they were largely developed in areas far removed from infectious disease and medicine.

A next generation sequencing (NGS) analysis to reveal genomic and proteomic mutation landscapes of SARS-CoV-2 in South Asia

•

48 SNPs were identified from the genome-wide analysis of 410 South Asian SARS-CoV-2 sequences.

•

About 85% SNPs are packaged in ORF1ab, spike protein, and nucleocapsid.

•

South Asian strains are highly related to the South American and European strains according to the phylogenetic analysis.

•

Unlike other countries, frequency of 1163A>T missense mutation is very high (78.80%) in Bangladeshi samples.

Counts for SARS-CoV-2 associated infections and fatalities are on the rise globally even in regions which contained the spread momentarily. The pattern of infections has been found to be controlled by the distinctive selection pressures exerted by fluctuating environmental nature and hosts. A total of 410 whole-genome sequences submitted by the South Asian countries were retrieved from the GISAID database and analyzed to assess the impact and pattern of mutations in this region. Most common and frequent mutations in the South Asian countries are 241C > T, 3037C > T, 14408C > T, and 23403A > G and about 85% SNPs are localized in ORF1ab, spike protein, and nucleocapsid. Among the identified mutations, the proportion of missense type (54.17%) was highest, followed by the synonymous (41.66%) and the non-coding types (4.17%). While analyzing transmission source in terms of geolocation, the largest clustered group from the South Asian countries was based on the G-clade (D614G) (81.7%; 335/410 samples), tracing the inception and transmission of SARS-CoV-2 infections in the South Asian countries from European regions. Phylogenetic analysis also revealed that the South Asian strains are highly related to the South American and European strains. We found that G-clade mutations are more prevalent (96.19%) in the samples of Bangladesh which were also prevalent in the European isolates. Surprisingly, one missense mutation (1163A > T) in ORF1ab gene became dominant only in Bangladesh (78.8%), which led to debates regarding effects on the pathogenicity and transmissibility of the virus. Overall, the findings of this study highlight the frequently mutated SARS-CoV-2 variants among the COVID-19 patients in the South Asian countries which might ease containment of the disease in this region through investigating the virulence reducing factors as the identified mutations are strongly correlated with low infection and mortality rate.

Jumper enables discontinuous transcript assembly in coronaviruses

Genes in SARS-CoV-2 and other viruses in the order of Nidovirales are expressed by a process of discontinuous transcription which is distinct from alternative splicing in eukaryotes and is mediated by the viral RNA-dependent RNA polymerase. Here, we introduce the DISCONTINUOUS TRANSCRIPT ASSEMBLYproblem of finding transcripts and their abundances given an alignment of paired-end short reads under a maximum likelihood model that accounts for varying transcript lengths. We show, using simulations, that our method, JUMPER, outperforms existing methods for classical transcript assembly. On short-read data of SARS-CoV-1, SARS-CoV-2 and MERS-CoV samples, we find that JUMPER not only identifies canonical transcripts that are part of the reference transcriptome, but also predicts expression of non-canonical transcripts that are supported by subsequent orthogonal analyses. Moreover, application of JUMPER on samples with and without treatment reveals viral drug response at the transcript level. As such, JUMPER enables detailed analyses of Nidovirales transcriptomes under varying conditions.

A common methodological phylogenomics framework for intra-patient heteroplasmies to infer SARS-CoV-2 sublineages and tumor clones

BACKGROUND

All diseases containing genetic material undergo genetic evolution and give rise to heterogeneity including cancer and infection. Although these illnesses are biologically very different, the ability for phylogenetic retrodiction based on the genomic reads is common between them and thus tree-based principles and assumptions are shared. Just as the different frequencies of tumor genomic variants presupposes the existence of multiple tumor clones and provides a handle to computationally infer them, we postulate that the different variant frequencies in viral reads offers the means to infer multiple co-infecting sublineages.

RESULTS

We present a common methodological framework to infer the phylogenomics from genomic data, be it reads of SARS-CoV-2 of multiple COVID-19 patients or bulk DNAseq of the tumor of a cancer patient. We describe the Concerti computational framework for inferring phylogenies in each of the two scenarios.To demonstrate the accuracy of the method, we reproduce some known results in both scenarios. We also make some additional discoveries.

CONCLUSIONS

Concerti successfully extracts and integrates information from multi-point samples, enabling the discovery of clinically plausible phylogenetic trees that capture the heterogeneity known to exist both spatially and temporally. These models can have direct therapeutic implications by highlighting "birth" of clones that may harbor resistance mechanisms to treatment, "death" of subclones with drug targets, and acquisition of functionally pertinent mutations in clones that may have seemed clinically irrelevant. Specifically in this paper we uncover new potential parallel mutations in the evolution of the SARS-CoV-2 virus. In the context of cancer, we identify new clones harboring resistant mutations to therapy.

Organ-specific genome diversity of replication-competent SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is not always confined to the respiratory system, as it impacts people on a broad clinical spectrum from asymptomatic to severe systemic manifestations resulting in death. Further, accumulation of intra-host single nucleotide variants during prolonged SARS-CoV-2 infection may lead to emergence of variants of concern (VOCs). Still, information on virus infectivity and intra-host evolution across organs is sparse. We report a detailed virological analysis of thirteen postmortem coronavirus disease 2019 (COVID-19) cases that provides proof of viremia and presence of replication-competent SARS-CoV-2 in extrapulmonary organs of immunocompromised patients, including heart, kidney, liver, and spleen (NCT04366882). In parallel, we identify organ-specific SARS-CoV-2 genome diversity and mutations of concern N501Y, T1027I, and Y453F, while the patient had died long before reported emergence of VOCs. These mutations appear in multiple organs and replicate in Vero E6 cells, highlighting their infectivity. Finally, we show two stages of fatal disease evolution based on disease duration and viral loads in lungs and plasma. Our results provide insights about the pathogenesis and intra-host evolution of SARS-CoV-2 and show that COVID-19 treatment and hygiene measures need to be tailored to specific needs of immunocompromised patients, even when respiratory symptoms cease.

Here the authors provide a detailed virological analysis of thirteen postmortem COVID-19 cases, including presence of replication-competent SARS-CoV-2 in extrapulmonary organs and tissue-specific patterns of SARS-CoV-2 genome diversity of an immunocompromised patient.

Historical Perspective on the Discovery of the Quasispecies Concept

Viral quasispecies are dynamic distributions of nonidentical but closely related mutant and recombinant viral genomes subjected to a continuous process of genetic variation, competition, and selection that may act as a unit of selection. The quasispecies concept owes its theoretical origins to a model for the origin of life as a collection of mutant RNA replicators. Independently, experimental evidence for the quasispecies concept was obtained from sampling of bacteriophage clones, which revealed that the viral populations consisted of many mutant genomes whose frequency varied with time of replication. Similar findings were made in animal and plant RNA viruses. Quasispecies became a theoretical framework to understand viral population dynamics and adaptability. The evidence came at a time when mutations were considered rare events in genetics, a perception that was to change dramatically in subsequent decades. Indeed, viral quasispecies was the conceptual forefront of a remarkable degree of biological diversity, now evident for cell populations and organisms, not only for viruses. Quasispecies dynamics unveiled complexities in the behavior of viral populations,with consequences for disease mechanisms and control strategies. This review addresses the origin of the quasispecies concept, its major implications on both viral evolution and antiviral strategies, and current and future prospects.

Intra-host evolution during SARS-CoV-2 prolonged infection

Long-term infection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) represents a challenge to virus dispersion and the control of coronavirus disease 2019 (COVID-19) pandemic. The reason why some people have prolonged infection and how the virus persists for so long are still not fully understood. Recent studies suggested that the accumulation of intra-host single nucleotide variants (iSNVs) over the course of the infection might play an important role in persistence as well as emergence of mutations of concern. For this reason, we aimed to investigate the intra-host evolution of SARS-CoV-2 during prolonged infection. Thirty-three patients who remained reverse transcription polymerase chain reaction (RT-PCR) positive in the nasopharynx for on average 18 days from the symptoms onset were included in this study. Whole-genome sequences were obtained for each patient at two different time points. Phylogenetic, populational, and computational analyses of viral sequences were consistent with prolonged infection without evidence of coinfection in our cohort. We observed an elevated within-host genomic diversity at the second time point samples positively correlated with cycle threshold (Ct) values (lower viral load). Direct transmission was also confirmed in a small cluster of healthcare professionals that shared the same workplace by the presence of common iSNVs. A differential accumulation of missense variants between the time points was detected targeting crucial structural and non-structural proteins such as Spike and helicase. Interestingly, longitudinal acquisition of iSNVs in Spike protein coincided in many cases with SARS-CoV-2 reactive and predicted T cell epitopes. We observed a distinguishing pattern of mutations over the course of the infection mainly driven by increasing A→U and decreasing G→A signatures. G→A mutations may be associated with RNA-editing enzyme activities; therefore, the mutational profiles observed in our analysis were suggestive of innate immune mechanisms of the host cell defense. Therefore, we unveiled a dynamic and complex landscape of host and pathogen interaction during prolonged infection of SARS-CoV-2, suggesting that the host’s innate immunity shapes the increase of intra-host diversity. Our findings may also shed light on possible mechanisms underlying the emergence and spread of new variants resistant to the host immune response as recently observed in COVID-19 pandemic.

Spread of Mink SARS-CoV-2 Variants in Humans: A Model of Sarbecovirus Interspecies Evolution

The rapid spread of SARS-CoV-2 variants has quickly spanned doubts and the fear about their ability escape vaccine protection. Some of these variants initially identified in caged were also found in humans. The claim that these variants exhibited lower susceptibility to antibody neutralization led to the slaughter of 17 million minks in Denmark. SARS-CoV-2 prevalence tests led to the discovery of infected farmed minks worldwide. In this study, we revisit the issue of the circulation of SARS-CoV-2 variants in minks as a model of sarbecovirus interspecies evolution by: (1) comparing human and mink angiotensin I converting enzyme 2 (ACE2) and neuropilin 1 (NRP-1) receptors; (2) comparing SARS-CoV-2 sequences from humans and minks; (3) analyzing the impact of mutations on the 3D structure of the spike protein; and (4) predicting linear epitope targets for immune response. Mink-selected SARS-CoV-2 variants carrying the Y453F/D614G mutations display an increased affinity for human ACE2 and can escape neutralization by one monoclonal antibody. However, they are unlikely to lose most of the major epitopes predicted to be targets for neutralizing antibodies. We discuss the consequences of these results for the rational use of SARS-CoV-2 vaccines.

Mortality in individuals treated with COVID-19 convalescent plasma varies with the geographic provenance of donors

Successful therapeutics and vaccines for coronavirus disease 2019 (COVID-19) have harnessed the immune response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Evidence that SARS-CoV-2 exists as locally evolving variants suggests that immunological differences may impact the effectiveness of antibody-based treatments such as convalescent plasma and vaccines. Considering that near-sourced convalescent plasma likely reflects the antigenic composition of local viral strains, we hypothesize that convalescent plasma has a higher efficacy, as defined by death within 30 days of transfusion, when the convalescent plasma donor and treated patient were in close geographic proximity. Results of a series of modeling techniques applied to approximately 28,000 patients from the Expanded Access to Convalescent Plasma program (ClinicalTrials.gov number: NCT04338360) support this hypothesis. This work has implications for the interpretation of clinical studies, the ability to develop effective COVID-19 treatments, and, potentially, for the effectiveness of COVID-19 vaccines as additional locally-evolving variants continue to emerge.

Regional differences in SARS-CoV-2 variants may affect treatment outcome. Here, the authors show that near-sourced convalescent plasma has higher efficacy, as defined by death within 30 days of transfusion, than plasma sourced more than 150 miles away.

SARS-CoV-2 Quasispecies Provides an Advantage Mutation Pool for the Epidemic Variants

The dynamics of quasispecies afford RNA viruses a great fitness on cell tropism and host range. To study the quasispecies features and the intra-host evolution of SARS-CoV-2, we collected nine confirmed patients and sequenced the haplotypes of spike gene using a single-molecule real-time platform. Fourteen samples were extracted from sputum, nasopharyngeal swabs, or stool, which in total produced 283,655 high-quality circular consensus sequences. We observed a stable quasispecies structure that one master mutant (mean abundance ∼0.70), followed by numerous minor mutants (mean abundance ∼1.21 × 10⁻³). Under high selective pressure, minor mutants may obtain a fitness advantage and become the master ones. The later predominant substitution D614G existed in the minor mutants of more than one early patient. An epidemic variant had a possibility to be independently originated from multiple hosts. The mutant spectrums covered ∼85% amino acid variations of public genomes (GISAID; frequency ≥ 0.1) and likely provided an advantage mutation pool for the current/future epidemic variants. Notably, 32 of 35 collected antibody escape substitutions were preexistent in the early quasispecies. Virus populations in different tissues/organs revealed potentially independent replications. The quasispecies complexity of sputum samples was significantly lower than that of nasopharyngeal swabs (P = 0.02). Evolution analysis revealed that three continuous S2 domains (HR1, CH, and CD) had undergone a positive selection. Cell fusion-related domains may play a crucial role in adapting to the intrahost immune system. Our findings suggested that future epidemiologic investigations and clinical interventions should consider the quasispecies information that has missed by routine single consensus genome.

IMPORTANCE RNA virus population in a host does not consist of a consensus single haplotype but rather an ensemble of related sequences termed quasispecies. The dynamics of quasispecies afford SARS-CoV-2 a great ability on genetic fitness during intrahost evolution. The process is likely achieved by changing the genetic characteristics of key functional genes, such as the spike glycoprotein. Previous studies have applied the next-generation sequencing (NGS) technology to evaluate the quasispecies of SARS-CoV-2, and results indicated a low genetic diversity of the spike gene. However, the NGS platform cannot directly obtain the full haplotypes without assembling, and it is also difficult to predict the extremely low-frequency variations. Therefore, we introduced a single-molecule real-time technology to directly obtain the haplotypes of the RNA population and further study the quasispecies features and intrahost evolution of the spike gene.

A Possible Role of Remdesivir and Plasma Therapy in the Selective Sweep and Emergence of New SARS-CoV-2 Variants

Since summer 2020, SARS-CoV-2 strains at the origin of the COVID-19 pandemic have suddenly been replaced by new SARS-CoV-2 variants, some of which are highly transmissible and spread at a high rate. These variants include the Marseille-4 lineage (Nextclade 20A.EU2) in Europe, the 20I/501Y.V1 variant first detected in the UK, the 20H/501Y.V2 variant first detected in South Africa, and the 20J/501Y.V3 variant first detected in Brazil. These variants are characterized by multiple mutations in the viral spike protein that is targeted by neutralizing antibodies elicited in response to infection or vaccine immunization. The usual coronavirus mutation rate through genetic drift alone cannot account for such rapid changes. Recent reports of the occurrence of such mutations in immunocompromised patients who received remdesivir and/or convalescent plasma or monoclonal antibodies to treat prolonged SARS-CoV-2 infections led us to hypothesize that experimental therapies that fail to cure the patients from COVID-19 could favor the emergence of immune escape SARS-CoV-2 variants. We review here the data that support this hypothesis and urge physicians and clinical trial promoters to systematically monitor viral mutations by whole-genome sequencing for patients who are administered these treatments.

Host- and Species-Dependent Quasispecies Divergence of Severe Acute Respiratory Syndrome Coronavirus-2 in Non-human Primate Models

Recently, newly emerging variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been continuously reported worldwide. However, the precise evaluation of SARS-CoV-2 microevolution in host is very limited because the exact genetic information of infected virus could not be acquired in human researches. In this report, we performed deep sequencing for seed virus and SARS-CoV-2 isolated in eight cynomolgus and rhesus macaques at 3 days postinoculation and evaluated single-nucleotide polymorphisms (SNPs) in SARS-CoV-2 by variant analysis. A total of 69 single-nucleotide variants (SNVs) were present in the 5′-untranslated region (UTR), 3′-UTR, ORF1ab, S, ORF3a, ORF8, and N genes of the seed virus passaged in VERO cells. Between those present on the seed virus and those on each SARS-CoV-2 isolated from the lungs of the macaques, a total of 29 variants was identified in 4 coding proteins (ORF1ab, S, ORF3a, and N) and non-coding regions (5′- and 3′-UTR). Variant number was significantly different according to individuals and ranged from 2 to 11. Moreover, the average major frequency variation was identified in six sites between the cynomolgus monkeys and rhesus macaques. As with diverse SNPs in SARS-CoV-2, the values of viral titers in lungs were significantly different according to individuals and species. Our study first revealed that the genomes of SARS-CoV-2 differ according to individuals and species despite infection of the identical virus in non-human primates (NHPs). These results are important for the interpretation of longitudinal studies evaluating the evolution of the SARS-CoV-2 in human beings and development of new diagnostics, vaccine, and therapeutics targeting SARS-CoV-2.

ACoRE: Accurate SARS-CoV-2 genome reconstruction for the characterization of intra-host and inter-host viral diversity in clinical samples and for the evaluation of re-infections

Sequencing the SARS-CoV-2 genome from clinical samples can be challenging, especially in specimens with low viral titer. Here we report Accurate SARS-CoV-2 genome Reconstruction (ACoRE), an amplicon-based viral genome sequencing workflow for the complete and accurate reconstruction of SARS-CoV-2 sequences from clinical samples, including suboptimal ones that would usually be excluded even if unique and irreplaceable. The protocol was optimized to improve flexibility and the combination of technical replicates was established as the central strategy to achieve accurate analysis of low-titer/suboptimal samples. We demonstrated the utility of the approach by achieving complete genome reconstruction and the identification of false-positive variants in >170 clinical samples, thus avoiding the generation of inaccurate and/or incomplete sequences. Most importantly, ACoRE was crucial to identify the correct viral strain responsible of a relapse case, that would be otherwise mis-classified as a re-infection due to missing or incorrect variant identification by a standard workflow.

SARS-CoV-2 Portrayed against HIV: Contrary Viral Strategies in Similar Disguise

SARS-CoV-2 and HIV are zoonotic viruses that rapidly reached pandemic scale, causing global losses and fear. The COVID-19 and AIDS pandemics ignited massive efforts worldwide to develop antiviral strategies and characterize viral architectures, biological and immunological properties, and clinical outcomes. Although both viruses have a comparable appearance as enveloped viruses with positive-stranded RNA and envelope spikes mediating cellular entry, the entry process, downstream biological and immunological pathways, clinical outcomes, and disease courses are strikingly different. This review provides a systemic comparison of both viruses' structural and functional characteristics, delineating their distinct strategies for efficient spread.

Tools and Techniques for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2)/COVID-19 Detection

The coronavirus disease 2019 (COVID-19) pandemic, caused by severe acute respiratory disease coronavirus 2 (SARS-CoV-2), has led to millions of confirmed cases and deaths worldwide. Efficient diagnostic tools are in high demand, as rapid and large-scale testing plays a pivotal role in patient management and decelerating disease spread. This paper reviews current technologies used to detect SARS-CoV-2 in clinical laboratories as well as advances made for molecular, antigen-based, and immunological point-of-care testing, including recent developments in sensor and biosensor devices. The importance of the timing and type of specimen collection is discussed, along with factors such as disease prevalence, setting, and methods. Details of the mechanisms of action of the various methodologies are presented, along with their application span and known performance characteristics. Diagnostic imaging techniques and biomarkers are also covered, with an emphasis on their use for assessing COVID-19 or monitoring disease severity or complications. While the SARS-CoV-2 literature is rapidly evolving, this review highlights topics of interest that have occurred during the pandemic and the lessons learned throughout. Exploring a broad armamentarium of techniques for detecting SARS-CoV-2 will ensure continued diagnostic support for clinicians, public health, and infection prevention and control for this pandemic and provide advice for future pandemic preparedness.

Characterization of SARS-CoV-2 different variants and related morbidity and mortality: a systematic review

INTRODUCTION

Coronavirus Disease-2019 (SARS-CoV-2) started its devastating trajectory into a global pandemic in Wuhan, China, in December 2019. Ever since, several variants of SARS-CoV-2 have been identified. In the present review, we aimed to characterize the different variants of SARS-CoV-2 and explore the related morbidity and mortality.

METHODS

A systematic review including the current evidence related to different variants of SARS-CoV-2 and the related morbidity and mortality was conducted through a systematic search utilizing the keywords in the online databases including Scopus, PubMed, Web of Science, and Science Direct; we retrieved all related papers and reports published in English from December 2019 to September 2020.

RESULTS

A review of identified articles has shown three main genomic variants, including type A, type B, and type C. we also identified three clades including S, V, and G. Studies have demonstrated that the C14408T and A23403G alterations in the Nsp12 and S proteins are the most prominent alterations in the world, leading to life-threatening mutations.The spike D614G amino acid change has become the most common variant since December 2019. From missense mutations found from Gujarat SARS-CoV-2 genomes, C28854T, deleterious mutation in the nucleocapsid (N) gene was significantly associated with patients' mortality. The other significant deleterious variant (G25563T) is found in patients located in Orf3a and has a potential role in viral pathogenesis.

CONCLUSION

Overall, researchers identified several SARS-CoV-2 variants changing clinical manifestations and increasing the transmissibility, morbidity, and mortality of COVID-19. This should be considered in current practice and interventions to combat the pandemic and prevent related morbidity and mortality.

Heterogeneity of SARS-CoV-2 virus produced in cell culture revealed by shotgun proteomics and supported by genome sequencing

COVID-19 is the most disturbing pandemic of the past hundred years. Its causative agent, the SARS-CoV-2 virus, has been the subject of an unprecedented investigation to characterize its molecular structure and intimate functioning. While markers for its detection have been proposed and several diagnostic methodologies developed, its propensity to evolve and evade diagnostic tools and the immune response is of great concern. The recent spread of new variants with increased infectivity requires even more attention. Here, we document how shotgun proteomics can be useful for rapidly monitoring the evolution of the SARS-CoV-2 virus. We evaluated the heterogeneity of purified SARS-CoV-2 virus obtained after culturing in the Vero E6 cell line. We found that cell culture induces significant changes that are translated at the protein level, such changes being detectable by tandem mass spectrometry. Production of viral particles requires careful quality control which can be easily performed by shotgun proteomics. Although considered relatively stable so far, the SARS-CoV-2 genome turns out to be prone to frequent variations. Therefore, the sequencing of SARS-CoV-2 variants from patients reporting only the consensus genome after its amplification would deserve more attention and could benefit from more in-depth analysis of low level but crystal-clear signals, as well as complementary and rapid analysis by shotgun proteomics.