Attenuation of replication by a 29 nucleotide deletion in SARS-coronavirus acquired during the early stages of human-to-human transmission

A glimpse into viral warfare: decoding the intriguing role of highly pathogenic coronavirus proteins in apoptosis regulation

Coronaviruses employ various strategies for survival, among which the activation of endogenous or exogenous apoptosis stands out, with viral proteins playing a pivotal role. Notably, highly pathogenic coronaviruses such as SARS-CoV-2, SARS-CoV, and MERS-CoV exhibit a greater array of non-structural proteins compared to low-pathogenic strains, facilitating their ability to induce apoptosis via multiple pathways. Moreover, these viral proteins are adept at dampening host immune responses, thereby bolstering viral replication and persistence. This review delves into the intricate interplay between highly pathogenic coronaviruses and apoptosis, systematically elucidating the molecular mechanisms underpinning apoptosis induction by viral proteins. Furthermore, it explores the potential therapeutic avenues stemming from apoptosis inhibition as antiviral agents and the utilization of apoptosis-inducing viral proteins as therapeutic modalities. These insights not only shed light on viral pathogenesis but also offer novel perspectives for cancer therapy.

SARS-CoV-2 ORF8 modulates lung inflammation and clinical disease progression

The virus severe acute respiratory syndrome coronavirus 2, SARS-CoV-2, is the causative agent of the current COVID-19 pandemic. It possesses a large 30 kilobase (kb) genome that encodes structural, non-structural, and accessory proteins. Although not necessary to cause disease, these accessory proteins are known to influence viral replication and pathogenesis. Through the synthesis of novel infectious clones of SARS-CoV-2 that lack one or more of the accessory proteins of the virus, we have found that one of these accessory proteins, ORF8, is critical for the modulation of the host inflammatory response. Mice infected with a SARS-CoV-2 virus lacking ORF8 exhibit increased weight loss and exacerbated macrophage infiltration into the lungs. Additionally, infection of mice with recombinant SARS-CoV-2 viruses encoding ORF8 mutations found in variants of concern reveal that naturally occurring mutations in this protein influence disease severity. Our studies with a virus lacking this ORF8 protein and viruses possessing naturally occurring point mutations in this protein demonstrate that this protein impacts pathogenesis.

Properties and Mechanisms of Deletions, Insertions, and Substitutions in the Evolutionary History of SARS-CoV-2

SARS-CoV-2 has accumulated many mutations since its emergence in late 2019. Nucleotide substitutions leading to amino acid replacements constitute the primary material for natural selection. Insertions, deletions, and substitutions appear to be critical for coronavirus's macro- and microevolution. Understanding the molecular mechanisms of mutations in the mutational hotspots (positions, loci with recurrent mutations, and nucleotide context) is important for disentangling roles of mutagenesis and selection. In the SARS-CoV-2 genome, deletions and insertions are frequently associated with repetitive sequences, whereas C>U substitutions are often surrounded by nucleotides resembling the APOBEC mutable motifs. We describe various approaches to mutation spectra analyses, including the context features of RNAs that are likely to be involved in the generation of recurrent mutations. We also discuss the interplay between mutations and natural selection as a complex evolutionary trend. The substantial variability and complexity of pipelines for the reconstruction of mutations and the huge number of genomic sequences are major problems for the analyses of mutations in the SARS-CoV-2 genome. As a solution, we advocate for the development of a centralized database of predicted mutations, which needs to be updated on a regular basis.

Hidden evolutionary constraints dictate the retention of coronavirus accessory genes

Coronaviruses exhibit many mechanisms of genetic innovation (1-5), including the acquisition of accessory genes that originate by capture of cellular genes or through duplication of existing viral genes (6,7). Accessory genes influence viral host range and cellular tropism, but little is known about how selection acts on these variable regions of virus genomes. We used experimental evolution of mouse hepatitis virus (MHV) encoding a cellular AKAP7 phosphodiesterase and an inactive native phosphodiesterase, NS2 (ref 8) to simulate the capture of a host gene and analyze its evolution. After courses of serial infection, the gene encoding inactive NS2, ORF2, unexpectedly remained intact, suggesting it is under cryptic constraint uncoupled from the function of NS2. In contrast, AKAP7 was retained under strong selection but rapidly lost under relaxed selection. Guided by the retention of ORF2 and similar patterns in related betacoronaviruses, we analyzed ORF8 of SARS-CoV-2, which arose via gene duplication6 and contains premature stop codons in several globally successful lineages. As with MHV ORF2, the coding-defective SARS-CoV-2 ORF8 gene remains largely intact, mirroring patterns observed during MHV experimental evolution, challenging assumptions on the dynamics of gene loss in virus genomes and extending these findings to viruses currently adapting to humans.

Unraveling the genetic variations underlying virulence disparities among SARS-CoV-2 strains across global regions: insights from Pakistan

Over the course of the COVID-19 pandemic, several SARS-CoV-2 variants have emerged that may exhibit different etiological effects such as enhanced transmissibility and infectivity. However, genetic variations that reduce virulence and deteriorate viral fitness have not yet been thoroughly investigated. The present study sought to evaluate the effects of viral genetic makeup on COVID-19 epidemiology in Pakistan, where the infectivity and mortality rate was comparatively lower than other countries during the first pandemic wave. For this purpose, we focused on the comparative analyses of 7096 amino-acid long polyprotein pp1ab. Comparative sequence analysis of 203 SARS-CoV-2 genomes, sampled from Pakistan during the first wave of the pandemic revealed 179 amino acid substitutions in pp1ab. Within this set, 38 substitutions were identified within the Nsp3 region of the pp1ab polyprotein. Structural and biophysical analysis of proteins revealed that amino acid variations within Nsp3's macrodomains induced conformational changes and modified protein-ligand interactions, consequently diminishing the virulence and fitness of SARS-CoV-2. Additionally, the epistatic effects resulting from evolutionary substitutions in SARS-CoV-2 proteins may have unnoticed implications for reducing disease burden. In light of these findings, further characterization of such deleterious SARS-CoV-2 mutations will not only aid in identifying potential therapeutic targets but will also provide a roadmap for maintaining vigilance against the genetic variability of diverse SARS-CoV-2 strains circulating globally. Furthermore, these insights empower us to more effectively manage and respond to potential viral-based pandemic outbreaks of a similar nature in the future.

SARS-CoV-2 and innate immunity: the good, the bad, and the "goldilocks"

An ancient conflict between hosts and pathogens has driven the innate and adaptive arms of immunity. Knowledge about this interplay can not only help us identify biological mechanisms but also reveal pathogen vulnerabilities that can be leveraged therapeutically. The humoral response to SARS-CoV-2 infection has been the focus of intense research, and the role of the innate immune system has received significantly less attention. Here, we review current knowledge of the innate immune response to SARS-CoV-2 infection and the various means SARS-CoV-2 employs to evade innate defense systems. We also consider the role of innate immunity in SARS-CoV-2 vaccines and in the phenomenon of long COVID.

Comparative Pathogenesis of Severe Acute Respiratory Syndrome Coronaviruses

Over the last two decades the world has witnessed the global spread of two genetically related highly pathogenic coronaviruses, severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2. However, the impact of these outbreaks differed significantly with respect to the hospitalizations and fatalities seen worldwide. While many studies have been performed recently on SARS-CoV-2, a comparative pathogenesis analysis with SARS-CoV may further provide critical insights into the mechanisms of disease that drive coronavirus-induced respiratory disease. In this review, we comprehensively describe clinical and experimental observations related to transmission and pathogenesis of SARS-CoV-2 in comparison with SARS-CoV, focusing on human, animal, and in vitro studies. By deciphering the similarities and disparities of SARS-CoV and SARS-CoV-2, in terms of transmission and pathogenesis mechanisms, we offer insights into the divergent characteristics of these two viruses. This information may also be relevant to assessing potential novel introductions of genetically related highly pathogenic coronaviruses.

Evolutionary deletions within the SARS-CoV-2 genome as signature trends for virus fitness and adaptation

Coronaviruses are large RNA viruses that can infect and spread among humans and animals. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for coronavirus disease 2019, has evolved since its first detection in December 2019. Deletions are a common occurrence in SARS-CoV-2 evolution, particularly in specific genomic sites, and may be associated with the emergence of highly competent lineages. While deletions typically have a negative impact on viral fitness, some persist and become fixed in viral populations, indicating that they may confer advantageous benefits for the virus's adaptive evolution. This work presents a literature review and data analysis on structural losses in the SARS-CoV-2 genome and the potential relevance of specific signatures for enhanced viral fitness and spread.

Distinct phenotype of SARS-CoV-2 Omicron BA.1 in human primary cells but no increased host range in cell lines of putative mammalian reservoir species

SARS-CoV-2's genetic plasticity has led to several variants of concern (VOCs). Here we studied replicative capacity for seven SARS-CoV-2 isolates (B.1, Alpha, Beta, Gamma, Delta, Zeta, and Omicron BA.1) in primary reconstituted airway epithelia (HAE) and lung-derived cell lines. Furthermore, to investigate the host range of Delta and Omicron compared to ancestral SARS-CoV-2, we assessed replication in 17 cell lines from 11 non-primate mammalian species, including bats, rodents, insectivores and carnivores. Only Omicron's phenotype differed in vitro, with rapid but short replication and efficient production of infectious virus in nasal HAEs, in contrast to other VOCs, but not in lung cell lines. No increased infection efficiency for other species was observed, but Delta and Omicron infection efficiency was increased in A549 cells. Notably replication in A549 and Calu3 cells was lower than in nasal HAE. Our results suggest better adaptation of VOCs towards humans, without an extended host range, and may be relevant to the search for the putative intermediate host and reservoirs prior to the pandemic.

ORF3c is expressed in SARS-CoV-2-infected cells and inhibits innate sensing by targeting MAVS

Most SARS-CoV-2 proteins are translated from subgenomic RNAs (sgRNAs). While the majority of these sgRNAs are monocistronic, some viral mRNAs encode more than one protein. One example is the ORF3a sgRNA that also encodes ORF3c, an enigmatic 41-amino-acid peptide. Here, we show that ORF3c is expressed in SARS-CoV-2-infected cells and suppresses RIG-I- and MDA5-mediated IFN-β induction. ORF3c interacts with the signaling adaptor MAVS, induces its C-terminal cleavage, and inhibits the interaction of RIG-I with MAVS. The immunosuppressive activity of ORF3c is conserved among members of the subgenus sarbecovirus, including SARS-CoV and coronaviruses isolated from bats. Notably, however, the SARS-CoV-2 delta and kappa variants harbor premature stop codons in ORF3c, demonstrating that this reading frame is not essential for efficient viral replication in vivo and is likely compensated by other viral proteins. In agreement with this, disruption of ORF3c does not significantly affect SARS-CoV-2 replication in CaCo-2, CaLu-3, or Rhinolophus alcyone cells. In summary, we here identify ORF3c as an immune evasion factor of SARS-CoV-2 that suppresses innate sensing in infected cells.

The 29-nucleotide deletion in SARS-CoV: truncated versions of ORF8 are under purifying selection

BACKGROUND

Accessory proteins have diverse roles in coronavirus pathobiology. One of them in SARS-CoV (the causative agent of the severe acute respiratory syndrome outbreak in 2002-2003) is encoded by the open reading frame 8 (ORF8). Among the most dramatic genomic changes observed in SARS-CoV isolated from patients during the peak of the pandemic in 2003 was the acquisition of a characteristic 29-nucleotide deletion in ORF8. This deletion cause splitting of ORF8 into two smaller ORFs, namely ORF8a and ORF8b. Functional consequences of this event are not entirely clear.

RESULTS

Here, we performed evolutionary analyses of ORF8a and ORF8b genes and documented that in both cases the frequency of synonymous mutations was greater than that of nonsynonymous ones. These results suggest that ORF8a and ORF8b are under purifying selection, thus proteins translated from these ORFs are likely to be functionally important. Comparisons with several other SARS-CoV genes revealed that another accessory gene, ORF7a, has a similar ratio of nonsynonymous to synonymous mutations suggesting that ORF8a, ORF8b, and ORF7a are under similar selection pressure.

CONCLUSIONS

Our results for SARS-CoV echo the known excess of deletions in the ORF7a-ORF7b-ORF8 complex of accessory genes in SARS-CoV-2. A high frequency of deletions in this gene complex might reflect recurrent searches in "functional space" of various accessory protein combinations that may eventually produce more advantageous configurations of accessory proteins similar to the fixed deletion in the SARS-CoV ORF8 gene.

Sarbecoviruses of British horseshoe bats; sequence variation and epidemiology

Horseshoe bats are the natural hosts of the Sarbecovirus subgenus that includes SARS-CoV and SARS-CoV- 2. Despite the devastating impact of the COVID-19 pandemic, there is still little known about the underlying epidemiology and virology of sarbecoviruses in their natural hosts, leaving large gaps in our pandemic preparedness. Here we describe the results of PCR testing for sarbecoviruses in the two horseshoe bat species (Rhinolophus hipposideros and R. ferrumequinum) present in Great Britain, collected in 2021-22 during the peak of COVID-19 pandemic. One hundred and ninety seven R. hipposideros samples from 33 roost sites and 277 R. ferrumequinum samples from 20 roost sites were tested. No coronaviruses were detected in any samples from R. ferrumequinum whereas 44 and 56 % of individual and pooled (respectively) faecal samples from R. hipposideros across multiple roost sites tested positive in a sarbecovirus-specific qPCR. Full genome sequences were generated from three of the positive samples (and partial genomes from two more) using Illumina RNAseq on unenriched samples. Phylogenetic analyses showed that the obtained sequences belong to the same monophyletic clade, with >95 % similarity to previously-reported European isolates from R. hipposideros. The sequences differed in the presence or absence of accessory genes ORF 7b, 9b and 10. All lacked the furin cleavage site of SARS-CoV-2 spike gene and are therefore unlikely to be infective for humans. These results demonstrate a lack, or at least low incidence, of SARS-CoV-2 spill over from humans to susceptible GB bats, and confirm that sarbecovirus infection is widespread in R. hipposideros. Despite frequently sharing roost sites with R. ferrumequinum, no evidence of cross-species transmission was found.

The relevant information about the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) using the five-question approach (when, where, what, why, and how) and its impact on the environment

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is regarded as a threat because it spreads quickly across the world without requiring a passport or establishing an identity. This tiny virus has wreaked havoc on people's lives, killed people, and created psychological problems all over the world. The viral spike protein (S) significantly contributes to host cell entry, and mutations associated with it, particularly in the receptor-binding protein (RBD), either facilitate the escape of virus from neutralizing antibodies or enhance its transmission by increasing the affinity for cell entry receptor, angiotensin-converting enzyme 2 (ACE2). The initial variants identified in Brazil, South Africa, and the UK have spread to various countries. On the other hand, new variants are being detected in India and the USA. The viral genome and proteome were applied for molecular detection techniques, and nanotechnology particles and materials were utilized in protection and prevention strategies. Consequently, the SARS-CoV-2 pandemic has resulted in extraordinary scientific community efforts to develop detection methods, diagnosis tools, and effective antiviral drugs and vaccines, where prevailing academic, governmental, and industrial institutions and organizations continue to engage themselves in large-scale screening of existing drugs, both in vitro and in vivo. In addition, COVID-19 pointed on the possible solutions for the environmental pollution globe problem. Therefore, this review aims to address SARS-CoV-2, its transmission, where it can be found, why it is severe in some people, how it can be stopped, its diagnosis and detection techniques, and its relationship with the environment.

ViReMa: a virus recombination mapper of next-generation sequencing data characterizes diverse recombinant viral nucleic acids

BACKGROUND

Genetic recombination is a tremendous source of intrahost diversity in viruses and is critical for their ability to rapidly adapt to new environments or fitness challenges. While viruses are routinely characterized using high-throughput sequencing techniques, characterizing the genetic products of recombination in next-generation sequencing data remains a challenge. Viral recombination events can be highly diverse and variable in nature, including simple duplications and deletions, or more complex events such as copy/snap-back recombination, intervirus or intersegment recombination, and insertions of host nucleic acids. Due to the variable mechanisms driving virus recombination and the different selection pressures acting on the progeny, recombination junctions rarely adhere to simple canonical sites or sequences. Furthermore, numerous different events may be present simultaneously in a viral population, yielding a complex mutational landscape.

FINDINGS

We have previously developed an algorithm called ViReMa (Virus Recombination Mapper) that bootstraps the bowtie short-read aligner to capture and annotate a wide range of recombinant species found within virus populations. Here, we have updated ViReMa to provide an "error density" function designed to accurately detect recombination events in the longer reads now routinely generated by the Illumina platforms and provide output reports for multiple types of recombinant species using standardized formats. We demonstrate the utility and flexibility of ViReMa in different settings to report deletion events in simulated data from Flock House virus, copy-back RNA species in Sendai viruses, short duplication events in HIV, and virus-to-host recombination in an archaeal DNA virus.

Unconventional secretion of unglycosylated ORF8 is critical for the cytokine storm during SARS-CoV-2 infection

Coronavirus disease 2019 is a respiratory infectious disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Evidence on the pathogenesis of SARS-CoV-2 is accumulating rapidly. In addition to structural proteins such as Spike and Envelope, the functional roles of non-structural and accessory proteins in regulating viral life cycle and host immune responses remain to be understood. Here, we show that open reading frame 8 (ORF8) acts as messenger for inter-cellular communication between alveolar epithelial cells and macrophages during SARS-CoV-2 infection. Mechanistically, ORF8 is a secretory protein that can be secreted by infected epithelial cells via both conventional and unconventional secretory pathways. Conventionally secreted ORF8 is glycosylated and loses the ability to recognize interleukin 17 receptor A of macrophages, possibly due to the steric hindrance imposed by N-glycosylation at Asn78. However, unconventionally secreted ORF8 does not undergo glycosylation without experiencing the ER-Golgi trafficking, thereby activating the downstream NF-κB signaling pathway and facilitating a burst of cytokine release. Furthermore, we show that ORF8 deletion in SARS-CoV-2 attenuates inflammation and yields less lung lesions in hamsters. Our data collectively highlights a role of ORF8 protein in the development of cytokine storms during SARS-CoV-2 infection.

Deletions across the SARS-CoV-2 Genome: Molecular Mechanisms and Putative Functional Consequences of Deletions in Accessory Genes

The analysis of deletions may reveal evolutionary trends and provide new insight into the surprising variability and rapidly spreading capability that SARS-CoV-2 has shown since its emergence. To understand the factors governing genomic stability, it is important to define the molecular mechanisms of deletions in the viral genome. In this work, we performed a statistical analysis of deletions. Specifically, we analyzed correlations between deletions in the SARS-CoV-2 genome and repetitive elements and documented a significant association of deletions with runs of identical (poly-) nucleotides and direct repeats. Our analyses of deletions in the accessory genes of SARS-CoV-2 suggested that there may be a hypervariability in ORF7A and ORF8 that is not associated with repetitive elements. Such recurrent search in a "sequence space" of accessory genes (that might be driven by natural selection) did not yet cause increased viability of the SARS-CoV-2 variants. However, deletions in the accessory genes may ultimately produce new variants that are more successful compared to the viral strains with the conventional architecture of the SARS-CoV-2 accessory genes.

The frequency of defective genomes in Omicron differs from that of the Alpha, Beta and Delta variants

The SARS-CoV-2 Omicron variant emerged showing higher transmissibility and possibly higher resistance to current COVID-19 vaccines than other variants dominating the global pandemic. In March 2020 we performed a study in clinical samples, where we found that a portion of genomes in the SARS-CoV-2 viral population accumulated deletions immediately before the S1/S2 cleavage site (furin-like cleavage site, PRRAR/S) of the spike gene, generating a frameshift and appearance of a premature stop codon. The main aim of this study was to determine the frequency of defective deletions in prevalent variants from the first to sixth pandemic waves in our setting and discuss whether the differences observed might support epidemiological proposals. The complete SARS-CoV-2 spike gene was deeply studied by next-generation sequencing using the MiSeq platform. More than 90 million reads were obtained from respiratory swab specimens of 78 COVID-19 patients with mild infection caused by the predominant variants circulating in the Barcelona city area during the six pandemic waves: B.1.5, B.1.1, B.1.177, Alpha, Beta, Delta, and Omicron. The frequency of defective genomes found in variants dominating the first and second waves was similar to that seen in Omicron, but differed from the frequencies seen in the Alpha, Beta and Delta variants. The changing pattern of mutations seen in the various SARS-CoV-2 variants driving the pandemic waves over time can affect viral transmission and immune escape. Here we discuss the putative biological effects of defective deletions naturally occurring before the S1/S2 cleavage site during adaption of the virus to human infection.

"Mutation blacklist" and "mutation whitelist" of SARS-CoV-2

Over the past two years, scientists throughout the world have completed more than 6 million SARS-CoV-2 genome sequences. Today, the number of SARS-CoV-2 genomes exceeds the total number of all other viral genomes. These genomes are a record of the evolution of SARS-CoV-2 in the human host, and provide information on the emergence of mutations. In this study, analysis of these sequenced genomes identified 296,728 de novo mutations (DNMs), and found that six types of base substitutions reached saturation in the sequenced genome population. Based on this analysis, a "mutation blacklist" of SARS-CoV-2 was compiled. The loci on the "mutation blacklist" are highly conserved, and these mutations likely have detrimental effects on virus survival, replication, and transmission. This information is valuable for SARS-CoV-2 research on gene function, vaccine design, and drug development. Through association analysis of DNMs and viral transmission rates, we identified 185 DNMs that positively correlated with the SARS-CoV-2 transmission rate, and these DNMs where classified as the "mutation whitelist" of SARS-CoV-2. The mutations on the "mutation whitelist" are beneficial for SARS-CoV-2 transmission and could therefore be used to evaluate the transmissibility of new variants. The occurrence of mutations and the evolution of viruses are dynamic processes. To more effectively monitor the mutations and variants of SARS-CoV-2, we built a SARS-CoV-2 mutation and variant monitoring and pre-warning system (MVMPS), which can monitor the occurrence and development of mutations and variants of SARS-CoV-2, as well as provide pre-warning for the prevention and control of SARS-CoV-2 (https://www.omicx.cn/). Additionally, this system could be used in real-time to update the "mutation whitelist" and "mutation blacklist" of SARS-CoV-2.

Emergence and spreading of the largest SARS-CoV-2 deletion in the Delta AY.20 lineage from Uruguay

The genetic variability of SARS-CoV-2 (genus Betacoronavirus, family Coronaviridae) has been scrutinized since its first detection in December 2019. Although the role of structural variants, particularly deletions, in virus evolution is little explored, these genome changes are extremely frequent. They are associated with relevant processes, including immune escape and attenuation. Deletions commonly occur in accessory ORFs and might even lead to the complete loss of one or more ORFs. This scenario poses an interesting question about the origin and spreading of extreme structural rearrangements that persist without compromising virus viability. Here, we analyze the genome of SARS-CoV-2 in late 2021 in Uruguay and identify a Delta lineage (AY.20) that experienced a large deletion (872 nucleotides according to the reference Wuhan strain) that removes the 7a, 7b, and 8 ORFs. Deleted viruses coexist with wild-type (without deletion) AY.20 and AY.43 strains. The Uruguayan deletion is like those identified in Delta strains from Poland and Japan but occurs in a different Delta clade. Besides providing proof of the circulation of this large deletion in America, we infer that the 872-deletion arises by the consecutive occurrence of a 6-nucleotide deletion, characteristic of delta strains, and an 866-nucleotide deletion that arose independently in the AY.20 Uruguayan lineage. The largest deletion occurs adjacent to transcription regulatory sequences needed to synthesize the nested set of subgenomic mRNAs that serve as templates for transcription. Our findings support the role of transcription sequences as a hotspot for copy-choice recombination and highlight the remarkable dynamic of SARS-CoV-2 genomes.

Extensive Recombination-driven Coronavirus Diversification Expands the Pool of Potential Pandemic Pathogens

The ongoing SARS-CoV-2 pandemic is the third zoonotic coronavirus identified in the last 20 years. Enzootic and epizootic coronaviruses of diverse lineages also pose a significant threat to livestock, as most recently observed for virulent strains of porcine epidemic diarrhea virus (PEDV) and swine acute diarrhea-associated coronavirus (SADS-CoV). Unique to RNA viruses, coronaviruses encode a proofreading exonuclease (ExoN) that lowers point mutation rates to increase the viability of large RNA virus genomes, which comes with the cost of limiting virus adaptation via point mutation. This limitation can be overcome by high rates of recombination that facilitate rapid increases in genetic diversification. To compare the dynamics of recombination between related sequences, we developed an open-source computational workflow (IDPlot) that bundles nucleotide identity, recombination, and phylogenetic analysis into a single pipeline. We analyzed recombination dynamics among three groups of coronaviruses with noteworthy impacts on human health and agriculture: SARSr-CoV, Betacoronavirus-1, and SADSr-CoV. We found that all three groups undergo recombination with highly diverged viruses from undersampled or unsampled lineages, including in typically highly conserved regions of the genome. In several cases, no parental origin of recombinant regions could be found in genetic databases, demonstrating our shallow characterization of coronavirus diversity and expanding the genetic pool that may contribute to future zoonotic events. Our results also illustrate the limitations of current sampling approaches for anticipating zoonotic threats to human and animal health.

SARS-CoV-2 disrupts host epigenetic regulation via histone mimicry

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged at the end of 2019 and caused the devastating global pandemic of coronavirus disease 2019 (COVID-19), in part because of its ability to effectively suppress host cell responses1-3. In rare cases, viral proteins dampen antiviral responses by mimicking critical regions of human histone proteins4-8, particularly those containing post-translational modifications required for transcriptional regulation9-11. Recent work has demonstrated that SARS-CoV-2 markedly disrupts host cell epigenetic regulation12-14. However, how SARS-CoV-2 controls the host cell epigenome and whether it uses histone mimicry to do so remain unclear. Here we show that the SARS-CoV-2 protein encoded by ORF8 (ORF8) functions as a histone mimic of the ARKS motifs in histone H3 to disrupt host cell epigenetic regulation. ORF8 is associated with chromatin, disrupts regulation of critical histone post-translational modifications and promotes chromatin compaction. Deletion of either the ORF8 gene or the histone mimic site attenuates the ability of SARS-CoV-2 to disrupt host cell chromatin, affects the transcriptional response to infection and attenuates viral genome copy number. These findings demonstrate a new function of ORF8 and a mechanism through which SARS-CoV-2 disrupts host cell epigenetic regulation. Further, this work provides a molecular basis for the finding that SARS-CoV-2 lacking ORF8 is associated with decreased severity of COVID-19.

SARS-COV-2/COVID-19: scenario, epidemiology, adaptive mutations, and environmental factors

The coronavirus pandemic of 2019 has already exerted an enormous impact. For over a year, the worldwide pandemic has ravaged the whole globe, with approximately 250 million verified human infection cases and a mortality rate surpassing 4 million. While the genetic makeup of the related pathogen (SARS-CoV-2) was identified, many unknown facets remain a mystery, comprising the virus's origin and evolutionary trend. There were many rumors that SARS-CoV-2 was human-borne and its evolution was predicted many years ago, but scientific investigation proved them wrong and concluded that bats might be the origin of SARS-CoV-2 and pangolins act as intermediary species to transmit the virus from bats to humans. Airborne droplets were found to be the leading cause of human-to-human transmission of this virus, but later studies showed that contaminated surfaces and other environmental factors are also involved in its transmission. The evolution of different SARS-CoV-2 variants worsens the condition and has become a challenge to overcome this pandemic. The emergence of COVID-19 is still a mystery, and scientists are unable to explain the exact origin of SARS-CoV-2. This review sheds light on the possible origin of SARS-CoV-2, its transmission, and the key factors that worsen the situation.

Genome Similarities between Human-Derived and Mink-Derived SARS-CoV-2 Make Mink a Potential Reservoir of the Virus

SARS-CoV-2 has RNA as the genome, which makes the virus more prone to mutations. Occasionally, mutations help a virus to cross the species barrier. SARS-CoV-2 infections in humans and minks (Neovison vison) are examples of zoonotic spillover. Many studies on the mutational analysis of human-derived SARS-CoV-2 have been published, but insight into the mink-derived SARS-CoV-2 genome of mutations is still required. Here, we performed a mutation analysis of the mink-derived SARS-CoV-2 genome sequences. We analyzed all available full-length mink-derived SARS-CoV-2 genome sequences on GISAID (214 genome sequences from the Netherlands and 133 genome sequences from Denmark). We found a striking resemblance between human-derived and mink-derived SARS-CoV-2. Our study showed that mutation patterns in the SARS-CoV-2 genome samples from the Netherlands and Denmark were different. Out of the 201 mutations we found, only 13 mutations were shared by the Netherlands' and Denmark's mink-derived samples. We found that six mutations were prevalent in the mink-derived SARS-CoV-2 genomes, and these six mutations are also known to be prevalent in human-derived SARS-CoV-2 variants. Our study reveals that the G27948T mutation in SARS-CoV-2 leads to truncation of ORF8, which was also reported in human-derived SARS-CoV-2, thus indicating that the virus can replicate without the full-length ORF8. These resemblances between mink-derived and human-derived SARS-CoV-2 enable the virus to cross the species barrier and suggest mink a potential reservoir for the virus.

Molecular adaptations during viral epidemics

In 1977, the world witnessed both the eradication of smallpox and the beginning of the modern age of genomics. Over the following half-century, 7 epidemic viruses of international concern galvanized virologists across the globe and led to increasingly extensive virus genome sequencing. These sequencing efforts exerted over periods of rapid adaptation of viruses to new hosts, in particular, humans provide insight into the molecular mechanisms underpinning virus evolution. Investment in virus genome sequencing was dramatically increased by the unprecedented support for phylogenomic analyses during the COVID-19 pandemic. In this review, we attempt to piece together comprehensive molecular histories of the adaptation of variola virus, HIV-1 M, SARS, H1N1-SIV, MERS, Ebola, Zika, and SARS-CoV-2 to the human host. Disruption of genes involved in virus-host interaction in animal hosts, recombination including genome segment reassortment, and adaptive mutations leading to amino acid replacements in virus proteins involved in host receptor binding and membrane fusion are identified as the key factors in the evolution of epidemic viruses.

Methodology-Centered Review of Molecular Modeling, Simulation, and Prediction of SARS-CoV-2

Despite tremendous efforts in the past two years, our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), virus-host interactions, immune response, virulence, transmission, and evolution is still very limited. This limitation calls for further in-depth investigation. Computational studies have become an indispensable component in combating coronavirus disease 2019 (COVID-19) due to their low cost, their efficiency, and the fact that they are free from safety and ethical constraints. Additionally, the mechanism that governs the global evolution and transmission of SARS-CoV-2 cannot be revealed from individual experiments and was discovered by integrating genotyping of massive viral sequences, biophysical modeling of protein-protein interactions, deep mutational data, deep learning, and advanced mathematics. There exists a tsunami of literature on the molecular modeling, simulations, and predictions of SARS-CoV-2 and related developments of drugs, vaccines, antibodies, and diagnostics. To provide readers with a quick update about this literature, we present a comprehensive and systematic methodology-centered review. Aspects such as molecular biophysics, bioinformatics, cheminformatics, machine learning, and mathematics are discussed. This review will be beneficial to researchers who are looking for ways to contribute to SARS-CoV-2 studies and those who are interested in the status of the field.

Attenuated Viral Replication of Avian Infectious Bronchitis Virus with a Novel 82-Nucleotide Deletion in the 5a Gene Indicates a Critical Role for 5a in Virus-Host Interactions

We previously found that a deletion in γ-coronavirus Infectious bronchitis virus (IBV) accessory gene 5a is critical for decreased viral pathogenicity in chickens. Here, we systematically analyzed IBV virus infection: invasion, genome replication, subgenomic mRNA (sgmRNA) synthesis, protein synthesis, and virion release. The ability of the mutant IBV strain rYN-Δ5a to invade susceptible cells was not significantly different from that of parental rYN. However, compared with rYN, the level of sgmRNA synthesis and genome replication after cell entry by rYN-Δ5a was significantly lower in the early stage, resulting in a significantly lower level of nucleoprotein (N) synthesis and a consequent significantly lower number of offspring viruses released into the supernatant. The detected 5a protein was diffusely distributed in the cytoplasm and perinuclear area. We identified 16 differentially expressed host proteins, 8 of which were found to be host nuclear and cytoplasmic transport-related proteins. Coimmunoprecipitation revealed an interaction between hemagglutinin (HA)-tagged TNPO1, TNPO3, XPO1, XPOT, RanBP1, and EIF2B4 proteins and Flag-tagged 5a protein, and laser confocal microscopy confirmed 5a protein colocalization with these proteins, indicating that 5a protein can cause changes in the host protein localization. These host proteins promote the nuclear localization of N proteins, so we believe that 5a protein can hijack host nucleoplasmic transport-related proteins to help N enter the nucleus. This may involve regulating the cell cycle to promote the optimal intracellular conditions for virus assembly or by participating in the regulation of nucleolar function as a strategy to optimize virus replication. IMPORTANCE Coronaviruses (CoVs) have a huge impact on humans and animals. It is important for the prevention and control of the viruses to assess the molecular mechanisms related to virulence attenuation. Here, we systematically analyzed a single cycle of virus infection by γ-CoV IBV lacking accessory protein 5a. We observed that a 5a deletion in the IBV genome affected virus replication and sgmRNA synthesis early in the virus life cycle, leading to decreases in protein synthesis, offspring virus assembly, and virion release in chicken embryonic kidney cells. IBV 5a protein was found to interact with multiple host nuclear and cytoplasmic transport- and translation-related proteins, which can also interact with IBV N and relocate it into the cell nucleus. These findings provide a comprehensive view regarding the importance of IBV accessory protein 5a and an important theoretical basis for studying the interaction between coronavirus and host cell proteins.

SARS-CoV-2 Accessory Protein ORF8 Decreases Antibody-Dependent Cellular Cytotoxicity

Viruses use many different strategies to evade host immune responses. In the case of SARS-CoV-2, its Spike mutates rapidly to escape from neutralizing antibodies. In addition to this strategy, ORF8, a small accessory protein encoded by SARS-CoV-2, helps immune evasion by reducing the susceptibility of SARS-CoV-2-infected cells to the cytotoxic CD8+ T cell response. Interestingly, among all accessory proteins, ORF8 is rapidly evolving and a deletion in this protein has been linked to milder disease. Here, we studied the effect of ORF8 on peripheral blood mononuclear cells (PBMC). Specifically, we found that ORF8 can bind monocytes as well as NK cells. Strikingly, ORF8 binds CD16a (FcγRIIIA) with nanomolar affinity and decreases the overall level of CD16 at the surface of monocytes and, to a lesser extent, NK cells. This decrease significantly reduces the capacity of PBMCs and particularly monocytes to mediate antibody-dependent cellular cytotoxicity (ADCC). Overall, our data identifies a new immune-evasion activity used by SARS-CoV-2 to escape humoral responses.

Genomic determinants of Furin cleavage in diverse European SARS-related bat coronaviruses

The furin cleavage site (FCS) in SARS-CoV-2 is unique within the Severe acute respiratory syndrome–related coronavirus (SrC) species. We re-assessed diverse SrC from European horseshoe bats and analyzed the spike-encoding genomic region harboring the FCS in SARS-CoV-2. We reveal molecular features in SrC such as purine richness and RNA secondary structures that resemble those required for FCS acquisition in avian influenza viruses. We discuss the potential acquisition of FCS through molecular mechanisms such as nucleotide substitution, insertion, or recombination, and show that a single nucleotide exchange in two European bat-associated SrC may suffice to enable furin cleavage. Furthermore, we show that FCS occurrence is variable in bat- and rodent-borne counterparts of human coronaviruses. Our results suggest that furin cleavage sites can be acquired in SrC via conserved molecular mechanisms known in other reservoir-bound RNA viruses and thus support a natural origin of SARS-CoV-2.

Genomic analyses of spike glycoprotein genes of European bat SARS-related coronaviruses suggest that furin cleavage sites can be acquired in the bat reservoir via conserved molecular mechanisms, supporting a natural origin of SARS-CoV-2.

Deep phylogenetic-based clustering analysis uncovers new and shared mutations in SARS-CoV-2 variants as a result of directional and convergent evolution

Nearly two decades after the last epidemic caused by a severe acute respiratory syndrome coronavirus (SARS-CoV), newly emerged SARS-CoV-2 quickly spread in 2020 and precipitated an ongoing global public health crisis. Both the continuous accumulation of point mutations, owed to the naturally imposed genomic plasticity of SARS-CoV-2 evolutionary processes, as well as viral spread over time, allow this RNA virus to gain new genetic identities, spawn novel variants and enhance its potential for immune evasion. Here, through an in-depth phylogenetic clustering analysis of upwards of 200,000 whole-genome sequences, we reveal the presence of previously unreported and hitherto unidentified mutations and recombination breakpoints in Variants of Concern (VOC) and Variants of Interest (VOI) from Brazil, India (Beta, Eta and Kappa) and the USA (Beta, Eta and Lambda). Additionally, we identify sites with shared mutations under directional evolution in the SARS-CoV-2 Spike-encoding protein of VOC and VOI, tracing a heretofore-undescribed correlation with viral spread in South America, India and the USA. Our evidence-based analysis provides well-supported evidence of similar pathways of evolution for such mutations in all SARS-CoV-2 variants and sub-lineages. This raises two pivotal points: (i) the co-circulation of variants and sub-lineages in close evolutionary environments, which sheds light onto their trajectories into convergent and directional evolution, and (ii) a linear perspective into the prospective vaccine efficacy against different SARS-CoV-2 strains.

ORF8 protein of SARS‐CoV‐2 reduces male fertility in mice

As one of the most rapidly evolving proteins of the genus Betacoronavirus, open reading frames (ORF8's) function and potential pathological consequence in vivo are still obscure. In this study, we show that the secretion of ORF8 is dependent on its N‐terminal signal peptide sequence and can be inhibited by reactive oxygen species scavenger and endoplasmic reticulum‐Golgi transportation inhibitor in cultured cells. To trace the effect of its possible in vivo secretion, we examined the plasma samples of coronavirus disease 2019 (COVID‐19) convalescent patients and found that the patients aged from 40 to 60 had higher antibody titers than those under 40. To explore ORF8's in vivo function, we administered the mice with ORF8 via tail‐vein injection to simulate the circulating ORF8 in the patient. Although no apparent difference in body weight, food intake, and vitality was detected between vehicle‐ and ORF8‐treated mice, the latter displayed morphological abnormalities of testes and epididymides, as indicated by the loss of the central ductal lumen accompanied by a decreased fertility in 5‐week‐old male mice. Furthermore, the analysis of gene expression in the testes between vehicle‐ and ORF8‐treated mice identified a decreased expression of Col1a1, the loss of which is known to be associated with mice's infertility. Although whether our observation in mice could be translated to humans remains unclear, our study provides a potential mouse model that can be used to investigate the impact of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) infection on the human reproductive system.

A simple reverse genetics method to generate recombinant coronaviruses

Engineering recombinant viruses is a pre‐eminent tool for deciphering the biology of emerging viral pathogens such as the severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2). However, the large size of coronavirus genomes renders the current reverse genetics methods challenging. Here, we describe a simple method based on “infectious subgenomic amplicons” (ISA) technology to generate recombinant infectious coronaviruses with no need for reconstruction of the complete genomic cDNA and apply this method to SARS‐CoV‐2 and also to the feline enteric coronavirus. In both cases we rescue wild‐type viruses with biological characteristics similar to original strains. Specific mutations and fluorescent red reporter genes can be readily incorporated into the SARS‐CoV‐2 genome enabling the generation of a genomic variants and fluorescent reporter strains for in vivo experiments, serological diagnosis, and antiviral assays. The swiftness and simplicity of the ISA method has the potential to facilitate the advance of coronavirus reverse genetics studies, to explore the molecular biological properties of the SARS‐CoV‐2 variants, and to accelerate the development of effective therapeutic reagents.

Probing the Immune System Dynamics of the COVID-19 Disease for Vaccine Designing and Drug Repurposing Using Bioinformatics Tools

The pathogenesis of COVID-19 is complicated by immune dysfunction. The impact of immune-based therapy in COVID-19 patients has been well documented, with some notable studies on the use of anti-cytokine medicines. However, the complexity of disease phenotypes, patient heterogeneity and the varying quality of evidence from immunotherapy studies provide problems in clinical decision-making. This review seeks to aid therapeutic decision-making by giving an overview of the immunological responses against COVID-19 disease that may contribute to the severity of the disease. We have extensively discussed theranostic methods for COVID-19 detection. With advancements in technology, bioinformatics has taken studies to a higher level. The paper also discusses the application of bioinformatics and machine learning tools for the diagnosis, vaccine design and drug repurposing against SARS-CoV-2.

Viroporins: Structure, function, and their role in the life cycle of SARS-CoV-2

Viroporins are indispensable for viral replication. As intracellular ion channels they disturb pH gradients of organelles and allow Ca²⁺ flux across ER membranes. Viroporins interact with numerous intracellular proteins and pathways and can trigger inflammatory responses. Thus, they are relevant targets in the search for antiviral drugs. Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) underlies the world-wide pandemic of COVID-19, where an effective therapy is still lacking despite impressive progress in the development of vaccines and vaccination campaigns. Among the 29 proteins of SARS-CoV-2, the E- and ORF3a proteins have been identified as viroporins that contribute to the massive release of inflammatory cytokines observed in COVID-19. Here, we describe structure and function of viroporins and their role in inflammasome activation and cellular processes during the virus replication cycle. Techniques to study viroporin function are presented, with a focus on cellular and electrophysiological assays. Contributions of SARS-CoV-2 viroporins to the viral life cycle are discussed with respect to their structure, channel function, binding partners, and their role in viral infection and virus replication. Viroporin sequences of new variants of concern (α–ο) of SARS-CoV-2 are briefly reviewed as they harbour changes in E and 3a proteins that may affect their function.

Next-Generation Sequencing for Confronting Virus Pandemics

Virus pandemics have happened, are happening and will happen again. In recent decades, the rate of zoonotic viral spillover into humans has accelerated, mirroring the expansion of our global footprint and travel network, including the expansion of viral vectors and the destruction of natural spaces, bringing humans closer to wild animals. Once viral cross-species transmission to humans occurs, transmission cannot be stopped by cement walls but by developing barriers based on knowledge that can prevent or reduce the effects of any pandemic. Controlling a local transmission affecting few individuals is more efficient that confronting a community outbreak in which infections cannot be traced. Genetic detection, identification, and characterization of infectious agents using next-generation sequencing (NGS) has been proven to be a powerful tool allowing for the development of fast PCR-based molecular assays, the rapid development of vaccines based on mRNA and DNA, the identification of outbreaks, transmission dynamics and spill-over events, the detection of new variants and treatment of vaccine resistance mutations, the development of direct-acting antiviral drugs, the discovery of relevant minority variants to improve knowledge of the viral life cycle, strengths and weaknesses, the potential for becoming dominant to take appropriate preventive measures, and the discovery of new routes of viral transmission.

SARS-CoV-2 accessory protein ORF8 is secreted extracellularly as a glycoprotein homodimer

ORF8 is an accessory protein encoded by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Consensus regarding the biological functions of ORF8 is lacking, largely because the fundamental characteristics of this protein in cells have not been determined. To clarify these features, we herein established an ORF8 expression system in 293T cells. Using this system, approximately 41% of the ORF8 expressed in 293T cells were secreted extracellularly as a glycoprotein homodimer with inter/intramolecular disulfide bonds. Intracellular ORF8 was sensitive to the glycosidase Endo H, whereas the secreted portion was Endo-H-resistant, suggesting that secretion occurs via a conventional pathway. Additionally, immunoblotting analysis showed that the total amounts of the major histocompatibility complex class Ι (MHC-I), angiotensin-converting enzyme 2 (ACE2), and SARS-CoV-2 spike (CoV-2 S) proteins coexpressed in cells were not changed by the increased ORF8 expression, although FACS analysis revealed that the expression of the cell surface MHC-I protein, but not that of ACE2 and CoV-2 S proteins, was reduced by ORF8 expression. Finally, we demonstrate by RNA-seq analysis that ORF8 had no significant stimulatory effects in human primary monocyte-derived macrophages (MDMs). Taken together, our results provide fundamental evidence that the ORF8 glycoprotein acts as a secreted homodimer, and its functions are likely associated with the intracellular transport and/or extracellular signaling in SARS-CoV-2 infection.

Robust Virus-Specific Adaptive Immunity in COVID-19 Patients with SARS-CoV-2 Δ382 Variant Infection

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) that have become dominant as the pandemic progresses bear the ORF8 mutation together with multiple spike mutations. A 382-nucleotide deletion (Δ382) in the ORF7b and ORF8 regions has been associated with milder disease phenotype and less systemic inflammation in COVID-19 patients. However, its impact on host immunity against SARS-CoV-2 remains undefined. Here, RNA-sequencing was performed to elucidate whole blood transcriptomic profiles and identify contrasting immune signatures between patients infected with either wildtype or Δ382 SARS-CoV-2 variant. Interestingly, the immune landscape of Δ382 SARS-CoV-2 infected patients featured an increased adaptive immune response, evidenced by enrichment of genes related to T cell functionality, a more robust SARS-CoV-2-specific T cell immunity, as well as a more rapid antibody response. At the molecular level, eukaryotic initiation factor 2 signaling was found to be upregulated in patients bearing Δ382, and its associated genes were correlated with systemic levels of T cell-associated and pro-inflammatory cytokines. This study provides more in-depth insight into the host-pathogen interactions of ORF8 with great promise as a therapeutic target to combat SARS-CoV-2 infection.

Metaviromic identification of discriminative genomic features in SARS-CoV-2 using machine learning

The COVID-19 pandemic caused by SARS-CoV-2 has become a major threat across the globe. Here, we developed machine learning approaches to identify key pathogenic regions in coronavirus genomes. We trained and evaluated 7,562,625 models on 3,665 genomes including SARS-CoV-2, MERS-CoV, SARS-CoV, and other coronaviruses of human and animal origins to return quantitative and biologically interpretable signatures at nucleotide and amino acid resolutions. We identified hotspots across the SARS-CoV-2 genome, including previously unappreciated features in spike, RdRp, and other proteins. Finally, we integrated pathogenicity genomic profiles with B cell and T cell epitope predictions for enrichment of sequence targets to help guide vaccine development. These results provide a systematic map of predicted pathogenicity in SARS-CoV-2 that incorporates sequence, structural, and immunologic features, providing an unbiased collection of genetic elements for functional studies. This metavirome-based framework can also be applied for rapid characterization of new coronavirus strains or emerging pathogenic viruses.

•

Machine learning identifies discriminative signatures in coronavirus genomes

•

Hotspots in key viral proteins have evolutionary and structural significance

•

Integration of hotspots with B cell and T cell epitopes identify joint features

•

Hotspots correlate with emerging variants of concern for mutation prioritization

Identifying which genomic regions of the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus are pathogenic remains a major challenge in COVID-19 research. However, there is currently a lack of systematic and unbiased methods for such functional characterization. In this study, we set up a machine learning-based approach to identify which genomic regions distinguish SARS-CoV-2 and other high case fatality rate coronaviruses from other coronaviruses. Discriminative scores were obtained for every nucleotide in the SARS-CoV-2 genome. We then performed a series of evolutionary and structural analyses of candidate hotspots, as well as integrative analyses with predicted B cell and T cell epitopes and emerging variants of concern. Our approach can be extended to other viral genomes or microbial pathogens to gain insights on which sequence features are pathogenic or immunogenic.

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.

Review: Development of SARS-CoV-2 immuno-enhanced COVID-19 vaccines with nano-platform

Vaccination is the most effective way to prevent coronavirus disease 2019 (COVID-19). Vaccine development approaches consist of viral vector vaccines, DNA vaccine, RNA vaccine, live attenuated virus, and recombinant proteins, which elicit a specific immune response. The use of nanoparticles displaying antigen is one of the alternative approaches to conventional vaccines. This is due to the fact that nano-based vaccines are stable, able to target, form images, and offer an opportunity to enhance the immune responses. The diameters of ultrafine nanoparticles are in the range of 1-100 nm. The application of nanotechnology on vaccine design provides precise fabrication of nanomaterials with desirable properties and ability to eliminate undesirable features. To be successful, nanomaterials must be uptaken into the cell, especially into the target and able to modulate cellular functions at the subcellular levels. The advantages of nano-based vaccines are the ability to protect a cargo such as RNA, DNA, protein, or synthesis substance and have enhanced stability in a broad range of pH, ambient temperatures, and humidity for long-term storage. Moreover, nano-based vaccines can be engineered to overcome biological barriers such as nonspecific distribution in order to elicit functions in antigen presenting cells. In this review, we will summarize on the developing COVID-19 vaccine strategies and how the nanotechnology can enhance antigen presentation and strong immunogenicity using advanced technology in nanocarrier to deliver antigens. The discussion about their safe, effective, and affordable vaccines to immunize against COVID-19 will be highlighted.

Possible Role of Accessory Proteins in the Viral Replication for the 20I/501Y.V1 (B.1.1.7) SARS CoV-2 Variant

The emergence of new severe acute respiratory syndrome coronavirus-2 (SARS CoV-2) has been a global concern. The B.1.1.7 variant of SARS CoV-2 is reported to cause higher transmission. The study investigates the replication cycle and transcriptional pattern of the B.1.1.7 to hypothesis the possible role of different genes in viral replication. It was observed that the B.1.1.7 variant required a longer maturation time. The transcriptional response demonstrated higher expression of ORF6 and ORF8 compared to nucleocapsid transcript till the eclipse period which might influence higher viral replication. The number of infectious viruses titer is higher in the B.1.1.7, despite a lesser copy number than B.1, indicating higher transmissibility. The experimental evidence published linked ORF6 and ORF8 to play important role in replication and we also observed their higher expression. This leads us to hypothesis the possible role of ORF6 and ORF8 in B.1.1.7 higher replication which causes higher transmission.

Host-dependent editing of SARS-CoV-2 in COVID-19 patients

A common trait among RNA viruses is their high capability to acquire genetic variability due to viral and host mechanisms. Next-generation sequencing (NGS) analysis enables the deep study of the viral quasispecies in samples from infected individuals. In this study, the viral quasispecies complexity and single nucleotide polymorphisms of the SARS-CoV-2 spike gene of coronavirus disease 2019 (COVID-19) patients with mild or severe disease were investigated using next-generation sequencing (Illumina platform). SARS-CoV-2 spike variability was higher in patients with long-lasting infection. Most substitutions found were present at frequencies lower than 1%, and had an A → G or T → C pattern, consistent with variants caused by adenosine deaminase acting on RNA-1 (ADAR1). ADAR1 affected a small fraction of replicating genomes, but produced multiple, mainly non-synonymous mutations. ADAR1 editing during replication rather than the RNA-dependent RNA polymerase (nsp12) was the predominant mechanism generating SARS-CoV-2 genetic variability. However, the mutations produced are not fixed in the infected human population, suggesting that ADAR1 may have an antiviral role, whereas nsp12-induced mutations occurring in patients with high viremia and persistent infection are the main source of new SARS-CoV-2 variants.

Characterization of SARS-CoV-2 ORF6 deletion variants detected in a nosocomial cluster during routine genomic surveillance, Lyon, France

During routine molecular surveillance of SARS-CoV-2 performed at the National Reference Center of Respiratory Viruses (Lyon, France) (n = 229 sequences collected February-April 2020), two frameshifting deletions were detected in the open reading frame 6, at the same position (27267). While a 26-nucleotide deletion variant (D26) was only found in one nasopharyngeal sample in March 2020, the 34-nucleotide deletion (D34) was found within a single geriatric hospital unit in 5/9 patients and one health care worker in April 2020. Phylogeny analysis strongly suggested a nosocomial transmission of D34, with potential fecal transmission, as also identified in a stool sample. No difference in disease severity was observed between patients hospitalized in the geriatric unit infected with WT or D34. In vitro D26 and D34 characterization revealed comparable replication kinetics with the wild-type (WT), but differential host immune responses. While interferon-stimulated genes were similarly upregulated after infection with WT and ORF6 variants, the latter specifically induced overexpression of 9 genes coding for inflammatory cytokines in the NF-kB pathway, including CCL2/MCP1, PTX3, and TNFα, for which high plasma levels have been associated with severe COVID-19. Our findings emphasize the need to monitor the occurrence of ORF6 deletions and assess their impact on the host immune response.

Expansion of a SARS-CoV-2 Delta variant with an 872 nt deletion encompassing ORF7a, ORF7b and ORF8, Poland, July to August 2021

Routine genomic surveillance on samples from COVID-19 patients collected in Poland during summer 2021 revealed the emergence of a SARS-CoV-2 Delta variant with a large 872 nt deletion. This change, confirmed by Sanger and deep sequencing, causes complete loss of ORF7a, ORF7b, and ORF8 genes. The index case carrying the deletion is unknown. The standard pipeline for sequencing may mask this deletion with a long stretch of N’s. Effects of this deletion on phenotype or immune evasion needs further study.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for next-generation sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called ‘Tiled-ClickSeq’, which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5’UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

Structural Analysis of SARS-CoV-2 ORF8 Protein: Pathogenic and Therapeutic Implications

Current therapeutic strategies and vaccines against SARS-CoV-2 are mainly focused on the Spike protein despite there are other viral proteins with important roles in COVID-19 pathogenicity. For example, ORF8 restructures vesicular trafficking in the host cell, impacts intracellular immunity through the IFN-I signaling, and growth pathways through the mitogen-activated protein kinases (MAPKs). In this mini-review, we analyze the main structural similarities of ORF8 with immunological molecules such as IL-1, contributing to the immunological deregulation observed in COVID-19. We also propose that the blockage of some effector functions of ORF8 with Rapamycin, such as the mTORC1 activation through MAPKs 40 pathway, with Rapamycin, can be a promising approach to reduce COVID-19 mortality.

Cloning of a Passage-Free SARS-CoV-2 Genome and Mutagenesis Using Red Recombination

The ongoing pandemic coronavirus (CoV) disease 2019 (COVID-19) by severe acute respiratory syndrome CoV-2 (SARS-CoV-2) has already caused substantial morbidity, mortality, and economic devastation. Reverse genetic approaches to generate recombinant viruses are a powerful tool to characterize and understand newly emerging viruses. To contribute to the global efforts for countermeasures to control the spread of SARS-CoV-2, we developed a passage-free SARS-CoV-2 clone based on a bacterial artificial chromosome (BAC). Moreover, using a Lambda-based Red recombination, we successfully generated different reporter and marker viruses, which replicated similar to a clinical isolate in a cell culture. Moreover, we designed a full-length reporter virus encoding an additional artificial open reading frame with wild-type-like replication features. The virus-encoded reporters were successfully applied to ease antiviral testing in cell culture models. Furthermore, we designed a new marker virus encoding 3xFLAG-tagged nucleocapsid that allows the detection of incoming viral particles and, in combination with bio-orthogonal labeling for the visualization of viral RNA synthesis via click chemistry, the spatiotemporal tracking of viral replication on the single-cell level. In summary, by applying BAC-based Red recombination, we developed a powerful, reliable, and convenient platform that will facilitate studies answering numerous questions concerning the biology of SARS-CoV-2.

Multiple Occurrences of a 168-Nucleotide Deletion in SARS-CoV-2 ORF8, Unnoticed by Standard Amplicon Sequencing and Variant Calling Pipelines

Genomic surveillance of the SARS-CoV-2 pandemic is crucial and mainly achieved by amplicon sequencing protocols. Overlapping tiled-amplicons are generated to establish contiguous SARS-CoV-2 genome sequences, which enable the precise resolution of infection chains and outbreaks. We investigated a SARS-CoV-2 outbreak in a local hospital and used nanopore sequencing with a modified ARTIC protocol employing 1200 bp long amplicons. We detected a long deletion of 168 nucleotides in the ORF8 gene in 76 samples from the hospital outbreak. This deletion is difficult to identify with the classical amplicon sequencing procedures since it removes two amplicon primer-binding sites. We analyzed public SARS-CoV-2 sequences and sequencing read data from ENA and identified the same deletion in over 100 genomes belonging to different lineages of SARS-CoV-2, pointing to a mutation hotspot or to positive selection. In almost all cases, the deletion was not represented in the virus genome sequence after consensus building. Additionally, further database searches point to other deletions in the ORF8 coding region that have never been reported by the standard data analysis pipelines. These findings and the fact that ORF8 is especially prone to deletions, make a clear case for the urgent necessity of public availability of the raw data for this and other large deletions that might change the physiology of the virus towards endemism.

Decomposing the sources of SARS-CoV-2 fitness variation in the United States

The fitness of a pathogen is a composite phenotype determined by many different factors influencing growth rates both within and between hosts. Determining what factors shape fitness at the host population-level is especially challenging because both intrinsic factors like pathogen genetics and extrinsic factors such as host behavior influence between-host transmission potential. This challenge has been highlighted by controversy surrounding the population-level fitness effects of mutations in the SARS-CoV-2 genome and their relative importance when compared against non-genetic factors shaping transmission dynamics. Building upon phylodynamic birth-death models, we develop a new framework to learn how hundreds of genetic and non-genetic factors have shaped the fitness of SARS-CoV-2. We estimate the fitness effects of all amino acid variants and several structural variants that have circulated in the United States between February 2020 and March 2021 from viral phylogenies. We also estimate how much fitness variation among pathogen lineages is attributable to genetic versus non-genetic factors such as spatial heterogeneity in transmission rates. Before September 2020, most fitness variation between lineages can be explained by background spatial heterogeneity in transmission rates across geographic regions. Starting in late 2020, genetic variation in fitness increased dramatically with the emergence of several new lineages including B.1.1.7, B.1.427, B.1.429 and B.1.526. Our analysis also indicates that genetic variants in less well-explored genomic regions outside of Spike may be contributing significantly to overall fitness variation in the viral population.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for Next-Generation Sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called 'Tiled-ClickSeq', which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5'UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

Molecular Level Dissection of Critical Spike Mutations in SARS-CoV-2 Variants of Concern (VOCs): A Simplified Review

SARS-CoV-2 virus during its spread in the last one and half year has picked up critical changes in its genetic code i.e. mutations, which have leads to deleterious epidemiological characteristics. Due to critical role of spike protein in cell entry and pathogenesis, mutations in spike regions have been reported to enhance transmissibility, disease severity, possible escape from vaccine-induced immune response and reduced diagnostic sensitivity/specificity. Considering the structure-function impact of mutations, understanding the molecular details of these key mutations of newly emerged variants/lineages is of urgent concern. In this review, we have explored the literature on key spike mutations harbored by alpha, beta, gamma and delta 'variants of concern' (VOCs) and discussed their molecular consequences in the context of resultant virus biology. Commonly all these VOCs i.e. B.1.1.7, B.1.351, P.1 and B.1.617.2 lineages have decisive mutation in Receptor Binding Motif (RBM) region and/or region around Furin cleavage site (FCS) of spike protein. In general, mutation induced disruption of intra-molecular interaction enhances molecular flexibility leading to exposure of spike protein surface in these lineages to make it accessible for inter-molecular interaction with hACE2. A disruption of spike antigen-antibody inter-molecular interactions in epitope region due to the chemical nature of substituting amino acid hampers the neutralization efficacy. Simplified surveillance of mutation induced changes and their consequences at molecular level can contribute in rationalizing mutation's impact on virus biology. It is believed that molecular level dissection of these key spike mutation will assist the future investigations for a more resilient outcome against severity of COVID-19.