A close shave: How SARS-CoV-2 induces the loss of cilia.

Wang et al. report in this issue (2022. J. Cell Biol.https://doi.org/10.1083/jcb.202108015) that the SARS-CoV-2 protein ORF10 increases the activity of the E3 ligase CUL2ZYG11B, leading to the degradation of multiple ciliary proteins. The resulting loss of cilia may facilitate the spread of SARS-CoV-2 in the respiratory tree.

Real-Time RT-PCR Allelic Discrimination Assay for Detection of N501Y Mutation in the Spike Protein of SARS-CoV-2 Associated with B.1.1.7 Variant of Concern.

The N501Y amino acid mutation caused by a single point substitution A23063T in the spike gene of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is possessed by three variants of concern (VOCs), B.1.1.7, B.1.351, and P.1. A rapid screening tool using this mutation is important for surveillance during the coronavirus disease 2019 (COVID-19) pandemic. We developed and validated a single nucleotide polymorphism real-time reverse transcription PCR assay using allelic discrimination of the spike gene N501Y mutation to screen for potential variants of concern and differentiate them from SARS-CoV-2 lineages without the N501Y mutation. A total of 160 clinical specimens positive for SARS-CoV-2 were characterized as mutant (N501Y) or N501 wild type by Sanger sequencing and were subsequently tested with the N501Y single nucleotide polymorphism real-time reverse transcriptase PCR assay. Our assay, compared to Sanger sequencing for single nucleotide polymorphism detection, demonstrated positive percent agreement of 100% for all 57 specimens displaying the N501Y mutation, which were confirmed by Sanger sequencing to be typed as A23063T, including one specimen with mixed signal for wild type and mutant. Negative percent agreement was 100% in all 103 specimens typed as N501 wild type, with A23063 identified as wild type by Sanger sequencing. The identification of circulating SARS-CoV-2 lineages carrying an N501Y mutation is critical for surveillance purposes. Current identification methods rely primarily on Sanger sequencing or whole-genome sequencing, which are time consuming, labor intensive, and costly. The assay described herein is an efficient tool for high-volume specimen screening for SARS-CoV-2 VOCs and for selecting specimens for confirmatory Sanger or whole-genome sequencing. IMPORTANCE During the coronavirus disease 2019 (COVID-19) pandemic, several variants of concern (VOCs) have been detected, for example, B.1.1.7, B.1.351, P.1, and B.1.617.2. The VOCs pose a threat to public health efforts to control the spread of the virus. As such, surveillance and monitoring of these VOCs is of the utmost importance. Our real-time RT-PCR assay helps with surveillance by providing an easy method to quickly survey SARS-CoV-2 specimens for VOCs carrying the N501Y single nucleotide polymorphism (SNP). Samples that test positive for the N501Y mutation in the spike gene with our assay can be sequenced to identify the lineage. Thus, our assay helps to focus surveillance efforts and decrease turnaround times.

Unmatched Virus Preservation Solution may Impede One-Step Reagent for Detection of 2019-nCoV RNA.

BACKGROUND AND METHODS: 2019 Corona Virus Disease (COVID-19) caused by SARS-CoV-2 is still pandemic now. RT-qPCR detection was the most common method for the diagnosis of SARS-CoV-2 infection, facilitated by amounts of nucleic acid testing kits. However, the accuracy of nucleic acid detection is affected by various factors such as specimen collection, specimen preparation, reagents deficiency, and personnel quality. RESULTS: In this study, we found that unmatched virus preservation solution will inhibit N gene and OFR-1ab gene (two independent genes of SARS-CoV-2) amplification in one-step detection reagent. CONCLUSIONS: Despite just being a particular phenomenon we found in our work to fight 2019-nCoV, we concluded that unmatched virus preservation solution may have an inhibitory effect on SARS-CoV-2 nucleic acid detection which may lead to incorrect clinical diagnosis.

Molecular features similarities between SARS-CoV-2, SARS, MERS and key human genes could favour the viral infections and trigger collateral effects.

In December 2019, rising pneumonia cases caused by a novel ß-coronavirus (SARS-CoV-2) occurred in Wuhan, China, which has rapidly spread worldwide, causing thousands of deaths. The WHO declared the SARS-CoV-2 outbreak as a public health emergency of international concern, since then several scientists are dedicated to its study. It has been observed that many human viruses have codon usage biases that match highly expressed proteins in the tissues they infect and depend on the host cell machinery for the replication and co-evolution. In this work, we analysed 91 molecular features and codon usage patterns for 339 viral genes and 463 human genes that consisted of 677,873 codon positions. Hereby, we selected the highly expressed genes from human lung tissue to perform computational studies that permit to compare their molecular features with those of SARS, SARS-CoV-2 and MERS genes. The integrated analysis of all the features revealed that certain viral genes and overexpressed human genes have similar codon usage patterns. The main pattern was the A/T bias that together with other features could propitiate the viral infection, enhanced by a host dependant specialization of the translation machinery of only some of the overexpressed genes. The envelope protein E, the membrane glycoprotein M and ORF7 could be further benefited. This could be the key for a facilitated translation and viral replication conducting to different comorbidities depending on the genetic variability of population due to the host translation machinery. This is the first codon usage approach that reveals which human genes could be potentially deregulated due to the codon usage similarities between the host and the viral genes when the virus is already inside the human cells of the lung tissues. Our work leaded to the identification of additional highly expressed human genes which are not the usual suspects but might play a role in the viral infection and settle the basis for further research in the field of human genetics associated with new viral infections. To identify the genes that could be deregulated under a viral infection is important to predict the collateral effects and determine which individuals would be more susceptible based on their genetic features and comorbidities associated.

Preliminary report of an outbreak of SARS-CoV-2 in mink and mink farmers associated with community spread, Denmark, June to November 2020.

In June-November 2020, SARS-CoV-2-infected mink were detected in 290 of 1,147 Danish mink farms. In North Denmark Region, 30% (324/1,092) of people found connected to mink farms tested SARS-CoV-2-PCR-positive and approximately 27% (95% confidence interval (CI): 25-30) of SARS-CoV-2-strains from humans in the community were mink-associated. Measures proved insufficient to mitigate spread. On 4 November, the government ordered culling of all Danish mink. Farmed mink constitute a potential virus reservoir challenging pandemic control.

Toolkit for Quickly Generating and Characterizing Molecular Probes Specific for SARS-CoV-2 Nucleocapsid as a Primer for Future Coronavirus Pandemic Preparedness.

Generating and characterizing immunoreagents to enable studies of novel emerging viruses is an area where ensembles of synthetic genes, recombinant antibody pipelines, and modular antibody-reporter fusion proteins can respond rapidly. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to spread through the global population causing widespread morbidity, mortality, and socioeconomic chaos. Using SARS-CoV-2 as our model and starting with a gBlocks encoded nucleocapsid (N) gene, we purified recombinant protein from E. coli, to serve as bait for selecting semisynthetic nanobodies from our Nomad single-pot library. Clones were isolated in days and first fused to Gaussia luciferase to determine EC50 in the tens of nM range, and second fused to the ascorbate peroxidase derivative APEX2 for sensitive detection of SARS-CoV-2 infected cells. To generate inherently fluorescent immunoreagents, we introduce novel periplasmic sdAb fusions made with mNeonGreen and mScarlet-I, which were produced at milligram amounts. The fluorescent fusion proteins enabled concise visualization of SARS-CoV-2 N in the cytoplasm but not in the nucleus 24 h post infection, akin to the distribution of SARS-CoV N, thereby validating these useful imaging tools. SdAb reactivity appeared specific to SARS-CoV-2 with very much weaker binding to SARS-CoV, and no noticeable cross-reactivity to a panel of overexpressed human codon optimized N proteins from other CoV. High periplasmic expression levels and in silico immortalization of the nanobody constructs guarantees a cost-effective and reliable source of SARS-CoV-2 immunoreagents. Our proof-of-principle study should be applicable to known and newly emerging CoV to broaden the tools available for their analysis and help safeguard human health in a more proactive than reactive manner.

Comparative genetic analyses of Korean bat coronaviruses with SARS-CoV and the newly emerged SARS-CoV-2.

BACKGROUND: Bats have been considered natural reservoirs for several pathogenic human coronaviruses (CoVs) in the last two decades. Recently, a bat CoV was detected in the Republic of Korea; its entire genome was sequenced and reported to be genetically similar to that of the severe acute respiratory syndrome CoV (SARS-CoV). OBJECTIVES: The objective of this study was to compare the genetic sequences of SARS-CoV, SARS-CoV-2, and the two Korean bat CoV strains 16BO133 and B15-21, to estimate the likelihood of an interaction between the Korean bat CoVs and the human angiotensin-converting enzyme 2 (ACE2) receptor. METHODS: The phylogenetic analysis was conducted with the maximum-likelihood (ML) method using MEGA 7 software. The Korean bat CoVs receptor binding domain (RBD) of the spike protein was analyzed by comparative homology modeling using the SWISS-MODEL server. The binding energies of the complexes were calculated using PRODIGY and MM/GBGA. RESULTS: Phylogenetic analyses of the entire RNA-dependent RNA polymerase, spike regions, and the complete genome revealed that the Korean CoVs, along with SARS-CoV and SARS-CoV-2, belong to the subgenus Sarbecovirus, within BetaCoVs. However, the two Korean CoVs were distinct from SARS-CoV-2. Specifically, the spike gene of the Korean CoVs, which is involved in host infection, differed from that of SARS-CoV-2, showing only 66.8%-67.0% nucleotide homology and presented deletions within the RBD, particularly within regions critical for cross-species transmission and that mediate interaction with ACE2. Binding free energy calculation revealed that the binding affinity of Korean bat CoV RBD to hACE2 was drastically lower than that of SARS-CoV and SARS-CoV-2. CONCLUSIONS: These results suggest that Korean bat CoVs are unlikely to bind to the human ACE2 receptor.

Rapid Spread of Mutant Alleles in Worldwide SARS-CoV-2 Strains Revealed by Genome-Wide Single Nucleotide Polymorphism and Variation Analysis.

The novel coronavirus (SARS-CoV-2) has become a pandemic and is threatening human health globally. Here, we report nine newly evolved SARS-CoV-2 single nucleotide polymorphism (SNP) alleles those underwent a rapid increase (seven cases) or decrease (two cases) in their frequency for 30-80% in the initial four months, which are further confirmed by intrahost single nucleotide variation analysis using raw sequence data including 8,217 samples. The nine SNPs are mostly (8/9) located in the coding region and are mainly (6/9) nonsynonymous substitutions. The nine SNPs show a complete linkage in SNP pairs and belong to three different linkage groups, named LG_1 to LG_3. Analyses in population genetics show signatures of adaptive selection toward the mutants in LG_1, but no signal of selection for LG_2. Population genetic analysis results on LG_3 show geological differentiation. Analyses on geographic COVID-19 cases and published clinical data provide evidence that the mutants in LG_1 and LG_3 benefit virus replication and those in LG_1 have a positive correlation with the disease severity in COVID-19-infected patients. The mutants in LG_2 show a bias toward mildness of the disease based on available public clinical data. Our findings may be instructive for epidemiological surveys and disease control of COVID-19 in the future.

Overcoming limitations in the availability of swabs systems used for SARS-CoV-2 laboratory diagnostics.

The diagnosis of COVID-19 relies on the direct detection of SARS-CoV-2 RNA in respiratory specimens by RT-PCR. The pandemic spread of the disease caused an imbalance between demand and supply of materials and reagents needed for diagnostic purposes including swab sets. In a comparative effectiveness study, we conducted serial follow-up swabs in hospitalized laboratory-confirmed COVID-19 patients. We assessed the diagnostic performance of an in-house system developed according to recommendations by the US CDC. In a total of 96 serial swabs, we found significant differences in the accuracy of the different swab systems to generate a positive result in SARS-CoV-2 RT-PCR, ranging from around 50 to 80%. Of note, an in-house swab system was superior to most commercially available sets as reflected by significantly lower Ct values of viral genes. Thus, a simple combination of broadly available materials may enable diagnostic laboratories to bypass global limitations in the supply of swab sets.

Network machine learning maps phytochemically rich "Hyperfoods" to fight COVID-19.

In this paper, we introduce a network machine learning method to identify potential bioactive anti-COVID-19 molecules in foods based on their capacity to target the SARS-CoV-2-host gene-gene (protein-protein) interactome. Our analyses were performed using a supercomputing DreamLab App platform, harnessing the idle computational power of thousands of smartphones. Machine learning models were initially calibrated by demonstrating that the proposed method can predict anti-COVID-19 candidates among experimental and clinically approved drugs (5658 in total) targeting COVID-19 interactomics with the balanced classification accuracy of 80-85% in 5-fold cross-validated settings. This identified the most promising drug candidates that can be potentially "repurposed" against COVID-19 including common drugs used to combat cardiovascular and metabolic disorders, such as simvastatin, atorvastatin and metformin. A database of 7694 bioactive food-based molecules was run through the calibrated machine learning algorithm, which identified 52 biologically active molecules, from varied chemical classes, including flavonoids, terpenoids, coumarins and indoles predicted to target SARS-CoV-2-host interactome networks. This in turn was used to construct a "food map" with the theoretical anti-COVID-19 potential of each ingredient estimated based on the diversity and relative levels of candidate compounds with antiviral properties. We expect this in silico predicted food map to play an important role in future clinical studies of precision nutrition interventions against COVID-19 and other viral diseases.

Spatial and temporal roles of SARS-CoV PLpro -A snapshot.

SARS-CoV and SARS-CoV-2 encode four structural and accessory proteins (spike, envelope, membrane and nucleocapsid proteins) and two polyproteins (pp1a and pp1ab). The polyproteins are further cleaved by 3C-like cysteine protease (3CLpro ) and papain-like protease (PLpro ) into 16 nonstructural proteins (nsps). PLpro is released from nsp3 through autocleavage, and then it cleaves the sites between nsp1/2, between nsp2/3 and between nsp3/4 with recognition motif of LXGG, and the sites in the C-terminus of ubiquitin and of protein interferon-stimulated gene 15 (ISG15) with recognition motif of RLRGG. Alone or together with SARS unique domain (SUD), PLpro can stabilize an E3 ubiquitin ligase, the ring-finger, and CHY zinc-finger domain-containing 1 (RCHY1), through domain interaction, and thus, promote RCHY1 to ubiquitinate its target proteins including p53. However, a dilemma appears in terms of PLpro roles. On the one hand, the ubiquitination of p53 is good for SARS-CoV because the ubiquitinated p53 cannot inhibit SARS-CoV replication. On the other hand, the ubiquitination of NF-κB inhibitor (IκBα), TNF receptor-associated factors (TRAFs), and stimulator of interferon gene (STING), and the ISGylation of targeted proteins are bad for SARS-CoV because these ubiquitination and ISGylation initiate the innate immune response and antiviral state. This mini-review analyzes the dilemma and provides a snapshot on how the viral PLpro smartly manages its roles to avoid its simultaneously contradictory actions, which could shed lights on possible strategies to deal with SARS-CoV-2 infections.

Spatial and temporal roles of SARS-CoV PLpro -A snapshot.

SARS-CoV and SARS-CoV-2 encode four structural and accessory proteins (spike, envelope, membrane and nucleocapsid proteins) and two polyproteins (pp1a and pp1ab). The polyproteins are further cleaved by 3C-like cysteine protease (3CLpro ) and papain-like protease (PLpro ) into 16 nonstructural proteins (nsps). PLpro is released from nsp3 through autocleavage, and then it cleaves the sites between nsp1/2, between nsp2/3 and between nsp3/4 with recognition motif of LXGG, and the sites in the C-terminus of ubiquitin and of protein interferon-stimulated gene 15 (ISG15) with recognition motif of RLRGG. Alone or together with SARS unique domain (SUD), PLpro can stabilize an E3 ubiquitin ligase, the ring-finger, and CHY zinc-finger domain-containing 1 (RCHY1), through domain interaction, and thus, promote RCHY1 to ubiquitinate its target proteins including p53. However, a dilemma appears in terms of PLpro roles. On the one hand, the ubiquitination of p53 is good for SARS-CoV because the ubiquitinated p53 cannot inhibit SARS-CoV replication. On the other hand, the ubiquitination of NF-κB inhibitor (IκBα), TNF receptor-associated factors (TRAFs), and stimulator of interferon gene (STING), and the ISGylation of targeted proteins are bad for SARS-CoV because these ubiquitination and ISGylation initiate the innate immune response and antiviral state. This mini-review analyzes the dilemma and provides a snapshot on how the viral PLpro smartly manages its roles to avoid its simultaneously contradictory actions, which could shed lights on possible strategies to deal with SARS-CoV-2 infections.

A review on viral data sources and search systems for perspective mitigation of COVID-19.

With the outbreak of the COVID-19 disease, the research community is producing unprecedented efforts dedicated to better understand and mitigate the effects of the pandemic. In this context, we review the data integration efforts required for accessing and searching genome sequences and metadata of SARS-CoV2, the virus responsible for the COVID-19 disease, which have been deposited into the most important repositories of viral sequences. Organizations that were already present in the virus domain are now dedicating special interest to the emergence of COVID-19 pandemics, by emphasizing specific SARS-CoV2 data and services. At the same time, novel organizations and resources were born in this critical period to serve specifically the purposes of COVID-19 mitigation while setting the research ground for contrasting possible future pandemics. Accessibility and integration of viral sequence data, possibly in conjunction with the human host genotype and clinical data, are paramount to better understand the COVID-19 disease and mitigate its effects. Few examples of host-pathogen integrated datasets exist so far, but we expect them to grow together with the knowledge of COVID-19 disease; once such datasets will be available, useful integrative surveillance mechanisms can be put in place by observing how common variants distribute in time and space, relating them to the phenotypic impact evidenced in the literature.

CoV2-ID, a MIQE-compliant sub-20-min 5-plex RT-PCR assay targeting SARS-CoV-2 for the diagnosis of COVID-19.

Accurate, reliable and rapid detection of SARS-CoV-2 is essential not only for correct diagnosis of individual COVID-19 disease but also for the development of a rational strategy aimed at lifting confinement restrictions and preparing for possible recurrent waves of viral infections. We have used the MIQE guidelines to develop two versions of a unique five plex RT-qPCR test, termed CoV2-ID, that allows the detection of three viral target genes, a human internal control for confirming the presence of human cells in a sample and a control artificial RNA for quality assessment and potential quantification. Viral targets can be detected either individually with separate fluorophores or jointly using the same fluorophore, thus increasing the test's reliability and sensitivity. It is robust, can consistently detect two copies of viral RNA, with a limit of detection of a single copy and can be completed in around 15 min. It was 100% sensitive and 100% specific when tested on 23 RNA samples extracted from COVID-19 positive patients and five COVID-19 negative patients. We also propose using multiple cycle fluorescence detection, rather than real-time PCR to reduce significantly the time taken to complete the assay as well as assuage the misunderstandings underlying the use of quantification cycles (Cq). Finally, we have designed an assay for the detection of the D614G mutation and show that all of the samples isolated in the Chelmsford, Essex area between mid-April and June 2020, have the mutant genotype whereas a sample originating in Australia was infected with the wild type genotype.

Simultaneous Dual-Gene Diagnosis of SARS-CoV-2 Based on CRISPR/Cas9-Mediated Lateral Flow Assay.

Few methods for the detection of SARS-CoV-2 currently have the capability to simultaneously detect two genes in a single test, which is a key measure to improve detection accuracy, as adopted by the gold standard RT-qPCR method. Developed here is a CRISPR/Cas9-mediated triple-line lateral flow assay (TL-LFA) combined with multiplex reverse transcription-recombinase polymerase amplification (RT-RPA) for rapid and simultaneous dual-gene detection of SARS-CoV-2 in a single strip test. This assay is characterized by the detection of envelope (E) and open reading frame 1ab (Orf1ab) genes from cell-cultured SARS-CoV-2 and SARS-CoV-2 viral RNA standards, showing a sensitivity of 100 RNA copies per reaction (25 µL). Furthermore, dual-gene analysis of 64 nasopharyngeal swab samples showed 100 % negative predictive agreement and 97.14 % positive predictive agreement. This platform will provide a more accurate and convenient pathway for diagnosis of COVID-19 or other infectious diseases in low-resource regions.

An Integrated Multi-Function Heterogeneous Biochemical Circuit for High-Resolution Electrochemistry-Based Genetic Analysis.

Modular construction of an autonomous and programmable multi-functional heterogeneous biochemical circuit that can identify, transform, translate, and amplify biological signals into physicochemical signals based on logic design principles can be a powerful means for the development of a variety of biotechnologies. To explore the conceptual validity, we design a CRISPR-array-mediated primer-exchange-reaction-based biochemical circuit cascade, which probes a specific biomolecular input, transform the input into a structurally accessible form for circuit wiring, translate the input information into an arbitrary sequence, and finally amplify the prescribed sequence through autonomous formation of a signaling concatemer. This upstream biochemical circuit is further wired with a downstream electrochemical interface, delivering an integrated bioanalytical platform. We program this platform to directly analyze the genome of SARS-CoV-2 in human cell lysate, demonstrating the capability and the utility of this unique integrated system.

Global analysis of more than 50,000 SARS-CoV-2 genomes reveals epistasis between eight viral genes.

Genome-wide epistasis analysis is a powerful tool to infer gene interactions, which can guide drug and vaccine development and lead to deeper understanding of microbial pathogenesis. We have considered all complete severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes deposited in the Global Initiative on Sharing All Influenza Data (GISAID) repository until four different cutoff dates, and used direct coupling analysis together with an assumption of quasi-linkage equilibrium to infer epistatic contributions to fitness from polymorphic loci. We find eight interactions, of which three are between pairs where one locus lies in gene ORF3a, both loci holding nonsynonymous mutations. We also find interactions between two loci in gene nsp13, both holding nonsynonymous mutations, and four interactions involving one locus holding a synonymous mutation. Altogether, we infer interactions between loci in viral genes ORF3a and nsp2, nsp12, and nsp6, between ORF8 and nsp4, and between loci in genes nsp2, nsp13, and nsp14. The paper opens the prospect to use prominent epistatically linked pairs as a starting point to search for combinatorial weaknesses of recombinant viral pathogens.

Compositional diversity and evolutionary pattern of coronavirus accessory proteins.

Accessory proteins play important roles in the interaction between coronaviruses and their hosts. Accordingly, a comprehensive study of the compositional diversity and evolutionary patterns of accessory proteins is critical to understanding the host adaptation and epidemic variation of coronaviruses. Here, we developed a standardized genome annotation tool for coronavirus (CoroAnnoter) by combining open reading frame prediction, transcription regulatory sequence recognition and homologous alignment. Using CoroAnnoter, we annotated 39 representative coronavirus strains to form a compositional profile for all of the accessary proteins. Large variations were observed in the number of accessory proteins of 1-10 for different coronaviruses, with SARS-CoV-2 and SARS-CoV having the most (9 and 10, respectively). The variation between SARS-CoV and SARS-CoV-2 accessory proteins could be traced back to related coronaviruses in other hosts. The genomic distribution of accessory proteins had significant intra-genus conservation and inter-genus diversity and could be grouped into 1, 4, 2 and 1 types for alpha-, beta-, gamma-, and delta-coronaviruses, respectively. Evolutionary analysis suggested that accessory proteins are more conservative locating before the N-terminal of proteins E and M (E-M), while they are more diverse after these proteins. Furthermore, comparison of virus-host interaction networks of SARS-CoV-2 and SARS-CoV accessory proteins showed that they share multiple antiviral signaling pathways, those involved in the apoptotic process, viral life cycle and response to oxidative stress. In summary, our study provides a tool for coronavirus genome annotation and builds a comprehensive profile for coronavirus accessory proteins covering their composition, classification, evolutionary pattern and host interaction.

Evolutionary dynamics of SARS-CoV-2 nucleocapsid protein and its consequences.

The emerged novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has created a global health crisis that warrants an accurate and detailed characterization of the rapidly evolving viral genome for understanding its epidemiology, pathogenesis, and containment. Here, we explored 61,485 sequences of the nucleocapsid (N) protein, a potent diagnostic and prophylactic target, for identifying the mutations to review their roles in real-time polymerase chain reaction based diagnosis and observe consequent impacts. Compared to the Wuhan reference strain, a total of 1034 unique nucleotide mutations were identified in the mutant strains (49.15%, n = 30,221) globally. Of these mutations, 367 occupy primer binding sites including the 3'-end mismatch to the primer-pair of 11 well-characterized primer sets. Noteworthily, CDC (USA) recommended the N2 primer set contained a lower mismatch than the other primer sets. Moreover, 684 amino acid (aa) substitutions were located across 317 (75.66% of total aa) unique positions including 82, 21, and 83 of those in the RNA binding N-terminal domain (NTD), SR-rich region, and C-terminal dimerization domain, respectively. Moreover, 11 in-frame deletions, mostly (n = 10) within the highly flexible linker region, were revealed, and the rest was within the NTD region. Furthermore, we predicted the possible consequence of high-frequency mutations (≥20) and deletions on the tertiary structure of the N protein. Remarkably, we observed that a high frequency (67.94% of mutated sequences) co-occuring mutations (R203K and G204R) destabilized and decreased overall structural flexibility. The N protein of SARS-CoV-2 comprises an average of 1.2 mutations per strain compared to 4.4 and 0.4 in Middle East respiratory syndrome-related coronavirus and SARS-CoV, respectively. Despite being proposed as the alternative target to spike protein for vaccine and therapeutics, the ongoing evolution of the N protein may challenge these endeavors, thus needing further immunoinformatics analyses. Therefore, continuous monitoring is required for tracing the ongoing evolution of the SARS-CoV-2 N protein in prophylactic and diagnostic interventions.

Dynamically evolving novel overlapping gene as a factor in the SARS-CoV-2 pandemic.

Understanding the emergence of novel viruses requires an accurate and comprehensive annotation of their genomes. Overlapping genes (OLGs) are common in viruses and have been associated with pandemics but are still widely overlooked. We identify and characterize ORF3d, a novel OLG in SARS-CoV-2 that is also present in Guangxi pangolin-CoVs but not other closely related pangolin-CoVs or bat-CoVs. We then document evidence of ORF3d translation, characterize its protein sequence, and conduct an evolutionary analysis at three levels: between taxa (21 members of Severe acute respiratory syndrome-related coronavirus), between human hosts (3978 SARS-CoV-2 consensus sequences), and within human hosts (401 deeply sequenced SARS-CoV-2 samples). ORF3d has been independently identified and shown to elicit a strong antibody response in COVID-19 patients. However, it has been misclassified as the unrelated gene ORF3b, leading to confusion. Our results liken ORF3d to other accessory genes in emerging viruses and highlight the importance of OLGs.