Continuous and Discontinuous RNA Synthesis in Coronaviruses

The Coronavirus helicase in replication

The coronavirus nonstructural protein (nsp) 13 encodes an RNA helicase (nsp13-HEL) with multiple enzymatic functions, including unwinding and nucleoside phosphatase (NTPase) activities. Attempts for enzymatic inactivation have defined the nsp13-HEL as a critical enzyme for viral replication and a high-priority target for antiviral development. Helicases have been shown to play numerous roles beyond their canonical ATPase and unwinding activities, though these functions are just beginning to be explored in coronavirus biology. Recent genetic and biochemical studies, as well as work in structurally-related helicases, have provided evidence that supports new hypotheses for the helicase's potential role in coronavirus replication. Here, we review several aspects of the coronavirus nsp13-HEL, including its reported and proposed functions in viral replication and highlight fundamental areas of research that may aid the development of helicase inhibitors.

Targeting sgRNA N secondary structure as a way of inhibiting SARS-CoV-2 replication

SARS-CoV-2 is a betacoronavirus that causes COVID-19, a global pandemic that has resulted in many infections, deaths, and socio-economic challenges. The virus has a large positive-sense, single-stranded RNA genome of ∼30 kb, which produces subgenomic RNAs (sgRNAs) through discontinuous transcription. The most abundant sgRNA is sgRNA N, which encodes the nucleocapsid (N) protein. In this study, we probed the secondary structure of sgRNA N and a shorter model without a 3' UTR in vitro, using the SHAPE (selective 2'-hydroxyl acylation analyzed by a primer extension) method and chemical mapping with dimethyl sulfate and 1-cyclohexyl-(2-morpholinoethyl) carbodiimide metho-p-toluene sulfonate. We revealed the secondary structure of sgRNA N and its shorter variant for the first time and compared them with the genomic RNA N structure. Based on the structural information, we designed gapmers, siRNAs and antisense oligonucleotides (ASOs) to target the N protein coding region of sgRNA N. We also generated eukaryotic expression vectors containing the complete sequence of sgRNA N and used them to screen for new SARS-CoV-2 gene N expression inhibitors. Our study provides novel insights into the structure and function of sgRNA N and potential therapeutic tools against SARS-CoV-2.

The coronavirus nsp15 endoribonuclease: A puzzling protein and pertinent antiviral drug target

The SARS-CoV-2 pandemic has bolstered unprecedented research efforts to better understand the pathogenesis of coronavirus (CoV) infections and develop effective therapeutics. We here focus on non-structural protein nsp15, a hexameric component of the viral replication-transcription complex (RTC). Nsp15 possesses uridine-specific endoribonuclease (EndoU) activity for which some specific cleavage sites were recently identified in viral RNA. By preventing accumulation of viral dsRNA, EndoU helps the virus to evade RNA sensors of the innate immune response. The immune-evading property of nsp15 was firmly established in several CoV animal models and makes it a pertinent target for antiviral therapy. The search for nsp15 inhibitors typically proceeds via compound screenings and is aided by the rapidly evolving insight in the protein structure of nsp15. In this overview, we broadly cover this fascinating protein, starting with its structure, biochemical properties and functions in CoV immune evasion. Next, we summarize the reported studies in which compound screening or a more rational method was used to identify suitable leads for nsp15 inhibitor development. In this way, we hope to raise awareness on the relevance and druggability of this unique CoV protein.

m6A Methylation of Transcription Leader Sequence of SARS-CoV-2 Impacts Discontinuous Transcription of Subgenomic mRNAs

The SARS-CoV-2 genome has been shown to be m6A methylated at several positions in vivo. Strikingly, a DRACH motif, the recognition motif for adenosine methylation, resides in the core of the transcriptional regulatory leader sequence (TRS-L) at position A74, which is highly conserved and essential for viral discontinuous transcription. Methylation at position A74 correlates with viral pathogenicity. Discontinuous transcription produces a set of subgenomic mRNAs that function as templates for translation of all structural and accessory proteins. A74 is base-paired in the short stem-loop structure 5'SL3 that opens during discontinuous transcription to form long-range RNA-RNA interactions with nascent (-)-strand transcripts at complementary TRS-body sequences. A74 can be methylated by the human METTL3/METTL14 complex in vitro. Here, we investigate its impact on the structural stability of 5'SL3 and the long-range TRS-leader:TRS-body duplex formation necessary for synthesis of subgenomic mRNAs of all four viral structural proteins. Methylation uniformly destabilizes 5'SL3 and long-range duplexes and alters their relative equilibrium populations, suggesting that the m6A74 modification acts as a regulator for the abundance of viral structural proteins due to this destabilization.

SARS-CoV-2 replication and drug discovery

The coronavirus disease 2019 (COVID-19) caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has killed millions of people and continues to wreak havoc across the globe. This sudden and deadly pandemic emphasizes the necessity for anti-viral drug development that can be rapidly administered to reduce morbidity, mortality, and virus propagation. Thus, lacking efficient anti-COVID-19 treatment, and especially given the lengthy drug development process as well as the critical death tool that has been associated with SARS-CoV-2 since its outbreak, drug repurposing (or repositioning) constitutes so far, the ideal and ready-to-go best approach in mitigating viral spread, containing the infection, and reducing the COVID-19-associated death rate. Indeed, based on the molecular similarity approach of SARS-CoV-2 with previous coronaviruses (CoVs), repurposed drugs have been reported to hamper SARS-CoV-2 replication. Therefore, understanding the inhibition mechanisms of viral replication by repurposed anti-viral drugs and chemicals known to block CoV and SARS-CoV-2 multiplication is crucial, and it opens the way for particular treatment options and COVID-19 therapeutics. In this review, we highlighted molecular basics underlying drug-repurposing strategies against SARS-CoV-2. Notably, we discussed inhibition mechanisms of viral replication, involving and including inhibition of SARS-CoV-2 proteases (3C-like protease, 3CLpro or Papain-like protease, PLpro) by protease inhibitors such as Carmofur, Ebselen, and GRL017, polymerases (RNA-dependent RNA-polymerase, RdRp) by drugs like Suramin, Remdesivir, or Favipiravir, and proteins/peptides inhibiting virus-cell fusion and host cell replication pathways, such as Disulfiram, GC376, and Molnupiravir. When applicable, comparisons with SARS-CoV inhibitors approved for clinical use were made to provide further insights to understand molecular basics in inhibiting SARS-CoV-2 replication and draw conclusions for future drug discovery research.

3D Puzzle at the Nanoscale-How do RNA Viruses Self-Assemble their Capsids into Perfectly Ordered Structures

The phenomenon of RNA virus self-organization, first observed in the mid-20th century in tobacco mosaic virus, is the subject of extensive research. Efforts to comprehend this process intensify due to its potential for producing vaccines or antiviral compounds as well as nanocarriers and nanotemplates. However, direct observation of the self-assembly is hindered by its prevalence within infected host cells. One of the approaches involves in vitro and in silico research using model viruses featuring a ssRNA(+) genome enclosed within a capsid made up of a single type protein. While various pathways are proposed based on these studies, their relevance in vivo remains uncertain. On the other hand, the development of advanced microscopic methods provide insights into the events within living cells, where following viral infection, specialized compartments form to facilitate the creation of nascent virions. Intriguingly, a growing body of evidence indicates that the primary function of packaging signals in viral RNA is to effectively initiate the virion self-assembly. This is in contrast to earlier opinions suggesting a role in marking RNA for encapsidation. Another noteworthy observation is that many viruses undergo self-assembly within membraneless liquid organelles, which are specifically induced by viral proteins.

Progress of research on coronaviruses and toroviruses in large domestic animals using reverse genetics systems

The generation of genetically engineered recombinant viruses from modified DNA/RNA is commonly referred to as reverse genetics, which allows the introduction of desired mutations into the viral genome. Reverse genetics systems (RGSs) are powerful tools for studying fundamental viral processes, mechanisms of infection, pathogenesis and vaccine development. However, establishing RGS for coronaviruses (CoVs) and toroviruses (ToVs), which have the largest genomes among vertebrate RNA viruses, is laborious and hampered by technical constraints. Hence, little research has focused on animal CoVs and ToVs using RGSs, especially in large domestic animals such as pigs and cattle. In the last decade, however, studies of porcine CoVs and bovine ToVs using RGSs have been reported. In addition, the coronavirus disease-2019 pandemic has prompted the development of new and simple CoV RGSs, which will accelerate RGS-based research on animal CoVs and ToVs. In this review, we summarise the general characteristics of CoVs and ToVs, the RGSs available for CoVs and ToVs and the progress made in the last decade in RGS-based research on porcine CoVs and bovine ToVs.

SARS-CoV-2 subgenomic RNA: formation process and rapid molecular diagnostic methods

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) which caused coronavirus disease-2019 (COVID-19) is spreading worldwide and posing enormous losses to human health and socio-economic. Due to the limitations of medical and health conditions, it is still a huge challenge to develop appropriate discharge standards for patients with COVID-19 and to use medical resources in a timely and effective manner. Similar to other coronaviruses, SARS-CoV-2 has a very complex discontinuous transcription process to generate subgenomic RNA (sgRNA). Some studies support that sgRNA of SARS-CoV-2 can only exist when the virus is active and is an indicator of virus replication. The results of sgRNA detection in patients can be used to evaluate the condition of hospitalized patients, which is expected to save medical resources, especially personal protective equipment. There have been numerous investigations using different methods, especially molecular methods to detect sgRNA. Here, we introduce the process of SARS-CoV-2 sgRNA formation and the commonly used molecular diagnostic methods to bring a new idea for clinical detection in the future.

The human microbiota is a beneficial reservoir for SARS-CoV-2 mutations

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mutations are rapidly emerging. In particular, beneficial mutations in the spike (S) protein, which can either make a person more infectious or enable immunological escape, are providing a significant obstacle to the prevention and treatment of pandemics. However, how the virus acquires a high number of beneficial mutations in a short time remains a mystery. We demonstrate here that variations of concern may be mutated due in part to the influence of the human microbiome. We searched the National Center for Biotechnology Information database for homologous fragments (HFs) after finding a mutation and the six neighboring amino acids in a viral mutation fragment. Among the approximate 8,000 HFs obtained, 61 mutations in S and other outer membrane proteins were found in bacteria, accounting for 62% of all mutation sources, which is 12-fold higher than the natural variable proportion. A significant proportion of these bacterial species-roughly 70%-come from the human microbiota, are mainly found in the lung or gut, and share a composition pattern with COVID-19 patients. Importantly, SARS-CoV-2 RNA-dependent RNA polymerase replicates corresponding bacterial mRNAs harboring mutations, producing chimeric RNAs. SARS-CoV-2 may collectively pick up mutations from the human microbiota that change the original virus's binding sites or antigenic determinants. Our study clarifies the evolving mutational mechanisms of SARS-CoV-2.

IMPORTANCE

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mutations are rapidly emerging, in particular advantageous mutations in the spike (S) protein, which either increase transmissibility or lead to immune escape and are posing a major challenge to pandemic prevention and treatment. However, how the virus acquires a high number of advantageous mutations in a short time remains a mystery. Here, we provide evidence that the human microbiota is a reservoir of advantageous mutations and aids mutational evolution and host adaptation of SARS-CoV-2. Our findings demonstrate a conceptual breakthrough on the mutational evolution mechanisms of SARS-CoV-2 for human adaptation. SARS-CoV-2 may grab advantageous mutations from the widely existing microorganisms in the host, which is undoubtedly an "efficient" manner. Our study might open a new perspective to understand the evolution of virus mutation, which has enormous implications for comprehending the trajectory of the COVID-19 pandemic.

De novo gene synthesis by an antiviral reverse transcriptase

Bacteria defend themselves from viral infection using diverse immune systems, many of which sense and target foreign nucleic acids. Defense-associated reverse transcriptase (DRT) systems provide an intriguing counterpoint to this immune strategy by instead leveraging DNA synthesis, but the identities and functions of their DNA products remain largely unknown. Here we show that DRT2 systems execute an unprecedented immunity mechanism that involves de novo gene synthesis via rolling-circle reverse transcription of a non-coding RNA (ncRNA). Unbiased profiling of RT-associated RNA and DNA ligands in DRT2-expressing cells revealed that reverse transcription generates concatenated cDNA repeats through programmed template jumping on the ncRNA. The presence of phage then triggers second-strand cDNA synthesis, leading to the production of long double-stranded DNA. Remarkably, this DNA product is efficiently transcribed, generating messenger RNAs that encode a stop codon-less, never-ending ORF (neo) whose translation causes potent growth arrest. Phylogenetic analyses and screening of diverse DRT2 homologs further revealed broad conservation of rolling-circle reverse transcription and Neo protein function. Our work highlights an elegant expansion of genome coding potential through RNA-templated gene creation, and challenges conventional paradigms of genetic information encoded along the one-dimensional axis of genomic DNA.

ADLM Guidance Document on Laboratory Diagnosis of Respiratory Viruses

Respiratory viral infections are among the most frequent infections experienced worldwide. The COVID-19 pandemic has highlighted the need for testing and currently several tests are available for the detection of a wide range of viruses. These tests vary widely in terms of the number of viral pathogens included, viral markers targeted, regulatory status, and turnaround time to results, as well as their analytical and clinical performance. Given these many variables, selection and interpretation of testing requires thoughtful consideration. The current guidance document is the authors' expert opinion based on the preponderance of available evidence to address key questions related to best practices for laboratory diagnosis of respiratory viral infections including who to test, when to test, and what tests to use. An algorithm is proposed to help laboratories decide on the most appropriate tests to use for the diagnosis of respiratory viral infections.

Role of marine natural products in the development of antiviral agents against SARS-CoV-2: potential and prospects

A novel coronavirus, known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has surfaced and caused global concern owing to its ferocity. SARS-CoV-2 is the causative agent of coronavirus disease 2019; however, it was only discovered at the end of the year and was considered a pandemic by the World Health Organization. Therefore, the development of novel potent inhibitors against SARS-CoV-2 and future outbreaks is urgently required. Numerous naturally occurring bioactive substances have been studied in the clinical setting for diverse disorders. The intricate infection and replication mechanism of SARS-CoV-2 offers diverse therapeutic drug targets for developing antiviral medicines by employing natural products that are safer than synthetic compounds. Marine natural products (MNPs) have received increased attention in the development of novel drugs owing to their high diversity and availability. Therefore, this review article investigates the infection and replication mechanisms, including the function of the SARS-CoV-2 genome and structure. Furthermore, we highlighted anti-SARS-CoV-2 therapeutic intervention efforts utilizing MNPs and predicted SARS-CoV-2 inhibitor design.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00215-9.

Prediction of Prospective Mutational Landscape of SARS-CoV-2 Spike ssRNA and Evolutionary Basis of Its Host Interaction

Booster doses are crucial against severe COVID-19, as rapid virus mutations and variant emergence prolong the pandemic crisis. The virus's quick evolution, short generation-time, and adaptive changes impact virulence and evolvability, helping predictions about variant of concerns' (VOCs') landscapes. Here, in this study, we used a new computational algorithm, to predict the mutational pattern in SARS-CoV-2 ssRNA, proteomics, structural identification, mutation stability, and functional correlation, as well as immune escape mechanisms. Interestingly, the sequence diversity of SARS Coronavirus-2 has demonstrated a predominance of G- > A and C- > U substitutions. The best validation statistics are explored here in seven homologous models of the expected mutant SARS-CoV-2 spike ssRNA and employed for hACE2 and IgG interactions. The interactome profile of SARS-CoV-2 spike with hACE2 and IgG revealed a strong correlation between phylogeny and divergence time. The systematic adaptation of SARS-CoV-2 spike ssRNA influences infectivity and immune escape. Data suggest higher propensity of Adenine rich sequence promotes MHC system avoidance, preferred by A-rich codons. Phylogenetic data revealed the evolution of SARS-CoV-2 lineages' epidemiology. Our findings may unveil processes governing the genesis of immune-resistant variants, prompting a critical reassessment of the coronavirus mutation rate and exploration of hypotheses beyond mechanical aspects.

Identification of Hammerhead-variant ribozyme sequences in SARS-CoV-2

The SARS-CoV-2 RNA virus and variants, responsible for the COVID-19 pandemic has become endemic, raised a need for further understanding of the viral genome and biology. Despite vast research on SARS-CoV-2, no ribozymes have been found in the virus genome. Here we report the identification of 39 Hammerhead-variant ribozyme sequences (CoV-HHRz) in SARS-CoV-2. These sequences are highly conserved within SARS-CoV-2 variants but show large diversity among other coronaviruses. In vitro CoV-HHRz sequences possess the characteristics of typical ribozymes; cleavage is pH and ion dependent, although their activity is relatively low and Mn2+ is required for cleavage. The cleavage sites of four CoV-HHRz coincide with the breakpoint of expressed subgenomic RNA (sgRNAs) in SARS-CoV-2 transcriptome data suggesting in vivo activity. The CoV-HHRz are involved in processing sgRNAs for ORF7b, ORF 10 and ORF1ab nsp13 which are essential for viral packaging and life cycle.

ER-export and ARFRP1/AP-1-dependent delivery of SARS-CoV-2 Envelope to lysosomes controls late stages of viral replication

The β-coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the global COVID-19 pandemic. Coronaviral Envelope (E) proteins are pentameric viroporins that play essential roles in assembly, release, and pathogenesis. We developed a nondisruptive tagging strategy for SARS-CoV-2 E and find that, at steady state, it localizes to the Golgi and to lysosomes. We identify sequences in E, conserved across Coronaviridae, responsible for endoplasmic reticulum-to-Golgi export, and relate this activity to interaction with COP-II via SEC24. Using proximity biotinylation, we identify an ADP ribosylation factor 1/adaptor protein-1 (ARFRP1/AP-1)-dependent pathway allowing Golgi-to-lysosome trafficking of E. We identify sequences in E that bind AP-1, are conserved across β-coronaviruses, and allow E to be trafficked from Golgi to lysosomes. We show that E acts to deacidify lysosomes and, by developing a trans-complementation assay for SARS-CoV-2 structural proteins, that lysosomal delivery of E and its viroporin activity is necessary for efficient viral replication and release.

Experimental and analytical pipeline for sub-genomic RNA landscape of coronavirus by Nanopore sequencer

Coronaviruses (CoVs), including severe acute respiratory syndrome coronavirus 2, can infect a variety of mammalian and avian hosts with significant medical and economic consequences. During the life cycle of CoV, a coordinated series of subgenomic RNAs, including canonical subgenomic messenger RNA and non-canonical defective viral genomes (DVGs), are generated with different biological implications. Studies that adopted the Nanopore sequencer (ONT) to investigate the landscape and dynamics of viral RNA subgenomic transcriptomes applied arbitrary bioinformatics parameters without justification or experimental validation. The current study used bovine coronavirus (BCoV), which can be performed under biosafety level 2 for library construction and experimental validation using traditional colony polymerase chain reaction and Sanger sequencing. Four different ONT protocols, including RNA direct and cDNA direct sequencing with or without exonuclease treatment, were used to generate RNA transcriptomic libraries from BCoV-infected cell lysates. Through rigorously examining the k-mer, gap size, segment size, and bin size, the optimal cutoffs for the bioinformatic pipeline were determined to remove the sequence noise while keeping the informative DVG reads. The sensitivity and specificity of identifying DVG reads using the proposed pipeline can reach 82.6% and 99.6% under the k-mer size cutoff of 15. Exonuclease treatment reduced the abundance of RNA transcripts; however, it was not necessary for future library preparation. Additional recovery of clipped BCoV nucleotide sequences with experimental validation expands the landscape of the CoV discontinuous RNA transcriptome, whose biological function requires future investigation. The results of this study provide the benchmarks for library construction and bioinformatic parameters for studying the discontinuous CoV RNA transcriptome.IMPORTANCEFunctional defective viral genomic RNA, containing all the cis-acting elements required for translation or replication, may play different roles in triggering cell innate immune signaling, interfering with the canonical subgenomic messenger RNA transcription/translation or assisting in establishing persistence infection. This study does not only provide benchmarks for library construction and bioinformatic parameters for studying the discontinuous coronavirus RNA transcriptome but also reveals the complexity of the bovine coronavirus transcriptome, whose functional assays will be critical in future studies.

SARS-CoV-2 biology and host interactions

The zoonotic emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the ensuing coronavirus disease 2019 (COVID-19) pandemic have profoundly affected our society. The rapid spread and continuous evolution of new SARS-CoV-2 variants continue to threaten global public health. Recent scientific advances have dissected many of the molecular and cellular mechanisms involved in coronavirus infections, and large-scale screens have uncovered novel host-cell factors that are vitally important for the virus life cycle. In this Review, we provide an updated summary of the SARS-CoV-2 life cycle, gene function and virus-host interactions, including recent landmark findings on general aspects of coronavirus biology and newly discovered host factors necessary for virus replication.

Quantitative and qualitative subgenomic RNA profiles of SARS-CoV-2 in respiratory samples: A comparison between Omicron BA.2 and non-VOC-D614G

The SARS-CoV-2 Omicron variants are notorious for their transmissibility, but little is known about their subgenomic RNA (sgRNA) expression. This study applied RNA-seq to delineate the quantitative and qualitative profiles of canonical sgRNA of 118 respiratory samples collected from patients infected with Omicron BA.2 and compared with 338 patients infected with non-variant of concern (non-VOC)-D614G. A unique characteristic profile depicted by the relative abundance of 9 canonical sgRNAs was reproduced by both BA.2 and non-VOC-D614G regardless of host gender, age and presence of pneumonia. Remarkably, such profile was lost in samples with low viral load, suggesting a potential application of sgRNA pattern to indicate viral activity of individual patient at a specific time point. A characteristic qualitative profile of canonical sgRNAs was also reproduced by both BA.2 and non-VOC-D614G. The presence of a full set of canonical sgRNAs carried a coherent correlation with crude viral load (AUC ​= ​0.91, 95% CI 0.88-0.94), and sgRNA ORF7b was identified to be the best surrogate marker allowing feasible routine application in characterizing the infection status of individual patient. Further potentials in using sgRNA as a target for vaccine and antiviral development are worth pursuing.

The leader RNA of SARS-CoV-2 sequesters polypyrimidine tract binding protein (PTBP1) and influences pre-mRNA splicing in infected cells

The large amount of viral RNA produced during infections has the potential to interact with and effectively sequester cellular RNA binding proteins, thereby influencing aspects of post-transcriptional gene regulation in the infected cell. Here we demonstrate that the abundant 5' leader RNA region of SARS-CoV-2 viral RNAs can interact with the cellular polypyrimidine tract binding protein (PTBP1). Interestingly, the effect of a knockdown of PTBP1 protein on cellular gene expression is also mimicked during SARS-CoV-2 infection, suggesting that this protein may be functionally sequestered by viral RNAs. Consistent with this model, the alternative splicing of mRNAs that is normally controlled by PTBP1 is dysregulated during SARS-CoV-2 infection. Collectively, these data suggest that the SARS-CoV-2 leader RNA sequesters the cellular PTBP1 protein during infection, resulting in significant impacts on the RNA biology of the host cell. These alterations in post-transcriptional gene regulation may play a role in SARS-CoV-2 mediated molecular pathogenesis.

Activity and inhibition of the SARS-CoV-2 Omicron nsp13 R392C variant using RNA duplex unwinding assays

SARS-CoV-2 nsp13 helicase is an essential enzyme for viral replication and a promising target for antiviral drug development. This study compares the double-stranded RNA (dsRNA) unwinding activity of nsp13 and the Omicron nsp13R392C variant, which is predominant in currently circulating lineages. Using in vitro gel- and fluorescence-based assays, we found that both nsp13 and nsp13R392C have dsRNA unwinding activity with equivalent kinetics. Furthermore, the R392C mutation had no effect on the efficiency of the nsp13-specific helicase inhibitor SSYA10-001. We additionally confirmed the activity of several other helicase inhibitors against nsp13, including punicalagin that inhibited dsRNA unwinding at nanomolar concentrations. Overall, this study reveals the utility of using dsRNA unwinding assays to screen small molecules for antiviral activity against nsp13 and the Omicron nsp13R392C variant. Continual monitoring of newly emergent variants will be essential for considering resistance profiles of lead compounds as they are advanced towards next-generation therapeutic development.

Coinfection and Nonrandom Recombination Drive the Evolution of Swine Enteric Coronaviruses

Coinfection with multiple viruses is a common phenomenon in clinical settings and is a crucial driver of viral evolution. Although numerous studies have demonstrated viral recombination arising from coinfections of different strains of a specific species, the role of coinfections of different species or genera during viral evolution is rarely investigated. Here, we analyzed coinfections of and recombination events between four different swine enteric coronaviruses that infect the jejunum and ileum in pigs, including porcine epidemic diarrhea virus (PEDV), transmissible gastroenteritis virus (TGEV), and swine acute diarrhea syndrome coronavirus (SADS-CoV), and a deltacoronavirus, porcine deltacoronavirus (PDCoV). Various coinfection patterns were observed in 4,468 fecal and intestinal tissue samples collected from pigs in a 4-year survey. PEDV/PDCoV was the most frequent coinfection. However, recombination analyses have only detected events involving PEDV/TGEV and SADS-CoV/TGEV, indicating that inter-species recombination among coronaviruses is most likely to occur within the same genus. We also analyzed recombination events within the newly identified genus Deltacoronavirus and found that sparrows have played a unique host role in the recombination history of the deltacoronaviruses. The emerging virus PDCoV, which can infect humans, has a different recombination history. In summary, our study demonstrates that swine enteric coronaviruses are a valuable model for investigating the relationship between viral coinfection and recombination, which provide new insights into both inter- and intraspecies recombination events among swine enteric coronaviruses, and extend our understanding of the relationship between coronavirus coinfection and recombination.

[Genetic factors associated with long COVID]

INTRODUCTION

The variability in expression and evolution of COVID is not completely explained by clinical factors. In fact, genetic factors play an important role. Moreover, it is unknown whether the genetic factor that contribute to susceptibility and severity are also involved in the onset and evolution of long-COVID. The objective of this review is to gather information from literature to understand which genetic factors are involved in the onset of persistent COVID.

MATERIAL AND METHODS

Systematic review in PubMed and bioRxiv and medRxiv repositories based on MeSH-descriptors and MeSH-terms related to COVID and genetic factors. Using these terms 2715 articles were pooled. An initial screening performed by authors independently, selected 205 articles of interest. A final deeper screening a total of 85 articles were chosen for complete reading and summarized in this review.

RESULTS

Although ACE2 and TMPSS6 are involved in COVID susceptibility, their involvement in long-COVID has not been found. On the other hand, the severity of the disease and the onset of long-COVID has been associated with different genes involved in the inflammatory and immune response. Particularly interesting has been the association found with the FOXP4 locus.

CONCLUSIONS

Although studies on long-COVID are insufficient to fully comprehend the cause, it is clear that the current identified genetic factors do not fully explain the progression and onset of long-COVID. Other factors such as polygenic action, pleiotropic genes, the microbiota and epigenetic changes must be considered and studied.

Classification, replication, and transcription of Nidovirales

Nidovirales is one order of RNA virus, with the largest single-stranded positive sense RNA genome enwrapped with membrane envelope. It comprises four families (Arterividae, Mesoniviridae, Roniviridae, and Coronaviridae) and has been circulating in humans and animals for almost one century, posing great threat to livestock and poultry，as well as to public health. Nidovirales shares similar life cycle: attachment to cell surface, entry, primary translation of replicases, viral RNA replication in cytoplasm, translation of viral proteins, virion assembly, budding, and release. The viral RNA synthesis is the critical step during infection, including genomic RNA (gRNA) replication and subgenomic mRNAs (sg mRNAs) transcription. gRNA replication requires the synthesis of a negative sense full-length RNA intermediate, while the sg mRNAs transcription involves the synthesis of a nested set of negative sense subgenomic intermediates by a discontinuous strategy. This RNA synthesis process is mediated by the viral replication/transcription complex (RTC), which consists of several enzymatic replicases derived from the polyprotein 1a and polyprotein 1ab and several cellular proteins. These replicases and host factors represent the optimal potential therapeutic targets. Hereby, we summarize the Nidovirales classification, associated diseases, "replication organelle," replication and transcription mechanisms, as well as related regulatory factors.

The most exposed regions of SARS-CoV-2 structural proteins are subject to strong positive selection and gene overlap may locally modify this behavior

The SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) pandemic that emerged in 2019 has been an unprecedented event in international science, as it has been possible to sequence millions of genomes, tracking their evolution very closely. This has enabled various types of secondary analyses of these genomes, including the measurement of their sequence selection pressure. In this work, we have been able to measure the selective pressure of all the described SARS-CoV-2 genes, even analyzed by sequence regions, and we show how this type of analysis allows us to separate the genes between those subject to positive selection (usually those that code for surface proteins or those exposed to the host immune system) and those subject to negative selection because they require greater conservation of their structure and function. We have also seen that when another gene with an overlapping reading frame appears within a gene sequence, the overlapping sequence between the two genes evolves under a stronger purifying selection than the average of the non-overlapping regions of the main gene. We propose this type of analysis as a useful tool for locating and analyzing all the genes of a viral genome when an adequate number of sequences are available.IMPORTANCEWe have analyzed the selection pressure of all severe acute respiratory syndrome coronavirus 2 genes by means of the nonsynonymous (Ka) to synonymous (Ks) substitution rate. We found that protein-coding genes are exposed to strong positive selection, especially in the regions of interaction with other molecules (host receptor and genome of the virus itself). However, overlapping coding regions are more protected and show negative selection. This suggests that this measure could be used to study viral gene function as well as overlapping genes.

Experimental and computational methods for studying the dynamics of RNA-RNA interactions in SARS-COV2 genomes

Long-range ribonucleic acid (RNA)-RNA interactions (RRI) are prevalent in positive-strand RNA viruses, including Beta-coronaviruses, and these take part in regulatory roles, including the regulation of sub-genomic RNA production rates. Crosslinking of interacting RNAs and short read-based deep sequencing of resulting RNA-RNA hybrids have shown that these long-range structures exist in severe acute respiratory syndrome coronavirus (SARS-CoV)-2 on both genomic and sub-genomic levels and in dynamic topologies. Furthermore, co-evolution of coronaviruses with their hosts is navigated by genetic variations made possible by its large genome, high recombination frequency and a high mutation rate. SARS-CoV-2's mutations are known to occur spontaneously during replication, and thousands of aggregate mutations have been reported since the emergence of the virus. Although many long-range RRIs have been experimentally identified using high-throughput methods for the wild-type SARS-CoV-2 strain, evolutionary trajectory of these RRIs across variants, impact of mutations on RRIs and interaction of SARS-CoV-2 RNAs with the host have been largely open questions in the field. In this review, we summarize recent computational tools and experimental methods that have been enabling the mapping of RRIs in viral genomes, with a specific focus on SARS-CoV-2. We also present available informatics resources to navigate the RRI maps and shed light on the impact of mutations on the RRI space in viral genomes. Investigating the evolution of long-range RNA interactions and that of virus-host interactions can contribute to the understanding of new and emerging variants as well as aid in developing improved RNA therapeutics critical for combating future outbreaks.

Defining Distinct RNA-Protein Interactomes of SARS-CoV-2 Genomic and Subgenomic RNAs

Host RNA binding proteins recognize viral RNA and play key roles in virus replication and antiviral mechanisms. SARS-CoV-2 generates a series of tiered subgenomic RNAs (sgRNAs), each encoding distinct viral protein(s) that regulate different aspects of viral replication. Here, for the first time, we demonstrate the successful isolation of SARS-CoV-2 genomic RNA and three distinct sgRNAs (N, S, and ORF8) from a single population of infected cells and characterize their protein interactomes. Over 500 protein interactors (including 260 previously unknown) were identified as associated with one or more target RNA. These included protein interactors unique to a single RNA pool and others present in multiple pools, highlighting our ability to discriminate between distinct viral RNA interactomes despite high sequence similarity. Individual interactomes indicated viral associations with cell response pathways, including regulation of cytoplasmic ribonucleoprotein granules and posttranscriptional gene silencing. We tested the significance of three protein interactors in these pathways (APOBEC3F, PPP1CC, and MSI2) using siRNA knockdowns, with several knockdowns affecting viral gene expression, most consistently PPP1CC. This study describes a new technology for high-resolution studies of SARS-CoV-2 RNA regulation and reveals a wealth of new viral RNA-associated host factors of potential functional significance to infection.

Identification of subgenomic mRNAs derived from the coronavirus 1a/1b protein gene: Implications for coronavirus transcription

Synthesis of coronavirus subgenomic mRNA (sgmRNA) is guided by the transcription regulatory sequence (TRS). sgmRNA derived from the body TRS (TRS-B) located at the 1a/1b protein gene is designated 1ab/sgmRNA. In the current study, we comprehensively identified the 1ab/sgmRNAs synthesized from TRS-Bs located at the 1a/1b protein genes of different coronavirus genera both in vitro and in vivo by RT‒PCR and sequencing. The results suggested that the degree of sequence homology between the leader TRS (TRS-L) and TRS-B may not be a decisive factor for 1ab/sgmRNA synthesis. This observation led us to revisit the coronavirus transcription mechanism and to propose that the disassociation of coronavirus polymerase from the viral genome may be a prerequisite for sgmRNA synthesis. Once the polymerase can disassociate at TRS-B, the sequence homology between TRS-L and TRS-B is important for sgmRNA synthesis. The study therefore extends our understanding of transcription mechanisms.

Intermolecular base-pairing interactions, a unique topology and exoribonuclease-resistant noncoding RNAs drive formation of viral chimeric RNAs in plants

In plants, exoribonuclease-resistant RNAs (xrRNAs) are produced by many viruses. Whereas xrRNAs contribute to the pathogenicity of these viruses, the role of xrRNAs in the virus infectious cycle remains elusive. Here, we show that xrRNAs produced by a benyvirus (a multipartite RNA virus with four genomic segments) in plants are involved in the formation of monocistronic coat protein (CP)-encoding chimeric RNAs. Naturally occurring chimeric RNAs, we discovered, are composed of 5'-end of RNA 2 and 3'-end of either RNA 3 or RNA 4 bearing conservative exoribonuclease-resistant 'coremin' region. Using computational tools and site-directed mutagenesis, we show that de novo formation of chimeric RNAs requires intermolecular base-pairing interaction between 'coremin' and 3'-proximal part of the CP gene of RNA 2 as well as a stem-loop structure immediately adjacent to the CP gene. Moreover, knockdown of the expression of the XRN4 gene, encoding 5'→3' exoribonuclease, inhibits biogenesis of both xrRNAs and chimeric RNAs. Our findings suggest a novel mechanism involving a unique tropology of the intermolecular base-pairing complex between xrRNAs and RNA2 to promote formation of chimeric RNAs in plants. XrRNAs, essential for chimeric RNA biogenesis, are generated through the action of cytoplasmic Xrn 4 5'→3' exoribonuclease conserved in all plant species.

Unveiling the role of PUS7-mediated pseudouridylation in host protein interactions specific for the SARS-CoV-2 RNA genome

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a positive single-stranded RNA virus, engages in complex interactions with host cell proteins throughout its life cycle. While these interactions enable the host to recognize and inhibit viral replication, they also facilitate essential viral processes such as transcription, translation, and replication. Many aspects of these virus-host interactions remain poorly understood. Here, we employed the catRAPID algorithm and utilized the RNA-protein interaction detection coupled with mass spectrometry technology to predict and validate the host proteins that specifically bind to the highly structured 5' and 3' terminal regions of the SARS-CoV-2 RNA. Among the interactions identified, we prioritized pseudouridine synthase PUS7, which binds to both ends of the viral RNA. Using nanopore direct RNA sequencing, we discovered that the viral RNA undergoes extensive post-transcriptional modifications. Modified consensus regions for PUS7 were identified at both terminal regions of the SARS-CoV-2 RNA, including one in the viral transcription regulatory sequence leader. Collectively, our findings offer insights into host protein interactions with the SARS-CoV-2 UTRs and highlight the likely significance of pseudouridine synthases and other post-transcriptional modifications in the viral life cycle. This new knowledge enhances our understanding of virus-host dynamics and could inform the development of targeted therapeutic strategies.

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

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

Three-Dimensional Remodeling of SARS-CoV2-Infected Cells Revealed by Cryogenic Soft X-ray Tomography

Plus-strand RNA viruses are proficient at remodeling host cell membranes for optimal viral genome replication and the production of infectious progeny. These ultrastructural alterations result in the formation of viral membranous organelles and may be observed by different imaging techniques, providing nanometric resolution. Guided by confocal and electron microscopy, this study describes the generation of wide-field volumes using cryogenic soft-X-ray tomography (cryo-SXT) on SARS-CoV-2-infected human lung adenocarcinoma cells. Confocal microscopy showed accumulation of double-stranded RNA (dsRNA) and nucleocapsid (N) protein in compact perinuclear structures, preferentially found around centrosomes at late stages of the infection. Transmission electron microscopy (TEM) showed accumulation of membranous structures in the vicinity of the infected cell nucleus, forming a viral replication organelle containing characteristic double-membrane vesicles and virus-like particles within larger vesicular structures. Cryo-SXT revealed viral replication organelles very similar to those observed by TEM but indicated that the vesicular organelle observed in TEM sections is indeed a vesiculo-tubular network that is enlarged and elongated at late stages of the infection. Overall, our data provide additional insight into the molecular architecture of the SARS-CoV-2 replication organelle.

SARS-CoV-2 and the DNA damage response

The recent coronavirus disease 2019 (COVID-19) pandemic was caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). COVID-19 is characterized by respiratory distress, multiorgan dysfunction and, in some cases, death. The virus is also responsible for post-COVID-19 condition (commonly referred to as 'long COVID'). SARS-CoV-2 is a single-stranded, positive-sense RNA virus with a genome of approximately 30 kb, which encodes 26 proteins. It has been reported to affect multiple pathways in infected cells, resulting, in many cases, in the induction of a 'cytokine storm' and cellular senescence. Perhaps because it is an RNA virus, replicating largely in the cytoplasm, the effect of SARS-Cov-2 on genome stability and DNA damage responses (DDRs) has received relatively little attention. However, it is now becoming clear that the virus causes damage to cellular DNA, as shown by the presence of micronuclei, DNA repair foci and increased comet tails in infected cells. This review considers recent evidence indicating how SARS-CoV-2 causes genome instability, deregulates the cell cycle and targets specific components of DDR pathways. The significance of the virus's ability to cause cellular senescence is also considered, as are the implications of genome instability for patients suffering from long COVID.

SARS-CoV-2 ORF8 accessory protein is a virulence factor

IMPORTANCE

The relevance of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) ORF8 in the pathogenesis of COVID-19 is unclear. Virus natural isolates with deletions in ORF8 were associated with wild milder disease, suggesting that ORF8 might contribute to SARS-CoV-2 virulence. This manuscript shows that ORF8 is involved in inflammation and in the activation of macrophages in two experimental systems: humanized K18-hACE2 transgenic mice and organoid-derived human airway cells. These results identify ORF8 protein as a potential target for COVID-19 therapies.

On a path toward a broad-spectrum anti-viral: inhibition of HIV-1 and coronavirus replication by SR kinase inhibitor harmine

IMPORTANCE

This study highlights the crucial role RNA processing plays in regulating viral gene expression and replication. By targeting SR kinases, we identified harmine as a potent inhibitor of HIV-1 as well as coronavirus (HCoV-229E and multiple SARS-CoV-2 variants) replication. Harmine inhibits HIV-1 protein expression and reduces accumulation of HIV-1 RNAs in both cell lines and primary CD4+ T cells. Harmine also suppresses coronavirus replication post-viral entry by preferentially reducing coronavirus sub-genomic RNA accumulation. By focusing on host factors rather than viral targets, our study offers a novel approach to combating viral infections that is effective against a range of unrelated viruses. Moreover, at doses required to inhibit virus replication, harmine had limited toxicity and minimal effect on the host transcriptome. These findings support the viability of targeting host cellular processes as a means of developing broad-spectrum anti-virals.

SND1 binds SARS-CoV-2 negative-sense RNA and promotes viral RNA synthesis through NSP9

Regulation of viral RNA biogenesis is fundamental to productive SARS-CoV-2 infection. To characterize host RNA-binding proteins (RBPs) involved in this process, we biochemically identified proteins bound to genomic and subgenomic SARS-CoV-2 RNAs. We find that the host protein SND1 binds the 5' end of negative-sense viral RNA and is required for SARS-CoV-2 RNA synthesis. SND1-depleted cells form smaller replication organelles and display diminished virus growth kinetics. We discover that NSP9, a viral RBP and direct SND1 interaction partner, is covalently linked to the 5' ends of positive- and negative-sense RNAs produced during infection. These linkages occur at replication-transcription initiation sites, consistent with NSP9 priming viral RNA synthesis. Mechanistically, SND1 remodels NSP9 occupancy and alters the covalent linkage of NSP9 to initiating nucleotides in viral RNA. Our findings implicate NSP9 in the initiation of SARS-CoV-2 RNA synthesis and unravel an unsuspected role of a cellular protein in orchestrating viral RNA production.

Chemical-guided SHAPE sequencing (cgSHAPE-seq) informs the binding site of RNA-degrading chimeras targeting SARS-CoV-2 5' untranslated region

One of the hallmarks of RNA viruses is highly structured untranslated regions (UTRs) in their genomes. These conserved RNA structures are often essential for viral replication, transcription, or translation. In this report, we discovered and optimized a new type of coumarin derivatives, such as C30 and C34, which bind to a four-way RNA helix called SL5 in the 5' UTR of the SARS-CoV-2 RNA genome. To locate the binding site, we developed a novel sequencing-based method namely cgSHAPE-seq, in which the acylating chemical probe was directed to crosslink with the 2'-OH groups of ribose at the ligand binding site. This crosslinked RNA could then create read-through mutations during reverse transcription (i.e., primer extension) at single-nucleotide resolution to uncover the acylation locations. cgSHAPE-seq unambiguously determined that a bulged G in SL5 was the primary binding site of C30 in the SARS-CoV-2 5' UTR, which was validated through mutagenesis and in vitro binding experiments. C30 was further used as a warhead in RNA-degrading chimeras to reduce viral RNA expression levels. We demonstrated that replacing the acylating moiety in the cgSHAPE probe with ribonuclease L recruiter (RLR) moieties yielded RNA degraders active in the in vitro RNase L degradation assay and SARS-CoV-2 5' UTR expressing cells. We further explored another RLR conjugation site on the E ring of C30/C34 and discovered improved RNA degradation activities in vitro and in cells. The optimized RNA-degrading chimera C64 inhibited live virus replication in lung epithelial carcinoma cells.

An Evolutionarily Conserved Strategy for Ribosome Binding and Host Translation Inhibition by β-coronavirus Non-structural Protein 1

An important pathogenicity factor of SARS-CoV-2 and related coronaviruses is Non-structural protein 1 (Nsp1), which suppresses host gene expression and stunts antiviral signaling. SARS-CoV-2 Nsp1 binds the ribosome to inhibit translation through mRNA displacement and induces degradation of host mRNAs. Here we show that Nsp1-dependent host shutoff is conserved in diverse coronaviruses, but only Nsp1 from β-Coronaviruses (β-CoV) inhibits translation through ribosome binding. The C-terminal domain (CTD) of all β-CoV Nsp1s confers high-affinity ribosome binding despite low sequence conservation. Modeling of interactions of four Nsp1s with the ribosome identified the few absolutely conserved amino acids that, together with an overall conservation in surface charge, form the β-CoV Nsp1 ribosome-binding domain. Contrary to previous models, the Nsp1 ribosome-binding domain is an inefficient translation inhibitor. Instead, the Nsp1-CTD likely functions by recruiting Nsp1's N-terminal "effector" domain. Finally, we show that a cis-acting viral RNA element has co-evolved to fine-tune SARS-CoV-2 Nsp1 function, but does not provide similar protection against Nsp1 from related viruses. Together, our work provides new insight into the diversity and conservation of ribosome-dependent host-shutoff functions of Nsp1, knowledge that could aid future efforts in pharmacological targeting of Nsp1 from SARS-CoV-2 and related human-pathogenic β-CoVs. Our study also exemplifies how comparing highly divergent Nsp1 variants can help to dissect the different modalities of this multi-functional viral protein.

Comprehensive Risk Assessment of Infection Induced by SARS-CoV-2

The pandemic COVID-19 (coronavirus disease 2019) is caused by the severe acute respiratory syndrome corona virus 2 (SARS-CoV-2), which devastated the global economy and healthcare system. The infection caused an unforeseen rise in COVID-19 patients and increased the mortality rate globally. This study gives an overall idea about host-pathogen interaction, immune responses to COVID-19, recovery status of infection, targeted organs and complications associated, and comparison of post-infection immunity in convalescent subjects and non-infected individuals. The emergence of the variants and episodes of COVID-19 infections made the situation worsen. The timely introduction of vaccines and precautionary measures helped control the infection's severity. Later, the population that recovered from COVID-19 grew significantly. However, understanding the impact of healthcare issues resulting after infection is paramount for improving an individual's health status. It is now recognised that COVID-19 infection affects multiple organs and exhibits a broad range of clinical manifestations. So, post COVID-19 infection creates a high risk in individuals with already prevailing health complications. The identification of post-COVID-19-related health issues and their appropriate management is of greater importance to improving patient's quality of life. The persistence, sequelae and other medical complications that normally last from weeks to months after the recovery of the initial infection are involved with COVID-19. A multi-disciplinary approach is necessary for the development of preventive measures, techniques for rehabilitation and strategies for clinical management when it comes to long-term care.

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

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

Unambiguous detection of SARS-CoV-2 subgenomic mRNAs with single-cell RNA sequencing

Single-cell RNA sequencing (scRNA-Seq) studies have provided critical insight into the pathogenesis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19). scRNA-Seq library preparation methods and data processing workflows are generally designed for the detection and quantification of eukaryotic host mRNAs and not viral RNAs. Here, we compare different scRNA-Seq library preparation methods for their ability to quantify and detect SARS-CoV-2 RNAs with a focus on subgenomic mRNAs (sgmRNAs). We show that compared to 10X Genomics Chromium Next GEM Single Cell 3' (10X 3') libraries or 10X Genomics Chromium Next GEM Single Cell V(D)J (10X 5') libraries sequenced with standard read configurations, 10X 5' libraries sequenced with an extended length read 1 (R1) that covers both cell barcode and transcript sequence (termed "10X 5' with extended R1") increase the number of unambiguous reads spanning leader-sgmRNA junction sites. We further present a data processing workflow, single-cell coronavirus sequencing (scCoVseq), which quantifies reads unambiguously assigned to viral sgmRNAs or viral genomic RNA (gRNA). We find that combining 10X 5' with extended R1 library preparation/sequencing and scCoVseq data processing maximizes the number of viral UMIs per cell quantified by scRNA-Seq. Corresponding sgmRNA expression levels are highly correlated with expression in matched bulk RNA-Seq data sets quantified with established tools for SARS-CoV-2 analysis. Using this scRNA-Seq approach, we find that SARS-CoV-2 gene expression is highly correlated across individual infected cells, which suggests that the proportion of viral sgmRNAs remains generally consistent throughout infection. Taken together, these results and corresponding data processing workflow enable robust quantification of coronavirus sgmRNA expression at single-cell resolution, thereby supporting high-resolution studies of viral RNA processes in individual cells. IMPORTANCE Single-cell RNA sequencing (scRNA-Seq) has emerged as a valuable tool to study host-virus interactions, especially for coronavirus disease 2019 (COVID-19). Here we compare the performance of different scRNA-Seq library preparation methods and sequencing strategies to detect SARS-CoV-2 RNAs and develop a data processing workflow to quantify unambiguous sequence reads derived from SARS-CoV-2 genomic RNA and subgenomic mRNAs. After establishing a workflow that maximizes the detection of SARS-CoV-2 subgenomic mRNAs, we explore patterns of SARS-CoV-2 gene expression across cells with variable levels of total viral RNA, assess host gene expression differences between infected and bystander cells, and identify non-canonical and lowly abundant SARS-CoV-2 RNAs. The sequencing and data processing strategies developed here can enhance studies of coronavirus RNA biology at single-cell resolution and thereby contribute to our understanding of viral pathogenesis.

Quantitative SARS-CoV-2 subgenomic RNA as a surrogate marker for viral infectivity: Comparison between culture isolation and direct sgRNA quantification

Detection of subgenomic (sg) SARS-CoV-2 RNAs are frequently used as a correlate of viral infectiousness, but few data about correlation between sg load and viable virus are available. Here, we defined concordance between culture isolation and E and N sgRNA quantification by ddPCR assays in 51 nasopharyngeal swabs collected from SARS-CoV-2 positive hospitalized patients. Among the 51 samples, 14 were SARS-CoV-2 culture-positive and 37 were negative. According to culture results, the sensitivity and specificity of E and N sgRNA assays were 100% and 100%, and 84% and 86%, respectively. ROC analysis showed that the best E and N cut-offs to predict positive culture isolation were 32 and 161 copies/mL respectively, with an AUC (95% CI) of 0.96 (0.91-1.00) and 0.96 (0.92-1.00), and a diagnostic accuracy of 88% and 92%, respectively. Even if no significant correlations were observed between sgRNA amount and clinical presentation, a higher number of moderate/severe cases and lower number of days from symptoms onset characterized patients with sgRNA equal to or higher than sgRNA cut-offs. Overall, this study suggests that SARS-CoV-2 sgRNA quantification could be helpful to estimate the replicative activity of SARS-CoV-2 and can represent a valid surrogate marker to efficiently recognize patients with active infection. The inclusion of this assay in available SARS-CoV-2 diagnostics procedure might help in optimizing fragile patients monitoring and management.

Genome analysis of SARS-CoV-2 isolates from a population reveals the rapid selective sweep of a haplotype carrying many pre-existing and new mutations

To understand the mechanism underlying the evolution of SARS-CoV-2 in a population, we sequenced 92 viral genomes from Assam, India. Analysis of these and database sequences revealed a complete selective sweep of a haplotype in Assam carrying 13 pre-existing variants, including a high leap in frequency of a variant on ORF8, which is involved in immune evasion. A comparative study between sequences of same lineage and similar time frames in and outside Assam showed that 10 of the 13 pre-existing variants had a frequency ranging from 96 to 99%, and the remaining 3 had a low frequency outside Assam. Using a phylogenetic approach to infer sequential occurrences of variants we found that the variant Phe120del on ORF8, which had a low frequency (1.75%) outside Assam, is at the base of the phylogenetic tree of variants and became totally fixed (100%) in Assam population. Based on this observation, we inferred that the variant on ORF8 had a selective advantage, so it carried the haplotype to reach the100% frequency. The haplotype also carried 32 pre-existing variants at a frequency from 1.00 to 80.00% outside Assam. Those of these variants that are more closely linked to the S-protein locus, which often carries advantageous mutations and is tightly linked to the ORF8 locus, retained higher frequencies, while the less tightly linked variants showed lower frequencies, likely due to recombination among co- circulating variants in Assam. The ratios of non-synonymous substitutions to synonymous substitutions suggested that some genes such as those coding for the S-protein and non-structural proteins underwent positive selection while others were subject to purifying selection during their evolution in Assam. Furthermore, we observed negative correlation of the Ct value of qRT-PCR of the patients with abundant ORF6 transcripts, suggesting that ORF6 can be used as a marker for estimating viral titer. In conclusion, our in-depth analysis of SARS-CoV-2 genomes in a regional population reveals the mechanism and dynamics of viral evolution.

A culture-free method for rapidly and accurately quantifying active SARS-CoV-2

Determining the quantity of active virus is the most important basis to judge the risk of virus infection, which usually relies on the virus median tissue culture infectious dose (TCID50) assay performed in a biosafety level 3 laboratory within 5-7 days. We have developed a culture-free method for rapid and accurate quantification of active severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by targeting subgenomic RNA (sgRNA) based on reverse transcription digital PCR (RT-dPCR). The dynamic range of quantitative assays for sgRNA-N and sgRNA-E by RT-dPCR was investigated, and the result showed that the limits of detection (LoD) and quantification (LoQ) were 2 copies/reaction and 10 copies/reaction, respectively. The delta strain (NMDC60042793) of SARS-CoV-2 was cultured at an average titer of 106.13 TCID50/mL and used to evaluate the developed quantification method. Copy number concentrations of the cultured SARS-CoV-2 sgRNA and genomic RNA (gRNA) gave excellent linearity (R2 = 0.9999) with SARS-CoV-2 titers in the range from 500 to 105 TCID50/mL. Validation of 63 positive clinical samples further proves that the quantification of sgRNA-N by RT-dPCR is more sensitive for active virus quantitative detection. It is notable that we can infer the active virus titer through quantification of SARS-CoV-2 sgRNA based on the linear relationship in a biosafety level 2 laboratory within 3 h. It can be used to timely and effectively identify infectious patients and reduce unnecessary isolation especially when a large number of COVID-19 infected people impose a burden on medical resources.

High-sensitivity profiling of SARS-CoV-2 noncoding region-host protein interactome reveals the potential regulatory role of negative-sense viral RNA

A deep understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-host interactions is crucial to developing effective therapeutics and addressing the threat of emerging coronaviruses. The role of noncoding regions of viral RNA (ncrRNAs) has yet to be systematically scrutinized. We developed a method using MS2 affinity purification coupled with liquid chromatography-mass spectrometry and designed a diverse set of bait ncrRNAs to systematically map the interactome of SARS-CoV-2 ncrRNA in Calu-3, Huh7, and HEK293T cells. Integration of the results defined the core ncrRNA-host protein interactomes among cell lines. The 5' UTR interactome is enriched with proteins in the small nuclear ribonucleoproteins family and is a target for the regulation of viral replication and transcription. The 3' UTR interactome is enriched with proteins involved in the stress granules and heterogeneous nuclear ribonucleoproteins family. Intriguingly, compared with the positive-sense ncrRNAs, the negative-sense ncrRNAs, especially the negative-sense of 3' UTR, interacted with a large array of host proteins across all cell lines. These proteins are involved in the regulation of the viral production process, host cell apoptosis, and immune response. Taken together, our study depicts the comprehensive landscape of the SARS-CoV-2 ncrRNA-host protein interactome and unveils the potential regulatory role of the negative-sense ncrRNAs, providing a new perspective on virus-host interactions and the design of future therapeutics. Given the highly conserved nature of UTRs in positive-strand viruses, the regulatory role of negative-sense ncrRNAs should not be exclusive to SARS-CoV-2. IMPORTANCE Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes COVID-19, a pandemic affecting millions of lives. During replication and transcription, noncoding regions of the viral RNA (ncrRNAs) may play an important role in the virus-host interactions. Understanding which and how these ncrRNAs interact with host proteins is crucial for understanding the mechanism of SARS-CoV-2 pathogenesis. We developed the MS2 affinity purification coupled with liquid chromatography-mass spectrometry method and designed a diverse set of ncrRNAs to identify the SARS-CoV-2 ncrRNA interactome comprehensively in different cell lines and found that the 5' UTR binds to proteins involved in U1 small nuclear ribonucleoprotein, while the 3' UTR interacts with proteins involved in stress granules and the heterogeneous nuclear ribonucleoprotein family. Interestingly, negative-sense ncrRNAs showed interactions with a large number of diverse host proteins, indicating a crucial role in infection. The results demonstrate that ncrRNAs could serve diverse regulatory functions.

A mutation in the coronavirus nsp13-helicase impairs enzymatic activity and confers partial remdesivir resistance

Coronaviruses (CoVs) encode nonstructural proteins 1-16 (nsps 1-16) which form replicase complexes that mediate viral RNA synthesis. Remdesivir (RDV) is an adenosine nucleoside analog antiviral that inhibits CoV RNA synthesis. RDV resistance mutations have been reported only in the nonstructural protein 12 RNA-dependent RNA polymerase (nsp12-RdRp). We here show that a substitution mutation in the nsp13-helicase (nsp13-HEL A335V) of the betacoronavirus murine hepatitis virus (MHV) that was selected during passage with the RDV parent compound confers partial RDV resistance independently and additively when expressed with co-selected RDV resistance mutations in the nsp12-RdRp. The MHV A335V substitution did not enhance replication or competitive fitness compared to WT MHV and remained sensitive to the active form of the cytidine nucleoside analog antiviral molnupiravir (MOV). Biochemical analysis of the SARS-CoV-2 helicase encoding the homologous substitution (A336V) demonstrates that the mutant protein retained the ability to associate with the core replication proteins nsps 7, 8, and 12 but had impaired helicase unwinding and ATPase activity. Together, these data identify a novel determinant of nsp13-HEL enzymatic activity, define a new genetic pathway for RDV resistance, and demonstrate the importance of surveillance for and testing of helicase mutations that arise in SARS-CoV-2 genomes. IMPORTANCE Despite the development of effective vaccines against COVID-19, the continued circulation and emergence of new variants support the need for antivirals such as RDV. Understanding pathways of antiviral resistance is essential for surveillance of emerging variants, development of combination therapies, and for identifying potential new targets for viral inhibition. We here show a novel RDV resistance mutation in the CoV helicase also impairs helicase functions, supporting the importance of studying the individual and cooperative functions of the replicase nonstructural proteins 7-16 during CoV RNA synthesis. The homologous nsp13-HEL mutation (A336V) has been reported in the GISAID database of SARS-CoV-2 genomes, highlighting the importance of surveillance of and genetic testing for nucleoside analog resistance in the helicase.

SARS-CoV-2 infection of human pluripotent stem cell-derived vascular cells reveals smooth muscle cells as key mediators of vascular pathology during infection

Although respiratory symptoms are the most prevalent disease manifestation of infection by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), nearly 20% of hospitalized patients are at risk for thromboembolic events 1 . This prothrombotic state is considered a key factor in the increased risk of stroke, which has been observed clinically during both acute infection and long after symptoms have cleared 2 . Here we developed a model of SARS-CoV-2 infection using human-induced pluripotent stem cell-derived endothelial cells, pericytes, and smooth muscle cells to recapitulate the vascular pathology associated with SARS-CoV-2 exposure. Our results demonstrate that perivascular cells, particularly smooth muscle cells (SMCs), are a specifically susceptible vascular target for SARS-CoV-2 infection. Utilizing RNA sequencing, we characterized the transcriptomic changes accompanying SARS-CoV-2 infection of SMCs, and endothelial cells (ECs). We observed that infected human SMCs shift to a pro-inflammatory state and increase the expression of key mediators of the coagulation cascade. Further, we showed human ECs exposed to the secretome of infected SMCs produce hemostatic factors that can contribute to vascular dysfunction, despite not being susceptible to direct infection. The findings here recapitulate observations from patient sera in human COVID-19 patients and provide mechanistic insight into the unique vascular implications of SARS-CoV-2 infection at a cellular level.

Interspecies transmission of SARS CoV-2 with special emphasis on viral mutations and ACE-2 receptor homology roles

COVID-19 outbreak was first reported in 2019, Wuhan, China. The spillover of the disease caused by Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), to a wide range of pet, zoo, wild, and farm animals has emphasized potential zoonotic and reverse zoonotic viral transmission. Furthermore, it has evoked inquiries about susceptibility of different animal species to SARS-CoV-2 infection and role of these animals as viral reservoirs. Therefore, studying susceptible and non-susceptible hosts for SARS-CoV-2 infection could give a better understanding for the virus and will help in preventing further outbreaks. Here, we review structural aspects of SARS-CoV-2 spike protein, the effect of the different mutations observed in the spike protein, and the impact of ACE2 receptor variations in different animal hosts on inter-species transmission. Moreover, the SARS-CoV-2 spillover chain was reviewed. Combination of SARS-CoV-2 high mutation rate and homology of cellular ACE2 receptors enable the virus to transcend species barriers and facilitate its transmission between humans and animals.

Advanced fluorescence microscopy in respiratory virus cell biology

Respiratory viruses are a major public health burden across all age groups around the globe, and are associated with high morbidity and mortality rates. They can be transmitted by multiple routes, including physical contact or droplets and aerosols, resulting in efficient spreading within the human population. Investigations of the cell biology of virus replication are thus of utmost importance to gain a better understanding of virus-induced pathogenicity and the development of antiviral countermeasures. Light and fluorescence microscopy techniques have revolutionized investigations of the cell biology of virus infection by allowing the study of the localization and dynamics of viral or cellular components directly in infected cells. Advanced microscopy including high- and super-resolution microscopy techniques available today can visualize biological processes at the single-virus and even single-molecule level, thus opening a unique view on virus infection. We will highlight how fluorescence microscopy has supported investigations on virus cell biology by focusing on three major respiratory viruses: respiratory syncytial virus (RSV), Influenza A virus (IAV) and SARS-CoV-2. We will review our current knowledge of virus replication and highlight how fluorescence microscopy has helped to improve our state of understanding. We will start by introducing major imaging and labeling modalities and conclude the chapter with a perspective discussion on remaining challenges and potential opportunities.

Comparative subgenomic mRNA profiles of SARS-CoV-2 Alpha, Delta and Omicron BA.1, BA.2 and BA.5 sub-lineages using Danish COVID-19 genomic surveillance data

BACKGROUND

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has rapidly spread worldwide in the population since it was first detected in late 2019. The transcription and replication of coronaviruses, although not fully understood, is characterised by the production of genomic length RNA and shorter subgenomic RNAs to make viral proteins and ultimately progeny virions. Observed levels of subgenomic RNAs differ between sub-lineages and open reading frames but their biological significance is presently unclear.

METHODS

Using a large and diverse panel of virus sequencing data produced as part of the Danish COVID-19 routine surveillance together with information in electronic health registries, we assessed the association of subgenomic RNA levels with demographic and clinical variables of the infected individuals.

FINDINGS

Our findings suggest no significant statistical relationship between levels of subgenomic RNAs and host-related factors.

INTERPRETATION

Differences between lineages and subgenomic ORFs may be related to differences in target cell tropism, early virus replication/transcription kinetics or sequence features.

FUNDING

The author(s) received no specific funding for this work.

Porcine epidemic diarrhea virus: an update overview of virus epidemiology, vaccines, and control strategies in South Korea

Porcine epidemic diarrhea virus (PEDV) has posed significant financial threats to the domestic pig industry over the last three decades in South Korea. PEDV infection will mostly result in endemic persistence in the affected farrow-to-finish (FTF) herds, leading to endemic porcine epidemic diarrhea (PED) followed by year-round recurrent outbreaks. This review aims to encourage collaboration among swine producers, veterinarians, and researchers to offer answers that strengthen our understanding of PEDV in efforts to prevent and control endemic PED and to prepare for the next epidemics or pandemics. We found that collaboratively implementing a PED risk assessment and customized four-pillar-based control measures is vital to interrupt the chain of endemic PED in affected herds: the former can identify on-farm risk factors while the latter aims to compensate for or improve weaknesses via herd immunity stabilization and virus elimination. Under endemic PED, long-term virus survival in slurry and asymptomatically infected gilts ("Trojan Pigs") that can transmit the virus to farrowing houses are key challenges for PEDV eradication in FTF farms and highlight the necessity for active monitoring and surveillance of the virus in herds and their environments. This paper underlines the current knowledge of molecular epidemiology and commercially available vaccines, as well as the risk assessment and customized strategies to control PEDV. The intervention measures for stabilizing herd immunity and eliminating virus circulation may be the cornerstone of establishing regional or national PED eradication programs.