Evolutionary dynamics of the SARS-CoV-2 ORF8 accessory gene

Structural and functional effects of the L84S mutant in the SARS-COV-2 ORF8 dimer based on microsecond molecular dynamics study

The L84S mutation has been observed frequently in the ORF8 protein of SARS-CoV-2, which is an accessory protein involved in various important functions such as virus propagation, pathogenesis, and evading the immune response. However, the specific effects of this mutation on the dimeric structure of ORF8 and its impacts on interactions with host components and immune responses are not well understood. In this study, we performed one microsecond molecular dynamics (MD) simulation and analyzed the dimeric behavior of the L84S and L84A mutants in comparison to the native protein. The MD simulations revealed that both mutations caused changes in the conformation of the ORF8 dimer, influenced protein folding mechanisms, and affected the overall structural stability. In particular, the 73YIDI76 motif has found to be significantly affected by the L84S mutation, leading to structural flexibility in the region connecting the C-terminal β4 and β5 strands. This flexibility might be responsible for virus immune modulation.  The free energy landscape (FEL) and principle component analysis (PCA) have also supported our investigation. Overall, the L84S and L84A mutations affect the ORF8 dimeric interfaces by reducing the frequency of protein-protein interacting residues (Arg52, Lys53, Arg98, Ile104, Arg115, Val117, Asp119, Phe120, and Ile121) in the ORF8 dimer.  Our findings provide detail insights for further research in designing structure-based therapeutics against the SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Abrogation of ORF8-IRF3 binding interface with Carbon nanotube derivatives to rescue the host immune system against SARS-CoV-2 by using molecular screening and simulation approaches

The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has led to over six million deaths worldwide. In human immune system, the type 1 interferon (IFN) pathway plays a crucial role in fighting viral infections. However, the ORF8 protein of the virus evade the immune system by interacting with IRF3, hindering its nuclear translocation and consequently downregulate the type I IFN signaling pathway. To block the binding of ORF8-IRF3 and inhibit viral pathogenesis a quick discovery of an inhibitor molecule is needed. Therefore, in the present study, the interface between the ORF8 and IRF3 was targeted on a high-affinity carbon nanotube by using computational tools. After analysis of 62 carbon nanotubes by multiple docking with the induced fit model, the top five compounds with high docking scores of - 7.94 kcal/mol, - 7.92 kcal/mol, - 7.28 kcal/mol, - 7.19 kcal/mol and - 7.09 kcal/mol (top hit1-5) were found to have inhibitory activity against the ORF8-IRF3 complex. Molecular dynamics analysis of the complexes revealed the high compactness of residues, stable binding, and strong hydrogen binding network among the ORF8-nanotubes complexes. Moreover, the total binding free energy for top hit1-5 was calculated to be - 43.21 ± 0.90 kcal/mol, - 41.17 ± 0.99 kcal/mol, - 48.85 ± 0.62 kcal/mol, - 43.49 ± 0.77 kcal/mol, and - 31.18 ± 0.78 kcal/mol respectively. These results strongly suggest that the identified top five nanotubes (hit1-5) possess significant potential for advancing and exploring innovative drug therapies. This underscores their suitability for subsequent in vivo and in vitro experiments, marking them as promising candidates worthy of further investigation.

Surfing the Waves of SARS-CoV-2: Analysis of Viral Genome Variants Using an NGS Survey in Verona, Italy

The availability of new technologies for deep sequencing, including next-generation sequencing (NGS), allows for the detection of viral genome variations. The epidemiological determination of SARS-CoV-2 viral genome changes during the pandemic waves displayed the genome evolution and subsequent onset of variants over time. These variants were often associated with a different impact on viral transmission and disease severity. We investigated, in a retrospective study, the trend of SARS-CoV-2-positive samples collected from the start of the Italian pandemic (January 2020) to June 2023. In addition, viral RNAs extracted from 938 nasopharyngeal swab samples were analyzed using NGS between February 2022 and June 2023. Sequences were analyzed with bioinformatic tools to identify lineages and mutations and for phylogenetic studies. Six pandemic waves were detected. In our samples, we predominantly detected BA.2, BQ.1, BA.5.1, BA.5.2, and, more recently, XBB.1 and its subvariants. The data describe the SARS-CoV-2 genome evolution involved in viral interactions with the host and the dynamics of specific genome mutations and deletions.

Prevalence of persistent SARS-CoV-2 in a large community surveillance study

Persistent SARS-CoV-2 infections may act as viral reservoirs that could seed future outbreaks1-5, give rise to highly divergent lineages6-8 and contribute to cases with post-acute COVID-19 sequelae (long COVID)9,10. However, the population prevalence of persistent infections, their viral load kinetics and evolutionary dynamics over the course of infections remain largely unknown. Here, using viral sequence data collected as part of a national infection survey, we identified 381 individuals with SARS-CoV-2 RNA at high titre persisting for at least 30 days, of which 54 had viral RNA persisting at least 60 days. We refer to these as 'persistent infections' as available evidence suggests that they represent ongoing viral replication, although the persistence of non-replicating RNA cannot be ruled out in all. Individuals with persistent infection had more than 50% higher odds of self-reporting long COVID than individuals with non-persistent infection. We estimate that 0.1-0.5% of infections may become persistent with typically rebounding high viral loads and last for at least 60 days. In some individuals, we identified many viral amino acid substitutions, indicating periods of strong positive selection, whereas others had no consensus change in the sequences for prolonged periods, consistent with weak selection. Substitutions included mutations that are lineage defining for SARS-CoV-2 variants, at target sites for monoclonal antibodies and/or are commonly found in immunocompromised people11-14. This work has profound implications for understanding and characterizing SARS-CoV-2 infection, epidemiology and evolution.

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.

SARS-CoV-2 ORF8 as a Modulator of Cytokine Induction: Evidence and Search for Molecular Mechanisms

Severe cases of SARS-CoV-2 infection are characterized by an immune response that leads to the overproduction of pro-inflammatory cytokines, resulting in lung damage, cardiovascular symptoms, hematologic symptoms, acute kidney injury and multiple organ failure that can lead to death. This remarkable increase in cytokines and other inflammatory molecules is primarily caused by viral proteins, and particular interest has been given to ORF8, a unique accessory protein specific to SARS-CoV-2. Despite plenty of research, the precise mechanisms by which ORF8 induces proinflammatory cytokines are not clear. Our investigations demonstrated that ORF8 augments production of IL-6 induced by Poly(I:C) in human embryonic kidney (HEK)-293 and monocyte-derived dendritic cells (mono-DCs). We discuss our findings and the multifaceted roles of ORF8 as a modulator of cytokine response, focusing on type I interferon and IL-6, a key component of the immune response to SARS-CoV-2. In addition, we explore the hypothesis that ORF8 may act through pattern recognition receptors of dsRNA such as TLRs.

Early monitoring-to-warning Internet of Things system for emerging infectious diseases via networking of light-triggered point-of-care testing devices

Early monitoring and warning arrangements are effective ways to distinguish infectious agents and control the spread of epidemic diseases. Current testing technologies, which cannot achieve rapid detection in the field, have a risk of slowing down the response time to the disease. In addition, there is still no epidemic surveillance system, implementing prevention and control measures is slow and inefficient. Motivated by these clinical needs, a sample-to-answer genetic diagnosis platform based on light-controlled capillary modified with a photocleavable linker is first developed, which could perform nucleic acid separation and release by light irradiation in less than 30 seconds. Then, on site polymerase chain reaction was performed in a handheld closed-loop convective system. Test reports are available within 20 min. Because this method is portable, rapid, and easy to operate, it has great potential for point-of-care testing. Additionally, through multiple device networking, a real-time artificial intelligence monitoring system for pathogens was developed on a cloud server. Through data reception, analysis, and visualization, the system can send early warning signals for disease control and prevention. Thus, anti-epidemic measures can be implemented effectively, and deploying and running this system can improve the capabilities for the prevention and control of infectious diseases.

Detection of SARS-CoV-2 Δ426 ORF8 Deletion Mutant Cluster in NGS Screening

Next-generation sequencing (NGS) from SARS-CoV-2-positive swabs collected during the last months of 2022 revealed a large deletion spanning ORF7b and ORF8 (426 nt) in six patients infected with the BA.5.1 Omicron variant. This extensive genome loss removed a large part of these two genes, maintaining in frame the first 22 aminoacids of ORF7b and the last three aminoacids of ORF8. Interestingly, the deleted region was flanked by two small repeats, which were likely involved in the formation of a hairpin structure. Similar rearrangements, comparable in size and location to the deletion, were also identified in 15 sequences in the NCBI database. In this group, seven out of 15 cases from the USA and Switzerland presented both the BA.5.1 variant and the same 426 nucleotides deletion. It is noteworthy that three out of six cases were detected in patients with immunodeficiency, and it is conceivable that this clinical condition could promote the replication and selection of these mutations.

A neonatal mouse model characterizes transmissibility of SARS-CoV-2 variants and reveals a role for ORF8

Small animal models have been a challenge for the study of SARS-CoV-2 transmission, with most investigators using golden hamsters or ferrets. Mice have the advantages of low cost, wide availability, less regulatory and husbandry challenges, and the existence of a versatile reagent and genetic toolbox. However, adult mice do not robustly transmit SARS-CoV-2. Here we establish a model based on neonatal mice that allows for transmission of clinical SARS-CoV-2 isolates. We characterize tropism, respiratory tract replication and transmission of ancestral WA-1 compared to variants Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2), Omicron BA.1 and Omicron BQ.1.1. We identify inter-variant differences in timing and magnitude of infectious particle shedding from index mice, both of which shape transmission to contact mice. Furthermore, we characterize two recombinant SARS-CoV-2 lacking either the ORF6 or ORF8 host antagonists. The removal of ORF8 shifts viral replication towards the lower respiratory tract, resulting in significantly delayed and reduced transmission in our model. Our results demonstrate the potential of our neonatal mouse model to characterize viral and host determinants of SARS-CoV-2 transmission, while revealing a role for an accessory protein in this context.

How SARS-CoV-2 Alters the Regulation of Gene Expression in Infected Cells

Nonstructural accessory proteins in viruses play a key role in hijacking the basic cellular mechanisms, which is essential to promote the virus survival and evasion of the immune system. The immonuglobulin-like open reading frame 8 (ORF8) protein expressed by SARS-CoV-2 accumulates in the nucleus and may influence the regulation of the gene expression in infected cells. In this contribution, by using microsecond time-scale all-atom molecular dynamics simulations, we unravel the structural bases behind the epigenetic action of ORF8. In particular, we highlight how the protein is able to form stable aggregates with DNA through a histone tail-like motif, and how this interaction is influenced by post-translational modifications, such as acetylation and methylation, which are known epigenetic markers in histones. Our work not only clarifies the molecular mechanisms behind the perturbation of the epigenetic regulation caused by the viral infection but also offers an unusual perspective which may foster the development of original antivirals.

SARS-CoV-2 ORF8: A Rapidly Evolving Immune and Viral Modulator in COVID-19

The COVID-19 pandemic has resulted in upwards of 6.8 million deaths over the past three years, and the frequent emergence of variants continues to strain global health. Although vaccines have greatly helped mitigate disease severity, SARS-CoV-2 is likely to remain endemic, making it critical to understand its viral mechanisms contributing to pathogenesis and discover new antiviral therapeutics. To efficiently infect, this virus uses a diverse set of strategies to evade host immunity, accounting for its high pathogenicity and rapid spread throughout the COVID-19 pandemic. Behind some of these critical host evasion strategies is the accessory protein Open Reading Frame 8 (ORF8), which has gained recognition in SARS-CoV-2 pathogenesis due to its hypervariability, secretory property, and unique structure. This review discusses the current knowledge on SARS-CoV-2 ORF8 and proposes actualized functional models describing its pivotal roles in both viral replication and immune evasion. A better understanding of ORF8's interactions with host and viral factors is expected to reveal essential pathogenic strategies utilized by SARS-CoV-2 and inspire the development of novel therapeutics to improve COVID-19 disease outcomes.

Correlated substitutions reveal SARS-like coronaviruses recombine frequently with a diverse set of structured gene pools

Quantifying SARS-like coronavirus (SL-CoV) evolution is critical to understanding the origins of SARS-CoV-2 and the molecular processes that could underlie future epidemic viruses. While genomic analyses suggest recombination was a factor in the emergence of SARS-CoV-2, few studies have quantified recombination rates among SL-CoVs. Here, we infer recombination rates of SL-CoVs from correlated substitutions in sequencing data using a coalescent model with recombination. Our computationally-efficient, non-phylogenetic method infers recombination parameters of both sampled sequences and the unsampled gene pools with which they recombine. We apply this approach to infer recombination parameters for a range of positive-sense RNA viruses. We then analyze a set of 191 SL-CoV sequences (including SARS-CoV-2) and find that ORF1ab and S genes frequently undergo recombination. We identify which SL-CoV sequence clusters have recombined with shared gene pools, and show that these pools have distinct structures and high recombination rates, with multiple recombination events occurring per synonymous substitution. We find that individual genes have recombined with different viral reservoirs. By decoupling contributions from mutation and recombination, we recover the phylogeny of non-recombined portions for many of these SL-CoVs, including the position of SARS-CoV-2 in this clonal phylogeny. Lastly, by analyzing >400,000 SARS-CoV-2 whole genome sequences, we show current diversity levels are insufficient to infer the within-population recombination rate of the virus since the pandemic began. Our work offers new methods for inferring recombination rates in RNA viruses with implications for understanding recombination in SARS-CoV-2 evolution and the structure of clonal relationships and gene pools shaping its origins.

The Spread of SARS-CoV-2 Omicron Variant in CALABRIA: A Spatio-Temporal Report of Viral Genome Evolution

We investigated the evolution of SARS-CoV-2 spread in Calabria, Southern Italy, in 2022. A total of 272 RNA isolates from nasopharyngeal swabs of individuals infected with SARS-CoV-2 were sequenced by whole genome sequencing (N = 172) and/or Sanger sequencing (N = 100). Analysis of diffusion of Omicron variants in Calabria revealed the prevalence of 10 different sub-lineages (recombinant BA.1/BA.2, BA.1, BA.1.1, BA.2, BA.2.9, BA.2.10, BA.2.12.1, BA.4, BA.5, BE.1). We observed that Omicron spread in Calabria presented a similar trend as in Italy, with some notable exceptions: BA.1 disappeared in April in Calabria but not in the rest of Italy; recombinant BA.1/BA.2 showed higher frequency in Calabria (13%) than in the rest of Italy (0.02%); BA.2.9, BA.4 and BA.5 emerged in Calabria later than in other Italian regions. In addition, Calabria Omicron presented 16 non-canonical mutations in the S protein and 151 non-canonical mutations in non-structural proteins. Most non-canonical mutations in the S protein occurred mainly in BA.5 whereas non-canonical mutations in non-structural or accessory proteins (ORF1ab, ORF3a, ORF8 and N) were identified in BA.2 and BA.5 sub-lineages. In conclusion, the data reported here underscore the importance of monitoring the entire SARS-CoV-2 genome.

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.

Non-uniform aspects of the SARS-CoV-2 intraspecies evolution reopen question of its origin

Several hypotheses have been presented on the origin of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) from its identification as the agent causing the current coronavirus disease 19 (COVID-19) pandemic. So far, no solid evidence has been found to support any hypothesis on the origin of this virus, and the issue continue to resurface over and over again. Here we have unfolded a pattern of distribution of several mutations in the SARS-CoV-2 proteins in 24 geo-locations across different continents. The results showed an evenly uneven distribution of the unique protein variants, distinct mutations, unique frequency of common conserved residues, and mutational residues across these 24 geo-locations. Furthermore, ample mutations were identified in the evolutionarily conserved invariant regions in the SARS-CoV-2 proteins across almost all geo-locations studied. This pattern of mutations potentially breaches the law of evolutionary conserved functional units of the beta-coronavirus genus. These mutations may lead to several novel SARS-CoV-2 variants with a high degree of transmissibility and virulence. A thorough investigation on the origin and characteristics of SARS-CoV-2 needs to be conducted in the interest of science and for the preparation of meeting the challenges of potential future pandemics.

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.

Temporal changes in the accessory protein mutations of SARS-CoV-2 variants and their predicted structural and functional effects

Emerging variants enable the continuous spread of SARS-CoV-2 in humans. The factors contributing to behavioral differences in variants remain elusive despite associations with several Spike protein mutations. Exploring accessory proteins may provide a wider understanding of these differences since these proteins may affect viral processes that occur beyond infection. Various bioinformatics tools were utilized to identify significant accessory protein mutations and determine their structural and functional effects over time. The ViruClust web application was used to retrieve accessory protein amino acid sequences and determine mutation frequencies in these sequences across time. The structural and functional effects of the mutations were determined using Missense3D and PROVEAN, respectively. The accessory and Spike protein mutations were compared using mutation densities. Q57H and T151I of ORF3a; T21I and W27L of ORF6; G38V, V82A, and T120I of ORF7a; S31P and T40I of ORF7b; and R52I, C61F, and I121L of ORF8 were highly frequent in most variants of concern and were within known functional domains. Thus, these are good candidates for further experimental evaluation. Among the accessory proteins, ORF6 and ORF8 were highlighted because of their strong and weak correlation with Spike protein mutations, respectively.

A neonatal mouse model characterizes transmissibility of SARS-CoV-2 variants and reveals a role for ORF8

Small animal models have been a challenge for the study of SARS-CoV-2 transmission, with most investigators using golden hamsters or ferrets ^1,2 . Mice have the advantages of low cost, wide availability, less regulatory and husbandry challenges, and the existence of a versatile reagent and genetic toolbox. However, adult mice do not transmit SARS-CoV-2 ³ . Here we establish a model based on neonatal mice that allows for transmission of clinical SARS-CoV-2 isolates. We characterize tropism, respiratory tract replication and transmission of ancestral WA-1 compared to variants alpha (B.1.1.7), beta (B.1.351), gamma (P.1), delta (B.1.617.2) and omicron (B.1.1.529). We identify inter-variant differences in timing and magnitude of infectious particle shedding from index mice, both of which shape transmission to contact mice. Furthermore, we characterize two recombinant SARS-CoV-2 lacking either the ORF6 or ORF8 host antagonists. The removal of ORF8 shifts viral replication towards the lower respiratory tract, resulting in significantly delayed and reduced transmission. Our results demonstrate the potential of our neonatal mouse model to characterize viral and host determinants of SARS-CoV-2 transmission, while revealing for the first time a role for an accessory protein this context.

SARS-CoV-2 ORF8: One protein, seemingly one structure, and many functions

SARS-CoV-2 is the virus responsible for the COVID-19 pandemic. The genome of SARS-CoV-2 encodes nine accessory proteins that are involved in host-pathogen interaction. ORF8 is unique among these accessory proteins. SARS-CoV-2 ORF8 shares a surprisingly low amino acid sequence similarity with SARS-COV ORF8 (30%), and it is presumed to have originated from bat. Studies have shown that ORF8 exerts multiple different functions that interfere with host immune responses, including the downregulation of MHC class I molecules. These functions may represent strategies of host immune evasion. The x-ray crystal structure of ORF8 revealed an immunoglobulin-like domain with several distinguishing features. To date, there are numerous unanswered questions about SARS-CoV-2 ORF8 protein and its structure-function relationship that we discuss in this mini-review. A better understanding of how ORF8 interacts with components of the immune system is needed for elucidating COVID-19 pathogenesis and to develop new avenues for the treatment of the disease.

Host genetic diversity and genetic variations of SARS-CoV-2 in COVID-19 pathogenesis and the effectiveness of vaccination

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), responsible for the outbreak of coronavirus disease 2019 (COVID-19), has shown a vast range of clinical manifestations from asymptomatic to life-threatening symptoms. To figure out the cause of this heterogeneity, studies demonstrated the trace of genetic diversities whether in the hosts or the virus itself. With this regard, this review provides a comprehensive overview of how host genetic such as those related to the entry of the virus, the immune-related genes, gender-related genes, disease-related genes, and also host epigenetic could influence the severity of COVID-19. Besides, the mutations in the genome of SARS-CoV-2 __leading to emerging of new variants__ per se affect the affinity of the virus to the host cells and enhance the immune escape capacity. The current review discusses these variants and also the latest data about vaccination effectiveness facing the most important variants.

SARS-CoV-2 modulation of RIG-I-MAVS signaling: Potential mechanisms of impairment on host antiviral immunity and therapeutic approaches

The coronavirus disease 2019 (COVID-19) is a global infectious disease aroused by RNA virus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Patients may suffer from severe respiratory failure or even die, posing a huge challenge to global public health. Retinoic acid-inducible gene I (RIG-I) is one of the major pattern recognition receptors, function to recognize RNA viruses and mediate the innate immune response. RIG-1 and melanoma differentiation-associated gene 5 contain an N-terminal caspase recruitment domain that is activated upon detection of viral RNA in the cytoplasm of virus-infected cells. Activated RIG-I and mitochondrial antiviral signaling (MAVS) protein trigger a series of corresponding immune responses such as the production of type I interferon against viral infection. In this review, we are summarizing the role of the structural, nonstructural, and accessory proteins from SARS-CoV-2 on the RIG-I-MAVS pathway, and exploring the potential mechanism how SARS-CoV-2 could evade the host antiviral response. We then proposed that modulation of the RIG-I-MAVS signaling pathway might be a novel and effective therapeutic strategy to against COVID-19 as well as the constantly mutating coronavirus.

A Putative long-range RNA-RNA interaction between ORF8 and Spike of SARS-CoV-2

SARS-CoV-2 has affected people worldwide as the causative agent of COVID-19. The virus is related to the highly lethal SARS-CoV-1 responsible for the 2002-2003 SARS outbreak in Asia. Research is ongoing to understand why both viruses have different spreading capacities and mortality rates. Like other beta coronaviruses, RNA-RNA interactions occur between different parts of the viral genomic RNA, resulting in discontinuous transcription and production of various sub-genomic RNAs. These sub-genomic RNAs are then translated into other viral proteins. In this work, we performed a comparative analysis for novel long-range RNA-RNA interactions that may involve the Spike region. Comparing in-silico fragment-based predictions between reference sequences of SARS-CoV-1 and SARS-CoV-2 revealed several predictions amongst which a thermodynamically stable long-range RNA-RNA interaction between (23660-23703 Spike) and (28025-28060 ORF8) unique to SARS-CoV-2 was observed. The patterns of sequence variation using data gathered worldwide further supported the predicted stability of the sub-interacting region (23679-23690 Spike) and (28031-28042 ORF8). Such RNA-RNA interactions can potentially impact viral life cycle including sub-genomic RNA production rates.

SARS-CoV-2 genome variations and evolution patterns in Egypt: a multi-center study

A serious global public health emergency emerged late November 2019 in Wuhan City, China, by a new highly pathogenic virus, SARS-CoV-2. The virus evolution spread has been tracked by three developing databases: GISAID, Nextstrain and PANGO to understand its circulating variants. In this study, 110 diagnosed positive COVID-19 patient's samples, were collected from Kasr Al-Aini Hospital and the Children Cancer Hospital Egypt 57357 between May 2020 and January 2021, with clinical severity ranging from mild to severe. The viral genomes were sequenced by next generation sequencing, and phylogenetic analysis was performed to understand viral transmission dynamics. According to Nextstrain clades, most of our sequenced samples belonged to clades 20A and 20D, which in addition to clade 20B were present from the beginning of sample collection in May 2020. Clades 19A and 19B, on the other hand, appeared in the mid and late 2020 respectively, followed by the disappearance of clade 20B at the end of 2020. We identified a relatively high prevalence of the D614G spike protein variant and novel patterns of mutations associated together and with different clades. We also identified four mutations, spike H49Y, ORF3a H78Y, ORF8 E64stop and nucleocapsid E378V, associated with higher disease severity. Altogether, our study contributes genetic, phylogenetic, and clinical correlation data about the spread of the SARS-CoV-2 pandemic in Egypt.

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.

Appraisal of SARS-CoV-2 mutations and their impact on vaccination efficacy: an overview

With the unexpected emergence of the novel 2019 Wuhan coronavirus, the world was faced with a sudden uproar that quickly shifted into a serious life-threatening pandemic. Affecting the lives of the global population and leaving drastic damage in various sections and systems, several measures have been constantly taken to tackle down this crisis. For instance, numerous vaccines have been developed in the past two years, some of which have been granted emergency use, thus providing sufficient immunity to the vaccinated individuals. However, the appearance of newly emerged SARS-CoV-2 variants with accelerated transmission and fatality has led the world towards another pandemic. Having undergone various mutations in genomic and/or amino acid profiles, some of the emerged variants of concern (VOCs) including Alpha, Beta, Gamma, and Delta have displayed immune evasion and pathogenicity even in the vaccinated population, hence raising concerns regarding the efficacy of current vaccines against new VOCs of COVID-19. Therefore, genomic investigations of SARS-CoV-2 mutations are expected to provide valuable insight into the evolution of SARS-CoV-2, while also determining the impact of different mutations on infection severity. This study was constructed with the aim of shining light on recent advances regarding mutations in major COVID-19 VOCs, as well as vaccination efficacy against those VOCs.

A Bioinformatics Approach to Investigate Structural and Non-Structural Proteins in Human Coronaviruses

Recent studies confirmed that people unexposed to SARS-CoV-2 have preexisting reactivity, probably due to previous exposure to widely circulating common cold coronaviruses. Such preexistent reactivity against SARS-CoV-2 comes from memory T cells that can specifically recognize a SARS-CoV-2 epitope of structural and non-structural proteins and the homologous epitopes from common cold coronaviruses. Therefore, it is important to understand the SARS-CoV-2 cross-reactivity by investigating these protein sequence similarities with those of different circulating coronaviruses. In addition, the emerging SARS-CoV-2 variants lead to an intense interest in whether mutations in proteins (especially in the spike) could potentially compromise vaccine effectiveness. Since it is not clear that the differences in clinical outcomes are caused by common cold coronaviruses, a deeper investigation on cross-reactive T-cell immunity to SARS-CoV-2 is crucial to examine the differential COVID-19 symptoms and vaccine performance. Therefore, the present study can be a starting point for further research on cross-reactive T cell recognition between circulating common cold coronaviruses and SARS-CoV-2, including the most recent variants Delta and Omicron. In the end, a deep learning approach, based on Siamese networks, is proposed to accurately and efficiently calculate a BLAST-like similarity score between protein sequences.

Insertion-and-Deletion Mutations between the Genomes of SARS-CoV, SARS-CoV-2, and Bat Coronavirus RaTG13

The evolutional process of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) development remains inconclusive. This study compared the genome sequences of severe acute respiratory syndrome coronavirus (SARS-CoV), bat coronavirus RaTG13, and SARS-CoV-2. In total, the genomes of SARS-CoV-2 and RaTG13 were 77.9% and 77.7% identical to the genome of SARS-CoV, respectively. A total of 3.6% (1,068 bases) of the SARS-CoV-2 genome was derived from insertion and/or deletion (indel) mutations, and 18.6% (5,548 bases) was from point mutations from the genome of SARS-CoV. At least 35 indel sites were confirmed in the genome of SARS-CoV-2, in which 17 were with ≥10 consecutive bases long. Ten of these relatively long indels were located in the spike (S) gene, five in nonstructural protein 3 (Nsp3) gene of open reading frame (ORF) 1a, and one in ORF8 and noncoding region. Seventeen (48.6%) of the 35 indels were based on insertion-and-deletion mutations with exchanged gene sequences of 7–325 consecutive bases. Almost the complete ORF8 gene was replaced by a single 325 consecutive base-long indel. The distribution of these indels was roughly in accordance with the distribution of the rate of point mutation rate around the indels. The genome sequence of SARS-CoV-2 was 96.0% identical to that of RaTG13. There was no long insertion-and-deletion mutation between the genomes of RaTG13 and SARS-CoV-2. The findings of the uneven distribution of multiple indels and the presence of multiple long insertion-and-deletion mutations with exchanged consecutive base sequences in the viral genome may provide insights into SARS-CoV-2 development.

IMPORTANCE The developmental mechanism of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains inconclusive. This study compared the base sequence one-by-one between severe acute respiratory syndrome coronavirus (SARS-CoV) or bat coronavirus RaTG13 and SARS-CoV-2. The genomes of SARS-CoV-2 and RaTG13 were 77.9% and 77.7% identical to the genome of SARS-CoV, respectively. Seventeen of the 35 sites with insertion and/or deletion mutations between SARS-CoV-2 and SARS-CoV were based on insertion-and-deletion mutations with the replacement of 7–325 consecutive bases. Most of these long insertion-and-deletion sites were concentrated in the nonstructural protein 3 (Nsp3) gene of open reading frame (ORF) 1a, S1 domain of the spike protein, and ORF8 genes. Such long insertion-and-deletion mutations were not observed between the genomes of RaTG13 and SARS-CoV-2. The presence of multiple long insertion-and-deletion mutations in the genome of SARS-CoV-2 and their uneven distributions may provide further insights into the development of the virus.

A bioinformatic approach of targeting SARS-CoV-2 replication by silencing a conserved alternative reserve of the orf8 gene using host miRNAs

The causative agent of the COVID-19 pandemic, the SARS-CoV-2 virus has yielded multiple relevant mutations, many of which have branched into major variants. The Omicron variant has a huge similarity with the original viral strain (first COVID-19 strain from Wuhan). Among different genes, the highly variable orf8 gene is responsible for crucial host interactions and has undergone multiple mutations and indels. The sequence of the orf8 gene of the Omicron variant is, however, identical with the gene sequence of the wild type. orf8 modulates the host immunity making it easier for the virus to conceal itself and remain undetected. Variants seem to be deleting this gene without affecting the viral replication. While analyzing, we came across the conserved orf7a gene in the viral genome which exhibits a partial sequence homology as well as functional similarity with the SARS-CoV-2 orf8. Hence, we have proposed here in our hypothesis that, orf7a might be an alternative reserve of orf8 present in the virus which was compensating for the lost gene. A computational approach was adopted where we screened various miRNAs targeted against the orf8 gene. These miRNAs were then docked onto the orf8 mRNA sequences. The same set of miRNAs was then used to check for their binding affinity with the orf7a reference mRNA. Results showed that miRNAs targeting the orf8 had favorable shape complementarity and successfully docked with the orf7a gene as well. These findings provide a basis for developing new therapeutic approaches where both orf8 and orf7a can be targeted simultaneously.

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.

Detection of SARS-CoV-2 Proteins in Wastewater Samples by Mass Spectrometry

The recent COVID-19 pandemic overwhelmed the health system worldwide, and there was a need to track outbreaks and try to use this information as an early warning system. Wastewater-based epidemiology (WBE) enabled detection of the SARS-CoV-2 virus in wastewater treatment plant influents. Until now, the most used technique for this detection has been the quantitative polymerase chain reaction (qPCR)-based quantification of SARS-CoV-2 RNA. This study proposes a mass spectrometry (MS)-based method that detected specific SARS-CoV-2 proteins in wastewater, 5 and 6 days ahead of the case data for two municipalities. We identified unique peptides of eight proteins related to the SARS-CoV-2 virus and COVID-19 infection. We detected the nonstructural protein (NSP) pp1ab (transcribed after host cell infection) most frequently in all of the samples. As a result, we suspect that in the active cases of COVID-19, the pp1ab protein is present in high abundance in the urine and feces and that this protein could be used as an alternative biomarker. These data were collected before mass vaccination occurred in the population.

An issue of concern: unique truncated ORF8 protein variants of SARS-CoV-2

Open reading frame 8 (ORF8) shows one of the highest levels of variability among accessory proteins in Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of Coronavirus Disease 2019 (COVID-19). It was previously reported that the ORF8 protein inhibits the presentation of viral antigens by the major histocompatibility complex class I (MHC-I), which interacts with host factors involved in pulmonary inflammation. The ORF8 protein assists SARS-CoV-2 in evading immunity and plays a role in SARS-CoV-2 replication. Among many contributing mutations, Q27STOP, a mutation in the ORF8 protein, defines the B.1.1.7 lineage of SARS-CoV-2, engendering the second wave of COVID-19. In the present study, 47 unique truncated ORF8 proteins (T-ORF8) with the Q27STOP mutations were identified among 49,055 available B.1.1.7 SARS-CoV-2 sequences. The results show that only one of the 47 T-ORF8 variants spread to over 57 geo-locations in North America, and other continents, which include Africa, Asia, Europe and South America. Based on various quantitative features, such as amino acid homology, polar/non-polar sequence homology, Shannon entropy conservation, and other physicochemical properties of all specific 47 T-ORF8 protein variants, nine possible T-ORF8 unique variants were defined. The question as to whether T-ORF8 variants function similarly to the wild type ORF8 is yet to be investigated. A positive response to the question could exacerbate future COVID-19 waves, necessitating severe containment measures.

Antibody escape and global spread of SARS-CoV-2 lineage A.27

In spring 2021, an increasing number of infections was observed caused by the hitherto rarely described SARS-CoV-2 variant A.27 in south-west Germany. From December 2020 to June 2021 this lineage has been detected in 31 countries. Phylogeographic analyses of A.27 sequences obtained from national and international databases reveal a global spread of this lineage through multiple introductions from its inferred origin in Western Africa. Variant A.27 is characterized by a mutational pattern in the spike gene that includes the L18F, L452R and N501Y spike amino acid substitutions found in various variants of concern but lacks the globally dominant D614G. Neutralization assays demonstrate an escape of A.27 from convalescent and vaccine-elicited antibody-mediated immunity. Moreover, the therapeutic monoclonal antibody Bamlanivimab and partially the REGN-COV2 cocktail fail to block infection by A.27. Our data emphasize the need for continued global monitoring of novel lineages because of the independent evolution of new escape mutations.

In silico Analysis of SARS-CoV-2 ORF8-Binding Proteins Reveals the Involvement of ORF8 in Acquired-Immune and Innate-Immune Systems

SARS-CoV-2 is the causative agent of a new type of coronavirus infection, COVID-19, which has rapidly spread worldwide. The overall genome sequence homology between SARS-CoV-2 and SARS-CoV is 79%. However, the homology of the ORF8 protein between these two coronaviruses is low, at ~26%. Previously, it has been suggested that infection by the ORF8-deleted variant of SARS-CoV-2 results in less severe symptoms than in the case of wild-type SARS-CoV-2. Although we found that ORF8 is involved in the proteasome autoimmunity system, the precise role of ORF8 in infection and pathology has not been fully clarified. In this study, we determined a new network of ORF8-interacting proteins by performing in silico analysis of the binding proteins against the previously described 47 ORF8-binding proteins. We used as a dataset 431 human protein candidates from Uniprot that physically interacted with 47 ORF8-binding proteins, as identified using STRING. Homology and phylogenetic profile analyses of the protein dataset were performed on 446 eukaryotic species whose genome sequences were available in KEGG OC. Based on the phylogenetic profile results, clustering analysis was performed using Ward's method. Our phylogenetic profiling showed that the interactors of the ORF8-interacting proteins were clustered into three classes that were conserved across chordates (Class 1: 152 proteins), metazoans (Class 2: 163 proteins), and eukaryotes (Class 3: 114 proteins). Following the KEGG pathway analysis, classification of cellular localization, tissue-specific expression analysis, and a literature study on each class of the phylogenetic profiling cluster tree, we predicted that the following: protein members in Class 1 could contribute to COVID-19 pathogenesis via complement and coagulation cascades and could promote sarcoidosis; the members of Class 1 and 2, together, may contribute to the downregulation of Interferon-β; and Class 3 proteins are associated with endoplasmic reticulum stress and the degradation of human leukocyte antigen.

Structural biology of SARS-CoV-2: open the door for novel therapies

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) is the causative agent of the pandemic disease COVID-19, which is so far without efficacious treatment. The discovery of therapy reagents for treating COVID-19 are urgently needed, and the structures of the potential drug-target proteins in the viral life cycle are particularly important. SARS-CoV-2, a member of the Orthocoronavirinae subfamily containing the largest RNA genome, encodes 29 proteins including nonstructural, structural and accessory proteins which are involved in viral adsorption, entry and uncoating, nucleic acid replication and transcription, assembly and release, etc. These proteins individually act as a partner of the replication machinery or involved in forming the complexes with host cellular factors to participate in the essential physiological activities. This review summarizes the representative structures and typically potential therapy agents that target SARS-CoV-2 or some critical proteins for viral pathogenesis, providing insights into the mechanisms underlying viral infection, prevention of infection, and treatment. Indeed, these studies open the door for COVID therapies, leading to ways to prevent and treat COVID-19, especially, treatment of the disease caused by the viral variants are imperative.

Genomic surveillance of SARS-CoV-2 in the state of Delaware reveals tremendous genomic diversity

The use of next generation sequencing is critical for the surveillance of severe acute respiratory syndrome coronavirus 2, SARS-CoV-2, transmission, as single base mutations have been identified with differences in infectivity. A total of 1,459 high quality samples were collected, sequenced, and analyzed in the state of Delaware, a location that offers a unique perspective on transmission given its proximity to large international airports on the east coast. Pangolin and Nextclade were used to classify these sequences into 16 unique clades and 88 lineages. A total of 411 samples belonging to the Alpha 20I/501Y.V1 (B.1.1.7) strain of concern were identified, as well as one sample belonging to Beta 20H/501.V2 (B.1.351), thirteen belonging to Epsilon 20C/S:452R (B.1.427/B.1.429), two belonging to Delta 20A/S:478K (B.1.617.2), and 15 belonging to Gamma 20J/501Y.V3 (p.1). A total of 2217 unique coding mutations were observed with an average of 17.7 coding mutations per genome. These data paired with continued sample collection and sequencing will give a deeper understanding of the spread of SARS-CoV-2 strains within Delaware and its surrounding areas.

A comprehensive evolutionary and epidemiological characterization of insertion and deletion mutations in SARS-CoV-2 genomes

SARS-CoV-2, which causes the current pandemic of respiratory illness, is evolving continuously and generating new variants. Nevertheless, most of the sequence analyses thus far focused on nucleotide substitutions despite the fact that insertions and deletions (indels) are equally important in the evolution of SARS-CoV-2. In this study, we analyzed 1,099,664 high-quality sequences of SARS-CoV-2 genomes to re-construct the evolutionary and epidemiological histories of indels. Our analysis revealed 289 circulating indel types (237 deletion and 52 insertion types, each represented by more than ten genomic sequences), among which eighteen were recurrent indel types, each represented by more than 500 genome sequences. Although indels were identified across the entire genome, most of them were identified in nsp6, S, ORF8, and N genes, among which ORF8 indel types had the highest frequencies of frameshift. Geographical and temporal analyses of these variants revealed a few alterations of dominant indel types, each accompanied by geographic expansion to different countries and continents, which resulted in the fixation of several types of indels in the field, including the current variants of concern. Evolutionary and structural analyses revealed that indels involving S N-terminal domain regions were linked to the 3/4 variants of concern, resulting in significantly altered S protein that might contribute to the selective advantage of the corresponding variant. In sum, our study highlights the important role of insertions and deletions in the evolution and spread of SARS-CoV-2.

An issue of concern: unique truncated ORF8 protein variants of SARS-CoV-2

Open reading frame 8 (ORF8) shows one of the highest levels of variability among accessory proteins in Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of Coronavirus Disease 2019 (COVID-19). It was previously reported that the ORF8 protein inhibits the presentation of viral antigens by the major histocompatibility complex class I (MHC-I), which interacts with host factors involved in pulmonary inflammation. The ORF8 protein assists SARS-CoV-2 in evading immunity and plays a role in SARS-CoV-2 replication. Among many contributing mutations, Q27STOP, a mutation in the ORF8 protein, defines the B.1.1.7 lineage of SARS-CoV-2, engendering the second wave of COVID-19. In the present study, 47 unique truncated ORF8 proteins (T-ORF8) with the Q27STOP mutations were identified among 49,055 available B.1.1.7 SARS-CoV-2 sequences. The results show that only one of the 47 T-ORF8 variants spread to over 57 geo-locations in North America, and other continents, which include Africa, Asia, Europe and South America. Based on various quantitative features, such as amino acid homology, polar/non-polar sequence homology, Shannon entropy conservation, and other physicochemical properties of all specific 47 T-ORF8 protein variants, nine possible T-ORF8 unique variants were defined. The question as to whether T-ORF8 variants function similarly to the wild type ORF8 is yet to be investigated. A positive response to the question could exacerbate future COVID-19 waves, necessitating severe containment measures.

In silico study on the effects of disulfide bonds in ORF8 of SARS-CoV-2

The COVID-19 epidemic, caused by virus SARS-CoV-2, has been a pandemic and threatening everyone's health in the past two years. In SARS-CoV-2, ORF8 is one of the most important accessory...

The importance of accessory protein variants in the pathogenicity of SARS-CoV-2

The coronavirus disease 2019 (COVID-19) is caused by the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS- CoV-2) with an estimated fatality rate of less than 1%. The SARS-CoV-2 accessory proteins ORF3a, ORF6, ORF7a, ORF7b, ORF8, and ORF10 possess putative functions to manipulate host immune mechanisms. These involve interferons, which appear as a consensus function, immune signaling receptor NLRP3 (NLR family pyrin domain-containing 3) inflammasome, and inflammatory cytokines such as interleukin 1β (IL-1β) and are critical in COVID-19 pathology. Outspread variations of each of the six accessory proteins were observed across six continents of all complete SARS-CoV-2 proteomes based on the data reported before November 2020. A decreasing order of percentage of unique variations in the accessory proteins was determined as ORF3a > ORF8 > ORF7a > ORF6 > ORF10 > ORF7b across all continents. The highest and lowest unique variations of ORF3a were observed in South America and Oceania, respectively. These findings suggest that the wide variations in accessory proteins seem to affect the pathogenicity of SARS-CoV-2.

Crystal Structures of Bat and Human Coronavirus ORF8 Protein Ig-Like Domain Provide Insights Into the Diversity of Immune Responses

ORF8 is a viral immunoglobulin-like (Ig-like) domain protein encoded by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA genome. It tends to evolve rapidly and interfere with immune responses. However, the structural characteristics of various coronavirus ORF8 proteins and their subsequent effects on biological functions remain unclear. Herein, we determined the crystal structures of SARS-CoV-2 ORF8 (S84) (one of the epidemic isoforms) and the bat coronavirus RaTG13 ORF8 variant at 1.62 Å and 1.76 Å resolution, respectively. Comparison of these ORF8 proteins demonstrates that the 62-77 residues in Ig-like domain of coronavirus ORF8 adopt different conformations. Combined with mutagenesis assays, the residue Cys20 of ORF8 is responsible for forming the covalent disulfide-linked dimer in crystal packing and in vitro biochemical conditions. Furthermore, immune cell-binding assays indicate that various ORF8 (SARS-CoV-2 ORF8 (L84), ORF8 (S84), and RaTG13 ORF8) proteins have different interaction capabilities with human CD14⁺ monocytes in human peripheral blood. These results provide new insights into the specific characteristics of various coronavirus ORF8 and suggest that ORF8 variants may influence disease-related immune responses.

Year-long COVID-19 infection reveals within-host evolution of SARS-CoV-2 in a patient with B cell depletion

Background

B-cell depleting therapies may lead to protracted disease and prolonged viral shedding in individuals infected with SARS-CoV-2. Viral persistence in the setting of immunosuppression raises concern for viral evolution.

Methods

Amplification of sub-genomic transcripts for the E gene (sgE) was done on nasopharyngeal samples over the course of 355 days in a patient infected with SARS-CoV-2 who had previously undergone CAR T cell therapy and had persistently positive SARS-CoV-2 nasopharyngeal swabs. Whole genome sequencing was performed on samples from the patient's original presentation and 10 months later.

Results

Over the course of almost a year, the virus accumulated a unique in-frame deletion in the amino-terminal domain of the spike protein, and complete deletion of ORF7b and ORF8, the first report of its kind in an immunocompromised patient. Also, minority variants that were identified in the early samples-reflecting the heterogeneity of the initial infection-were found to be fixed late in the infection. Remdesivir and high-titer convalescent plasma treatment were given, and the infection was eventually cleared after 335 days of infection.

Conclusions

The unique viral mutations found in this study highlight the importance of analyzing viral evolution in protracted SARS-CoV-2 infection, especially in immunosuppressed hosts, and the implication of these mutations in the emergence of viral variants.

Summary

We report an immunocompromised patient with persistent symptomatic SARS-CoV-2 infection for 335 days. During this time, the virus accumulated a unique in-frame deletion in the spike, and a complete deletion of ORF7b and ORF8 which is the first report of its kind in an immunocompromised patient.

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.

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.

Whole-genome analysis of SARS-CoV-2 samples indicate no tissue specific genetic adaptation of the virus in COVID-19 patients' upper and lower respiratory tract

Sample panels of SARS-CoV-2 cases were retrospectively whole-genome sequenced. In three individuals, samples of upper and lower respiratory tract resulted in identical sequences suggesting virus stability including the spike protein cleavage site. In a fourth case, low-level intra-host genomic evolution and a unique 5-nucleotide deletion was observed.

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.

Molecular modelling on SARS-CoV-2 papain-like protease: an integrated study with homology modelling, molecular docking, and molecular dynamics simulations

SARS-CoV-2 PLpro was investigated as a therapeutic target for potent antiviral drugs due to its essential role in not only viral replication but also in regulating the inborn immune response. Several computational approaches, including homology modelling, molecular docking, and molecular dynamics (MD) studies, were employed to search for promising drugs in treating SARS-CoV-2. Eighty-one compounds, sub-structurally similar to the antiviral drug, were used as potential inhibitors of PLpro. From our results, three complexes containing the ligands with Pubchem IDs: 153012995, 12149203, and 123608715 showed lower binding energies than the control (Ritonavir), indicating that they may become promising inhibitors for PLpro. MD was performed in a water solvent to validate the stability of the three complexes. All complexes achieved stable structure during the simulation as no significant fluctuations were observed in the validation parameters. Moreover, the binding energy for each complex was estimated using the MM-GBSA method. Complex 1 was the most stable structure based on the lowest binding energy score and its structure remained in a similar cavity with the docket snapshot. Based on our studies, three ligands were assumed to be potential inhibitors. The ligand of complex 1 may become the most promising antiviral drug against SARS-CoV-2 targeting PLpro.