RNA structure analysis of alphacoronavirus terminal genome regions

Alpha- and betacoronavirus cis-acting RNA elements

Coronaviruses have exceptionally large RNA genomes and employ multiprotein replication/transcription complexes to orchestrate specific steps of viral RNA genome replication and expression. Most of these processes involve viral cis-acting RNA elements that are engaged in vital RNA-RNA and/or RNA-protein interactions. Over the past years, a large number of studies provided interesting new insight into the structures and, to a lesser extent, functions of specific RNA elements for representative coronaviruses, and there is evidence to suggest that (a majority of) these RNA elements are conserved across genetically divergent coronavirus genera. It is becoming increasingly clear that at least some of these elements do not function in isolation but operate through complex and highly dynamic RNA-RNA interactions. This article reviews structural and functional aspects of cis-acting RNA elements conserved in alpha- and betacoronavirus 5'- and 3'-terminal genome regions, focusing on their critical roles in viral RNA synthesis and gene expression.

Tertiary folds of the SL5 RNA from the 5' proximal region of SARS-CoV-2 and related coronaviruses

Coronavirus genomes sequester their start codons within stem-loop 5 (SL5), a structured, 5' genomic RNA element. In most alpha- and betacoronaviruses, the secondary structure of SL5 is predicted to contain a four-way junction of helical stems, some of which are capped with UUYYGU hexaloops. Here, using cryogenic electron microscopy (cryo-EM) and computational modeling with biochemically determined secondary structures, we present three-dimensional structures of SL5 from six coronaviruses. The SL5 domain of betacoronavirus severe-acute-respiratory-syndrome-related coronavirus 2 (SARS-CoV-2), resolved at 4.7 Å resolution, exhibits a T-shaped structure, with its UUYYGU hexaloops at opposing ends of a coaxial stack, the T's "arms." Further analysis of SL5 domains from SARS-CoV-1 and MERS (7.1 and 6.4 to 6.9 Å resolution, respectively) indicate that the junction geometry and inter-hexaloop distances are conserved features across these human-infecting betacoronaviruses. The MERS SL5 domain displays an additional tertiary interaction, which is also observed in the non-human-infecting betacoronavirus BtCoV-HKU5 (5.9 to 8.0 Å resolution). SL5s from human-infecting alphacoronaviruses, HCoV-229E and HCoV-NL63 (6.5 and 8.4 to 9.0 Å resolution, respectively), exhibit the same coaxial stacks, including the UUYYGU-capped arms, but with a phylogenetically distinct crossing angle, an X-shape. As such, all SL5 domains studied herein fold into stable tertiary structures with cross-genus similarities and notable differences, with implications for potential protein-binding modes and therapeutic targets.

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.

Tertiary folds of the SL5 RNA from the 5' proximal region of SARS-CoV-2 and related coronaviruses

Coronavirus genomes sequester their start codons within stem-loop 5 (SL5), a structured, 5' genomic RNA element. In most alpha- and betacoronaviruses, the secondary structure of SL5 is predicted to contain a four-way junction of helical stems, some of which are capped with UUYYGU hexaloops. Here, using cryogenic electron microscopy (cryo-EM) and computational modeling with biochemically-determined secondary structures, we present three-dimensional structures of SL5 from six coronaviruses. The SL5 domain of betacoronavirus SARS-CoV-2, resolved at 4.7 Å resolution, exhibits a T-shaped structure, with its UUYYGU hexaloops at opposing ends of a coaxial stack, the T's "arms." Further analysis of SL5 domains from SARS-CoV-1 and MERS (7.1 and 6.4-6.9 Å resolution, respectively) indicate that the junction geometry and inter-hexaloop distances are conserved features across the studied human-infecting betacoronaviruses. The MERS SL5 domain displays an additional tertiary interaction, which is also observed in the non-human-infecting betacoronavirus BtCoV-HKU5 (5.9-8.0 Å resolution). SL5s from human-infecting alphacoronaviruses, HCoV-229E and HCoV-NL63 (6.5 and 8.4-9.0 Å resolution, respectively), exhibit the same coaxial stacks, including the UUYYGU-capped arms, but with a phylogenetically distinct crossing angle, an X-shape. As such, all SL5 domains studied herein fold into stable tertiary structures with cross-genus similarities, with implications for potential protein-binding modes and therapeutic targets.

Targeting the conserved coronavirus octamer motif GGAAGAGC is a strategy for the development of coronavirus vaccine

BACKGROUND

Coronaviruses are pathogens of humans and animals that cause widespread and costly diseases. The development of effective strategies to combat the threat of coronaviruses is therefore a top priority. The conserved coronavirus octamer motif 5'GGAAGAGC3' exists in the 3' untranslated region of all identified coronaviruses. In the current study, we aimed to examine whether targeting the coronavirus octamer motif GGAAGAGC is a promising approach to develop coronavirus vaccine.

METHODS

Plaque assays were used to determine the titers of mouse hepatitis virus (MHV)-A59 octamer mutant (MHVoctm) and wild-type (wt) MHV-A59 (MHVwt). Western blotting was used for the determination of translation efficiency of MHVoctm and MHVwt. Plaque assays and RT-qPCR were employed to examine whether MHVoctm was more sensitive to interferon treatment than MHVwt. Weight loss, clinical signs, survival rate, viral RNA detection and histopathological examination were used to evaluate whether MHVoctm was a vaccine candidate against MHVwt infection in BALB/c mice.

RESULTS

In this study, we showed that (i) the MHVoctm with mutation of coronavirus octamer was able to grow to high titers but attenuated in mice, (ii) with the reduced multiplicity of infection (MOI), the difference in gene expression between MHVoctm and MHVwt became more evident in cultured cells, (iii) MHVoctm was more sensitive to interferon treatment than MHVwt and (iv) mice inoculated with MHVoctm were protected from MHVwt infection.

CONCLUSIONS

Based on the results obtained from cultured cells, it was suggested that the synergistic effects of octamer mutation, multiplicity of infection and immune response may be a mechanism explaining the distinct phenotypes of octamer-mutated coronavirus in cell culture and mice. In addition, targeting the conserved coronavirus octamer motif is a strategy for development of coronavirus vaccine. Since the conserved octamer exists in all coronaviruses, this strategy of targeting the conserved octamer motif can also be applied to other human and animal coronaviruses for the development of coronavirus vaccines, especially the emergence of novel coronaviruses such as SARS-CoV-2, saving time and cost for vaccine development and disease control.

A viral pan-end RNA element and host complex define a SARS-CoV-2 regulon

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, generates multiple protein-coding, subgenomic RNAs (sgRNAs) from a longer genomic RNA, all bearing identical termini with poorly understood roles in regulating viral gene expression. Insulin and interferon-gamma, two host-derived, stress-related agents, and virus spike protein, induce binding of glutamyl-prolyl-tRNA synthetase (EPRS1), within an unconventional, tetra-aminoacyl-tRNA synthetase complex, to the sgRNA 3'-end thereby enhancing sgRNA expression. We identify an EPRS1-binding sarbecoviral pan-end activating RNA (SPEAR) element in the 3'-end of viral RNAs driving agonist-induction. Translation of another co-terminal 3'-end feature, ORF10, is necessary for SPEAR-mediated induction, independent of Orf10 protein expression. The SPEAR element enhances viral programmed ribosomal frameshifting, thereby expanding its functionality. By co-opting noncanonical activities of a family of essential host proteins, the virus establishes a post-transcriptional regulon stimulating global viral RNA translation. A SPEAR-targeting strategy markedly reduces SARS-CoV-2 titer, suggesting a pan-sarbecoviral therapeutic modality.

The preference signature of the SARS-CoV-2 Nucleocapsid NTD for its 5'-genomic RNA elements

The nucleocapsid protein (N) of SARS-CoV-2 plays a pivotal role during the viral life cycle. It is involved in RNA transcription and accounts for packaging of the large genome into virus particles. N manages the enigmatic balance of bulk RNA-coating versus precise RNA-binding to designated cis-regulatory elements. Numerous studies report the involvement of its disordered segments in non-selective RNA-recognition, but how N organizes the inevitable recognition of specific motifs remains unanswered. We here use NMR spectroscopy to systematically analyze the interactions of N's N-terminal RNA-binding domain (NTD) with individual cis RNA elements clustering in the SARS-CoV-2 regulatory 5'-genomic end. Supported by broad solution-based biophysical data, we unravel the NTD RNA-binding preferences in the natural genome context. We show that the domain's flexible regions read the intrinsic signature of preferred RNA elements for selective and stable complex formation within the large pool of available motifs.

Magnipore: Prediction of differential single nucleotide changes in the Oxford Nanopore Technologies sequencing signal of SARS-CoV-2 samples

Oxford Nanopore Technologies (ONT) allows direct sequencing of ribonucleic acids (RNA) and, in addition, detection of possible RNA modifications due to deviations from the expected ONT signal. The software available so far for this purpose can only detect a small number of modifications. Alternatively, two samples can be compared for different RNA modifications. We present Magnipore, a novel tool to search for significant signal shifts between samples of Oxford Nanopore data from similar or related species. Magnipore classifies them into mutations and potential modifications. We use Magnipore to compare SARS-CoV-2 samples. Included were representatives of the early 2020s Pango lineages (n=6), samples from Pango lineages B.1.1.7 (n=2, Alpha), B.1.617.2 (n=1, Delta), and B.1.529 (n=7, Omicron). Magnipore utilizes position-wise Gaussian distribution models and a comprehensible significance threshold to find differential signals. In the case of Alpha and Delta, Magnipore identifies 55 detected mutations and 15 sites that hint at differential modifications. We predicted potential virus-variant and variant-group-specific differential modifications. Magnipore contributes to advancing RNA modification analysis in the context of viruses and virus variants.

Stimulation of interferon-β responses by aberrant SARS-CoV-2 small viral RNAs acting as retinoic acid-inducible gene-I agonists

Patients with severe COVID-19 exhibit a cytokine storm characterized by greatly elevated levels of cytokines. Despite this, the interferon (IFN) response is delayed, contributing to disease progression. Here, we report that SARS-CoV-2 excessively generates small viral RNAs (svRNAs) encoding exact 5' ends of positive-sense genes in human cells in vitro and ex vivo, whereas endemic human coronaviruses (OC43 and 229E) produce significantly fewer similar svRNAs. SARS-CoV-2 5' end svRNAs are RIG-I agonists and induce the IFN-β response in the later stages of infection. The first 60-nt ends bearing duplex structures and 5'-triphosphates are responsible for immune-stimulation. We propose that RIG-I activation by accumulated SARS-CoV-2 5' end svRNAs may contribute to later drive over-exuberant IFN production. Additionally, the differences in the amounts of svRNAs produced and the corresponding IFN response among CoV strains suggest that lower svRNA production during replication may correlate with the weaker immune response seen in less pathogenic CoVs.

Genome Structure, Life Cycle, and Taxonomy of Coronaviruses and the Evolution of SARS-CoV-2

Coronaviruses have a broad host range and exhibit high zoonotic potential. In this chapter, we describe their genomic organization in terms of encoded proteins and provide an introduction to the peculiar discontinuous transcription mechanism. Further, we present evolutionary conserved genomic RNA secondary structure features, which are involved in the complex replication mechanism. With a focus on computational methods, we review the emergence of SARS-CoV-2 starting with the 2019 strains. In that context, we also discuss the debated hypothesis of whether SARS-CoV-2 was created in a laboratory. We focus on the molecular evolution and the epidemiological dynamics of this recently emerged pathogen and we explain how variants of concern are detected and characterised. COVID-19, the disease caused by SARS-CoV-2, can spread through different transmission routes and also depends on a number of risk factors. We describe how current computational models of viral epidemiology, or more specifically, phylodynamics, have facilitated and will continue to enable a better understanding of the epidemic dynamics of SARS-CoV-2.

Women in the European Virus Bioinformatics Center

Viruses are the cause of a considerable burden to human, animal and plant health, while on the other hand playing an important role in regulating entire ecosystems. The power of new sequencing technologies combined with new tools for processing "Big Data" offers unprecedented opportunities to answer fundamental questions in virology. Virologists have an urgent need for virus-specific bioinformatics tools. These developments have led to the formation of the European Virus Bioinformatics Center, a network of experts in virology and bioinformatics who are joining forces to enable extensive exchange and collaboration between these research areas. The EVBC strives to provide talented researchers with a supportive environment free of gender bias, but the gender gap in science, especially in math-intensive fields such as computer science, persists. To bring more talented women into research and keep them there, we need to highlight role models to spark their interest, and we need to ensure that female scientists are not kept at lower levels but are given the opportunity to lead the field. Here we showcase the work of the EVBC and highlight the achievements of some outstanding women experts in virology and viral bioinformatics.

1H, 13C and 15N chemical shift assignment of the stem-loops 5b + c from the 5'-UTR of SARS-CoV-2

The ongoing pandemic of the respiratory disease COVID-19 is caused by the SARS-CoV-2 (SCoV2) virus. SCoV2 is a member of the Betacoronavirus genus. The 30 kb positive sense, single stranded RNA genome of SCoV2 features 5'- and 3'-genomic ends that are highly conserved among Betacoronaviruses. These genomic ends contain structured cis-acting RNA elements, which are involved in the regulation of viral replication and translation. Structural information about these potential antiviral drug targets supports the development of novel classes of therapeutics against COVID-19. The highly conserved branched stem-loop 5 (SL5) found within the 5'-untranslated region (5'-UTR) consists of a basal stem and three stem-loops, namely SL5a, SL5b and SL5c. Both, SL5a and SL5b feature a 5'-UUUCGU-3' hexaloop that is also found among Alphacoronaviruses. Here, we report the extensive ¹H, ¹³C and ¹⁵N resonance assignment of the 37 nucleotides (nts) long sequence spanning SL5b and SL5c (SL5b + c), as basis for further in-depth structural studies by solution NMR spectroscopy.

Automated Flow Synthesis of Peptide-PNA Conjugates

Antisense peptide nucleic acids (PNAs) have yet to translate to the clinic because of poor cellular uptake, limited solubility, and rapid elimination. Cell-penetrating peptides (CPPs) covalently attached to PNAs may facilitate clinical development by improving uptake into cells. We report an efficient technology that utilizes a fully automated fast-flow instrument to manufacture CPP-conjugated PNAs (PPNAs) in a single shot. The machine is rapid, with each amide bond being formed in 10 s. Anti-IVS2-654 PPNA synthesized with this instrument presented threefold activity compared to transfected PNA in a splice-correction assay. We demonstrated the utility of this approach by chemically synthesizing eight anti-SARS-CoV-2 PPNAs in 1 day. A PPNA targeting the 5' untranslated region of SARS-CoV-2 genomic RNA reduced the viral titer by over 95% in a live virus infection assay (IC₅₀ = 0.8 μM). Our technology can deliver PPNA candidates to further investigate their potential as antiviral agents.

Duplex formation between the template and the nascent strand in the transcription-regulating sequences is associated with the site of template switching in SARS - CoV-2

Recently published transcriptomic data of the SARS-CoV-2 coronavirus show that there is a large variation in the frequency and steady state levels of subgenomic mRNA sequences. This variation is derived from discontinuous subgenomic RNA synthesis, where the polymerase switches template from a 3' proximal genome body sequence to a 5' untranslated leader sequence. This leads to a fusion between the common 5' leader sequence and a 3' proximal body sequence in the RNA product. This process revolves around a common core sequence (CS) that is present at both the template sites that make up the fusion junction. Base-pairing between the leader CS and the nascent complementary minus strand body CS, and flanking regions (together called the transcription regulating sequence, TRS) is vital for this template switching event. However, various factors can influence the site of template switching within the same TRS duplex. Here, we model the duplexes formed between the leader and complementary body TRS regions, hypothesizing the role of the stability of the TRS duplex in determining the major sites of template switching for the most abundant mRNAs. We indicate that the stability of secondary structures and the speed of transcription play key roles in determining the probability of template switching in the production of subgenomic RNAs. We speculate on the effect of reported variant nucleotide substitutions on our models.

SARS-CoV-2: from its discovery to genome structure, transcription, and replication

SARS-CoV-2 is an extremely contagious respiratory virus causing adult atypical pneumonia COVID-19 with severe acute respiratory syndrome (SARS). SARS-CoV-2 has a single-stranded, positive-sense RNA (+RNA) genome of ~ 29.9 kb and exhibits significant genetic shift from different isolates. After entering the susceptible cells expressing both ACE2 and TMPRSS2, the SARS-CoV-2 genome directly functions as an mRNA to translate two polyproteins from the ORF1a and ORF1b region, which are cleaved by two viral proteases into sixteen non-structural proteins (nsp1-16) to initiate viral genome replication and transcription. The SARS-CoV-2 genome also encodes four structural (S, E, M and N) and up to six accessory (3a, 6, 7a, 7b, 8, and 9b) proteins, but their translation requires newly synthesized individual subgenomic RNAs (sgRNA) in the infected cells. Synthesis of the full-length viral genomic RNA (gRNA) and sgRNAs are conducted inside double-membrane vesicles (DMVs) by the viral replication and transcription complex (RTC), which comprises nsp7, nsp8, nsp9, nsp12, nsp13 and a short RNA primer. To produce sgRNAs, RTC starts RNA synthesis from the highly structured gRNA 3' end and switches template at various transcription regulatory sequence (TRSB) sites along the gRNA body probably mediated by a long-distance RNA-RNA interaction. The TRS motif in the gRNA 5' leader (TRSL) is responsible for the RNA-RNA interaction with the TRSB upstream of each ORF and skipping of the viral genome in between them to produce individual sgRNAs. Abundance of individual sgRNAs and viral gRNA synthesized in the infected cells depend on the location and read-through efficiency of each TRSB. Although more studies are needed, the unprecedented COVID-19 pandemic has taught the world a painful lesson that is to invest and proactively prepare future emergence of other types of coronaviruses and any other possible biological horrors.

The SARS-CoV-2 subgenome landscape and its novel regulatory features

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is currently a global pandemic. CoVs are known to generate negative subgenomes (subgenomic RNAs [sgRNAs]) through transcription-regulating sequence (TRS)-dependent template switching, but the global dynamic landscapes of coronaviral subgenomes and regulatory rules remain unclear. Here, using next-generation sequencing (NGS) short-read and Nanopore long-read poly(A) RNA sequencing in two cell types at multiple time points after infection with SARS-CoV-2, we identified hundreds of template switches and constructed the dynamic landscapes of SARS-CoV-2 subgenomes. Interestingly, template switching could occur in a bidirectional manner, with diverse SARS-CoV-2 subgenomes generated from successive template-switching events. The majority of template switches result from RNA-RNA interactions, including seed and compensatory modes, with terminal pairing status as a key determinant. Two TRS-independent template switch modes are also responsible for subgenome biogenesis. Our findings reveal the subgenome landscape of SARS-CoV-2 and its regulatory features, providing a molecular basis for understanding subgenome biogenesis and developing novel anti-viral strategies.

Analysis of Emerging Variants in Structured Regions of the SARS-CoV-2 Genome

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has motivated a widespread effort to understand its epidemiology and pathogenic mechanisms. Modern high-throughput sequencing technology has led to the deposition of vast numbers of SARS-CoV-2 genome sequences in curated repositories, which have been useful in mapping the spread of the virus around the globe. They also provide a unique opportunity to observe virus evolution in real time. Here, we evaluate two sets of SARS-CoV-2 genomic sequences to identify emerging variants within structured cis-regulatory elements of the SARS-CoV-2 genome. Overall, 20 variants are present at a minor allele frequency of at least 0.5%. Several enhance the stability of Stem Loop 1 in the 5' untranslated region (UTR), including a group of co-occurring variants that extend its length. One appears to modulate the stability of the frameshifting pseudoknot between ORF1a and ORF1b, and another perturbs a bi-ss molecular switch in the 3'UTR. Finally, 5 variants destabilize structured elements within the 3'UTR hypervariable region, including the S2M (stem loop 2 m) selfish genetic element, raising questions as to the functional relevance of these structures in viral replication. Two of the most abundant variants appear to be caused by RNA editing, suggesting host-viral defense contributes to SARS-CoV-2 genome heterogeneity. Our analysis has implications for the development of therapeutics that target viral cis-regulatory RNA structures or sequences.

Structure and regulation of coronavirus genomes: state-of-the-art and novel insights from SARS-CoV-2 studies

Coronaviruses (CoV) are positive-sense single-stranded RNA viruses, harboring the largest viral RNA genomes known to date. Apart from the primary sequence encoding for all the viral proteins needed for the generation of new viral particles, certain regions of CoV genomes are known to fold into stable structures, controlling several aspects of CoV life cycle, from the regulation of the discontinuous transcription of subgenomic mRNAs, to the packaging of the genome into new virions. Here we review the current knowledge on CoV RNA structures, discussing it in light of the most recent discoveries made possible by analyses of the SARS-CoV-2 genome.

Global Geographic and Temporal Analysis of SARS-CoV-2 Haplotypes Normalized by COVID-19 Cases During the Pandemic

Since the identification of SARS-CoV-2, a large number of genomes have been sequenced with unprecedented speed around the world. This marks a unique opportunity to analyze virus spreading and evolution in a worldwide context. Currently, there is not a useful haplotype description to help to track important and globally scattered mutations. Also, differences in the number of sequenced genomes between countries and/or months make it difficult to identify the emergence of haplotypes in regions where few genomes are sequenced but a large number of cases are reported. We propose an approach based on the normalization by COVID-19 cases of relative frequencies of mutations using all the available data to identify major haplotypes. Furthermore, we can use a similar normalization approach to tracking the temporal and geographic distribution of haplotypes in the world. Using 171,461 genomes, we identify five major haplotypes or operational taxonomic units (OTUs) based on nine high-frequency mutations. OTU_3 characterized by mutations R203K and G204R is currently the most frequent haplotype circulating in four of the six continents analyzed (South America, North America, Europe, Asia, Africa, and Oceania). On the other hand, during almost all months analyzed, OTU_5 characterized by the mutation T85I in nsp2 is the most frequent in North America. Recently (since September), OTU_2 has been established as the most frequent in Europe. OTU_1, the ancestor haplotype, is near to extinction showed by its low number of isolations since May. Also, we analyzed whether age, gender, or patient status is more related to a specific OTU. We did not find OTU's preference for any age group, gender, or patient status. Finally, we discuss structural and functional hypotheses in the most frequently identified mutations, none of those mutations show a clear effect on the transmissibility or pathogenicity.

Inhibitors of Eukaryotic Translational Machinery as Therapeutic Agents

Inhibiting eukaryotic protein translation with small molecules is emerging as a powerful therapeutic strategy. The advantage of targeting cellular translational machinery is that it is required for the highly proliferative state of many neoplastic cells, replication of certain viruses, and ultimately the expression of a wide variety of protein targets. Although, this approach has been exploited to develop clinical agents, such as homoharringtonine (HHT, 1), used to treat chronic myeloid leukemia (CML), inhibiting components of the translational machinery is often associated with cytotoxic phenotypes. However, recent studies have demonstrated that certain small molecules can inhibit the translation of specific subsets of proteins, leading to lower cytotoxicity, and opening-up therapeutic opportunities for translation inhibitors to be deployed in indications beyond oncology and infectious disease. This review summarizes efforts to develop inhibitors of the eukaryotic translational machinery as therapeutic agents and highlights emerging opportunities for translation inhibitors in the future.

Rfam 14: expanded coverage of metagenomic, viral and microRNA families

Rfam is a database of RNA families where each of the 3444 families is represented by a multiple sequence alignment of known RNA sequences and a covariance model that can be used to search for additional members of the family. Recent developments have involved expert collaborations to improve the quality and coverage of Rfam data, focusing on microRNAs, viral and bacterial RNAs. We have completed the first phase of synchronising microRNA families in Rfam and miRBase, creating 356 new Rfam families and updating 40. We established a procedure for comprehensive annotation of viral RNA families starting with Flavivirus and Coronaviridae RNAs. We have also increased the coverage of bacterial and metagenome-based RNA families from the ZWD database. These developments have enabled a significant growth of the database, with the addition of 759 new families in Rfam 14. To facilitate further community contribution to Rfam, expert users are now able to build and submit new families using the newly developed Rfam Cloud family curation system. New Rfam website features include a new sequence similarity search powered by RNAcentral, as well as search and visualisation of families with pseudoknots. Rfam is freely available at https://rfam.org.

RNA genome conservation and secondary structure in SARS-CoV-2 and SARS-related viruses: a first look

As the COVID-19 outbreak spreads, there is a growing need for a compilation of conserved RNA genome regions in the SARS-CoV-2 virus along with their structural propensities to guide development of antivirals and diagnostics. Here we present a first look at RNA sequence conservation and structural propensities in the SARS-CoV-2 genome. Using sequence alignments spanning a range of betacoronaviruses, we rank genomic regions by RNA sequence conservation, identifying 79 regions of length at least 15 nt as exactly conserved over SARS-related complete genome sequences available near the beginning of the COVID-19 outbreak. We then confirm the conservation of the majority of these genome regions across 739 SARS-CoV-2 sequences subsequently reported from the COVID-19 outbreak, and we present a curated list of 30 "SARS-related-conserved" regions. We find that known RNA structured elements curated as Rfam families and in prior literature are enriched in these conserved genome regions, and we predict additional conserved, stable secondary structures across the viral genome. We provide 106 "SARS-CoV-2-conserved-structured" regions as potential targets for antivirals that bind to structured RNA. We further provide detailed secondary structure models for the extended 5' UTR, frameshifting stimulation element, and 3' UTR. Lastly, we predict regions of the SARS-CoV-2 viral genome that have low propensity for RNA secondary structure and are conserved within SARS-CoV-2 strains. These 59 "SARS-CoV-2-conserved-unstructured" genomic regions may be most easily accessible by hybridization in primer-based diagnostic strategies.

In silico identification of conserved cis-acting RNA elements in the SARS-CoV-2 genome

Aim:

The aim of this study was to computationally predict conserved RNA sequences and structures known as cis-acting RNA elements (CREs) in the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome.

Materials & methods:

Bioinformatics tools were used to analyze and predict CREs by obtaining viral sequences from available databases.

Results:

Computational analysis revealed the presence of RNA stem-loop structures within the 3′ end of the ORF1ab region analogous to previously identified SARS-CoV genomic packaging signals. Alignment-based RNA secondary structure predictions of the 5′ end of the SARS-CoV-2 genome also identified conserved CREs.

Conclusion:

These CREs may be potential vaccine and/or antiviral therapeutic targets; however, further studies are warranted to confirm their roles in the SARS-CoV-2 life cycle.

Analysis of Rapidly Emerging Variants in Structured Regions of the SARS-CoV-2 Genome

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has motivated a widespread effort to understand its epidemiology and pathogenic mechanisms. Modern high-throughput sequencing technology has led to the deposition of vast numbers of SARS-CoV-2 genome sequences in curated repositories, which have been useful in mapping the spread of the virus around the globe. They also provide a unique opportunity to observe virus evolution in real time. Here, we evaluate two cohorts of SARS-CoV-2 genomic sequences to identify rapidly emerging variants within structured cis-regulatory elements of the SARS-CoV-2 genome. Overall, twenty variants are present at a minor allele frequency of at least 0.5%. Several enhance the stability of Stem Loop 1 in the 5'UTR, including a set of co-occurring variants that extend its length. One appears to modulate the stability of the frameshifting pseudoknot between ORF1a and ORF1b, and another perturbs a bi-stable molecular switch in the 3'UTR. Finally, five variants destabilize structured elements within the 3'UTR hypervariable region, including the S2M stem loop, raising questions as to the functional relevance of these structures in viral replication. Two of the most abundant variants appear to be caused by RNA editing, suggesting host-viral defense contributes to SARS-CoV-2 genome heterogeneity. This analysis has implications for the development of therapeutics that target viral cis-regulatory RNA structures or sequences, as rapidly emerging variations in these regions could lead to drug resistance.

RNA genome conservation and secondary structure in SARS-CoV-2 and SARS-related viruses

As the COVID-19 outbreak spreads, there is a growing need for a compilation of conserved RNA genome regions in the SARS-CoV-2 virus along with their structural propensities to guide development of antivirals and diagnostics. Using sequence alignments spanning a range of betacoronaviruses, we rank genomic regions by RNA sequence conservation, identifying 79 regions of length at least 15 nucleotides as exactly conserved over SARS-related complete genome sequences available near the beginning of the COVID-19 outbreak. We then confirm the conservation of the majority of these genome regions across 739 SARS-CoV-2 sequences reported to date from the current COVID-19 outbreak, and we present a curated list of 30 'SARS-related-conserved' regions. We find that known RNA structured elements curated as Rfam families and in prior literature are enriched in these conserved genome regions, and we predict additional conserved, stable secondary structures across the viral genome. We provide 106 'SARS-CoV-2-conserved-structured' regions as potential targets for antivirals that bind to structured RNA. We further provide detailed secondary structure models for the 5´ UTR, frame-shifting element, and 3´ UTR. Last, we predict regions of the SARS-CoV-2 viral genome have low propensity for RNA secondary structure and are conserved within SARS-CoV-2 strains. These 59 'SARS-CoV-2-conserved-unstructured' genomic regions may be most easily targeted in primer-based diagnostic and oligonucleotide-based therapeutic strategies.

Identification of Diverse Bat Alphacoronaviruses and Betacoronaviruses in China Provides New Insights Into the Evolution and Origin of Coronavirus-Related Diseases

Outbreaks of severe acute respiratory syndrome (SARS) in 2002, Middle East respiratory syndrome in 2012 and fatal swine acute diarrhea syndrome in 2017 caused serious infectious diseases in humans and in livestock, resulting in serious public health threats and huge economic losses. All such coronaviruses (CoVs) were confirmed to originate from bats. To continuously monitor the epidemic-related CoVs in bats, virome analysis was used to classify CoVs from 831 bats of 15 species in Yunnan, Guangxi, and Sichuan Provinces between August 2016 and May 2017. We identified 11 CoV strains from 22 individual samples of four bat species. Identification of four alpha-CoVs from Scotophilus kuhlii in Guangxi, which was closely related to a previously reported bat CoV and porcine epidemic diarrhea virus (PEDV), revealed a bat-swine lineage under the genus Alphacoronavirus. A recombinant CoV showed that the PEDV probably originated from the CoV of S. kuhlii. Another alpha-CoV, α-YN2018, from Rhinolophus sinicus in Yunnan, suggested that this alpha-CoV lineage had multiple host origins, and α-YN2018 had recombined with CoVs of other bat species over time. We identified five SARS-related CoVs (SARSr-CoVs) in Rhinolophus bats from Sichuan and Yunnan and confirmed that angiotensin-converting enzyme 2 usable SARSr-CoVs were continuously circulating in Rhinolophus spp. in Yunnan. The other beta-CoV, strain β-GX2018, found in Cynopterus sphinx of Guangxi, represented an independently evolved lineage different from known CoVs of Rousettus and Eonycteris bats. The identification of diverse CoVs here provides new genetic data for understanding the distribution and source of pathogenic CoVs in China.

Gene Variations in Cis-Acting Elements between the Taiwan and Prototype Strains of Porcine Epidemic Diarrhea Virus Alter Viral Gene Expression

In 2013, the outbreak of porcine epidemic diarrhea (PED) in Taiwan caused serious economic losses. In this study, we examined whether the variations of the cis-acting elements between the porcine epidemic diarrhea virus (PEDV) Taiwan (TW) strain and the prototype strain CV777 alter gene expression. For this aim, we analyzed the variations of the cis-acting elements in the 5' and 3' untranslated regions (UTRs) between the PEDV TW, CV777, and other reference strains. We also determined the previously unidentified transcription regulatory sequence (TRS), a sequence motif required for coronavirus transcription, and found that a nucleotide deletion in the TW strain, in comparison with CV777 strain, immediately downstream of the leader core sequence alters the identity between the leader TRS and the body TRS. Functional analyses using coronavirus defective interfering (DI) RNA revealed that such variations in cis-acting elements for the TW strain compared with the CV777 strain have an influence on the efficiency of gene expression. The current data show for the first time the evolution of PEDV in terms of cis-acting elements and their effects on gene expression, and thus may contribute to our understanding of recent PED outbreaks worldwide.

Full-length genome sequences of porcine epidemic diarrhoea virus strain CV777; Use of NGS to analyse genomic and sub-genomic RNAs

Porcine epidemic diarrhoea virus, strain CV777, was initially characterized in 1978 as the causative agent of a disease first identified in the UK in 1971. This coronavirus has been widely distributed among laboratories and has been passaged both within pigs and in cell culture. To determine the variability between different stocks of the PEDV strain CV777, sequencing of the full-length genome (ca. 28kb) has been performed in 6 different laboratories, using different protocols. Not surprisingly, each of the different full genome sequences were distinct from each other and from the reference sequence (Accession number AF353511) but they are >99% identical. Unique and shared differences between sequences were identified. The coding region for the surface-exposed spike protein showed the highest proportion of variability including both point mutations and small deletions. The predicted expression of the ORF3 gene product was more dramatically affected in three different variants of this virus through either loss of the initiation codon or gain of a premature termination codon. The genome of one isolate had a substantially rearranged 5´-terminal sequence. This rearrangement was validated through the analysis of sub-genomic mRNAs from infected cells. It is clearly important to know the features of the specific sample of CV777 being used for experimental studies.

Structural and functional conservation of cis-acting RNA elements in coronavirus 5'-terminal genome regions

Structure predictions suggest a partial conservation of RNA structure elements in coronavirus terminal genome regions. Here, we determined the structures of stem-loops (SL) 1 and 2 of two alphacoronaviruses, human coronavirus (HCoV) 229E and NL63, by RNA structure probing and studied the functional relevance of these putative cis-acting elements. HCoV-229E SL1 and SL2 mutants generated by reverse genetics were used to study the effects on viral replication of single-nucleotide substitutions predicted to destabilize the SL1 and SL2 structures. The data provide conclusive evidence for the critical role of SL1 and SL2 in HCoV-229E replication and, in some cases, revealed parallels with previously characterized betacoronavirus SL1 and SL2 elements. Also, we were able to rescue viable HCoV-229E mutants carrying replacements of SL2 with equivalent betacoronavirus structural elements. The data obtained in this study reveal a remarkable degree of structural and functional conservation of 5′-terminal RNA structural elements across coronavirus genus boundaries.

•

Conservation of 5'-terminal SL1 and SL2 elements in alpha- and betacoronaviruses.

•

HCoV-229E reverse genetics data suggest critical role for SL1/2 in viral replication.

•

Replacement of SL2 in HCoV-229E with betacoronavirus equivalents is tolerated.

Software Dedicated to Virus Sequence Analysis "Bioinformatics Goes Viral"

Computer-assisted technologies of the genomic structure, biological function, and evolution of viruses remain a largely neglected area of research. The attention of bioinformaticians to this challenging field is currently unsatisfying in respect to its medical and biological importance. The power of new genome sequencing technologies, associated with new tools to handle "big data", provides unprecedented opportunities to address fundamental questions in virology. Here, we present an overview of the current technologies, challenges, and advantages of Next-Generation Sequencing (NGS) in relation to the field of virology. We present how viral sequences can be detected de novo out of current short-read NGS data. Furthermore, we discuss the challenges and applications of viral quasispecies and how secondary structures, commonly shaped by RNA viruses, can be computationally predicted. The phylogenetic analysis of viruses, as another ubiquitous field in virology, forms an essential element of describing viral epidemics and challenges current algorithms. Recently, the first specialized virus-bioinformatic organizations have been established. We need to bring together virologists and bioinformaticians and provide a platform for the implementation of interdisciplinary collaborative projects at local and international scales. Above all, there is an urgent need for dedicated software tools to tackle various challenges in virology.

The Nonstructural Proteins Directing Coronavirus RNA Synthesis and Processing

Coronaviruses are animal and human pathogens that can cause lethal zoonotic infections like SARS and MERS. They have polycistronic plus-stranded RNA genomes and belong to the order Nidovirales, a diverse group of viruses for which common ancestry was inferred from the common principles underlying their genome organization and expression, and from the conservation of an array of core replicase domains, including key RNA-synthesizing enzymes. Coronavirus genomes (~ 26–32 kilobases) are the largest RNA genomes known to date and their expansion was likely enabled by acquiring enzyme functions that counter the commonly high error frequency of viral RNA polymerases. The primary functions that direct coronavirus RNA synthesis and processing reside in nonstructural protein (nsp) 7 to nsp16, which are cleavage products of two large replicase polyproteins translated from the coronavirus genome. Significant progress has now been made regarding their structural and functional characterization, stimulated by technical advances like improved methods for bioinformatics and structural biology, in vitro enzyme characterization, and site-directed mutagenesis of coronavirus genomes. Coronavirus replicase functions include more or less universal activities of plus-stranded RNA viruses, like an RNA polymerase (nsp12) and helicase (nsp13), but also a number of rare or even unique domains involved in mRNA capping (nsp14, nsp16) and fidelity control (nsp14). Several smaller subunits (nsp7–nsp10) act as crucial cofactors of these enzymes and contribute to the emerging “nsp interactome.” Understanding the structure, function, and interactions of the RNA-synthesizing machinery of coronaviruses will be key to rationalizing their evolutionary success and the development of improved control strategies.

Coronavirus cis-Acting RNA Elements

Coronaviruses have exceptionally large RNA genomes of approximately 30 kilobases. Genome replication and transcription is mediated by a multisubunit protein complex comprised of more than a dozen virus-encoded proteins. The protein complex is thought to bind specific cis-acting RNA elements primarily located in the 5′- and 3′-terminal genome regions and upstream of the open reading frames located in the 3′-proximal one-third of the genome. Here, we review our current understanding of coronavirus cis-acting RNA elements, focusing on elements required for genome replication and packaging. Recent bioinformatic, biochemical, and genetic studies suggest a previously unknown level of conservation of cis-acting RNA structures among different coronavirus genera and, in some cases, even beyond genus boundaries. Also, there is increasing evidence to suggest that individual cis-acting elements may be part of higher-order RNA structures involving long-range and dynamic RNA–RNA interactions between RNA structural elements separated by thousands of nucleotides in the viral genome. We discuss the structural and functional features of these cis-acting RNA elements and their specific functions in coronavirus RNA synthesis.

Continuous and Discontinuous RNA Synthesis in Coronaviruses

Replication of the coronavirus genome requires continuous RNA synthesis, whereas transcription is a discontinuous process unique among RNA viruses. Transcription includes a template switch during the synthesis of subgenomic negative-strand RNAs to add a copy of the leader sequence. Coronavirus transcription is regulated by multiple factors, including the extent of base-pairing between transcription-regulating sequences of positive and negative polarity, viral and cell protein-RNA binding, and high-order RNA-RNA interactions. Coronavirus RNA synthesis is performed by a replication-transcription complex that includes viral and cell proteins that recognize cis-acting RNA elements mainly located in the highly structured 5' and 3' untranslated regions. In addition to many viral nonstructural proteins, the presence of cell nuclear proteins and the viral nucleocapsid protein increases virus amplification efficacy. Coronavirus RNA synthesis is connected with the formation of double-membrane vesicles and convoluted membranes. Coronaviruses encode proofreading machinery, unique in the RNA virus world, to ensure the maintenance of their large genome size.

Dissection of amino-terminal functional domains of murine coronavirus nonstructural protein 3

Coronaviruses, the largest RNA viruses, have a complex program of RNA synthesis that entails genome replication and transcription of subgenomic mRNAs. RNA synthesis by the prototype coronavirus mouse hepatitis virus (MHV) is carried out by a replicase-transcriptase composed of 16 nonstructural protein (nsp) subunits. Among these, nsp3 is the largest and the first to be inserted into the endoplasmic reticulum. nsp3 comprises multiple structural domains, including two papain-like proteases (PLPs) and a highly conserved ADP-ribose-1″-phosphatase (ADRP) macrodomain. We have previously shown that the ubiquitin-like domain at the amino terminus of nsp3 is essential and participates in a critical interaction with the viral nucleocapsid protein early in infection. In the current study, we exploited atypical expression schemes to uncouple PLP1 from the processing of nsp1 and nsp2 in order to investigate the requirements of nsp3 domains for viral RNA synthesis. In the first strategy, a mutant was created in which replicase polyprotein translation initiated with nsp3, thereby establishing that complete elimination of nsp1 and nsp2 does not abolish MHV viability. In the second strategy, a picornavirus autoprocessing element was used to separate a truncated nsp1 from nsp3. This provided a platform for further dissection of amino-terminal domains of nsp3. From this, we found that catalytic mutation of PLP1 or complete deletion of PLP1 and the adjacent ADRP domain was tolerated by the virus. These results showed that neither the PLP1 domain nor the ADRP domain of nsp3 provides integral activities essential for coronavirus genomic or subgenomic RNA synthesis.

IMPORTANCE

The largest component of the coronavirus replicase-transcriptase complex, nsp3, contains multiple modules, many of which do not have clearly defined functions in genome replication or transcription. These domains may play direct roles in RNA synthesis, or they may have evolved for other purposes, such as to combat host innate immunity. We initiated a dissection of MHV nsp3 aimed at identifying those activities or structures in this huge molecule that are essential to replicase activity. We found that both PLP1 and ADRP could be entirely deleted, provided that the requirement for proteolytic processing by PLP1 was offset by an alternative mechanism. This demonstrated that neither PLP1 nor ADRP plays an essential role in coronavirus RNA synthesis.