Nucleocapsid mutations in SARS-CoV-2 augment replication and pathogenesis

A specific phosphorylation-dependent conformational switch in SARS-CoV-2 nucleocapsid protein inhibits RNA binding

The nucleocapsid protein of severe acute respiratory syndrome coronavirus 2 encapsidates the viral genome and is essential for viral function. The central disordered domain comprises a serine-arginine-rich (SR) region that is hyperphosphorylated in infected cells. This modification regulates function, although mechanistic details remain unknown. We use nuclear magnetic resonance to follow structural changes occurring during hyperphosphorylation by serine arginine protein kinase 1, glycogen synthase kinase 3, and casein kinase 1, that abolishes interaction with RNA. When eight approximately uniformly distributed sites have been phosphorylated, the SR domain binds the same interface as single-stranded RNA, resulting in complete inhibition of RNA binding. Phosphorylation by protein kinase A does not prevent RNA binding, indicating that the pattern resulting from physiologically relevant kinases is specific for inhibition. Long-range contacts between the RNA binding, linker, and dimerization domains are abrogated, phenomena possibly related to genome packaging and unpackaging. This study provides insight into the recruitment of specific host kinases to regulate viral function.

COVID-19 in a common woolly monkey (Lagothrix lagothricha): First evidence of fatal outcome in a nonhuman primate after natural SARS-CoV-2 infection

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was declared a pandemic by the World Health Organization in March 2020. Since then, viral spread from humans to animals has occurred worldwide. Nonhuman primates (NHPs) have been found to be susceptible to reverse-zoonosis transmission of SARS-CoV-2, but initial research suggested that platyrrhine primates are less susceptible than catarrhine primates. Here we report the natural SARS-CoV-2 infection of a common woolly monkey (Lagothrix lagothricha) from a wildlife rehabilitation center in Ecuador. The course of the disease, the eventual death of the specimen, and the pathological findings are described. Our results show the susceptibility of a new platyrrhine species to SARS-CoV-2 and provide evidence for the first time of a COVID-19-associated death in a naturally infected NHP. The putative route of transmission from humans, and implications for captive NHPs management, are also discussed. Given that common woolly monkeys are at risk of extinction in Ecuador, further understanding of the potential threat of SARS-CoV-2 to their health should be a conservation priority. A One Health approach is the best way to protect NHPs from a new virus in the same way that we would protect the human population.

SARS-CoV-2-related bat viruses evade human intrinsic immunity but lack efficient transmission capacity

Circulating bat coronaviruses represent a pandemic threat. However, our understanding of bat coronavirus pathogenesis and transmission potential is limited by the lack of phenotypically characterized strains. We created molecular clones for the two closest known relatives of SARS-CoV-2, BANAL-52 and BANAL-236. We demonstrated that BANAL-CoVs and SARS-CoV-2 have similar replication kinetics in human bronchial epithelial cells. However, BANAL-CoVs have impaired replication in human nasal epithelial cells and in the upper airway of mice. We also observed reduced pathogenesis in mice and diminished transmission in hamsters. Further, we observed that diverse bat coronaviruses evade interferon and downregulate major histocompatibility complex class I. Collectively, our study demonstrates that despite high genetic similarity across bat coronaviruses, prediction of pandemic potential of a virus necessitates functional characterization. Finally, the restriction of bat coronavirus replication in the upper airway highlights that transmission potential and innate immune restriction can be uncoupled in this high-risk family of emerging viruses.

Comparative Proteomics and Interactome Analysis of the SARS-CoV-2 Nucleocapsid Protein in Human and Bat Cell Lines

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of COVID-19 and responsible for the global coronavirus pandemic which started in 2019. Despite exhaustive efforts to trace its origins, including potential links with pangolins and bats, the precise origins of the virus remain unclear. Bats have been recognized as natural hosts for various coronaviruses, including the Middle East respiratory coronavirus (MERS-CoV) and the SARS-CoV. This study presents a comparative analysis of the SARS-CoV-2 nucleocapsid protein (N) interactome in human and bat cell lines. We identified approximately 168 cellular proteins as interacting partners of SARS-CoV-2 N in human cells and 196 cellular proteins as interacting partners with this protein in bat cells. The results highlight pathways and events that are both common and unique to either bat or human cells. Understanding these interactions is crucial to comprehend the reasons behind the remarkable resilience of bats to viral infections. This study provides a foundation for a deeper understanding of host-virus interactions in different reservoirs.

Modulation of biophysical properties of nucleocapsid protein in the mutant spectrum of SARS-CoV-2

Genetic diversity is a hallmark of RNA viruses and the basis for their evolutionary success. Taking advantage of the uniquely large genomic database of SARS-CoV-2, we examine the impact of mutations across the spectrum of viable amino acid sequences on the biophysical phenotypes of the highly expressed and multifunctional nucleocapsid protein. We find variation in the physicochemical parameters of its extended intrinsically disordered regions (IDRs) sufficient to allow local plasticity, but also observe functional constraints that similarly occur in related coronaviruses. In biophysical experiments with several N-protein species carrying mutations associated with major variants, we find that point mutations in the IDRs can have nonlocal impact and modulate thermodynamic stability, secondary structure, protein oligomeric state, particle formation, and liquid-liquid phase separation. In the Omicron variant, distant mutations in different IDRs have compensatory effects in shifting a delicate balance of interactions controlling protein assembly properties, and include the creation of a new protein-protein interaction interface in the N-terminal IDR through the defining P13L mutation. A picture emerges where genetic diversity is accompanied by significant variation in biophysical characteristics of functional N-protein species, in particular in the IDRs.

Assembly of SARS-CoV-2 nucleocapsid protein with nucleic acid

The viral genome of SARS-CoV-2 is packaged by the nucleocapsid (N-)protein into ribonucleoprotein particles (RNPs), 38 ± 10 of which are contained in each virion. Their architecture has remained unclear due to the pleomorphism of RNPs, the high flexibility of N-protein intrinsically disordered regions, and highly multivalent interactions between viral RNA and N-protein binding sites in both N-terminal (NTD) and C-terminal domain (CTD). Here we explore critical interaction motifs of RNPs by applying a combination of biophysical techniques to ancestral and mutant proteins binding different nucleic acids in an in vitro assay for RNP formation, and by examining nucleocapsid protein variants in a viral assembly assay. We find that nucleic acid-bound N-protein dimers oligomerize via a recently described protein-protein interface presented by a transient helix in its long disordered linker region between NTD and CTD. The resulting hexameric complexes are stabilized by multivalent protein-nucleic acid interactions that establish crosslinks between dimeric subunits. Assemblies are stabilized by the dimeric CTD of N-protein offering more than one binding site for stem-loop RNA. Our study suggests a model for RNP assembly where N-protein scaffolding at high density on viral RNA is followed by cooperative multimerization through protein-protein interactions in the disordered linker.

Differential host responses to COVID-19: Unraveling the complexity

These diverse outcomes of Covid-19 are influenced by various factors including age, gender, underlying health conditions, immune responses, viral variants, external factors, and overall quality of life. Demographic analysis of patients aged 0-18 years experienced mild to moderate cases, above 55 years with co-morbidities, were more severely affected.COVID-19 incidence was higher in males (58 %) & (42 %) in females. The reduced expression of Toll-like receptors (TLR) in severe and critical patients is a crucial determinant. This reduced TLR expression is primarily attributed to the dominance of the PLpro viral protein of COVID-19. Disease enrichment analysis highlights the long-term impact of COVID-19, which can lead to post-recovery complications such as hypertension, diabetes, cardiac diseases, and brain ischemia in Covid-19 patients. In conclusion, a comprehensive strategy targeting key factors like PLpro, TLR, and inflammatory cytokines such as IL-1 and IL-6 could offer an effective approach to mitigate the devastating effects of COVID-19.

A longitudinal study on SARS-CoV-2 seroconversion, reinfection and neutralisation spanning several variant waves and vaccination campaigns, Heinsberg, Germany, April 2020 to November 2022

BackgroundSince its emergence in December 2019, over 700 million people worldwide have been infected with SARS-CoV-2 up to May 2024. While early rollout of mRNA vaccines against COVID-19 has saved many lives, there was increasing immune escape of new virus variants. Longitudinal monitoring of population-wide SARS-CoV-2 antibody responses from regular sample collection irrespective of symptoms provides representative data on infection and seroconversion/seroreversion rates.AimTo examine adaptive and cellular immune responses of a German SARS-CoV-2 outbreak cohort through several waves of infection with different virus variants.MethodsUtilising a 31-month longitudinal seroepidemiological study (n = 1,446; mean age: 50 years, range: 2-103) initiated during the first SARS-CoV-2 superspreading event (February 2020) in Heinsberg, Germany, we analysed acute infection, seroconversion and virus neutralisation at five follow-up visits between October 2020 and November 2022; cellular and cross-protective immunity against SARS-CoV-2 Omicron variants were also examined.ResultsSARS-CoV-2 spike (S)-specific IgAs decreased shortly after infection, while IgGs remained stable. Both increased significantly after vaccination. We predict an 18-month half-life of S IgGs upon infection. Nucleocapsid (N)-specific responses declined over 12 months post-infection but increased (p < 0.0001) during Omicron. Frequencies of SARS-CoV-2-specific TNF-alpha+/IFN-gamma+ CD4+  T-cells declined over 12 months after infection (p < 0.01). SARS-CoV-2 S antibodies and neutralisation titres were highest in triple-vaccinated participants infected between April 2021 and November 2022 compared with infections between April 2020 and January 2021. Cross neutralisation against Omicron BQ.1.18 and XBB.1.5 was very low in all groups.ConclusionInfection and/or vaccination did not provide the population with cross-protection against Omicron variants.

Early, rapidly progressive vasculopathy in a transplanted heart: A possible complication of COVID-19

The epidemic caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has had a significant global impact, especially on immunosuppressed populations such as heart transplant recipients. While SARS-CoV-2 initially infects the respiratory system, cardiovascular complications induced by coronavirus disease 2019 (COVID-19) include cardiac arrest, myocardial infarction, heart failure, myocarditis, arrhythmia, acute myocyte injury, thrombotic events, and cardiogenic shock. Here, we present a case of a 45-year-old African American male who tested positive for COVID-19 infection six months after receiving a heart transplant. The patient was asymptomatic initially, but two weeks later he developed dyspnea, early satiety, and abdominal bloating. The patient was admitted to the hospital for acute renal failure and subsequently diagnosed with moderate acute T cell-mediated allograft rejection (Grade 2R) by endomyocardial biopsy. Three months after testing positive for COVID-19, the patient suffered a sudden cardiac death. At autopsy, the epicardium was diffusely edematous and showed vascular congestion. The coronary arteries showed a striking concentric narrowing of lumens and diffusely thickened arterial walls of all major extramural arteries deemed consistent with a rapidly progressive form of cardiac allograft vasculopathy (CAV). SARS-CoV-2 nucleocapsid protein was localized by immunohistochemistry (IHC) in endothelial cells of venules and capillaries within the epicardium. Our localization of SARS-CoV-2 in coronary vessel endothelial cells by IHC suggests that endothelial cell infection, endotheliitis, and immune-related inflammation may be a primary mechanism of vascular injury. The present case represents an early onset rapidly progressive form of CAV. This case may be the first case of post-transplant arteriopathy occurring in such a short time that includes corresponding autopsy, surgical pathology, and IHC data.

SARS-CoV-2 variant of concern fitness and adaptation in primary human airway epithelia

The severe acute respiratory syndrome coronavirus 2 pandemic is characterized by the emergence of novel variants of concern (VOCs) that replace ancestral strains. Here, we dissect the complex selective pressures by evaluating variant fitness and adaptation in human respiratory tissues. We evaluate viral properties and host responses to reconstruct forces behind D614G through Omicron (BA.1) emergence. We observe differential replication in airway epithelia, differences in cellular tropism, and virus-induced cytotoxicity. D614G accumulates the most mutations after infection, supporting zoonosis and adaptation to the human airway. We perform head-to-head competitions and observe the highest fitness for Gamma and Delta. Under these conditions, RNA recombination favors variants encoding the B.1.617.1 lineage 3' end. Based on viral growth kinetics, Alpha, Gamma, and Delta exhibit increased fitness compared to D614G. In contrast, the global success of Omicron likely derives from increased transmission and antigenic variation. Our data provide molecular evidence to support epidemiological observations of VOC emergence.

Molecular Insights into the Variability in Infection and Immune Evasion Capabilities of SARS-CoV-2 Variants: A Sequence and Structural Investigation of the RBD Domain

As severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants continuously emerge, an increasing number of mutations are accumulating in the Spike protein receptor-binding domain (RBD) region. Through sequence analysis of various Variants of Concern (VOC), we identified that they predominantly fall within the ο lineage although recent variants introduce any novel mutations in the RBD. Molecular dynamics simulations were employed to compute the binding free energy of these variants with human Angiotensin-converting enzyme 2 (ACE2). Structurally, the binding interface of the ο RBD displays a strong positive charge, complementing the negatively charged binding interface of ACE2, resulting in a significant enhancement in the electrostatic potential energy for the ο variants. Although the increased potential energy is partially offset by the rise in polar solvation free energy, enhanced electrostatic potential contributes to the long-range recognition between the ο variant's RBD and ACE2. We also conducted simulations of glycosylated ACE2-RBD proteins. The newly emerged ο (JN.1) variant has added a glycosylation site at N-354@RBD, which significantly weakened its binding affinity with ACE2. Further, our interaction studies with three monoclonal antibodies across multiple SARS-CoV-2 strains revealed a diminished neutralization efficacy against the ο variants, primarily attributed to the electrostatic repulsion between the antibodies and RBD interface. Considering the characteristics of the ο variant and the trajectory of emerging strains, we propose that newly developed antibodies against SARS-CoV-2 RBD should have surfaces rich in negative potential and, postbinding, exhibit strong van der Waals interactions. These findings provide invaluable guidance for the formulation of future therapeutic strategies.

Comparison of SARS-CoV-2 variants of concern in primary human nasal cultures demonstrates Delta as most cytopathic and Omicron as fastest replicating

The SARS-CoV-2 pandemic was marked with emerging viral variants, some of which were designated as variants of concern (VOCs) due to selection and rapid circulation in the human population. Here, we elucidate functional features of each VOC linked to variations in replication rate. Patient-derived primary nasal cultures grown at air-liquid interface were used to model upper respiratory infection and compared to cell lines derived from human lung epithelia. All VOCs replicated to higher titers than the ancestral virus, suggesting a selection for replication efficiency. In primary nasal cultures, Omicron replicated to the highest titers at early time points, followed by Delta, paralleling comparative studies of population sampling. All SARS-CoV-2 viruses entered the cell primarily via a transmembrane serine protease 2 (TMPRSS2)-dependent pathway, and Omicron was more likely to use an endosomal route of entry. All VOCs activated and overcame dsRNA-induced cellular responses, including interferon (IFN) signaling, oligoadenylate ribonuclease L degradation, and protein kinase R activation. Among the VOCs, Omicron infection induced expression of the most IFN and IFN-stimulated genes. Infections in nasal cultures resulted in cellular damage, including a compromise of cell barrier integrity and loss of nasal cilia and ciliary beating function, especially during Delta infection. Overall, Omicron was optimized for replication in the upper respiratory tract and least favorable in the lower respiratory cell line, and Delta was the most cytopathic for both upper and lower respiratory cells. Our findings highlight the functional differences among VOCs at the cellular level and imply distinct mechanisms of pathogenesis in infected individuals.

IMPORTANCE

Comparative analysis of infections by SARS-CoV-2 ancestral virus and variants of concern, including Alpha, Beta, Delta, and Omicron, indicated that variants were selected for efficiency in replication. In infections of patient-derived primary nasal cultures grown at air-liquid interface to model upper respiratory infection, Omicron reached the highest titers at early time points, a finding that was confirmed by parallel population sampling studies. While all infections overcame dsRNA-mediated host responses, infections with Omicron induced the strongest interferon and interferon-stimulated gene response. In both primary nasal cultures and lower respiratory cell line, infections by Delta were most damaging to the cells as indicated by syncytia formation, loss of cell barrier integrity, and nasal ciliary function.

Resistance mechanisms of SARS-CoV-2 3CLpro to the non-covalent inhibitor WU-04

Drug resistance poses a significant challenge in the development of effective therapies against SARS-CoV-2. Here, we identified two double mutations, M49K/M165V and M49K/S301P, in the 3C-like protease (3CLpro) that confer resistance to a novel non-covalent inhibitor, WU-04, which is currently in phase III clinical trials (NCT06197217). Crystallographic analysis indicates that the M49K mutation destabilizes the WU-04-binding pocket, impacting the binding of WU-04 more significantly than the binding of 3CLpro substrates. The M165V mutation directly interferes with WU-04 binding. The S301P mutation, which is far from the WU-04-binding pocket, indirectly affects WU-04 binding by restricting the rotation of 3CLpro's C-terminal tail and impeding 3CLpro dimerization. We further explored 3CLpro mutations that confer resistance to two clinically used inhibitors: ensitrelvir and nirmatrelvir, and revealed a trade-off between the catalytic activity, thermostability, and drug resistance of 3CLpro. We found that mutations at the same residue (M49) can have distinct effects on the 3CLpro inhibitors, highlighting the importance of developing multiple antiviral agents with different skeletons for fighting SARS-CoV-2. These findings enhance our understanding of SARS-CoV-2 resistance mechanisms and inform the development of effective therapeutics.

Sequence signatures within the genome of SARS-CoV-2 can be used to predict host source

We conducted an in silico analysis to better understand the potential factors impacting host adaptation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in white-tailed deer, humans, and mink due to the strong evidence of sustained transmission within these hosts. Classification models trained on single nucleotide and amino acid differences between samples effectively identified white-tailed deer-, human-, and mink-derived SARS-CoV-2. For example, the balanced accuracy score of Extremely Randomized Trees classifiers was 0.984 ± 0.006. Eighty-eight commonly identified predictive mutations are found at sites under strong positive and negative selective pressure. A large fraction of sites under selection (86.9%) or identified by machine learning (87.1%) are found in genes other than the spike. Some locations encoded by these gene regions are predicted to be B- and T-cell epitopes or are implicated in modulating the immune response suggesting that host adaptation may involve the evasion of the host immune system, modulation of the class-I major-histocompatibility complex, and the diminished recognition of immune epitopes by CD8+ T cells. Our selection and machine learning analysis also identified that silent mutations, such as C7303T and C9430T, play an important role in discriminating deer-derived samples across multiple clades. Finally, our investigation into the origin of the B.1.641 lineage from white-tailed deer in Canada discovered an additional human sequence from Michigan related to the B.1.641 lineage sampled near the emergence of this lineage. These findings demonstrate that machine-learning approaches can be used in combination with evolutionary genomics to identify factors possibly involved in the cross-species transmission of viruses and the emergence of novel viral lineages.IMPORTANCESevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a highly transmissible virus capable of infecting and establishing itself in human and wildlife populations, such as white-tailed deer. This fact highlights the importance of developing novel ways to identify genetic factors that contribute to its spread and adaptation to new host species. This is especially important since these populations can serve as reservoirs that potentially facilitate the re-introduction of new variants into human populations. In this study, we apply machine learning and phylogenetic methods to uncover biomarkers of SARS-CoV-2 adaptation in mink and white-tailed deer. We find evidence demonstrating that both non-synonymous and silent mutations can be used to differentiate animal-derived sequences from human-derived ones and each other. This evidence also suggests that host adaptation involves the evasion of the immune system and the suppression of antigen presentation. Finally, the methods developed here are general and can be used to investigate host adaptation in viruses other than SARS-CoV-2.

SARS-CoV-2 outbreak: role of viral proteins and genomic diversity in virus infection and COVID-19 progression

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection is the cause of coronavirus disease 2019 (COVID-19); a severe respiratory distress that has emerged from the city of Wuhan, Hubei province, China during December 2019. COVID-19 is currently the major global health problem and the disease has now spread to most countries in the world. COVID-19 has profoundly impacted human health and activities worldwide. Genetic mutation is one of the essential characteristics of viruses. They do so to adapt to their host or to move to another one. Viral genetic mutations have a high potentiality to impact human health as these mutations grant viruses unique unpredicted characteristics. The difficulty in predicting viral genetic mutations is a significant obstacle in the field. Evidence indicates that SARS-CoV-2 has a variety of genetic mutations and genomic diversity with obvious clinical consequences and implications. In this review, we comprehensively summarized and discussed the currently available knowledge regarding SARS-CoV-2 outbreaks with a fundamental focus on the role of the viral proteins and their mutations in viral infection and COVID-19 progression. We also summarized the clinical implications of SARS-CoV-2 variants and how they affect the disease severity and hinder vaccine development. Finally, we provided a massive phylogenetic analysis of the spike gene of 214 SARS-CoV-2 isolates from different geographical regions all over the world and their associated clinical implications.

Interaction between host G3BP and viral nucleocapsid protein regulates SARS-CoV-2 replication and pathogenicity

G3BP1/2 are paralogous proteins that promote stress granule formation in response to cellular stresses, including viral infection. The nucleocapsid (N) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) inhibits stress granule assembly and interacts with G3BP1/2 via an ITFG motif, including residue F17, in the N protein. Prior studies examining the impact of the G3PB1-N interaction on SARS-CoV-2 replication have produced inconsistent findings, and the role of this interaction in pathogenesis is unknown. Here, we use structural and biochemical analyses to define the residues required for G3BP1-N interaction and structure-guided mutagenesis to selectively disrupt this interaction. We find that N-F17A mutation causes highly specific loss of interaction with G3BP1/2. SARS-CoV-2 N-F17A fails to inhibit stress granule assembly in cells, has decreased viral replication, and causes decreased pathology in vivo. Further mechanistic studies indicate that the N-F17-mediated G3BP1-N interaction promotes infection by limiting sequestration of viral genomic RNA (gRNA) into stress granules.

Modulation of Biophysical Properties of Nucleocapsid Protein in the Mutant Spectrum of SARS-CoV-2

Genetic diversity is a hallmark of RNA viruses and the basis for their evolutionary success. Taking advantage of the uniquely large genomic database of SARS-CoV-2, we examine the impact of mutations across the spectrum of viable amino acid sequences on the biophysical phenotypes of the highly expressed and multifunctional nucleocapsid protein. We find variation in the physicochemical parameters of its extended intrinsically disordered regions (IDRs) sufficient to allow local plasticity, but also exhibiting functional constraints that similarly occur in related coronaviruses. In biophysical experiments with several N-protein species carrying mutations associated with major variants, we find that point mutations in the IDRs can have nonlocal impact and modulate thermodynamic stability, secondary structure, protein oligomeric state, particle formation, and liquid-liquid phase separation. In the Omicron variant, distant mutations in different IDRs have compensatory effects in shifting a delicate balance of interactions controlling protein assembly properties, and include the creation of a new protein-protein interaction interface in the N-terminal IDR through the defining P13L mutation. A picture emerges where genetic diversity is accompanied by significant variation in biophysical characteristics of functional N-protein species, in particular in the IDRs.

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.

Impact of mutations on the stability of SARS-CoV-2 nucleocapsid protein structure

The nucleocapsid (N) protein of SARS-CoV-2 is known to participate in various host cellular processes, including interferon inhibition, RNA interference, apoptosis, and regulation of virus life cycles. Additionally, it has potential as a diagnostic antigen and/or immunogen. Our research focuses on examining structural changes caused by mutations in the N protein. We have modeled the complete tertiary structure of native and mutated forms of the N protein using Alphafold2. Notably, the N protein contains 3 disordered regions. The focus was on investigating the impact of mutations on the stability of the protein's dimeric structure based on binding free energy calculations (MM-PB/GB-SA) and RMSD fluctuations after MD simulations. The results demonstrated that 28 mutations out of 37 selected mutations analyzed, compared with wild-type N protein, resulted in a stable dimeric structure, while 9 mutations led to destabilization. Our results are important to understand the tertiary structure of the N protein dimer of SARS-CoV-2 and the effect of mutations on it, their behavior in the host cell, as well as for the research of other viruses belonging to the same genus additionally, to anticipate potential strategies for addressing this viral illness․.

Variant mutation in SARS-CoV-2 nucleocapsid enhances viral infection via altered genomic encapsidation

The evolution of SARS-CoV-2 variants and their respective phenotypes represents an important set of tools to understand basic coronavirus biology as well as the public health implications of individual mutations in variants of concern. While mutations outside of Spike are not well studied, the entire viral genome is undergoing evolutionary selection, particularly the central disordered linker region of the nucleocapsid (N) protein. Here, we identify a mutation (G215C), characteristic of the Delta variant, that introduces a novel cysteine into this linker domain, which results in the formation of a disulfide bond and a stable N-N dimer. Using reverse genetics, we determined that this cysteine residue is necessary and sufficient for stable dimer formation in a WA1 SARS-CoV-2 background, where it results in significantly increased viral growth both in vitro and in vivo. Finally, we demonstrate that the N:G215C virus packages more nucleocapsid per virion and that individual virions are larger, with elongated morphologies.

Evolution of a Distinct SARS-CoV-2 Lineage Identified during an Investigation of a Hospital Outbreak

The SARS-CoV-2 virus steadily evolves, and numerous antigenically distinct variants have emerged over the past three years. Tracking the evolution of the virus would help us understand the process that generates the diverse variants and predict the future evolutionary trajectory of SARS-CoV-2. Here, we report the evolutionary trajectory of a unique Omicron lineage identified during an outbreak investigation that occurred in a residence unit in the healthcare system. The new lineage had four distinct non-synonymous and two distinct synonymous mutations apart from its parental lineage. Since this lineage of virus was exclusively found during the outbreak, we were able to track the detailed evolutionary history of the entire lineage along the transmission path. Furthermore, we estimated the evolutionary rate of the SARS-CoV-2 Omicron variant from the analysis of the evolution of the lineage. This new Omicron sub-lineage acquired 3 mutations in a 12-day period, and the evolutionary rate was estimated as 3.05 × 10-3 subs/site/year. This study provides more insight into an ever-evolving virus.

Epitopes recognition of SARS-CoV-2 nucleocapsid RNA binding domain by human monoclonal antibodies

Coronavirus nucleocapsid protein (NP) of SARS-CoV-2 plays a central role in many functions important for virus proliferation including packaging and protecting genomic RNA. The protein shares sequence, structure, and architecture with nucleocapsid proteins from betacoronaviruses. The N-terminal domain (NPRBD) binds RNA and the C-terminal domain is responsible for dimerization. After infection, NP is highly expressed and triggers robust host immune response. The anti-NP antibodies are not protective and not neutralizing but can effectively detect viral proliferation soon after infection. Two structures of SARS-CoV-2 NPRBD were determined providing a continuous model from residue 48 to 173, including RNA binding region and key epitopes. Five structures of NPRBD complexes with human mAbs were isolated using an antigen-bait sorting. Complexes revealed a distinct complement-determining regions and unique sets of epitope recognition. This may assist in the early detection of pathogens and designing peptide-based vaccines. Mutations that significantly increase viral load were mapped on developed, full length NP model, likely impacting interactions with host proteins and viral RNA.

Genomic Analyses of Major SARS-CoV-2 Variants Predicting Multiple Regions of Pathogenic and Transmissive Importance

The rapid evolution of SARS-CoV-2 has fueled its global proliferation since its discovery in 2019, with several notable variants having been responsible for increases in cases of coronavirus disease 2019 (COVID-19). Analyses of codon bias and usage in these variants between phylogenetic clades or lineages may grant insights into the evolution of SARS-CoV-2 and identify target codons indicative of evolutionary or mutative trends that may prove useful in tracking or defending oneself against emerging strains. We processed a cohort of 120 SARS-CoV-2 genome sequences through a statistical and bioinformatic pipeline to identify codons presenting evidence of selective pressure as well as codon coevolution. We report the identification of two codon sites in the orf8 and N genes demonstrating such evidence with real-world impacts on pathogenicity and transmissivity.

Analysis of SARS-CoV-2 genome evolutionary patterns

The spread of SARS-CoV-2 virus accompanied by public availability of abundant sequence data provides a window for the determination of viral evolutionary patterns. In this study, SARS-CoV-2 genome sequences were collected from seven countries in the period January 2020-December 2022. The sequences were classified into three phases, namely, pre-vaccination, post-vaccination, and recent period. Comparison was performed between these phases based on parameters like mutation rates, selection pressure (dN/dS ratio), and transition to transversion ratios (Ti/Tv). Similar comparisons were performed among SARS-CoV-2 variants. Statistical significance was tested using Graphpad unpaired t-test. The analysis showed an increase in the percent genomic mutation rates post-vaccination and in recent periods across all countries from the pre-vaccination sequences. Mutation rates were highest in NSP3, S, N, and NSP12b before and increased further after vaccination. NSP4 showed the largest change in mutation rates after vaccination. The dN/dS ratios showed purifying selection that shifted toward neutral selection after vaccination. N, ORF8, ORF3a, and ORF10 were under highest positive selection before vaccination. Shift toward neutral selection was driven by E, NSP3, and ORF7a in the after vaccination set. In recent sequences, the largest dN/dS change was observed in E, NSP1, and NSP13. The Ti/Tv ratios decreased with time. C→U and G→U were the most frequent transitions and transversions. However, U→G was the most frequent transversion in recent period. The Omicron variant had the highest genomic mutation rates, while Delta showed the highest dN/dS ratio. Protein-wise dN/dS ratio was also seen to vary across the different variants.IMPORTANCETo the best of our knowledge, there exists no other large-scale study of the genomic and protein-wise mutation patterns during the time course of evolution in different countries. Analyzing the SARS-CoV-2 evolutionary patterns in view of the varying spatial, temporal, and biological signals is important for diagnostics, therapeutics, and pharmacovigilance of SARS-CoV-2.

The NSP3 protein of SARS-CoV-2 binds fragile X mental retardation proteins to disrupt UBAP2L interactions

Viruses interact with numerous host factors to facilitate viral replication and to dampen antiviral defense mechanisms. We currently have a limited mechanistic understanding of how SARS-CoV-2 binds host factors and the functional role of these interactions. Here, we uncover a novel interaction between the viral NSP3 protein and the fragile X mental retardation proteins (FMRPs: FMR1, FXR1-2). SARS-CoV-2 NSP3 mutant viruses preventing FMRP binding have attenuated replication in vitro and reduced levels of viral antigen in lungs during the early stages of infection. We show that a unique peptide motif in NSP3 binds directly to the two central KH domains of FMRPs and that this interaction is disrupted by the I304N mutation found in a patient with fragile X syndrome. NSP3 binding to FMRPs disrupts their interaction with the stress granule component UBAP2L through direct competition with a peptide motif in UBAP2L to prevent FMRP incorporation into stress granules. Collectively, our results provide novel insight into how SARS-CoV-2 hijacks host cell proteins and provides molecular insight into the possible underlying molecular defects in fragile X syndrome.

Evolution of the SARS-CoV-2 Omicron Variants: Genetic Impact on Viral Fitness

Over the last three years, the pandemic of COVID-19 has had a significant impact on people's lives and the global economy. The incessant emergence of variant strains has compounded the challenges associated with the management of COVID-19. As the predominant variant from late 2021 to the present, Omicron and its sublineages, through continuous evolution, have demonstrated iterative viral fitness. The comprehensive elucidation of the biological implications that catalyzed this evolution remains incomplete. In accordance with extant research evidence, we provide a comprehensive review of subvariants of Omicron, delineating alterations in immune evasion, cellular infectivity, and the cross-species transmission potential. This review seeks to clarify the underpinnings of biology within the evolution of SARS-CoV-2, thereby providing a foundation for strategic considerations in the post-pandemic era of COVID-19.

Within-host evolutionary dynamics and tissue compartmentalization during acute SARS-CoV-2 infection

The global evolution of SARS-CoV-2 depends in part upon the evolutionary dynamics within individual hosts with varying immune histories. To characterize the within-host evolution of acute SARS-CoV-2 infection, we sequenced saliva and nasal samples collected daily from vaccinated and unvaccinated individuals early during infection. We show that longitudinal sampling facilitates high-confidence genetic variant detection and reveals evolutionary dynamics missed by less-frequent sampling strategies. Within-host dynamics in both unvaccinated and vaccinated individuals appeared largely stochastic; however, in rare cases, minor genetic variants emerged to frequencies sufficient for forward transmission. Finally, we detected significant genetic compartmentalization of viral variants between saliva and nasal swab sample sites in many individuals. Altogether, these data provide a high-resolution profile of within-host SARS-CoV-2 evolutionary dynamics.IMPORTANCEWe detail the within-host evolutionary dynamics of SARS-CoV-2 during acute infection in 31 individuals using daily longitudinal sampling. We characterized patterns of mutational accumulation for unvaccinated and vaccinated individuals, and observed that temporal variant dynamics in both groups were largely stochastic. Comparison of paired nasal and saliva samples also revealed significant genetic compartmentalization between tissue environments in multiple individuals. Our results demonstrate how selection, genetic drift, and spatial compartmentalization all play important roles in shaping the within-host evolution of SARS-CoV-2 populations during acute infection.

Machine learning-based approach KEVOLVE efficiently identifies SARS-CoV-2 variant-specific genomic signatures

Machine learning was shown to be effective at identifying distinctive genomic signatures among viral sequences. These signatures are defined as pervasive motifs in the viral genome that allow discrimination between species or variants. In the context of SARS-CoV-2, the identification of these signatures can assist in taxonomic and phylogenetic studies, improve in the recognition and definition of emerging variants, and aid in the characterization of functional properties of polymorphic gene products. In this paper, we assess KEVOLVE, an approach based on a genetic algorithm with a machine-learning kernel, to identify multiple genomic signatures based on minimal sets of k-mers. In a comparative study, in which we analyzed large SARS-CoV-2 genome dataset, KEVOLVE was more effective at identifying variant-discriminative signatures than several gold-standard statistical tools. Subsequently, these signatures were characterized using a new extension of KEVOLVE (KANALYZER) to highlight variations of the discriminative signatures among different classes of variants, their genomic location, and the mutations involved. The majority of identified signatures were associated with known mutations among the different variants, in terms of functional and pathological impact based on available literature. Here we showed that KEVOLVE is a robust machine learning approach to identify discriminative signatures among SARS-CoV-2 variants, which are frequently also biologically relevant, while bypassing multiple sequence alignments. The source code of the method and additional resources are available at: https://github.com/bioinfoUQAM/KEVOLVE.

Targeted amplification and genetic sequencing of the severe acute respiratory syndrome coronavirus 2 surface glycoprotein

IMPORTANCE

The COVID-19 pandemic was accompanied by an unprecedented surveillance effort. The resulting data were and will continue to be critical for surveillance and control of SARS-CoV-2. However, some genomic surveillance methods experienced challenges as the virus evolved, resulting in incomplete and poor quality data. Complete and quality coverage, especially of the S-gene, is important for supporting the selection of vaccine candidates. As such, we developed a robust method to target the S-gene for amplification and sequencing. By focusing on the S-gene and imposing strict coverage and quality metrics, we hope to increase the quality of surveillance data for this continually evolving gene. Our technique is currently being deployed globally to partner laboratories, and public health representatives from 79 countries have received hands-on training and support. Expanding access to quality surveillance methods will undoubtedly lead to earlier detection of novel variants and better inform vaccine strain selection.

TRIM28-mediated nucleocapsid protein SUMOylation enhances SARS-CoV-2 virulence

Viruses, as opportunistic intracellular parasites, hijack the cellular machinery of host cells to support their survival and propagation. Numerous viral proteins are subjected to host-mediated post-translational modifications. Here, we demonstrate that the SARS-CoV-2 nucleocapsid protein (SARS2-NP) is SUMOylated on the lysine 65 residue, which efficiently mediates SARS2-NP's ability in homo-oligomerization, RNA association, liquid-liquid phase separation (LLPS). Thereby the innate antiviral immune response is suppressed robustly. These roles can be achieved through intermolecular association between SUMO conjugation and a newly identified SUMO-interacting motif in SARS2-NP. Importantly, the widespread SARS2-NP R203K mutation gains a novel site of SUMOylation which further increases SARS2-NP's LLPS and immunosuppression. Notably, the SUMO E3 ligase TRIM28 is responsible for catalyzing SARS2-NP SUMOylation. An interfering peptide targeting the TRIM28 and SARS2-NP interaction was screened out to block SARS2-NP SUMOylation and LLPS, and consequently inhibit SARS-CoV-2 replication and rescue innate antiviral immunity. Collectively, these data support SARS2-NP SUMOylation is critical for SARS-CoV-2 virulence, and therefore provide a strategy to antagonize SARS-CoV-2.

Tracking SARS-CoV-2 variants during the 2023 flu season and beyond in Lebanon

BACKGROUND

Early SARS-CoV-2 variant detection relies on testing and genomic surveillance. The Omicron variant (B.1.1.529) has quickly become the dominant type among the previous circulating variants worldwide. Several subvariants have emerged exhibiting greater infectivity and immune evasion. In this study we aimed at studying the prevalence of the Omicron subvariants during the flu season and beyond in Lebanon through genomic screening and at determining the overall standing and trajectory of the pandemic in the country.

METHODS

A total of 155 SARS-CoV-2 RNA samples were sequenced, using Nanopore sequencing technology.

RESULTS

Nanopore sequencing of 155 genomes revealed their distribution over 39 Omicron variants. XBB.1.5 (23.29 %) was the most common, followed by XBB.1.9.1 (10.96 %) and XBB.1.42 (7.5 %). The first batch collected between September and November 2022, included the BA.2.75.2, BA.5.2, BA.5.2.20, BA.5.2.25 and BQ.1.1.5 lineages. Between December 2022 and January 2023, those lineages were replaced by BA.2.75.5, BN.1, BN.1.4, BQ.1, BQ.1.1, BQ.1.1.23, CH.1.1, CM.4 and XBK. Starting February 2023, we observed a gradual emergence and dominance of the recombinant XBB and its sub-lineages (XBB.1, XBB.1.5, XBB.1.5.2, XBB.1.5.3, XBB.1.9, XBB.1.9.1, XBB.1.9.2, XBB.1.16, XBB.1.22 and XBB.1.42).

CONCLUSIONS

The timely detection and characterization of SARS-CoV-2 variants is important to reduce transmission through established disease control measures and to avoid introductions into animal populations that could lead to serious public health implications.

Comparison of SARS-CoV-2 variants of concern in primary human nasal cultures demonstrates Delta as most cytopathic and Omicron as fastest replicating

The SARS-CoV-2 pandemic was marked with emerging viral variants, some of which were designated as variants of concern (VOCs) due to selection and rapid circulation in the human population. Here we elucidate functional features of each VOC linked to variations in replication rate. Patient-derived primary nasal cultures grown at air-liquid-interface (ALI) were used to model upper-respiratory infection and human lung epithelial cell lines used to model lower-respiratory infection. All VOCs replicated to higher titers than the ancestral virus, suggesting a selection for replication efficiency. In primary nasal cultures, Omicron replicated to the highest titers at early time points, followed by Delta, paralleling comparative studies of population sampling. All SARS-CoV-2 viruses entered the cell primarily via a transmembrane serine protease 2 (TMPRSS2)-dependent pathway, and Omicron was more likely to use an endosomal route of entry. All VOCs activated and overcame dsRNA-induced cellular responses including interferon (IFN) signaling, oligoadenylate ribonuclease L degradation and protein kinase R activation. Among the VOCs, Omicron infection induced expression of the most IFN and IFN stimulated genes. Infections in nasal cultures resulted in cellular damage, including a compromise of cell-barrier integrity and loss of nasal cilia and ciliary beating function, especially during Delta infection. Overall, Omicron was optimized for replication in the upper-respiratory system and least-favorable in the lower-respiratory cell line; and Delta was the most cytopathic for both upper and lower respiratory cells. Our findings highlight the functional differences among VOCs at the cellular level and imply distinct mechanisms of pathogenesis in infected individuals.

Introduction, Dispersal, and Predominance of SARS-CoV-2 Delta Variant in Rio Grande do Sul, Brazil: A Retrospective Analysis

Mutations in the SARS-CoV-2 genome can alter the virus' fitness, leading to the emergence of variants of concern (VOC). In Brazil, the Gamma variant dominated the pandemic in the first half of 2021, and from June onwards, the first cases of Delta infection were documented. Here, we investigate the introduction and dispersal of the Delta variant in the RS state by sequencing 1077 SARS-CoV-2-positive samples from June to October 2021. Of these samples, 34.7% were identified as Gamma and 65.3% as Delta. Notably, 99.2% of Delta sequences were clustered within the 21J lineage, forming a significant Brazilian clade. The estimated clock rate was 5.97 × 10-4 substitutions per site per year. The Delta variant was first reported on 17 June in the Vinhedos Basalto microregion and rapidly spread, accounting for over 70% of cases within nine weeks. Despite this, the number of cases and deaths remained stable, possibly due to vaccination, prior infections, and the continued mandatory mask use. In conclusion, our study provides insights into the Delta variant circulating in the RS state, highlighting the importance of genomic surveillance for monitoring viral evolution, even when the impact of new variants may be less severe in a given region.

Overview of Nucleocapsid-Targeting Vaccines against COVID-19

The new SARS-CoV-2 coronavirus, which emerged in late 2019, is a highly variable causative agent of COVID-19, a contagious respiratory disease with potentially severe complications. Vaccination is considered the most effective measure to prevent the spread and complications of this infection. Spike (S) protein-based vaccines were very successful in preventing COVID-19 caused by the ancestral SARS-CoV-2 strain; however, their efficacy was significantly reduced when coronavirus variants antigenically different from the original strain emerged in circulation. This is due to the high variability of this major viral antigen caused by escape from the immunity caused by the infection or vaccination with spike-targeting vaccines. The nucleocapsid protein (N) is a much more conserved SARS-CoV-2 antigen than the spike protein and has therefore attracted the attention of scientists as a promising target for broad-spectrum vaccine development. Here, we summarized the current data on various N-based COVID-19 vaccines that have been tested in animal challenge models or clinical trials. Despite the high conservatism of the N protein, escape mutations gradually occurring in the N sequence can affect its protective properties. During the three years of the pandemic, at least 12 mutations have arisen in the N sequence, affecting more than 40 known immunogenic T-cell epitopes, so the antigenicity of the N protein of recent SARS-CoV-2 variants may be altered. This fact should be taken into account as a limitation in the development of cross-reactive vaccines based on N-protein.

Assembly of SARS-CoV-2 ribonucleosomes by truncated N∗ variant of the nucleocapsid protein

The nucleocapsid (N) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) compacts the RNA genome into viral ribonucleoprotein (vRNP) complexes within virions. Assembly of vRNPs is inhibited by phosphorylation of the N protein serine/arginine (SR) region. Several SARS-CoV-2 variants of concern carry N protein mutations that reduce phosphorylation and enhance the efficiency of viral packaging. Variants of the dominant B.1.1 viral lineage also encode a truncated N protein, termed N∗ or Δ(1-209), that mediates genome packaging despite lacking the N-terminal RNA-binding domain and SR region. Here, we use mass photometry and negative stain electron microscopy to show that purified Δ(1-209) and viral RNA assemble into vRNPs that are remarkably similar in size and shape to those formed with full-length N protein. We show that assembly of Δ(1-209) vRNPs requires the leucine-rich helix of the central disordered region and that this helix promotes N protein oligomerization. We also find that fusion of a phosphomimetic SR region to Δ(1-209) inhibits RNA binding and vRNP assembly. Our results provide new insights into the mechanisms by which RNA binding promotes N protein self-association and vRNP assembly, and how this process is modulated by phosphorylation.

Cross-neutralization and viral fitness of SARS-CoV-2 Omicron sublineages

The rapid evolution of SARS-CoV-2 Omicron sublineages mandates a better understanding of viral replication and cross-neutralization among these sublineages. Here we used K18-hACE2 mice and primary human airway cultures to examine the viral fitness and antigenic relationship among Omicron sublineages. In both K18-hACE2 mice and human airway cultures, Omicron sublineages exhibited a replication order of BA.5 ≥ BA.2 ≥ BA.2.12.1 > BA.1; no difference in body weight loss was observed among different sublineage-infected mice. The BA.1-, BA.2-, BA.2.12.1-, and BA.5-infected mice developed distinguishable cross-neutralizations against Omicron sublineages, but exhibited little neutralization against the index virus (i.e. USA-WA1/2020) or the Delta variant. Surprisingly, the BA.5-infected mice developed higher neutralization activity against heterologous BA.2 and BA.2.12.1 than that against homologous BA.5; serum neutralizing titres did not always correlate with viral replication levels in infected animals. Our results revealed a distinct antigenic cartography of Omicron sublineages and support the bivalent vaccine approach.

Less neutralization evasion of SARS-CoV-2 BA.2.86 than XBB sublineages and CH.1.1

The highly mutated BA.2.86, with over 30 spike protein mutations in comparison to Omicron BA.2 and XBB.1.5 variants, has raised concerns about its potential to evade COVID-19 vaccination or prior SARS-CoV-2 infection-elicited immunity. In this study, we employ a live SARS-CoV-2 neutralization assay to compare the neutralization evasion ability of BA.2.86 with other emerged SARS-CoV-2 subvariants, including BA.2-derived CH.1.1, Delta-Omicron recombinant XBC.1.6, and XBB descendants XBB.1.5, XBB.1.16, XBB.2.3, EG.5.1 and FL.1.5.1. Our results show that BA.2.86 is less neutralization evasive than XBB sublineages. XBB descendants XBB.1.16, EG.5.1, and FL.1.5.1 continue to significantly evade neutralization induced by the parental COVID-19 mRNA vaccine and a BA.5 Bivalent booster. Notably, when compared to XBB.1.5, the more recent XBB descendants, particularly EG.5.1, display increased resistance to neutralization. Among all the tested variants, CH.1.1 exhibits the greatest neutralization evasion. In contrast, XBC.1.6 shows a slight reduction but remains comparably sensitive to neutralization when compared to BA.5. Furthermore, a recent XBB.1.5-breakthrough infection significantly enhances the breadth and potency of cross-neutralization. These findings reinforce the expectation that the upcoming XBB.1.5 mRNA vaccine would likely boost the neutralization of currently circulating variants, while also underscoring the critical importance of ongoing surveillance to monitor the evolution and immune evasion potential of SARS-CoV-2 variants.

Reverse genetics systems for SARS-CoV-2: Development and applications

The recent emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) caused serious harm to human health and struck a blow to global economic development. Research on SARS-CoV-2 has greatly benefited from the use of reverse genetics systems, which have been established to artificially manipulate the viral genome, generating recombinant and reporter infectious viruses or biosafety level 2 (BSL-2)-adapted non-infectious replicons with desired modifications. These tools have been instrumental in studying the molecular biological characteristics of the virus, investigating antiviral therapeutics, and facilitating the development of attenuated vaccine candidates. Here, we review the construction strategies, development, and applications of reverse genetics systems for SARS-CoV-2, which may be applied to other CoVs as well.

Assembly reactions of SARS-CoV-2 nucleocapsid protein with nucleic acid

The viral genome of SARS-CoV-2 is packaged by the nucleocapsid (N-) protein into ribonucleoprotein particles (RNPs), 38±10 of which are contained in each virion. Their architecture has remained unclear due to the pleomorphism of RNPs, the high flexibility of N-protein intrinsically disordered regions, and highly multivalent interactions between viral RNA and N-protein binding sites in both N-terminal (NTD) and C-terminal domain (CTD). Here we explore critical interaction motifs of RNPs by applying a combination of biophysical techniques to mutant proteins binding different nucleic acids in an in vitro assay for RNP formation, and by examining mutant proteins in a viral assembly assay. We find that nucleic acid-bound N-protein dimers oligomerize via a recently described protein-protein interface presented by a transient helix in its long disordered linker region between NTD and CTD. The resulting hexameric complexes are stabilized by multi-valent protein-nucleic acid interactions that establish crosslinks between dimeric subunits. Assemblies are stabilized by the dimeric CTD of N-protein offering more than one binding site for stem-loop RNA. Our study suggests a model for RNP assembly where N-protein scaffolding at high density on viral RNA is followed by cooperative multimerization through protein-protein interactions in the disordered linker.

SARS-CoV-2 variants evolve convergent strategies to remodel the host response

SARS-CoV-2 variants of concern (VOCs) emerged during the COVID-19 pandemic. Here, we used unbiased systems approaches to study the host-selective forces driving VOC evolution. We discovered that VOCs evolved convergent strategies to remodel the host by modulating viral RNA and protein levels, altering viral and host protein phosphorylation, and rewiring virus-host protein-protein interactions. Integrative computational analyses revealed that although Alpha, Beta, Gamma, and Delta ultimately converged to suppress interferon-stimulated genes (ISGs), Omicron BA.1 did not. ISG suppression correlated with the expression of viral innate immune antagonist proteins, including Orf6, N, and Orf9b, which we mapped to specific mutations. Later Omicron subvariants BA.4 and BA.5 more potently suppressed innate immunity than early subvariant BA.1, which correlated with Orf6 levels, although muted in BA.4 by a mutation that disrupts the Orf6-nuclear pore interaction. Our findings suggest that SARS-CoV-2 convergent evolution overcame human adaptive and innate immune barriers, laying the groundwork to tackle future pandemics.

Clinical Utility of Sero-Immunological Responses Against SARS-CoV-2 Nucleocapsid Protein During Subsequent Prevalence of Wild-Type, Delta Variant, and Omicron Variant

As nucleocapsid protein of severe acute respiratory syndrome coronavirus 2 is immunogenic but not targeted in vaccines, it could be useful in distinguishing natural infection from vaccination. We aimed to investigate the clinical utility of sero-immunological responses against the nucleocapsid protein. Nucleocapsid antibody immunoassay study with 302 coronavirus disease 2019 (COVID-19) patients showed lower titers in immunocompromised patients (P < 0.001), higher titers in higher severity (P = 0.031), and different seroconversion rates and titers according to variants of concern. Longitudinal evaluation of nucleocapsid antibodies using 513 samples from 291 COVID-19 patients revealed that it could persist up to 556 days from symptom onset. Interferon gamma release assay against the nucleocapsid protein showed poor response, precluding the deduction of a cut-off for the nucleocapsid protein. In conclusion, nucleocapsid antibody provides instructive clues about the immunogenicity of nucleocapsid proteins by different seroconversion rates and titers according to the severity of infection, host immune status, and different variants of concern.

SARS-CoV-2 hijacks fragile X mental retardation proteins for efficient infection

Viruses interact with numerous host factors to facilitate viral replication and to dampen antiviral defense mechanisms. We currently have a limited mechanistic understanding of how SARS-CoV-2 binds host factors and the functional role of these interactions. Here, we uncover a novel interaction between the viral NSP3 protein and the fragile X mental retardation proteins (FMRPs: FMR1 and FXR1-2). SARS-CoV-2 NSP3 mutant viruses preventing FMRP binding have attenuated replication in vitro and have delayed disease onset in vivo. We show that a unique peptide motif in NSP3 binds directly to the two central KH domains of FMRPs and that this interaction is disrupted by the I304N mutation found in a patient with fragile X syndrome. NSP3 binding to FMRPs disrupts their interaction with the stress granule component UBAP2L through direct competition with a peptide motif in UBAP2L to prevent FMRP incorporation into stress granules. Collectively, our results provide novel insight into how SARS-CoV-2 hijacks host cell proteins for efficient infection and provides molecular insight to the possible underlying molecular defects in fragile X syndrome.

Evolution of transient RNA structure-RNA polymerase interactions in respiratory RNA virus genomes

RNA viruses are important human pathogens that cause seasonal epidemics and occasional pandemics. Examples are influenza A viruses (IAV) and coronaviruses (CoV). When emerging IAV and CoV spill over to humans, they adapt to evade immune responses and optimize their replication and spread in human cells. In IAV, adaptation occurs in all viral proteins, including the viral ribonucleoprotein (RNP) complex. RNPs consist of a copy of the viral RNA polymerase, a double-helical coil of nucleoprotein, and one of the eight segments of the IAV RNA genome. The RNA segments and their transcripts are partially structured to coordinate the packaging of the viral genome and modulate viral mRNA translation. In addition, RNA structures can affect the efficiency of viral RNA synthesis and the activation of host innate immune response. Here, we investigated if RNA structures that modulate IAV replication processivity, so-called template loops (t-loops), vary during the adaptation of pandemic and emerging IAV to humans. Using cell culture-based replication assays and in silico sequence analyses, we find that the sensitivity of the IAV H3N2 RNA polymerase to t-loops increased between isolates from 1968 and 2017, whereas the total free energy of t-loops in the IAV H3N2 genome was reduced. This reduction is particularly prominent in the PB1 gene. In H1N1 IAV, we find two separate reductions in t-loop free energy, one following the 1918 pandemic and one following the 2009 pandemic. No destabilization of t-loops is observed in the influenza B virus genome, whereas analysis of SARS-CoV-2 isolates reveals destabilization of viral RNA structures. Overall, we propose that a loss of free energy in the RNA genome of emerging respiratory RNA viruses may contribute to the adaption of these viruses to the human population.

Total and Subgenomic RNA Viral Load in Patients Infected With SARS-CoV-2 Alpha, Delta, and Omicron Variants

BACKGROUND

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomic and subgenomic RNA levels are frequently used as a correlate of infectiousness. The impact of host factors and SARS-CoV-2 lineage on RNA viral load is unclear.

METHODS

Total nucleocapsid (N) and subgenomic N (sgN) RNA levels were measured by quantitative reverse transcription polymerase chain reaction (RT-qPCR) in specimens from 3204 individuals hospitalized with coronavirus disease 2019 (COVID-19) at 21 hospitals. RT-qPCR cycle threshold (Ct) values were used to estimate RNA viral load. The impact of time of sampling, SARS-CoV-2 variant, age, comorbidities, vaccination, and immune status on N and sgN Ct values were evaluated using multiple linear regression.

RESULTS

Mean Ct values at presentation for N were 24.14 (SD 4.53) for non-variants of concern, 25.15 (SD 4.33) for Alpha, 25.31 (SD 4.50) for Delta, and 26.26 (SD 4.42) for Omicron. N and sgN RNA levels varied with time since symptom onset and infecting variant but not with age, comorbidity, immune status, or vaccination. When normalized to total N RNA, sgN levels were similar across all variants.

CONCLUSIONS

RNA viral loads were similar among hospitalized adults, irrespective of infecting variant and known risk factors for severe COVID-19. Total N and subgenomic RNA N viral loads were highly correlated, suggesting that subgenomic RNA measurements add little information for the purposes of estimating infectivity.

Evolution of transient RNA structure-RNA polymerase interactions in respiratory RNA virus genomes

RNA viruses are important human pathogens that cause seasonal epidemics and occasional pandemics. Examples are influenza A viruses (IAV) and coronaviruses (CoV). When emerging IAV and CoV spill over to humans, they adapt to evade immune responses and optimize their replication and spread in human cells. In IAV, adaptation occurs in all viral proteins, including the viral ribonucleoprotein (RNP) complex. RNPs consists of a copy of the viral RNA polymerase, a double-helical coil of nucleoprotein, and one of the eight segments of the IAV RNA genome. The RNA segments and their transcripts are partially structured to coordinate the packaging of the viral genome and modulate viral mRNA translation. In addition, RNA structures can affect the efficiency of viral RNA synthesis and the activation of host innate immune response. Here, we investigated if RNA structures that modulate IAV replication processivity, so called template loops (t-loops), vary during the adaptation of pandemic and emerging IAV to humans. Using cell culture-based replication assays and in silico sequence analyses, we find that the sensitivity of the IAV H3N2 RNA polymerase to t-loops increased between isolates from 1968 and 2017, whereas the total free energy of t-loops in the IAV H3N2 genome was reduced. This reduction is particularly prominent in the PB1 gene. In H1N1 IAV, we find two separate reductions in t-loop free energy, one following the 1918 pandemic and one following the 2009 pandemic. No destabilization of t-loops is observed in the IBV genome, whereas analysis of SARS-CoV-2 isolates reveals destabilization of viral RNA structures. Overall, we propose that a loss of free energy in the RNA genome of emerging respiratory RNA viruses may contribute to the adaption of these viruses to the human population.

Loss-of-function mutation in Omicron variants reduces spike protein expression and attenuates SARS-CoV-2 infection

SARS-CoV-2 Omicron variants emerged in 2022 with >30 novel amino acid mutations in the spike protein alone. While most studies focus on receptor binding domain changes, mutations in the C-terminus of S1 (CTS1), adjacent to the furin cleavage site, have largely been ignored. In this study, we examined three Omicron mutations in CTS1: H655Y, N679K, and P681H. Generating a SARS-CoV-2 triple mutant (YKH), we found that the mutant increased spike processing, consistent with prior reports for H655Y and P681H individually. Next, we generated a single N679K mutant, finding reduced viral replication in vitro and less disease in vivo. Mechanistically, the N679K mutant had reduced spike protein in purified virions compared to wild-type; spike protein decreases were further exacerbated in infected cell lysates. Importantly, exogenous spike expression also revealed that N679K reduced overall spike protein yield independent of infection. Although a loss-of-function mutation, transmission competition demonstrated that N679K had a replication advantage in the upper airway over wild-type SARS-CoV-2 in hamsters, potentially impacting transmissibility. Together, the data show that N679K reduces overall spike protein levels during Omicron infection, which has important implications for infection, immunity, and transmission.

An improved method for determining frequency of multiple variants of SARS-CoV-2 in wastewater using qPCR assays

Wastewater-based surveillance has become an effective tool around the globe for indirect monitoring of COVID-19 in communities. Variants of Concern (VOCs) have been detected in wastewater by use of reverse transcription polymerase chain reaction (RT-PCR) or whole genome sequencing (WGS). Rapid, reliable RT-PCR assays continue to be needed to determine the relative frequencies of VOCs and sub-lineages in wastewater-based surveillance programs. The presence of multiple mutations in a single region of the N-gene allowed for the design of a single amplicon, multiple probe assay, that can distinguish among several VOCs in wastewater RNA extracts. This approach which multiplexes probes designed to target mutations associated with specific VOC's along with an intra-amplicon universal probe (non-mutated region) was validated in singleplex and multiplex. The prevalence of each mutation (i.e. VOC) is estimated by comparing the abundance of the targeted mutation with a non-mutated and highly conserved region within the same amplicon. This is advantageous for the accurate and rapid estimation of variant frequencies in wastewater. The N200 assay was applied to monitor frequencies of VOCs in wastewater extracts from several communities in Ontario, Canada in near real time from November 28, 2021 to January 4, 2022. This includes the period of the rapid replacement of the Delta variant with the introduction of the Omicron variant in these Ontario communities in early December 2021. The frequency estimates using this assay were highly reflective of clinical WGS estimates for the same communities. This style of qPCR assay, which simultaneously measures signal from a non-mutated comparator probe and multiple mutation-specific probes contained within a single qPCR amplicon, can be applied to future assay development for rapid and accurate estimations of variant frequencies.

Interaction between host G3BP and viral nucleocapsid protein regulates SARS-CoV-2 replication

G3BP1/2 are paralogous proteins that promote stress granule formation in response to cellular stresses, including viral infection. G3BP1/2 are prominent interactors of the nucleocapsid (N) protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, the functional consequences of the G3BP1-N interaction in the context of viral infection remain unclear. Here we used structural and biochemical analyses to define the residues required for G3BP1-N interaction, followed by structure-guided mutagenesis of G3BP1 and N to selectively and reciprocally disrupt their interaction. We found that mutation of F17 within the N protein led to selective loss of interaction with G3BP1 and consequent failure of the N protein to disrupt stress granule assembly. Introduction of SARS-CoV-2 bearing an F17A mutation resulted in a significant decrease in viral replication and pathogenesis in vivo, indicating that the G3BP1-N interaction promotes infection by suppressing the ability of G3BP1 to form stress granules.

SARS-CoV-2 Inhibits NRF2-Mediated Antioxidant Responses in Airway Epithelial Cells and in the Lung of a Murine Model of Infection

Several viruses have been shown to modulate the transcription factor nuclear factor erythroid 2-related factor 2 (NRF2), the master regulator of redox homeostasis. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for the COVID-19 pandemic, also seems to disrupt the balance between oxidants and antioxidants, which likely contributes to lung damage. Using in vitro and in vivo models of infection, we investigated how SARS-CoV-2 modulates the transcription factor NRF2 and its dependent genes, as well as the role of NRF2 during SARS-CoV-2 infection. We found that SARS-CoV-2 infection downregulates NRF2 protein levels and NRF2-dependent gene expression in human airway epithelial cells and in lungs of BALB/c mice. Reductions in cellular levels of NRF2 seem to be independent of proteasomal degradation and the interferon/promyelocytic leukemia (IFN/PML) pathway. Furthermore, lack of the Nrf2 gene in SARS-CoV-2-infected mice exacerbates clinical disease, increases lung inflammation, and is associated with a trend toward increased lung viral titers, indicating that NRF2 has a protective role during this viral infection. In summary, our results suggest that SARS-CoV-2 infection alters the cellular redox balance by downregulating NRF2 and its dependent genes, which exacerbates lung inflammation and disease, therefore, suggesting that the activation of NRF2 could be explored as therapeutic approach during SARS-CoV-2 infection. IMPORTANCE The antioxidant defense system plays a major function in protecting the organism against oxidative damage caused by free radicals. COVID-19 patients often present with biochemical characteristics of uncontrolled pro-oxidative responses in the respiratory tract. We show herein that SARS-CoV-2 variants, including Omicron, are potent inhibitors of cellular and lung nuclear factor erythroid 2-related factor 2 (NRF2), the master transcription factor that controls the expression of antioxidant and cytoprotective enzymes. Moreover, we show that mice lacking the Nrf2 gene show increased clinical signs of disease and lung pathology when infected with a mouse-adapted strain of SARS-CoV-2. Overall, this study provides a mechanistic explanation for the observed unbalanced pro-oxidative response in SARS-CoV-2 infections and suggests that therapeutic strategies for COVID-19 may consider the use of pharmacologic agents that are known to boost the expression levels of cellular NRF2.

Highly Networked SARS-CoV-2 Peptides Elicit T Cell Responses with Enhanced Specificity

Identifying SARS-CoV-2-specific T cell epitope-derived peptides is critical for the development of effective vaccines and measuring the duration of specific SARS-CoV-2 cellular immunity. In this regard, we previously identified T cell epitope-derived peptides within topologically and structurally essential regions of SARS-CoV-2 spike and nucleocapsid proteins by applying an immunoinformatics pipeline. In this study, we selected 30 spike- and nucleocapsid-derived peptides and assessed whether these peptides induce T cell responses and avoid major mutations found in SARS-CoV-2 variants of concern. Our peptide pool was highly specific, with only a single peptide driving cross-reactivity in people unexposed to SARS-COV-2, and immunogenic, inducing a polyfunctional response in CD4+ and CD8+ T cells from COVID-19 recovered individuals. All peptides were immunogenic and individuals recognized broad and diverse peptide repertoires. Moreover, our peptides avoided most mutations/deletions associated with all four SARS-CoV-2 variants of concern while retaining their physicochemical properties even when genetic changes are introduced. This study contributes to an evolving definition of individual CD4+ and CD8+ T cell epitopes that can be used for specific diagnostic tools for SARS-CoV-2 T cell responses and is relevant to the development of variant-resistant and durable T cell-stimulating vaccines.