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

Geospatial and demographic patterns of SARS-CoV-2 spread in Massachusetts from over 130,000 genomes

Despite intensive study, gaps remain in our understanding of SARS-CoV-2 transmission patterns during the COVID-19 pandemic, in part due to limited contextual metadata accompanying most large genomic surveillance datasets. We analyzed over 130,000 SARS-CoV-2 genomes, over 85,000 with matched epidemiological data, collected in Massachusetts from November 2021 to January 2023, to investigate viral transmission dynamics at high resolution. The data were drawn from diagnostic testing at >600 facilities representing schools, workplaces, public testing, and other sectors, and encompass the emergence of six major viral lineages, each representing a new outbreak. We found urban areas as key hubs for new lineage introduction and interurban transmission as facilitating spread throughout the state. Young adults, especially those on college campuses, served as early indicators of emerging lineage dominance. Resident-aged populations in college campuses and nursing homes exhibited a higher likelihood of being linked to within-facility transmission, while staff-aged at those facilities were more linked to their surrounding community. Individuals with recent vaccine doses, including boosters, had a lower likelihood of initiating transmission. This dataset shows the value of linking genomic and epidemiologic data at scale for higher resolution insights into viral dynamics and their implication for public health strategy.

In silico analysis of balsaminol as anti-viral agents targeting SARS-CoV-2 main protease, spike receptor binding domain and papain-like protease receptors

Plant-derived phytochemicals from medicinal plants are becoming increasingly attractive natural sources of antimicrobial and antiviral agents due to their therapeutic value, mechanism of action, level of toxicity and bioavailability. The continued emergence of more immune-evasive strains and the rate of resistance to current antiviral drugs have created a need to identify new antiviral agents against SARS-CoV-2. This study investigated the antiviral potential of balsaminol, a bioactive compound from Momordica balsamina, and its inhibitory activities against SARS-CoV-2 receptor proteins. In this study, three Food and Drug Administration (FDA) COVID-19 approved drugs namely; nirmatrelvir, ritonavir and remdesivir were used as positive control. Molecular docking was performed to determine the predominant binding mode (most negative Gibbs free energy of binding/ΔG) and inhibitory activity of balsaminol against SARS-CoV-2 receptor proteins. The pharmacokinetics, toxicity, physicochemical and drug-like properties of balsaminol were evaluated to determine its potential as an active oral drug candidate as well as its non-toxicity in humans. The results show that balsaminol E has the highest binding affinity to the SARS CoV-2 papain-like protease (7CMD) with a free binding energy of - 8.7 kcal/mol, followed by balsaminol A interacting with the spike receptor binding domain (6VW1) with - 8.5 kcal/mol and balsaminol C had a binding energy of - 8.1 kcal/mol with the main protease (6LU7) comparable to the standard drugs namely ritonavir, nirmatrelvir and remdesivir. However, the ADMET and drug-like profile of balsaminol F favours it as a better potential drug candidate and inhibitor of the docked SARS-CoV-2 receptor proteins. Further preclinical studies are therefore recommended.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-024-00241-0.

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

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

The Chameleon Strategy-A Recipe for Effective Ligand Screening for Viral Targets Based on Four Novel Structure-Binding Strength Indices

The RNA viruses SARS-CoV, SARS-CoV-2 and MERS-CoV encode the non-structural Nsp16 (2'-O-methyltransferase) that catalyzes the transfer of a methyl group from S-adenosylmethionine (SAM) to the first ribonucleotide in mRNA. Recently, it has been found that breaking the bond between Nsp16 and SAM substrate results in the cessation of mRNA virus replication. To date, only a limited number of such inhibitors have been identified, which can be attributed to a lack of an effective "recipe". The aim of our study was to propose and verify a rapid and effective screening protocol dedicated to such purposes. We proposed four new indices describing structure-binding strength (structure-binding affinity, structure-hydrogen bonding, structure-steric and structure-protein-ligand indices) were then applied and shown to be extremely helpful in determining the degree of increase or decrease in binding affinity in response to a relatively small change in the ligand structure. After initial pre-selection, based on similarity to SAM, we limited the study to 967 compounds, so-called molecular chameleons. They were then docked in the Nsp16 protein pocket, and 10 candidate ligands were selected using the novel structure-binding affinity index. Subsequently the selected 10 candidate ligands and 8 known inhibitors and were docked to Nsp16 pockets from SARS-CoV-2, MERS-CoV and SARS-CoV. Based on the four new indices, the best ligands were selected and a new one was designed by tuning them. Finally, ADMET profiling and molecular dynamics simulations were performed for the best ligands. The new structure-binding strength indices can be successfully applied not only to screen and tune ligands, but also to determine the effectiveness of the ligand in response to changes in the target viral entity, which is particularly useful for assessing drug effectiveness in the case of alterations in viral proteins. The developed approach, the so-called chameleon strategy, has the capacity to introduce a novel universal paradigm to the field of drugs design, including RNA antivirals.

Mapping of long stretches of highly conserved sequences in over 6 million SARS-CoV-2 genomes

We identified 11 conserved stretches in over 6.3 million SARS-CoV-2 genomes including all the major variants of concerns. Each conserved stretch is ≥100 nucleotides in length with ≥99.9% conservation at each nucleotide position. Interestingly, six of the eight conserved stretches in ORF1ab overlapped significantly with well-folded experimentally verified RNA secondary structures. Furthermore, two of the conserved stretches were mapped to regions within the S2-subunit that undergo dynamic structural rearrangements during viral fusion. In addition, the conserved stretches were significantly depleted for zinc-finger antiviral protein (ZAP) binding sites, which facilitated the recognition and degradation of viral RNA. These highly conserved stretches in the SARS-CoV-2 genome were poorly conserved at the nucleotide level among closely related β-coronaviruses, thus representing ideal targets for highly specific and discriminatory diagnostic assays. Our findings highlight the role of structural constraints at both RNA and protein levels that contribute to the sequence conservation of specific genomic regions in SARS-CoV-2.

Genes involved in the limited spread of SARS-CoV-2 in the lower respiratory airways of hamsters may be associated with adaptive evolution

Novel respiratory viruses can cause a pandemic and then evolve to coexist with humans. The Omicron strain of severe acute respiratory syndrome coronavirus 2 has spread worldwide since its emergence in late 2021, and its sub-lineages are now established in human society. Compared to previous strains, Omicron is markedly less invasive in the lungs and causes less severe disease. One reason for this is that humans are acquiring immunity through previous infection and vaccination, but the nature of the virus itself is also changing. Using our newly established low-volume inoculation system, which reflects natural human infection, we show that the Omicron strain spreads less efficiently into the lungs of hamsters compared with an earlier Wuhan strain. Furthermore, by characterizing chimeric viruses with the Omicron gene in the Wuhan strain genetic background and vice versa, we found that viral genes downstream of ORF3a, but not the S gene, were responsible for the limited spread of the Omicron strain in the lower airways of the virus-infected hamsters. Moreover, molecular evolutionary analysis of SARS-CoV-2 revealed a positive selection of genes downstream of ORF3a (M and E genes). Our findings provide insight into the adaptive evolution of the virus in humans during the pandemic convergence phase.IMPORTANCEThe severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant has spread worldwide since its emergence in late 2021, and its sub-lineages are established in human society. Compared to previous strains, the Omicron strain is less invasive in the lower respiratory tract, including the lungs, and causes less severe disease; however, the mechanistic basis for its restricted replication in the lower airways is poorly understood. In this study, using a newly established low-volume inoculation system that reflects natural human infection, we demonstrated that the Omicron strain spreads less efficiently into the lungs of hamsters compared with an earlier Wuhan strain and found that viral genes downstream of ORF3a are responsible for replication restriction in the lower respiratory tract of Omicron-infected hamsters. Furthermore, we detected a positive selection of genes downstream of ORF3a (especially the M and E genes) in SARS-CoV-2, suggesting that these genes may undergo adaptive changes in humans.

What makes SARS-CoV-2 unique? Focusing on the spike protein

Severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) seriously threatens public health and safety. Genetic variants determine the expression of SARS-CoV-2 structural proteins, which are associated with enhanced transmissibility, enhanced virulence, and immune escape. Vaccination is encouraged as a public health intervention, and different types of vaccines are used worldwide. However, new variants continue to emerge, especially the Omicron complex, and the neutralizing antibody responses are diminished significantly. In this review, we outlined the uniqueness of SARS-CoV-2 from three perspectives. First, we described the detailed structure of the spike (S) protein, which is highly susceptible to mutations and contributes to the distinct infection cycle of the virus. Second, we systematically summarized the immunoglobulin G epitopes of SARS-CoV-2 and highlighted the central role of the nonconserved regions of the S protein in adaptive immune escape. Third, we provided an overview of the vaccines targeting the S protein and discussed the impact of the nonconserved regions on vaccine effectiveness. The characterization and identification of the structure and genomic organization of SARS-CoV-2 will help elucidate its mechanisms of viral mutation and infection and provide a basis for the selection of optimal treatments. The leaps in advancements regarding improved diagnosis, targeted vaccines and therapeutic remedies provide sound evidence showing that scientific understanding, research, and technology evolved at the pace of the pandemic.

Improved detection of low-frequency within-host variants from deep sequencing: A case study with human papillomavirus

High-coverage sequencing allows the study of variants occurring at low frequencies within samples, but is susceptible to false-positives caused by sequencing error. Ion Torrent has a very low single nucleotide variant (SNV) error rate and has been employed for the majority of human papillomavirus (HPV) whole genome sequences. However, benchmarking of intrahost SNVs (iSNVs) has been challenging, partly due to limitations imposed by the HPV life cycle. We address this problem by deep sequencing three replicates for each of 31 samples of HPV type 18 (HPV18). Errors, defined as iSNVs observed in only one of three replicates, are dominated by C→T (G→A) changes, independently of trinucleotide context. True iSNVs, defined as those observed in all three replicates, instead show a more diverse SNV type distribution, with particularly elevated C→T rates in CCG context (CCG→CTG; CGG→CAG) and C→A rates in ACG context (ACG→AAG; CGT→CTT). Characterization of true iSNVs allowed us to develop two methods for detecting true variants: (1) VCFgenie, a dynamic binomial filtering tool which uses each variant's allele count and coverage instead of fixed frequency cut-offs; and (2) a machine learning binary classifier which trains eXtreme Gradient Boosting models on variant features such as quality and trinucleotide context. Each approach outperforms fixed-cut-off filtering of iSNVs, and performance is enhanced when both are used together. Our results provide improved methods for identifying true iSNVs in within-host applications across sequencing platforms, specifically using HPV18 as a case study.

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

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

The SARS-CoV-2 protein ORF3c is a mitochondrial modulator of innate immunity

The SARS-CoV-2 genome encodes a multitude of accessory proteins. Using comparative genomic approaches, an additional accessory protein, ORF3c, has been predicted to be encoded within the ORF3a sgmRNA. Expression of ORF3c during infection has been confirmed independently by ribosome profiling. Despite ORF3c also being present in the 2002-2003 SARS-CoV, its function has remained unexplored. Here we show that ORF3c localizes to mitochondria, where it inhibits innate immunity by restricting IFN-β production, but not NF-κB activation or JAK-STAT signaling downstream of type I IFN stimulation. We find that ORF3c is inhibitory after stimulation with cytoplasmic RNA helicases RIG-I or MDA5 or adaptor protein MAVS, but not after TRIF, TBK1 or phospho-IRF3 stimulation. ORF3c co-immunoprecipitates with the antiviral proteins MAVS and PGAM5 and induces MAVS cleavage by caspase-3. Together, these data provide insight into an uncharacterized mechanism of innate immune evasion by this important human pathogen.

Alternative Reading Frames are an Underappreciated Source of Protein Sequence Novelty

Protein-coding DNA sequences can be translated into completely different amino acid sequences if the nucleotide triplets used are shifted by a non-triplet amount on the same DNA strand or by translating codons from the opposite strand. Such "alternative reading frames" of protein-coding genes are a major contributor to the evolution of novel protein products. Recent studies demonstrating this include examples across the three domains of cellular life and in viruses. These sequences increase the number of trials potentially available for the evolutionary invention of new genes and also have unusual properties which may facilitate gene origin. There is evidence that the structure of the standard genetic code contributes to the features and gene-likeness of some alternative frame sequences. These findings have important implications across diverse areas of molecular biology, including for genome annotation, structural biology, and evolutionary genomics.

From Alpha to Omicron: How Different Variants of Concern of the SARS-Coronavirus-2 Impacted the World

SARS-CoV-2, the virus that causes COVID-19, is prone to mutations and the generation of genetic variants. Since its first outbreak in 2019, SARS-CoV-2 has continually evolved, resulting in the emergence of several lineages and variants of concern (VOC) that have gained more efficient transmission, severity, and immune evasion properties. The World Health Organization has given these variants names according to the letters of the Greek Alphabet, starting with the Alpha (B.1.1.7) variant, which emerged in 2020, followed by the Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2), and Omicron (B.1.1.529) variants. This review explores the genetic variation among different VOCs of SARS-CoV-2 and how the emergence of variants made a global impact on the pandemic.

Selection for immune evasion in SARS-CoV-2 revealed by high-resolution epitope mapping and sequence analysis

Here, we exploit a deep serological profiling strategy coupled with an integrated, computational framework for the analysis of SARS-CoV-2 humoral immune responses. Applying a high-density peptide array (HDPA) spanning the entire proteomes of SARS-CoV-2 and endemic human coronaviruses allowed identification of B cell epitopes and relate them to their evolutionary and structural properties. We identify hotspots of pre-existing immunity and identify cross-reactive epitopes that contribute to increasing the overall humoral immune response to SARS-CoV-2. Using a public dataset of over 38,000 viral genomes from the early phase of the pandemic, capturing both inter- and within-host genetic viral diversity, we determined the evolutionary profile of epitopes and the differences across proteins, waves, and SARS-CoV-2 variants. Lastly, we show that mutations in spike and nucleocapsid epitopes are under stronger selection between than within patients, suggesting that most of the selective pressure for immune evasion occurs upon transmission between hosts.

Timeline of SARS-CoV-2 Transmission in Sabah, Malaysia: Tracking the Molecular Evolution

BACKGROUND

The COVID-19 pandemic poses an unprecedented public health challenge in Malaysia. The impact of COVID-19 varies between countries, including geographically divided states within a country. The deadly transmission of COVID-19 has taken a heavy toll in Sabah, Malaysia's third most populous state, contributing nearly 10% to the recorded national death toll as of 31 December 2022. Although several SARS-CoV-2 genome sequences have been analysed in Malaysia, molecular epidemiology data from Sabah focusing on the diversity and evolution of SARS-CoV-2 variants are still lacking. This study examines the major SARS-CoV-2 variants and emerging mutations from Sabah, the Malaysian Borneo, which is geographically divided from West Malaysia by the South China Sea.

METHODS

A total of 583 COVID-19 samples were subjected to whole genome sequencing and analysed with an additional 1123 Sabah COVID-19 sequences retrieved from the GISAID EpiCoV consortium. Nextclade and Pangolin were used to classify these sequences according to the clades and lineages. To determine the molecular evolutionary characteristics, Bayesian time-scaled phylogenetic analysis employing the maximum likelihood algorithm was performed on selected SARS-CoV-2 genome sequences, using the Wuhan-Hu-1 sequence as a reference.

RESULTS

Sabah was affected starting from the second COVID-19 wave in Malaysia, and the early sequences were classified under the O clade. The clade was gradually replaced during subsequent waves by G, GH, GK and GRA, with the latter being dominant as of December 2022. Phylogenetically, the Delta isolates in this study belong to the three main subclades 21A, 21J and 21I, while Omicron isolates belong to 21M, 21L and 22B. The time-scaled phylogeny suggested that SARS-CoV-2 introduced into Sabah originated from Peninsular Malaysia in early March 2020, and phylodynamic analysis indicated that increased viral spread was observed in early March and declined in late April, followed by an evolutionary stationary phase in June 2020.

CONCLUSION

Continuous molecular epidemiology of SARS-CoV-2 in Sabah will provide a deeper understanding of the emergence and dominance of each variant in the locality, thus facilitating public health intervention measures.

Natural selection shapes the evolution of SARS-CoV-2 Omicron in Bangladesh

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved to give rise to a highly transmissive and immune-escaping variant of concern, known as Omicron. Many aspects of the evolution of SARS-CoV-2 and the driving forces behind the ongoing Omicron outbreaks remain unclear. Substitution at the receptor-binding domain (RBD) in the spike protein is one of the primary strategies of SARS-CoV-2 Omicron to hinder recognition by the host angiotensin-converting enzyme 2 (ACE2) receptor and avoid antibody-dependent defense activation. Here, we scanned for adaptive evolution within the SARS-CoV-2 Omicron genomes reported from Bangladesh in the public database GISAID (www.gisaid.org; dated 2 April 2023). The ratio of the non-synonymous (Ka) to synonymous (Ks) nucleotide substitution rate, denoted as ω, is an indicator of the selection pressure acting on protein-coding genes. A higher proportion of non-synonymous to synonymous substitutions (Ka/Ks or ω > 1) indicates positive selection, while Ka/Ks or ω near zero indicates purifying selection. An equal amount of non-synonymous and synonymous substitutions (Ka/Ks or ω = 1) refers to neutrally evolving sites. We found evidence of adaptive evolution within the spike (S) gene of SARS-CoV-2 Omicron isolated from Bangladesh. In total, 22 codon sites of the S gene displayed a signature of positive selection. The data also highlighted that the receptor-binding motif within the RBD of the spike glycoprotein is a hotspot of adaptive evolution, where many of the codons had ω > 1. Some of these adaptive sites at the RBD of the spike protein are known to be associated with increased viral fitness. The M gene and ORF6 have also experienced positive selection. These results suggest that although purifying selection is the dominant evolutionary force, positive Darwinian selection also plays a vital role in shaping the evolution of SARS-CoV-2 Omicron in Bangladesh.

Plasticity, Paralogy, and Pseudogenization: Rhabdoviruses of Freshwater Mussels Elucidate Mechanisms of Viral Genome Diversification and the Evolution of the Finfish-Infecting Rhabdoviral Genera

Viruses in the family Rhabdoviridae display remarkable genomic variation and ecological diversity. This plasticity occurs despite the fact that, as negative sense RNA viruses, rhabdoviruses rarely if ever recombine. Here, we describe nonrecombinatorial evolutionary processes leading to genomic diversification in the Rhabdoviridae inferred from two novel rhabdoviruses of freshwater mussels (Mollusca: Bivalvia: Unionida). Killamcar virus 1 (KILLV-1) from a plain pocketbook (Lampsilis cardium) is closely related phylogenetically and transcriptionally to finfish-infecting viruses in the subfamily Alpharhabdovirinae. KILLV-1 offers a novel example of glycoprotein gene duplication, differing from previous examples in that the paralogs overlap. Evolutionary analyses reveal a clear pattern of relaxed selection due to subfunctionalization in rhabdoviral glycoprotein paralogs, which has not previously been described in RNA viruses. Chemarfal virus 1 (CHMFV-1) from a western pearlshell (Margaritifera falcata) is closely related phylogenetically and transcriptionally to viruses in the genus Novirhabdovirus, the sole recognized genus in the subfamily Gammarhabdovirinae, representing the first known gammarhabdovirus of a host other than finfish. The CHMFV-1 G-L noncoding region contains a nontranscribed remnant gene of precisely the same length as the NV gene of most novirhabdoviruses, offering a compelling example of pseudogenization. The unique reproductive strategy of freshwater mussels involves an obligate parasitic stage in which larvae encyst in the tissues of finfish, offering a plausible ecological mechanism for viral host-switching. IMPORTANCE Viruses in the family Rhabdoviridae infect a variety of hosts, including vertebrates, invertebrates, plants and fungi, with important consequences for health and agriculture. This study describes two newly discovered viruses of freshwater mussels from the United States. One virus from a plain pocketbook (Lampsilis cardium) is closely related to fish-infecting viruses in the subfamily Alpharhabdovirinae. The other virus from a western pearlshell (Margaritifera falcata) is closely related to viruses in the subfamily Gammarhabdovirinae, which until now were only known to infect finfish. Genome features of both viruses provide new evidence of how rhabdoviruses evolved their extraordinary variability. Freshwater mussel larvae attach to fish and feed on tissues and blood, which may explain how rhabdoviruses originally jumped between mussels and fish. The significance of this research is that it improves our understanding of rhabdovirus ecology and evolution, shedding new light on these important viruses and the diseases they cause.

SARS-CoV-2 Evasion of the Interferon System: Can We Restore Its Effectiveness?

Type I and III Interferons (IFNs) are the first lines of defense in microbial infections. They critically block early animal virus infection, replication, spread, and tropism to promote the adaptive immune response. Type I IFNs induce a systemic response that impacts nearly every cell in the host, while type III IFNs' susceptibility is restricted to anatomic barriers and selected immune cells. Both IFN types are critical cytokines for the antiviral response against epithelium-tropic viruses being effectors of innate immunity and regulators of the development of the adaptive immune response. Indeed, the innate antiviral immune response is essential to limit virus replication at the early stages of infection, thus reducing viral spread and pathogenesis. However, many animal viruses have evolved strategies to evade the antiviral immune response. The Coronaviridae are viruses with the largest genome among the RNA viruses. Severe Acute Respiratory Syndrome-Coronavirus-2 (SARS-CoV-2) caused the coronavirus disease 2019 (COVID-19) pandemic. The virus has evolved numerous strategies to contrast the IFN system immunity. We intend to describe the virus-mediated evasion of the IFN responses by going through the main phases: First, the molecular mechanisms involved; second, the role of the genetic background of IFN production during SARS-CoV-2 infection; and third, the potential novel approaches to contrast viral pathogenesis by restoring endogenous type I and III IFNs production and sensitivity at the sites of infection.

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

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

Translation landscape of SARS-CoV-2 noncanonical subgenomic RNAs

The ongoing COVID-19 pandemic is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with a positive-stranded RNA genome. Current proteomic studies of SARS-CoV-2 mainly focus on the proteins encoded by its genomic RNA (gRNA) or canonical subgenomic RNAs (sgRNAs). Here, we systematically investigated the translation landscape of SARS-CoV-2, especially its noncanonical sgRNAs. We first constructed a strict pipeline, named vipep, for identifying reliable peptides derived from RNA viruses using RNA-seq and mass spectrometry data. We applied vipep to analyze 24 sets of mass spectrometry data related to SARS-CoV-2 infection. In addition to known canonical proteins, we identified many noncanonical sgRNA-derived peptides, which stably increase after viral infection. Furthermore, we explored the potential functions of those proteins encoded by noncanonical sgRNAs and found that they can bind to viral RNAs and may have immunogenic activity. The generalized vipep pipeline is applicable to any RNA viruses and these results have expanded the SARS-CoV-2 translation map, providing new insights for understanding the functions of SARS-CoV-2 sgRNAs.

Druggable targets and therapeutic development for COVID-19

Coronavirus disease (COVID-19), which is caused by SARS-CoV-2, is the biggest challenge to the global public health and economy in recent years. Until now, only limited therapeutic regimens have been available for COVID-19 patients, sparking unprecedented efforts to study coronavirus biology. The genome of SARS-CoV-2 encodes 16 non-structural, four structural, and nine accessory proteins, which mediate the viral life cycle, including viral entry, RNA replication and transcription, virion assembly and release. These processes depend on the interactions between viral polypeptides and host proteins, both of which could be potential therapeutic targets for COVID-19. Here, we will discuss the potential medicinal value of essential proteins of SARS-CoV-2 and key host factors. We summarize the most updated therapeutic interventions for COVID-19 patients, including those approved clinically or in clinical trials.

SARS-CoV-2 ORF6 disrupts innate immune signalling by inhibiting cellular mRNA export

SARS-CoV-2 is a betacoronavirus and the etiological agent of COVID-19, a devastating infectious disease. Due to its far-reaching effect on human health, there is an urgent and growing need to understand the viral molecular biology of SARS-CoV-2 and its interaction with the host cell. SARS-CoV-2 encodes 9 predicted accessory proteins, which are presumed to be dispensable for in vitro replication, most likely having a role in modulating the host cell environment to aid viral replication. Here we show that the ORF6 accessory protein interacts with cellular Rae1 to inhibit cellular protein production by blocking mRNA export. We utilised cell fractionation coupled with mRNAseq to explore which cellular mRNA species are affected by ORF6 expression and show that ORF6 can inhibit the export of many mRNA including those encoding antiviral factors such as IRF1 and RIG-I. We also show that export of these mRNA is blocked in the context of SARS-CoV-2 infection. Together, our studies identify a novel mechanism by which SARS-CoV-2 can manipulate the host cell environment to supress antiviral responses, providing further understanding to the replication strategies of a virus that has caused an unprecedented global health crisis.

Sequence similarity of SARS-CoV-2 and humans: Implications for SARS-CoV-2 detection

Detecting severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) needs human samples, which inevitably contain trace human DNA and RNA. Sequence similarity may cause invalid detection results; however, there is still a lack of gene similarity analysis of SARS-CoV-2 and humans. All publicly reported complete genome assemblies in the Entrez genome database were collected for multiple sequence alignment, similarity and phylogenetic analysis. The complete genomes showed high similarity (>99.88% sequence identity). Phylogenetic analysis divided these viruses into three major clades with significant geographic group effects. Viruses from the United States showed considerable variability. Sequence similarity analysis revealed that SARS-CoV-2 has 612 similar sequences with the human genome and 100 similar sequences with the human transcriptome. The sequence characteristics and genome distribution of these similar sequences were confirmed. The sequence similarity and evolutionary mutations provide indispensable references for dynamic updates of SARS-CoV-2 detection primers and methods.

Multiple Novel Mosquito-Borne Zoonotic Viruses Revealed in Pangolin Virome

Swab samples were collected from 34 pangolins in Guangxi Province, China. Metavirome sequencing and bioinformatics approaches were undertaken to determine the abundant viral sequences in the viromes. The results showed that the viral sequences belong to 24 virus taxonomic families. To verify the results, PCR combined with phylogenetic analysis was conducted. Some viral sequences including Japanese encephalitis virus (JEV), Getah virus (GETV), and chikungunya virus (CHIKV) were detected. On the basis of the metavirome analysis, seven segments belonging to JEV were further identified through PCR amplification. Sequence comparison showed that, among seven sequences, JEV-China/P2020E-1 displayed the highest nucleotide (80.6%), with the JEV isolated in South Korea, 1988, and all of which belonging to genotype III. Seven CHIKV sequences were detected, with the highest homology (80.6%) to the Aedes africanus in Côte d’Ivoire, 1993. Moreover, passage from BHK-21 to Vero cells makes the newly isolated CHIKV-China/P2020-1 more contagious. In addition, the newly verified GETV sequences shared 86.4% identity with the 1955 GETV isolated from Malaysia. Some sudden and recurrent viruses have also been observed from the virome of pangolin in Guangxi Province, China; hence, dissemination tests will be implemented in the future.

Molecular Virology of SARS-CoV-2 and Related Coronaviruses

Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The global COVID-19 pandemic continues to threaten the lives of hundreds of millions of people, with a severe negative impact on the global economy. Although several COVID-19 vaccines are currently being administered, none of them is 100% effective. Moreover, SARS-CoV-2 variants remain an important worldwide public health issue. Hence, the accelerated development of efficacious antiviral agents is urgently needed. Coronavirus depends on various host cell factors for replication. An ongoing research objective is the identification of host factors that could be exploited as targets for drugs and compounds effective against SARS-CoV-2. In the present review, we discuss the molecular mechanisms of SARS-CoV-2 and related coronaviruses, focusing on the host factors or pathways involved in SARS-CoV-2 replication that have been identified by genome-wide CRISPR screening.

Systems biology of malaria explored with nonhuman primates

"The Primate Malarias" book has been a uniquely important resource for multiple generations of scientists, since its debut in 1971, and remains pertinent to the present day. Indeed, nonhuman primates (NHPs) have been instrumental for major breakthroughs in basic and pre-clinical research on malaria for over 50 years. Research involving NHPs have provided critical insights and data that have been essential for malaria research on many parasite species, drugs, vaccines, pathogenesis, and transmission, leading to improved clinical care and advancing research goals for malaria control, elimination, and eradication. Whilst most malaria scientists over the decades have been studying Plasmodium falciparum, with NHP infections, in clinical studies with humans, or using in vitro culture or rodent model systems, others have been dedicated to advancing research on Plasmodium vivax, as well as on phylogenetically related simian species, including Plasmodium cynomolgi, Plasmodium coatneyi, and Plasmodium knowlesi. In-depth study of these four phylogenetically related species over the years has spawned the design of NHP longitudinal infection strategies for gathering information about ongoing infections, which can be related to human infections. These Plasmodium-NHP infection model systems are reviewed here, with emphasis on modern systems biological approaches to studying longitudinal infections, pathogenesis, immunity, and vaccines. Recent discoveries capitalizing on NHP longitudinal infections include an advanced understanding of chronic infections, relapses, anaemia, and immune memory. With quickly emerging new technological advances, more in-depth research and mechanistic discoveries can be anticipated on these and additional critical topics, including hypnozoite biology, antigenic variation, gametocyte transmission, bone marrow dysfunction, and loss of uninfected RBCs. New strategies and insights published by the Malaria Host-Pathogen Interaction Center (MaHPIC) are recapped here along with a vision that stresses the importance of educating future experts well trained in utilizing NHP infection model systems for the pursuit of innovative, effective interventions against malaria.

Homology-based classification of accessory proteins in coronavirus genomes uncovers extremely dynamic evolution of gene content

Coronaviruses (CoVs) have complex genomes that encode a fixed array of structural and nonstructural components, as well as a variety of accessory proteins that differ even among closely related viruses. Accessory proteins often play a role in the suppression of immune responses and may represent virulence factors. Despite their relevance for CoV phenotypic variability, information on accessory proteins is fragmentary. We applied a systematic approach based on homology detection to create a comprehensive catalogue of accessory proteins encoded by CoVs. Our analyses grouped accessory proteins into 379 orthogroups and 12 super‐groups. No orthogroup was shared by the four CoV genera and very few were present in all or most viruses in the same genus, reflecting the dynamic evolution of CoV genomes. We observed differences in the distribution of accessory proteins in CoV genera. Alphacoronaviruses harboured the largest diversity of accessory open reading frames (ORFs), deltacoronaviruses the smallest. However, the average number of accessory proteins per genome was highest in betacoronaviruses. Analysis of the evolutionary history of some orthogroups indicated that the different CoV genera adopted similar evolutionary strategies. Thus, alphacoronaviruses and betacoronaviruses acquired phosphodiesterases and spike‐like accessory proteins independently, whereas horizontal gene transfer from reoviruses endowed betacoronaviruses and deltacoronaviruses with fusion‐associated small transmembrane (FAST) proteins. Finally, analysis of accessory ORFs in annotated CoV genomes indicated ambiguity in their naming. This complicates cross‐communication among researchers and hinders automated searches of large data sets (e.g., PubMed, GenBank). We suggest that orthogroup membership is used together with a naming system to provide information on protein function.

Functional and structural segregation of overlapping helices in HIV-1

Overlapping coding regions balance selective forces between multiple genes. One possible division of nucleotide sequence is that the predominant selective force on a particular nucleotide can be attributed to just one gene. While this arrangement has been observed in regions in which one gene is structured and the other is disordered, we sought to explore how overlapping genes balance constraints when both protein products are structured over the same sequence. We use a combination of sequence analysis, functional assays, and selection experiments to examine an overlapped region in HIV-1 that encodes helical regions in both Env and Rev. We find that functional segregation occurs even in this overlap, with each protein spacing its functional residues in a manner that allows a mutable non-binding face of one helix to encode important functional residues on a charged face in the other helix. Additionally, our experiments reveal novel and critical functional residues in Env and have implications for the therapeutic targeting of HIV-1.

Prevalence and circulation patterns of SARS-CoV-2 variants in European sewage mirror clinical data of 54 European cities

For community-level monitoring, the European Commission under the EU Sewage Sentinel System recommends wastewater-based SARS-CoV-2 surveillance. Tracking SARS-CoV-2 variants in a community is pivotal for appropriate public health response. Genome sequencing of SARS-CoV-2 in wastewater samples for tracking variants is challenging, often resulting in low coverage genome sequences, thereby impeding the detection of the SARS-CoV-2 mutations. Therefore, we aimed at high-coverage SARS-CoV-2 genome sequences from sewage samples which we successfully accomplished. This first pan-European surveillance compared the mutation profiles associated with the variants of concerns: B.1.1.7, P.1, B.1.351 and B.1.617.2 across 20 European countries, including 54 municipalities. The results highlight that SARS-CoV-2 variants detected in the wastewater samples mirror the variants profiles reported in clinical data. This study demonstrated that >98% coverage of SARS-CoV-2 genomic sequences is possible and can be used to track SARS-CoV-2 mutations in wastewater to support identifying variants circulating in a city at the community level.

SARS-CoV-2 Nsp14 mediates the effects of viral infection on the host cell transcriptome

Viral infection involves complex set of events orchestrated by multiple viral proteins. To identify functions of SARS-CoV-2 proteins, we performed transcriptomic analyses of cells expressing individual viral proteins. Expression of Nsp14, a protein involved in viral RNA replication, provoked a dramatic remodeling of the transcriptome that strongly resembled that observed following SARS-CoV-2 infection. Moreover, Nsp14 expression altered the splicing of more than 1000 genes and resulted in a dramatic increase in the number of circRNAs, which are linked to innate immunity. These effects were independent of the Nsp14 exonuclease activity and required the N7-guanine-methyltransferase domain of the protein. Activation of the NFkB pathway and increased expression of CXCL8 occurred early upon Nsp14 expression. We identified IMPDH2, which catalyzes the rate-limiting step of guanine nucleotides biosynthesis, as a key mediator of these effects. Nsp14 expression caused an increase in GTP cellular levels, and the effect of Nsp14 was strongly decreased in the presence of IMPDH2 inhibitors. Together, our data demonstrate an unknown role for Nsp14 with implications for therapy.

Viruses are parasites, relying on the cells they infect to make more of themselves. In doing so they change how an infected cell turns its genes on and off, forcing it to build new virus particles and turning off the immune surveillance that would allow the body to intervene. This is how SARS-CoV-2, the virus that causes COVID, survives with a genome that carries instructions to make just 29 proteins.

One of these proteins, known as Nsp14, is involved in both virus reproduction and immune escape. Previous work has shown that it interacts with IMPDH2, the cellular enzyme that controls the production of the building blocks of the genetic code. The impact of this interaction is not clear.

To find out more, Zaffagni et al. introduced 26 of the SARS-CoV-2 proteins into human cells one at a time. Nsp14 had the most dramatic effect, dialing around 4,000 genes up or down and changing how the cell interprets over 1,000 genes. Despite being just one protein, it mimicked the genetic changes seen during real SARS-CoV-2 infection. Blocking IMPDH2 partially reversed the effects, which suggests that the interaction of Nsp14 with the enzyme might be responsible for the effects of SARS-CoV-2 on the genes of the cell.

Understanding how viral proteins affect cells can explain what happens during infection. This could lead to the discovery of new treatments designed to counteract the effects of the virus. Further work could investigate whether interfering with Nsp14 helps cells to overcome infection.

A BioID-Derived Proximity Interactome for SARS-CoV-2 Proteins

The novel coronavirus SARS-CoV-2 is responsible for the ongoing COVID-19 pandemic and has caused a major health and economic burden worldwide. Understanding how SARS-CoV-2 viral proteins behave in host cells can reveal underlying mechanisms of pathogenesis and assist in development of antiviral therapies. Here, the cellular impact of expressing SARS-CoV-2 viral proteins was studied by global proteomic analysis, and proximity biotinylation (BioID) was used to map the SARS-CoV-2 virus–host interactome in human lung cancer-derived cells. Functional enrichment analyses revealed previously reported and unreported cellular pathways that are associated with SARS-CoV-2 proteins. We have established a website to host the proteomic data to allow for public access and continued analysis of host–viral protein associations and whole-cell proteomes of cells expressing the viral–BioID fusion proteins. Furthermore, we identified 66 high-confidence interactions by comparing this study with previous reports, providing a strong foundation for future follow-up studies. Finally, we cross-referenced candidate interactors with the CLUE drug library to identify potential therapeutics for drug-repurposing efforts. Collectively, these studies provide a valuable resource to uncover novel SARS-CoV-2 biology and inform development of antivirals.

Unheeded SARS-CoV-2 proteins? A deep look into negative-sense RNA

SARS-CoV-2 is a novel positive-sense single-stranded RNA virus from the Coronaviridae family (genus Betacoronavirus), which has been established as causing the COVID-19 pandemic. The genome of SARS-CoV-2 is one of the largest among known RNA viruses, comprising of at least 26 known protein-coding loci. Studies thus far have outlined the coding capacity of the positive-sense strand of the SARS-CoV-2 genome, which can be used directly for protein translation. However, it has been recently shown that transcribed negative-sense viral RNA intermediates that arise during viral genome replication from positive-sense viruses can also code for proteins. No studies have yet explored the potential for negative-sense SARS-CoV-2 RNA intermediates to contain protein-coding loci. Thus, using sequence and structure-based bioinformatics methodologies, we have investigated the presence and validity of putative negative-sense ORFs (nsORFs) in the SARS-CoV-2 genome. Nine nsORFs were discovered to contain strong eukaryotic translation initiation signals and high codon adaptability scores, and several of the nsORFs were predicted to interact with RNA-binding proteins. Evolutionary conservation analyses indicated that some of the nsORFs are deeply conserved among related coronaviruses. Three-dimensional protein modeling revealed the presence of higher order folding among all putative SARS-CoV-2 nsORFs, and subsequent structural mimicry analyses suggest similarity of the nsORFs to DNA/RNA-binding proteins and proteins involved in immune signaling pathways. Altogether, these results suggest the potential existence of still undescribed SARS-CoV-2 proteins, which may play an important role in the viral lifecycle and COVID-19 pathogenesis.

Spotlight on alternative frame coding: Two long overlapping genes in Pseudomonas aeruginosa are translated and under purifying selection

The existence of overlapping genes (OLGs) with significant coding overlaps revolutionizes our understanding of genomic complexity. We report two exceptionally long (957 nt and 1536 nt), evolutionarily novel, translated antisense open reading frames (ORFs) embedded within annotated genes in the pathogenic Gram-negative bacterium Pseudomonas aeruginosa. Both OLG pairs show sequence features consistent with being genes and transcriptional signals in RNA sequencing. Translation of both OLGs was confirmed by ribosome profiling and mass spectrometry. Quantitative proteomics of samples taken during different phases of growth revealed regulation of protein abundances, implying biological functionality. Both OLGs are taxonomically restricted, and likely arose by overprinting within the genus. Evidence for purifying selection further supports functionality. The OLGs reported here, designated olg1 and olg2, are the longest yet proposed in prokaryotes and are among the best attested in terms of translation and evolutionary constraint. These results highlight a potentially large unexplored dimension of prokaryotic genomes.

Are There Hidden Genes in DNA/RNA Vaccines?

Due to the fast global spreading of the Severe Acute Respiratory Syndrome Coronavirus - 2 (SARS-CoV-2), prevention and treatment options are direly needed in order to control infection-related morbidity, mortality, and economic losses. Although drug and inactivated and attenuated virus vaccine development can require significant amounts of time and resources, DNA and RNA vaccines offer a quick, simple, and cheap treatment alternative, even when produced on a large scale. The spike protein, which has been shown as the most antigenic SARS-CoV-2 protein, has been widely selected as the target of choice for DNA/RNA vaccines. Vaccination campaigns have reported high vaccination rates and protection, but numerous unintended effects, ranging from muscle pain to death, have led to concerns about the safety of RNA/DNA vaccines. In parallel to these studies, several open reading frames (ORFs) have been found to be overlapping SARS-CoV-2 accessory genes, two of which, ORF2b and ORF-Sh, overlap the spike protein sequence. Thus, the presence of these, and potentially other ORFs on SARS-CoV-2 DNA/RNA vaccines, could lead to the translation of undesired proteins during vaccination. Herein, we discuss the translation of overlapping genes in connection with DNA/RNA vaccines. Two mRNA vaccine spike protein sequences, which have been made publicly-available, were compared to the wild-type sequence in order to uncover possible differences in putative overlapping ORFs. Notably, the Moderna mRNA-1273 vaccine sequence is predicted to contain no frameshifted ORFs on the positive sense strand, which highlights the utility of codon optimization in DNA/RNA vaccine design to remove undesired overlapping ORFs. Since little information is available on ORF2b or ORF-Sh, we use structural bioinformatics techniques to investigate the structure-function relationship of these proteins. The presence of putative ORFs on DNA/RNA vaccine candidates implies that overlapping genes may contribute to the translation of smaller peptides, potentially leading to unintended clinical outcomes, and that the protein-coding potential of DNA/RNA vaccines should be rigorously examined prior to administration.

Gene Overlapping as a Modulator of Begomovirus Evolution

In RNA viruses, which have high mutation—and fast evolutionary— rates, gene overlapping (i.e., genomic regions that encode more than one protein) is a major factor controlling mutational load and therefore the virus evolvability. Although DNA viruses use host high-fidelity polymerases for their replication, and therefore should have lower mutation rates, it has been shown that some of them have evolutionary rates comparable to those of RNA viruses. Notably, these viruses have large proportions of their genes with at least one overlapping instance. Hence, gene overlapping could be a modulator of virus evolution beyond the RNA world. To test this hypothesis, we use the genus Begomovirus of plant viruses as a model. Through comparative genomic approaches, we show that terminal gene overlapping decreases the rate of virus evolution, which is associated with lower frequency of both synonymous and nonsynonymous mutations. In contrast, terminal overlapping has little effect on the pace of virus evolution. Overall, our analyses support a role for gene overlapping in the evolution of begomoviruses and provide novel information on the factors that shape their genetic diversity.

Overlapping genes in natural and engineered genomes

Modern genome-scale methods that identify new genes, such as proteogenomics and ribosome profiling, have revealed, to the surprise of many, that overlap in genes, open reading frames and even coding sequences is widespread and functionally integrated into prokaryotic, eukaryotic and viral genomes. In parallel, the constraints that overlapping regions place on genome sequences and their evolution can be harnessed in bioengineering to build more robust synthetic strains and constructs. With a focus on overlapping protein-coding and RNA-coding genes, this Review examines their discovery, topology and biogenesis in the context of their genome biology. We highlight exciting new uses for sequence overlap to control translation, compress synthetic genetic constructs, and protect against mutation.

Genome Evolution from Random Ligation of RNAs of Autocatalytic Sets

The evolutionary origin of the genome remains elusive. Here, I hypothesize that its first iteration, the protogenome, was a multi-ribozyme RNA. It evolved, likely within liposomes (the protocells) forming in dry-wet cycling environments, through the random fusion of ribozymes by a ligase and was amplified by a polymerase. The protogenome thereby linked, in one molecule, the information required to seed the protometabolism (a combination of RNA-based autocatalytic sets) in newly forming protocells. If this combination of autocatalytic sets was evolutionarily advantageous, the protogenome would have amplified in a population of multiplying protocells. It likely was a quasispecies with redundant information, e.g., multiple copies of one ribozyme. As such, new functionalities could evolve, including a genetic code. Once one or more components of the protometabolism were templated by the protogenome (e.g., when a ribozyme was replaced by a protein enzyme), and/or addiction modules evolved, the protometabolism became dependent on the protogenome. Along with increasing fidelity of the RNA polymerase, the protogenome could grow, e.g., by incorporating additional ribozyme domains. Finally, the protogenome could have evolved into a DNA genome with increased stability and storage capacity. I will provide suggestions for experiments to test some aspects of this hypothesis, such as evaluating the ability of ribozyme RNA polymerases to generate random ligation products and testing the catalytic activity of linked ribozyme domains.

Biological factors in the synthetic construction of overlapping genes

BACKGROUND

Overlapping genes (OLGs) with long protein-coding overlapping sequences are disallowed by standard genome annotation programs, outside of viruses. Recently however they have been discovered in Archaea, diverse Bacteria, and Mammals. The biological factors underlying life's ability to create overlapping genes require more study, and may have important applications in understanding evolution and in biotechnology. A previous study claimed that protein domains from viruses were much better suited to forming overlaps than those from other cellular organisms - in this study we assessed this claim, in order to discover what might underlie taxonomic differences in the creation of gene overlaps.

RESULTS

After overlapping arbitrary Pfam domain pairs and evaluating them with Hidden Markov Models we find OLG construction to be much less constrained than expected. For instance, close to 10% of the constructed sequences cannot be distinguished from typical sequences in their protein family. Most are also indistinguishable from natural protein sequences regarding identity and secondary structure. Surprisingly, contrary to a previous study, virus domains were much less suitable for designing OLGs than bacterial or eukaryotic domains were. In general, the amount of amino acid change required to force a domain to overlap is approximately equal to the variation observed within a typical domain family. The resulting high similarity between natural sequences and those altered so as to overlap is mostly due to the combination of high redundancy in the genetic code and the evolutionary exchangeability of many amino acids.

CONCLUSIONS

Synthetic overlapping genes which closely resemble natural gene sequences, as measured by HMM profiles, are remarkably easy to construct, and most arbitrary domain pairs can be altered so as to overlap while retaining high similarity to the original sequences. Future work however will need to assess important factors not considered such as intragenic interactions which affect protein folding. While the analysis here is not sufficient to guarantee functional folding proteins, further analysis of constructed OLGs will improve our understanding of the origin of these remarkable genetic elements across life and opens up exciting possibilities for synthetic biology.

Computational methods for inferring location and genealogy of overlapping genes in virus genomes: approaches and applications

Viruses may evolve to increase the amount of encoded genetic information by means of overlapping genes, which utilize several reading frames. Such overlapping genes may be especially impactful for genomes of small size, often serving a source of novel accessory proteins, some of which play a crucial role in viral pathogenicity or in promoting the systemic spread of virus. Diverse genome-based metrics were proposed to facilitate recognition of overlapping genes that otherwise may be overlooked during genome annotation. They can detect the atypical codon bias associated with the overlap (e.g. a statistically significant reduction in variability at synonymous sites) or other sequence-composition features peculiar to overlapping genes. In this review, I compare nine computational methods, discuss their strengths and limitations, and survey how they were applied to detect candidate overlapping genes in the genome of SARS-CoV-2, the etiological agent of COVID-19 pandemic.

ADAR Editing in Viruses: An Evolutionary Force to Reckon with

Adenosine Deaminases that Act on RNA (ADARs) are RNA editing enzymes that play a dynamic and nuanced role in regulating transcriptome and proteome diversity. This editing can be highly selective, affecting a specific site within a transcript, or nonselective, resulting in hyperediting. ADAR editing is important for regulating neural functions and autoimmunity, and has a key role in the innate immune response to viral infections, where editing can have a range of pro- or antiviral effects and can contribute to viral evolution. Here we examine the role of ADAR editing across a broad range of viral groups. We propose that the effect of ADAR editing on viral replication, whether pro- or antiviral, is better viewed as an axis rather than a binary, and that the specific position of a given virus on this axis is highly dependent on virus- and host-specific factors, and can change over the course of infection. However, more research needs to be devoted to understanding these dynamic factors and how they affect virus-ADAR interactions and viral evolution. Another area that warrants significant attention is the effect of virus-ADAR interactions on host-ADAR interactions, particularly in light of the crucial role of ADAR in regulating neural functions. Answering these questions will be essential to developing our understanding of the relationship between ADAR editing and viral infection. In turn, this will further our understanding of the effects of viruses such as SARS-CoV-2, as well as many others, and thereby influence our approach to treating these deadly diseases.

SARS-CoV-2 transgressing LncRNAs uncovers the known unknowns

SARS-CoV-2 harbors many known unknown regions in the form of hypothetical open reading frames (ORFs). Although the mechanisms underlying the disease pathogenesis are not clearly understood, molecules such as long noncoding RNAs (lncRNAs) play a key regulatory role in the viral pathogenesis from endocytosis. We asked whether or not the lncRNAs in the host are associated with the viral proteins and argue that lncRNA-mRNAs molecules related to viral infection may regulate SARS-CoV-2 pathogenesis. Toward the end of the perspective, we provide challenges and insights into investigating these transgression pathways.

Prediction of two novel overlapping ORFs in the genome of SARS-CoV-2

Six candidate overlapping genes have been detected in SARS-CoV-2, yet current methods struggle to detect overlapping genes that recently originated. However, such genes might encode proteins beneficial to the virus, and provide a model system to understand gene birth. To complement existing detection methods, I first demonstrated that selection pressure to avoid stop codons in alternative reading frames is a driving force in the origin and retention of overlapping genes. I then built a detection method, CodScr, based on this selection pressure. Finally, I combined CodScr with methods that detect other properties of overlapping genes, such as a biased nucleotide and amino acid composition. I detected two novel ORFs (ORF-Sh and ORF-Mh), overlapping the spike and membrane genes respectively, which are under selection pressure and may be beneficial to SARS-CoV-2. ORF-Sh and ORF-Mh are present, as ORF uninterrupted by stop codons, in 100% and 95% of the SARS-CoV-2 genomes, respectively.

Simple Classification of RNA Sequences of Respiratory-Related Coronaviruses

A very simple, fast, and efficient approach to analyze and identify respiratory-related virus sequences based on machine learning is proposed. Such schemes are very important in identifying viruses, especially in view of spreading pandemics. The method is based on genetic code rules and the open reading frame (ORF). Data from the respiratory-related coronaviruses are collected and features are extracted based on reoccurring nucleobase 3-tuples in the RNA. Our methodology is simply based on counting nucleobase triplets, normalizing the count to the length of the sequence, and applying principal component analysis (PCA) techniques. The triplet counting can be further used for classification purposes. DNA sequences from the herpes virus family can be considered as the first step towards a complete and accurate classification including more complex factors, such as mutations. The proposed classification scheme is simply based on "counting" biological information. It can serve as the first fast detection method, widely accessible and portable to a variety of distinct architectures for fast and on-the-fly detection. We provide an approach that can be further optimized and combined with supervised techniques to allow for more accurate detection and read out of the exact virus type or sequence. We discuss the relevance of this scheme in identifying differences in similar viruses and their impact on biochemical analysis.

Science's Response to CoVID-19

CoVID-19 is a multi-symptomatic disease which has made a global impact due to its ability to spread rapidly, and its relatively high mortality rate. Beyond the heroic efforts to develop vaccines, which we do not discuss herein, the response of scientists and clinicians to this complex problem has reflected the need to detect CoVID-19 rapidly, to diagnose patients likely to show adverse symptoms, and to treat severe and critical CoVID-19. Here we aim to encapsulate these varied and sometimes conflicting approaches and the resulting data in terms of chemistry and biology. In the process we highlight emerging concepts, and potential future applications that may arise out of this immense effort.

Coronavirus, the King Who Wanted More Than a Crown: From Common to the Highly Pathogenic SARS-CoV-2, Is the Key in the Accessory Genes?

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), that emerged in late 2019, is the etiologic agent of the current "coronavirus disease 2019" (COVID-19) pandemic, which has serious health implications and a significant global economic impact. Of the seven human coronaviruses, all of which have a zoonotic origin, the pandemic SARS-CoV-2, is the third emerging coronavirus, in the 21st century, highly pathogenic to the human population. Previous human coronavirus outbreaks (SARS-CoV-1 and MERS-CoV) have already provided several valuable information on some of the common molecular and cellular mechanisms of coronavirus infections as well as their origin. However, to meet the new challenge caused by the SARS-CoV-2, a detailed understanding of the biological specificities, as well as knowledge of the origin are crucial to provide information on viral pathogenicity, transmission and epidemiology, and to enable strategies for therapeutic interventions and drug discovery. Therefore, in this review, we summarize the current advances in SARS-CoV-2 knowledges, in light of pre-existing information of other recently emerging coronaviruses. We depict the specificity of the immune response of wild bats and discuss current knowledge of the genetic diversity of bat-hosted coronaviruses that promotes viral genome expansion (accessory gene acquisition). In addition, we describe the basic virology of coronaviruses with a special focus SARS-CoV-2. Finally, we highlight, in detail, the current knowledge of genes and accessory proteins which we postulate to be the major keys to promote virus adaptation to specific hosts (bat and human), to contribute to the suppression of immune responses, as well as to pathogenicity.

SARS-CoV-2 Accessory Proteins in Viral Pathogenesis: Knowns and Unknowns

There are still many unanswered questions concerning viral SARS-CoV-2 pathogenesis in COVID-19. Accessory proteins in SARS-CoV-2 consist of eleven viral proteins whose roles during infection are still not completely understood. Here, a review on the current knowledge of SARS-CoV-2 accessory proteins is summarized updating new research that could be critical in understanding SARS-CoV-2 interaction with the host. Some accessory proteins such as ORF3b, ORF6, ORF7a and ORF8 have been shown to be important IFN-I antagonists inducing an impairment in the host immune response. In addition, ORF3a is involved in apoptosis whereas others like ORF9b and ORF9c interact with cellular organelles leading to suppression of the antiviral response in infected cells. However, possible roles of ORF7b and ORF10 are still awaiting to be described. Also, ORF3d has been reassigned. Relevant information on the knowns and the unknowns in these proteins is analyzed, which could be crucial for further understanding of SARS-CoV-2 pathogenesis and to design strategies counteracting their actions evading immune responses in COVID-19.

Genetic Diversity and Evolution of the Biological Features of the Pandemic SARS-CoV-2

The new coronavirus infection (COVID-19) represents a challenge for global health. Since the outbreak began, the number of confirmed cases has exceeded 117 million, with more than 2.6 million deaths worldwide. With public health measures aimed at containing the spread of the disease, several countries have faced a crisis in the availability of intensive care units. Currently, a large-scale effort is underway to identify the nucleotide sequences of the SARS-CoV-2 coronavirus that is an etiological agent of COVID-19. Global sequencing of thousands of viral genomes has revealed many common genetic variants, which enables the monitoring of the evolution of SARS-CoV-2 and the tracking of its spread over time. Understanding the current evolution of SARS-CoV-2 is necessary not only for a retrospective analysis of the new coronavirus infection spread, but also for the development of approaches to the therapy and prophylaxis of COVID-19. In this review, we have focused on the general characteristics of SARS-CoV-2 and COVID-19. Also, we have analyzed available publications on the genetic diversity of the virus and the relationship between the diversity and the biological properties of SARS-CoV-2, such as virulence and contagiousness.

A Novel SARS-CoV-2 Viral Sequence Bioinformatic Pipeline Has Found Genetic Evidence That the Viral 3' Untranslated Region (UTR) Is Evolving and Generating Increased Viral Diversity

An unprecedented amount of SARS-CoV-2 sequencing has been performed, however, novel bioinformatic tools to cope with and process these large datasets is needed. Here, we have devised a bioinformatic pipeline that inputs SARS-CoV-2 genome sequencing in FASTA/FASTQ format and outputs a single Variant Calling Format file that can be processed to obtain variant annotations and perform downstream population genetic testing. As proof of concept, we have analyzed over 229,000 SARS-CoV-2 viral sequences up until November 30, 2020. We have identified over 39,000 variants worldwide with increased polymorphisms, spanning the ORF3a gene as well as the 3' untranslated (UTR) regions, specifically in the conserved stem loop region of SARS-CoV-2 which is accumulating greater observed viral diversity relative to chance variation. Our analysis pipeline has also discovered the existence of SARS-CoV-2 hypermutation with low frequency (less than in 2% of genomes) likely arising through host immune responses and not due to sequencing errors. Among annotated non-sense variants with a population frequency over 1%, recurrent inactivation of the ORF8 gene was found. This was found to be present in the newly identified B.1.1.7 SARS-CoV-2 lineage that originated in the United Kingdom. Almost all VOC-containing genomes possess one stop codon in ORF8 gene (Q27∗), however, 13% of these genomes also contains another stop codon (K68∗), suggesting that ORF8 loss does not interfere with SARS-CoV-2 spread and may play a role in its increased virulence. We have developed this computational pipeline to assist researchers in the rapid analysis and characterization of SARS-CoV-2 variation.

The First Molecular Characterization of Serbian SARS-CoV-2 Isolates From a Unique Early Second Wave in Europe

March 6, 2020 is considered as the official date of the beginning of the COVID-19 epidemic in Serbia. In late spring and early summer 2020, Europe recorded a decline in the rate of SARS-CoV-2 infection and subsiding of the first wave. This trend lasted until the fall, when the second wave of the epidemic began to appear. Unlike the rest of Europe, Serbia was hit by the second wave of the epidemic a few months earlier. Already in June 2020, newly confirmed cases had risen exponentially. As the COVID-19 pandemic is the first pandemic in which there has been instant sharing of genomic information on isolates around the world, the aim of this study was to analyze whole SARS-CoV-2 viral genomes from Serbia, to identify circulating variants/clade/lineages, and to explore site-specific mutational patterns in the unique early second wave of the European epidemic. This analysis of Serbian isolates represents the first publication from Balkan countries, which demonstrates the importance of specificities of local transmission especially when preventive measures differ among countries. One hundred forty-eight different genome variants among 41 Serbian isolates were detected in this study. One unique and seven extremely rare mutations were identified, with locally specific continuous dominance of the 20D clade. At the same time, amino acid substitutions of newly identified variants of concern were found in our isolates from October 2020. Future research should be focused on functional characterization of novel mutations in order to understand the exact role of these variations.

Functions of Coronavirus Accessory Proteins: Overview of the State of the Art

Coronavirus accessory proteins are a unique set of proteins whose genes are interspersed among or within the genes encoding structural proteins. Different coronavirus genera, or even different species within the same coronavirus genus, encode varying amounts of accessory proteins, leading to genus- or species-specificity. Though accessory proteins are dispensable for the replication of coronavirus in vitro, they play important roles in regulating innate immunity, viral proliferation, and pathogenicity. The function of accessory proteins on virus infection and pathogenesis is an area of particular interest. In this review, we summarize the current knowledge on accessory proteins of several representative coronaviruses that infect humans or animals, including the emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), with an emphasis on their roles in interaction between virus and host, mainly involving stress response, innate immunity, autophagy, and apoptosis. The cross-talking among these pathways is also discussed.