Causes and Consequences of Purifying Selection on SARS-CoV-2

Use of a risk assessment tool to determine the origin of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)

The origin of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is contentious. Most studies have focused on a zoonotic origin, but definitive evidence such as an intermediary animal host is lacking. We used an established risk analysis tool for differentiating natural and unnatural epidemics, the modified Grunow-Finke assessment tool (mGFT) to study the origin of SARS-COV-2. The mGFT scores 11 criteria to provide a likelihood of natural or unnatural origin. Using published literature and publicly available sources of information, we applied the mGFT to the origin of SARS-CoV-2. The mGFT scored 41/60 points (68%), with high inter-rater reliability (100%), indicating a greater likelihood of an unnatural than natural origin of SARS-CoV-2. This risk assessment cannot prove the origin of SARS-CoV-2 but shows that the possibility of a laboratory origin cannot be easily dismissed.

Selection on synonymous sites: the unwanted transcript hypothesis

Although translational selection to favour codons that match the most abundant tRNAs is not readily observed in humans, there is nonetheless selection in humans on synonymous mutations. We hypothesize that much of this synonymous site selection can be explained in terms of protection against unwanted RNAs - spurious transcripts, mis-spliced forms or RNAs derived from transposable elements or viruses. We propose not only that selection on synonymous sites functions to reduce the rate of creation of unwanted transcripts (for example, through selection on exonic splice enhancers and cryptic splice sites) but also that high-GC content (but low-CpG content), together with intron presence and position, is both particular to functional native mRNAs and used to recognize transcripts as native. In support of this hypothesis, transcription, nuclear export, liquid phase condensation and RNA degradation have all recently been shown to promote GC-rich transcripts and suppress AU/CpG-rich ones. With such 'traps' being set against AU/CpG-rich transcripts, the codon usage of native genes has, in turn, evolved to avoid such suppression. That parallel filters against AU/CpG-rich transcripts also affect the endosomal import of RNAs further supports the unwanted transcript hypothesis of synonymous site selection and explains the similar design rules that have enabled the successful use of transgenes and RNA vaccines.

Wastewater-based surveillance as a tool for public health action: SARS-CoV-2 and beyond

Wastewater-based surveillance (WBS) has undergone dramatic advancement in the context of the coronavirus disease 2019 (COVID-19) pandemic. The power and potential of this platform technology were rapidly realized when it became evident that not only did WBS-measured SARS-CoV-2 RNA correlate strongly with COVID-19 clinical disease within monitored populations but also, in fact, it functioned as a leading indicator. Teams from across the globe rapidly innovated novel approaches by which wastewater could be collected from diverse sewersheds ranging from wastewater treatment plants (enabling community-level surveillance) to more granular locations including individual neighborhoods and high-risk buildings such as long-term care facilities (LTCF). Efficient processes enabled SARS-CoV-2 RNA extraction and concentration from the highly dilute wastewater matrix. Molecular and genomic tools to identify, quantify, and characterize SARS-CoV-2 and its various variants were adapted from clinical programs and applied to these mixed environmental systems. Novel data-sharing tools allowed this information to be mobilized and made immediately available to public health and government decision-makers and even the public, enabling evidence-informed decision-making based on local disease dynamics. WBS has since been recognized as a tool of transformative potential, providing near-real-time cost-effective, objective, comprehensive, and inclusive data on the changing prevalence of measured analytes across space and time in populations. However, as a consequence of rapid innovation from hundreds of teams simultaneously, tremendous heterogeneity currently exists in the SARS-CoV-2 WBS literature. This manuscript provides a state-of-the-art review of WBS as established with SARS-CoV-2 and details the current work underway expanding its scope to other infectious disease targets.

Introduction and effect of natural selection analysis at common mutations of SARS-CoV-2 spike gene in Iran

The epidemic of coronavirus disease 2019 (COVID-19) was caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The spike (S) protein of SARS-Cov-2 is composed of two subunits, S1 and S2. This study aimed to describe SARS-CoV-2 haplotypes in Iranians based on the S gene, which plays a key role in the receptor recognition and cell membrane fusion proses. 95 positive saliva samples for SARS-CoV-2 were amplified and sequenced for the S gene. The sequences were classified into 35 haplotypes, which 11 haplotypes were new (H1, H2, H3, H4, H6, H7, H11, H13, H15, H16, H25) and have not been reported so far. Amino acid substitutions were found at 40 positions that 23 were located at S1 subunit and 16 were at S2 subunit and one was at cleavage loop (P681H/R), thus polymorphisms at S1 subunit were found to be higher than S2. The neutrality index (NI) analyses showed a negative departure from the neutral substitution patterns (NI > 1) for S1 and S2 subunit in the studied sequences. The co-occurrence of B-cell epitopes and mutation sites were found in seven positions with more probably to be exposed the immune system pressure. In conclusion, the results provide the significant data to design an effective vaccine based on this protein.

Demographic and Viral-Genetic Analyses of COVID-19 Severity in Bahrain Identify Local Risk Factors and a Protective Effect of Polymerase Mutations

A multitude of demographic, health, and genetic factors are associated with the risk of developing severe COVID-19 following infection by the SARS-CoV-2. There is a need to perform studies across human societies and to investigate the full spectrum of genetic variation of the virus. Using data from 869 COVID-19 patients in Bahrain between March 2020 and March 2021, we analyzed paired viral sequencing and non-genetic host data to understand host and viral determinants of severe COVID-19. We estimated the effects of demographic variables specific to the Bahrain population and found that the impact of health factors are largely consistent with other populations. To extend beyond the common variants of concern in the Spike protein analyzed by previous studies, we used a viral burden approach and detected a protective effect of low-frequency missense viral mutations in the RNA-dependent RNA polymerase (Pol) gene on disease severity. Our results contribute to the survey of severe COVID-19 in diverse populations and highlight the benefits of studying rare viral mutations.

Genome analysis of SARS-CoV-2 isolates from a population reveals the rapid selective sweep of a haplotype carrying many pre-existing and new mutations

To understand the mechanism underlying the evolution of SARS-CoV-2 in a population, we sequenced 92 viral genomes from Assam, India. Analysis of these and database sequences revealed a complete selective sweep of a haplotype in Assam carrying 13 pre-existing variants, including a high leap in frequency of a variant on ORF8, which is involved in immune evasion. A comparative study between sequences of same lineage and similar time frames in and outside Assam showed that 10 of the 13 pre-existing variants had a frequency ranging from 96 to 99%, and the remaining 3 had a low frequency outside Assam. Using a phylogenetic approach to infer sequential occurrences of variants we found that the variant Phe120del on ORF8, which had a low frequency (1.75%) outside Assam, is at the base of the phylogenetic tree of variants and became totally fixed (100%) in Assam population. Based on this observation, we inferred that the variant on ORF8 had a selective advantage, so it carried the haplotype to reach the100% frequency. The haplotype also carried 32 pre-existing variants at a frequency from 1.00 to 80.00% outside Assam. Those of these variants that are more closely linked to the S-protein locus, which often carries advantageous mutations and is tightly linked to the ORF8 locus, retained higher frequencies, while the less tightly linked variants showed lower frequencies, likely due to recombination among co- circulating variants in Assam. The ratios of non-synonymous substitutions to synonymous substitutions suggested that some genes such as those coding for the S-protein and non-structural proteins underwent positive selection while others were subject to purifying selection during their evolution in Assam. Furthermore, we observed negative correlation of the Ct value of qRT-PCR of the patients with abundant ORF6 transcripts, suggesting that ORF6 can be used as a marker for estimating viral titer. In conclusion, our in-depth analysis of SARS-CoV-2 genomes in a regional population reveals the mechanism and dynamics of viral evolution.

Distinguishing COVID-19 from seasonal influenza in patients under age 65 years-a retrospective observational cohort study comparing the 2009 influenza A (H1N1) and 2022 SARS-CoV-2 pandemics

Introduction

This study explored the differences in clinical characteristics between the 2009 pandemic influenza A (H1N1) and SARS-CoV-2 BA.2 variant (Omicron) infections in patients younger than age 65 years, to improve identification of these diseases and better respond to the current epidemic.

Methods

Data from 127 patients with the 2009 pandemic influenza A (H1N1) diagnosed between May and July of 2009 and 3,265 patients with Omicron diagnosed between March and May of 2022 were collected. Using a 1:2 match based on age (difference <2 years), sex, and underlying diseases, data from 115 patients with the 2009 pandemic influenza A (H1N1) infection (H1N1 group) and 230 patients with SARS-CoV-2 Omicron BA.2 infection (Omicron group) were analyzed. The clinical manifestations were compared between the groups, logistic regression was performed to identify possible independent risk factors for each group, and multiple linear regression was used to analyze the factors predicting time for nucleic acid negativization (NAN).

Results

The median [interquartile range] age of the two groups was 21 [11, 26] years. Compared with the H1N1 group, the Omicron group had: lower white blood cell counts and C-reactive protein levels; less fever, nasal congestion, sore throat, cough, sputum, and headache; and more olfactory loss, muscle soreness, and lactate dehydrogenase (LDH) abnormalities. Patients in the Omicron group used fewer antibiotics and antiviral drugs, and the time for NAN was longer (17 [14,20] VS 4 [3,5] days, P<0.001). Logistic regression showed that fever, cough, headache, and increased white blood cell count were more strongly correlated with the H1N1 group, while muscle soreness and LDH abnormalities were more strongly correlated with the Omicron group. Fever (B 1.529, 95% confidence interval [0.149,2.909], P=0.030) significantly predicted a longer time for NAN in patients with Omicron.

Discussion

There are significant differences in clinical characteristics between SARS-CoV-2 Omicron infection and the 2009 pandemic influenza A (H1N1) infection. Recognition of these differences has important implications for clinical practice.

Intra-Host Co-Existing Strains of SARS-CoV-2 Reference Genome Uncovered by Exhaustive Computational Search

The COVID-19 pandemic caused by SARS-CoV-2 has had a severe impact on people worldwide. The reference genome of the virus has been widely used as a template for designing mRNA vaccines to combat the disease. In this study, we present a computational method aimed at identifying co-existing intra-host strains of the virus from RNA-sequencing data of short reads that were used to assemble the original reference genome. Our method consisted of five key steps: extraction of relevant reads, error correction for the reads, identification of within-host diversity, phylogenetic study, and protein binding affinity analysis. Our study revealed that multiple strains of SARS-CoV-2 can coexist in both the viral sample used to produce the reference sequence and a wastewater sample from California. Additionally, our workflow demonstrated its capability to identify within-host diversity in foot-and-mouth disease virus (FMDV). Through our research, we were able to shed light on the binding affinity and phylogenetic relationships of these strains with the published SARS-CoV-2 reference genome, SARS-CoV, variants of concern (VOC) of SARS-CoV-2, and some closely related coronaviruses. These insights have important implications for future research efforts aimed at identifying within-host diversity, understanding the evolution and spread of these viruses, as well as the development of effective treatments and vaccines against them.

Developing an appropriate evolutionary baseline model for the study of SARS-CoV-2 patient samples

Over the past 3 years, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spread through human populations in several waves, resulting in a global health crisis. In response, genomic surveillance efforts have proliferated in the hopes of tracking and anticipating the evolution of this virus, resulting in millions of patient isolates now being available in public databases. Yet, while there is a tremendous focus on identifying newly emerging adaptive viral variants, this quantification is far from trivial. Specifically, multiple co-occurring and interacting evolutionary processes are constantly in operation and must be jointly considered and modeled in order to perform accurate inference. We here outline critical individual components of such an evolutionary baseline model-mutation rates, recombination rates, the distribution of fitness effects, infection dynamics, and compartmentalization-and describe the current state of knowledge pertaining to the related parameters of each in SARS-CoV-2. We close with a series of recommendations for future clinical sampling, model construction, and statistical analysis.

Immunosuppression as a Hub for SARS-CoV-2 Mutational Drift

The clinical course of coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), is largely determined by host factors, with a wide range of outcomes. Despite an extensive vaccination campaign and high rates of infection worldwide, the pandemic persists, adapting to overcome antiviral immunity acquired through prior exposure. The source of many such major adaptations is variants of concern (VOCs), novel SARS-CoV-2 variants produced by extraordinary evolutionary leaps whose origins remain mostly unknown. In this study, we tested the influence of factors on the evolutionary course of SARS-CoV-2. Electronic health records of individuals infected with SARS-CoV-2 were paired to viral whole-genome sequences to assess the effects of host clinical parameters and immunity on the intra-host evolution of SARS-CoV-2. We found slight, albeit significant, differences in SARS-CoV-2 intra-host diversity, which depended on host parameters such as vaccination status and smoking. Only one viral genome had significant alterations as a result of host parameters; it was found in an immunocompromised, chronically infected woman in her 70s. We highlight the unusual viral genome obtained from this woman, which had an accelerated mutational rate and an excess of rare mutations, including near-complete truncating of the accessory protein ORF3a. Our findings suggest that the evolutionary capacity of SARS-CoV-2 during acute infection is limited and mostly unaffected by host characteristics. Significant viral evolution is seemingly exclusive to a small subset of COVID-19 cases, which typically prolong infections in immunocompromised patients. In these rare cases, SARS-CoV-2 genomes accumulate many impactful and potentially adaptive mutations; however, the transmissibility of such viruses remains unclear.

Variation in synonymous evolutionary rates in the SARS-CoV-2 genome

Introduction

Coronavirus disease 2019 is an infectious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Influential variants and mutants of this virus continue to emerge, and more effective virus-related information is urgently required for identifying and predicting new mutants. According to earlier reports, synonymous substitutions were considered phenotypically silent; thus, such mutations were frequently ignored in studies of viral mutations because they did not directly cause amino acid changes. However, recent studies have shown that synonymous substitutions are not completely silent, and their patterns and potential functional correlations should thus be delineated for better control of the pandemic.

Methods

In this study, we estimated the synonymous evolutionary rate (SER) across the SARS-CoV-2 genome and used it to infer the relationship between the viral RNA and host protein. We also assessed the patterns of characteristic mutations found in different viral lineages.

Results

We found that the SER varies across the genome and that the variation is primarily influenced by codon-related factors. Moreover, the conserved motifs identified based on the SER were found to be related to host RNA transport and regulation. Importantly, the majority of the existing fixed-characteristic mutations for five important virus lineages (Alpha, Beta, Gamma, Delta, and Omicron) were significantly enriched in partially constrained regions.

Discussion

Taken together, our results provide unique information on the evolutionary and functional dynamics of SARS-CoV-2 based on synonymous mutations and offer potentially useful information for better control of the SARS-CoV-2 pandemic.

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

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

Attenuation and Degeneration of SARS-CoV-2 Despite Adaptive Evolution

The evolution of severe acute respiratory syndrome coronavirus 2 (SARS‑CoV‑2) has followed similar trends as other RNA viruses, such as human immunodeficiency virus type 1 and the influenza A virus. Rapid initial diversification was followed by strong competition and a rapid succession of dominant variants. Host-initiated RNA editing has been the primary mechanism for introducing mutations. A significant number of mutations detrimental to viral replication have been quickly purged. Fixed mutations are mostly diversifying mutations selected for host adaptation and immune evasion, with the latter accounting for the majority of the mutations. However, immune evasion often comes at the cost of functionality, and thus, optimal functionality is still far from being accomplished. Instead, selection for antibody-escaping variants and accumulation of near-neutral mutations have led to suboptimal codon usage and reduced replicative capacity, as demonstrated in non-respiratory cell lines. Beneficial adaptation of the virus includes reduced infectivity in lung tissues and increased tropism for the upper airway, resulting in shorter incubation periods, milder diseases, and more efficient transmission between people.

Mutational Patterns Observed in SARS-CoV-2 Genomes Sampled From Successive Epochs Delimited by Major Public Health Events in Ontario, Canada: Genomic Surveillance Study

BACKGROUND

The emergence of SARS-CoV-2 variants with mutations associated with increased transmissibility and virulence is a public health concern in Ontario, Canada. Characterizing how the mutational patterns of the SARS-CoV-2 genome have changed over time can shed light on the driving factors, including selection for increased fitness and host immune response, that may contribute to the emergence of novel variants. Moreover, the study of SARS-CoV-2 in the microcosm of Ontario, Canada can reveal how different province-specific public health policies over time may be associated with observed mutational patterns as a model system.

OBJECTIVE

This study aimed to perform a comprehensive analysis of single base substitution (SBS) types, counts, and genomic locations observed in SARS-CoV-2 genomic sequences sampled in Ontario, Canada. Comparisons of mutational patterns were conducted between sequences sampled during 4 different epochs delimited by major public health events to track the evolution of the SARS-CoV-2 mutational landscape over 2 years.

METHODS

In total, 24,244 SARS-CoV-2 genomic sequences and associated metadata sampled in Ontario, Canada from January 1, 2020, to December 31, 2021, were retrieved from the Global Initiative on Sharing All Influenza Data database. Sequences were assigned to 4 epochs delimited by major public health events based on the sampling date. SBSs from each SARS-CoV-2 sequence were identified relative to the MN996528.1 reference genome. Catalogues of SBS types and counts were generated to estimate the impact of selection in each open reading frame, and identify mutation clusters. The estimation of mutational fitness over time was performed using the Augur pipeline.

RESULTS

The biases in SBS types and proportions observed support previous reports of host antiviral defense activity involving the SARS-CoV-2 genome. There was an increase in U>C substitutions associated with adenosine deaminase acting on RNA (ADAR) activity uniquely observed during Epoch 4. The burden of novel SBSs observed in SARS-CoV-2 genomic sequences was the greatest in Epoch 2 (median 5), followed by Epoch 3 (median 4). Clusters of SBSs were observed in the spike protein open reading frame, ORF1a, and ORF3a. The high proportion of nonsynonymous SBSs and increasing dN/dS metric (ratio of nonsynonymous to synonymous mutations in a given open reading frame) to above 1 in Epoch 4 indicate positive selection of the spike protein open reading frame.

CONCLUSIONS

Quantitative analysis of the mutational patterns of the SARS-CoV-2 genome in the microcosm of Ontario, Canada within early consecutive epochs of the pandemic tracked the mutational dynamics in the context of public health events that instigate significant shifts in selection and mutagenesis. Continued genomic surveillance of emergent variants will be useful for the design of public health policies in response to the evolving COVID-19 pandemic.

Monoclonal antibodies for the treatment of COVID-19 infection in children

INTRODUCTION

Monoclonal antibodies (mAbs) have been authorized for the treatment of COVID-19 in pediatric populations, however, there is a lack of evidence for their use in these populations.

AREAS COVERED

We outline the evidence of mAbs for COVID-19, discuss their use in the treatment of COVID-19 infection for pediatric patients, and consider alternative treatment options and challenges to COVID-19 drug approvals.

EXPERT OPINION

Limited evidence exists for the safety and efficacy of mAbs to treat COVID-19 in children as new variants emerge. In rare pediatric outpatient settings, such as profound immunodeficiency or severe pulmonary disease, the benefits of antiviral treatment for COVID-19 likely outweigh the relatively small risks. However, for the great majority of pediatric patients, mAb treatment is likely not indicated. Small molecule antiviral therapies are another potential treatment for COVID-19 in children in an outpatient setting, though neither mAb nor small molecule antiviral treatments have significant supporting evidence in children and developing a strong evidence base for these decisions will be challenging if not impractical. Ultimately, these decisions are likely to be made at the level of individual cases using expert opinion as the primary guiding principle.

Identification of mutations in SARS-CoV-2 PCR primer regions

Due to the constantly increasing number of mutations in the SARS-CoV-2 genome, concerns have emerged over the possibility of decreased diagnostic accuracy of reverse transcription-polymerase chain reaction (RT-PCR), the gold standard diagnostic test for SARS-CoV-2. We propose an analysis pipeline to discover genomic variations overlapping the target regions of commonly used PCR primer sets. We provide the list of these mutations in a publicly available format based on a dataset of more than 1.2 million SARS-CoV-2 samples. Our approach distinguishes among mutations possibly having a damaging impact on PCR efficiency and ones anticipated to be neutral in this sense. Samples are categorized as "prone to misclassification" vs. "likely to be correctly detected" by a given PCR primer set based on the estimated effect of mutations present. Samples susceptible to misclassification are generally present at a daily rate of 2% or lower, although particular primer sets seem to have compromised performance when detecting Omicron samples. As different variant strains may temporarily gain dominance in the worldwide SARS-CoV-2 viral population, the efficiency of a particular PCR primer set may change over time, therefore constant monitoring of variations in primer target regions is highly recommended.

The health benefit of physical exercise on COVID-19 pandemic: Evidence from mainland China

Objectives

Our study aims to investigate the health benefit of regular physical exercise participation on a series of COVID-19 outcomes including COVID-19 morbidity, mortality, and cure rate.

Methods

Prefecture-level panel data related to physical exercise and the COVID-19 pandemic in China were collected from January 1 to March 17, 2020, (N = 21379). Multiple linear regression was conducted, and the ordinary least squares technique was used to estimate the coefficient.

Results

It was shown that regular sports participation significantly negatively affected COVID-19 morbidity (estimate = -1.1061, p<0.01) and mortality (estimate = -0.3836, p<0.01), and positively affected cure rate (estimate = 0.0448, p<0.01), implying that engaging in physical exercise regularly does have a significant positive effect on COVID-19 outcomes. Then, we explored the heterogeneity of the effect of physical exercise on areas with different risk levels and it was revealed that the effect of physical exercise was more pronounced in high-risk areas in terms of morbidity (estimate = -1.8776, p<0.01 in high-risk areas; estimate = -0.0037, p<0.01 in low-risk areas), mortality (estimate = -0.3982, p<0.01 in high-risk areas; estimate = -0.3492, p<0.01 in low-risk areas), and cure rate (estimate = 0.0807, p<0.01 in high-risk areas; 0.0193 = -0.0037, p<0.05 in low-risk areas).

Conclusions

Our results suggest that regularly engaging in physical exercise before the pandemic has positive health effects, especially in the case of a more severe epidemic. Therefore, we urge readers to actively engage in physical exercise so that we can reduce the risks in the event of a pandemic.

SARS-CoV-2 omicron variant: the black swan of microbiology

The sudden emergence of the SARS-CoV-2 Omicron variant is causing major global concern due to its high number of mutations compared to previous variants, which is a relatively rare but significant event that can change the course of viral evolution, the occurrence of which might have huge consequences for the natural evolution of species in general, prompting us to rethink our knowledge on evolution.

The current status of gene expression profilings in COVID‐19 patients

Background

The global pandemic of coronavirus disease 2019 (COVID‐19) caused by severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) has swept through every part of the world. Because of its impact, international efforts have been underway to identify the variants of SARS‐CoV‐2 by genome sequencing and to understand the gene expression changes in COVID‐19 patients compared to healthy donors using RNA sequencing (RNA‐seq) assay. Within the last two and half years since the emergence of SARS‐CoV‐2, a large number of OMICS data of COVID‐19 patients have accumulated. Yet, we are still far from understanding the disease mechanism. Further, many people suffer from long‐term effects of COVID‐19; calling for a more systematic way to data mine the generated OMICS data, especially RNA‐seq data.

Methods

By searching gene expression omnibus (GEO) using the key terms, COVID‐19 and RNA‐seq, 108 GEO entries were identified. Each of these studies was manually examined to categorize the studies into bulk or single‐cell RNA‐seq (scRNA‐seq) followed by an inspection of their original articles.

Results

The currently available RNA‐seq data were generated from various types of patients’ samples, and COVID‐19 related sample materials have been sequenced at the level of RNA, including whole blood, different components of blood [e.g., plasma, peripheral blood mononuclear cells (PBMCs), leukocytes, lymphocytes, monocytes, T cells], nasal swabs, and autopsy samples (e.g., lung, heart, liver, kidney). Of these, RNA‐seq studies using whole blood, PBMCs, nasal swabs and autopsy/biopsy samples were reviewed to highlight the major findings from RNA‐seq data analysis.

Conclusions

Based on the bulk and scRNA‐seq data analysis, severe COVID‐19 patients display shifts in cell populations, especially those of leukocytes and monocytes, possibly leading to cytokine storms and immune silence. These RNA‐seq data form the foundation for further gene expression analysis using samples from individuals suffering from long COVID.

RNA polymerase inaccuracy underlies SARS-CoV-2 variants and vaccine heterogeneity

Both the SARS-CoV-2 virus and its mRNA vaccines depend on RNA polymerases (RNAP)1,2; however, these enzymes are inherently error-prone and can introduce variants into the RNA3. To understand SARS-CoV-2 evolution and vaccine efficacy, it is critical to identify the extent and distribution of errors introduced by the RNAPs involved in each process. Current methods lack the sensitivity and specificity to measure de novo RNA variants in low input samples like viral isolates3. Here, we determine the frequency and nature of RNA errors in both SARS-CoV-2 and its vaccine using a targeted Accurate RNA Consensus sequencing method (tARC-seq). We found that the viral RNA-dependent RNAP (RdRp) makes ~1 error every 10,000 nucleotides - higher than previous estimates4. We also observed that RNA variants are not randomly distributed across the genome but are associated with certain genomic features and genes, such as S (Spike). tARC-seq captured a number of large insertions, deletions and complex mutations that can be modeled through non-programmed RdRp template switching. This template switching feature of RdRp explains many key genetic changes observed during the evolution of different lineages worldwide, including Omicron. Further sequencing of the Pfizer-BioNTech COVID-19 vaccine revealed an RNA variant frequency of ~1 in 5,000, meaning most of the vaccine transcripts produced in vitro by T7 phage RNAP harbor a variant. These results demonstrate the extraordinary genetic diversity of viral populations and the heterogeneous nature of an mRNA vaccine fueled by RNAP inaccuracy. Along with functional studies and pandemic data, tARC-seq variant spectra can inform models to predict how SARS-CoV-2 may evolve. Finally, our results may help improve future vaccine development and study design as mRNA therapies continue to gain traction.

How to Deal With Vaccine Breakthrough Infection With SARS-CoV-2 Variants

Novel Coronary Pneumonia is the most infectious disease with the highest number of morbidity and mortality in 100 years. Despite aggressive and effective COVID-19 prevention and control measures, countries have been unable to stop its outbreaks. With the widespread use of vaccines, the occurrence of COVID-19 has declined markedly. April 21, 2021, New York scholars reported Vaccine Breakthrough Infections with SARS-CoV-2 Variants, which immediately attracted widespread attention. In this mini-review, we focus on the characteristics of SARS-CoV-2 and its mutant strains and vaccine breakthrough infections. We have found that outbreaks of vaccine-breaking SARS-CoV-2 Delta infections in many countries are primarily the result of declining vaccine-generated antibody titers and relaxed outbreak management measures. For this reason, we believe that the main response to vaccine-breaking infections with the SARS-CoV-2 variant is to implement a rigorous outbreak defense policy and vaccine application. Only by intensifying the current vaccination intensity, gradually improving the vaccine and its application methods, and strengthening non-pharmaceutical measures such as travel restrictions, social distancing, masking and hand hygiene, can the COVID-19 outbreak be fully controlled at an early date.

Emergence of SARS-CoV-2 Variants in the World: How Could This Happen?

The COVID-19 pandemic has had a significant global impact, with more than 280,000,000 people infected and 5,400,000 deaths. The use of personal protective equipment and the anti-SARS-CoV-2 vaccination campaigns have reduced infection and death rates worldwide. However, a recent increase in infection rates has been observed associated with the appearance of SARS-CoV-2 variants, including the more recently described lineage B.1.617.2 (Delta variant) and lineage B.1.1.529/BA.1 (Omicron variant). These new variants put the effectiveness of international vaccination at risk, with the appearance of new outbreaks of COVID-19 throughout the world. This emergence of new variants has been due to multiple predisposing factors, including molecular characteristics of the virus, geographic and environmental conditions, and the impact of social determinants of health that favor the genetic diversification of SARS-CoV-2. We present a literature review on the most recent information available on the emergence of new variants of SARS-CoV-2 in the world. We analyzed the biological, geographical, and sociocultural factors that favor the development of these variants. Finally, we evaluate the surveillance strategies for the early detection of new variants and prevent their distribution outside these regions.

Aggregation of high-frequency RBD mutations of SARS-CoV-2 with three VOCs did not cause significant antigenic drift

Variants of SARS‐CoV‐2 continue to emerge, posing great challenges in outbreak prevention and control. It is important to understand in advance the impact of possible variants of concern (VOCs) on infectivity and antigenicity. Here, we constructed one or more of the 15 high‐frequency naturally occurring amino acid changes in the receptor‐binding domain (RBD) of Alpha, Beta, and Gamma variants. A single mutant of A520S, V367F, and S494P in the above three VOCs enhanced infectivity in ACE2‐overexpressing 293T cells of different species, LLC‐MK2 and Vero cells. Aggregation of multiple RBD mutations significantly reduces the infectivity of the possible three VOCs. Regarding neutralization, it is noteworthy that E484K, N501Y, K417N, and N439K predispose to monoclonal antibodies (mAbs) protection failure in the 15 high‐frequency mutations. Most importantly, almost all possible VOCs (single RBD mutation or aggregation of multiple mutations) showed no more than a fourfold decrease in neutralizing activity with convalescent sera, vaccine sera, and immune sera of guinea pigs with different immunogens, and no significant antigenic drift was formed. In conclusion, our pseudovirus results could reduce the concern that the aggregation of multiple high‐frequency mutations in the RBD of the spike protein of the three VOCs would lead to severe antigenic drift, and this would provide value for vaccine development strategies.

Infectivity increased by adding three VOCs of V367F, S494P, or A520S.

The infectivity of the three VOCs with multiple high‐frequency mutations decreased.

Almost all of the possible variants of the three VOCs did not show severe antigenic drift.

A global analysis of conservative and non-conservative mutations in SARS-CoV-2 detected in the first year of the COVID-19 world-wide diffusion

The ability of SARS-CoV-2 to rapidly mutate represents a remarkable complicancy. Quantitative evaluations of the effects that these mutations have on the virus structure/function is of great relevance and the availability of a large number of SARS-CoV-2 sequences since the early phases of the pandemic represents a unique opportunity to follow the adaptation of the virus to humans. Here, we evaluated the SARS-CoV-2 amino acid mutations and their progression by analyzing publicly available viral genomes at three stages of the pandemic (2020 March 15th and October 7th, 2021 February 7th). Mutations were classified in conservative and non-conservative based on the probability to be accepted during the evolution according to the Point Accepted Mutation substitution matrices and on the similarity of the encoding codons. We found that the most frequent substitutions are T > I, L > F, and A > V and we observe accumulation of hydrophobic residues. These findings are consistent among the three stages analyzed. We also found that non-conservative mutations are less frequent than conservative ones. This finding may be ascribed to a progressive adaptation of the virus to the host. In conclusion, the present study provides indications of the early evolution of the virus and tools for the global and genome-specific evaluation of the possible impact of mutations on the structure/function of SARS-CoV-2 variants.