The emergence, spread and vanishing of a French SARS-CoV-2 variant exemplifies the fate of RNA virus epidemics and obeys the Mistigri rule

Role of SARS-CoV-2 mutations in the evolution of the COVID-19 pandemic

OBJECTIVES

The SARS-CoV-2 pandemic and large-scale genomic surveillance provided an exceptional opportunity to analyze mutations that appeared over three years in viral genomes. Here we studied mutations and their epidemic consequences for SARS-CoV-2 genomes from our center.

METHODS

We analyzed 61,397 SARS-CoV-2 genomes we sequenced from respiratory samples for genomic surveillance. Mutations frequencies were calculated using Nextclade, Microsoft Excel, and an in-house Python script.

RESULTS

A total of 22,225 nucleotide mutations were identified, 220 (1.0%) being each at the root of ≥836 genomes, classifying mutations as 'hyperfertile'. Two seeded the European pandemic: P323L in RNA polymerase, associated with an increased mutation rate, and D614G in spike that improved fitness. Most 'hyperfertile' mutations occurred in areas not predicted with increased virulence. Their mean number was 8±6 (0-22) per 1000 nucleotides per gene. They were 3.7-times more frequent in accessory than informational genes (13.8 versus 3.7/1000 nucleotides). Particularly, they were 4.1-times more frequent in ORF8 than in the RNA polymerase gene. Interestingly, stop codons were present in 97 positions, almost only in accessory genes, including ORF8 (21/100 codons).

CONCLUSIONS

most 'hyperfertile' mutations did not predict emergence of a new epidemic, and some were stop codons indicating the existence of so-named 'non-virulence' genes.

Spike-Seq: An amplicon-based high-throughput sequencing approach for the sensitive detection and characterization of SARS-CoV-2 genetic variations in environmental samples

Ever since the outbreak of COVID-19 disease in Wuhan, China, different variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been identified. Wastewater-based epidemiology (WBE), an approach that has been successfully applied in numerous case studies worldwide, offers a cost-effective and rapid way for monitoring trends of SARS-Cov-2 in the community level without selection bias. Despite being a gold-standard procedure, WBE is a challenging approach due to the sample instability and the moderate efficiency of SARS-CoV-2 concentration in wastewater. In the present study, we introduce Spike-Seq, a custom amplicon-based approach for the S gene sequencing of SARS-CoV-2 in wastewater samples, which enables not only the accurate identification of the existing Spike-related genetic markers, but also the estimation of their frequency in the investigated samples. The implementation of Spike-Seq involves the combination of nested PCR-based assays that efficiently amplify the entire nucleotide sequence of the S gene and next-generation sequencing, which enables the variant detection and the estimation of their frequency. In the framework of the current work, Spike-Seq was performed to investigate the mutational profile of SARS-CoV-2 in samples from the Wastewater Treatment Plant (WWTP) of Athens, Greece, which originated from multiple timepoints, ranging from March 2021 until July 2022. Our findings demonstrate that Spike-Seq efficiently detected major genetic markers of B.1.1.7 (Alpha), B.1.617.2 (Delta) as well as B.1.1.529 (Omicron) variants in wastewater samples and provided their frequency levels, showing similar variant distributions with the published clinical data from the National Public Health organization. The presented approach can prove to be a useful tool for the detection of SARS-CoV-2 in challenging wastewater samples and the identification of the existing genetic variants of S gene.

Resistance-associated mutations to the anti-SARS-CoV-2 agent nirmatrelvir: Selection not induction

Mutations associated with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) resistance to antiprotease nirmatrelvir were reported. We aimed to detect them in SARS-CoV-2 genomes and quasispecies retrieved in our institute before drug availability in January 2022 and to analyze the impact of mutations on protease (3CLpro) structure. We sought for 38 3CLpro nirmatrelvir resistance mutations in a set of 62 673 SARS-CoV-2 genomes obtained in our institute from respiratory samples collected between 2020 and 2023 and for these mutations in SARS-CoV-2 quasispecies for 90 samples collected in 2020, using Python. SARS-CoV-2 protease with major mutation E166V was generated with Swiss Pdb Viewer and Molegro Molecular Viewer. We detected 22 (58%) of the resistance-associated mutations in 417 (0.67%) of the genomes analyzed; 325 (78%) of these genomes had been obtained from samples collected in 2020-2021. APOBEC signatures were found for 12/22 mutations. We also detected among viral quasispecies from 90 samples some minority reads harboring any of 15 nirmatrelvir resistance mutations, including E166V. Also, we predicted that E166V has a very limited effect on 3CLpro structure but may prevent drug attachment. Thus, we evidenced that mutations associated with nirmatrelvir resistance pre-existed in SARS-CoV-2 before drug availability. These findings further warrant SARS-CoV-2 genomic surveillance and SARS-CoV-2 quasispecies characterization.

Assessing the Potential Role of Cats (Felis catus) as Generators of Relevant SARS-CoV-2 Lineages during the Pandemic

Several questions regarding the evolution of SARS-CoV-2 remain poorly elucidated. One of these questions is the possible evolutionary impact of SARS-CoV-2 after the infection in domestic animals. In this study, we aimed to evaluate the potential role of cats as generators of relevant SARS-CoV-2 lineages during the pandemic. A total of 105 full-length genome viral sequences obtained from naturally infected cats during the pandemic were evaluated by distinct evolutionary algorithms. Analyses were enhanced, including a set of highly related SARS-CoV-2 sequences recovered from human populations. Our results showed the apparent high susceptibility of cats to the infection SARS-CoV-2 compared with other animal species. Evolutionary analyses indicated that the phylogenomic characteristics displayed by cat populations were influenced by the dominance of specific SARS-CoV-2 genetic groups affecting human populations. However, disparate dN/dS rates at some genes between populations recovered from cats and humans suggested that infection in these two species may suggest a different evolutionary constraint for SARS-CoV-2. Interestingly, the branch selection analysis showed evidence of the potential role of natural selection in the emergence of five distinct cat lineages during the pandemic. Although these lineages were apparently irrelevant to public health during the pandemic, our results suggested that additional studies are needed to understand the role of other animal species in the evolution of SARS-CoV-2 during the pandemic.

The return of the "Mistigri" (virus adaptative gain by gene loss) through the SARS-CoV-2 XBB.1.5 chimera that predominated in 2023

Severe acute respiratory syndrome coronavirus 2 XBB.1.5 is the first recombinant lineage to predominate at the country and global scales. Very interestingly, like the Marseille-4B subvariant (or B.1.160) and the pandemic variant B.1.1.7 (or Alpha) previously, it has its ORF8 gene inactivated by a stop codon. We aimed here to study the distribution of stop codons in ORF8 of XBB.1.5 and non-XBB.1.5 genomes. We identified that a stop codon was present at 89 (74%) ORF8 codons in ≥1 of 15 222 404 genomes available in GISAID. The mean proportion of genomes with a stop codon per codon was 0.11% (range, 0%-7.8%). In addition, a stop codon was detected at 15 (12%) codons in at least 1000 genomes. These 15 codons are notably located on seven stem-loop hairpin regions and in the signal peptide region for the case of the XBB.1.5 lineage (codon 8). Thus, it is very likely that stop codons in ORF8 gene contributed on at least three occasions and independently during the pandemic to the evolutionary success of a lineage that became transiently predominant. Such association of gene loss with evolutionary success, which suits the recently described Mistigri rule, is an important biological phenomenon very unknown in virology while largely described in cellular organisms.

Mutations in SARS-CoV-2 variant nsp6 enhance type-I interferon antagonism

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to evolve after its emergence. Given its importance in viral infection and vaccine development, mutations in the viral Spike gene have been studied extensively; however, the impact of mutations outside the Spike gene are poorly understood. Here, we report that a triple deletion (ΔSGF or ΔLSG) in nonstructural protein 6 (nsp6) independently acquired in Alpha and Omicron sublineages of SARS-CoV-2 augments nsp6-mediated antagonism of type-I interferon (IFN-I) signaling. Specifically, these triple deletions enhance the ability of mutant nsp6 to suppress phosphorylation of STAT1 and STAT2. A parental SARS-CoV-2 USA-WA1/2020 strain containing the nsp6 ΔSGF deletion (ΔSGF-WA1) shows reduced susceptibility to IFN-I treatment in vitro, outcompetes the parental strain in human primary airway cultures, and increases virulence in mice; however, the ΔSGF-WA1 virus is less virulent than the Alpha variant (which has the nsp6 ΔSGF deletion and additional mutations in other genes). Analyses of host responses from ΔSGF-WA1-infected mice and primary airway cultures reveal activation of pathways indicative of a cytokine storm. These results provide evidence that mutations outside the Spike protein affect virus-host interactions and may alter pathogenesis of SARS-CoV-2 variants in humans.

Characterisation of SARS-CoV-2 variants in Beijing during 2022: an epidemiological and phylogenetic analysis

BACKGROUND

Due to the national dynamic zero-COVID strategy in China, there were no persistent local transmissions of SARS-CoV-2 in Beijing before December, 2022. However, imported cases have been frequently detected over the past 3 years. With soaring growth in the number of COVID-19 cases in China recently, there are concerns that there might be an emergence of novel SARS-CoV-2 variants. Routine surveillance of viral genomes has been carried out in Beijing over the last 3 years. Spatiotemporal analyses of recent viral genome sequences compared with that of global pooled and local data are crucial for the global response to the ongoing COVID-19 pandemic.

METHODS

We routinely collected respiratory samples covering both imported and local cases in Beijing for the last 3 years (of which the present study pertains to samples collected between January and December, 2022), and then randomly selected samples for analysis. Next-generation sequencing was used to generate the SARS-CoV-2 genomes. Phylogenetic and population dynamic analyses were performed using high-quality complete sequences in this study.

FINDINGS

We obtained a total of 2994 complete SARS-CoV-2 genome sequences in this study, among which 2881 were high quality and were used for further analysis. From Nov 14 to Dec 20, we sequenced 413 new samples, including 350 local cases and 63 imported cases. All of these genomes belong to the existing 123 Pango lineages, showing there are no persistently dominant variants or novel lineages. Nevertheless, BA.5.2 and BF.7 are currently dominant in Beijing, accounting for 90% of local cases since Nov 14 (315 of 350 local cases sequenced in this study). The effective population size for both BA.5.2 and BF.7 in Beijing increased after Nov 14, 2022.

INTERPRETATION

The co-circulation of BF.7 and BA.5.2 has been present in the current outbreak since Nov 14, 2022 in Beijing, and there is no evidence that novel variants emerged. Although our data were only from Beijing, the results could be considered a snapshot of China, due to the frequent population exchange and the presence of circulating strains with high transmissibility.

FUNDING

National Key Research and Development Program of China and Strategic Priority Research Program of the Chinese Academy of Sciences.

TRANSLATION

For the Chinese translation of the abstract see Supplementary Materials section.

Ribosomal control in RNA virus-infected cells

Viruses are strictly intracellular parasites requiring host cellular functions to complete their reproduction cycle involving virus infection of host cell, viral genome replication, viral protein translation, and virion release. Ribosomes are protein synthesis factories in cells, and viruses need to manipulate ribosomes to complete their protein synthesis. Viruses use translation initiation factors through their own RNA structures or cap structures, thereby inducing ribosomes to synthesize viral proteins. Viruses also affect ribosome production and the assembly of mature ribosomes, and regulate the recognition of mRNA by ribosomes, thereby promoting viral protein synthesis and inhibiting the synthesis of host antiviral immune proteins. Here, we review the remarkable mechanisms used by RNA viruses to regulate ribosomes, in particular, the mechanisms by which RNA viruses induce the formation of specific heterogeneous ribosomes required for viral protein translation. This review provides valuable insights into the control of viral infection and diseases from the perspective of viral protein synthesis.

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.